Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Understanding Proton Movement in [Fe-Fe] Hydrogenases

Nature uses specialized metalloenzymes to carry out small molecule activation reactions, including CO 2 fixation, O 2 activation, and proton reduction, with unparalleled efficiency, rates, and selectivity. The latter reactions are performed by hydrogenases, protein metallo-complexes that interconverts H 2 to protons and electrons (H 2 oxidation) and the reverse reaction (H 2 production) with incredibly low energy input and amazingly fast kinetics. Reproducing both the activity and efficiency of metalloenzymes in sustainable anthropogenic systems remains one of the “holy grails” of inorganic chemistry. However, identifying the precise molecular components responsible for these desirable properties has been challenging in the natural metalloenzymes, hindering efforts to develop analogous processes in synthetic compounds. Considering the inherent complexity of a metalloenzyme and the many interactions, both strong and weak, that contribute to the function of an enzyme, we have elected to model natural metalloenzymes on a biochemical platform. Towards this end, we have a developed structural, functional, and mechanistic mimic of the [Ni-Fe] hydrogenases within a robust protein scaffold, rubredoxin, to understand the influence of the secondary coordination environment on the metal center. This involved preparing a series of three rubredoxin constructs containing a single point mutation at Val positions and making the NMR chemical shift assignments for the paramagnetic (nickel-substituted) and non-paramagnetic (zinc-substituted) form. These physical studies were complemented with in silico molecular dynamic studies on proteins with [Fe-S] clusters to engineer new proton channels to test in vitro . To assist our search for new natural metalloprotein scaffolds within the vast number of sequenced genomes, we created a neural network-based program to identify proteins with specific metal-binding sites. These computational and physical studies with metalloenzymes provide direct insight into the fundamental chemical principles driving the natural systems and offer design principles for developing catalysts that utilize analogous principles.

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

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

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

Optimization of Magnesium-Doped Lithium Metal Anode for High Performance Lithium Metal Batteries through Modeling and Experiment

Lithium (Li)-magnesium (Mg) alloy with limited Mg amount, which can also be called Mg-doped Li (Li-Mg), has been considered as a potential alternative anode for high energy density rechargeable Li metal batteries. However, the optimum doping-content of Mg in Li-Mg anode and the mechanism of the improved performance are not well understood. In this study, density functional theory (DFT) calculations are used to investigate the effect of Mg amount in Li-Mg anode. The Li-Mg with about 5 wt. % Mg (abbreviated as Li-Mg5) has the lowest absorption energy of Li, thus all the surface area can be “controlled” by Mg atoms, leading to the smooth and continuous deposition of Li on the surface around the Mg center. A localized high concentration electrolyte enables Li-Mg5 to exhibit the best cycling stability in Li metal batteries with high-loading cathode and lean electrolyte under 4.4 V high-voltage, which is approaching the demand of practical application. This electrolyte also helps generate an inorganic-rich solid electrolyte interphase, which leads to smooth, compact and less corrosion layer on the Li-Mg5 surface. Both theoretical simulations and experimental results prove that Li-Mg5 has optimum Mg content and gives best battery cycling performance.

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

Adsorption, Orientation, and Speciation of Amino Acids at Air–Aqueous Interfaces for the Direct Air Capture of CO 2

Amino acids make up a promising family of molecules capable of direct air capture (DAC) of CO 2 from the atmosphere. Under alkaline conditions, CO 2 reacts with the anionic form of an amino acid to produce carbamates and deactivated zwitterionic amino acids. The presence of the various species of amino acids and reactive intermediates can have a significant effect on DAC chemistry, the role of which is poorly understood. In this study, all-atom molecular dynamics (MD) based computational simulations and vibrational sum frequency generation (vSFG) spectroscopy studies were conducted to understand the role of competitive interactions at the air–aqueous interface in the context of DAC. We find that the presence of potassium bicarbonate ions, in combination with the anionic and zwitterionic forms of amino acids, induces concentration and charge gradients at the interface, generating a layered molecular arrangement that changes under pre- and post-DAC conditions. In parallel, an enhancement in the surface activity of both anionic and zwitterionic forms of amino acids is observed, which is attributed to enhanced interfacial stability and favorable intermolecular interactions between the adsorbed amino acids in their anionic and zwitterionic forms. The collective influence of these competitive interactions, along with the resulting interfacial heterogeneity, may in turn affect subsequent capture reactions and associated rates. Finally, these effects underscore the need to consider dynamic changes in interfacial chemical makeup to enhance DAC efficiency and to develop successful negative emission and storage technologies.

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

Diffusion-Limited Kinetics of Isovalent Cation Exchange in III–V Nanocrystals Dispersed in Molten Salt Reaction Media

The goal of this work is to determine the kinetic factors that govern isovalent cation exchange in III–V colloidal quantum dots using molten salts as the solvent and cation source. We focus on the reactions of InP + GaI 3 → In 1 – x Ga x P and InAs + GaI 3 → In 1–x Ga x As to create technologically important ternary III–V phases. We find that the molten salt reaction medium causes the transformation of nearly spherical InP nanocrystals to tetrahedron-shaped In 1–x Ga x P nanocrystals. Furthermore, we determine that the activation energy for the cation exchange reaction is 0.9 eV for incorporation of Ga into InP and 1.2 eV for incorporation of Ga into InAs, both much lower than the measured values in bulk semiconductors. Next, we use powder XRD simulations to constrain our understanding of the structure of the In 1–x Ga x P nanocrystals. Together our results reveal several important features of molten salt-mediated cation exchange and provide guidance for future development of these materials.

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

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗
DOE OSTI2022

Ion Transport in 2D Nanostructured $\pi$-Conjugated Thieno[3,2- b ]thiophene-Based Liquid Crystal

Leveraging the self-assembling behavior of liquid crystals designed for controlling ion transport is of both fundamental and technological significance. Here, we have designed and prepared a liquid crystal that contains (2,5-bis(5-(2,5,8,11-tetraoxatridecan-13-yl)thiophen-2yl)thieno[3,2-b]thiophene (BTTT) as mesogenic core and conjugated segment and symmetric tetra(ethylene oxide) (EO4) as polar side chains for ionic conducting regions. Driven by the crystallization of the BTTT cores, BTTT/dEO4 exhibits well-ordered smectic phases below 71.5 °C as confirmed by differential scanning calorimetry, polarized optical microscope, temperature dependent wide-angle X-ray scattering and grazing incidence wide-angle X-ray scattering (GIWAXS). We adopted a combination of experimental GIWAXS and molecular dynamics (MD) simulations to better understand the molecular packing of BTTT/dEO4 films, particularly when loaded with the ion conducting salt, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI). Ionic conduction of BTTT/dEO4 is realized by the addition of LiTFSI, with the material able to maintain smectic phases up to r = [Li + ]/[EO] = 0.1. In this work, The highest ionic conductivity of 8 × 10 -3 S/cm was attained at an intermedium salt concentration of r = 0.05. It was also found that ion conduction in BTTT/dEO4 is enhanced by forming a smectic layered structure with irregular interfaces between the BTTT and EO4 layers and by the lateral film expansion upon salt addition. This can be explained by the enhancement of the misalignment and configurational entropy of the side chains, which increase their local mobility and that of the solvated ions. Our molecular design thus illustrates how, beyond the favorable energetic interactions that drive the assembly and ion solvating domains, modulation of entropic effects can also be favorably harnessed to improve ion conduction.

thin film electrochemical impedance spectroscopy↗
DOE OSTI2024

Detecting underscreening and generalized Kirkwood transitions in aqueous electrolytes

We establish the connection between the measured small angle x-ray scattering signal and the charge–charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a 0 , and the inverse periodicity length, Q 0 . Above the KTs, we find a universal scaling of a 0 ∝ c – $\sqrt{ζ/3}$ and Q 0 ∝ c 1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor.

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

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

Here, in this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

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

Modeling the Photoelectrochemical Evolution of Lead-Based, Mixed-Halide Perovskites Due to Photosegregation

Lead-based, mixed-halide perovskites such as methylammonium lead iodide-bromide [MAPb(I 1-x Br x ) 3 ] undergo anion photosegregation under illumination. This is observed as low band gap photoluminescence from photogenerated iodine-rich domains due to favorable band offsets that induce carrier funneling into them. Unfortunately, theoretical rationalizations of mixed-halide photosegregation are complicated by biases inherent to photoluminescence-based observations. Recent compositionally-weighted X-ray diffraction (XRD) measurements now reveal broad distributions of photosegregated stoichiometries not captured by existing photosegregation models. To better bridge experiment and theory, we perform kinetic Monte Carlo (KMC) simulations of photosegregation within the context of a band gap-based thermodynamic model, which has previously accounted for numerous experimental observations. Our KMC simulations are modified to consider high carrier density Fermi-Dirac statistics that result from carrier funneling and accumulation within photosegregated I-rich domains. Obtained KMC results reproduce broad x terminal distributions seen experimentally and illustrate how they are characterized by a central, heavily I-enriched stoichiometry. I-rich domain “drifting” during photosegregation rationalizes the long photosegregation timescales seen experimentally with drifting simultaneously producing a wake of variable stoichiometry I-rich inclusions that form the lion’s share of stoichiometric heterogeneities seen in compositionally-weighted XRD measurements. These simulations and accompanying rationalizations further reveal a general criterion for realizing favorable free energies to induce demixing. Central to the criterion is the statistical occupation of low gap inclusions in the parent alloy by excitations. Furthermore, the resulting model thus provides a general framework for conceptualizing mixed-halide perovskite light and temperature sensitivities, mediated by photocarriers.

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

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Analytical solution for the cumulative wake of wind turbines in wind farms

This paper solves an approximate form of conservation of mass and momentum for a turbine in a wind farm array. The solution is a fairly simple explicit relationship that predicts the streamwise velocity distribution within a wind farm with an arbitrary layout. As this model is obtained by solving flow-governing equations directly for a turbine that is subject to upwind turbine wakes, no ad hoc superposition technique is needed to predict wind farm flows. A suite of large-eddy simulations (LES) of wind farm arrays is used to examine self-similarity as well as validity of the so-called conservation of momentum deficit for turbine wakes in wind farms. The simulations are performed with and without the presence of some specific turbines in the wind farm. This allows us to systematically study some of the assumptions made to develop the analytical model. A modified version of the conservation of momentum deficit is also proposed to provide slightly better results at short downwind distances, as well as in the far wake of turbines deep inside a wind farm. Model predictions are validated against the LES data for turbines in both full-wake and partial-wake conditions. While our results highlight the limitation in capturing the flow speed-up between adjacent turbine columns, the model is overall able to acceptably predict flow distributions for a moderately sized wind farm. Finally, the paper employs the new model to provide insights on the accuracy of common wake superposition methods.

17 WIND ENERGY↗
DOE OSTI2021

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗
DOE OSTI2021

WaterTAP3 (The Water Technoeconomic Assessment Pipe-Parity Platform)

The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience.

Miara, Ariel↗
DOE OSTI2021

Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3)

The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience. A publicly available graphical user interface is currently under development.

Miara, Ariel↗