USGS
South Florida Information Access
SOFIA home
Help
Projects
by Title
by Investigator
by Region
by Topic
by Program
Results
Publications
Meetings
South Florida Restoration Science Forum
Synthesis
Information
Personnel
About SOFIA
USGS Science Strategy
DOI Greater Everglades Science Plan
Education
Upcoming Events
Data
Data Exchange
Metadata
publications > paper > simulation of integrated surface-water/ground-water flow and salinity... > numerical implementation

Simulation of integrated surface-water/ground-water flow and salinity for a coastal wetland and adjacent estuary

3. Numerical implementation

Home
Governing Equations
>Numerical Implement.
- Coupling Procedure
- Leakage Calculation
- Surface/Subsurface Solute Exchange
- Rainfall & ET
- Rewetting
Application of Model
Discussion
Summary & Conclusions
Acknowledgments & References
Tables, Figures, Eqs.
PDF Version
To solve the coupled surface-water and ground-water equations presented in the previous section, the finite-difference programs, SWIFT2D and SEAWAT, were modified to run timesteps sequentially under the control of a master program called FTLOADDS (Flow and Transport in a Linked Overland/Aquifer Density Dependent System). The SWIFT2D two-dimensional hydrodynamic flow and solute-transport code was originally developed for bays and shallow estuaries (Leendertse and Gritton, 1971; Leendertse 1987). The code has been applied to Jamaica Bay, New York (Leendertse, 1972), to Delta Works, The Netherlands (Leendertse et al., 1981), to Tampa, Florida (Goodwin, 1987; Goodwin 1991), to Pamlico River Estuary, North Carolina (Bales and Robbins, 1995), to Charlotte Harbor, Florida (Goodwin, 1996) and to the Neuse River Estuary, North Carolina (Robbins and Bales, 1995). The SWIFT2D program was later modified by Swain et al. (2004) to represent overland flow in coastal wetlands and to include the effects of spatially distributed rainfall and evapotranspiration. SWIFT2D uses a finite-difference approximation to solve the surface-water equations (Eqs. 1-5). SEAWAT, a combined version of MODFLOW (McDonald and Harbaugh, 1988) and MT3DMS (Zheng and Wang, 1999), was designed to solve the three-dimensional variable-density ground-water flow and solute-transport equations (Eqs. 5, 6, 7) using finite-difference methods (Guo and Bennett, 1998; Guo and Langevin, 2002; Langevin et al., 2003). Examples of SEAWAT applications include simulation of submarine ground-water discharge to a marine estuary (Langevin, 2001, 2003) and intercode comparisons (Bakker, 2003; Bakker et al., 2004).

SWIFT2D uses an alternating-direction implicit (ADI) method and a space- and time-staggered grid to solve the governing equations, such that each surface-water timestep is divided into two half timesteps-one half timestep for flow and transport in the x-direction and the other for the y-direction. In each of the two phases of the ADI method, the continuity equation and one of the components of the momentum equations are solved with local storage (and corresponding transport term of the continuity equation), local acceleration, pressure gradient, and the frictional term of the momentum equation treated implicitly. The last three terms on the left-hand side of Eq. 1 (the source and sink terms) are not included in the finite-difference solution, but are separately added to, or subtracted from, the cell volume. SEAWAT uses an implicit finite-difference approximation to solve the ground-water flow equation (Eq. (6)), and contains several alternative methods for solving the solute-transport equation (Eq. (7)), including implicit and explicit finite-difference methods with various weighting options and the method of characteristics.

The integrated code for SWIFT2D and SEAWAT requires cells that coincide and are identical in size. The integrated code was designed such that the domains of the two models need not be identical, provided that leakage fluxes can be neglected in areas where the two models do not overlap. Although not used for the Everglades application, this feature may prove useful for certain applications where the extension of the model domain is necessary in only one of the two systems.

< Previous: Governing Equations | Next: Coupling Procedure >



| Disclaimer | Privacy Statement | Accessibility |

U.S. Department of the Interior, U.S. Geological Survey
This page is: http://sofia.usgs.gov/publications/papers/swgwflow_salinity/numerical.html
Comments and suggestions? Contact: Heather Henkel - Webmaster
Last updated: 19 October, 2005 @ 10:03 AM(KP)