SEARCH · Engineering Papers

Results for “surface complexation 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗

DOE OSTI2022

Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution

The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Improvement of Atmospheric Objective Analysis Over Sloping Terrain and Its Impact on Shallow-Cumulus Clouds in Large-Eddy Simulations

Surface topography strongly impacts the regional atmospheric circulation dynamically and thermodynamically. Over the Great Plains in the United States, the gently tilted slope is an important factor that impacts clouds, convection and regional circulations. This study enhances an atmospheric constrained variational analysis by using a terrain-following sigma vertical coordinate and applies to the data collected at the Southern Great Plain (SGP) site of the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program. Sensitivity studies are performed to examine the impact of the terrain effects on the derived large-scale atmospheric forcing fields and the simulated shallow-cumulus clouds in large-eddy simulation (LES) driven by the forcing. We found that the terrain impacts on the large-scale forcing fields are mainly at lower levels and are strongly controlled by up/downslope winds. The response of the derived forcing to the slope of the terrain is monotonic, but the response of the simulated shallow-cumulus clouds is more complex and depends on several factors. Overall, the terrain impact is small over SGP due to the small slope angle. However, the flat-surface assumption may cause larger biases in the large-scale forcing fields at locations with steeper terrain. The new sigma-coordinate algorithm, with its consideration of surface slope, should be more suitable to derive large-scale objective analysis over regions with steep terrain for application to force single-column models, cloud resolving models and LES models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗

DOE OSTI2021

Summertime low clouds mediate the impact of the large-scale circulation on Arctic sea ice

The rapid Arctic sea ice retreat in the early 21 st century is believed to be driven by several dynamic and thermodynamic feedbacks, such as ice-albedo feedback and water vapor feedback. However, the role of clouds in these feedbacks remains unclear since the causality between clouds and these processes is complex. Here, we use NASA CERES satellite products and NCAR CESM model simulations to suggest that summertime low clouds have played an important role in driving sea ice melt by amplifying the adiabatic warming induced by a stronger anticyclonic circulation aloft. The upper-level high pressure regulates low clouds through stronger downward motion and increasing lower troposphere relative humidity. The increased low clouds favor more sea ice melt via emitting stronger longwave radiation. Then decreased surface albedo triggers a positive ice-albedo feedback, which further enhances sea ice melt. Considering the importance of summertime low clouds, accurate simulation of this process is a prerequisite for climate models to produce reliable future projections of Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗

DOE OSTI2020

WHONDRS-GUI: a web application for global survey of surface water metabolites

Background The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) is a consortium that aims to understand complex hydrologic, biogeochemical, and microbial connections within river corridors experiencing perturbations such as dam operations, floods, and droughts. For one ongoing WHONDRS sampling campaign, surface water metabolite and microbiome samples are collected through a global survey to generate knowledge across diverse river corridors. Metabolomics analysis and a suite of geochemical analyses have been performed for collected samples through the Environmental Molecular Sciences Laboratory (EMSL). The obtained knowledge and data package inform mechanistic and data-driven models to enhance predictions of outcomes of hydrologic perturbations and watershed function, one of the most critical components in model-data integration. To support efforts of the multi-domain integration and make the ever-growing data package more accessible for researchers across the world, a Shiny/R Graphical User Interface (GUI) called WHONDRS-GUI was created. Results The web application can be run on any modern web browser without any programming or operational system requirements, thus providing an open, well-structured, discoverable dataset for WHONDRS. Together with a context-aware dynamic user interface, the WHONDRS-GUI has functionality for searching, compiling, integrating, visualizing and exporting different data types that can easily be used by the community. The web application and data package are available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1484811 , which enables users to simultaneously obtain access to the data and code and to subsequently run the web app locally. The WHONDRS-GUI is also available for online use at Shiny Server ( https://xmlin.shinyapps.io/whondrs/ ).

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1995

Studies of Two-Phase Gas-Liquid Flow in Microgravity

Two-phase gas-liquid flows are expected to occur in many future space operations. Due to a lack of buoyancy in the microgravity environment, two-phase flows are known to behave differently than those in earth gravity. Despite these concerns, little research has been conducted on microgravity two-phase flow and the current understanding is poor. This dissertation describes an experimental and modeling study of the characteristics of two-phase flows in microgravity. An experiment was operated onboard NASA aircraft capable of producing short periods of microgravity. In addition to high speed photographs of the flows, electronic measurements of void fraction, liquid film thickness, bubble and wave velocity, pressure drop and wall shear stress were made for a wide range of liquid and gas flow rates. The effects of liquid viscosity, surface tension and tube diameter on the behavior of these flows were also assessed. From the data collected, maps showing the occurrence of various flow patterns as a function of gas and liquid flow rates were constructed. Earth gravity two-phase flow models were compared to the results of the microgravity experiments and in some cases modified. Models were developed to predict the transitions on the flow pattern maps. Three flow patterns, bubble, slug and annular flow, were observed in microgravity. These patterns were found to occur in distinct regions of the gas-liquid flow rate parameter space. The effect of liquid viscosity, surface tension and tube diameter on the location of the boundaries of these regions was small. Void fraction and Weber number transition criteria both produced reasonable transition models. Void fraction and bubble velocity for bubble and slug flows were found to be well described by the Drift-Flux model used to describe such flows in earth gravity. Pressure drop modeling by the homogeneous flow model was inconclusive for bubble and slug flows. Annular flows were found to be complex systems of ring-like waves and a substrate film. Pressure drop was best fitted with the Lockhart- Martinelli model. Force balances suggest that droplet entrainment may be a large component of the total pressure drop.

Bousman, William Scott↗

NASA NTRS2011

GIS Data Based Automatic High-Fidelity 3D Road Network Modeling

3D road models are widely used in many computer applications such as racing games and driving simulations_ However, almost all high-fidelity 3D road models were generated manually by professional artists at the expense of intensive labor. There are very few existing methods for automatically generating 3D high-fidelity road networks, especially those existing in the real world. This paper presents a novel approach thai can automatically produce 3D high-fidelity road network models from real 2D road GIS data that mainly contain road. centerline in formation. The proposed method first builds parametric representations of the road centerlines through segmentation and fitting . A basic set of civil engineering rules (e.g., cross slope, superelevation, grade) for road design are then selected in order to generate realistic road surfaces in compliance with these rules. While the proposed method applies to any types of roads, this paper mainly addresses automatic generation of complex traffic interchanges and intersections which are the most sophisticated elements in the road networks

Wang, Jie↗

NASA NTRS2014

NMME Monthly / Seasonal Forecasts for NASA SERVIR Applications Science

This work details use of the North American Multi-Model Ensemble (NMME) experimental forecasts as drivers for Decision Support Systems (DSSs) in the NASA / USAID initiative, SERVIR (a Spanish acronym meaning "to serve"). SERVIR integrates satellite observations, ground-based data and forecast models to monitor and forecast environmental changes and to improve response to natural disasters. Through the use of DSSs whose "front ends" are physically based models, the SERVIR activity provides a natural testbed to determine the extent to which NMME monthly to seasonal projections enable scientists, educators, project managers and policy implementers in developing countries to better use probabilistic outlooks of seasonal hydrologic anomalies in assessing agricultural / food security impacts, water availability, and risk to societal infrastructure. The multi-model NMME framework provides a "best practices" approach to probabilistic forecasting. The NMME forecasts are generated at resolution more coarse than that required to support DSS models; downscaling in both space and time is necessary. The methodology adopted here applied model output statistics where we use NMME ensemble monthly projections of sea-surface temperature (SST) and precipitation from 30 years of hindcasts with observations of precipitation and temperature for target regions. Since raw model forecasts are well-known to have structural biases, a cross-validated multivariate regression methodology (CCA) is used to link the model projected states as predictors to the predictands of the target region. The target regions include a number of basins in East and South Africa as well as the Ganges / Baramaputra / Meghna basin complex. The MOS approach used address spatial downscaling. Temporal disaggregation of monthly seasonal forecasts is achieved through use of a tercile bootstrapping approach. We interpret the results of these studies, the levels of skill by several metrics, and key uncertainties.

Robertson, Franklin R.↗

NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

Mesh motion is the process by which a computational domain is updated in time to reflect physical changes in the material the domain represents. Such a technique is needed in the study of the thermal response of ablative materials, which erode when strong heating is applied to the boundary. Traditionally, the thermal solver is coupled with a linear elastic or biharmonic system whose sole purpose is to update mesh node locations in response to altering boundary heating. Simple mesh motion algorithms rely on boundary surface normals. In such schemes, evolution in time will eventually cause the mesh to intersect and "tangle" with itself, causing failure. Furthermore, such schemes are greatly limited in the problems geometries on which they will be successful. This paper presents a comprehensive and sophisticated scheme that tailors the directions of motion based on context. By choosing directions for each node smartly, the inevitable tangle can be completely avoided and mesh motion on complex geometries can be modeled accurately.

Droba, Justin↗

DOE OSTI2023

Systematic modification of functionality in disordered elastic networks through free energy surface tailoring

A combined machine learning–physics–based approach is explored for molecular and materials engineering. Specifically, collective variables, akin to those used in enhanced sampled simulations, are constructed using a machine learning model trained on data gathered from a single system. Through the constructed collective variables, it becomes possible to identify critical molecular interactions in the considered system, the modulation of which enables a systematic tailoring of the system’s free energy landscape. To explore the efficacy of the proposed approach, we use it to engineer allosteric regulation and uniaxial strain fluctuations in a complex disordered elastic network. Its successful application in these two cases provides insights regarding how functionality is governed in systems characterized by extensive connectivity and points to its potential for design of complex molecular systems.

42 ENGINEERING↗

DOE OSTI2020

Understanding asymmetries using integrated simulations of capsule implosions in low gas-fill hohlraums at the National Ignition Facility

Current capsule implosions at the National Ignition Facility (NIF) using high-density-carbon ablators and laser energies close to 2 MJ have shown neutron yields in excess of 50 kJ. Improving on this performance requires understanding of the different degradation mechanisms. For many NIF implosions, nuclear diagnostic signatures have inferred residual hot-spot velocities that correlate with fuel areal density variations consistent with a low-order mode-1 asymmetry. A current working hypothesis attributes these asymmetries to a combination of beam to beam variations in the laser delivery and possibly coupling to target features, such as the diagnostic holes needed for x-ray imaging. Recently, a new source has been identified, thickness variation in the ablator shell. To gain better understanding and eventually mitigate the causes of <ρR> asymmetries, 3D integrated simulations using the actual delivered laser powers are needed. To capture the effect of the diagnostic holes (DHs) using direct numerical simulation would require significantly large computational resources. Instead, our 3D simulations make use of a subgrid model developed using highly-resolved 2D simulations that include several details of the DH engineering complexity. Simulations of NIF shots using hohlraums without DHs, to isolate the effect of beam-to-beam variations, reproduce fairly well the observed nuclear diagnostic signatures. Similarly, reasonable agreement between data and simulations is also obtained in the presence of diagnostic holes. To account for the remaining discrepancies a sensitivity study of ablator thickness variation showed that 1% thickness asymmetries are comparable in effect to 1% peak drive mode-1 asymmetries. Furthermore, this study identified sensitivity to variations in the imbedded doped layer (needed to shield the DT ice from the high energy x-rays generated in the hohlraum) thickness even when the inner and outer surface of the ablator are perfectly spherical.

radiation hydrodynamics↗

NASA NTRSDate not supplied

How Well Satellite Remote Sensing Can Inform Surface Level Ozone Production Sensitivity to Precursor Trace Gas Emissions?

Surface-level ozone (O3) pollution mitigation strategies rely on a clear understanding of the atmospheric chemical processes involving O3, NOx (nitric oxide (NO) + nitrogen dioxide (NO2)) and volatile organic compounds (VOCs). A robust spatiotemporal classification of O3 photochemical regimes with relative chemical sensitivity of local O3 formation to emission reductions of NOx (NOx-limited regime) versus VOCs (radical-limited regime) is required to design anthropogenic emission reduction policies in order to improve surface air quality. Numerous studies found that the ratio of satellite-retrieved Vertical Column Densities (VCDs) of formaldehyde (HCHO) to NO2 can diagnose O3 production sensitivity. However, uncertainty still remains about the accuracy of such estimates using satellite-based retrievals. This study conducts an investigation to identify how well satellite data can inform surface-level O3 sensitivity by validating satellite-based estimates against aircraft-based retrievals taken during the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) and Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaigns. For this study, we use HCHO and NO2 retrievals from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and high spatiotemporal resolution TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite. The Community Multiscale Air Quality (CMAQ) modelling system is used to provide trace gas vertical profiles to recalculate air mass factor for deriving consistent satellite-based VCDs. To demonstrate the ability of satellites to diagnose O3 production regimes, we validate satellite-derived O3 sensitivity regimes, along with the OMI and TROPOMI products applied, using HCHO and NO2 retrievals from the aircraft-based instruments: Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) spectrometer and Geostationary Coastal and Air Pollution Events Airborne Simulator (GCAS) instruments deployed during the LISTOS-2018 and OWLETS-2 campaigns. The results of this work will advance our fundamental knowledge of the spatiotemporal evolution of the complex O3-NOx-VOC chemistry, and will advance our understanding about the quality of remote-sensing, suborbital, and satellite data for accurately observing surface-level O3 production sensitivity.

Satellite↗

DOE OSTI2024

Influence of porous aluminosilicate grain size materials in experimental and modelling Cs + adsorption kinetics and wastewater column process

This paper focuses on the influence of the grain size of a geopolymer based adsorbent on its Cs + adsorption performances both in batch and fixed-bed process. The geopolymer phase was used as a binder to support NaY zeolite particle in a 20 wt% charged porous composite with 160 m 2 .g –1 of porous surface area. These samples were shaped with three grain sizes (50 /100/500 µm) to remove 80–90 mg/g of Cs + in batch and column operations. After their microstructural and porous characterizations, their efficiency and adsorption characteristics were investigated through adsorption isotherms and kinetic in the two processes. While the grain size has no influence on the maximal extraction capacity of the adsorbent, it strongly affects the sorption kinetic. By coupling experimental data and a modelling approach, the complex sorption mechanism was highlighted, suggesting a new insight of the contaminant sorption kinetic. Then, comparison of batch and column adsorption experiments illustrates the detailed explanation of various process parameters for column study. The results show challenges for fixed-bed column utilization by the choice of the appropriate grain size as a compromise between the material sorption kinetic and hydrodynamic considerations. Furthermore, this is of high importance to more accurately optimize the design of column adsorption to assess the transport of Cs+ in multi-porous tailored grain size materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Hourly PM 2.5 Estimates across California from 2018 to 2023

This study presents a new data set of hourly PM 2.5 concentrations across California from 2018 to 2023 at a three-kilometer resolution. This data set was developed by assimilating observations from PurpleAir and the U.S. EPA Air Quality System monitors into wildfire smoke forecasts from the High-Resolution Rapid Refresh Smoke (HRRR-Smoke) model using the Gridpoint Statistical Interpolation (GSI) three-dimensional variational data assimilation framework. Archived forecasts of modeled wildfire smoke PM 2.5 from HRRR-Smoke create the background field for assimilation, which is then corrected using surface observations of total PM 2.5 . The resulting reanalysis from GSI provides an estimate of total PM 2.5 that minimizes error from both the observational and the model data. Validation results indicate strong performance, with monthly R 2 values ranging from 0.73 to 0.91 across the six-year data set, comparable to other PM 2.5 data sets. Case studies are presented for three major fire events, the 2018 Camp Fire, 2019 Kincade Fire, and 2020 Lightning Complex Fires to demonstrate the data set’s fidelity in resolving plume dynamics and local exposure patterns. Root-mean-squared error averaged over each month scales with average PM 2.5 concentrations, resulting in a low error under typical conditions but higher absolute errors during extreme smoke events. This is the first long-term, hourly PM 2.5 data set of its kind for California and enables the generation of subdaily exposure metrics, such as peak hourly concentrations, exceedance durations, and time-of-day exposure peaks. The novelty and strong validation of this data set make it a compelling resource for future studies on the impact and significance of subdaily PM 2.5 exposure.

PM2.5↗

DOE OSTI2022

The influence of recent and future climate change on spring Arctic cyclones

In recent decades, the Arctic has experienced rapid atmospheric warming and sea ice loss, with an ice-free Arctic projected by the end of this century. Cyclones are synoptic weather events that transport heat and moisture into the Arctic, and have complex impacts on sea ice, and the local and global climate. However, the effect of a changing climate on Arctic cyclone behavior remains poorly understood. This study uses high resolution (4 km), regional modeling techniques and downscaled global climate reconstructions and projections to examine how recent and future climatic changes alter cyclone behavior. Results suggest that recent climate change has not yet had an appreciable effect on Arctic cyclone characteristics. However, future sea ice loss and increasing surface temperatures drive large increases in the near-surface temperature gradient, sensible and latent heat fluxes, and convection during cyclones. The future climate can alter cyclone trajectories and increase and prolong intensity with greatly augmented wind speeds, temperatures, and precipitation. Such changes in cyclone characteristics could exacerbate sea ice loss and Arctic warming through positive feedbacks. The increasing extreme nature of these weather events has implications for local ecosystems, communities, and socio-economic activities.

54 ENVIRONMENTAL SCIENCES↗