This tutorial document is based on a steady-state problem from Thomann & Mueller (1987) “Principles of Surface Water Quality Modeling and Control” Sample Problem 6.10 (which was itself adapted from Driscoll et al., 1981).
1. Tasks to be performed
This document will take the user through the following steps:
- building a grid
- setting hydraulic parameters
- setting initial conditions
- setting boundary conditions
- setting a conservative and non-conservative dye tracer
- creating data files for calibration data
- configuring calibration plots
2. Description of the Problem
The problem is taken from Thomann and Mueller (1987, p.356). A city of approximately 60,000 people discharges its wastewater into a relatively small river with an average annual flow of about 250 cfs. The city’s wastewater is presently treated by a trickling filter plant which provides about 85% BOD removal and has reached its design capacity of 7.5MGD. The population is expected to increase by 50% to 92,000 people by the year 2000. Expansion of the treatment plan to the capacity of 11.5 MGD and provision of an activated sludge system for more efficient secondary treatment is proposed.
A calibration analysis at Qu = 100 cfs is done first to establish kinetic rated based on observed in stream data.
2.2 Basic Data and Data Processing
2.2.1 Data Requirements
- Upstream and tributary boundary inflows of water and pollutants
- Initial conditions of water depth and pollutants in time and space
- Physical transport of water (velocity, travel time)
- Waterbody geometry (length, depth, area, volume)
- Pollutant fate reactions and interactions
- Sources/loading rates of flow & pollutants
2.2.2 River Characteristics
Relatively constant spatially overreach of interest. From three field surveys of different flows, the following area flow and depth flow relationships were derived:
- A = 19.5Q0.6 ; A (ft2), Q (cfs)
- H= 0.312Q0.5; H (ft), Q (cfs)
- W = A/H
- USGS gaging station upstream of WWTP
- Long term average summer flow = 100 cfs
- 7Q10 = 30 cfs (for min 7-consecutive day flow with a return period of once in 10 years)
Figure 1. 1-D River schematic.
- From the time of travel studies at various flows, U = 0.0513Q0.4; U (fps), Q (cfs)
- 25 oC during the low flow summer survey, 27 oC maximum monthly average temperature (August) to be used with 7Q10 to check conformance with standards.
Stream @ x = 0
Q = 100 cfs
Qe = 7.5 MGD (11.6 cfs)
Q = 111.6 cfs
CBOD5 = 1.0 mg/l
CBOD5 = 40 mg/l
CBOD5 = 5.05 mg/l
NH3-N = 0.2 mg/l
NH3-N = 15 mg/l
NH3-N = 1.74 mg/l
DO = 8.3 mg/l (sat.)
DO = 8.3 mg/l (sat.)
DO = 8.3 mg/l
CBODU = 2.0 x 5.05 = 10.1 mg/l
NBOD = 4.57 x 1.74 = 7.8 mg/l
H = 0.312 x (111.5)0.5 = 3.3 ft
U = 0.0513 x (111.6)0.4 = 0.34 fps = 5.6 mpd = 10.32 cm/s
Ka = [12.9 x 0.340.5 / (3.3)1.5] x (1.024)25-20 = 1.42 /day @25 oC
2.3 Model Input Data
- Located at x = 227,030 m
- Upstream inflow Qu= 2.83 cms (100 cfs)
- Upstream bottom elevation = 303 m
East Branch inflow (note that this inflow can be ignored and is not in the original example)
- Located at x = 216,705.5 m
Waste Water Treatment Plant
- Located at x = 213,528.7 m
- Qw = 0.33 cms (7.5 mgd, 11.6 cfs)
- Located at x = 165,083 m
- Downstream bottom elevation = 300 m
- Flow below WWTP = 3.16 cms (100 + 11.6 =111.6 cfs)
Area-Depth-Velocity Relationships as f (flow, cfs)
- Area = 19.5*Q0.6 =330.1 ft2 = 30.67 m2
- Depth = 0.312*Q0.5 = 3.3 ft = 1.006 m
- Velocity = 0.0513*Q0.4 =0.34 ft/s = 0.1036 m/s
- Width = Area/Depth = 100.15 ft ~ 30.0 m
- Grid cell length DY=1588.4 m
- Total length of 1D river = 61,947.0 m
- Number of Grid Cells = 39
- Constant head for downstream elevation control = 301 m.
User experience and trial and error method to determine appropriate slope and bottom roughness
- Bottom roughness: Z0 = 0.06 m
- Channel slope = - 0.5 e-4 m/m
- Cell depths computed by EFDC model from channel slope and upstream & downstream bottom elevations
3. Generating a New Model
3.1 Generate a New Model
3.1.1 Generating a Grid
In this step, we will generate a simple model grid to present the 1D river in our problem.
1. Open EFDC_Explorer.
2. From the empty main form shown in Figure 2 select the New Model button or press Ctrl + N from the keyboard.
Figure 2. EFDC_Explorer main form.
The form is shown in Figure 3 should now be displayed.
Figure 3. Generate the EFDC model.
3. Select Generate Uniform Grid for this case
4. We will now enter values for Lower Left and Upper Right of the model domain. These values correspond to the US and DS boundary conditions outlined in Model Input Data above
5. Set Lower-Left (m): (0; 165,083) in X and Y direction respectively
6. Set Upper-Right (m): (30; 227,030) in X and Y direction respectively
7. Cell Size (m): 30 and 1588.4 for X and Y respectively
8. Number of Cells: Press the Calculator symbol to update the number of Cell in X and Y directions.
9. Leave Rotation Angle = 0 as we won’t rotate the grid then
10. Press Generate button as shown in Figure 4. Note that EE creates two extra cells outside of each edge of the active cells, so the grid is (2+1+2) x (2+39+2) = 5 x 43
11. Press OK button, ignore the messenger box about UTM Zone as this is not in an actual physical location.
Figure 4. Generate EFDC model – with data.
12. In the Main Form, the General Information of the model is now displayed (Figure 5)
3.1.2 Naming and Saving the Model
1. The user should enter the name of the project by RMC on Descriptions in General tab. These will become the default names in the View Plan and Time Series titles. Run ID, Run Title, Project Title as shown in Figure 6. The user also can keep track of changes from model to model by RMC on Change Logs and updating information in the form as shown in Figure 7.
Figure 6. Model Description form – General tab.
Figure 7. Change Logs form – General tab.
2. Save the model by selecting Save button from main menu toolbars and create a new directory and then, press OK button. (See Figure 8)
3.1.3 Viewing and Checking the Grid
Once the model grid is created and saved, we can display and check the grid.
1. Click the 2DH button on the main toolbar , or the user can go to 2DH View on Main Menu and select New 2DH View, to see the model in plan view.
2. Enable Bottom Elevation by turning on the light bulb icon from the View Layer Control menu as shown in Figure 9.
This now displays the bottom elevation of the model domain. As the model is very long and narrow, it is hard to see clearly. To increase the width of the model, the user must use horizontal exaggeration.
Figure 9. 2DH View.
5. RMC on the legend to bring up the 2DView Settings as shown in Figure 10. From this form, the user can set the Easting value to 100 in the E/N Scales Modifiers frame. This will stretch the model by a factor of 100 in the Y direction and is easier to view.
Figure 10. 2DH View – 2DH View Settings.
Now we want to better understand the cell indices of our model.
8. Add Cell Indices layer by LMC on this button or RMC on the blank area in Layer Control and select Add new layer ... It will show View Options form, in Primary Group select Model Grid, in Parameters select Cell Indices, in Cell Indices Options choose Show I J's and set the value of I,J Skip to "1", then click on Add button to finish and close the windows (Figure 11)
Figure 11. Add Cell Indices layer to 2DH View.
9. Click on any of the cells in the model domain to display an information box (yellow) that describes the cell L number, I and J number, the X, Y coordinates, bottom elevation, and depth.
Figure 12 2DH View – cell indices.
3.2 Assigning Initial Conditions
In this step, we will set up the initial conditions for bottom and water surface elevations.
3.2.1 Assigning Bathymetry
1. Return to the 2DH View and enable edit for Bottom Elevation by clicking on the pen button shown in Figure 13.
Figure 13. Enable edit for Bottom Elevation.
2. LMC on Grid Cell Tool button and click on Set IJ Slopes (Figure 14) to open the Assign Slope form.
3. In Assign Slope windows, set Start Value = 300 and J Slope = - 0.00005 then press OK. Note that this is a negative slope to assign from the DS to US end of the model, starting at 300 m at the DS end. (Figure 15). The model's bathymetry is updated and shown in Figure 15.1 after clicking on OK button on the Figure 15.
Figure 14. Open Assign Slope windows.
Figure 15. Assign data for model bathymetry.
Figure 15.1 Model bathymetry after data Assignment.
3.2.2 Assigning Water Surface Elevation
Return to the Main Form and select the Initial Conditions tab to set water surface elevations as shown in Figure 16.
1. RMC on Water Depth/Elevation then select the Assign Depth to display the form shown in Figure 17.
2. Click Assign Depth button
3. Use Constant option and set the constant = 1.
4. Click Apply Defined Conditions
5. Click OK button
Once the Model Bathymetry and Water surface elevation are assigned, we can check the data using the viewing options (2D plan view or the longitudinal profiles) by following the steps:
To check the data, return to the Main Form and select the Longitudinal Profile Plot button The Data Extraction form is shown in Figure 18 should now be in view. Using this form the user may select a profile along a pre-defined polyline (Drape file) or along an I or J entered by the user.
In this case, choose Use I and set Specific I = 3, then move to Define Parameter to Plot
To add Bottom Elevation to the plot, in Primary Group choose Bottom, in Parameter choose Bottom Elevation and press Add to plot button
To add WSEL to the plot, in Primary Group choose Water level, in Parameter choose Water Surface Elevation then press Add to plot and OK button. This displays the plot shown in Figure 19.
Figure 18. Data Extraction for Longitudinal Profile Plots.
Figure 19. Longitudinal Profile.
3.2.3 Assigning Bottom Roughness
1. Similarly with WSEL, RMC on Roughness in Initial Conditions to open Bottom Roughness windows. LMC on Assign to display the form to assign the value (Figure 20)
2. Choose Constant and set the value of Bottom Roughness = 0.06
3. Click on Interpolate to assign and Ok to close the form
3.2.5 Placing Labels in the 2DH View
We now want to create labels for this model.
- Open a text editor and create a simple text file with three columns, space delimited. X, Y, Z, and Label are the four columns.
X and Y correspond to the locations that we want our boundary conditions as shown in Figure 21.
Figure 21. Creating a labels file.
2. Save the file with the suffix .dat or .txt
3. Proceed to the ViewPlan 2D form
4. Go to Import an External Overlay Layer button at the bottom Layer Control windows then select the label file and remember to choose the type of import file as Label Files (*.lbl, *.dat, *.txt),
5. The label file is loaded and show as a layer on the 2D view.
When displayed the labels are as shown in Figure 23 and have not yet been correctly aligned and symbols have not yet been selected.
Figure 22. Loading Label File.
Figure 23. 2DH View: labels.
7. RMC on the label layer just loaded and select Properties to correctly align the labels and set the symbols as shown in Figure 24 and Figure 25.
Figure 24. 2DH View: Edit Labels.
Figure 25. 2DH View: symbol properties.
8. Update the symbols and labels as desired such as shown in Figure 26.
Figure 26. 2DH View with labels and makers.
9. Select Save
It is assumed for the purposes of this tutorial that three data recording stations are providing measured data. The user should create a file to record the locations of data recording stations called 1DRiverStations.dat (Figure 27)
Import this file in the same way as above to have a new layer for this data recording stations label file (Figure 28).
Figure 27. 1DRiverStations dat file.
Figure 28. Loading data posting file.
Figure 29. Overlay of data posting file on labels.
4 Boundary Conditions
The boundary conditions for this model have been taken from the Thomman textbook example and are summarized in the table below. In this case, dye is used as a proxy for the contaminant transport.
To enter these BCs into the model we will first create the time series and then set the BC locations and link them to the time series.
Proceed to the External Forcing Data Tab shown in Figure 30
From this form, we can create the Flow boundary time series and the Dye boundary time series
Figure 30. Main Form – boundary condition settings.
4.1. Flow boundary time series
1. RMC on Flow and select Add New Data Series. This will display the form shown in Figure 31.
2. Set the number of flow series to 3 (since Downstream will be set as a WSEL or pressure boundary) by clicking on Add New button 3 times
3. Modified the value of each flow series by choosing them in Select Series or type the number of series base on the corresponding series on the beside text box
4. Set the title to “Upstream” and type or copy and paste into the form for the time and flow data. The value of "200" represents the start time of the data series in Julian time. At least two times must be entered for each series. The values are entered shown in Figure 32
5. Repeat this process for the WWTP series in Figure 33
6. Repeat this process for the East Branch series, in Figure 34.
7. Click on OK button to finish editing Flow time series for boundaries.
Figure 31. Boundary condition settings: flow data series.
Figure 32. Boundary Condition: US flow.
Figure 33. Boundary Condition: WWTP flow.
4.2. Dye boundary time series
1. Return for the Main Form shown in Figure 35 and enable Dye module by RMC on Modules and ticks on the checkbox Dye then press OK to return again to Main Form
Figure 35. Active Dye Module.
The Dye module will be added to the Modules, RMC on the Dye then select Settings, the form Dye Module will appear as shown in Figure 35.1 from this form, click on the drop-down list in the Type column then select
Figure 35.1 Dye Module – Conservative Dye setting.
2. Proceed to the External Forcing Data Tab. Now, you can see Dye is already shown in the drop-down menu
3. RMC on Dye and select Add New Data Series.
4. Repeat the above process for dye concentrations for the Upstream (US) as shown in Figure 36
5. Create time series for dye concentration for EastBranch (EB) with the same values as for “Upstream”, ie 0.2 (Figure 37)
6. Create time series for “WWTP” boundaries as shown in Figure 38.
Figure 36. Boundary Condition: upstream dye concentration.
Figure 37. Boundary Condition: East Branch dye concentration.
Figure 38. Boundary Condition: WWTP dye concentration.
4.3. Linking BC locations to the time series.
- In the Layer Control, turn on the View / Edit function of Boundaries, LMC on light, mouse, pen buttons
- RMC on the cell most upstream cell and select Add New Boundary Group to create a new boundary condition.
- Name the Boundary Group Name with ID “Upstream” and select Boundary Group Type = Flow Boundary as shown in Figure 39.
Figure 39. 2DH View – Assign Boundary Conditions: Upstream
The Flow Boundary Conditions will now be displayed. This form is used to link the flow table to the BC.
6. Select Time-Varying option and use the drop-down table under the Forcing Assignment | Flow Table to select the “US” table.
7. In the Boundary Group Conditions | Time-Varying Concentration | Constituent | Dye, select US table from the drop-down in the Data Series column (Figure 40)
8. Click Edit (E) to check which series have been defined if necessary.
Apply similar steps for EastBranch and WWTP boundaries as shown in Figure 41 and Figure 42.
Now we will set the Downstream BC. This will use a fixed WSEL to set the pressure boundary (clamped boundary). Therefore, we set a constant head at the downstream boundary.
- RMC on the bottom-most cell (I=3, J=3)
- Name the Boundary group Name “Downstream”. Click OK
- Select South Open Boundary for Boundary Group Type as shown in Figure 43. Click OK
- After click OK button, the Open Boundary Conditions form appears as shown in Figure 44, select Time-Varying in Forcing Approach, in the Level Data click the Edit button, the Boundary Data Series will now be displayed as shown in Figure 45.
- Click on Add New button to create new boundary data series
- Give the name as “DS” for Series Name then press Enter
- Enter time and WSEL as shown in Figure 45 and click OK.
- Returning to the Open Boundary Conditions form. Select the new “DS” series from the drop-down menu then click OK as shown in Figure 46. Note that dye is not set to the Downstream boundary
All boundary conditions have now been set as shown in Figure 47.
Figure 43. Downstream Boundary Conditions.
Figure 44. Set BC for Downstream.
Figure 45. Create WSEL for Downstream Boundary.
Figure 46. Link WSEL time series to the downstream open boundary.
Figure 47. 2DH View: Boundaries.
5 Running the Model
We want to set the model to run from day 200 for 20 periods of 24 hours at a time step of 100 seconds as shown in Figure 48.
- Save the model
- Select the Timing/Linkage tab and RMC on Timing
- Set the Model start time to 200;
- Set # of Reference Periods to 20;
- Set Duration of reference periods to 24;
- Set Time step (seconds) to 100.
Figure 48. Main Form – Time / Linkage: model run timing.
7. Click Run EFDC button from the main toolbar, the EFDC+ Run Options form appears, navigate to the General tab and set the Output Interval (min) to 15 minutes, select the version of EFDC+ executable file by browsing to the file locates in the default installation folder (e.g C:\Program Files\DSI\EEMS10.2) as shown in Figure 49.
8. Click the Run EFDC+ button on the form to start run the model as shown in Figure 50.
Figure 49. EFDC+ Run Options: Set Output Interval.
Figure 50. EFDC+ Running window.
6 View Model Result
To View dye in 2DH View:
- When the model has finished running. Click on the Refresh Model Output button ( ). Click the 2DH button from the main toolbar to view the model results in the 2DH View.
- Add New Model Layer from Layer Control Select Primary Group: Water Column; Parameter: Dye; Option: Depth Average
- To see the change of dye, move the timing bar or using short-key Ctrl+G to go to a specific date then click OK button
- Click on a grid cell, the dye concentration will be shown in the yellow rectangle at a specific time as shown in Figure 51.
To view a time series of dye from selected grid cells:
- Select a cell (or multiple cells) first
- Click on the Time Series Plot button ( ) from the main toolbar, the Extract Time Series Data form appears as shown in Figure 52, click OK button, a time series of the dye at those cells will be displayed as shown in Figure 53.
Figure 51. View dye concentration in the 2DH View.
Figure 52. View dye time series of selected cells (1).
Figure 53. View dye time series of selected cells (2).
7 Model Result vs. Measured Data
7.1 Conservative Dye Case
The longitudinal profile comparison will be selected for this case. A profile plot may be set up using the Longitudinal Profile button
- Click the Longitudinal Profile button ( ) on the main toolbar, the Data Extraction for Longitudinal Profile form appears, select Use I option then enter 3 to Specify I.
- Select Primary Group: Water Column, Parameter: Dye, Options: Depth Average for Layer Settings then click Arrow button to add to Parameters to Plots, then click OK button as shown in Figure 54.
- After click OK button, the plot will be displayed, use Ctrl+G to go to t=218 then click OK button as shown in Figure 55.
- In order to compare to measured data, load the data file by clicking the Import button ( ) from main toolbar (data file name: DyeCX_T218_readMe.dat)
- RMC on the legend to set the series line options for the data to dots using Line Formatting options as shown in Figure 56. As a result, the comparison plot is shown as Figure 57.
Figure 54. Define to plot longitudinal profile.
Figure 55. Longitudinal profile for I=3.
Figure 56. Plotting style .
Figure 57. Dye longitudinal profile compared with measured data at Time = 218.
7.2 Non-Conservative Dye Case
To compare with the theoretical results it is necessary to reset the dye as a non-conservative substance as shown below.
Figure 58. Decay of Non-conservative substances in Cartesian coordinates and Logarithmic form.
Base on the conservative dye model, we will set a new model for the non-conservative dye model by following steps:
- RMC on the Dye module
- Then select Settings, Dye Module will appear, from this form click the drop-down list in the Type column then select the Non-conservative option
- Set Dye decay Rate to 0.15 in the 1st-order Decay Rate (1/day) column as shown in Figure 59
- Save as model to a new model (e.g Run02)
- Run the new model again
The resulting longitudinal profile of dye for I=3 at time = 218 will be shown in Figure 60
Figure 59. Dye Module – Non-conservative Dye setting.
Figure 60. Longitudinal Profile: K=0.15 Non-conservative dye.
Now we use the data from Figure 61 (provided in the file “River1D-Analytical Solution.dat”) as longitudinal calibration data for this non-conservative model.
- Click Import Data Series button from the main toolbar shown in Figure 62.
- The Open form will appear, browse to folder contains the file “River1D-Analytical Solution.dat” then select that file
- Select “2” for the value for the analytical solution then click OK button as shown in Figure 63
- Use CTRL + G to go to time = 218, the longitudinal plot for dye comparison is shown in Figure 64
Figure 61. Longitudinal calibration Data.
Figure 62. Longitudinal calibration Data.
Figure 63. Longitudinal calibration Data.
Figure 64. Longitudinal plot: Data vs. modeled.
Time Series Comparison
We will now compare measured data at three stations, RIV01, RIV02, and RIV03 described earlier (Figure 29). Data for these stations is as shown below in Figure 65 and
- Return to Main Form | Model Analysis | Calibration Time Series (Figure 66)
- RMC to Calibration Time Series then select Define Calibration Series
Figure 65. RIV01 calibration data.
Figure 66. Model Analysis: Calibration Time Series.
3. Set the Number of time series to 3 then press enter, three rows will be added in the below table of the form as shown in Figure 67
4. To fill in the grid table, the user can enter the values to each column or RMC on any cell on the table, the Station Information form will appear, from this form the user fills in the station's information as coordinates, parameter to compare (e.g Dye) and data file path (Figure 68). The user can click on the Get I & J button to see the station locates in which cell of the model. Finally, click OK button.
5. Enter group number and tick on Check box for all stations in the User column (unchecked box will be ignored to compare) as shown in Figure 69. Click OK button.
Figure 68. Model Analysis: Station Information.
Figure 69. Update information of the three measured stations completely and checked boxes to show plots.
6. Click on the form of Figure 69 to show the plot, the user can edit the time series of data from line to point and reset Y scale of the plot. Figure 70 shows modeled compared to data at RIV01. To go to the next calibration plot, the user can click on the Next Calibration Plotting button on the top main toolbar as shown in Figure 70
Figure 71 and Figure 72 show modeled compared to data at RIV02 and RIV03.
To save the calibration stations configured in the model, the user should save the model then select Save Profile File Only (EFDC.EE) option as shown in Figure 73. So next time when the user loads the model again, the calibration configuration is still existing in the model.
Figure 71. Time series comparisons: RIV02.
Figure 73. Save Model to save calibration configuration.