Search NASA⌕ Search

SEARCH · Search NASA

Results for “WATER FLOW”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

Deep Learning for Intelligent Bubble Size Detection in the Spallation Neutron Source Visual Target

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) will undergo proton power upgrade (PPU), increasing the proton beam power from 1.4 MW to 2.8 MW. From 2.8 MW, 2.0 MW will go to the current First Target Station and the rest will go to the future Second Target Station (STS). The First Target Station uses a liquid mercury target that is contained in a 316L stainless steel vessel. The proton beam is pulsed at 60 Hz, with a pulse of about 0.7μs. When the proton beam hits the target, the intense energy deposition leads to a rapid rise in temperature in the mercury. This temperature rise creates pressure waves that propagate through the mercury and cause cavitation erosion. The power upgrade will cause stronger pressure waves that will further increase damage because of cavitation. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at 1.4 MW. However, at higher power, additional mitigation is necessary. Therefore, the 2 MW target vessel will be equipped with swirl bubblers and an additional gas injection port near the nose to inject more gas in the target. To develop a gas injection strategy and design, flow visualization in water with a transparent prototypical target (“visual target”) was performed. Bubble sizes and their spatial distribution in the flow loop are crucial to understanding the effectiveness of the bubbles in mitigating pressure waves. Bubbles were generated in the visual target under varied conditions of input pressures with helium and air. Images were captured using a high-speed camera at varied frame rates at different positions away from the swirl bubbler and different depths in the flow loop under varying lighting conditions. Initially, methods such as circular Hough transforms were applied after a series of images processing to obtain a general distribution of bubble sizes. Bubbles smaller than 500 μm are preferred to effectively mitigate the effect of pressure waves, which demands an accurate bubble detection and sizing system. Intelligent detection and identification of bubble sizes alleviate misdetection and improves accuracies. Employing neural networks, intelligent detection of bubble sizes and their distribution was developed and provides a robust alternative to traditional techniques. Human intervention was employed to label in-focus and out-of-focus bubbles in the set of training images. An object detection network using a pretrained convolutional neural network was created that extracted the features from the training images. Data augmentation was used to improve network accuracy through a random transformation of the original data.

Rasheed, Fayaz↗

Scoping Study of Integrated Geothermal CO 2 Heat Pump and Water Well System in Cold Climate Region

This study evaluates the feasibility of using ground-source heat pumps (GSHPs) integrated with water wells for residential heating demands in cold regions of the United States. Four different heat pump systems were evaluated and compared: a CO 2 air-source heat pump (ASHP), a R-410A ASHP, a CO 2 GSHP integrated with well water, and a R-410A GSHP integrated with well water. Simulations were conducted at both equipment and building integration levels, followed by a nationwide analysis across 10 selected cities. Results indicate that integrating GSHP technology with existing water wells (e.g., artesian, bedrock, drilled, and collector wells) is a feasible heating solution because of their high flow rates, good water quality, stable temperatures, and low environmental impact. Although R-410A systems demonstrate higher efficiency, CO 2 systems are more suitable for heating from contamination-free and energy security perspective. The CO 2 GSHPs maintain better performance in colder climates compared with that of ASHPs, providing higher coefficient of performance and environmental benefits because CO 2 refrigerant leakage would not contaminate well water. A national analysis confirms that CO 2 GSHPs outperform furnaces, especially in cold climates, making them a promising solution for residential heating.

99 GENERAL AND MISCELLANEOUS↗

Scoping Study of Integrated Geothermal CO 2 Heat Pump and Water Well System in Cold Climate Region (CRADA Final Report)

This study evaluates the feasibility of using ground-source heat pumps (GSHPs) integrated with water wells for residential heating demands in cold regions of the United States. Four different heat pump systems were evaluated and compared: a CO 2 air-source heat pump (ASHP), a R-410A ASHP, a CO 2 GSHP integrated with well water, and a R-410A GSHP integrated with well water. Simulations were conducted at both equipment and building integration levels, followed by a nationwide analysis across 10 selected cities. Results indicate that integrating GSHP technology with existing water wells (e.g., artesian, bedrock, drilled, and collector wells) is a feasible heating solution because of their high flow rates, good water quality, stable temperatures, and low environmental impact. Although R-410A systems demonstrate higher efficiency, CO 2 systems are more suitable for heating from contamination-free and energy security perspective. The CO 2 GSHPs maintain better performance in colder climates compared with that of ASHPs, providing higher coefficient of performance and environmental benefits because CO 2 refrigerant leakage would not contaminate well water. A national analysis confirms that CO 2 GSHPs outperform furnaces, especially in cold climates, making them a promising solution for residential heating.

15 GEOTHERMAL ENERGY↗

System and method for high efficiency electrochemical desalination

The present disclosure relates to a capacitive deionization (CDI) system for desalinating salt water. The system may have a capacitor formed by spaced apart first and second electrodes, which enable a fluid flow containing salt water to pass either between them or through them. An input electrical power source is configured to generate an electrical forcing signal between the two electrodes. The electrical forcing signal represents a periodic signal including at least one of voltage or current, and which can be represented as a Fourier series. One component of the Fourier series is a constant, and a second component of the Fourier series is a sinusoidal wave of non-zero frequency which has the highest amplitude of the additive components of the Fourier series. The amplitude of the sinusoidal wave component is between 0.85 and 1.25 times the amplitude of the periodic signal.

Hawks, Steven↗

Digitally manufactured air plasma-on-water reactor for nitrate production

The sustainable production of food to support the increasing world population is one of humanity’s most pressing challenges. Plasma activated water, produced using renewable energy, can help fulfill plants’ needs in sustainable agriculture approaches. The design, implementation, and characterization of a digitally manufactured air plasma-on-water reactor (POWR) for the synthesis of nitrate as green nitrogen fertilizer is presented. The interaction of air plasma-generated reactive oxygen and nitrogen species with water produces nitrate (NO 3 - ) and related species, which are the main nitrogen-containing nutrients for plants. The mild conditions of the operation of the POWR opens the possibility to use plastics, particularly through digital manufacturing strategies such as 3D-printing, for its fabrication. A pin-to-plate reactor configuration powered by high-voltage alternating power is chosen due to its simplicity and efficacy. A computational thermal-fluid model is used to evaluate the design and attain expected operational characteristics. The experimental characterization of the POWR encompassed design and operation parameters, namely electrode-water spacing, air flow rate, and voltage level. A machine learning approach is implemented to extract and quantify characteristic features of the plasma–water interaction, such plasma volume and plasma–water interface area. Experimental results revealed that the nitrate production rate varies linearly with dimensionless plasma volume. The design, fabrication, and characterization methods presented can be adapted to other POWRs and help enable on-demand nitrogen fertilizer production at low environmental and economic cost.

Physics↗

Numerical and Physical Simulation of the Fluid Flow in a Beam Blank Mold Fed Through the Center of the Web

Fluid flow inside a beam blank mold fed through a three-port SEN (two lateral ports and one bottom port), positioned at the center of the mold, has been investigated. Literature survey shows that this kind of configuration is not frequently used. It had been shown that it is possible to get a symmetrical flow with this configuration. The downward inclination of the ports should influence the transient fluid flow. The greater the downward inclination, the greater is the instability of the jet and of the interface between the immiscible fluids that simulate the slag/metal behavior. The flow characterization was made by dye dispersion, PIV technique, measurement of the meniscus oscillation using ultrasonic sensors, and CFD simulations. The slag/metal behavior was simulated using water and oils with different physical properties. The results from PIV as well as from observations of the water–oil interface have been used to validate the CFD simulations. Oils with density close to water resulted in more instability at the interface with entrainment starting from fluid flow rate of 125 L/min (equivalent to casting speed of 1 m/min). Decreasing the viscosity of slag (oils) and increasing the casting velocity (water flow rate) result in reduction of the interfacial stability.

36 MATERIALS SCIENCE↗

Financial Analysis of the High Flow Experiment conducted at the Glen Canyon Dam during Water Year 2023

The Glen Canyon Dam (GCD) is a Colorado River Storage Project (CRSP) power resource that is a component of the Salt Lake City Area Integrated Projects (SLCA/IP). The 2016 record of decision (ROD) for the GCD long-term experimental and management plan (LTEMP) final Environmental Impact Statement (EIS) specified criteria for GCD monthly water releases, daily and hourly operating limits, and experimental releases. This report examines the financial implications of the high flow experiment (HFE) conducted at GCD during the spring of Water Year (WY) 2023 as required by the LTEMP HFE Protocol. This report is part of a series of reports that describe the financial costs of LTEMP experimental releases since the 2016 ROD was adopted in January 2017. Previous reports analyzed the impact of several past HFEs and Bug Flow Experiments. This report focuses on the HFE conducted in April 2023. For this experimental release, financial costs of approximately $1.33 million were incurred because the HFE required sustained water releases exceeding the power plant’s maximum turbine flow rate. In addition, during the experiment, operators were not allowed to shape GCD power production, either to follow Firm Electric Service (FES) customer day-ahead energy deliveries or to respond to market prices. This study identifies the main factors contributing to the HFE costs and examines the interdependencies among these factors. It applies an integrated set of tools to estimate Western Area Power Administration (WAPA) financial impacts by simulating GCD under two types of cases; namely, (1) a “With Experiment” case that mimics the operations that actually occurred and (2) a “Without Experiment” case that simulates operations under the assumption that the HFE did not occur. The “With Experiment” case mimics operations during the HFE and the entire month the HFE occurred. It complies with LTEMP hourly and daily operating criteria. The “Without Experiment” case assumes that the HFE did not occur. The monthly water release volume is assumed to be identical under both cases. The Colorado River Storage Project Python-based model (CRiSPPy) model was the main modeling tool used to simulate the dispatch of the GCD hydropower plant and associated water releases from Lake Powell. In the modeling process, the research team used extensive data sets and historical information on SLCA/IP power plant characteristics, hydrologic conditions, and WAPA’s power purchases and sales prices. In addition to estimating the financial impact of the HFE, the team used the CRiSPPy model to gain insights into the interplay among ROD operating criteria, exceptions made to criteria to accommodate the HFE, and WAPA operating practices.

13 HYDRO ENERGY↗

Peatland hydrology across scale: a probabilistic framework for confronting variability, heterogeneity, and uncertainty (Final Report)

Our project aimed to guide the representations of peatland hydrology and water-carbon feedbacks within Earth System Models, by improving process understanding and developing parsimonious models that elucidate the effects of hydroclimatic variability and spatial heterogeneity. Our research questions centered on the three key drivers of peatland hydrology: (1) seasonal and interannual hydroclimatic fluctuations, (2) spatial heterogeneity of peatland microtopography, and (3) hydrological connectivity across landscape units within a peatland watershed. Our analyses and model testing used existing long-term datasets within the Marcell Experimental Forest (MEF), with new data collected whenever necessary. To date, we have demonstrated the strong influence of water table elevations on the temperature sensitivity of CH 4 emissions, using a newly developed eddy covariance dataset spanning eleven years at MEF. Specifically, higher water tables dampen the increase in CH 4 emissions in the spring as well as their decrease in the fall, resulting in hysteresis. These results imply that any hydroclimatological changes in peatlands that shift seasonal water availability from winter to summer will increase annual CH 4 emissions, even if soil temperature remains unchanged. To further investigate water table elevation change within the peatland microtopography, we installed automated water table gauges across four bog-forest transects (the “lagg”) in two peatland watersheds at MEF, which are hotspots of intense biogeochemical activity. These measurements will give us new information about the extent of lagg expansion and contraction during snowmelt and high intensity rainfall events, as well as the directionality of water and nutrient flow into the bog from the surrounding forests. Finally, we investigated the hydrological connectivity across the peatland watershed complex during early spring – a critical period where water table is elevated from snowmelt – using extended hydrological records from MEF (e.g., snow and frost depths, water table elevations, streamflow). The importance of thermal and wintertime processes within the Energy Exascale Earth System Model (E3SM) was determined through sensitivity analyses. Results from data analyses show that (i) streamflow has decreased over decades likely due to increased evapotranspiration rates (despite no detectable trends in precipitation), (ii) spring streamflow generation is controlled by a sequence of fill-and-spill pathways from the snowpack to the stream through peatland water stores, and (iii) frost depth is a key explanatory variable for the timing and magnitude of streamflow. These results suggest that frost plays an important role in connecting surface water storage in peatlands to stream outlets, and that the hydrological connectivity across the peatland watershed complex can mediate the sensitivity of hydrological responses to climate variations.

54 ENVIRONMENTAL SCIENCES↗

Conceptual and Mathematical Foundation for the FE/NETL CO 2 Prophet Model for Simulating CO 2 Enhanced Oil Recovery, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document describes the mathematical foundation for the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production and Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

User’s Manual for StrmtbGen, the Stream Tube Generating Part of the FE/NETL CO 2 Prophet Model, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document is a User’s Manual for StrmtbGen, which is one part of the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Export of Ice Sheet Meltwater from Upernavik Fjord, West Greenland

Meltwater from Greenland is an important freshwater source for the North Atlantic Ocean, released into the ocean at the head of fjords in the form of runoff, submarine melt, and icebergs. The meltwater release gives rise to complex in-fjord transformations that result in its dilution through mixing with other water masses. The transformed waters, which contain the meltwater, are exported from the fjords as a new water mass Glacially Modified Water (GMW). Here we use summer hydrographic data collected from 2013 to 2019 in Upernavik, a major glacial fjord in northwest Greenland, to describe the water masses that flow into the fjord from the shelf and the exported GMWs. Using an optimum multi-parameter technique across multiple years we then show that GMW is composed of 57.8% ± 8.1% Atlantic Water (AW), 41.0% ± 8.3% Polar Water (PW), 1.0% ± 0.1% subglacial discharge, and 0.2% ± 0.2% submarine meltwater. We show that the GMW fractional composition cannot be described by buoyant plume theory alone since it includes lateral mixing within the upper layers of the fjord not accounted for by buoyant plume dynamics. Consistent with its composition, we find that changes in GMW properties reflect changes in the AW and PW source waters. Using the obtained dilution ratios, this study suggests that the exchange across the fjord mouth during summer is on the order of 50 mSv (1 Sv ≡ 10 6 m 3 s –1 ) (compared to a freshwater input of 0.5 mSv). This study provides a first-order parameterization for the exchange at the mouth of glacial fjords for large-scale ocean models.

58 GEOSCIENCES↗

Managed aquifer recharge as a strategy to redistribute excess surface flow to baseflow in snowmelt hydrologic regimes

Water management in snowmelt hydrologic regimes, characterized by large annual fluctuations in stream flow driven by seasonal snow melt, faces the challenge of highly variable supply that often does not align with timing of demand. Climate change may exacerbate management challenges by significantly reducing snowpack or shifting snow melt earlier. Here, managed aquifer recharge (MAR) is evaluated as a potential strategy to reallocate excess early-season stream flow to time periods when less surface water is available. This strategy differs from traditional MAR, where the goal is to minimize loss to surface water. We assess how to site MAR operations such that groundwater recharge flows back to the surface water system in a lagged manor to benefit water management objectives, which we term “enhanced baseflow.” We use a regional groundwater model for the Treasure Valley aquifer located in southwestern Idaho, United States to demonstrate a generalizable approach using regional groundwater models as tools to identify favorable baseflow enhancement locations. Hypothetical MAR is simulated at 197 candidate locations, which are then evaluated for how effectively they meet potential management objectives. In addition to demonstrating the modeling and evaluation approach, we discuss lessons learned from applying a pre-existing regional groundwater model to MAR for enhanced baseflow and also describe important considerations, such as the physical and institutional availability of surface flows and specific management objectives, when assessing regional and site-specific suitability of MAR for enhanced baseflow as a potential management strategy.

54 ENVIRONMENTAL SCIENCES↗

Effect of PVD-coated chromium on the subcooled flow boiling performance of nuclear reactor cladding materials

Here we elucidate the separate effect of a thin Cr coating deposited by physical vapor deposition (PVD) on the subcooled flow boiling performance of zircaloy-4. First, we run flow boiling experiments on prototypical zircaloy-4 surfaces mimicking the scratch pattern and surface roughness of nuclear reactor claddings. Then, we PVD-coat a 0.3 µm thick chromium layer on the same exact surface and repeat the same flow boiling investigations. All experiments are run using deionized water at atmospheric pressure, flowing on a 1 × 3 cm 2 rectangular cross section channel at a rate of 1000 kg/m 2 /s and a subcooling of 10 K. We measure the average temperature of the boiling surface at increasing surface heat fluxes, covering a wide range of heat transfer regimes, from single-phase forced convection to the boiling crisis. We also record high-speed videos of the boiling process, which we postprocess to measure bubble nucleation site density, growth time, departure diameter and frequency. The surface analysis reveals that, while the chromium coating does not seem change the surface roughness and morphology, it improves surface wettability. However, it decreases the critical heat flux. The chromium coating causes an increase of nucleation temperature, bubble departure diameter and growth time, and a reduction of the nucleation site density. The concurrence of these observations indicates that a size reduction of the nucleation sites, conformally covered by the chromium coating, may be the cause of the boiling performance deterioration. We confirm this hypothesis repeating the same analysis on a FeCrAl sample prepared and tested using the same protocol as the zircaloy-4 sample, but with a different initial surface texture.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Single-phase, natural circulation annular flow measurements for cartridge loop irradiation experiments

The nuclear industry is increasingly considering cartridge-style experiments for irradiation testing of advanced reactor fuels and materials under flowing conditions. Cartridge loops do not require the extensive support infrastructure that are necessary for external flow loops and minimize the possibility of coolant solidification over the long distance from the reactor to the external facilities. However, there is a general lack of quality flow data for internally heated fluids in an annular configuration representative of a cartridge-type irradiation experiment, particularly one with natural circulation. To address this data need, a series of experiments was conducted to measure the natural circulation flow rates of pressurized water in a sealed, internally heated vessel with annular flow conditions that represent a molten salt or sodium cartridge loop. Temperatures and flow rates were measured under steady-state and transient conditions. Here, this paper describes the facility, methods, and results of the experiments, including the determination of nondimensional parameters. A simple 1D model of the natural convection flow rates agrees well with the experimental results. Applying this model to simulate a liquid salt cartridge experiment predicts that natural circulation flow might be able to provide liquid salt Reynolds numbers similar to those of some molten salt reactor concepts at relevant power densities.

42 ENGINEERING↗

Desalinating a Real Hyper-Saline Pre-Treated Produced Water via Direct-Heat Vacuum Membrane Distillation

Membrane distillation (MD) is an emerging thermal desalination technology capable of desalinating waters of any salinity. During typical MD processes, the saline feedwater is heated and acts as the thermal energy carrier; however, temperature polarization (as well as thermal energy loss) contributes to low distillate fluxes, low single-pass water recovery and poor thermal efficiency. An alternative approach is to integrate an extra thermal energy carrier as part of the membrane and/or module assembly, which can channel externally provided heat directly to the membrane-feedwater interface and/or along the feed channel length. This direct-heat delivery has been demonstrated to increase single-pass water recovery and enhance the overall thermal efficiency. We developed a bench-scale direct-heated vacuum MD (DHVMD) process to desalinate pre-treated oil and gas "produced water" with an initial total dissolved solids of 115,500 ppm at a feed temperature ranging between 24 and 32 degrees C. We evaluated both water flux and specific energy consumption (SEC) as a function of water recovery. The system achieved a 50% water recovery without significant scaling, with an average flux >6 kg m-2 hr-1 and a SEC as low as 2,530 kJ kg-1. The major species of mineral scales (i.e., NaCl, CaSO4, and SrSO4) that limited the water recovery to 68% were modeled in terms of thermodynamics and identified by scanning electron microscopy and energy-dispersive X-ray spectroscopy. In addition, we further developed and employed a physics-based process model to estimate temperature, salinity, water transport and energy flows for full-scale vacuum MD and DHVMD modules. Model results show that a direct-heat input rate of 3,600 W can increase single-pass water recovery from 2.1% to 3.1% while lowering the thermal SEC from 7,800 kJ kg-1 to 6,517 kJ kg-1 in an unoptimized module. Finally, the scaling up potential of DHVMD process is briefly discussed.

brine↗

Geochemical and microbial characteristics of seepage water and mineral precipitates in a radwaste disposal facility impacted by seawater intrusion and high alkalinity

The construction of an underground facility can dramatically change the quality, flow direction, and level of groundwater. It may also impact subsurface microbial composition and activity. Groundwater quality was monitored over eight years in two observational wells near an underground disposal facility on the east coast of South Korea. The results showed dramatic increases in dissolved ions such as O 2 , Na, Ca, Mg, and SO 4 during facility construction. Seepage water samples downgradient from the silos and tunnels, and precipitates deposited along the seepage water flow path were collected to determine the impact inside the disposal facility. X-ray analysis (powder X-ray diffraction (pXRD) and X-ray absorption fine structure (XAFS)) were used to characterize the mineral precipitates. Microbial community composition was determined by 16S rRNA gene sequencing. The seepage water composition was of two types: Ca-Cl and Ca-Na-HCO 3 . The ratio of Cl and δ 18 O showed that the Ca-Cl type seepage water was influenced by groundwater mixed with seawater ranging from 2.7% to 15.1%. Various sulfate-reducing bacteria were identified in the Ca-Cl type seepage water, exhibiting relatively high sulfate content from seawater intrusion. Samples from the Ca-Na-HCO 3 type seepage water had an extremely high pH (>10) and abundance of Hydrogenophaga. The precipitates observed along the flow path of the seepage water included calcite, ferrihydrite, green rust, and siderite, depending on seepage water chemistry and microbial activity. This study suggests that the construction of underground structures creates distinct, localized geochemical conditions (e.g., high alkalinity, high salinity, and oxic conditions), which may impact microbial communities. These biogeochemical changes may have undesirable large-scale impacts such as water pump clogging. Overall, an understanding of the process and long-term monitoring are essential to assess the safety of underground facilities.

54 ENVIRONMENTAL SCIENCES↗

Data used in manuscript "Quantifying Subsurface Flow and Solute Transport in a Snowmelt-recharged Hillslope with Multiyear Water Balance", submitted to Water Resources Research on 5-26-2022

Data in this package were used to calculate subsurface flow and solute (specifically dissolved nitrogen) transport in the Pumphouse Lower Montane (PLM) hillslope. These data are included in separate CSV files associated with each figure presented in the manuscript. The files describe site locations and elevations, weather data (air temperatures and precipitation), modeled evapotranspiration, time dependent groundwater elevations, calculated annual subsurface flow, hydraulic conductivity values, time trends in subsurface flow within different zones, pore water concentrations of total dissolve nitrogen (TDN), rates of TDN transport, soil water matric potentials following snowmelt, soil water contents during drainage, and evapotranspiration calculated by soil water mass balance versus by CLM4 model.

54 ENVIRONMENTAL SCIENCES↗

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗