GMS:Conductance

From XMS Wiki
Revision as of 14:51, 10 June 2011 by Ksmith (talk | contribs) (→‎Conductance in conceptual models)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigationJump to search
MODFLOW
Pre-processing
MODFLOW Commands
Building a MODFLOW Model
Map to MODFLOW
Calibration
Packages Supported in GMS
Saving a MODFLOW Simulation
Importing MODFLOW Files
Unsupported MODFLOW Features
Run MODFLOW
Post-processing
MODFLOW Display Options
MODFLOW Post-Processing Viewing Options
Reading a MODFLOW Simulation
Tutorials
Packages
Flow: BCF6, HUF, LPF, UPW
Solvers:

DE4, GMG, NWT, PCG,

PCGN, LMG, SIP, SOR,

SMS
Other:

BAS6, BFH, CHD1, CLN,

DRN1, DRT1, EVT1, ETS1,

GAGE, GHB1, GNC, HFB1,

HUF, LAK3, MNW1, MNW2,

OUT1, RCH1, RIV1, SFR2,

STR1, SUB1, SWI2, WEL1,

UZF1

Many MODFLOW boundary condition types (ie general head, rivers, streams, drains) include a conductance parameter. MODFLOW uses the conductance to determine the amount of water that flows in or out of the model due to the boundary condition stresses.

When using a conceptual model, the manner in which the conductance term should be computed and entered depends on whether the feature object is a polygon, arc or point. Before explaining this fully, a short review of the definition of conductance is appropriate.

Conductance formulation

Darcy's law states:

File:Eq cond1.gif

where Q is the flow rate, k is the hydraulic conductivity, i represents the hydraulic gradient, and A represents the gross cross-sectional area of flow. Darcy's law can also be expressed as:

File:Eq cond2.gif

where DH represents the head loss and L represents the length of flow. Since the unknown on the right side is the head, it is convenient to group all of the other terms together and call them conductance:

File:Eq cond3.gif

This results in the following general definition for conductance:

File:Eq cond4.gif

This may be represented more specifically in the following form.

File:Eq cond5.gif

Where t represents the thickness of the material in the direction of flow, and lw represents the cross-sectional area perpendicular to the flow direction.

Conductance in conceptual models

In the case of a river boundary condition, the conductance is defined in MODFLOW as the hydraulic conductivity of the river bed materials divided by the vertical thickness (length of travel based on vertical flow) of the river bed materials, multiplied by the area (width times the length) of the river in the cell. The last term, area, is the hardest parameter to determine by hand since it varies from cell to cell.

Arcs

Fortunately, GMS can automatically calculate the lengths of arcs and areas of polygons. Therefore, when a conductance is entered for an arc, it should be entered in terms of conductance per unit length. For example, in the case of rivers, conductance should be entered as:

File:Eq cond6.gif

Where t is the thickness of the material and w is the width of the material along the length of the arc. When GMS applies the boundary condition from the arc to the grid cell, it automatically multiplies the entered value of conductance by the length of the arc that intersects the cell to create an accurate conductance value for the cell.

Polygons

For polygons, conductance should be entered in a conductance per unit area form:

File:Eq cond7.gif

Where t is the thickness of the material. When GMS converts the stress from a polygon to a grid cell, it automatically multiplies the entered value of conductance by the area of the cell that is covered by the polygon to create an appropriate conductance value for the cell. This restores the dimensional accuracy to the expression for conductance.

Points

When a general head, river, stream or drain attribute is assigned to an individual point, the conductance should be entered as a normal conductance value. This conductance is then directly assigned to the cell containing the point.

Refine Points

Refine attributes are assigned to points or nodes and are used to automatically increase the grid density around a point when the grid is constructed. Although refine attributes may be associated with any point or node, they are usually assigned in conjunction with wells.