Gaussian Plume Simulator Back
Environmental Engineering Simulator

Gaussian Plume Simulator — Atmospheric Dispersion Model

Compute downwind dispersion of pollutants from a point source using the Gaussian plume model (Pasquill stability class D). Sliders for emission rate, wind speed, stack height and observation distance update centerline ground level concentration and spreads in real time.

Emission Conditions
Emission Q
g/s
Wind speed u
m/s
Stack height H
m
Distance x
km
Fixed Constants
Pasquill class = D (neutral)
sigma_y(x) = 0.10 x^0.92
sigma_z(x) = 0.08 x^0.92
Ground reflection at z = 0
Results
Ground concentration C
Spread sigma_y
Spread sigma_z
Ratio H/sigma_z
Plume side view — stack, dispersion, observation point
Downwind distance x vs. ground concentration C
Theory & Key Formulas

Centerline ground concentration (z = 0, with ground reflection): $$C(x,0,0) = \frac{Q}{\pi\,u\,\sigma_y\,\sigma_z}\exp\!\left(-\frac{H^2}{2\sigma_z^2}\right)$$

Pasquill D class spread (simplified, x in m): $$\sigma_y(x)=0.10\,x^{0.92},\qquad \sigma_z(x)=0.08\,x^{0.92}$$

$Q$=emission rate [kg/s], $u$=wind speed [m/s], $H$=effective stack height [m], $\sigma_y,\sigma_z$=horizontal and vertical spread [m]. The denominator is $\pi$ rather than $2\pi$ because ground reflection has been folded in.

About the Gaussian Plume Simulator

🙋
So the Gaussian plume model is how we predict where smoke from a stack will go? Where exactly does the Gaussian shape come from?
🎓
Good question. If you slice the plume across the wind, the concentration profile is assumed to be a normal (Gaussian) distribution in both the horizontal (y) and vertical (z) directions. The standard deviations are sigma_y and sigma_z, and they grow with downwind distance x. This tool uses the simplified Pasquill D class form sigma = a * x^0.92. At x = 2 km you get sigma_y ≈ 109 m and sigma_z ≈ 87.4 m. Watch the side view on the right and you can literally see the plume widen.
🙋
OK, and what happens if I raise the stack height H? How much does ground concentration drop?
🎓
That is exactly the point of stack design. The formula has an exp(-H^2/(2*sigma_z^2)) factor. When H exceeds sigma_z this factor falls very quickly. With the defaults H = 50 m and sigma_z = 87.4 m we have H/sigma_z = 0.572 and the factor is about 0.85. Push H up to 120 m and the ratio becomes 1.37 and the factor drops to 0.39. Drag the H slider down to 10 m and you can watch the ground concentration jump up dramatically.
🙋
And the wind speed u? Stronger wind should mean lower concentration, right?
🎓
Yes, your intuition is right. Concentration is inversely proportional to u, so doubling the wind speed roughly halves the concentration. The caveat is the calm-wind regime. When u drops below 1 m/s the basic Gaussian assumption fails. The slider lower bound is 0.5 m/s but for real low-wind cases you should switch to a puff model or AERMOD low-wind treatment. Press the Distance sweep button to scan x from 0.1 to 20 km and watch the ground concentration peak shift.

Physical Model and Key Equations

The Gaussian plume model assumes a steady point source whose pollutant is advected by the mean wind in the x direction and spreads laterally and vertically following Gaussian profiles. With ground reflection at z = 0 the centerline concentration is $C(x,0,0)=\dfrac{Q}{\pi u\sigma_y\sigma_z}\exp\!\left(-\dfrac{H^2}{2\sigma_z^2}\right)$. This tool uses the Pasquill D class simplified spread $\sigma_y(x)=0.10\,x^{0.92}$ and $\sigma_z(x)=0.08\,x^{0.92}$. With the defaults Q = 10 g/s, u = 5 m/s, H = 50 m and x = 2 km the result is about 56.6 μg/m³. Note that the denominator is $\pi$ rather than $2\pi$ because the ground reflection has been folded into the exponential.

Real World Applications

Environmental impact assessment (EIA): When designing new stacks for power plants or factories, regulators require that ground level concentration stays below standards across the full Pasquill A to F range. This tool gives quick D class sensitivity, the standard reference case.

Effective stack height design: Geometric stack height plus plume rise (Briggs formula) gives the effective H. Engineers minimise capital cost by finding the lowest stack that still meets WHO or national air quality limits.

Emergency response: After a chemical plant leak, a Gaussian plume estimate of the downwind impact zone can be produced in seconds, long before a full CFD run is available. It supports evacuation planning.

Urban air quality monitoring: Networks of point sources (roads, plants) are summed using Gaussian plume kernels to estimate annual mean and 1 hour peak concentrations. It remains the workhorse formula behind many regulatory codes.

Common Misconceptions and Caveats

First, the Gaussian plume model assumes a steady state and a time-averaged concentration. Short bursts and rapidly varying wind directions cannot be reproduced. When wind direction variance is large, switch to a puff model or a CFD simulation for higher fidelity.

Second, the $\sigma=a\,x^{0.92}$ form used here is a simplified Pasquill D representation. Production codes such as AERMOD and CALPUFF use class-specific coefficients and functions (Briggs open country, ASME formulation, etc.). Treat the absolute value as illustrative and the trend as the takeaway.

Third, terrain and buildings are ignored. A tall building near the stack causes downwash and the model will underestimate the near-field ground concentration. In practice apply the GEP rule (stack at least 2.5 times the tallest nearby building) and add building wake corrections.

Frequently Asked Questions

The cross-wind concentration profile (in the y and z directions) is assumed to follow a normal (Gaussian) distribution. The standard deviations sigma_y and sigma_z grow with downwind distance x, which captures the gradual dilution as the plume spreads.
5 m/s is a representative mid-latitude annual mean wind speed and a typical value for Pasquill class D (neutral). For actual assessments use site-specific wind rose data weighted by stability class frequency.
The concentration drops by exp(-y^2/(2*sigma_y^2)) as y increases. At y = sigma_y ≈ 109 m the value falls to about 0.61 of the centerline. This tool plots only the centerline (y = 0) but the lateral fall-off is part of the same Gaussian shape.
Theoretically the maximum centerline ground concentration occurs where sigma_z = H/sqrt(2). For the defaults H = 50 m under D class this is around x ≈ 0.6 km. Watch the chart on the right to see how the peak position moves when you change H.