03_EFDC_Model_Overview

EFDC Model TrainingDrew Ackerman, Cardno ENTRIXSeptember X, 2013
EFDC Model Background•Solves the three-dimensional, vertically hydrostatic, free surface, turbulent averaged equations of motions for a variable density fluid•Sigma coordinates and Cartesian or curvilinear, orthogonal grid•Assumes incompressible flow and hydrostatic pressure distribution with dynamically coupled salinity and temperature transport
EFDC Model Background, Continued•Simulates:•Wetting and drying•Simulation of controlled flow structures•Vegetation resistance•Wave-current boundary layers and wave-induced currents•Embedded single port buoyant jet module
EFDC Model Linkage•Provides dynamics for sediment, water quality and toxics fate and transport•Sediment simulations–Cohesive/non-cohesive, deposition/resuspension, sediment load and bed loads•Water quality–Eutrophication model, sediment diagenesis•Toxics–Trace metals/organic hydrocarbons and sediment interactions
EFDC Versions•Tetra Tech version•EPA version•Original EFDC model by Hamrick•Last updated in 2002•Sandia Labs•Research version•Dynamic Solutions•Commercial version with GUI•Additional capabilities and support•Last updated in 2013
EFDC Application•Original model•EFDC (EPA) version•Expanded model•Dynamic Solutions (without GUI)•DOS window•Text Editor
Model Configuration•Model developed to support hydrodynamic and water quality simulations•Stage, temperature, salinity•Nutrient transformations, dissolved oxygen kinetics, and eutrophication•Focus of this training is temperature and salinity
Model Files
File Types•Model control file•Describes model and sets up its output•Model grid files•Defines model grid orientation, elevations•Model boundary/input files•Characterizes what is entering/leaving model domain
Model FilesEFDC.inp Control primary simulation control fileSHOW.inpFile controlling screen printWQ3DWC.inpWater quality simulation control file
EFDC Model
Model Grid FilesCELL.inp Grid cell typeCELLLT.inp Grid cell typeCORNERS.inpCorner coordinates for each cellDXDY.inp Grid cell dimensionsLXLY.inp Grid cell orientation
Initial Condition FilesBEDBDN.inp Sediment bed bulk densityBEDDDN.inp Sediment bed porosityBEDLAY.inp Sediment bed thicknessSALT.inp Salinity concentrationSEDB.inp Cohesive sediment in sediment bedSEDW.inp Cohesive sediment in water columnTEMP.inpTemperature
Boundary/Input FilesCWQSERxx.inpWater quality input filesGWMAP.inpDefines groundwater seepage zonesGWSEEP.inpDefines groundwater flow ratePSER.inp Time series tide boundary conditionQSER.inp Time series flowSDSER.inp Time series cohesive sedimentSNSER.inp Time series non-cohesive sedimentSSER.inp Time series salinityTSER.inpTime series temperature
Boundary/Input Files, ContinuedWQALGG.inpAlgal dynamicsWQBNENFLX.inpTime series benthic fluxWQBENMAP.inpMapping info for spatial benthic fluxesWQSETL.inpWater quality settling ratesWQWCMAP.inpWater quality kinetic zonesWSER.inp Time series wind
General File Comments•Discuss relevant files to hydrodynamic simulation and model grid•Free format•Sensitive to integer/real numbers•View in Notepad++
Model Grid
Model Grid Overview•Need detailed information describing model cells•Cell orientation and where it is in the world•Current grid has one open boundary to east defining oceanic conditions
Grid Development•Three options•GRIDGEN•Dynamic Solutions •Delft3d
GEFDC•Original model grid generation program•Define boundary points•Define model grid type•Specify grid relaxation parameters•Cartesian, curvilinear •DOS and text based•Outputs image of grid in DXF format
Dynamic Solutions•More flexible grid generation•Cartesian, curvilinear •Use geographic files to help define grid boundary•Can import other grids to EFDC grids•Deflt RGFGrid, Grid95, SEAGRID
Delft3d RGFGrid•Intuitive GUI •Good online support•Active user community•Online training videos
LTPR EFDC Model Grid•Lower Tar Pamlico River •6 sigma coordinate layer models•593 horizontal cells
LTPREFDC Model Grid
Grid Files
•Model consists of a series of quadrilateral cells•Orthogonal and curvilinear•Files define cells•Location•Orientation•ShapeEFDC Grid FilesCELL.inp Grid cell typeCELLLT.inp Grid cell typeCORNERS.inpCoordinates for each cellDXDY.inp Grid cell dimensionsLXLY.inp Grid cell orientation
CELL.inp
i1234567891011121234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000911000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000099999999999999999999100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000959959955555555555559099999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999995555555555555555555809555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555555570999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999555556000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000009555555555500000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000999999999940000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000093000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001314151617181234567890123456789012345678901234567890123456789012345678903500000000000000000000000000000000000000000000000000000000000034000000000000000000000000000000000000000000000000000000000000330000000000000000000000000000000000000009999900000000000000003200000000000000000000000000000000000000095559000000000000000031000000000000000000000000000000999999999959990000000000000000300000000000000000000000000000999555555555590000000000000000002900000000000000000000000000009555555555555900000000000000000028000000000000000000000000000095555555555599000000000000000000270000000000000000000000000000995559999999900000000000000000002600000000000000000000000000000995555900000000000000000000000025000000000000000000000000000000955559000000000000000000000000240000000000000000000000000000009555590000000000000000000000002300000000000000000000000000000095555900000000000000000000000022000000000000000000000000000000955559000000000000000000000000210000000000000000000000000000009555599000000000000000000000002000000000000000000000000000000095555590000000000000000000000019000000000000000000000000000000955555900000000000000000000000180000000000000000000000000000009555559900000000000000000000001700000000000000000000000000000095555559900000000000000000000016000000000000000000000000000000955555559900000000000000000000150000000000009990000099999999999555555559990000000000000000001400000009999995999999955555555555555555555900000000000000000013009999995555555555555555555555555555555559000000000000000000129995555555555555555555555555555555555555590000000000000000001155555555555555555555555555555555555555555900000000000000000010555555555555555555555555555555555555555559000000000000000000955555555555555555555555555555555555555555900000000000000000085555555555555555555555555595555555555555590000000000000000007555555555555555555555555555555555555555559000000000000000000655555995555995555555559995595999999999999900000000000000000055555999999999999559999909999990000000000000000000000000000004955590000000000999900000000000000000000000000000000000000000399599000000000000000000000000000000000000000000000000000000020999000000000000000000000000000000000000000000000000000000001000000000000000000000000000000000000000000000000000000000000
Grid Numbering
CORNERS.inp
DXDY.inp
LXLY.inp
Initial Conditions
•Specify model conditions at start of the simulation•Focus on salinity and temperatureInitial Condition FilesBEDBDN.inp Sediment bed bulk densityBEDDDN.inp Sediment bed porosityBEDLAY.inp Sediment bed thicknessSALT.inp Salinity concentrationSEDB.inp Cohesive sediment in sediment bedSEDW.inp Cohesive sediment in water columnTEMP.inpTemperature
Model Layers•Sigma (σ-stretched) layers•Layers expand and contract•Like an accordion•6 layers defined for the model•Lowest number on bottom123456
SALT.inp
TEMP.inp
Boundary Files
Boundary/Input FilesASER.inpAtmospheric condition fileCWQSERxx.inpWater quality input filesGWMAP.inpDefines groundwater seepage zonesGWSEEP.inpDefines groundwater flow ratePSER.inp Time series tide boundary conditionQSER.inp Time series flowSDSER.inp Time series cohesive sedimentSNSER.inp Time series non-coheisve sedimentSSER.inp Time series salinityTSER.inpTime series temperatureWQALGG.inpAlgal dynamicsWQBNENFLX.inpTime series benthic fluxWQBENMAP.inpMapping info for spatial benthic fluxesWQSETL.inpWater quality settling ratesWQWCMAP.inpWater quality kinetic zonesWSER.inp Time series wind
Atmospheric Conditions
Atmospheric Conditions•Describe atmospheric conditions•Processes affected•Heating•Dissolved oxygen•Rainfall•Wet deposition of nutrients•Photosynthesis
Atmospheric condition input file: ASER.inp•Define environmental conditions•Atmospheric pressure•Dry air temperature•Relative humidity•Rainfall•Evaporation•Solar radiation•Fraction cloud cover
ASER.inp
Wind input file: WSER.inp•Wind conditions•Speed•Direction
WSER.inp
Tidal Boundary
Pressure time series: PSER.inp•Define seaward tidal condition at the “open boundary”•Three options•Define elevations (current model configuration for calibration and validation years)•Define model tidal harmonics•Use regression models based on stage at Washington (used to set up 2007 and 2008)
Define Elevations•Current model configuration•Can capture storm surges•Need good local monitoring information for simulation period•Need time series for entire simulation period•Utilize monitoring data to define absolute or relative water surface elevation
PSER.inpMPSERTCPSERTAPSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changingthe input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion
Define Harmonics•Modeled data from other sources•JTides, ACRIC model•Estimate tidal harmonics for different constituents•M2, S2, N2, K2, O1, K1, Q1, M4, M6•http://www.unc.edu/ims/ccats/tides/tides.htm•Extrapolate observed stage at Washington to open boundary condition at the estuary mouth •Based on lag time and distance (Xu et al, 2008)•Can model any period once defined•More computationally efficient
Model Card 14/15
Model Card 17/18
Surface Flow Inputs
QSER.inp•Flow time series•Define flows into/out of model•River/stream flows•Point sources•Return flows•19 identified sources in the model
QSER FlowsNS = 1, Greenville + GUC WTP IntakeNS = 2, Chicod Creek flow (cfs)NS = 3, Grindle Creek flow (cfs)NS = 4, Tranters Creek flow (cfs)NS = 5, Greenville WWTP Flow (cfs)NS = 6, Washington WWTP Flow (cfs)NS = 7, Belhaven WWTP Flow (cfs)NS = 8, PCPS WWTP Flow (cfs)NS = 9, ChocotowingNS = 10, BlountsNS = 11, DurhamNS = 12, SouthNS = 13, GooseNS = 14, BathNS = 15, Pungo CreekNS = 16, PantegoNS = 17, Pungo CanalNS = 18, GUC WTP Raw Water Withdrawal Rate (cms)NS = 19, GUC WTP Return Rate (cms)
Sources
QSER.inpMQSERTCQSERTAQSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion
Groundwater
Groundwater Flows•Different flows into model grid•Define zones and constant inflow rates•4 zones with 2 inflow rates (m/d)
GWMAP.inp
GWSEEP.inp
Temperature and Salinity
SSER.inpMCSERTCCSERTACSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion
TSER.inpMCSERTCCSERTACSERRMULADJADDADJNumber of data pointsMultiplying conversion factor changing the input time units to secondsAdditive time adjustmentMultiplying conversionAdditive conversion
EFDC Input Files
EFDC Control File•Main control file: EFDC.inp•Controls•Grid definitions•Inputs•Time steps•Output•Model calibration parameters•Arranged in “cards”
Model Execution
Simulations on Citrix Server•Run EFDC model on appserver.ncwater.org server•Same log in information as with OASIS modeling
Virtual DesktopLogging onto the Citrix server allows for virtual desktop access to the model
Virtual Desktop
Model SimulationRun model from icon on desktop•Runs base case model•D:\Tar_users\dackerman\EFDCBaseCase
Scenario Development
Example Scenarios and Required Input Changes•Alteration of conditions for a pre-developed year•Development of a new simulation year•Simulation of sea level rise•Simulation of potential impacts of climate change
Alteration of Conditions for a Pre-developed Year •Potential model years include 2001, 2003, 2007, 2008•Key parameters/factors could be altered•Flow•Temperature•Salinity•Files to change•QSER.INP •TSER.INP •SSER.INP
Development of a New Simulation Year•Key parameters/factors•Flow, temperature, salinity, boundary conditions, atmospheric conditions, initial conditions•Files to change•QSER.INP, TSER.INP, SSER.INP, ASER.INP, WSER.INP, PSER.INP, TEMP.INP, SALT.INP, DXDY.INP•It will be important to process all data to develop accurate representation of new conditions•This will be a significant effort to undertake
Simulation of Sea Level Rise•Key parameters/factors•Boundary conditions, initial conditions•Files to change•TSER.INP, SSER.INP, TEMP.INP, SALT.INP, DXDY.INP, PSER.INP, EFDC.INP•Requires changing seawater elevation but also salinity and temperature
Simulation of Potential Impacts of Climate Change•Potentially hotter/cooler or dryer/wetter conditions •Key parameters/factors•Flow, temperature, salinity, boundary conditions, atmospheric conditions, initial conditions•List files to change•QSER.INP, TSER.INP, SSER.INP, ASER.INP, WSER.INP, PSER.INP, TEMP.INP, SALT.INP, DXDY.INP, EFDC.INP•Climate change will impact terrestrial and oceanic conditions•There will be a large degree of uncertainty and care should be taken with these simulations
Model Scenario Runs
Scenario Simulations on ServerCopy EFDCBaseCase directory into a new directoryModify input files in new directoryRun simulation in new directory
Model Output
EFDC Output•EFDC example output file•Copy post-processing files (PostProc.exe) into directory where simulation was run•Copy post-processing Excel file (PostProcessing.xlsx) onto your computer
Your local computerCitrix Server
Data from EFDC Output•Post-process results on server•Will speed up download times (output files are large)•Copy “PostProc.exe” into scenario folder on the server•From \EFDC_PostProcessing
Data from EFDC Output•Extract data from salinity and temperature output files•16 stations specified in EFDC.INP file–Card 87•SALTSxx.out and TEMTSxx.out•Run the executable file to write a file that averages the top, middle, and bottom two layers (temperature and salinity)•Double click on PostProc.exe•To generate TempSalt.csv
Example model outputadfdfafafaf>werwreqwrewreafdafasfdasfdas>adffasfdasfafdsa
Calculate Percentiles•Steps to calculate percentiles1.Double click “PostProc.exe”2.Open PostProcessing.xlsx on your computer3.Copy TempSalt.csv from the scenario folder to your computer4.Open TempSalt.csv in EXCEL5.Select columns B through J and copy6.Paste into Cell B1 in the “Data Import” tab7.Results are shown in “Temp Salt Percentiles”
PostProcessing.xlsxOpen TempSalt.csv in EXCEL and paste those columns into Columns “B” through “J”Percentiles for the top two layers, middle two layers, and bottom two layers are calculated in the “Temp Salt Percentiles” tab
Temp Salt Percentiles
Additional Model Post Processing
Advanced Model Post Processing•Requires model user to write post processing code•FORTRAN, C++, BASIC, Python, Matlab•Will require more knowledge of model output structure than has been presented here
Issues with direct linkage to the OASIS model
Constraints to Linking the LTPR Model to the OASIS Model•LTPR model has been set up for four specific years •Headwater flows, open boundary condition, atmospheric conditions, withdrawals, and discharges•OASIS model is a time-series model developed for a longer period•Changes to the LTPR model must be hard wired into input files •Cannot read output files from OASIS directly•Running additional years with the LTPR model will require reconfiguring the model boundary and forcing files
Questions?