Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures

Performance improvements in biomolecular simulations based on molecular dynamics (MD) codes are widely desired. Unfortunately, the factors, which allowed past performance improvements, particularly the microprocessor clock frequencies, are no longer increasing. Hence, novel software and hardware solutions are being explored for accelerating performance of widely used MD codes. In this paper, we describe our efforts on porting, optimizing and tuning of Large-scale Atomic/Molecular Massively Parallel Simulator, a popular MD framework, on heterogeneous architectures: multi-core processors with graphical processing unit (GPU) accelerators. Our implementation is based on accelerating the most computationally expensive non-bonded interaction terms on the GPUs and overlapping the computation on the CPU and GPUs. This functionality is built on top of message passing interface that allows multi-level parallelism to be extracted even at the workstation level with the multi-core CPUs and allows extension of the implementation on GPU-enabled clusters. We hypothesize that the optimal benefit of heterogeneous architectures for applications will come by utilizing all possible resources (for example, CPU-cores and GPU devices on GPU-enabled clusters). Benchmarks for a range of biomolecular system sizes are provided, and an analysis is performed on four generations of NVIDIA's GPU devices. On GPU-enabled Linux clusters, by overlapping and pipelining computation and communication, we observe up to 10-folds application acceleration in multi-core and multi-GPU environments illustrating significant performance improvements. Detailed analysis of the implementation is presented that allows identification of bottlenecks in algorithm, indicating that code optimization and improvements on GPUs could allow microsecond scale simulation throughput on workstations and inexpensive GPU clusters, putting widely desired biologically relevant simulation time-scales within reach of a large user community. In order to systematically optimize simulation throughput and to enable performance prediction, we have developed a parameterized performance model that will allow developers and users to explore the performance potential of future heterogeneous systems for biological simulations.

Agarwal, Pratul K.↗
DOE OSTI2021

Precipitation and nitrogen application stimulate soil nitrous oxide emission

Precipitation and nitrogen (N) fertilization are the two most important drivers for soil nitrous oxide (N 2 O) emission. However, the effects of changes in N fertilization and precipitation patterns (i.e., precipitation intensity and frequency) on N 2 O emissions in agricultural fields are still unclear. In this study, we simulated soil N 2 O emission under different precipitation patterns (6 precipitation intensities, and 12 precipitation frequencies by either merging or splitting precipitation events) and N fertilization rates (low, typical, and high N fertilization) in a cornfield using the DeNitrification-DeComposition model. The model was parameterized and validated using meteorological data and N experimental measurements in Nashville, Tennessee, USA. Results showed that soil water filled pore space (WFPS) and simulated soil N 2 O emission increased as precipitation intensity increased. Less frequent but high intensity precipitation treatments reduced the soil WFPS by 25.2% and stimulated soil N 2 O emission by 45.3%, while more frequent but low intensity precipitation treatments increased soil WFPS by 9.0% and reduced soil N 2 O emission by 23.9%. Compared to typical N fertilization, the sensitivity of soil N 2 O emission to precipitation was higher under high N than low N fertilization treatments, and the response ratios were 50.0% and 40.1%, respectively. There was significant interactive effect of precipitation intensity and N fertilization on soil N 2 O emission. Furthermore, these findings improved our understanding of precipitation and N impacts on soil N 2 O emissions and provided useful knowledge for irrigation and N fertilizer management in agriculture to mitigate greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A ModEx Framework for Watershed Subsurface Investigation With Limited Geophysical Data Using Machine Learning and Hydrologic Modeling

Abstract Subsurface heterogeneity influences watershed hydrology strongly but remains difficult to characterize at catchment scales with sparse and costly field data. Geophysical surveys such as electromagnetic induction (EMI) provide local spatial subsurface images yet scaling them to watershed scales and converting EMI‐derived resistivity into hydraulic properties remains a challenge. We present a Model–Experiment (ModEx) framework that integrates limited EMI data with machine learning (ML) and hydrologic modeling to improve process representation and guide field investigations. Sparse EMI surveys were scaled to the catchment scale using a Random Forest model, and the resulting resistivity fields were combined with nearby borehole constraints to parameterize a hydrologic model. The EMI‐informed hydrological simulations improved predictions of streamflow sustained by subsurface flow and shallow saturation patterns. By combining EMI data and ML with hydrologic modeling, the ModEx framework guides future subsurface surveys, providing a transferable and efficient strategy for data–model integration across diverse watersheds. Plain Language Summary Mapping the underground network of soil and rock that controls water is essential for predicting floods and droughts, but seeing underground is difficult and expensive. We cannot drill everywhere, so scientists use geophysical tools to scan broad areas. There are two key challenges: these geophysical scans are often sparse across the whole watershed, and the geophysical data is hard to translate into water‐related properties. We used artificial intelligence to solve these problems. We taught a computer to find patterns linking the limited geophysical data to the land surface properties. This allowed it to fill in the gaps and create a complete, useful subsurface map for the entire watershed. This new map improves hydrologic simulations, leading to more accurate predictions of water movement in the watershed. It also helps scientists build better models with less data and generates a priority map showing where to measure next, making future investigations more efficient. Key Points Limited EMI scaled with ML improves catchment‐scale subsurface parameterization for hydrologic models The framework integrates hydrologic modeling with limited geophysical data to support subsurface investigation design ModEx framework offers a transferable data–model integration strategy that quantifies and reduces uncertainty guiding watershed studies

Chen, Hang↗
DOE OSTI2015

Quantifying sources and sinks of reactive gases in the lower atmosphere using airborne flux observations

Atmospheric composition is governed by the interplay of emissions, chemistry, deposition, and transport. Substantial questions surround each of these processes, especially in forested environments with strong biogenic emissions. Utilizing aircraft observations acquired over a forest in the southeast U.S., we calculate eddy covariance fluxes for a suite of reactive gases and apply the synergistic information derived from this analysis to quantify emission and deposition fluxes, oxidant concentrations, aerosol uptake coefficients, and other key parameters. Evaluation of results against state-of-the-science models and parameterizations provides insight into our current understanding of this system and frames future observational priorities. As a near-direct measurement of fundamental process rates, airborne fluxes offer a new tool to improve biogenic and anthropogenic emissions inventories, photochemical mechanisms, and deposition parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A new process sensitivity index to identify important system processes under process model and parametric uncertainty

Hydrological models are always composed of multiple components that represent processes key to intended model applications. When a process can be simulated by multiple conceptual-mathematical models (process models), model uncertainty in representing the process arises. While global sensitivity analysis methods have been widely used for identifying important processes in hydrologic modeling, the existing methods consider only parametric uncertainty but ignore the model uncertainty for process representation. To address this problem, this study develops a new method to probe multimodel process sensitivity by integrating the model averaging methods into the framework of variance-based global sensitivity analysis, given that the model averaging methods quantify both parametric and model uncertainty. A new process sensitivity index is derived as a metric of relative process importance, and the index includes variance in model outputs caused by uncertainty in both process models and model parameters. For demonstration, the new index is used to evaluate the processes of recharge and geology in a synthetic study of groundwater reactive transport modeling. The recharge process is simulated by two models that converting precipitation to recharge, and the geology process is also simulated by two models of different parameterizations of hydraulic conductivity; each process model has its own random parameters. The new process sensitivity index is mathematically general, and can be applied to a wide range of problems in hydrology and beyond.

Dai, Heng↗
DOE OSTI2024

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GW-PINN: A deep learning algorithm for solving groundwater flow equations

Machine learning methods provide new perspective for more convenient and efficient prediction of groundwater flow. In this study, a deep learning method “GW-PINN” without labeled data for solving groundwater flow equations with wells was proposed. GW-PINN takes the physics inform neural network (PINN) as the backbone and uses either the hard or soft constraint in the loss function for training. A locally refined sampling strategy (LRS) is adopted to generate the consistent spatial sampling points for problems with strong hydraulic head change, and then combined with an appropriate temporal sampling scheme to obtain the final spatial-temporal sampling points. A snowball-style two-stage training strategy by dividing the temporal domain into two subdomains is designed to decrease the sampling points. Five cases were designed to test the training performance of GW-PINN under different sampling strategies and two constraints. The predicted results of GW-PINN were compared with MODFLOW and the analytical solution. The results demonstrate that GW-PINN possesses strong ability in capturing the hydraulic head change for both confined and un-confined aquifers. The hard constraint owns more robust learning ability than the soft constraint. The LRS strategy can generate more accurate results with much fewer sampling points than traditional sampling strategies, and the snowball-style two-stage training strategy is significantly efficient for problems with the drastic change of hydraulic head. Additionally, the application of GW-PINN as a surrogate model for parameterized groundwater flow equations is illustrated. This study provides an option tool for efficient groundwater flow simulation, especially for those with local refinements are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

New perspectives on superparameterization for geophysical turbulence

This is a research expository paper regarding superparameterization, a class of multi-scale numerical methods designed to cope with the intermittent multi-scale effects of inhomogeneous geophysical turbulence where energy often inverse-cascades from the unresolved scales to the large scales through the effects of waves, jets, vortices, and latent heat release from moist processes. Original as well as sparse space–time superparameterization algorithms are discussed for the important case of moist atmospheric convection including the role of multi-scale asymptotic methods in providing self-consistent constraints on superparameterization algorithms and related deterministic and stochastic multi-cloud parameterizations. Test models for the statistical numerical analysis of superparameterization algorithms are discussed both to elucidate the performance of the basic algorithms and to test their potential role in efficient multi-scale data assimilation. The very recent development of grid-free seamless stochastic superparameterization methods for geophysical turbulence appropriate for “eddy-permitting” mesoscale ocean turbulence is presented here including a general formulation and illustrative applications to two-layer quasigeostrophic turbulence, and another difficult test case involving one-dimensional models of dispersive wave turbulence. This last test case has randomly generated solitons as coherent structures which collapse and radiate wave energy back to the larger scales, resulting in strong direct and inverse turbulent energy cascades.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Maximum likelihood Bayesian model averaging and its predictive analysis for groundwater reactive transport models

While Bayesian model averaging (BMA) has been widely used in groundwater modeling, it is infrequently applied to groundwater reactive transport modeling because of multiple sources of uncertainty in the coupled hydrogeochemical processes and because of the long execution time of each model run. To resolve these problems, this study analyzed different levels of uncertainty in a hierarchical way, and used the maximum likelihood version of BMA, i.e., MLBMA, to improve the computational efficiency. Our study demonstrates the applicability of MLBMA to groundwater reactive transport modeling in a synthetic case in which twenty-seven reactive transport models were designed to predict the reactive transport of hexavalent uranium (U(VI)) based on observations at a former uranium mill site near Naturita, CO. Moreover, these reactive transport models contain three uncertain model components, i.e., parameterization of hydraulic conductivity, configuration of model boundary, and surface complexation reactions that simulate U(VI) adsorption. These uncertain model components were aggregated into the alternative models by integrating a hierarchical structure into MLBMA. The modeling results of the individual models and MLBMA were analyzed to investigate their predictive performance. The predictive logscore results show that MLBMA generally outperforms the best model, suggesting that using MLBMA is a sound strategy to achieve more robust model predictions relative to a single model. MLBMA works best when the alternative models are structurally distinct and have diverse model predictions. When correlation in model structure exists, two strategies were used to improve predictive performance by retaining structurally distinct models or assigning smaller prior model probabilities to correlated models. Since the synthetic models were designed using data from the Naturita site, the results of this study are expected to provide guidance for real-world modeling. Finally, limitations of applying MLBMA to the synthetic study and future real-world modeling are discussed.

58 GEOSCIENCES↗
DOE OSTI2023

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Dislocation dynamics in hexagonal close-packed crystals

Extensions of the dislocation dynamics methodology necessary to enable accurate simulations of crystal plasticity in hexagonal close-packed (HCP) metals are presented. They concern the introduction of dislocation motion in HCP crystals through linear and non-linear mobility laws, as well as the treatment of composite dislocation physics. Formation, stability and dissociation of and other dislocations with large Burgers vectors defined as composite dislocations are examined and a new topological operation is proposed to enable their dissociation. Furthermore, the results of our simulations suggest that composite dislocations are omnipresent and may play important roles both in specific dislocation mechanisms and in bulk crystal plasticity in HCP materials. While fully microscopic, our bulk DD simulations provide wealth of data that can be used to develop and parameterize constitutive models of crystal plasticity at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Seasonal and Morphological Controls on Nitrate Retention in Arctic Deltas

Abstract Estimates of nitrate loading to the Arctic Ocean are limited by the lack of field observations within deltas partly due to logistical constraints. To overcome this limitation, we use a remote sensing framework to estimate retention of nitrate in Arctic deltas. We achieve this by coupling hydrological and biogeochemical process models at the network scale for five major Arctic deltas. Binary masks of delta channels were used to simulate flow direction and magnitude through networks. Models were parameterized using historical and seasonal observations. Simulated nitrate retention ranged from 2.9% to 15% of the incoming load. Retention rates were largest during winter but smallest during spring conditions when increased discharges export large nitrate masses to the coast. Under future climate scenarios, retention rates fall by ∼1%–10%. Arctic deltas have an important effect on the magnitude of nitrate entering Arctic seas and the inclusion of processing in deltas can improve flux estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uptake of N 2 O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗