Advection-diffusion constitutive law for 3D solid eulerian elements (ADVE3)

This law is only used for linear pollutant transport in isotropic solids by upwind methods.
This constitutive law takes into account the advection‑dispersion in the moving fluid but also degradation, adsorption on the solid matrix and immobile fluid effect as linear phenomena. This law is used for three‑dimensional flow.

Prepro: LADV3D.F

1. Upwind parameters:
   • Optimum formulation in FRENET axes
2. Physical parameters:
   • RDM=$1 + \theta_s \rho_s p K_d / \theta_m$
   • RDIM=$1 + \theta_s \rho_s (1-p) K_d / \theta_{im}$
   • AM =$k_s.(RDM -1).+ k_m + ALM$
   • AIM=$k_s.(RDIM -1).+ k_{im} + ALIM$
   • ALM =$\alpha_\alpha/\theta_m$
   • ALIM=$\alpha_\alpha/\theta_{im}$
     With:
     • $\theta_s$=1-total porosity
     • $\theta_{m}$=effective porosity
     • $\theta_{im}$=non effective porosity
     • $k_{s,m,im}$=degradation coefficients in solid, moving and non-moving fluids
     • $K_{d}$=linear adsorption coefficient
     • $p$=part of solid surface in contact with the moving fluid
     • $\rho_s$=solid density
     • $\alpha_d$=exchange coefficient between mobile and immobile fluids

10

5