13.07.2015 Views

Phthalates and Nonylphenols in Roskilde Fjord

Phthalates and Nonylphenols in Roskilde Fjord

Phthalates and Nonylphenols in Roskilde Fjord

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.

The vertical flux is a comb<strong>in</strong>ation of sedimentation of particles, extractionof molecular substance through sediment suction, atmosphericdeposition <strong>and</strong> freshwater sources.∂ CsN wz = - D sz ⋅ - S ⋅ R ⋅ C + C∂ zz=0q f,narrq p,narr+ ⋅ C tot,f + ⋅ Cbbs w,narr tot,deptot,p (17)The first term <strong>in</strong> Equation 17 is the boundary condition at the sedimentsurface caused by sediment “suction” of substance. The total mass balancefor the water compartment furthermore <strong>in</strong>cludes 1 st order degradation<strong>and</strong> thus becomes( q + q )∂ Cw,narr ∂= Q Cnarr w,narr f,narr p,narr⋅ +⋅ C w,narr +∂ t b ⋅ h ∂ xb ⋅ h2∂ Cw,narrDsz∂ CsD wx ⋅ - ⋅ +2∂ x R w ⋅ h ∂ zz=0CdepR h - S ⋅ R s⋅ C w,narr- k 1⋅ Cw,narrw ⋅ R w ⋅ h R wq f,narrq p,narr⋅ C tot,f +⋅ Ctot,p(18)b ⋅ h ⋅ Rb ⋅ h ⋅ RwBoundary conditions:C w,narr = C w,bred for x = 0C w,narr = C w,sea = 0 for x = L narr + 5 kmwSediment:The horizontal flux is negligible. The vertical flux arises from moleculardiffusion <strong>and</strong> advection from sediment build-upN = - Dszsz∂ Cs⋅ + S ⋅ R s ⋅ C∂ zw(19)The total mass balance for the sediment compartment is∂ C∂ ts= D Rszs2∂ Cs∂ Cs - S - k 1⋅ ⋅ ⋅ C2s(20)∂ z ∂ z RsIn the upper 5 cm the aerobic degradation rate, k 1N , is used <strong>and</strong> furtherdown the anoxic degradation rate, k 1D , is used.The two coupled 2 nd order differential Equations 18 <strong>and</strong> 20 each hav<strong>in</strong>gtwo <strong>in</strong>dependent variables, x <strong>and</strong> t for the water model <strong>and</strong> z <strong>and</strong> t for thesediment model respectively, can be solved numerically us<strong>in</strong>g a grid consist<strong>in</strong>gof discrete nodes for the (x,t) <strong>and</strong> the (z,t) system respectively,<strong>and</strong> by employ<strong>in</strong>g a forward time central space scheme (cf. Figure 36).70

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

Saved successfully!

Ooh no, something went wrong!