| 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 (i.e., 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:
where:
- Q = flow rate [L^3/T]
- k = hydraulic conductivity [L/T]
- i = hydraulic gradient (unitless)
- A = gross cross-sectional area of flow [L^2].
Darcy's law can also be expressed as:
where:
- ΔH = the head loss [L]
- L = the length of flow [L]
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:
where:
- C = conductance [L^2/T]
This results in the following general definition for conductance:
This may be represented more specifically in the following form.
where:
- t = the thickness of the material in the direction of flow [L]
- lw = the cross-sectional area perpendicular to the flow direction [L^2].
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:
where:
- Carc = conductance per unit length [(L^2/T)/L] or [L/T]
- t = the thickness of the material [L]
- w = the width of the material along the length of the arc [L]
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:
where:
- Cpoly = conductance per unit area [(L^2/T)/L^2] or [1/T]
- t = the thickness of the material [L]
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.
Parameter Factors
Beginning with GMS version 7.0, the Well, Drain, River, General head, and Stream boundary conditions have a parameter factor field that is associated with each boundary condition. The parameter factor is used with the parameter value to compute the final value of conductance for a given boundary condition (or in the case of the well package the final Q value). The parameter factor is automatically set by GMS when doing the Map→MODFLOW command. This value will be set to the length of the arc in the cell or the area of the polygon in the cell associated with the boundary condition; for a point feature the factor is set to 1.0.
See also
