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

| main | Tutorials | Functions | website |