Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DYNAMICS”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Thermal Interface Conductance Between Aluminum and Silicon by Molecular Dynamics Simulations

The thermal interface conductance between Al and Si was simulated by a non-equilibrium molecular dynamics method. In the simulations, the coupling between electrons and phonons in Al are considered by using a stochastic force. The findings demonstrate the size dependence of the interface thermal conductance and the effect of electron–phonon coupling on the interface thermal conductance. To understand the mechanism of interface resistance, the vibration power spectra are calculated. We postulate that the atomic level disorder near the interface is an important aspect of interfacial phonon transport, which leads to a modification of the phonon states near the interface. There, the vibrational spectrum near the interface greatly differs from the bulk. This change in the vibrational spectrum affects the results predicted by AMM and DMM theories and indicates new physics is involved with phonon transport across interfaces.

42 ENGINEERING↗
DOE OSTI2015

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models. Current status and future directions

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and loss from soil accounts for a large pro portion of land-atmosphere C exchange. Due to large pool size and variable residence time from years to millennia, even small changes in soil organic C(SOC) have substantial effects on the terrestrial C budget, thereby affecting atmospheric carbon dioxide (CO2)concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain and identifying major driving forces controlling soil C storage and fluxes remains a key research challenge his study has compiled century-long (1901-2010)estimates of SOC storage and heterotrophic respiration (Rh) from ten terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) and two observation based datasets. The ten-TBM ensemble shows that global SOC estimate range from 4 to 2111 Pg C (1 Pg = 10 15 g) with a median value of 1158 Pg C33 in 2010. Modeling approach estimates a broad range of Rh from 35 to 69 Pg C yr -1 with a median value of 51Pg C yr -1 during 200–2010. The largest uncertainty in SOC stocks exists in the 40–65°N latitude band while Rh differences are the largest in the tropics. All the models agreed that climate and land use changes have decreased SOC stocks while elevated CO 2 and atmospheric nitrogen deposition have increased SOC stocks though the response varied significantly among models. Model representations of temperature and moisture sensitivity,nutrient limitation and land use partially explain the divergent estimates of global SOC stocks and soil fluxes in this study. In addition, major sources of uncertainty from model estimation include exclusion of SOC storage in wetlands and peatlands as well as C storage in deep soil layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Intramolecular Hydrogen Bonding Restricts Gd-Aqua-Ligand Dynamics [The Day the Water Stood Still: Intramolecular Hydrogen Bonding to Restrict Gd-Aqua Ligand Dynamics]

Aqua ligands can undergo rapid internal rotation about the M-O bond. For magnetic resonance contrast agents, this rotation results in diminished relaxivity. Herein, we show that an intramolecular hydrogen bond to the aqua ligand can reduce this internal rotation and increase relaxivity. Molecular modeling was used to design a series of four Gd complexes capable of forming an intramolecular H-bond to the coordinated water ligand, and these complexes had anomalously high relaxivities compared to similar complexes lacking a H-bond acceptor. Molecular dynamics simulations supported the formation of a stable intramolecular H-bond, while alternative hypotheses that could explain the higher relaxivity were systematically ruled out. Finally, intramolecular H-bonding represents a useful strategy to limit internal water rotational motion and increase relaxivity of Gd complexes.

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

Interface and embedded joint methods for modeling dynamic fracture opening by explosive products [Two methods for modeling dynamic fracture opening by explosive products]

Two computational approaches are proposed in the paper to model dynamic fracture opening by explosive products. The first method assumes that the fractures may be modeled using flow elements embedded along the mesh lines. This method models crack opening in a straightforward way by splitting the nodes of the computational grid. It can account for crack branching; however, the crack directions are constrained by existing mesh faces, which may lead to mesh dependence. Also, the stress in flow elements is calculated explicitly separate from the surrounding solid elements that can impose additional limits on the time step stability condition for explicit integration. The second approach uses embedded flow elements to model the cracks. Typical thickness of the cracks is much smaller than the element size. Therefore, gas pressure in the cracks is assumed to be in stress equilibrium with the element stress. To achieve this, the crack thickness and the state of the gas is updated simultaneously with the state of the solid element which contains the crack. Therefore, the time step is controlled by the explicit solver applied for the solid and does not depend on the thickness of the crack. In conclusion, the main disadvantage of the second approach is due to the complexity of modeling multiple intersecting cracks, which go through the same element. We discuss the areas of possible applications of these 2 methods and the ways to improve and enhance them for future practical applications.

58 GEOSCIENCES↗
DOE OSTI2017

Error modeling for surrogates of dynamical systems using machine learning: Machine-learning-based error model for surrogates of dynamical systems

A machine learning–based framework for modeling the error introduced by surrogate models of parameterized dynamical systems is proposed. The framework entails the use of high-dimensional regression techniques (eg, random forests, and LASSO) to map a large set of inexpensively computed “error indicators” (ie, features) produced by the surrogate model at a given time instance to a prediction of the surrogate-model error in a quantity of interest (QoI). This eliminates the need for the user to hand-select a small number of informative features. The methodology requires a training set of parameter instances at which the time-dependent surrogate-model error is computed by simulating both the high-fidelity and surrogate models. Using these training data, the method first determines regression-model locality (via classification or clustering) and subsequently constructs a “local” regression model to predict the time-instantaneous error within each identified region of feature space. We consider 2 uses for the resulting error model: (1) as a correction to the surrogate-model QoI prediction at each time instance and (2) as a way to statistically model arbitrary functions of the time-dependent surrogate-model error (eg, time-integrated errors). We then apply the proposed framework to model errors in reduced-order models of nonlinear oil-water subsurface flow simulations, with time-varying well-control (bottom-hole pressure) parameters. The reduced-order models used in this work entail application of trajectory piecewise linearization in conjunction with proper orthogonal decomposition. Moreover, when the first use of the method is considered, numerical experiments demonstrate consistent improvement in accuracy in the time-instantaneous QoI prediction relative to the original surrogate model, across a large number of test cases. When the second use is considered, results show that the proposed method provides accurate statistical predictions of the time- and well-averaged errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hydrodynamics of countercurrent flows in a structured packed column: effects of initial wetting and dynamic contact angle

The structured packed column has gained considerable attention as a prospective technology to mitigate carbon dioxide emissions from thermal power plants. The computational flow investigation of such systems is a multiscale problem, and local hydrodynamics plays a key role in overall column efficiency. Accordingly, multiphase flow studies in the representative elementary unit (REU) of a packed column can offer insights into the local flow behavior such as flow pattern, interfacial area, etc. The interfacial area plays a key role in the overall interphase mass transfer and thereby in separation efficiency. The predicted effective area matches well with corresponding ones in experiments for Mellapak 250.Y packings. The contact angle (i.e., solid substrate behavior) is one of the critical factors that dictates wettings, and thereby the interfacial area. The dynamic contact angle (DCA) was considered in order to explore the effects of contact angle hysteresis on the interfacial area. DCA has more pronounced impact on the interfacial area for a solvent possessing a higher value of surface tension and equilibrium contact angle. The interfacial area shows temporal undulation and does not achieve a pseudo steady state due to the capillary force. In contrary, the interfacial area gets a net value for a solvent having low surface tension value. Further, specification of the static contact angle at the solid substrate shows the least value of the interfacial area. We also extensively studied the effect of the initial sheet condition (dry vs wet) on the interfacial area. The initially wetted sheets show slightly higher value of the interfacial area as compared to the initially dry sheets at a fixed flow rate for a solvent having high surface tension value. Wetting hysteresis decreases with increasing liquid loads. On the other hand, wetting hysteresis does not occur for a solvent having value of low surface tension and equilibrium contact angle.

structured packings, dynamic contact angle, interf↗
DOE OSTI2018

Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution

Hydrodynamics, oxygen transfer, and oxygen uptake in bubble-column and airlift bioreactors were studied using multiphase Euler-Euler computational fluid dynamics (CFD) simulations. Interphase mass transfer of oxygen was modeled with a liquid mass-transfer coefficient and sustained driving force provided by a phenomenological model for microbial oxygen uptake rate (OUR). Laboratory-scale reactor simulations of a bubble-column (0.15 m diameter and 0.75 m initial liquid height) showed reasonable agreement with gas holdup and mass transfer experiments reported in the literature. Commercial-scale bubble-column and airlift reactor simulations (5 m diameter and 25 m initial liquid height) were simulated using this validated model. Similar trends for the variation of overall gas holdup and oxygen concentrations with superficial velocity were observed for both of the commercial-scale reactors. However, the simulation results indicate differences in hydrodynamics, such as better recirculation in the airlift reactors that enables favorable distribution of oxygen in the reactor.

09 BIOMASS FUELS↗
DOE OSTI2022

Membrane synthesis via in-situ pore formation in silica gels through dynamic miscibility with soybean oil

The dynamic immiscibility between vegetable oil and organic silica precursors during sol-gel synthesis was used to tailor the microstructure of gels made with tetraethoxysilane (TEOS) or tetraethoxysilane+methyltrimethoxysilane (MTMS). Vegetable oil was investigated as an environmentally friendly and cheaper alternative to many other porogens. Oil concentrations of 9.09 v/v% and 33.33 v/v% were studied. For both TEOS and TEOS+MTMS gels, the miscibility of the oil and precursor decreases as the silica gel network forms, leading to oil-rich regions within the silica matrix. Gels containing MTMS maintained partial miscibility with the oil during gel formation due to the non-hydrolysable methyl groups on the MTMS. The pore structure of MTMS gels appeared more interconnected compared to the isolated pores observed in gels made with only TEOS. To create membranes with hierarchical pore structures, solvent exchanges were performed on gels or dried gels (i.e., xerogels) were calcined. The effect of the methyl group was also observed in the results of the time sweep rheology, as the addition of oil to the TEOS+MTMS gels caused an increase in the time to form the spanning gel. The trend was the opposite for oil additions to TEOS gels, which only contained polar groups. Xerogel membranes made with TEOS+MTMS had higher specific surface areas (SSAs) than those made with TEOS. For oil concentrations of 33.33 v/v%, xerogels made from TEOS+MTMS had an average SSA of 1068 + 113 m2 g-1, an order of magnitude higher than the average SSA of 103 + 6 m2 g-1 for xerogels made with TEOS.

sol-gel, dynamic immiscibility, oil-alkoxide inter↗
DOE OSTI2025

Mechanisms and stability of Li dynamics in amorphous Li-Ti-P-S-based mixed ionic–electronic conductors: A machine learning molecular dynamics study

Mixed ionic–electronic conductors (MIECs) exhibit both high ionic and electronic conductivity to improve the battery performance. In this work, we investigate the mechanism and stability of transport channels in our recently developed MIEC material, amorphous Ti-doped lithium phosphorus sulfide (LPS), using molecular dynamics (MD) simulations with a 99% accurate machine-learning force field (MLFF) trained on ab initio MD data. The achieved MLFF helps efficient large-scale MD simulations on LPS with three Ti concentrations (10%, 20%, and 30%) and six temperatures (25°C to 225°C) to calculate ionic conductivity, activation energy, Li-ion transport mechanism, and configurational entropy. Results show that ionic conductivities and activation energies are consistent with our recent experimental values. Moreover, Li-ion transport occurs via free-volume diffusion facilitated by the formation of disordered Li-S polyhedra. The enhanced stability of transport channels at 10% and 20% Ti doping, compared to 0% and 30%, is observed by analyzing the vibrational and configurational entropy of these disordered Li-S polyhedra. Overall, this study highlights the utility of MLFF-based large-scale MD simulations in explaining the transport mechanism and the stability of Li-ion in Ti-doped LPS electrolyte with significant computational efficiency.

And configuration entropy↗
DOE OSTI2015

Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model—the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4-fold for a two-level implementation, and can be increased up to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of r{sub s}=0.912.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗