SEAWAT 4 Tutorial
Contents
1. xl Settings Values Conjugate Gradient Solver fig Run Translate Layer Property Flow Packa Y Run Translate Run Translate Solver Basic Package Layer Property Package Well Package Oc Package E Specified Head Package Translate CUNIT 20 Extension CHD LUNIT 30 r lv Run Translate Run Translate Run Translate Run Translate ha Run Translate Discretization Package FTL Package Calb Package Nodry Package Run Translate Run Translate Basic Package Advection Package Dispersion Package Run Translate Sink Sources Package Run Translate Chemical Reactions Package Run Translate GCG Solver Package Run Translate Viscosity Package Run Translate hha Basic Transport Package Reset OK Cancel 22 SEAWAT 4 Tutorial Example Problem Beside Layer Property Package Layer Property Flow Package from the drop down list box r OK E Translate amp Run The VMEngine dialog will open which will track the progress of the simulation Tyres Ce oes and eel Belo ek aga ar banal vet wl 00 yenannples Twelve ea T ata titel Fils Fun an Hip E E e Fa teh to active enges SEAWAT H7 T Shea panad 1 Timon sings 14 BA HRALA pak fannie See Sre raO O ji Max head change v eration number l Properties Valu P Has utes Berations2. On your Windows desktop you will see an icon for Visual MODFLOW Visual MODFLOW to start the program Se EF To open the model er File from the top menu bar or Open An Open Visual MODFLOW file window will appear Navigate to the following directory on your computer C My Documents Visual MODFLOWTutorials Seawat_Heat i Seawat4_Tutorial vmf s Open Next you will take a quick look at some aspects of the model that have already been prepared for you After that we will create a new SEAWAT variant for simulating heat and salt transport Exploring the Model Grid and Properties In this section you will explore some aspects of the model that have already been defined for you such as the model grid and properties This information can be viewed in the Input window Input from the Visual MODFLOW main menu Exploring the Model Grid and Properties 3 File Grid Wells Properties Boundaries Particles ZBud Tools Help View Column View Layer ae Grid s7 nactive Cells gt Edit Extents Smoothing mport Elevation Assign Elevation Contouring Z Row I Column J F1 Fa Eg FY FB 5 FE 5 Eg F8 verti E2 Queri F i Layer K f Save nap and A outl A Pan Kjerag Pay E nen Grid edit module As mentioned in the description the grid used in this model is 2 dimensional The grid and cell dimensions are shown in the following table Input Parameter Value Number of Columns NCOL Number of Rows NROW Table 1 Mod
3. Defining Model Boundary Conditions 13 File Help Description Constant Head 0 pa ae amp Assign to appropriate layer J Linear gradient Edit selected rowfs or column Start Time day Stop Time day Start Time Head Stop Time Head Density Options Density kg m3 m m 0 400000 0 OK Cancel In the Stop Time day field type 200000 In the Start Time Head field type 0 In the Stop Time Head field type 0 In the Density Options drop down list box i No Conversion In the Density field type 0 a OK to close the dialog box Now that the flow conditions have been assigned at each boundary you can now proceed to define the point source boundary conditions We will start at the right boundary as you should already be zoomed into that portion of the model grid To assign a point source boundary condition G Boundaries gt Point Source from the menu bar er Yes at the Save warning message s Assign gt Single from the side bar menu The following dialog box will appear on your screen 14 SEAWAT 4 Tutorial Example Problem Assign Point Source Code fi m Time Schedule Copy from flow C Define new Start Time day Stop Time day Salt mg L Temperature C With the dialog still opened s Each cell in Column 100 You will notice that when a cell is clicked the color of the cell turns pink Once every cell in column 1 has been sele
4. MITER a Mae miner irar seen te to DC En gee 1 es change chen HELO Roina aban ROLOSE FELAK 3 _Weper boured of cette HPE Freed erbeste MPRUL 2 Lecter WAM 1 Upper oi Sice sring 1AE uf ne pater a Te pilles kemar ET er am fee le ie t ic alg le E re vi gi ith pite annie T tar te a Ai tere de es eee et E jE 1 rie Bi 5 i att HE i tee j E rant mn Eh ais atte fing ita zig a elit A d te g 7 eter ve ir B lettre T ee we ee tH iad Be i iE a dd i 14137 14157 14147 1415 Tata 1 attr ERE a ts as erabon iunie NOTE SEAWAT requires longer run times as the complexity of the model increases On a machine with a PentiumlI processor this model requires approximately five minutes to run NOTE A SEAWAT model takes about 2 4 times longer to run than a MODFLOW model of same dimensions due to the coupled flow and transport Once the simulation has finished running ee Close Viewing Output To view the output er Output from the Main Menu The following view will be presented Viewing Output 23 Fils Maps Grapte Tools He AD y ao aa 2h aL rip rrit a m E anini meie mr B iB Gi 4 MARAEA Tegght enfell qnl bens wrlh mpopmbms odp maas ch e Maps gt Contouring gt Concentration to activate the concentration contours Both Salt and Heat transport were simulated in the model By defau
5. Options er Solve Implicitly Use GCG Solver to unselect the checkbox For the Advection drop down list box TVD option ee Under Dispersion Reaction and Sink Sources s Solve Implicitly Use GCG Solver option For the Maximum Step Size field type 0 For the Multipler field type 1 For the Min sat thickness field type 5e 02 For the Courant number field type 1 To accept all modifications and close the Simulation Scenario window OK button ee Defining Simulation Run Settings Setting Flow Time Steps From main menu s SEAWAT gt Flow gt Time Steps The Time Step Options dialog will appear on your screen For the Time Steps field type 10 default For the Multiplier field type 1 xf aooo pooo Start day Stop day Time steps Multiplier Steady state 0 200000 10 1 OK Cancel e OK button Configuring Layer Settings From main menu s SEAWAT gt Flow gt Layers The Layer Settings dialog will appear on your screen For each layer change the layer type from 3 Confined Unconfined variable S T to 0 Confined Constant S T To save time you can select all the rows by clicking and dragging your mouse along the left side of the grid Once all the rows are selected you can select 0 Confined Constant S T from the drop down list box at the top of the window 20 SEAWAT 4 Tutorial Example Problem 00 oo 00 00 00 ne 00 O0 Harmonic mean 2 a a oe a O Con
6. SEAWAT 4 Tutorial Example Problem Introduction About SEAWAT 4 SEAWAT 4 is the latest release of SEAWAT computer program for simulation of three dimensional variable density transient ground water flow in porous media SEAWAT was originally designed by combining a modified version of MODFLOW 2000 and MT3DMS into a single computer program SEAWAT 4 includes all of the capabilities of previous versions with the addition of some new features including The ability to simulate heat and solute transport simultaneously The ability to simulate the effect that fluid viscosity variations have on resistance to ground water flow using the Viscosity VSC Package Simulate fluid density as a function of temperature and one or more MT3DMS species and optionally a pressure term You can also specify one of several non linear equations to represent the dependence of viscosity on temperature For the Time Variant Constant Head CHD package auxiliary variables can now be used to designate the fluid density value associated with the prescribed head Sample Applications Some examples of applications for SEAWAT 4 include Background Introduction Simulating the transport of heat and salinity in coastal aquifers and thermally driven convection in deep aquifers ASR Storage of freshwater in saline aquifers simulate buoyancy effects and the affect of groundwater temperature on the fluid viscosity Geothermal investigations The heat balan
7. ature front For each time step the temperature concentration contours and color shading should look similar to the ones shown below SEAWAT 4 Tutorial Example Problem Viewing Output PHIO1 5 000 Days IRU UUE BE WAS SIOROUO DENA yat kia A tr 29 30 ORON BES at s AR art 200 000 Days The contours should appear similar to the results published in the USGS SEAWAT 4 User s Manual shown below time 30 000 days time 5 000 days SEAWAT 4 Tutorial Example Problem Please note that your output results may vary slightly from the results produced in the USGS USer s Manual This is because some parts of this model have been modified to allow for shorter run times Viewing Output 31 32 SEAWAT 4 Tutorial Example Problem
8. can proceed to translate and run the SEAWAT model However first you must configure the appropriate simulation settings er File gt Main Menu eo Run from the Visual MODFLOW main menu bar e SEAWAT gt Simulation Scenario VDF Settings First you will define the settings for the Variable Density Flow VDF package Simulation Scenario i x VDF Settings vsc Settings Flow Settings Transport Settings Coupling Option of YDF and IMT Density scenario E densiy from sat only Explicitly coupled YDF and IMT Mode E Calculate density from multiple species Implicitl VDF IMT C Uncoupled simulation using pre calculated density Mode D mpiciNy couped VDE and Mic Mote EJ M Density specification IDENSE e density a T Maximum number of implicitly coupling iterations for YDF and IMT 2 0 Fluid density convergence criterion R i e 4 Kg m 3 Use convergence criterion NSWTCF Intermodal density calculation algorithm MFNADVFD Centralin space algorithm i Density fram concentration UEN file Upstream weighted algorithm Variable density water table correction I WTABLE C Applied Not applied Initial Time Step Initial Time Step fo 01 day Internodal density calculation algorithm MFNADVFD Bact Next gt OK Cancel Help For the VDF Package most of the default values are acceptable except for the following Unde
9. ce calculation and the impact of groundwater heat pumps on the groundwater temperature Affect of volcanic activities on groundwater and the water table Fate of wastewater injected into deep wells and impact on the aquifer This tutorial assumes that you are already familiar with the Visual MODFLOW interface and with the process of building a groundwater model If you are not familiar with Visual MODFLOW it is recommended that you work through the Visual MODFLOW tutorial first This tutorial illustrates a common application for SEAWAT and is based on the Example Problem Case 7 published in User s Guide to SEAWAT 4 Please note that output results may vary slightly from published results as this model has been scaled down to allow for shorter run times The problem consists of a two dimensional cross section of a confined coastal aquifer initially saturated with relatively cold seawater at a temperature of 5 C Warmer freshwater with a temperature of 25 C is injected into the coastal aquifer along the left boundary to represent flow from inland areas The warmer freshwater flows to the right where it discharges into a vertical ocean boundary The ocean boundary is represented with hydrostatic conditions based on a fluid density calculated from seawater salinities at 5 C No flow conditions are assigned to the top and bottom boundaries Reference Langevin C D Thorne D T Jr Dausman A M Sukop M C and Guo Weixing 2008 SEAWAT V
10. convenience the well boundary has already been assigned for you To view the well boundary follow the steps below If you are currently zoomed in e Zoom Out located at the bottom tool bar When you see the full grid on your screen er Wells gt Pumping from the menu bar G Yes at the Save message At the left boundary you will see that a well boundary condition has been defined From the left side bar e Database button The Edit Well window will appear on your screen Here you can see for PWOOI a pumping rate injection of 10 m d has been assigned to the entire column resulting in arate 0 1 m day at each individual cell G OK to close the Edit Well window Along the right boundary we will assign a Constant Head flow boundary condition To make it easier to assign this boundary condition you may want to zoom into the right boundary portion of the model grid You can do this by using the Zoom In button located at the bottom toolbar The boundary condition will be assigned to all cells in column 100 Once this column is visible on your screen proceed with the steps below e Boundaries gt Constant Head from the menu bar e Assign gt Line from the side menu bar F Inside cell Row 1 Column 100 Layer 1 s Inside cell Row 1 Column 100 Layer 100 e Right click in cell Row 1 Column 100 Layer 100 Upon right clicking all cells in column 100 will be colored pink and the following dialog box will appear on your screen
11. cted colored pink type 35000 in the Salt field of the Assign Point Source dialog box e OK in the Point Source dialog box Next we will assign a point source boundary condition to the left boundary G Zoom Out from the bottom toolbar To make it easier to assign this boundary condition you may want to zoom into the left boundary portion of the model grid You can do this by using the Zoom In button located at the bottom toolbar The boundary condition will be assigned to all cells in column 1 Once this column is completely visible on your screen proceed with the steps below P Boundaries gt Point Source from the menu bar s Yes at the Save warning message s Assign gt Single from the side bar menu The following dialog box will appear on your screen Assign Point Source Code fi m Time Schedule Copy from flow C Define new Start Time day Stop Time day Salt mg L Temperature C of 20000 0 25 Import OK Cancel With the dialog still opened Defining Model Boundary Conditions 15 er Each cell in Column 1 You will notice that when a cell is clicked the color of the cell will turn pink Once every cell in column 200 has been clicked colored pink type 25 in the Temperature field of the Assign Point Source dialog box e OK in the Point Source dialog box Defining Simulation Run Settings Now that the grid properties and boundaries have been defined you
12. el grid and cell dimensions The current view is set to layer view planar To view the model grid from a cross section view er View Row from the side menu On the grid SEAWAT 4 Tutorial Example Problem G Row 1 the only row Due to the ratio of length to width of the model grid you may need to zoom in on the model to successfully select the row Once selected your view will look similar to the one shown in the image shown below Next we will look at the assigned property values using the Cell Inspector tool e Tools from the main menu Cell Inspector The Cell Inspector window will appear on your screen er Options tab Properties node expand this node Kx Ky Kz Specific Yield 47 4 4 4 4 Specific Storage oe Effective Porosity The Cell Inspector window should look similar to the image below Exploring the Model Grid and Properties 5 a Cell Values Options Grid Position Q Model Position Q world Position r Specific Yield Specific Storage amp Effective Porosity Initial Head Initial Concentration Longitudinal dispersivity Boundary Conditions Q Output AllOn All Off e Cell Values tab While leaving the Cell Inspector window opened move your mouse across the cells in the grid You will see the assigned property values populate in the Cell Inspector window Also you will notice that all properties are uni
13. ersion 4 A Computer Program for Simulation of Multi Species Solute and Heat Transport U S Geological Survey Techniques and Methods Book 6 Chapter A22 39 p Terms and Notations For the purposes of this tutorial the following terms and notations will be used Type type in the given word or value lt gt press the Tab key on your keyboard lt Enter gt press the Enter key on your keyboard click the left mouse button where indicated eS double click the left mouse button where indicated The bold faced type indicates menu or window items to click on or values to type in denotes a button to click on either in a window or in the side or bottom menu bars Opening the Visual MODFLOW Model For your convenience part of the model has already been prepared for you The Visual MODFLOW project file is located on your computer and can be opened by following the directions below A full version of the model following this step by step tutorial has already been created and is located on the FTP site To download the full version of the model please visit the following ftp address ftp ftp flowpath com Software VMOD 2010 1 Tutorials 2 SEAWAT 4 Tutorial Example Problem NOTE Some features described in this tutorial are only available in a Standard SEAWAT Pro or Premium version If you are using a standard version of VMOD and wish to upgrade your software please contact sales waterloohydrogeologic com
14. es for various model parameters In the corresponding fields modify the following parameters according to the values in the table below Setting Up the SEAWAT Flow Engine 9 Parameter Value RHOB Bulk Density kg m 3 1760 DENSEMIN Minimum Fluid Density oc DESNMAX Maximum Fluid Density Ey DENSESLP Density Concentration Slope VISCMIN Minimum Fluid Viscosity oO VISCMAX Maximum Fluid Viscosity Table 5 Constant Variant Model Parameters The default values for the remaining parameters are acceptable Now that all the parameters have been defined we can create the variant a OK button e OK button at the Edit Engines window Editing Dispersion Properties we Input from the main menu s Properties gt Dispersion s Layer Options from the side menu The Dispersion Package window will appear on your screen In the empty field above the Vert Long Dispersivity column type 0 1 10 SEAWAT 4 Tutorial Example Problem Dispersion Packar Edit columns er OK button e Database from the side menu The Dispersion window will appear on your screen In the DI m field type 1 BED DI m Distribution Array Longitudinal dispersivity Value 1 EERIE Reset Order ean Up Standard lt lt OF Cancel er OK button Setting Up the SEAWAT Flow Engine 11 Defining Model Boundary Conditions Constant Concentration Along the right side of the grid a constant concentration boundary is used to h
15. fined constant 5 OK Cancel Multiple assigning layers 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 2 er OK button Setting Transport Time Steps From the main menu e SEAWAT gt Transport gt Output Time Steps The Output and Time Step Control dialog box will appear on your screen In the Simulation Time field type 200000 default In the Max of Transport Steps type 1000000 Under Save Simulation Results at s Specified Time day s Under Output Time s day type 5000 lt Enter gt type 10000 lt Enter gt type 30000 lt Enter gt type 60000 lt Enter gt type 200000 The Output and Time Step Control dialog box should appear similar to the one shown below Defining Simulation Run Settings 21 Simulation time 200000 day s Output Timels Max transport Steps 1000000 a s Specify max stepsize fo day s Save simulation results at C End of simulation only Every Nth time step N fi 0 Specified times day s IV Save CBM CCM file for Mass Balance output OK Running SEAWAT To run the model simulation G Run from the toolbar er In the check box beside SEAWAT Because the Variable Density Flow VDF package is only compatible with the Layer Property Flow Package you have to set the model to run with the VDF package s Advanced gt gt a User Defined under the Settings column for SEAWAT e button
16. form across the entire model grid The values for each property are shown in the following table Property Specific Yield Specific Storage le 05 1 m Effective Porosity 0 35 Table 2 Assigned property values To close the Cell Inspector s X top right corner of window s File gt Main Menu from the Input main menu The next section will guide you through the steps of setting up a SEAWAT flow engine including the creation of a transport variant that contains both salt and temperature species SEAWAT 4 Tutorial Example Problem Setting Up the SEAWAT Flow Engine To configure the SEAWAT flow engine er Setup gt Engines from the Visual MODFLOW main menu The Edit Engines dialog box will appear on your screen a USGS SEAWAT from the flow engine drop down list box e Yes at the Confirmation message a OK at the Information message x M Time Options Start Date 21 01 2009 Start Time 12 00 00 AM Run Type Transient Flow Steady State Simulation Time 1 lda Flow Options Flow Engine ERIE Rg evn Simulation Type Groundwater flow Compatible Transport Variant Transport Engine NONE Total Number of Species Sorption P Number of Mobile Species oo Reactions oOo O O Reaction Parameters Po Description Delete Copy New Edit OK Cancel Help In order to support all of the available options for the multi species reactive transport prog
17. lt the salt concentration contours are shown in red To make the contours more visible turn off the other layers shown on the display e Overlay from the bottom menu er Uncheck C O Temperature e Uncheck C O Head Equipotentials a OK To further visually enhance the concentration output you can enable color shading E Options button directly above the Time button Concentration Options The following dialog will load 24 SEAWAT 4 Tutorial Example Problem Contours Color Shading JV Automatic reset min max and interval values Contour Lines V Automatic contour levels Minimum 5000 mg L Custom contour levels Maximum 40000 mg L Interval 5000 mg L Color m Maroon gt Width 1 Contour Label IV Visible 4 2 Decimals 1 Ti Size E Add Delete Insert T fi Width al Contouring resolution Speed Highest Slowest s Color Shading tab The following dialog will load C O Salt Contours Color Shading I7 Use Color Shading Range of Values Cut off Levels Minimum 5000 mg L I Lower Level 0 mg L Maximum 40000 mg L Upper Level D mg L Ranges to colo Transparency 30 B z Min I 5000 058 823 058 250 Reset Apply OK Cancel or Use Color Shading check box Leave the rest of the settings as default Viewing Output er OK to apply the changes and close the dialog s Nex
18. old the temperature constant at 5 degrees allowing for heat conduction to occur to or from the seawater boundary To create this boundary condition er Boundaries gt Constant Concentration e Yes to save changes To make it easier to assign this boundary condition you may want to zoom into the right boundary portion of the model grid You can do this by using the Zoom In button located at the bottom toolbar The boundary condition will be assigned to all cells in column 99 Once this column is visible on your screen proceed with the steps below e Assign gt Line from the side menu e Inside cell Row 1 Column 99 Layer 1 er right click Inside cell Row 1 Column 99 Layer 100 After right clicking all the cells in column 99 will be colored blue and the following dialog will appear on your screen one H f H a New Multiply values by ho mz start Time Stop Time ssi tmgvL Temperature a i 7 oon NA Eopy schedule tron In the Stop Time field type 200000 In the Temperature field type 5 OK button 12 SEAWAT 4 Tutorial Example Problem Point Source The next step is to assign point source concentrations at both the left and right boundaries However because point source concentrations must be associated with a flow boundary you must assign the appropriate flow boundaries first On the left boundary the well package is used to define an inflow of water from the boundary For your
19. r Density Scenario s Calculate density from multiple species option 16 SEAWAT 4 Tutorial Example Problem Selecting this option will ensure that both temperature and salinity are included when calculating density values Under Internodal density calculation algorithm MFNADVFD s Upstream weighted algorithm option Under Initial Time Step type 0 01 e VSC Settings tab VSC Settings The VSC Settings tab contains various options for the Viscosity VSC package x4 VDF Settings VSC Settings Flow Settings Transport Settings Viscosity Options MT3DMUFLG Calcualte viscosity from single salt only species no temperature dependence Calcualte viscosity from temperature and multiple species Use specified viscosity distribution Viscosity Temperature Dependence Options MUTEMPOPT Simple linear relation Voss 1984 parametrization C Pawlowski 1991 parametrization Guo and Zhao 2005 parametrization not applicable for water Specified Viscosity Options INVISC f Values for the viscosity will be set ta YVISCREF for 1 st stress period For the next period viscosity will be reused from the previous Bead specified viscosity from UCN file Read concentration from UCN file and convert to viscosity using linear expression Viscosity UEN file name EA Concentration UCN file name lt Back Next gt OK Cancel Help Under Viscosity Option MT3DMUFLG e Calculate viscosit
20. rams Visual MODFLOW requires you to setup the initial conditions for the transport scenario e g number of species names of each species initial concentrations decay rates partitioning coefficients etc Each scenario is referred to as a Transport Variant and you can have more than one variant for a given flow model To create a new Transport Variant e New button The Variant Parameters dialog box will appear on your screen Enter the following information into the available fields For the Variant Title field type VAR002 For the Description field type Simultaneous heat and salt transport example Setting Up the SEAWAT Flow Engine 7 For the Sorption drop down list box s Linear Isotherm equilibrium controlled For the Reactions drop down list P No kinetic reactions leave as is The top half of the Variant Parameters dialog box show appear similar to the image below x Variant Title vAR002 Mar002 Transport Engine USGS SEAWAT Description Example el Sorption Linear isotherm equilibrium controlled x Total Number of species fi Reactions No kinetic reactions Number of mobile species fi Are the reaction parameters constant or spatially variable Constant Next you will specify the appropriate parameters for the variant parameters In the bottom half of the dialog you will see that by default the Salt species already exists However in order to simulate heat transport simultaneousl
21. s ITER1 type 25 For Head change criterion HCLOSE type le 5 For Residual criterion RCLOSE type 1 For Printout interval IPRPCG type 5 along the left side Finally under package options top right corner of dialog User Defined Specified Head Package Run 4 4 4 Transport Settings tab node to unselect the checkbox SEAWAT 4 Tutorial Example Problem Transport Settings Simulation Scenario VDF Settings VSC Settings Flow Settings Transport Settings Solution Options Solve Implicitly Use GCG solver Advection GG Options C Solve Explicitly Package options Setting C Automatic UserDefined Reset Chemical Reactions Package GCG Solver Package M Time Steps and GCG options Viscosity Package TvD x Particle Options Basic Package Run Trans Advection Package Run Trans M Dispersion Reaction and Sink Sources Dispersion Package Run Trans Solve Implicitl Use GCG solver Sink Sources Package Run Trans Bun Trans JRun Trans Run Trans Initial step size DTO 0 day Maximum Step Size 1 day Multiplier 1 M Porosity Options Use Effective Porosity Use Total Porosity Min sat thickness fo 01 fraction of cell thick Courant number 0 75 Dispersion Reaction and Sources Sink Under Solution
22. t in the side menu beside Time to advance the contours to the next time period 26 oe Concentration eS ees Seet H r Watch as the contours trace the advancement of the salt concentration front from left to right For each time step the salt concentration contours and color shading should look similar to the ones shown below 5000 Days nut 0 Eier eoon Sietle iraa Ei PULL L RUU OE BEM ES Eii roo U Eii SEAWAT 4 Tutorial Example Problem Viewing Output 1200 1200 30 000 Days LOR OOO DENE 200 000 Days 21 The contours should appear similar to the ones published in the USGS SEAWAT 4 User s Manual shown below time 60 000 days time 10 000 days time 30 000 days time 5 000 days EXPLANATION 5 percent seawater line for Case 7 50 percent seawater line for Case 6 You can do the same for the Temperature species e Species gt Temperature from the side menu Before turning on the temperature contours disable the salt concentration contour overlay e Overlay from the bottom menu s Uncheck C O Salt s Check C O Temperature E OK You can enable color shading for the temperature species by following the steps described previously s Next in the side menu beside Time to advance the contours to the next time period Watch as the contours trace the advancement of the temper
23. y from temperature and multiple species When this option is selected you can select from various parametrization options for calculating the temperature term in the viscosity equation In this case we will retain the default setting Voss 1984 Parametrization eo Flow Settings tab Defining Simulation Run Settings 17 Flow Settings 18 Simulation Scenario i x VDF Settings VSC Settings Flow Settings Transport Settings Solver m Package options iterati 100 Max outer iterations MITER Pei a E T Aon PCG Max inner iterations ITER1 25 l Setting Value Head change criterion HCLOSE 1e5 Basic Package Run Translate ee Layer Property Package Block Cente Residual criterion RCLOSE 1 Well Package Run Translate D iad pooo Oc Package Run Translate SIP S E Specified Head Package Translate Printout interval IPRPCG CUNIT 20 LUNIT 30 M Pre conditioning method NPCOND Extension CHD Cholesky NPCOND 1 TARU E ES C SOR Translate Iv C Polynomial NPCOND 2 Discretization Package Run Translate FTL Package Translate Zz Calb Package Run Translate Relaxation parameter RELAX L Nodry Package Run Translate C WHS J Calculate estimate NPBGL 1 C AMG C GMG Printout interval For the Solver s PCG option For Max outer iterations MXITER type 100 For Max inner iteration
24. y with solute salt transport you need to add a Temperature species To add a Temperature species er Temperature button Species Model Params Species Params Desgnaion Component Deseipton Mabie SCONCTng I sPTIVngAl OAHODE oner 0 0 0 Salt Density controlling species Temperature Degrees Celsius 0 0 0 New Species Delete Species Temperature OK Cancel Help For both the Salt and Temperature species you can now enter the appropriate parameter values For the Salt species row 1 type the following parameter values into the appropriate fields under the Species or Species Params tab Note The Species Params tab contains descriptions of each parameter SEAWAT 4 Tutorial Example Problem Parameter SCONC mg L SP1 1 mg L DRHODC CREF mg L MDCOEF m7 Day DMUDC m7 Day Table 3 Variant parameter values for salt species Once these values have been entered under the Species tab set the Temperature species to mobile s Yes under the Mobile column for the Temperature species For the Temperature species type the following parameter values into the appropriate fields under the Species or Species Params tab Parameter SCONC mg L SP1 1 mg L DRHODC CREF mg L MDCOEF m7 Day DMUDC m Day Table 4 Variant parameter values for temperature The next step is to enter the appropriate model parameters eo Model Param tab Visual MODFLOW provides typical valu
Download Pdf Manuals
Related Search