SEARCH · Engineering Papers
Results for “MODELING”
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.
Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data
Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.
Modelling the sulfate capacity of simulated radioactive waste borosilicate glasses
The capacity of simulated high-level radioactive waste borosilicate glasses to incorporate sulfate has been studied as a function of glass composition. Combined Raman, 57 Fe Mössbauer and literature evidence supports the attribution of coordination numbers and oxidation states of constituent cations for the purposes of modelling, and results confirm the validity of correlating sulfate incorporation in multicomponent borosilicate radioactive waste glasses with different models. A strong compositional dependency is observed and this can be described by an inverse linear relationship between incorporated sulfate (mol% SO 4 2- ) and total cation field strength index of the glass, Σ(z/a 2 ), with a high goodness-of-fit (R 2 ≈ 0.950). Similar relationships are also obtained if theoretical optical basicity, Λ th (R 2 ≈ 0.930) or non-bridging oxygen per tetrahedron ratio, NBO/T (R 2 ≈ 0.919), are used. Results support the application of these models, and in particular Σ(z/a 2 ), as predictive tools to aid the development of new glass compositions with enhanced sulfate capacities.
A staggered-grid finite-difference scheme optimized in the time–space domain for modeling scalar-wave propagation in geophysical problems
For modeling scalar-wave propagation in geophysical problems using finite-difference schemes, optimizing the coefficients of the finite-difference operators can reduce numerical dispersion. Most optimized finite-difference schemes for modeling seismic-wave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finite-difference scheme for numerical scalar-wave modeling to control dispersion errors not only in space but also in time. Our optimized scheme is based on a new stencil that contains a few more grid points than the standard stencil. We design an objective function for minimizing relative errors of phase velocities of waves propagating in all directions within a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finite-difference scheme is computationally up to 2.5 times faster than the optimized schemes using the standard stencil to achieve the similar modeling accuracy for a given 2D or 3D problem. Compared with the high-order finite-difference scheme using the same new stencil, our optimized scheme reduces 50 percent of the computational cost to achieve the similar modeling accuracy. This new optimized finite-difference scheme is particularly useful for large-scale 3D scalar-wave modeling and inversion.
A Bayesian framework for adaptive selection, calibration, and validation of coarse-grained models of atomistic systems
A general adaptive modeling algorithm for selection and validation of coarse-grained models of atomistic systems is presented. A Bayesian framework is developed to address uncertainties in parameters, data, and model selection. Algorithms for computing output sensitivities to parameter variances, model evidence and posterior model plausibilities for given data, and for computing what are referred to as Occam Categories in reference to a rough measure of model simplicity, make up components of the overall approach. Computational results are provided for representative applications.
Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows
We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number RE = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.
Effect of grain morphology on gas bubble swelling in UMo fuels – A 3D microstructure dependent Booth model
Not Available
Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model
Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.
Modeling of flame extinction/re-ignition in oxygen-reduced environments
Not Available
Towards robust and scalable dispatch modeling of long-duration energy storage
Not Available
Machine Learning Applied to a Variable Charge Atomistic Model for Cu/Hf Binary Alloy Oxide Heterostructures [Machine Learning a Variable Charge Atomistic Model for Cu/Hf Binary Alloy Oxide Heterostructures]
Alloy oxide heterostructures are essential for a wide variety of energy applications, including anticorrosion coatings, bifunctional catalysis, energy storage, as well as emerging platforms for multilevel nonvolatile memory and neuromorphic devices. The key role played by oxide composition, density, and stoichiometry in governing electrochemical reactions, interface structure, and ionic transport in alloy oxides remains largely unknown. This is primarily due to the lack of variable charge interatomic potential models that can adequately describe the various atomic/ionic interactions in alloy oxides within a unified framework. Here, we introduce a charge transfer ionic potential (CTIP) model for Cu/Hf/O alloy system and demonstrate its ability to accurately capture the complex potential energy surface owing to dynamically varying interactions, including metallic (Cu/Hf), ionic (Cu/Hf/Ox), mixed environment (interfaces). We leverage supervised machine learning methods powered by genetic algorithms coupled with local simplex optimization to efficiently navigate through a high-dimensional parameter space and identify an optimum set of an independent set of CTIP parameters. We train our model against an extensive first-principles based data set that includes lattice constants, cohesive energies, equations of state, and elastic constants of various experimentally observed polymorphs of hafnia, copper oxide, hafnium, as well as copper. Our machine-learned CTIP model captures the structure, elastic properties, thermodynamics, and energetic ordering of various polymorphs. It also accurately predicts the surface properties of both oxides and their metals. To demonstrate the suitability of our CTIP model for investigating dynamic processes, we employed it to identify the atomistic scale mechanisms associated with the initial nanoscale oxidation and surface oxide growth on Cu-doped hafnia surfaces. Here, our machine-learned CTIP model can be used to probe the dynamic response of Cu/Hf/O alloy interfaces subjected to external stimuli (e.g., electric field, pressure, temperature, strain, etc.) and a variety of atomistic phenomena including the dynamics of switching in emerging neuromorphic platforms.
A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes
To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.
Modeling the Impacts of Hydrogen Extraction on Anaerobic Fermentation Systems Using a Modified ADM1 Model
Here, this paper presents a modified ADM1 (mADM1) model for use in assessing the influence of active dissolved H 2 extraction on biogas production from anaerobic fermentation reactors, specifically in the context of high-strength brewery wastewater as the feed. The modified model adds regulation factors to account for the varying product stoichiometry of glucose degradation reactions as a function of H 2 concentration in anaerobic digestion, and incorporates the use of both lactate and ethanol as intermediates. The changes to the model are expected to better predict the effects of H 2 removal on the production of H 2 in fermentation processes. Both ADM1 and mADM1 were calibrated against experimental data from a fluidized-bed reactor containing encapsulated acidogenic-acetogenic microbial communities treating real brewery wastewater. The calibrated models were validated against another fluidized-bed reactor identical to the calibration reactor, with the addition of a hollow-fiber membrane module that was used to reduce dissolved H 2 concentrations in the validation reactor. It was found that ADM1 over-estimated the H 2 production rates of the validation reactor (117 mL/day simulated vs. 16.48 mL/day experimental), whereas mADM1 was able to successfully predict fermentation, and specifically the H 2 production rate (17 mL/day). The presented model is a first step towards a model that can predict the effects of dissolved H 2 removal on fermentation processes.
Predictive Model for Particle Residence Time Distributions in Riser Reactors. Part 1: Model Development and Validation
Here in this computational study, we model the mixing of biomass pyrolysis vapor with solid catalyst in circulating riser reactors with a focus on the determination of solid catalyst residence time distributions (RTDs). A comprehensive set of 2D and 3D simulations were conducted for a pilot-scale riser using the Eulerian-Eulerian two-fluid modeling framework with and without sub-grid-scale models for the gas-solids interaction. A validation test case was also simulated and compared to experiments, showing agreement in the pressure gradient and RTD mean and spread. For simulation cases, it was found that for accurate RTD prediction, the Johnson and Jackson partial slip solids boundary condition was required for all models and a sub-grid model is useful so that ultra high resolutions grids that are very computationally intensive are not required. Finally, we discovered a 2/3 scaling relation for the RTD mean and spread when comparing resolved 2D simulations to validated unresolved 3D sub-grid-scale model simulations.
Improving Land Surface Temperature Simulation in CoLM Over the Tibetan Plateau through Fractional Vegetation Cover Derived from a Remotely Sensed Clumping Index and Model-Simulated Leaf Area Index
Parameterizations of fractional vegetation cover (FVC) in land surface models have important effects on simulations of surface energy budget, especially in arid and semiarid regions. This study uses a FVC scheme in which FVC is derived from leaf area index and a remotely sensed clumping index. The performance of the new scheme (SMFVC) is evaluated against Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) and in situ soil temperature observations, together with two other FVC schemes, a general FVC scheme (CTL) based on land cover map and a climatology-based FVC scheme (RSFVC) that uses long-term remotely sensed Normalized Difference Vegetation Index of MODIS. The three FVC schemes were implemented in the Common Land Model (CoLM) and applied in the Tibetan Plateau using the same forcing data and default parameters. Our results demonstrate that FVC schemes have significant influence on the CoLM performance: (1) the RSFVC and SMFVC schemes significantly reduce the LST biases found in CTL, particularly in grassland and during summer; (2) soil temperature evaluation by in situ observations from three networks on the Tibetan Plateau corroborates the LST results; and (3) the improvements are mainly related to representing temporal (seasonal) variability and subgrid heterogeneity of FVC, which improves surface albedo and surface energy balance. In other words, by including more vegetation characteristics, such as using a clumping index, land surface models may better simulate surface vegetation condition and further better represent the land surface energy budget over the Tibetan Plateau.
Effects of groundwater pumping on ground surface temperature: A regional modeling study in the North China Plain
Over-exploitation of groundwater (GW) in the North China Plain (NCP) since the 1960s has many environmental consequences. However, understanding of the dominant mechanisms remains limited, particularly at regional scale. In this study, the coupled ParFlow.CLM model representing subsurface and land-surface processes and their interactions was applied in the NCP at high spatio-temporal resolutions. The model was validated using the water and energy fluxes reported in previous studies and from the JRA-55 reanalysis. Numerical experiments were designed to examine the relative impacts of GW pumping and irrigation on the ground surface temperature (GST). Results show significant effects of GW pumping on GST in the NCP. Generally, the subsurface acts as a buffer to temporal variations in heat fluxes at the land-surface, but long-term pumping can gradually weaken this buffer, resulting in increases in the spatio-temporal variability of GST, as exemplified by hotter summers and colder winters. Considering that changes of water table depth (WTD) can significantly affect land surface heat fluxes when WTD ranges roughly between 0.01–10 m, the 0.5 m/year increase of WTD simulated by the model due to pumping can continue to increase the average WTD and hence, disturb GST, for about 20 years from the pre-pumping WTD in the NCP, before the WTD exceeds 10 m. The variations of GST are expected to be faster initially and gradually slow down due to the nonlinear behaviors of GST with WTD. The findings from this study in the NCP may also have implications for other regions with GW depletion.
An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis
Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.
Mechanistic Modeling of Microtopographic Impacts on CO 2 and CH 4 Fluxes in an Alaskan Tundra Ecosystem Using the CLM-Microbe Model
Spatial heterogeneities in soil hydrology have been confirmed as a key control on CO 2 and CH 4 fluxes in the Arctic tundra ecosystem. In this study, we applied a mechanistic ecosystem model, CLM-Microbe, to examine the microtopographic impacts on CO 2 and CH 4 fluxes across seven landscape types in Utqiaġvik, Alaska: trough, low-centered polygon (LCP) center, LCP transition, LCP rim, high-centered polygon (HCP) center, HCP transition, and HCP rim. We first validated the CLM-Microbe model against static-chamber measured CO 2 and CH 4 fluxes in 2013 for three landscape types: trough, LCP center, and LCP rim. Model application showed that low-elevation and thus wetter landscape types (i.e., trough, transitions, and LCP center) had larger CH 4 emissions rates with greater seasonal variations than high-elevation and drier landscape types (rims and HCP center). Sensitivity analysis indicated that substrate availability for methanogenesis (acetate, CO 2 + H 2 ) is the most important factor determining CH 4 emission, and vegetation physiological properties largely affect the net ecosystem carbon exchange and ecosystem respiration in Arctic tundra ecosystems. Modeled CH 4 emissions for different microtopographic features were upscaled to the eddy covariance (EC) domain with an area-weighted approach before validation against EC-measured CH 4 fluxes. The model underestimated the EC-measured CH 4 flux by 20% and 25% at daily and hourly time steps, suggesting the importance of the time step in reporting CH 4 flux. The strong microtopographic impacts on CO 2 and CH 4 fluxes call for a model-data integration framework for better understanding and predicting carbon flux in the highly heterogeneous Arctic landscape.