ChemFlux Tutorial Manual
Contents
1. 9 Select Model gt FEM Options from the menu to open the FEM Options dialog 10 Click on the Advanced button 11 Click on the Mesh Generation Controls tab 12 Set the NODELIMIT to 10000 13 Press OK to close the FEM Options dialog The next step indefining the model is to enter the Material Properties for the single material that will be used in the model Only one material is used for the model with these properties Longitudinal Dispersivity 0 Im Transverse Dispersivity o Im Diffusion Coefficient D Om day 1 Open the Materials Manager dialog by selecting Model gt Materials gt Manager from the menu 2 Click the New Material button to create a material type in a name for the material as 3D Tutorial Soil and click OK The Material Properties dialog will open automatically 3 Move to the Dispersion tab 4 Refer to the data provided above Enter the Longitudinal Dispersivity a 1m 5 Enter the Transverse Dispersivity o Im 6 The Diffusion option is set to Constant as the gradient file specified does not contain volumetric water content which is required to define a diffusion curve 7 Enter the Diffusion Coefficient D Om day 8 Close the Material Properties and Material Manager dialogs SOILVISION systems LTD A Three DimensionalExample Model 30 of 37 Each region will cut through all the layers in a model creating a separate block on each layer Each block can be assigned a mate2. Boundary Condition Zero Flux Boundary Name 3 Close the dialog The newly specified boundary condition will be displayed with symbols on the CAD window d Apply Material Properties The next step in defining the model is to specify the settings that will be used for the model The Settings dialog will contain information about the current model System Units Time and contaminant transport processes 1 To open the Settings dialog select Model gt Settings in the menu 2 Puta check mark in the Advection and Dispersion boxes in the Processes section under the General tab if they are not already there Move to the Time tab Enter a Start Time of 0 an Initial Increment of 50 days and an End Time of 400 days Select the Advection tab N ME oS Choose Import from the Advection Control option SOILVISION SYSTEMS LTD A Three Dimensional Example Model 29 of 37 7 Click Browse 8 Specify the file ChemfluxInput_Reservoir3D_l trn that was generated by SVFLUX This file can be found in the following directory C SVS ModelFilesSVSlope Tutorial 3d SteadyState Reservoir3D NOTE It is very important that the TRN file and the geometry are obtained from the same SVFLUX model In order to improve solution time for the purposes of this tutorial certain finite element options will be set The finite element mesh node limit and grid spacing will be set to generate a simpler mesh that will reduce the solution time
3. CHEMFLUX 2D 3D Contaminant Transport Modeling Software Tutorial Manual Written by Robert Thode B Sc G E Edited by Murray Fredlund Ph D SoilVision Systems Ltd Saskatoon Saskatchewan Canada Software License The software described in this manual is furnished under a license agreement The software may be used or copied only in accordance with the terms of the agreement Software Support Support for the software is furnished under the terms of a support agreement Copyright Information contained within this User s Manual is copyrighted and all rights are reserved by Soil Vision Systems Ltd The CHEMFLUX software is a proprietary product and trade secret of Soil Vision Systems The User s Manual may be reproduced or copied in whole or in part by the software licensee for use with running the software The User s Manual may not be reproduced or copied in any form or by any means for the purpose of selling the copies Disclaimer of Warranty SoilVision Systems Ltd reserves the right to make periodic modifications of this product without obligation to notify any person of such revision SoilVision does not guarantee warrant or make any representation regarding the use of or the results of the programs in terms of correctness accuracy reliability currentness or otherwise the user is expected to make the final evaluation in the context of his her own problems Trademarks Windows is a registered trade
4. Below is a description of the seepage model solved by SVFLUX SOILVISION SYSTEMS UTD A Three Dimensional Example Model 23 of Project Ponds Model Reservoir3D Minimum authorization required FULL e Model Dimensions 2 1 R 24 The data points for the surface grids can be found Appendix B Enter these points to set up the SVFLUX model geometry Boundary Conditions 37 SOILVISION systems LTD A Three DimensionalExample Model 24 of 37 Head 10 5m ak The steady state seepage model is set up to simulate a pond or reservoir a certain distance from ariver channel The water levels in the reservoir and river channel are set using head boundary conditions The level of water in the reservoir is set using a Head Expression 10 5m set on surface 2 for the reservoir region The level of water in the river channel is set using a Head Expression 7m set on the line segment extending from point 14 0 to 14 27 on surface 1 Material Properties There is only one material in the saturated 3D example model Two regions have been implemented in this model in order to apply the necessary boundary conditions The material in the model has a hydraulic conductivity ksat 2 17e 01 m d e FlowRegime SOILVISION systems LTD A Three Dimensional Example Model 25 of 37 Flow lines show that groundwater is flowing from the reservoir toward the adjacent river channel The presence of unsaturated material n
5. Gradients VRY t from 0 by 5 to 20 Solution Mesh att 20 Mesh att 0 Add New Plot l Delete Properties Multiple Update 2 l Plot Settings Custom Plots OUTPUT MANAGER There are many output file types that can be specified to export the results of the model One will be SOILVISION SYSTEMS LTD Sudicky Model 18 of 37 generated for this tutorial example model afile to transfer the results to ACUMESH 1 Open the Output Manager dialog by selecting Model gt Reporting gt Output Manager from the menu 2 The toolbar at the bottom left corner of the dialog contains a button for each output file type Click on the ACUMESH button to begin adding the output file The Output File Properties dialog will open 3 Enter the title ACUMESH 4 Click OK to close the dialog and add the output file to the list Specify Time Steps The Output File Properties dialog also allows the user todefine timesteps for the current model This can be accessed by using the Update Method tab on the Output File Properties dialog 1 Enter a Start Time of 0 a Time Increment of 1 yr and an End Time of 20 yr 2 Click OK to close the Output File Properties dialog and return to the workspace 3 Click OK to close the Output Manager and return to the workspace f Run Model The current model may be run by selecting the Solve gt Analyze menu option g Visualize Results The flow vectors for the current model
6. The World Coordinate Settings and the View Settings may need to be set up again 1 Access the World Coordinate System dialog by selecting View gt World Coordinate System from the menu 2 Enter the World Coordinates System coordinates shown below into the dialog x minimum 10 y minimum 2 x maximum 260 y maximum 8 3 Click OK toclose the dialog For betterviewing results set the magnification factor in the View Settings SOILVISION SYSTEMS LTD Sudicky Model 15 of 37 1 Select View gt Settings from the menu 2 Set the aspect ratio near 1 10 This will magnify the model s height 3 Click OK toclose the dialog c Specify Boundary Conditions In general flow models must have a defined entry and exit point for water to flow The boundary conditions shown at the start of this model may be entered through the following steps 6 Boundaries Region Name Region 1 Select Shape Index 585083471 Y Boundary Condition Expression or Data Units Tm 250 0 Continue 250 5 375 Concentration Expression 0 afm 3 175 55 Continue 125 6 333 Continue 80 6 393 Concentration Expression ift lt 5 then 1 else 0 g m 3 40 6 447 Concentration Expression 0 g m 3 0 6 5 Flux Expression 0 g vr m Update Selected Segment Boundary Condition Flux Expression v E Lenath Expression 0 l NOTE boundary conditions defined at a point remain in effect until re defined at a subsequent point 1 Select Region 1 Region 1 m
7. Y Z 0 0 11 14 0 11 0 10 11 14 10 11 0 11 10 14 11 11 0 16 10 14 16 11 0 17 11 14 17 11 0 27 11 14 27 11 2 0 11 21 0 4 2 10 11 21 10 4 2 11 10 21 11 4 2 16 10 21 16 4 2 17 11 21 17 4 2 27 11 21 27 4 3 0 11 24 0 4 3 10 11 24 10 4 3 11 11 24 11 4 3 16 11 24 16 4 3 17 11 24 17 4 3 27 11 24 27 4 37 of 37 SOILVISION SYSTEMS LTD References 38 of 37 4 References FlexPDE 5 x Reference Manual 2007 PDE Solutions Inc Spokane Valley WA 99206 Fredlund D G and Xing A 1994 Equations for the soil water characteristic curve Canadian Geotechnical Journal Vol 31 No 3 pp 521 532 Freeze R Allan and Cherry John A 1979 Groundwater Prentice Hall Inc Englewood Cliffs New Jersey Sudicky E A 1989 The Laplace transform Galerkin technique A time continuous finite element theory and application to mass transport in groundwater Water Resources Research Volume 25 Issue 8 p 1833 1846 Zheng C and Wang P 1999 MT3DMS Documentation and User s Guide Report to the US Army Corps of Engineers Waterways Experiment Station available at http hydro geo ua edu
8. gradient vectors SOILVISION SYSTEMS LTD Sudicky Model 11 of 37 1 Open the Plot Manager dialog by selecting Model gt Reporting gt Plot Manager from the menu Plot Manager Plots Point Piezometer Area Volume Flux Sections Boundary Flux Climate Review Boundary Other Title Variable Restriction PLOT MO KO Pressure u Yes zhead h Yes zpressure u Yes Head h vector gradx grady Mesh Add New Plot Copy Delete Properties l Multiple Update 2 The toolbar at the bottom left of the dialog contains a button for each plot type Click on the Contour button tobegin adding the first contour plot The Plot Properties dialog will open Enter the title Pressure Select u as the variable to plot from the drop down Click the Output Options tab and ensure that only Plot is checked off Click OK to close the dialog and add the plot to the list sl 2 ta A Repeat steps 2 to 6 to create the plots shown in the above dialog The zoomed plots are not necessary they are used to closely examine key zones in the problem 8 Click OK to close the Plot Manager and return to the workspace Alternatively the user may press the Add Default Plots button and typical plots will be added to the plot list OUTPUTFILES There are four output file types that can be specified to export the results of the model One will be generated for this tutorial example model a plot to transfer the results to AC
9. may be visualized through the following steps 1 Open ACUMESH View gt ACUMESH menu option 222 Results and Discussion After the model has finished solving the results will be displayed in the dialog of thumbnail plots within the CHEMFLUX solver Right click the mouse and select Maximize to enlarge any of the thumbnail plots The following is a short summary of plots illustrating the movement of the plume through the model for times of 8 years 12 years and 20 years SOILVISION SYSTEMS LTD Sudicky Model 19 of e Time 8 years 7 00 The source has been shut off for 3 years e Time 12 years 7 00 6 00 5 00 4 00 3 00 2 00 1 00 0 00 SOILVISION SYSTEMS LTD Sudicky Model 20 of e Time 20 years 2 3 Appendix A Region 1 Region 2 Region 3 X Y X Y X Y 0 0 0 2 180 2 250 0 120 2 250 2 250 5 375 120 4 250 4 175 5 5 0 4 180 4 125 6 333 80 6 393 40 6 447 0 6 5 SOILVISION SYSTEMS UTD Sudicky Model Boundary Conditions X Y Boundary Condition 0 0 Gradient Expresion 0 250 0 Continue 250 5 375 Concentration Expression 0 175 5 5 Continue 125 6 333 Continue 80 6 393 Concentration Expression if t lt 5 then 1 else 0 40 6 447 Concentration Expression 0 0 6 5 Gradient Expresion 0 21 of 37 SOILVISION SYSTEMS LTD A Three Dimension
10. Fit button 6 Enter the Hydraulic Conductivity data Choose the Hydraulic Conductivity Tab and enter the Ksat and Ky ratio values The dialog can be closed once material properties are entered The ChemFlux Soil2 material s properties are available in Appendix A and can be entered in the same manner as ChemFlux Soil 7 Apply toregions The material properties can be applied to regions by opening the Regions dialog Model gt Geometry gt Regions and selecting the appropriate materials from the drop down boxes ChemFlux Soill should be applied to Region 1 and ChemFlux Soil2 should be applied to Regions 2 and 3 e Specify Model Output Two levels of output may be specified i output graphs contour plots fluxes etc which are displayed during model solution and ii output which is written to a standard finite element file for viewing with ACUMESH software Output is specified in the following twodialogs in the software i Plot Manager Output displayed during model solution ii Output Manager Standard finite element files written out for visualization in ACUMESH or for initial condition input to other finite element packages PLOT MANAGER The Plot Manager dialog is first opened to display appropriate solver graphs There are many plot types that can be specified to visualize the results of the model A few will be generated for this tutorial example model including a plot of the solution mesh pressure contours head contours and
11. Open the SVOFFICE Manager dialog 2 Select the projectcalled UserTutorial 3 Press the New button under the Models heading and enter User CHEMFLUX3D as the model title SOILVISION SYSTEMS LTD A Three Dimensional Example Model 27 of 37 4 Select the following Application CHEMFLUX System 3D Vertical Type Transient Units Metric Time Units Days 5 Click the OK button to save the model and close the New Model dialog 6 The new model will automatically be opened in the workspace b Enter Geometry The geometry for the model must be imported from SVFLUX before any other modeling can be done in CHEMFLUX 1 Select the Model gt Geometry gt Import Geometry gt From Existing Model menu 2 The Import Geometries menu will pop up Select the appropriate project name Tutorials 3 Select the Reservoir3D model 4 Press the Import button 5 A pop up message will appear stating current surfaces geometry features art objects flux sections and plots referencing a specific region to be deleted Do you wish to continue Click on Yes 6 A pop up message will appear asking if you want to copy material properties and assignments Click on No The import includes any regions region shapes surfaces surface grids and elevations These parts of the model definition are fixed in CHEMFLUX World Coordinate System settings and features are also imported if present but may be edited in CHEMFLUX c Specify Boundary Co
12. UMESH 1 Open the Output Manager dialog by selecting Model gt Reporting gt Output Manager from the menu 2 The toolbar at the bottom left corner of the dialog contains a button for each output file type Click on the ACUMESH button to add the output file with the default variables Press the CHEMFLUX button to add the output file with the default variables Click OK to close the dialog and add the output file to the list Click the Settings button on the Output Manager dialog to open the Output Settings dialog N ow Pp amp Ensure the Region Separation checkbox is checked press OK to close the dialog SOILVISION SYSTEMS LTD Sudicky Model 12 of 37 7 Click OK to close the Output Manager and return to the workspace f Run Model The current model can be run by selecting the Solve gt Analyze menu option g Visualize Results The flow vectors for the current model can be visualized through the following steps 1 Open ACUMESH View gt ACUMESH menu option 2 Plot Flow Lines Plot gt Flow Lines 2 1 2 Results and Discussion 7 6 E 5 c 4 Y 3 D i 2 i 100 200 Distance Along Flow Direction m 2 2 CHEMFLUX Model Now that the steady state flow hydraulic head gradients have been established in the SVFLUX software the focus turns to solving for the chemical concentrations with time for the solution domain In order to solve this model the user needs to perform the following steps 1 Create a
13. al Example Model 22 of 37 3 A Three Dimensional Example Model The following example will introduce you to the three dimensional model in CHEMFLUX The model will be used to investigate if contaminant from areservoir will travel to ariver channel due to advection and dispersion processes within a 400 day time period The 400 day time period was chosen as the time necessary to install a pumping well between the river channel and the reservoir The well will be used to pump contaminant from the ground to ensure the plume will notreach the river channel The example model begins with a brief description of the steady state seepage analysis completed to provide CHEMFLUX with computed seepage gradients Next adetailed set of instructions guides the user through the creation of the 3D contaminant transport model Project Ponds Model ResevoirChemFlux Minimum authorization required FULL Model Description and Geometry It is important to note that you will be analyzing the SVFLUX model before the CHEMFLUX model is completed 3 1 Steady State SVFLUX Solution Advection is known as the process by which solutes are transported by the bulk motion of the flowing groundwater Freeze and Cherry 1979 The bulk motion of the flowing groundwater or seepage gradients are solved using SVFLUX SVFLUX calculates the seepage gradients and writes them to a text file The CHEMFLUX solver then reads this text file when calculating the contaminant transport solution
14. al forRegion 1 e repeat for Region 2 andRegion 3 click OK to close the dialog Next the settings that will be used for the model must be specified To open the Settings dialog select Model gt Settings in the workspace menu The Settings dialog will contain information about the current model System Units Time and contaminant transport processes 1 To open the Settings dialog select Model gt Settings in the workspace menu 2 Check Advection and Dispersion in the Processes box under the General tab 3 Choose the Time tab Enter a Start Time of 0 a Time Increment of 1 yr and an End Time of 20 YT 4 Select the Advection tab 5 Choose Import from the Advection Control option 6 Click Browse 7 Specify the gradient file Examples_ChemFlux2D trn that was generated by SVFLUX in the previous example 8 Press OK to close the Settings dialog NOTE It is very important that the TRN file and the geometry are obtained from the same SVFLUX model e Specify Model Output Two levels of output may be specified i output graphs contour plots fluxes etc which are displayed during model solution and ii output which is written to a standard finite element file for viewing with ACUMESH software Output is specified in the following twodialogs in the software SOILVISION SYSTEMS LTD Sudicky Model 17 of 37 i Plot Manager Output displayed during model solution ii Output Manager Standard finite element file
15. at flow is from right to left towards the river in this view withhigher gradients near the reservoir The following is a short summary of plots created in ACUMESH illustrating the movement of the plume through the model for times of 50 days 100 days and 400 days Note that the plume does not reach the river channel in within the 400 day time period The below diagram was created in ACUMESH by plotting concentration contours and varying time 1 Open ACUMESH byselecting Window gt ACUMESH from the menu 2 Select Plots gt Contours from the menu 3 Select c from the Variable Name drop down 4 Click OK toclose the Contours dialog 5 Select the desired timestep from the Time drop down on the toolbar SOILVISION systems LTD A Three Dimensional Example Model e Time 50 days e Time 100 days of 37 SOILVISION SYSTEMS LTD A Three DimensionalExample Model e Time 400 days 36 of 37 SOILVISION SYSTEMS LTD A Three Dimensional Example Model 3 4 Appendix B e Surface 1 Grid X Y Z X Y Z 0 0 0 14 0 0 0 10 0 14 10 0 0 11 0 14 11 0 0 16 0 14 16 0 0 17 0 14 17 0 0 27 0 14 27 0 2 0 0 21 0 0 2 10 0 21 10 0 2 11 0 21 11 0 2 16 0 21 16 0 2 17 0 21 17 0 2 27 0 21 27 0 3 0 0 24 0 0 3 10 0 24 10 0 3 11 0 24 11 0 3 16 0 24 16 0 3 17 0 24 17 0 3 27 0 24 2T 0 e Surface 2 Grid X Y Z X
16. button next to the list of projects Fe a Create a new model called User_Vanderheijde by pressing the New button next to the list of models The new model will be automatically added under the recently created UserTutorial project Use the settings below whencreating a new model 5 Select the following Application CHEMFLUX System 2D Vertical Type Steady State Units Metric Time Units Seconds s 6 Clickon OK Before entering any model geometry it is best to set the World Coordinate System to ensure that the model will fit into the drawing space 1 Access the World Coordinate System tab by selecting the World Coordinate System tab on the Create New Model dialog 2 Enter the World Coordinates System coordinates shown below into the dialog X minimum 10 y minimum 2 x maximum 260 y maximum 8 3 Click OK toclose the dialog For betterviewing results set the magnification factor in the View Settings 1 Select View gt Settings from the menu 2 Set the aspect ratio near 1 10 This will magnify the model s height 3 Click OK toclose the dialog b Enter Geometry The geometry must be defined for the SVFLUX model SOILVISION SYSTEMS LTD Sudicky Model 9 of 37 Add Region The model is created with one default region Another region can be added under the Model gt Geometry gt Regions dialog by pressing the New button andclosing the dialog SelectRegion The user must select the region they would like to dra
17. e user must decide the project under which the CHEMFLUX SOILVISION SYSTEMS LTD Sudicky Model 14 of 37 model is going to be organized If the project is not yet included in the Projects section of the SVOFFICE Manager you must add the project before proceeding withcreating the model In this case the model is placed under the projectcalled UserTutorial To add a model 1 Open the SVOFFICE Manager dialog Select the project called UserTutorial Press the New button under the Models heading Select CHEMFLUX for the Application Enter User_Example2D in the Model Name box Select 2D for System Transient for Type Metric for Units and Years for Time Units Click the OK button to save the model and close the New Model dialog CO ae ON a es le The new model will automatically added be added to the Models list NOTE You will notice that there is no distinction between steady state and transient state in CHEMFLUX This is because all CHEMFLUX models are considered to be transient state b Enter Geometry The geometry for the model can be obtained in the spreadsheet located here Entering the geometry into the newly created SVFLUX model can be accomplished through the following steps 1 Select the Model gt Geometry gt Import Geometry gt From Existing Model menu 2 The Import Geometries dialog will pop up Select the appropriate project name Tutorial 3 Select the Vanderhiejde model 4 Press the Import button
18. ear the surface of the model is causing water to first flow down to the saturated zone and then move toward the river channel 3 2 CHEMFLUX Model Setup Once the gradients have been calculated in the SVFLUX software the focus may be directed towards the calculation of contaminant movement in the CHEMFLUX software This part of the tutorial involves setting up the CHEMFLUX model which will use the gradients calculated in SVFLUX as well as the diffusion process to determine the location of the resulting contaminant contours Project Ponds Model ResevoirChemFlux Minimum authorization required FULL SOILVISION systems LTD A Three Dimensional Example Model 26 of 37 30 30 e CHEMFLUX Material Properties Please note the SVFLUX Solution shown in the diagrams above are a result of the SVFLUX Reservoir3D Tutorial In order to set up the CHEMFLUX model described for this tutorial the following steps will be required The steps forcreating a model fall under the general categories of a Create model b Enter geometry c Specify boundary conditions d Apply material properties e Specify model output f Run model g Visualize results a Create Model The first step indefining a model is to decide the project under which the model is going to be organized If the project is not yet included you must add the project before proceeding with the model In this case the model is placed under the projectcalled Tutorial To add amodel 1
19. mark of Microsoft Corporation Soil Vision is a registered trademark of Soil Vision Systems Ltd SVOFFICE is a trademark of SoilVision Systems Ltd SVFLUX TM is a trademark of Soil Vision Systems Ltd CHEMFLUX is a trademark of Soil Vision Systems Ltd SVAIRFLOW is a trademark of SoilVision Systems Ltd SVHEAT is a trademark of SoilVision Systems Ltd SVSOLID TM is a trademark of Soil Vision Systems Ltd SVSLOPE TM is a trademark of Soil Vision Systems Ltd ACUMESH TM is a trademark of Soil Vision Systems Ltd FlexPDE is a registered trademark of PDE Solutions Inc Copyright 2008 by Soil Vision Systems Ltd Saskatoon Saskatchewan Canada ALLRIGHTS RESERVED Printed in Canada SOILVISION SYSTEMS LTD Table of Contents Pee T Nirod UCT ON EE ER ER Hee eel eh ei hie ist Naver aad Pe 2 Sudicky Mleto Ce PEE ee ee ee ee EE ee ee ee Ee ee ee ee ee ee ee ee ee ee ee ee AN 2 1 Steady State SVFLUX Model iese ee ek ee Ge Re ee ee ke ee 2 1 1 Model Setup Ge EER EWER Ge 2 2 CHEMPLUX Model isusswii Ese ses Ee EVER GER ese eek esse Gee ee Geek ee ee ew ke 2 2 1 Model Setup sd d Appendi Ai eo cet tenet lords 3 A Three Dimensional Example Model esse ed ee ee ee ded Steady State SVFLUX Solution ee ee ee ke ee Ge ee ee 3 2 CHEMFLUX Model SetuP sesse ese ee eek ees ee ee ee ee ee be ee 343 Results and DisCuSSION ee ee ee ee ee AR ee GR Ke ee AR be ME ie AE OE ME EE OO T SOILVISION SYSTEMS LTD Intr
20. n transient models where water storage will change with time The initial positive slope on the SWCC applies for the low suction range SOILVISION SYSTEMS LTD Sudicky Model 7 of 37 e SVFLUX Flow Regime 7 6 E 5 c 4 3 D i 2 i 100 200 Distance Along Flow Direction m 2 1 1 Model Setup In order to set up the model described in the preceding section the following steps or categories will be required a Create model b Enter geometry c Specify boundary conditions d Apply material properties e Specify model output f Run model g Visualize results a Create Model Since FULL authorization is required for this tutorial the user must perform the following steps to ensure full authorization is activated 1 Plug in the USB security key 2 Go to the File gt Authorization dialog on the SVOFFICE Manager 3 Software should display full authorization of Standard or Professional If not it means that the security codes provided by SoilVision Systems at the time of purchase have not yet been entered Please see the the Authorization section of the SVOFFICE User s Manual for SOILVISION SYSTEMS LTD Sudicky Model 8 of 37 instructions on entering these codes The following steps are required to create the model 1 Open the SVOFFICE Manager dialog Select ALL under the Applications combo box and ALL for the Model Origin combo box Create a new projectcalled UserTutorial by pressing the New
21. nditions In general flow models must have a defined entry and exit point for water to flow The boundary conditions shown at the start of this model may be entered through the following steps 1 Select Region 1 Slope and Surface 1 must be selected by clicking on the Region amp Surface dialogues at the top of the screen Slope m Fa GP surface 1 2 Enter Boundary Conditions 1 The Boundary Conditions dialog may be displayed under the Model gt Boundaries gt Boundary Conditions menu option Once in the dialog the user needs to Select the starting node point 14 0 SOILVISION SYSTEMS LTD A Three Dimensional Example Model 28 of 37 Then select a Concentration expression from the combo box e Enter a value of 0 1 g m3 Select the node point 14 27 and specify a zero flux boundary condition Please note although the above boundary conditions may appear to be entered already by default from the geometry import you performed you still need to enter the above conditions for CHEMFLUX analysis to occur 6 Boundary Conditions Location Region Name Slope Surface Surface 1 Segment Boundary Conditions Surface Boundary Conditions x ii Boundary Condition Description Units Boundary Name 0 000 0 000 Zero Flux eae 14 000 0 000 Concentration Expression 1 g m 3 24 000 0 000 Continue 24 000 27 000 Continue 14 000 27 000 Zero Flux Update Selected Segment length 14 000 m
22. new CHEMFLUX model 2 Apply appropriate boundary conditions in the CHEMFLUX model 3 Apply appropriate material properties The methodology for setting up the model is detailed in the following sections SOILVISION SYSTEMS LTD Sudicky Model 13 of 37 Model Description and Geometry Concentration 0 if t lt 5 then Concentration 1 Concentration 0 else Concentration 0 Zero Flux gt _ _ Zero Flux gt Zero Flux Material Properties The Material Properties for the numerical model are as follows Longitudinal Dispersivity o 0 5m Transverse Dispersivity Or 0 005m Diffusion Coefficient D 0 0423m7 yr 2 2 1 Model Setup In order to set up the model described in the preceding section the following steps are required The steps fall under the general categories of a Create model b Enter geometry c Specify boundary conditions d Apply material properties e Specify model output f Run model g Visualize results a Create Model A gradient file generated by SVFLUX is required for this example The seepage model described above has been included in the model files distributed with the SVFLUX software This file was generated previously in the Steady State SVFLUX model example When the solution for the model is finished a gradient file will be automatically created in the solution folder by SVFLUX The file is called gradient trn Th
23. nts VRVZ Y 14 t from 0 by 100 to 400 Yes SP Solution Mesh Y 14 att 400 Yes Mesh 3D att 0 Yes Add New Plot 2 Plot Settings Custom Plots 2 The toolbar at the bottom left corner of the dialog contains a button for each plot type Click on the Contour button tobegin adding the first contour plot The Plot Properties dialog will open SOILVISION SYSTEMS LTD A Three Dimensional Example Model 31 of 37 OO SA ER APE oO 11 12 13 14 15 Enter the title Concentration Select c as the variable to plot from the drop down Move to the Update Method tab If not already entered enter in the following values for Start 0 Increment 50 End 400 Move to the Projection tab Select Plane Projection option Select Y from the Coordinate Direction drop down Enter 14 in the Coordinate field This will generate a 2D slice at Y 14m on which the concentration contours will be plotted Move to the Output Options tab Select the PLOT output option Click OK toclose the dialog and add the plot to the list Repeat these steps 2 to 12 tocreate the plots shown above Note that the Mesh plot does not require entry of a variable under the Description Tab Also note that the Solution Mesh amp Mesh should only have an entry value for Start under the Update Method tab Time Steps Click OK to close the Plot Manager and return to the workspace There are a few output file types that can be specified to export
24. oduction 4 of 37 1 Introduction The Tutorial Manual serves a special role in guiding the first time users of the CHEMFLUX software through a typicalexample problem The example is typical in the sense that it is not too rigorous on one hand and not too simple on the other hand The Tutorial Manual serves as a guide by assisting the user with the input of data necessary to solve the boundary value problem ii explaining the relevance of the solution from anengineering standpoint and iii assisting with the visualization of the computer output An attempt has been made to ascertain and respond to questions most likely to be asked by first time users of CHEMFLUX SOILVISION SYSTEMS LTD Sudicky Model 5 of 37 2 Sudicky Model Sudicky 1989 developed the following example The model considers flow and solute transport in a heterogeneous cross section with a highly irregular flow field dispersion parameters that are small compared with the spatial discretization and alarge contrast between longitudinal and transverse dispersivities Zheng and Wang 1999 Project ContaminantPlumes Model VanderHeijdeSS VanderHeijde Minimum authorization required FULL 1 10 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 X It is important to note that you will be analyzing the SVFLUX model before the CHEMFLUX model is completed 2 1 Steady State SVFLUX Model The completed model is present
25. og The user must then e select a Head Expression boundary condition type e enter avalue of 5 375 m Close the dialog The newly specified boundary condition will be displayed with symbols on the CAD window A summary of the boundary conditions for this model can be found in Appendix A d Apply Material Properties Material Properties must now be entered and applied to specific regions in the model The following steps are SOILVISION SYSTEMS LTD Sudicky Model 10 of 37 required in order to properly apply material properties 1 Open Materials Manager Model gt Material gt Manager 2 Click New This will open a new material record Enter ChemFlux Soil1 as the material name and click OK 3 Enter Properties The material properties for ChemFlux Soill must be entered Click the Properties button on the Material Manager dialog 4 Laboratory SWCC data Choose the SWCC tab and click the Data button Enter the material properties as found in Table 1 Click OK to accept the data entered and close the SWCC Data dialog 5 Fitting Laboratory data can be best fit with the Fredlund amp Xing 1994 soil water characteristic curve equation Fitting the curve can be accomplished through the following steps e Select Fredlund amp Xing as the fitting method from the SWCC drop down e Enter 0 351 in the field for Saturated VWC e Click the Properties button to set the properties of the Fredlund amp Xing fit e Click the Apply
26. ots This section will give abrief analysis foreach plot that was generated 0 The Mesh plot displays the finite element mesh generated by the solver The mesh is automatically refined in critical areas Right click on the plot and select Rotate to enable the rotate window SOILVISION SYSTEMS LTD A Three Dimensional Example Model ChemFlux Model Example3D lengthunit m timeunit day L Concentration on Y 14 33 of 10 19 30 11 9 07 FlexPDE 5 0 19 BR porao amao mn rreman ie 1 00 1 00 0 95 0 90 0 85 0 80 0 75 0 70 0 65 0 60 0 55 0 50 0 45 0 40 0 35 0 30 0 25 0 20 0 15 0 10 0 05 le 3 Tutorial_Example3D Cycle 27 Time 400 00 dt 48 350 p2 Nodes 5657 Cells 3458 RMS Err 6 2e 4 Integral 116 3422 37 In this contour plot it can be seen the concentration is equal to 1 at the reservoir and decreases to 0 at the river SOILVISION SYSTEMS LTD A Three Dimensional Example Model on Y 14 34 of 10 19 30 11 9 07 Water Gradients FlexPDE 5 0 19 SEE Saas ae SoS 10 SSasa Steet ieee Nant ttt tte 0 20 0 19 0 18 0 17 0 16 0 15 0 14 0 13 0 12 0 11 0 10 0 09 0 08 0 07 0 06 0 05 0 04 0 03 0 02 0 01 0 00 Tutorial Example3D Cycle 27 Time 400 00 dt 48 350 p2 Nodes 5657 Cells 3458 RMS Err 6 2e d 37 Gradient Vectors show both the direction and the magnitude of the flow at specific points in the model Vectors illustrate th
27. rial or be left as void A void area is essentially air space In this model all blocks will be assigned a material 1 Select Slope in the Region Selector 2 Select Model gt Materials gt Material Layers from the menu to open the Material Layers dialog 3 Select the 3D Tutorial Soil material from the drop down for Layer 1 4 Close the dialog using the OK button 5 Select Reservoir in the Region Selector 6 Select Model gt Materials gt Material Layerss from the menu to open the Material Layers dialog 7 Select the 3D Tutorial Soil material from the drop down for Layer 1 8 Close the dialog using the OK button e Specify Model Output The next step is to specify the data which will be generated by the finite element solver Both the graphs displayed by the FlexPDE solver as well as the output generated for the subsequent CHEMFLUX analysis must be specified There are many plot types that can be specified to visualize the results of the model A plots few will be generated for this tutorial example model including a plot of the concentration contours solution mesh and water gradient vectors 1 Open the Plot Manager dialog by selecting Model gt Reporting gt Plot Manager from the menu 6 Plot Manager Plots Paint Area Volume Flux Sections Surface Flux Other Title Variable Restriction Projection Update Method PLOT MONITOR KG Concentration c Y 14 t from 0 by 50 to 400 Yes Y Water Gradie
28. s written out for visualization in ACUMESH or for inputting to other finite element packages PLOT MANAGER The Plot Manager dialog is first opened to display appropriate solver graphs The next step in model setup is to specify the data which will be generated by the finite element solver Both the graphs displayed by the FlexPDE solver as well as the output generated for the subsequent CHEMFLUX analyses must be specified 1 2 N da 10 Open the Plot Managerdialog by selecting Model gt Reporting gt Plot Manager from the menu The toolbar at the bottom left corner of the dialog contains a button foreach plot type Click on the Contour button tobegin adding the first contour plot The Plot Properties dialog will open Enter the title Concentration Select c as the variable to plot from the drop down Select the Plot from the Output Option tab Move to the Update Method tab and enter a Start Time 0 a Time Increment 1 and an End Time 20 Move to the Zoom tab and enter X 100 Y 0 1 Width 100 and Height 6 6 Click OK toclose the dialog and add the plot to the list Repeat steps 2 to 8 tocreate the remaining plots Note that the Mesh plot does not require entry of avariable Click OK toclose the Plot Manager and return to the workspace 6 Plot Manager Plots Point Area Volume Flux Sections Boundary Flux Other Title Variable Restriction Update Method KG Concentration c t from 0 by 1 to 20 Z Water
29. the results of the model One will be generated for this tutorial example model a file to output the results to ACUMESH for advanced visualization 1 Open the Output Manager dialog by selecting Model gt Reporting gt Plot Manager from the menu The toolbar at the bottom left corner of the dialog contains a button for each output file type Click on the ACUMESH button to begin adding the output file The Output File Properties dialog will open The title will be entered automatically as User_Example3D_ACUMESH The variables c vx vy and vz appear automatically as the default Type in the following values under the Update Method Time Steps Start 0 Increment 40 End 400 Click OK toclose the dialog and add the output file to the list Click OK toclose the Output Manager and return to the workspace f Run Model The current model may be run by selecting the Solve gt Analyze menu option g Visualize Results The flow vectors for the current model may be visualized through the following steps SOILVISION SYSTEMS LTD A Three DimensionalExample Model 32 of 37 1 Open ACUMESH View gt ACUMESH menu option 2 Plot Contours Plot gt Contours 3 Model State States toolbar drop down 3 3 Results and Discussion After the model has finished solving the results will be displayed in the dialog of thumbnail plots within the CHEMFLUX solver Right click the mouse and select Maximize toenlarge any of the thumbnail pl
30. under the project and model listed below The user may open this model to run and display the results of the analysis The user can also recreate this model under a separate model file The seepage model shown below gives the model dimensions boundary conditions material properties and the final flow regime This is followed by step by step instructions on how to enter and solve the contaminant transport model SOILVISION SYSTEMS LTD Sudicky Model 6 of 37 Model Description Geometry Concentration 0 if t lt 5 then 1 else 0 Concentration 0 Normal Flux Expression 0 1 m yr gt j 40 6 447 80 6 393 0 _ gt 0 5 375 m gt h Gradient Gradient Zero Flux gt o Ui Yo Is 250 gt Zero Flux gt lt Gradient 0 gt Material Properties Material Properties used for the SVFLUX steady state model are as follows Ka 158 m yr Ka 3156 m yr k ratio 1 0 k ratio 1 0 Volumetric water content 0 35 Volumetric water content 0 35 The soil water characteristic curve data was used for both materials in the model Table 1 Soil Suction kPa Volumetric Water Content Ratio 0 0001 0 351 100 0 35 1000 0 349 NOTE Steady state seepage solutions do notrequire that the soil water characteristic curves have an initial positive slope An initial positive slope is only required i
31. ust be selected by clicking on the region 2 Enter Boundary Conditions The Boundary Conditions dialog may be displayed under the Model gt Boundary Conditions gt Properties menu option Once in the dialog the user needs to e selectthe point 0 0 from the list e from the Boundary Condition drop down select a Flux Expression boundary condition equal to 0 e repeat these steps to define the boundary conditions for the remaining Region 1 segment as shown in the diagram and in the screen shot above be sure to define a Flux Expression boundary condition equal to 0 for the last point in the list 3 Close the dialog The newly specified boundary condition will be displayed with symbols on the CAD window d Apply Material Properties Material Properties must now be entered and applied to specific regions in the model The following steps are SOILVISION SYSTEMS LTD Sudicky Model 16 of 37 required in order to properly apply material properties 1 2 Open Materials Manager Model gt Material gt Manager Click New This will open a new material record Enter Properties Move to the Dispersion tab Enter the Longitudinal Dispersivity aL 0 5 m Enter the Transverse Dispersivity aT 0 005 m Select Constant as the Diffusion option Enter the Diffusion Coefficient D 0 0423 m2 yr Apply toregions open the regions dialog selecting Model gt Geometry gt Regions from the menu e select Material 1 from the drop down as the materi
32. w using the Draw gt Region Polygon from the menu Draw Region 1 The firstregion can be extended by drawing the geometry on the CAD window using the Draw gt Region Polygoncommand Alternatively the Region Properties dialog can be opened Model gt Geometry gt Region Properties and the region polyline points cut and pasted in from the provided spreadsheet The points can also be pasted into the dialog See Appendix A for the geometry of eachregion SelectRegion 2 The user must selectRegion 2 as the currentregion on which to draw Draw Region 2 The second region can be entered in a manner similar to that explained for Region 1 Repeat for Region 3 c Specify Boundary Conditions Flow models must generally have a defined entry and exit point for water to flow The boundary conditions shown at the start of this model can be entered using the following steps 1 2 Select Region 1 Region 1 must be selected by clicking on the Region Enter Boundary Conditions 1 The Boundary Conditions dialog can be displayed under the Model gt Boundary Conditions gt Properties menu option Once in the dialog the user needs to e select the node point 250 5 375 e then select anormal flux expression from the combo box e enter avalue of 0 1 m yr e select the point 0 6 5 select the Zero Flux expression from the boundary condition drop down box Enter Boundary Condition 2 The node 250 0 must be selected in the Boundary Conditions dial
Download Pdf Manuals
Related Search