1-D Bounded, Distributed Diffusion of Contaminant

[edge333]

Homework Statement

I'm having some difficulty deriving the equation for a concentration of CO2 as a function of length and time. Ultimately I end up with an equation that includes the summation of two error function terms that appear to have incorrect signs.

Given:

A cylinder of infinite length has a mass input, M₀ of 1.0 g CO² added over a width, L, of 5 cm at t = 0. The cylinder has a diameter of 5 cm. Assume a fluid temperature of 25°C.

Find:

The time to reach a concentration of 2 ppm at x = 100 cm for:

1.) Molecular diffusion in air
2.) Molecular diffusion in water
3.) If the fluid was turbulent, how would that affect your answers above? Substitute a value of 1 cm2/s for D as a reflection of uniformly generated (isotropic and homogenous) turbulence in water. How does this time compare to part 2.) above?

Assumptions:
Symmetry (d/dz = d/dy = 0)
No flow (u = v = w = 0)
Conservative tracer ( r = 0)

Homework Equations

At x = 0, the mass is evenly distributed between -L/2 and L/2

Governing Equation:

\frac{\delta c}{\delta t}=D \frac{\delta^{2} c}{\delta t^{2}}

Initial conditions:

C = 0 for x < -L/2
C = 0 for x > +L/2
C = C_{0} for -L/2 < x < +L/2

Boundary conditions:

C → 0 as x → ±∞

The Attempt at a Solution

Based on a dirac delta function:

c \left( x, t \right) = \left[ \frac{M^{0}}{A \sqrt{4 \pi D t}} \right] exp \left[ - \frac{x}{4 D t} \right]

Using superposition:

dc = \frac{C_{0} A dx_{1}}{A \sqrt{4 \pi D t}} exp \left( - \frac { \left( x - x_{1} \right )^{2} }{4 D t} \right) dt

where A is the cross-sectional area of the cylinder, x₁ is the location of the infinitesimal mass dm₁

Integrating the sources (of contaminants) from x₁→-L/2 to x₁→+L/2:

c \left( x, t \right) = \int ^{L/2}_{-L/2} \frac{C_{0} }{ \sqrt{4 \pi D t}} exp \left( - \frac { \left( x - x_{1} \right )^{2} }{4 D t} \right) dx_{1}

c \left( x, t \right) = \frac{C_{0}}{\sqrt{\pi} } \left[ <br /> \int ^{L/2}_{-\infty} \frac{1}{ \sqrt{4 D t}} exp \left( - \frac { \left( x - x_{1} \right )^{2} }{4 D t} \right) dx_{1} - \int ^{-L/2}_{-\infty} \frac{1}{ \sqrt{4 D t}} exp \left( - \frac { \left( x - x_{1} \right )^{2} }{4 D t} \right) dx_{1} \right]

Via substitution:

where \eta = \frac{ \left( x - x_{1} \right) }{ \sqrt{4 D t}} and d \eta = - \frac{dx_{1} }{ \sqrt{4 D t} }

c \left( x, t \right) = \frac{C_{0}}{\sqrt{\pi} } \left[ <br /> \int ^{\infty}_{\frac{x-L/2}{\sqrt{4 D t}}} - exp \left( - \eta^{2} \right) d \eta - \int ^{\infty}_{\frac{x+L/2}{\sqrt{4 D t}}} - exp \left( - \eta^{2} \right) d \eta \right]

= \frac{C_{0}}{ 2 } \left[ - \left( 1 - erf \left( \frac{x-L/2}{\sqrt{4 D t}} \right) \right) + \left( 1 - erf \left( \frac{x+L/2}{\sqrt{4 D t}} \right) \right) \right]

= \frac{C_{0}}{ 2 } \left[ erf \left( \frac{x-L/2}{\sqrt{4 D t}} \right) - erf \left( \frac{x+L/2}{\sqrt{4 D t}} \right) \right]

From what I have seen, this is not correct though. The two error function terms should be switched such that:

= \frac{C_{0}}{ 2 } \left[ erf \left( \frac{x+L/2}{\sqrt{4 D t}} \right) - erf \left( \frac{x-L/2}{\sqrt{4 D t}} \right) \right]

I'm wondering if I switched the signs somewhere or I made a mistake in the integration, perhaps the domain.

 

[Chestermiller]

The reference result you gave (with the signs switched) is correct. You must have made a mistake somewhere. It is much simpler to do the integration in a little different way than you have done. When you split it into two integrals, take one integral from -L/2 to 0, and the other integral from 0 to +L/2. The error functions with their correct signs will emerge virtually immediately from this.