28.01.2013 Views

Three-dimensional Lagrangian Tracer Modelling in Wadden Sea ...

Three-dimensional Lagrangian Tracer Modelling in Wadden Sea ...

Three-dimensional Lagrangian Tracer Modelling in Wadden Sea ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CHAPTER 2. THEORY 19<br />

side which correspond to the transport of C by turbulent fluctuations. These<br />

additional terms can be rewritten together with the diffusive transport D∆C<br />

as<br />

∂<br />

∂x<br />

�<br />

D ∂C<br />

∂x − u′ C ′<br />

�<br />

+ ∂<br />

�<br />

D<br />

∂y<br />

∂C<br />

∂y − v′ C ′<br />

�<br />

+ ∂<br />

�<br />

D<br />

∂z<br />

∂C<br />

∂z − w′ C ′<br />

�<br />

. (2.2.35)<br />

All terms <strong>in</strong> the parentheses represent mass transport. The first of these<br />

terms stands for the transport due to molecular diffusion (Fick’s law) and<br />

the second is a turbulent flux that arises by virtue of the correlation between<br />

�u ′ and C ′<br />

where �u ′ = (u ′ , v ′ , w ′ ). S<strong>in</strong>ce the molecular diffusion coefficient is<br />

usually a very small quantity (i.e. �u ′ C ′ ≫ D ∂C ), the molecular diffusion<br />

∂x<br />

terms can be neglected compared to turbulent flux. The terms denot<strong>in</strong>g<br />

the turbulent transport represent unknown quantities. Assum<strong>in</strong>g that the<br />

turbulent diffusion is a Fickian process, this closure problem is solved by<br />

<strong>in</strong>troduc<strong>in</strong>g a coefficient of eddy diffusivity for each direction. With this<br />

assumption and the notation Ax, Ay for the horizontal eddy diffusivities and<br />

ν ′ for the vertical eddy diffusivity, the closure scheme is<br />

u ′ C ′ ∂C<br />

= −Ax<br />

∂x<br />

v ′ C ′ ∂C<br />

= −Ay<br />

∂y<br />

w ′ C ′ = −ν<br />

′ ∂C<br />

. (2.2.36)<br />

∂z<br />

The coefficients Ax, Ay, ν ′ are strongly flow dependent and vary with<strong>in</strong> the<br />

flow field. With the assumption that turbulent mix<strong>in</strong>g is much stronger than<br />

molecular diffusion the transport equation for turbulent flow is<br />

∂C<br />

∂t<br />

+ ∇ · (�u C) = ∂<br />

∂x<br />

�<br />

Ax<br />

∂C<br />

∂x<br />

�<br />

+ ∂<br />

∂y<br />

�<br />

Ay<br />

∂C<br />

∂y<br />

�<br />

+ ∂<br />

∂z<br />

�<br />

ν<br />

2.2.4 Eulerian vs. <strong>Lagrangian</strong> perspective<br />

′ ∂C<br />

∂z<br />

�<br />

. (2.2.37)<br />

Fluid motion and any constituent (temperature, sal<strong>in</strong>ity, concentration of an<br />

arbitrary substance) transported by fluid motion can be described from two<br />

frames of reference, from a stationary frame (Eulerian) and from one which<br />

is mov<strong>in</strong>g along with the flow (<strong>Lagrangian</strong>). In the Eulerian perspective, the<br />

flow and its constituents are described with respect to fixed spatial positions<br />

�x = (x, y, z) and time t. Thus, they are written as e.g. �u(�x, t) and C(�x, t).<br />

The <strong>Lagrangian</strong> perspective follows the flow and traces the history of <strong>in</strong>dividual<br />

fluid particles. Unlike <strong>in</strong> the Eulerian description, spatial position is<br />

not a fixed reference but another variable of the particle. The flow variables<br />

are written with respect to time t and to a s<strong>in</strong>gle, <strong>in</strong>itial reference position,<br />

e.g. �x0 the particle position at t = 0. In this case all variables are recorded<br />

as e.g. �x(�x0, t) and C(�x0, t).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!