Hwai-Ping Cheng

Integrated Modeling of the Battlespace Environment

The goal of the Battlespace Environments Institute (BEI) is to integrate Earth and space modeling capabilities into a seamless, whole-Earth common modeling infrastructure that facilitates interservice development of multiple, mission-specific environmental simulations and supports battlefield decisions, improves interoperability, and reduces operating costs.

Modeling surface and subsurface hydrologic interactions in a south Florida watershed near the Biscayne Bay

Restoration of the South Florida ecosystem is a major undertaking for the U.S. Army Corps of Engineers and the South Florida Water management District. The Biscayne Bay Coastal Wetlands (BBCW) Project is one component of more than 60 restoration plans and has a goal to restore or enhance freshwater wetlands, tidal wetlands, and near shore bay habitat. The primary purpose of the BBCW project is to redistribute runoff from the watershed into the Biscayne Bay, away from the canal discharges that exist today and provide a more natural and historical overland flow through the existing and/or improved coastal wetlands. In an effort to restore wetlands, several structures, and management plans and scenarios are considered. One of the plans is to deliver fresh water from the existing canals through a shallow spreader swale system that is to distribute fresh water through wetlands into the Biscayne Bay. To achieve this, a tool is needed to design this complicated shallow spreader swale system. This paper presents how a spreader swale system, which includes 1D canal network routing, 2D overland flow, 3D subsurface flow, and flow through the interface of any two sub-domains of the spreader system, is simulated with the WASH123D computer code. Brief physics-based mathematical statements and numerical strategies of the model will be given. A hypothetical example that uses topographic data for the project area will be provided to demonstrate how WASH123D can help the design of a spreader swale system. Some issues that concern the numerical convergence of the coupled flow model will also be discussed in this paper.

Parallelization of the WASH123D Code—Phase III: 1-Dimensional Channel, 2-Dimensional Overland, and 3-Dimensional Subsurface Flows

Watershed models are used to simulate and predict major hydrological processes, such as surface or subsurface flows, which may occur on different spatial domains and temporal scales. A key feature of watershed models is the ability to model interactions among different processes and domains. Such interactions can be strongly or weakly coupled depending on the relevant time scales for each process. WASH123D is a first-principles, physics-based model for simulating a coupled system of channel flow, overland flow, and subsurface flow. A parallel version, pWASH, has been used for the calibration, validation, and evaluation of proposed alternatives of the Biscayne Bay Coastal Wetlands project. The goal is to rehydrate wetlands using the best alternative based on the simulation results. In the pWASH code, channel flow is modeled as a one-dimensional (1-D) channel network, overland flow as a 2-D process, and subsurface flow as a 3-D process. Different algorithms are implemented to account for the interactions between these different domains. The pWASH code is designed to tackle large watershed problems on parallel high performance computers. A software-engineering approach was used to efficiently parallelize the complex coupling algorithms. The resulting software toolkit encapsulates the parallel data structures and message passing required for multi- domain/process interactions. In this paper, the authors briefly describe the numerical methods in WASH123D, the parallelization of pWASH, and the scalability of coupled flow problems running on parallel high performance computers.

A parallel software development for watershed simulations

A watershed software application is designed to model a coupled system of multiple physics on multiple domains. Tremendous computational resources are required to integrate the system equations on large spatial domains with multiple temporal scales among them. Supported by the Department of Defense Common High Performance Computing Software Initiative, the parallel WASH123D software development aims to efficiently simulate one aspect (i.e., soil and land) of the battlespace environment. Currently, the coupled two-dimensional overland and three-dimensional subsurface flows have been completed. Different numerical approaches are implemented to solve different components of the coupled system. The parallelization of such a complex system is developed on an IT-based approach—modular, hierarchical model construction, portable, scalable, and embedded parallel computational tools development and integration. Experimental results are presented to demonstrate the successful implementation of the parallel algorithms. Detailed profiling is also provided to show the imposed light-weight communication overhead.

Parallelization of the WASH123D Code—Phase II: Coupled Two-Dimensional Overland and Three-Dimensional Subsurface Flows

The parallel WASH123D is designed to simulate watershed systems using a coupled system of one-dimensional (1-D) channels, 2-D overland areas, and 3- D subsurface media on parallel scalable computers. The U.S. Department of Defense High Performance Computing Modernization Program funds the parallelization of the watershed model through the Common High Performance Computing Software Initiative in order to solve one aspect of the battlespace environment, which in- cludes space, weather, ocean, and soil, to develop a complete coupled battlespace environment. Tasks in this project include parallel algorithm development, software toolkit development, and performance studies. Currently, the parallel version, which includes coupled 2-D overland and 3-D subsurface flows, is completed and packaged for application projects, such as calibration and validation of the coupled 2-D/3-D Biscayne Bay Coastal Wetland inland flow model presented in this paper, on various computer architectures.

Parallelization of the WASH123D code—Phase I: 2-dimensional overland and 3-dimensional subsurface flows

The parallel WASH123D, which is supported by the DoD CHSSI (Department of Defense Common High Performance Computing Software Support Initiative), is designed to solve watershed problems on scalable computing systems. WASH123D is a first-principle, physics-based model to compute water flow and/or contaminant and sediment transport within a watershed system. In the WASH123D model, a watershed is conceptualized as a coupled system of one-dimensional (1-D) channel network, two-dimensional (2-D) overland regime, and three-dimensional (3-D) subsurface media. It aims to address the environmental issues concerning both water quantity and quality. To reach numerical solutions with reasonable and tolerable computer time for simulations that embrace large meshes, numerical algorithm improvement and code parallelization are two essential tasks. Mathematically, 1-D channel flow and 2-D overland flow are described with the St. Venant equations, which are solved with either the Semi-Lagrangian or the Eulerian finite element method. The 3-D subsurface flow is governed by the modified Richards equation, which is solved with the Eulerian finite element method. The contaminant transport and sediment transport equations, which are solved with the Lagrangian-Eulerian finite element method, are derived based on the mass conservation principle. A parallel in-element particle-tracking algorithm for unsteady flow is applied to backtrack fictitious particles from global nodes to determine the so-called Lagrangian values when the Semi-Lagrangian or the Lagrangian-Eulerian method is used. This paper addresses the parallelization of such a complex numerical model. In phase I, tasks including data structure and software design, software tool development, as well as tool integration are accomplished. The 2- and 3-D flow modules are expected to be employed in the production stage.

Numerical Strategies to Model Surface and Groundwater Interactions for the Biscayne Bay Coastal Wetlands Project Alternatives

WASH123D is a first-principle, physics-based numerical model that computes flow and transport in a watershed system that is conceptualized as a combination of 1-D channel network, 2-D overland regimes, and 3-D subsurface media. It has been selected as the tool to help evaluate the proposed alternatives of the Biscayne Bay Coastal Wetlands project that is one of the 40 projects included in the Florida Comprehensive Everglade Restoration Plan. In order to best rehydrate wetlands and reduce point source discharge to Biscayne Bay, rule-controlled coastal canal structures, rule-controlled pump stations, spreader swales, stormwater treatment areas, flowways, levees, culverts, roads, and backfilling canals are included in these project alternatives. Specified target freshwater flows for Biscayne Bay and the wetlands within the redistribution system are to be computed in each alternative, which will be used in performance measurement to determine the most adequate alternative for further investigation. In this paper, the numerical strategies to incorporate all the aforementioned hydrological features and processes included in the alternatives are presented. An example alternative will be used for demonstration.

Using the Parallel WASH123D Code to Simulate Overland-Subsurface Interactions

WASH123D model is a first-principle, physics-based model, where water flow and/or contaminant and sediment transport within a watershed system are computed. In the WASH123D model, a watershed is conceptualized as a coupled system of 1-D canal /stream network, 2-D overland regime, and 3-D subsurface media. It is designed to answer the environmental issues concerning both water quantity and quality. To reach numerical solutions with reasonable and tolerable computer time for a regional scale watershed simulation, numerical algorithm improvement and code parallelization are two essential tasks when a distributed numerical model, WASH123D, is used. This paper presents the code parallelization approach, followe d by demonstrating its scalability performance. The test problem of a large mesh uses the topographic data in the C -111 Spreader Canal Project. Background The C-111 Spreader Canal Project is one component of the more than 60 restoration plans under the Comprehensive Environmental Restoration Plan (CERP) and has a goal to provide water deliveries that will enhance the connection between the natural areas in the Southern Glades and Model Lands area of South Dade County. The spreader canal

Using WASH123D to Design Spreader Canals for Water Management in Watersheds

In the federally approved Comprehensive Everglades Restoration Plan (CERP, http://www.evergladesplan.org/), the restoration of the South Florida ecosystem is a major task for the U.S. Army Corps of Engineers and the South Florida Water Management District. Many competing entities have an interest in the restoration process which will probably include physical changes to the land surface and adjustments to water deliveries. The Biscayne Bay Coastal Wetlands (BBCW) Project is one component of more than 60 restoration plans and has a goal to restore or enhance freshwater wetlands, tidal wetlands, and near shore bay habitat. In an effort to restore wetlands, several structures, and management plans and scenarios are considered. One of the plans is to deliver fresh water from the existing canals through a shallow spreader canal system that is to distribute fresh water through wetlands into the Biscayne Bay. To achieve this, a tool is needed to design this complicated shallow spreader canal system. This paper presents how a spreader canal system, which includes 1D canal network routing, 2D overland flow, 3D subsurface flow, and flow through the interface of any two sub domains of the spreader canal system, is simulated with the WASH123D computer code. A brief description of the model will be given. A hypothetical example that includes two cases and uses topographic data and a high -resolution computational mesh for one project area will be considered to demonstrate how WASH123D can help design a spreader canal system. A couple of issues concerning run time and numerical convergence of the coupled flow model will also be discussed in this paper. Background The Biscayne Bay Coastal Wetlands (BBCW) Project is one component of the more than 60 restoration plans and has a goal to restore the coastal wetlands area in Central and South Biscayne Bay along its western shoreline. In the existing condition, fresh water

Modeling Surface and Subsurface Hydrologic Interactions in the Biscayne Bay Coastal Wetlands

Restoration of the South Florida ecosystem is a major undertaking for the U.S. Army Corps of Engineers and the South Florida Water Management District in the federally approved Comprehensive Environmental Restoration Plan (CERP). Many competing entities have an interest in the restoration process which will probably include physical changes to the land surface and adjustments to water deliveries. Restoration plans will be developed using numerical hydrologic models that optimize the benefits among the various stakeholders. Since the health of the South Florida ecosystem is impacted by subtle exchanges of surface and subsurface water, there is a need for great accuracy in the models that will identify the best restoration plans. The accuracy needed for design level simulations of the various projects approaches +/0.1 ft in the water elevations. Models with insufficient physical processes and/or highly schematized models with large grid or mesh induced error will not be sufficient to make design level simulation of South Florida hydrologic processes. Additionally, since a large number of models will be set up and analyzed, effective and efficient graphical user interfaces are needed for the models. This paper presents the results of testing WASH123D code to simulate surface water and subsurface interactions in the groundwater system.