Air Quality & Dispersion formula solvers

Barometric Pressure with Altitude

P=P0 e−Mgz/RTP = P_0 \, e^{-Mgz/RT}

Air Quality & DispersionAtmospheric pressure at a height above a reference level, for an isothermal layer. The exponential decay that gives the atmosphere a scale height of about 8.4 km.

Briggs Buoyancy Flux

F=g vs d2(Ts−Ta)4 TsF = \frac{g \, v_s \, d^2 (T_s - T_a)}{4 \, T_s}

Air Quality & DispersionThe buoyancy flux parameter of a hot stack plume, from exit velocity, stack diameter and the temperature difference against ambient. The input to every Briggs plume-rise equation.

Briggs Plume Rise (Neutral and Unstable)

Δh=1.6 F1/3x2/3u\Delta h = \frac{1.6 \, F^{1/3} x^{2/3}}{u}

Air Quality & DispersionHow far a buoyant plume climbs above the stack at a given downwind distance, in neutral or unstable air — the Briggs two-thirds law that underlies most regulatory dispersion modelling.

Clear-Sky Temperature (Berdahl-Martin)

Tsky=Tair[ 0.711+0.56(tdp100)+0.73(tdp100)2]1/4T_{sky} = T_{air}\left[\,0.711 + 0.56\left(\tfrac{t_{dp}}{100}\right) + 0.73\left(\tfrac{t_{dp}}{100}\right)^{2}\right]^{1/4}

Air Quality & DispersionHeat TransferEffective radiating temperature of a cloudless sky, from air temperature and dew point, using the Berdahl-Martin (1984) clear-sky emissivity. This is why frost forms on a clear night and not a cloudy one: a clear sky radiates like a body tens of degrees colder than the air.

Coriolis Parameter from Latitude

f=2 Ωsin⁡φf = 2\,\Omega \sin\varphi

Air Quality & DispersionThe Coriolis parameter is not an abstract constant — it is a place on the Earth. Twice the planet's rotation rate times the sine of the latitude, zero on the equator and maximum at the poles, and it is the f that the Rossby and Ekman numbers both divide by.

Dew Point from Humidity Ratio

Td=243.04 γ17.625−γ,γ=ln⁡ ⁣pv610.94,pv=W p0.62198+WT_d = \frac{243.04\,\gamma}{17.625 - \gamma}, \quad \gamma = \ln\!\frac{p_v}{610.94}, \quad p_v = \frac{W\,p}{0.62198 + W}

HVAC & HydronicsThermodynamicsAir Quality & DispersionThe temperature at which air of a given humidity ratio and barometric pressure starts to condense, by running the vapour-pressure chain backwards through the Magnus fit.

Effective Stack Height

H=hs+ΔhH = h_s + \Delta h

Air Quality & DispersionThe height dispersion actually starts from: the physical stack plus the plume rise. Because ground-level concentration falls with the square of this height, plume rise is worth more than steel.

Ekman Number

Ek=νf L2\mathrm{Ek} = \frac{\nu}{f \, L^{2}}

Air Quality & DispersionViscosity against the Coriolis force. It is tiny almost everywhere in the atmosphere and ocean, which is why friction only matters in thin boundary layers at the top and bottom — and the thickness of those layers is exactly what this number sets.

Emission Correction to Reference Oxygen

Ccorr=Cmeas 20.9−O2,ref20.9−O2,measC_{corr} = C_{meas} \, \frac{20.9 - O_{2,ref}}{20.9 - O_{2,meas}}

Air Quality & DispersionRestates a measured stack concentration at the reference oxygen a limit is written at, so that adding dilution air can no longer make an emission look cleaner than it is. The US EPA convention, built on 20.9 % oxygen in ambient air.

Emission Rate from Stack Concentration

E=C QvE = C \, Q_v

Air Quality & DispersionThe mass of pollutant leaving a stack per unit time, from the measured concentration and the volumetric flow. This is the number a permit limit is written against and the source strength every dispersion model asks for.

Environmental Lapse Rate

Γ=T1−T2z2−z1\Gamma = \frac{T_1 - T_2}{z_2 - z_1}

Air Quality & DispersionHow fast the air cools with height, from two temperatures at two altitudes. Comparing it against the dry adiabatic rate of 9.8 °C/km is what decides whether the atmosphere is stable or unstable.

Excess Air from Flue Gas Oxygen

EA=O220.9−O2EA = \frac{O_2}{20.9 - O_2}

Air Quality & DispersionHVAC & HydronicsHow much air beyond stoichiometric is passing through a burner, read straight off the oxygen in the flue gas. The single most useful number a combustion analyser gives you, because everything about efficiency follows from it.

Gaussian Plume Ground-Level Concentration

C=Qπσyσzu e−H2/(2σz2)C = \frac{Q}{\pi \sigma_y \sigma_z u} \, e^{-H^{2}/(2\sigma_z^{2})}

Air Quality & DispersionGround-level concentration directly downwind of an elevated point source, on the plume centreline, with full reflection from the ground. The core equation of regulatory air-quality modelling.

Good Engineering Practice Stack Height

HGEP=hb+1.5LH_{GEP} = h_b + 1.5 L

Air Quality & DispersionThe stack height that just clears the turbulent wake of a nearby building, by the US EPA formula: the building height plus one and a half times its lesser dimension. It is a ceiling on credit, not a permit to build tall.

Holland Plume Rise

Δh=vsdu(1.5+2.68×10−3Pd Ts−TaTs)\Delta h = \frac{v_s d}{u}\left(1.5 + 2.68\times10^{-3} P d \, \frac{T_s - T_a}{T_s}\right)

Air Quality & DispersionPlume rise by Holland's 1953 equation, combining momentum and buoyancy in one step. Simpler than Briggs and generally conservative, which is why it survives as a screening estimate.

Humidity Ratio from Vapour Pressure

W=0.62198 pvp−pvW = 0.62198\,\frac{p_v}{p - p_v}

HVAC & HydronicsThermodynamicsAir Quality & DispersionKilograms of water vapour per kilogram of dry air, from the vapour pressure and the barometric pressure — the one humidity variable that does not move when you heat the air.

Isokinetic Sampling Rate

Qn=vsAnQ_n = v_s A_n

Air Quality & DispersionThe flow a sampling train must draw so that gas enters the nozzle at exactly the stack velocity. Get this wrong and the particulate result is biased in a direction you can predict from the sign of the error.

Maximum Ground-Level Concentration

Cmax=2QeπuH2⋅σzσyC_{max} = \frac{2Q}{e \pi u H^{2}} \cdot \frac{\sigma_z}{\sigma_y}

Air Quality & DispersionThe worst ground-level concentration a plume ever produces, wherever downwind it occurs. Falls with the square of effective stack height, which is the whole argument for building tall.

Moist Air Density at Altitude (and the 1.08 Correction)

ρ=pz (1+W)Rda T (1+1.6078 W),pz=101 325 (1−2.25577×10−5z)5.25588\rho = \frac{p_z\,(1 + W)}{R_{da}\,T\,(1 + 1.6078\,W)}, \quad p_z = 101\,325\,(1 - 2.25577 \times 10^{-5} z)^{5.25588}

HVAC & HydronicsThermodynamicsAir Quality & DispersionDensity of moist air from altitude, temperature and humidity ratio, and the correction factor it forces on the 1.08, 0.68 and 4.5 rules — which are sea-level, standard-air numbers and nothing else.

Particulate Collection Efficiency

η=Cin−CoutCin\eta = \frac{C_{in} - C_{out}}{C_{in}}

Air Quality & DispersionThe fraction of the dust entering a control device that does not leave it, from inlet and outlet loadings. The number on every baghouse, cyclone and precipitator datasheet, and the one most often quoted with too many nines.

Pasquill–Gifford Dispersion Coefficient

σ=a xb\sigma = a \, x^{b}

Air Quality & DispersionThe crosswind or vertical spread of a plume at a downwind distance, in the power-law form fitted to the Pasquill stability classes. The width term the Gaussian model needs.

ppm to mg/m³ Conversion

C=ppm⋅M24.45C = \frac{ppm \cdot M}{24.45}

Air Quality & DispersionConverts a gas concentration by volume into one by mass, using the molar volume of an ideal gas at 25 °C and 101.325 kPa. The bridge between an analyser that reads in ppm and a limit written in mg/m³.

Relative Humidity from a Sling Psychrometer

φ=pws(twb)−A p (tdb−twb)pws(tdb)\varphi = \frac{p_{ws}(t_{wb}) - A\,p\,(t_{db} - t_{wb})}{p_{ws}(t_{db})}

HVAC & HydronicsThermodynamicsAir Quality & DispersionRelative humidity from a dry-bulb and a wet-bulb reading, using the psychrometric equation with A = 6.66×10⁻⁴ K⁻¹ — the calculation behind every paper psychrometric slide rule.

Relative Humidity from Vapour Pressure

φ=pvpws\varphi = \frac{p_v}{p_{ws}}

HVAC & HydronicsThermodynamicsAir Quality & DispersionRelative humidity as what it actually is: the vapour pressure present divided by the saturation vapour pressure at the same temperature — a ratio to a target that moves whenever the air is heated.

Rossby Number

Ro=vf L\mathrm{Ro} = \frac{v}{f \, L}

Air Quality & DispersionInertia against the Coriolis force. Small Rossby numbers mean rotation runs the flow, which is why weather systems circle their low-pressure centres instead of filling them; large Rossby numbers mean rotation is irrelevant, which is why your bath does not care which hemisphere it is in.

Saturation Vapour Pressure (Magnus / Alduchov–Eskridge)

pws=610.94exp⁡ ⁣(17.625 tt+243.04)p_{ws} = 610.94 \exp\!\left(\frac{17.625\,t}{t + 243.04}\right)

HVAC & HydronicsThermodynamicsAir Quality & DispersionThe vapour pressure of water at saturation, from temperature alone, by the Alduchov–Eskridge Magnus fit — the ceiling every other psychrometric quantity is measured against.

Stack Draft Pressure (Chimney Effect)

Δp=h g (ρa−ρs)\Delta p = h \, g \, (\rho_a - \rho_s)

Air Quality & DispersionHVAC & HydronicsThe pressure a chimney generates on its own, from the height of the column and the density difference between cold outside air and hot flue gas. The same equation explains why a tall building's lobby doors are hard to open in January.

Stack Exit Velocity

vs=4Qvπd2v_s = \frac{4 Q_v}{\pi d^2}

Air Quality & DispersionHVAC & HydronicsThe speed exhaust leaves a round stack, from the volumetric flow and the inside diameter. Every plume-rise calculation starts here, and so does the check against stack-tip downwash.

Stokes Number (particle inertia)

Stk=ρp d2 v18 μ L\mathrm{Stk} = \frac{\rho_p \, d^{2} \, v}{18 \, \mu \, L}

Air Quality & DispersionThe ratio of a particle's stopping time to the time the flow takes to go round an obstacle. It is the number that decides whether a particle follows the air or hits the thing the air is flowing past — the whole of cyclone capture, impactor sampling, filtration and whether a droplet lands on your windscreen or slides around it.

Surface Temperature Depression Under a Clear Sky

Ts=Tair−qnethcT_s = T_{air} - \frac{q_{net}}{h_c}

Heat TransferAir Quality & DispersionHow far a surface settles below air temperature when it loses net radiation to a clear sky and gains heat back only by convection. This is why frost appears on a windscreen while the porch thermometer reads a few degrees above freezing.

Wind Speed at Height (Power Law)

u2=u1(z2z1)pu_2 = u_1 \left(\frac{z_2}{z_1}\right)^{p}

Air Quality & DispersionExtrapolates a wind measurement to another height using the power-law profile. Anemometers sit at 10 m and stacks do not, so this step precedes every plume-rise calculation.