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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics: A DOE-NIH Joint Workshop Executive Summary

The DOE-NIH Joint Workshop on Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics, held on September 27, 2024, brought together experts from government, academia, and industry to address critical challenges in radionuclide production and clinical translation. Here, the workshop emphasized interdisciplinary collaboration, particularly between the Department of Energy (DOE) and the National Institutes of Health (NIH), to strengthen the domestic isotope supply, streamline regulatory pathways, and further integrate computational tools into radiopharmaceutical therapy (RPT). Key discussions explored the role of AI-driven modeling, machine learning, and digital twin technologies in optimizing dosimetry, dynamically personalizing treatments, and reducing time to clinical adoption. Advances in predictive computational modeling were highlighted as essential for improving radionuclide yield, purity, and synthesis efficiency. Regulatory considerations and equitable access were central themes, with participants advocating for harmonized global standards, adaptive trial designs, and expanded infrastructure for clinical implementation. DOE computational and production infrastructure was emphasized. Future priorities identified include increased investment in radionuclide production infrastructure, expanded workforce development in radiopharmaceutical sciences and computational modeling, and the creation of robust public-private partnerships. The workshop concluded that continued strategic collaboration and sustained resources will be vital for advancing next-generation radiotheranostics, ensuring safe and effective therapies accessible to all patients.

digital twins↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Toward “optimal” integration of terrestrial biosphere models

Multimodel ensembles (MME) are commonplace in Earth system modeling. Here we perform MME integration using a 10-member ensemble of terrestrial biosphere models (TBMs) from the Multiscale synthesis and Terrestrial Model Intercomparison Project (MsTMIP). We contrast optimal (skill based for present-day carbon cycling) versus naive (one model-one vote) integration. MsTMIP optimal and naive mean land sink strength estimates (-1.16 versus -1.15 Pg C per annum respectively) are statistically indistinguishable. This holds also for grid cell values and extends to gross uptake, biomass, and net ecosystem productivity. TBM skill is similarly indistinguishable. The added complexity of skill-based integration does not materially change MME values. This suggests that carbon metabolism has predictability limits and/or that all models and references are misspecified. Finally, resolving this issue requires addressing specific uncertainty types (initial conditions, structure, and references) and a change in model development paradigms currently dominant in the TBM community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Hybrid Simulation of Low Temperature Plasmas: A Brief Tutorial

Modeling and simulation of low temperature plasmas is widely accepted as a useful tool from both physics and plasma reactor design points of view. Of the modeling approaches, hybrid models strive to combine the predictability of kinetic models with the lower computational burden of fluid models. This work presents a brief tutorial of hybrid plasma modeling and representative simulation results using hybrid models. Whenever possible, these results are compared to kinetic and/or fluid simulations as well as experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spatial patterns in occupancy and density of larval lampreys in freshwater habitats restored to a Stage 0 condition

Abstract We examined occupancy and density of larval lampreys ( Entosphenus tridentatus and Lampetra spp.) in two rivers in Oregon (USA) restored to a Stage 0 condition 1–5 years prior, using a multiscale occupancy model and a zero‐inflated Poisson mixture model. We sampled lampreys using backpack electrofishing in randomly distributed, paired, 1‐m 2 quadrats and recorded environmental data. Probabilities of occupancy and density were higher when water velocity was low, the substrate was noncompacted, and sediment was dominated by fines (<4 mm). At mean water depth (0.34 m) and velocity (0.09 m/s), estimated densities in occupied quadrats were 4.8 lampreys/m 2 (95%: 3.4–6.9) when the substrate was compacted, and fines were not dominant, and 21.1 lampreys/m 2 (95%: 17.7–25.3) when the substrate was noncompacted and fines were dominant. Probabilities of detecting occupancy in a 1‐m 2 quadrat sampled by backpack electrofishing were 0.76 (95%: 0.64–0.87) when captured after visual observation and 0.80 (95%: 0.71–0.88) with blind sweeps (i.e., constantly moving the net regardless of observation). The probability of capturing a single lamprey in a quadrat sampled by blind sweeps was 0.32 (95%: 0.27–0.37). Sampling in paired 1‐m 2 quadrats facilitated concurrent examination of patterns in occupancy and density while accounting for capture probability, which could aid temporal monitoring of restored habitats. To the best of our knowledge, this is the first study to document occupancy and estimate densities of larval lampreys in habitats that underwent valley floor restoration to Stage 0. We observed both lamprey genera within 5 years of restoration. Aquatic restoration that increases low‐velocity, noncompacted, fine sediment habitats could benefit lampreys.

Harris, Julianne E.↗
DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Mathematical modeling in cancer nanomedicine: a review

Cancer continues to be among the leading healthcare problems worldwide, and efforts continue not just to find better drugs, but also better drug delivery methods. The need for delivering cytotoxic agents selectively to cancerous cells, for improved safety and efficacy, has triggered the application of nanotechnology in medicine. This effort has provided drug delivery systems that can potentially revolutionize cancer treatment. Nanocarriers, due to their capacity for targeted drug delivery, can shift the balance of cytotoxicity from healthy to cancerous cells. The field of cancer nanomedicine has made significant progress, but challenges remain that impede its clinical translation. Several biophysical barriers to the transport of nanocarriers to the tumor exist, and a much deeper understanding of nano-bio interactions is necessary to change the status quo. Mathematical modeling has been instrumental in improving our understanding of the physicochemical and physiological underpinnings of nanomaterial behavior in biological systems. Here, we present a comprehensive review of literature on mathematical modeling works that have been and are being employed towards a better understanding of nano-bio interactions for improved tumor delivery efficacy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
DOE OSTI2026

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the dynamic evolution of atomically dispersed Cu catalyst for CO 2 electrochemical conversion using integrated XANES analysis and mechanistic studies

Direct electrochemical conversion of CO 2 to ethanol (CH 3 CH 2 OH) offers a promising strategy to lower CO 2 emission while storing energy from renewable electricity. Our recent study reported a carbon-supported atomically dispersed Cu catalyst that achieved the highest reported selectivity for CH 3 CH 2 OH formation (91%) at a relatively low potential (-0.6 V), however, the active site structure that is responsible for such high activity and selectivity has yet to be understood. Here, in this paper, we demonstrate a computational investigation combining X-ray absorption near edge structure (XANES) simulations and a mechanistic study via density functional theory (DFT) to understand the catalyst structures of this Cu catalyst during electrocatalysis and the corresponding reaction mechanisms of the key products. An integrated computational and experimental XANES analysis depicted the dynamic evolution of the catalytic site during electrocatalysis. The as-prepared, atomically dispersed Cu catalyst aggregates and forms metallic clusters/nanoparticles under electrochemical condition, which then break down to smaller oxidized clusters after electrocatalysis. The formed Cu clusters/nanoparticles showed distinct catalytic activity and selectivity as a function of particle size based on the mechanistic investigation using DFT, which is consistent with experimental observations for catalyst samples with different Cu loadings. This comprehensive study which combines experimental and computational XANES investigation, mechanistic study via DFT calculations, and experimental performance of the catalysts, provides unprecedented dynamic and mechanistic insights into the supported atomically dispersed metal catalysts for CO 2 reduction. Such strategy and details gained can further guide discovery of novel catalyst materials for CO 2 electrochemical reduction.

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

Supervised paralled-in-time algorithm for long-time Lagrangian simulations of stochastic dynamics: Application to hydrodynamics

Lagrangian particle methods based on detailed atomic and molecular and nanoscale systems. However, the maximum time step is limited by the smallest oscillation period of the fastest atomic motion, rendering long-time simulations very expensive. To resolve this bottleneck, we propose a supervised parallel-in-time algorithm for stochastic dynamics to accelerate long-time Lagrangian particle simulations. Our method is inspired by bottom-up coarse-graining projections that yield mean-field hydrodynamic behavior in the continuum limit.

Multiscale modeling, Lagrangian method, Parallel-i↗