Assemble_LHS_std
| main | Tutorials | Functions | website |
Assembles the global advection-dispersion and sorption matrix of the advection dispersion equation for steady state flow field.
Version : 1.0
Author : George Kourakos
email: giorgk@gmail.com
web : http://groundwater.ucdavis.edu/msim
Date 28-Mar-2014
Department of Land Air and Water
University of California Davis
Contents
Usage
[Dglo Mglo c]= Assemble_LHS_std(p, MSH, A, V, rho_b, K_d, lambda, theta, Dm, CC, opt)
Input
p: coordinates of the Mesh nodes
MSH: Mesh structure
A: dispertivity [aL, aTH, aTV].
V: velocity
rho_b: bulk density
K_d: equilibrium distribution coefficient
lambda: decay constant
theta: porosity
Dm: Molecular diffusion coefficient
CC: Dirichlet boundary conditions
opt: option structure with the following fields : assemblemode : 'vect' for vectorize or 'nested' (just dont used the latter) capacmode : 'lumped' or 'consistent' int_ord : Gaussian intergration order. If the field does not exist a default value will be used according to element order and the dimension of the problem el_order : 'linear' Currently this is the only option for transport dim : Dimension of the problem. el_type : this is the element type.
Output
Dglo: Global advection dispersion matrix Mglo: Global sorption matrix c: The unkown concentration vector with nan values for unknown nodes scalar values for the dirichlet nodes