Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Water model”

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 217 records · Page 12

A simple hydrologically based model of land surface water and energy fluxes for general circulation models

A generalization of the single soil layer variable infiltration capacity (VIC) land surface hydrological model previously implemented in the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model (GCM) is described. The new model is comprised of a two-layer characterization of the soil column, and uses an aerodynamic representation of the latent and sensible heat fluxes at the land surface. The infiltration algorithm for the upper layer is essentially the same as for the single layer VIC model, while the lower layer drainage formulation is of the form previously implemented in the Max-Planck-Institut GCM. The model partitions the area of interest (e.g., grid cell) into multiple land surface cover types; for each land cover type the fraction of roots in the upper and lower zone is specified. Evapotranspiration consists of three components: canopy evaporation, evaporation from bare soils, and transpiration, which is represented using a canopy and architectural resistance formulation. Once the latent heat flux has been computed, the surface energy balance is iterated to solve for the land surface temperature at each time step. The model was tested using long-term hydrologic and climatological data for Kings Creek, Kansas to estimate and validate the hydrological parameters, and surface flux data from three First International Satellite Land Surface Climatology Project Field Experiment (FIFE) intensive field campaigns in the summer-fall of 1987 to validate the surface energy fluxes.

Liang, XU↗

Assimilation of Remotely Sensed Leaf Area Index into the Noah-MP Land Surface Model: Impacts on Water and Carbon Fluxes and States over the Continental U.S.

Accurate representation of vegetation states is required for the modeling of terrestrial water-energy-carbon exchanges and the characterization of the impacts of natural and anthropogenic vegetation changes on the land surface. This study presents a comprehensive evaluation of the impact of assimilating remote sensing-based Leaf Area Index (LAI) retrievals over the Continental U.S. in the Noah-MP land surface model, during a time period of 2000 to 2017. The results demonstrate that the assimilation has a beneficial impact on the simulation of key water budget terms such as soil moisture, evapotranspiration, snow depth, terrestrial water storage and streamflow, when compared with a large suite of reference datasets. In addition, the assimilation of LAI is also found to improve the carbon fluxes of Gross Primary Production (GPP) and Net Ecosystem Exchange (NEE). Most prominent improvements in the water and carbon variables are observed over the agricultural areas of the U.S., where assimilation improves the representation of vegetation seasonality impacted by cropping schedules. The systematic, added improvements from assimilation in a configuration that employs high quality boundary conditions highlight the significant utility of LAI data assimilation in capturing the impacts of vegetation changes.

anthropogenic vegetation↗

Modeling of Dense Water Production and Salt Transport from Alaskan Coastal Polynyas

The main significance of this paper is that a realistic, three-dimensional, high-resolution primitive equation model has been developed to study the effects of dense water formation in Arctic coastal polynyas. The model includes realistic ambient stratification, realistic bottom topography, and is forced by time-variant surface heat flux, surface salt flux, and time-dependent coastal flow. The salt and heat fluxes, and the surface ice drift, are derived from satellite observations (SSM/I and NSCAT sensors). The model is used to study the stratification, salt transport, and circulation in the vicinity of Barrow Canyon during the 1996/97 winter season. The coastal flow (Alaska coastal current), which is an extension of the Bering Sea throughflow, is formulated in the model using the wind-transport regression. The results show that for the 1996/97 winter the northeastward coastal current exports 13% to 26% of the salt produced by coastal polynyas upstream of Barrow Canyon in 20 to 30 days. The salt export occurs more rapidly during less persistent polynyas. The inclusion of ice-water stress in the model makes the coastal current slightly weaker and much wider due to the combined effects of surface drag and offshore Ekman transport.

Signorini, Sergio R.↗

The desorptivity model of bulk soil-water evaporation

Available models of bulk evaporation from a bare-surfaced soil are difficult to apply to field conditions where evaporation is complicated by two main factors: rate-limiting climatic conditions and redistribution of soil moisture following infiltration. Both factors are included in the "desorptivity model', wherein the evaporation rate during the second stage (the soil-limiting stage) of evaporation is related to the desorptivity parameter, A. Analytical approximations for A are presented. The approximations are independent of the surface soil moisture. However, calculations using the approximations indicate that both soil texture and soil moisture content at depth significantly affect A. Because the moisture content at depth decreases in time during redistribution, it follows that the A parameter also changes with time. Consequently, a method to calculate a representative value of A was developed. When applied to field data, the desorptivity model estimated cumulative evaporation well. The model is easy to calculate, but its usefulness is limited because it requires an independent estimate of the time of transition between the first and second stages of evaporation. The model shows that bulk evaporation after the transition to the second stage is largely independent of climatic conditions.

Clapp, R. B.↗

Numerical modeling of runback water on ice protected aircraft surfaces

A numerical simulation for 'running wet' aircraft anti-icing systems is developed. The model includes breakup of the water film, which exists in regions of direct impingement, into individual rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted in the numerical solution procedure. The solid wall is modeled as a multilayer structure and the anti-icing system used is of the thermal type utilizing hot air and/or electrical heating elements embedded with the layers. Details of the calculation procedure and the methods used are presented.

Al-Khalil, Kamel M.↗

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water↗

Biophysical model of eelgrass and water quality in Coos Bay, OR shows greater mitigation potential for ocean acidification than hypoxia

Seagrass beds provide important ecosystem services and are valued, in part, for their potential to mediate stressors such as ocean acidification and hypoxia (OAH) for sensitive species. However, the susceptibility of seagrasses to anthropogenic impacts and recent declines motivate the need to better understand the drivers of seagrass and the water quality consequences that occur with variation in seagrass abundance. To meet this need, we leveraged existing monitoring data (water quality and seagrass), hydrodynamic circulation model, and biogeochemical model framework with seagrass submodel, to produce a biophysical model of Coos Bay estuary, Oregon, U.S. The model includes biogeochemical processes involving water quality, plankton, seagrass, and sediment-water interactions. Ecosystem models like this are useful for evaluating complex estuarine systems because they allow us to extend our understanding of system dynamics beyond existing observations and perform experiments to identify the processes driving observed patterns. We used the biophysical model of Coos Bay to evaluate the dynamics of water quality and native eelgrass (Zostera marina) under three eelgrass abundance scenarios (zero eelgrass, current extent, and maximum observed extent) to elucidate the relationship between eelgrass and OAH. Including eelgrass in the Coos Bay model produced results that more closely resembled water quality observations - dissolved oxygen (DO) and pH were more dynamic in simulations with eelgrass, often having both higher highs and lower lows. While there were some areas of the estuary where DO improved with the addition of eelgrass to the model there was overall a small net increase in harmful DO conditions (based on a salmon physiological threshold). In contrast, ocean acidification conditions, pH and calcium carbonate saturation state for aragonite (Ω), were improved (based on oyster requirements) with the addition of eelgrass - although the magnitude of improvement differed seasonally and spatially. Our new model represents a useful tool - one which accounts for and controls the relevant physical and biogeochemical processes - to evaluate conditions that confer resilience or enhance vulnerability to OAH in an important Pacific Northwest coastal estuary and results can inform the OAH-related dynamics occurring in other eastern boundary current estuaries.

FVCOM-ICM↗

Molecular dynamics study of ion capture from water by a model ionophore, tetraprotonated cryptand SC24

The molecular dynamics of chloride capture from water by the tetraprotonated cryptand SC24 has been studied for the cases of 19 distances between the criptand and the chloride. The chloride capture is found to be characterized by a rapid cooperative change in the conformation of the cryptand when the Cl(-) begins to enter the ligand and just as it encounters the energy barrier. The conformational transition is associated with a shift of three N-H bonds from the pure endo orientation, such that they point toward the chloride.

Owenson, Brian↗

Dynamic Model Investigation of Water Pressures and Accelerations Encountered During Landings of the Apollo Spacecraft

An experimental investigation was made to determine impact water pressures, accelerations, and landing dynamics of a 1/4-scale dynamic model of the command module of the Apollo spacecraft. A scaled-stiffness aft heat shield was used on the model to simulate the structural deflections of the full-scale heat shield. Tests were made on water to obtain impact pressure data at a simulated parachute letdown (vertical) velocity component of approximately 30 ft/sec (9.1 m/sec) full scale. Additional tests were made on water, sand, and hard clay-gravel landing surfaces at simulated vertical velocity components of 23 ft/sec (7.0 m/sec) full scale. Horizontal velocity components investigated ranged from 0 to 50 ft/sec (15 m/sec) full scale and the pitch attitudes ranged from -40 degrees to 29 degrees. Roll attitudes were O degrees, 90 degrees, and 180 degrees, and the yaw attitude was 0 degrees.

Tests↗

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗

The Liquid-Water Oscillation in Modeling Boundary-Layer Cumuli with Third-Order Turbulence Closure Models

A hierarchy of third-order turbulence closure models are used to simulate boundary-layer cumuli in this study. An unrealistically strong liquid-water oscillation (LWO) is found in the fully prognostic model, which predicts all third moments. The LWO propagates from cloud base to cloud top with a speed of 1 m/s. The period of the oscillation is about 1000 s. Liquid-water buoyancy terms in the third-moment equations contribute to the LWO. The LWO mainly affects the vertical profiles of cloud fraction, mean liquid-water mixing ratio and the fluxes of liquid-water potential temperature and total water, but has less impact on the vertical profiles of other second-moments and third-moments. In order to minimize the LWO, a moderate large diffusion coefficient and a large turbulent dissipation at its originating level are needed. However, this approach distorts the vertical distributions of cloud fraction and liquid-water mixing ratio. A better approach is to parameterize liquid-water buoyancy more reasonably. A minimally prognostic model, which diagnoses all third moments except for vertical velocity, is shown to produce better results, compared to a fully prognostic model.

Cheng, A.↗

Development of an improved model for runback water on aircraft surfaces

A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.

Al-Khalil, Kamel M.↗

Accounting for Static and Dynamic Open Water in the Modeling of SMAP Brightness Temperatures over Peatlands

Hydrological change in peatlands due to anthropogenic disturbance and global warming can release enormous amounts of greenhouse gas emissions. Passive microwave satellite observations are an opportunity to globally monitor these changes. Abundant static and dynamic open water surfaces in peatlands strongly affect observed brightness temperatures (Tb). Here, we account for these contributions in radiative transfer modeling using NASA's Goddard Earth Observing System Model version 5 (GEOS-5) static open water mask and, for the dynamic open water fraction, the simulated inundated area using a version of the GEOS-5 Catchment land surface model that has been modified for peatland areas (PEAT-CLSM). Modeled Tb is compared against two years of SMAP L-band Tb. Preliminary results indicate: (i) a bias reduction when including the static open water fraction in a simple RTM mixing model, and ii) significantly improved correlation between modeled and observed Tb when using land surface output from PEAT-CLSM instead of the operational CLSM.

passive microwave↗

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA'S ISRU Technology

Introduction: In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s) range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant. System Model: In FY24, the System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team developed a lunar water processing system model using the Mission Analysis and Integration Tool (MAIT) to estimate the M/V/P for ISRU subsystems operating under a wide range of Hydrogen (H2) and Oxygen (O2) production targets for the Space Technology Mission Directorate (STMD) [1]. Based on Japan Aerospace Exploration Agency’s (JAXA) surface operational requirements, this system architecture was modified to include the ability to excavate consolidated icy regolith (versus granular ice excavation using Kennedy Space Center’s (KSC) ISRU Pilot Excavator, IPEx) and explore the feasibility of processing the lunar water both inside and outside of the PSR. For the consolidated icy regolith case study, excavation was performed via a mobility transport chassis outfitted with The Regolith Ice Drill for Exploring New Terrain (TRIDENT) for drilling [2] and the Cold Operable Lunar Deployable Arm (COLDArm) [3] for regolith transfer. The system model determines the required rover and payload. M/V/P to handle the required regolith processing rates. The regolith is then sorted and heated to sublimate the ice (via an auger dryer). The exiting high temperature, low pressure vapor is cleaned of volatiles (via cold trap) and electrolyzed to produce H2 and O2. These products are then dried, liquified with 20 K and 90 K cryocoolers (for H2 and O2, respectively), and stored in cylindrical tanks. Study Goals: Due to the different ConOps options of regolith transport to the ridge for processing versus processing it directly inside the PSR, as well as the unknown regolith/water-ice composition, new excavation techniques and their power configurations are being evaluated within a ISRU system architecture for production targets less than NASA’s pilot plant (1 mT). This analysis investigates the feasibility of numerous excavation techniques, power architectures, and logistical operations and determines an optimal system configuration with regards to M/V/P. It aims to investigate which parameters, both locally and globally, have the greatest effect on each subsystem within the plant. This can be used to identify the most critical components of the plant, and guide future decisions on allocating funding for research and development. The results from this study may provide subsystem developers with appropriate interfaces with excavation subsystems and downstream processes, and assessing the overall feasibility of each excavation technique, power architecture, and logistical timeframe. References: [1] Carlson, A. et. al. (2024) ICES. [2] Zacny, K., et. al. (2024) “ASCE Earth and Space”. [3] McCormick, R., et. Al. (2024) IEEE Xplore.

ISRU↗

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗