Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High temperature gamma radiation-induced chromium redox chemistry via in situ spectroscopic measurements

Chromium ions can make their way into the primary coolant of nuclear power reactors from the corrosion of stainless-steel reactor components, decreasing the material's corrosion resistance and resulting in increased transport of further corrosion products. Despite these potential effects, the radiation-induced redox speciation of chromium ions in aqueous solution is not well understood, especially at the elevated temperatures experienced by reactor coolants. Here, in the present work, we report new experimental results demonstrating that in aerated aqueous solution, the radiolytic oxidation of Cr(III) to Cr(VI) occurs at pH 4, while the reduction of Cr(VI) to Cr(III) occurs at pH 2. The oxidation of Cr(III) is primarily attributed to the reaction of the hydroxyl radical (˙OH) with the Cr(OH) 2+ species, while the reduction of Cr(VI) is attributed to reactions involving the hydrated electron (e aq - ) and hydrogen atom (H˙). Additionally, the steady-state equilibrium yield of Cr(VI) from the gamma irradiation of pH 4 Cr(III) solutions decreased with increasing temperature (over a range of 37–195 °C). This observation indicates that the activation energy of the Cr(VI) reduction reactions is higher than that for the Cr(III) oxidation reactions, such that it becomes relatively more favorable at higher temperatures. Overall, these data are important for the development of complementary multiscale models for the prediction of metal ion speciation in high temperature radiation environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2014

A novel method for calculating the energy barriers for carbon diffusion in ferrite under heterogeneous stress

A novel method for accurate and efficient evaluation of the change in energy barriers for carbon diffusion in ferrite under heterogeneous stress is introduced. This method, called Linear Combination of Stress States, is based on the knowledge of the effects of simple stresses (uniaxial or shear) on these diffusion barriers. Then, it is assumed that the change in energy barriers under a complex stress can be expressed as a linear combination of these already known simple stress effects. The modifications of energy barriers by either uniaxial traction/compression and shear stress are determined by means of atomistic simulations with the Climbing Image-Nudge Elastic Band method and are stored as a set of functions. The results of this method are compared to the predictions of anisotropic elasticity theory. It is shown that, linear anisotropic elasticity fails to predict the correct energy barrier variation with stress (especially with shear stress) whereas the proposed method provides correct energy barrier variation for stresses up to ∼3 GPa. This study provides a basis for the development of multiscale models of diffusion under non-uniform stress.

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

Switching times of nanoscale FePt: Finite size effects on the linear reversal mechanism

The linear reversal mechanism in FePt grains ranging from 2.316 nm to 5.404 nm has been simulated using atomistic spin dynamics, parametrized from ab-initio calculations. The Curie temperature and the critical temperature (T{sup *}), at which the linear reversal mechanism occurs, are observed to decrease with system size whilst the temperature window T{sup *}<T<T{sub C} increases. The reversal paths close to the Curie temperature have been calculated, showing that for decreasing system size the reversal path becomes more elliptic at lower temperatures, consistent with the decrease in the Curie temperature arising from finite size effects. Calculations of the minimum pulse duration show faster switching in small grains and are qualitatively described by the Landau-Lifshitz-Bloch equation with finite size atomistic parameterization, which suggests that multiscale modeling of FePt down to a grain size of ≈3.5 nm is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Hexagonal boron nitride and water interaction parameters

The study of hexagonal boron nitride (hBN) in microfluidic and nanofluidic applications at the atomic level requires accurate force field parameters to describe the water-hBN interaction. In this work, we begin with benchmark quality first principles quantum Monte Carlo calculations on the interaction energy between water and hBN, which are used to validate random phase approximation (RPA) calculations. We then proceed with RPA to derive force field parameters, which are used to simulate water contact angle on bulk hBN, attaining a value within the experimental uncertainties. This paper demonstrates that end-to-end multiscale modeling, starting at detailed many-body quantum mechanics and ending with macroscopic properties, with the approximations controlled along the way, is feasible for these systems.

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

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

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

Intrinsic anion diffusivity in lead halide perovskites is facilitated by a soft lattice

Significance Understanding and controlling the ion migration in a solid-state semiconductor is essential for device performance in modern photonics, electronics, and optoelectronics. Direct microstructural observation and atomic-scale mechanism of the intrinsic lattice diffusion are still highly challenging. Here, we report the direct visualization of halide anion interdiffusion in CsPbCl 3 –CsPbBr 3 single crystalline nanowire heterojunctions via nanoscale imaging, allowing for the direct measuring of intrinsic lattice diffusivities. The intrinsic halide diffusivities were found to be several orders of magnitudes lower than those reported in polycrystalline thin films. Together with multiscale modeling, we further demonstrate how anionic diffusion in lead halides is determined by low barriers to vacancy hopping and facilitated by the soft perovskite lattice.

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

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗
DOE OSTI2022

Mesoscale structural gradients in human tooth enamel

The outstanding mechanical and chemical properties of dental enamel emerge from its complex hierarchical architecture. An accurate, detailed multiscale model of the structure and composition of enamel is important for understanding lesion formation in tooth decay (dental caries), enamel development (amelogenesis) and associated pathologies (e.g., amelogenesis imperfecta or molar hypomineralization), and minimally invasive dentistry. Although features at length scales smaller than 100 nm (individual crystallites) and greater than 50 µm (multiple rods) are well understood, competing field of view and sampling considerations have hindered exploration of mesoscale features, i.e., at the level of single enamel rods and the interrod enamel (1 to 10 µm). Here, we combine synchrotron X-ray diffraction at submicrometer resolution, analysis of crystallite orientation distribution, and unsupervised machine learning to show that crystallographic parameters differ between rod head and rod tail/interrod enamel. This variation strongly suggests that crystallites in different microarchitectural domains also differ in their composition. Thus, we use a dilute linear model to predict the concentrations of minority ions in hydroxylapatite (Mg 2+ and CO 3 2- /Na + ) that plausibly explain the observed lattice parameter variations. While differences within samples are highly significant and of similar magnitude, absolute values and the sign of the effect for some crystallographic parameters show interindividual variation that warrants further investigation. In conclusion, by revealing additional complexity at the rod/interrod level of human enamel and leaving open the possibility of modulation across larger length scales, these results inform future investigations into mechanisms governing amelogenesis and introduce another feature to consider when modeling the mechanical and chemical performance of enamel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The case for digital twins in metal additive manufacturing

The digital twin (DT) is a relatively new concept that is finding increased acceptance in industry. A DT is generally considered as comprising a physical entity, its virtual replica, and two-way digital data communications in-between. Its primary purpose is to leverage the process intelligence captured within digital models—or usually their faster-solving surrogates—towards generating increased value from the physical entities. The surrogate models are created using machine learning based on data obtained from the field, experiments and digital models, which may be physics-based or statistics-based. Anomaly detection and correction, and diagnostic closed-loop process control are examples of how a process DT can be deployed. In the manufacturing industry, its use can achieve improvements in product quality and process productivity. Metal additive manufacturing (AM) stands to gain tremendously from the use of DTs. This is because the AM process is inherently chaotic, resulting in poor repeatability. However, a DT acting in a supervisory role can inject certainty into the process by actively keeping it within bounds through real-time control commands. Closed-loop feedforward control is achieved by observing the process through sensors that monitor critical parameters and, if there are any deviations from their respective optimal ranges, suitable corrective actions are triggered. The type of corrective action (e.g. a change in laser power or a modification to the scanning speed) and its magnitude are determined by interrogating the surrogate models. Because of their artificial intelligence (AI)-endowed predictive capabilities, which allow them to foresee a future state of the physical twin (e.g. the AM process), DTs proactively take context-sensitive preventative steps, whereas traditional closed-loop feedback control is usually reactive. Apart from assisting a build process in real-time, a DT can help with planning the build of a part by pinpointing the optimum processing window relevant to the desired outcome. Again, the surrogate models are consulted to obtain the required information. In this article, we explain how the application of DTs to the metal AM process can significantly widen its application space by making the process more repeatable (through quality assurance) and cheaper (by getting builds right the first time).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Irradiation damage concurrent challenges with RAFM and ODS steels for fusion reactor first-wall/blanket: a review

Reduced activation ferritic martensitic (RAFM) and oxide dispersion strengthened (ODS) steels are the most promising candidates for fusion first-wall/blanket (FW/B) structures. The performance of these steels will deteriorate during service due to neutron damage and transmutation-induced gases, such as helium/hydrogen, at elevated operating temperatures. Here, after highlighting the operating conditions of fusion reactor concepts and a brief overview, the main irradiation-induced degradation challenges associated with RAFM/ODS steels are discussed. Their long-term degradation scenarios such as (a) low-temperature hardening embrittlement (LTHE)—including dose-temperature dependent yield stress, tensile elongations, necking ductility, test temperature effect on hardening, Charpy impact ductile-to-brittle transition temperature and fracture toughness, (b) intermediate temperature cavity swelling, (c) the effect of helium on LTHE and cavity swelling, (d) irradiation creep and (e) tritium management issues are reviewed. The potential causes of LTHE are discussed, which highlights the need for advanced characterisation techniques. The mechanical properties, including the tensile/Charpy impact of RAFM and ODS steels, are compared to show that the current generation of ODS steels also suffers from LTHE, and shows irradiation hardening up to high temperatures of ~400 °C–500 °C. To minimise this, future ODS steel development for FW/B-specific application should target materials with a lower Cr concentration (to minimise α'), and minimise other elements that could form embrittling phases under irradiation. RAFM steel-designing activities targeting improvements in creep and LTHE are reviewed. The need to better understand the synergistic effects of helium on the thermo-mechanical properties in the entire temperature range of FW/B is highlighted. Because fusion operating conditions will be complex, including stresses due to the magnetic field, primary loads like coolant pressure, secondary loads from thermal gradients, and due to spatial variation in damage levels and gas production rates, an experimentally validated multiscale modelling approach is suggested as a pathway to future reactor component designing such as for the fusion neutron science facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Toward durable stacks: glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems

Protonic ceramic electrochemical cells (PCECs) are emerging as promising technologies for efficient energy conversion and hydrogen production because they operate at intermediate temperatures with improved efficiency and durability compared with conventional solid oxide electrochemical cells. However, the long-term reliability and commercialization of PCEC stacks remain strongly limited by the performance of sealants, which are required to maintain gas tightness, electrical insulation, and mechanical integrity under harsh thermal and chemical environments. Among various sealing approaches, glass-ceramic sealants are considered the most practical and scalable due to their excellent wettability, chemical tunability, and strong interfacial adhesion. This review provides a comprehensive overview of recent advances in glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems. The fundamental design principles of sealant compositions are first discussed, followed by recent developments in deposition methods, sintering strategies, surface treatments, and degradation monitoring techniques. Particular attention is given to the unique challenges associated with PCEC operating conditions, including hydrothermal degradation, interfacial reactions with barium-containing electrolytes, and thermal mismatch. Finally, future opportunities involving sustainable materials, multiscale modeling, additive manufacturing, and artificial intelligence-assisted sealant optimization are highlighted.

glass–ceramic sealants↗
DOE OSTI2021

Multiscale Microtubule Dynamics in Active Nematics

In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multimotor and multifilament interactions present in a dense 2D active nematic. Furthermore, they also provide quantitative data that is essential for building multiscale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Correlated anharmonicity and dynamic disorder control carrier transport in halide perovskites

Halide pervoskites are an important class of semiconducting materials that hold great promise for optoelectronic applications. In this work we investigate the relationship between vibrational anharmonicity and dynamic disorder in this class of solids. Via a multiscale model parametrized from first-principles calculations, we demonstrate that the non-Gaussian lattice motion in halide perovskites is microscopically connected to the dynamic disorder of overlap fluctuations among electronic states. This connection allows us to rationalize the emergent differences in temperature-dependent mobilities of prototypical MAPbI 3 and MAPbBr 3 compounds across structural phase transitions, in agreement with experimental findings. In conclusion, our analysis suggests that the details of vibrational anharmonicity and dynamic disorder can complement known predictors of electronic conductivity and can provide structure-property guidelines for the tuning of carrier transport characteristics in anharmonic semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature dependence of vacancy/self-interstitial recombination volumes in copper

We use molecular dynamics to calculate rate coefficients for recombination of vacancies with self-interstitial atoms (SIAs) in Cu at temperatures from 300 to 700 K. From these results, we calculate vacancy/SIA recombination volumes and find that they decrease from around 290 (where is one atomic volume) at 300 K to 160 at 500 K and above. By counting the number of distinct pathways by which a stable SIA may migrate to a site of spontaneous recombination with a nearby vacancy, we find a lower bound estimate of 168 for the recombination volume. We furthermore rationalize its temperature dependence based on differences between the activation energies for recombination and SIA migration. Furthermore, our work sheds light on the fundamental nature of the recombination volume and provides information that may be incorporated into multiscale models of radiation response in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Community Organizations: Changing the Culture in Which Research Software Is Developed and Sustained

Software is the key crosscutting technology that enables advances in mathematics, computer science, and domain-specific science and engineering to achieve robust simulations and analysis for science, engineering, and other research fields. However, software itself has not traditionally received focused attention from research communities; rather, software has evolved organically and inconsistently, with its development largely as by-products of other initiatives. Moreover, challenges in scientific software are expanding due to disruptive changes in computer hardware, increasing scale and complexity of data, and demands for more complex simulations involving multiphysics, multiscale modeling and outer-loop analysis. In recent years, community members have established a range of grass-roots organizations and projects to address these growing technical and social challenges in software productivity, quality, reproducibility, and sustainability. This article provides an overview of such groups and discusses opportunities to leverage their synergistic activities while nurturing work toward emerging software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Numbering-Up of Microscale Devices for Megawatt-Scale Supercritical Carbon Dioxide Concentrating Solar Power Receivers

Concentrated solar power (CSP) plants have the potential to reduce the consumption of nonrenewable resources and greenhouse gas emissions in electricity production. In CSP systems, a field of heliostats focuses solar radiation on a central receiver, and energy is then transferred to a thermal power plant at high temperature. However, maximum receiver surface fluxes are low (30–100 W cm -2 ) with high thermal losses, which has contributed to the limited market penetration of CSP systems. Recently, small (~4 cm 2 ), laminated micro pin-fin devices have shown potential to achieve concentrated surface fluxes over 100 W cm -2 using supercritical CO 2 as the working fluid. The present study explores the feasibility of using these microscale unit cells as building blocks for a megawatt-scale (250 MW thermal) open solar receiver through a numbering-up approach, where multiple microscale unit cell devices are connected in parallel. A multiscale model of the full-scale central receiver is developed. The model consists of interconnected unit cell and module level (i.e., multiple unit cells in parallel) submodels which predict local performance of the central receiver. Each full-scale receiver consists of 3000 micro pin-fin unit cells divided into 250 modules. The performance of three different full-scale receivers is simulated under representative operating conditions. The results show that the microscale unit cells have the potential to be numbered up to megawatt applications while providing high heat flux and thermal efficiency. At the design incident flux and surface emissivity, a global receiver efficiency of approximately 90% when heating sCO 2 from 550 °C to 650 °C at an average incident flux of 110 W cm -2 can be achieved.

14 SOLAR ENERGY↗
DOE OSTI2017

Nano- and Mesoscale Ion and Water Transport in Perfluorosulfonic-Acid Membranes

Water and aqueous cations transport along multiple length scales in perfluorosulfonic-acid membranes. Molecular interactions in hydrophilic domains dictate nanoscale resistances while the connectivity of domains controls mesoscale transport. The developed multiscale models probes the synergy of these length scales. Concentrated solution theory and electrokinetics are used to predict transport in the aqueous domains as a function of hydration. A resistor network upscales the nanoscale properties to predict effective membrane ion and water transport and their coupling. The nature of macroscopic and nanoscale properties differs drastically because the mesoscale network mediates transport. Moreover, the effective tortuosity and connectivity is not the same for water and ion transport. The methodology and findings highlight improvement opportunities for membrane performance.

25 ENERGY STORAGE↗
DOE OSTI2018

NimbleSM

NimbleSM is a Lagrangian finite-element code for solid mechanics. Its primary application is the solution of mechanics problems on nonuniform, three-dimensional meshes using either explicit transient dynamic or implicit quasi-static time integration. Additional features include multiscale modeling and the simulation of contact. The NimbleSM code base is designed for performance portability across varying hardware architectures.

Shelton, Timothy↗