β—€ Back to Note

Tony's Intertwingled Memex β€’ Reference Manual

Underground Contaminant Dispersion Plume Solution

Formulas & Definitions

In environmental and hydrological engineering, subterranean chemical spills threaten regional drinking water aquifers. When an industrial solvent leaks through porous underground gravel and fractured limestone, environmental remediation engineers construct 3D spatial dispersion models to estimate pollutant concentrations, design containment grouting walls, and drill monitoring recovery wells.

An environmental engineering firm in Nevada models the chemical pollutant concentration CCC (in parts per million, ppm\text{ppm}ppm) migrating from a buried ruptured storage pipe using the 3D rational function:

C(x,y,z)=100x2+y2+z2βˆ’9C(x, y, z) = \frac{100}{x^2 + y^2 + z^2 - 9}C(x,y,z)=x2+y2+z2βˆ’9100​
where x,y,zx, y, zx,y,z represent 3D spatial coordinates (measured in meters) relative to an injection test wellhead (0,0,0)(0, 0, 0)(0,0,0).

A deep monitoring sensor is installed at the coordinate position (2,1,2)(2, 1, 2)(2,1,2).

  1. State the algebraic domain condition for C(x,y,z)C(x, y, z)C(x,y,z) and write the mathematical domain DDD in set-builder notation. Describe the 3D surface that is physically and mathematically excluded from the domain.
  2. Determine the mathematical range RRR of the concentration function over its valid domain, analyzing the behavior both inside and outside the excluded boundary surface.
  3. Evaluate the contaminant concentration at the sensor location (2,1,2)(2, 1, 2)(2,1,2). Explain what this value means for environmental safety zones and groundwater monitoring.

Step-by-Step Solution:

  1. State the Denominator Condition and Identify Excluded Surface: State the Three-Variable Non-Zero Denominator Condition:

    x2+y2+z2βˆ’9β‰ 0β€…β€ŠβŸΉβ€…β€Šx2+y2+z2β‰ 9x^2 + y^2 + z^2 - 9 \ne 0 \implies x^2 + y^2 + z^2 \ne 9x2+y2+z2βˆ’9ξ€ =0⟹x2+y2+z2ξ€ =9

    • Set-Builder Notation:
      D={(x,y,z)∈R3∣x2+y2+z2β‰ 9}D = \{(x, y, z) \in \mathbb{R}^3 \mid x^2 + y^2 + z^2 \ne 9\}D={(x,y,z)∈R3∣x2+y2+z2ξ€ =9}
    • Geometric Description: The mathematical domain consists of all points in 3D space except for those lying exactly on a hollow spherical shell of radius ρ=9=3Β meters\rho = \sqrt{9} = 3\text{ meters}ρ=9​=3Β meters centered at the wellhead (0,0,0)(0, 0, 0)(0,0,0). (Defined in Domain and Range of Multivariable Functions.)
  2. Determine the Concentration Range: State the Radial Distance Coordinate Substitution: Let ρ2=x2+y2+z2\rho^2 = x^2 + y^2 + z^2ρ2=x2+y2+z2, where ρβ‰₯0\rho \ge 0ρβ‰₯0. The denominator is d=ρ2βˆ’9d = \rho^2 - 9d=ρ2βˆ’9.

    • Inside the Excluded Shell (ρ∈[0,3)\rho \in [0, 3)ρ∈[0,3)): As ρ2\rho^2ρ2 ranges from 000 to 999, the denominator d∈[βˆ’9,0)d \in [-9, 0)d∈[βˆ’9,0).
      • At the center wellhead (0,0,0)(0, 0, 0)(0,0,0), C=100βˆ’9β‰ˆβˆ’11.11Β ppmC = \frac{100}{-9} \approx -11.11\text{ ppm}C=βˆ’9100β€‹β‰ˆβˆ’11.11Β ppm.
      • As ρ→3βˆ’\rho \to 3^-ρ→3βˆ’, dβ†’0βˆ’β€…β€ŠβŸΉβ€…β€ŠCβ†’βˆ’βˆžd \to 0^- \implies C \to -\inftydβ†’0βˆ’βŸΉCβ†’βˆ’βˆž.
      • Inner range component: (βˆ’βˆž,βˆ’100/9](-\infty, -100/9](βˆ’βˆž,βˆ’100/9].
    • Outside the Excluded Shell (ρ∈(3,∞)\rho \in (3, \infty)ρ∈(3,∞)): As ρ2\rho^2ρ2 ranges from 999 to ∞\infty∞, the denominator d∈(0,∞)d \in (0, \infty)d∈(0,∞).
      • As ρ→3+\rho \to 3^+ρ→3+, dβ†’0+β€…β€ŠβŸΉβ€…β€ŠCβ†’+∞d \to 0^+ \implies C \to +\inftydβ†’0+⟹Cβ†’+∞.
      • As Οβ†’βˆž\rho \to \inftyΟβ†’βˆž, dβ†’βˆžβ€…β€ŠβŸΉβ€…β€ŠCβ†’0+d \to \infty \implies C \to 0^+dβ†’βˆžβŸΉCβ†’0+.
      • Outer range component: (0,∞)(0, \infty)(0,∞).
    • Full Mathematical Range:
      R=(βˆ’βˆž,βˆ’1009]βˆͺ(0,∞)R = \left(-\infty, -\frac{100}{9}\right] \cup (0, \infty)R=(βˆ’βˆž,βˆ’9100​]βˆͺ(0,∞)
  3. Evaluate at Sensor Coordinate (2,1,2)(2, 1, 2)(2,1,2): Compute the squared distance from the wellhead:

    ρ2=22+12+22=4+1+4=9\rho^2 = 2^2 + 1^2 + 2^2 = 4 + 1 + 4 = 9ρ2=22+12+22=4+1+4=9
    Notice that ρ2=9β€…β€ŠβŸΉβ€…β€ŠΟ=3Β meters\rho^2 = 9 \implies \rho = 3\text{ meters}ρ2=9⟹ρ=3Β meters! Substitute into the concentration expression:
    C(2,1,2)=1009βˆ’9=1000β€…β€ŠβŸΉβ€…β€ŠUndefinedΒ (Singularity)C(2, 1, 2) = \frac{100}{9 - 9} = \frac{100}{0} \implies \textbf{Undefined (Singularity)}C(2,1,2)=9βˆ’9100​=0100β€‹βŸΉUndefinedΒ (Singularity)

Meaning of the Answer & Real-Life Application:

The sensor at (2,1,2)(2, 1, 2)(2,1,2) sits precisely on the mathematical boundary sphere ρ=3Β m\rho = 3\text{ m}ρ=3Β m. In physical hydrology, this radius corresponds to the outer perimeter of an impermeable compacted grout containment barrier where the fluid flow model exhibits an asymptotic boundary singularity. Because the theoretical equation diverges at ρ=3Β m\rho = 3\text{ m}ρ=3Β m, field engineers must reposition monitoring sensors slightly beyond the containment shell (e.g., at ρ=4Β m\rho = 4\text{ m}ρ=4Β m, where C=10016βˆ’9=14.29Β ppmC = \frac{100}{16 - 9} = 14.29\text{ ppm}C=16βˆ’9100​=14.29Β ppm) to obtain reliable, finite sensor readings without digital telemetry overflow.