pdm: Flow-normal cross-sectional area PDM
In CJBarry/DNAPL: DNAPL Source Term Framework

integral of Cs*erfc(z/(2*sqrt(ax*wp0))) between 0 and hp, or else lots of smaller pools if Np isn't 1

1
2
3

C0_pool_XsectA(Cs, hp, hL)

C0_pool_Hunt88(Cs, wp0, hL, aV = wp0/1000, Np = 1L)

`Cs`	numeric []; solubility
`hp`	numeric []; pool height
`hL`	numeric []; layer thickness
`wp0`	numeric []; historic maximum of pool width
`aV`	numeric [1]; vertical dispersivity, in units of length
`Np`	integer [1]; number of pools within the layer; more pools will result in a greater average effluent concentration as no point will be so far from a pool

numeric []; C0, the historic maximum of effluent concentration averaged through the layer

It is assumed that dispersion only occurs above the pool, as pools rest on top of aquitards, so hydrodynamic dispersion below the pool should be fairly small, instead just molecular diffusion should occur into the underlying aquitard. If you want to dispense with this assumption, then double the value for Np, as the two-way dispersion from one pool at the bottom of the layer will result in the same average concentration as the one-way dispersion from two pools, one at the bottom and one in the middle.

CJBarry/DNAPL documentation built on May 6, 2019, 9:25 a.m.