GMS:MF6 MDT Package: Difference between revisions

From XMS Wiki
Jump to navigationJump to search
 
(37 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{MODFLOW 6 Beta}}
{{MODFLOW 6 Beta}}
{{MODFLOW 6 Links}}
{{MODFLOW 6 Links}}
[[File:MDTPackage.png|thumb|350 px|Example of the ''MDT Package'' dialog.]]
[[File:GMS MDT Matrix Diffusion.png|thumb|450 px|Example of the ''MDT Package'' dialog.]]
The Matrix Diffusion Transport package (MDT) works with MODFLOW-USG Transport. The MDT package allow existing flow and chemical transport models of sites to include a full accounting of matrix diffusion effects. The package works with the BCT process.
The Matrix Diffusion Transport package (MDT) works with MODFLOW 6. The MDT package allow existing flow and chemical transport models of sites to include a full accounting of matrix diffusion effects. The package works with the BCT process.
:'''''Note:''''' the MDT package is not an official package released by USGS. The MDT package for MODFLOW 6 provided with GMS is still being developed and may change in the future. Newer versions of the package may not be compatible with older projects.


The ''MDT Package'' dialog is divided into two parts: one on the left which shows the main sections, and one on the right which shows the inputs for the current section.
The ''MDT Package'' dialog is divided into two parts: one on the left which shows the main sections, and one on the right which shows the inputs for the current section.


==Comments==
==COMMENTS==
Comment lines can be added and will appear at the top of the MDT file. Use the '''Add Row''' [[File:Insert Row Icon.svg|16 px]] button to create comments. Comments can be removed using the '''Delete Row''' [[File:Delete Row Icon.svg|16 px]] button.
Comment lines can be added and will appear at the top of the MDT file.  


==Variables==
==OPTIONS==
*''Cell-by-cell mass flux unit (IMDTCB)'' – The mass flux between the fracture and the matrix for every cell.
<!--*''Print_Input'' &ndash; Print input to listing file.
*''Immobile domain concentration unit (IMDTCF)'' &ndash; enter a value to activate the immobile domain concentration unit.
*''Print_Concentration'' &ndash; Print concentration to listing file.-->
*''SAVE_FLOWS'' &ndash; Save flows to budget file.
*''PRINT_FLOWS'' &ndash; Print flows to listing file.
*''SORPTION'' &ndash; Activates sorption. Requires values entered for bulk density and distance coefficient.
*''FIRST ORDER DECAY'' &ndash; first order decay will occur if this option is on.
* ''TSHIFTMD'' &ndash; Enter the time elapsed since matrix diffusion began.
* ''MD FILEOUT'' &ndash; Enter a file name for writing AIOLD1MD and AIOLD2MD coefficients to a binary file.
<!--*''Concentration Fileout'' &ndash; Save concentration to file.
*''Budget Fileout'' &ndash; Save budget to file.
*''Hydraulic conductivity for fracture only (FRANHK)'' &ndash; Indicates that the hydraulic conductivity and storage terms (for transient simulations) which are input are only for the fracture (mobile) domain. Otherwise, it is assumed that the hydraulic conductivity and storage input for the flow simulation are effective hydraulic conductivity and storage terms for the medium, representative of the entire volume.
*''Hydraulic conductivity for fracture only (FRANHK)'' &ndash; Indicates that the hydraulic conductivity and storage terms (for transient simulations) which are input are only for the fracture (mobile) domain. Otherwise, it is assumed that the hydraulic conductivity and storage input for the flow simulation are effective hydraulic conductivity and storage terms for the medium, representative of the entire volume.
*''Darcy flux hydraulic conductivity for fracture only (FRADARCY)'' &ndash; Indicates that the Darcy flux should be computed only for the fracture domain by use of this hydraulic conductivity value. Thus, the input hydraulic conductivity is divided by the fracture volume fraction. Also, the storage terms (for transient simulations) which are input are only for the fracture (mobile) domain.
*''Darcy flux hydraulic conductivity for fracture only (FRADARCY)'' &ndash; Indicates that the Darcy flux should be computed only for the fracture domain by use of this hydraulic conductivity value. Thus, the input hydraulic conductivity is divided by the fracture volume fraction. Also, the storage terms (for transient simulations) which are input are only for the fracture (mobile) domain.
*''Shift matrix diffusion start time (TSHIFTMD)'' &ndash; indicates that the MODFLOW start time (t = 0) is not the starting time of matrix diffusion, and the subsequent real variable ''tshift'' is the value of time already elapsed since matrix diffusion began.
*''Shift matrix diffusion start time (TSHIFTMD)'' &ndash; indicates that the MODFLOW start time (t = 0) is not the starting time of matrix diffusion, and the subsequent real variable.
*''Time elapsed since matrix diffusion began (tshift)'' &ndash; Enter the value of time already elapsed since matrix diffusion began.
-->


==Aquifer Properties==
==GRIDDATA==
*''Property Name'' &ndash; List of available properties.
*''Griddata'' &ndash; List of available properties.
**''MDFLAG'' — The matrix diffusion type flag at any location.
**''MD_TYPE_FLAG'' — The matrix diffusion type flag at any location.
**''VOLFRACMD'' &ndash;  The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
**''MD_FRACTION'' &ndash;  The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
**''PORMD'' &ndash;  The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
**''MD_POROSITY'' &ndash;  The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
**''RHOBMD'' &ndash; The dry bulk density of the porous matrix of the matrix domain at any location
**''BULK_DENSITY'' &ndash; The dry bulk density of the porous matrix of the matrix domain at any location
**''DIFFLENMD'' &ndash;  The diffusion length for diffusive transport within the matrix domain at any location (only used for embedded heterogeneity).  
**''MD_DIFF_LENGTH'' &ndash;  The diffusion length for diffusive transport within the matrix domain at any location (only used for embedded heterogeneity).  
**''TORTMD'' &ndash;  The tortuosity of the matrix domain at any location.
**''MD_TORTUOSITY'' &ndash;  The tortuosity of the matrix domain at any location.
*''Use Constant'' &ndash; Turn on to enter a constant value.
<!--**''KDMD'' &ndash; The adsorption coefficient (Kd) of a contaminant species in the immobile domain at any location. -->
*''Constant Value'' &ndash; Enter a constant value for the property.
**''MD_DIST_COEFF'' &ndash; The distribution coefficient of each species component within the matrix domain at any location.
*''Array Editor'' &ndash; Clicking the '''Edit''' will open an [[GMS:MODFLOW Array Editor|array editor dialog]] for the specific property.
**''MD_DECAY'' &ndash; The first-order decay coefficient in water and on soil (1/time) in the matrix domain at any location.
 
<!--**''YEILDMD'' &ndash; The yield coefficient for chain decay at any location. The value indicates mass of species generated by decay of the parent component. -->
==Species Properties==
**''MD_DIFF_COEFF'' &ndash; The diffusion coefficient of each species component within the matrix domain at any location.
*''Species ID''
*''Property Name'' &ndash; List of available properties.
**''KDMD'' &ndash; The adsorption coefficient (Kd) of a contaminant species in the immobile domain at any location.  
**''DECAYMD'' &ndash; The first-order decay coefficient in water and on soil (1/time) in the matrix domain at any location.
**''YEILDMD'' &ndash; The yield coefficient for chain decay at any location. The value indicates mass of species generated by decay of the parent component.  
**''DIFFMD'' &ndash; The diffusion coefficient of each species component within the matrix domain at any location.
**''AIOLD1MD'' &ndash; The concentration integral (equation 8) of each species in an element associated with matrix diffusion into an infinite acting adjacent aquitard (MDFLAG = 1, 3, 4, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.  
**''AIOLD1MD'' &ndash; The concentration integral (equation 8) of each species in an element associated with matrix diffusion into an infinite acting adjacent aquitard (MDFLAG = 1, 3, 4, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.  
**''AIOLD2MD'' &ndash; The concentration integral (equation 8) of each species in an element associated with matrix diffusion in low permeability zones that are embedded in the volume element (MDFLAG = 2, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.
**''AIOLD2MD'' &ndash; The concentration integral (equation 8) of each species in an element associated with matrix diffusion in low permeability zones that are embedded in the volume element (MDFLAG = 2, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.
*''Layered'' &ndash; Turn on and enter the layer value to apply values to specific layers.
*''Use Constant'' &ndash; Turn on to enter a constant value.
*''Use Constant'' &ndash; Turn on to enter a constant value.
*''Constant Value'' &ndash; Enter a constant value for the property.
*''Constant Value'' &ndash; Enter a constant value for the property.
*''Array Editor'' &ndash; Clicking the '''Edit''' will open an [[GMS:MODFLOW Array Editor|array editor dialog]] for the specific property.
*'''Dataset to Array''' &ndash; Opens a ''Select Dataset'' dialog where a dataset can be selected for the array.
 
==Conversion from MODFLOW-USG Transport==
When converting a MODFLOW-USG Transport model to MODFLOW 6, use the following table when converting the MDT package:
{|- border="1"
|-
|colspan="2" width="250" style="background: #E1E1E1;"|'''MODFLOW-USG Transport'''
|colspan="2" width="250" style="background: #E1E1E1;"|'''MODFLOW 6'''
|-
|style="background: #E1E1E1;"|'''Package'''
|style="background: #E1E1E1;"|'''Variable'''
|style="background: #E1E1E1;"|'''Package'''
|style="background: #E1E1E1;"|'''Variable'''
|-
|PCB
|
|CNC
|
|-
|BCT
|DIFFNC
|DSP
|DIFFC
|-
|BCT
|PRSITY
|MST
|POROSITY
|-
|BCT
|BULKD
|MST
|BULK_DENSITY
|-
|BCT
|ADSORB
|MST
|DISTCOEF
|-
|MDT
|MDFLAG
|MDT
|MD_TYPE_FLAG
|-
|MDT
|VOLFRACMD
|MDT
|MD_FRACTION
|-
|MDT
|PORMD
|MDT
|MD_POROSITY
|-
|MDT
|RHOBMD
|MDT
|BULK_DENSITY
|-
|MDT
|DIFFLENMD
|MDT
|MD_DIFF_LENGTH
|-
|MDT
|TORTMD
|MDT
|MD_TORTUOSITY
|-
|MDT
|DIFFMD
|MDT
|MD_DIFF_COEFF
|}


==Related Topics==
==Related Topics==
Line 48: Line 125:
{{Navbox GMS}}
{{Navbox GMS}}
[[Category:MODFLOW 6|MDT]]
[[Category:MODFLOW 6|MDT]]
[[Category:MODFLOW 6 Dialogs|MDT]]

Latest revision as of 18:35, 23 August 2023

MODFLOW 6
Models & Tools
GWF Model
GWT Model
Cell Properties Dialog
Zone File
Packages
Flow: GNC, HFB, NPF
GWF: BUY, CHD, CSUB,
DRN, EVT, GHB,
LAK, MAW, MVR,
OBS, RCH,RIV,
SFR, STO, UZF,
WEL
GWT: ADV, CNC, DSP
FMI, IC, IST, LKT,
MDT, MST, MVT,
MWT, OBS, SFT,
SRC,SSM, UZT
Other DIS, DISU, DISV,
IMS, OC, TDIS,
PEST

Example of the MDT Package dialog.

The Matrix Diffusion Transport package (MDT) works with MODFLOW 6. The MDT package allow existing flow and chemical transport models of sites to include a full accounting of matrix diffusion effects. The package works with the BCT process.

Note: the MDT package is not an official package released by USGS. The MDT package for MODFLOW 6 provided with GMS is still being developed and may change in the future. Newer versions of the package may not be compatible with older projects.

The MDT Package dialog is divided into two parts: one on the left which shows the main sections, and one on the right which shows the inputs for the current section.

COMMENTS

Comment lines can be added and will appear at the top of the MDT file.

OPTIONS

  • SAVE_FLOWS – Save flows to budget file.
  • PRINT_FLOWS – Print flows to listing file.
  • SORPTION – Activates sorption. Requires values entered for bulk density and distance coefficient.
  • FIRST ORDER DECAY – first order decay will occur if this option is on.
  • TSHIFTMD – Enter the time elapsed since matrix diffusion began.
  • MD FILEOUT – Enter a file name for writing AIOLD1MD and AIOLD2MD coefficients to a binary file.

GRIDDATA

  • Griddata – List of available properties.
    • MD_TYPE_FLAG — The matrix diffusion type flag at any location.
    • MD_FRACTION – The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
    • MD_POROSITY – The mobile fraction at any location, i.e., the fraction of the total space that is occupied by the fracture/high permeability domain.
    • BULK_DENSITY – The dry bulk density of the porous matrix of the matrix domain at any location
    • MD_DIFF_LENGTH – The diffusion length for diffusive transport within the matrix domain at any location (only used for embedded heterogeneity).
    • MD_TORTUOSITY – The tortuosity of the matrix domain at any location.
    • MD_DIST_COEFF – The distribution coefficient of each species component within the matrix domain at any location.
    • MD_DECAY – The first-order decay coefficient in water and on soil (1/time) in the matrix domain at any location.
    • MD_DIFF_COEFF – The diffusion coefficient of each species component within the matrix domain at any location.
    • AIOLD1MD – The concentration integral (equation 8) of each species in an element associated with matrix diffusion into an infinite acting adjacent aquitard (MDFLAG = 1, 3, 4, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.
    • AIOLD2MD – The concentration integral (equation 8) of each species in an element associated with matrix diffusion in low permeability zones that are embedded in the volume element (MDFLAG = 2, 5, 6, or 7). This value is used to restart a simulation with prior matrix diffusion.
  • Layered – Turn on and enter the layer value to apply values to specific layers.
  • Use Constant – Turn on to enter a constant value.
  • Constant Value – Enter a constant value for the property.
  • Dataset to Array – Opens a Select Dataset dialog where a dataset can be selected for the array.

Conversion from MODFLOW-USG Transport

When converting a MODFLOW-USG Transport model to MODFLOW 6, use the following table when converting the MDT package:

MODFLOW-USG Transport MODFLOW 6
Package Variable Package Variable
PCB CNC
BCT DIFFNC DSP DIFFC
BCT PRSITY MST POROSITY
BCT BULKD MST BULK_DENSITY
BCT ADSORB MST DISTCOEF
MDT MDFLAG MDT MD_TYPE_FLAG
MDT VOLFRACMD MDT MD_FRACTION
MDT PORMD MDT MD_POROSITY
MDT RHOBMD MDT BULK_DENSITY
MDT DIFFLENMD MDT MD_DIFF_LENGTH
MDT TORTMD MDT MD_TORTUOSITY
MDT DIFFMD MDT MD_DIFF_COEFF

Related Topics