- #1

edge333

- 16

- 0

## 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

_{0}of 1.0 g CO

^{2}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:__

[itex] \frac{\delta c}{\delta t}=D \frac{\delta^{2} c}{\delta t^{2}}[/itex]

__Initial conditions:__

C = 0 for x < -L/2

C = 0 for x > +L/2

C = [itex]C_{0}[/itex] for -L/2 < x < +L/2

__Boundary conditions:__

C → 0 as x → ±∞

## The Attempt at a Solution

Based on a dirac delta function:

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

Using superposition:

[itex] 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 [/itex]

where A is the cross-sectional area of the cylinder, x

_{1}is the location of the infinitesimal mass dm

_{1}

Integrating the sources (of contaminants) from x

_{1}→-L/2 to x

_{1}→+L/2:

[itex] 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} [/itex]

[itex] c \left( x, t \right) = \frac{C_{0}}{\sqrt{\pi} } \left[

\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] [/itex]

Via substitution:

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

[itex] c \left( x, t \right) = \frac{C_{0}}{\sqrt{\pi} } \left[

\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] [/itex]

[itex] = \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] [/itex]

[itex] = \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] [/itex]

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

[itex] = \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] [/itex]

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