Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics - Plasma physics , Mathematics and Computing, Physics”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ensuring Σ s Y s = 1 in transport of species mass fractions

When transporting species mass fractions in reacting flow simulations, there are physical constraints that must be met. Unfortunately, nonlinear transport schemes such as weighted essentially non-oscillatory (WENO) schemes do not ensure that the sum of mass fractions equals 1. In detonation simulations, errors in the sum of mass fractions are observed to worsen over time when the standard WENO scheme is used. To prevent this, typically one species is forced to absorb all or most of the error in the sum of mass fractions. Here, this work presents an alternative method for correcting the WENO-interpolated mass fractions such that inert species do not change. The method is demonstrated for both argon and nitrogen-diluted hydrogen detonations in one dimension. Using the modified WENO interpolations, the error in the sum is reduced significantly. In addition, the new approach is better able to capture the physical instability expected for nitrogen-diluted detonations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Validation of electron backscattering in the MCNP6.3 code

A new electron backscattering validation suite has been developed for the MCNP® code. The calculations in this suite cover both condensed history and single-event transport methods, seven different elemental materials ranging from beryllium to uranium, and incident electron energies from 200 eV to 14.1 MeV. In general, the accuracy of the MCNP code for electron backscattering calculations is energy-dependent. For incident electron energies below 2 keV, the condensed history method is not applied while the single-event method shows poor agreement with experimental measurements due to poor accuracy of the underlying atomic data. For incident energies between 2 keV and 256 keV, the single-event method typically shows better agreement with experiments, while for incident energies above 256 keV the condensed history method gives better accuracy. Additionally, agreement with experimental data is generally worse for low- Z materials such as beryllium and carbon, which may reflect deficiencies in the data or physics underlying electron transport through light elements in the MCNP code. The new electron backscattering validation suite, along with the previously introduced electron stopping power and Lockwood energy deposition validation suites, provide MCNP users with a quantified understanding of the uncertainties in electron transport calculations. Future improvements to electron transport in the MCNP code are necessary for low energies (≤ 2 keV) and for low- Z materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three dimensional cluster analysis for atom probe tomography using Ripley’s K-function and machine learning

The size and structure of spatial molecular and atomic clustering can significantly impact material properties and is therefore important to accurately quantify. Ripley’s K-function (K(r)), a measure of spatial correlation, can be used to perform such quantification when the material system of interest can be represented as a marked point pattern. This work demonstrates how machine learning models based on K (r)-derived metrics can accurately estimate cluster size and intra-cluster density in simulated three dimensional (3D) point patterns containing spherical clusters of varying size; over 90% of model estimates for cluster size and intra-cluster density fall within 11% and 18% error of the true values, respectively. These K (r)-based size and density estimates are then applied to an experimental APT reconstruction to characterize MgZn clusters in a 7000 series aluminum alloy. Here we find that the estimates are more accurate, consistent, and robust to user interaction than estimates from the popular maximum separation algorithm. Using K (r) and machine learning to measure clustering is an accurate and repeatable way to quantify this important material attribute.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗
DOE OSTI2025

Confusion-Driven Machine Learning of Structural Phases of a Flexible, Magnetic Stockmayer Polymer

We use a semisupervised, neural-network-based machine learning technique, the confusion method, to investigate structural transitions in magnetic polymers, which we model as chains of magnetic colloidal nanoparticles characterized by dipole–dipole and Lennard-Jones interactions. As input for the neural network, we use the particle positions and magnetic dipole moments of equilibrium polymer configurations, which we generate via replica-exchange Wang–Landau simulations. We demonstrate that by measuring the classification accuracy of neural networks, we can effectively identify transition points between multiple structural phases without any prior knowledge of their existence or location. We corroborate our findings by investigating relevant conventional order parameters. Our study furthermore examines previously unexplored low-temperature regions of the phase diagram, where we find new structural transitions between highly ordered helicoidal polymer configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles-Based Study of the Decomposition of Phenol and Hydroquinone on Pt(111) Combined with Quantitative Information from XPS Spectra to Address the Impact of Coverage and Number of Hydroxyl Functional Groups

A combined first-principles-based and experimental X-ray photoelectron spectroscopy approach was used to investigate the thermal decomposition of two model biofuel compounds, phenol and hydroquinone, on Pt(111) at both low and high coverages. The DFT-based approach yields adsorption geometries and energies, activation barriers and core-level binding energy shifts for C 1s and O 1s. Increasing the coverage in the theoretical model leads to slight shifts in core-level binding energies─toward higher values for C 1s and lower values for O 1s. It also alters the energy profiles of the decomposition reaction pathway, resulting in weaker adsorption energies and changes in both reaction and activation barriers. At low temperatures, we observe a multilayer for phenol and hydroquinone upon adsorption, with desorption occurring at 200 and 270 K, respectively. Following desorption of the multilayer, decomposition proceeds via initial O–H bond scission, followed by two parallel pathways involving either C–H or C–C bond scission, whereby in the case of phenol C–H bond scission occurs first. Here, we further provide characteristic core level binding energies by theoretical calculations that are subsequently used in experimental analyses, establishing a reference database for key spectra of phenolic functionalities applicable to a range of catalytic reactions.

09 BIOMASS FUELS↗
DOE OSTI2025

How Solvent Mediates the Catalytic Fate of 2-Methoxyphenol in Hydrodeoxygenation Catalysis?

Kinetic and isotopic assessments combined with theoretical calculations shed light on the mechanistic differences in terms of the number of hydrogen addition steps and the electronic charges of reactive hydrogen in the catalytic sequence leading to the Ar−OCH 3 scission of 2-methoxyphenol occurring at various interfaces of vapor−Ru, cyclohexane−Ru, and water−Ru. At vapor−Ru and cyclohexane−Ru interfaces, the initial hydrogen adatom (H*) additions on the phenyl ring disrupt its aromaticity, as a step required to occur before kinetically relevant Ar−OCH 3 scission on an ensemble of Ru sites. The number of H* addition events prerequisite to the formation of the kinetically relevant Ar−OCH 3 scission transition state, however, is larger at the cyclohexane−Ru interface than those at the vapor−Ru and water−Ru interfaces. The preferential solvation of 2-methoxyphenol derived reaction precursors and H* weakens their adsorption on Ru sites, which leads to a decrease in the activation barrier for the third H* addition step and in turn a larger number of H* additions. At the water−Ru interface, water as a polar protic solvent promotes the Ar−OCH 3 scission by enabling a mechanistically distinct catalytic route an adjacent water molecule from the solvent layer assists with the intramolecular proton shuttling from the hydroxyl group (−OH) to the vicinal methoxy group (−OCH 3 ) of 2-methoxyphenol, evolving a highly charged transition state that is further stabilized by the water solvent layer, thereby leading to Ar−OCH 3 scission rate constants that are much higher than those at the vapor−Ru interface and, when comparing at the same H 2 chemical potential, higher Ar−OCH 3 scission turnovers. Furthermore, these distinct mechanistic details of the Ar−OCH 3 scission at various reaction interfaces, inferred from experiment and theory, illustrate how a solvent either solvates reaction intermediates and transition states or directly participates in the formation of kinetically relevant transition states, leading to changes in the free energy landscape, thus rewriting the catalytic fates in hydrodeoxygenation reactions.

09 BIOMASS FUELS↗
DOE OSTI2025

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗
DOE OSTI2021

Correspondence between neuroevolution and gradient descent

Abstract We show analytically that training a neural network by conditioned stochastic mutation or neuroevolution of its weights is equivalent, in the limit of small mutations, to gradient descent on the loss function in the presence of Gaussian white noise. Averaged over independent realizations of the learning process, neuroevolution is equivalent to gradient descent on the loss function. We use numerical simulation to show that this correspondence can be observed for finite mutations, for shallow and deep neural networks. Our results provide a connection between two families of neural-network training methods that are usually considered to be fundamentally different.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Learning in continuous action space for developing high dimensional potential energy models

Reinforcement learning (RL) approaches that combine a tree search with deep learning have found remarkable success in searching exorbitantly large, albeit discrete action spaces, as in chess, Shogi and Go. Many real-world materials discovery and design applications, however, involve multi-dimensional search problems and learning domains that have continuous action spaces. Exploring high-dimensional potential energy models of materials is an example. Traditionally, these searches are time consuming (often several years for a single bulk system) and driven by human intuition and/or expertise and more recently by global/local optimization searches that have issues with convergence and/or do not scale well with the search dimensionality. Here, in a departure from discrete action and other gradient-based approaches, we introduce a RL strategy based on decision trees that incorporates modified rewards for improved exploration, efficient sampling during playouts and a “window scaling scheme" for enhanced exploitation, to enable efficient and scalable search for continuous action space problems. Using high-dimensional artificial landscapes and control RL problems, we successfully benchmark our approach against popular global optimization schemes and state of the art policy gradient methods, respectively. We demonstrate its efficacy to parameterize potential models (physics based and high-dimensional neural networks) for 54 different elemental systems across the periodic table as well as alloys. We analyze error trends across different elements in the latent space and trace their origin to elemental structural diversity and the smoothness of the element energy surface. Broadly, our RL strategy will be applicable to many other physical science problems involving search over continuous action spaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phonon thermal Hall effect in a metallic spin ice

It has become common knowledge that phonons can generate thermal Hall effect in a wide variety of materials, although the underlying mechanism is still controversial. We study longitudinal κ xx and transverse κ xy thermal conductivity in Pr 2 Ir 2 O 7 , which is a metallic analog of spin ice. Despite the presence of mobile charge carriers, we find that both κ xx and κ xy are dominated by phonons. A T/H scaling of κ xx unambiguously reveals that longitudinal heat current is substantially impeded by resonant scattering of phonons on paramagnetic spins. Upon cooling, the resonant scattering is strongly affected by a development of spin ice correlation and κ xx deviates from the scaling in an anisotropic way with respect to field directions. Strikingly, a set of the κ xx and κ xy data clearly shows that κ xy correlates with κ xx in its response to magnetic field including a success of the T/H scaling and its failure at low temperature. This remarkable correlation provides solid evidence that an indispensable role is played by spin-phonon scattering not only for hindering the longitudinal heat conduction, but also for generating the transverse response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗