phydrus.model.Model.add_solute_transport¶

Model.add_solute_transport(model=0, epsi=0.5, lupw=False, lartd=False, ltdep=False, ctola=0, ctolr=0, maxit=0, pecr=2, ltort=True, lwatdep=False, top_bc=- 1, bot_bc=0, dsurf=None, catm=None, tpulse=1)[source]¶

Method to add solute transport to the model.

Parameters

model (int, optional) – Code describing type of nonequilibrium for solute transport: 0 = equilibrium solute transport (Default) 1 = one-site sorption model (chemical nonequilibrium) 2 = two-site sorption model (chemical nonequilibrium) 3 = two kinetic sorption sites model (attachment/detachment; chemical nonequilibrium). This model is often used for particle (viruses, colloids, bacteria) transport. 4 = two kinetic sorption sites model (attachment/detachment) ( chemical nonequilibrium). Attachment coefficients are calculated using filtration theory. 5 = dual-porosity model (mobile-immobile regions; physical nonequilibrium). 6 = dual-porosity model (mobile-immobile regions) with two-site sorption in the mobile zone (physical and chemical nonequilibrium). 7 = dual-permeability model (physical nonequilibrium). 8 = dual-permeability model with either an immobile region in the matrix domain (physical nonequilibrium) or with two-site sorption in both domains (physical and chemical nonequilibrium).
epsi (float, optional) – Temporal weighing coefficient. 0.0 for an explicit scheme ( Default). 0.5 for a Crank-Nicholson implicit scheme. =1.0 for a fully implicit scheme.
lupw (bool, optional) – True if upstream weighing formulation is to be used. False if the original Galerkin formulation is to be used.
lartd (bool, optional) – True if artificial dispersion is to be added in order to fulfill the stability criterion PeCr (see Section 8.4.4), else False.
ltdep (bool, optional) – True if at least one transport or reaction coefficient (ChPar) is temperature dependent, else False. If ltdep=True, then all values of ChPar(i,M) should be specified at a reference temperature Tr=20 degrees celsius.
ctola (float, optional) – Absolute concentration tolerance [ML-3], the value is dependent on the units used (set equal to zero if nonlinear adsorption is not considered).
ctolr (float, optional) – Relative concentration tolerance [-] (set equal to zero if nonlinear adsorption is not considered).
maxit (int, optional) – Maximum number of iterations allowed during any time step for solute transport - usually 20 (set equal to zero if nonlinear adsorption is not considered).
pecr (float optional) – Stability criteria. Set to zero when lUpW=True.
ltort (bool, optional) – True if the tortuosity factor [Millington and Quirk, 1961] is to be used. False if the tortuosity factor is assumed to be equal to one.
lwatdep (bool, optional) – True if at least one degradation coefficient (ChPar) is water content dependent.
top_bc (int, optional) – Code which specifies the type of upper boundary condition 1 = Dirichlet boundary condition, -1 = Cauchy boundary condition. -2 = a special type of boundary condition for volatile solutes as described by equation (3.46).
bot_bc (int, optional) – Code which specifies the type of lower boundary condition: 1 = Dirichlet boundary condition, 0 = continuous concentration profile, -1 = Cauchy boundary condition.
dsurf (float, optional) – Thickness of the stagnant boundary layer, d [L]. Only when kTopCh=-2.
catm (float, optional) – Concentration above the stagnant boundary layer, g_atm [ML-3]. Only when kTopCh=-2.
tpulse (float, optional) – Time duration of the concentration pulse [T].