Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗
DOE OSTI2024

Evolution and engineering of pathways for aromatic O -demethylation in Pseudomonas putida KT2440

In this study, biological conversion of lignin from biomass offers a promising strategy for sustainable production of fuels and chemicals. However, aromatic compounds derived from lignin commonly contain methoxy groups, and O-demethylation of these substrates is often a rate-limiting reaction that influences catabolic efficiency. Several enzyme families catalyze aromatic O-demethylation, but they are rarely compared in vivo to determine an optimal biocatalytic strategy. Here, two pathways for aromatic O-demethylation were compared in Pseudomonas putida KT2440. The native Rieske non-heme iron monooxygenase (VanAB) and, separately, a heterologous tetrahydrofolate-dependent demethylase (LigM) were constitutively expressed in P. putida, and the strains were optimized via adaptive laboratory evolution (ALE) with vanillate as a model substrate. All evolved strains displayed improved growth phenotypes, with the evolved strains harboring the native VanAB pathway exhibiting growth rates ~1.8x faster than those harboring the heterologous LigM pathway. Enzyme kinetics and transcriptomics studies investigated the contribution of selected mutations toward enhanced utilization of vanillate. The VanAB-overexpressing strains contained the most impactful mutations, including those in VanB, the reductase for vanillate O-demethylase, PP_3494, a global regulator of vanillate catabolism, and fghA, involved in formaldehyde detoxification. These three mutations were combined into a single strain, which exhibited approximately 5x faster vanillate consumption than the wild-type strain in the first 8 h of cultivation. Overall, this study illuminates the details of vanillate catabolism in the context of two distinct enzymatic mechanisms, yielding a platform strain for efficient O-demethylation of lignin-related aromatic compounds to value-added products.

09 BIOMASS FUELS↗
DOE OSTI2026

Multi-resolution enhancement for full-spectrum neural representations

Scientific data acquisition continues to outpace storage and analysis capabilities, making voxel-basedrepresentations increasingly intractable. Implicit neural representations (INRs) offer a promising solutionby encoding signals through coordinate-based neural networks, serving as surrogates of data, withcomputational and storage requirements scaling with network complexity rather than data dimensionality.However, smaller INRs struggle to faithfully represent multiscale structures, high-frequency informationand fine textures that constitute a large proportion of scientific measurements. We propose WIEN-INR, atheoretically guided hierarchical INR framework that distributes modelling across resolution scales andenables improved representation capacity through a novel enhancement network to recover subtle details.This multiscale architecture allows smaller networks to retain the full spatial-frequency content of thesignal as well as preserve training efficiency and lower storage cost. Evaluated on distinct raw experimentalmeasurements across scales and complexities, WIEN-INR represents a practical step towards a broaderadoption of neural representations in scientific workflows, delivering compact, robust and high-fidelityrepresentations.

Ni, Yuan [SLAC National Accelerator Laboratory (SL↗
DOE OSTI2020

New Physics from the Neutrino Portal: Early Universe Implications and Detection Prospects

Despite their ubiquity, the nature of Dark Matter (DM) and neutrinos remains mysterious. Both are examples of particles which are neutral, stable, and require physics beyond the Standard Model (SM) of particle physics. It is therefore natural to investigate theories which relate them. If DM has more than gravitational interactions, it must not couple very sizably to the Standard Mode (SM) in order to have escaped detection so far. The nature of this coupling (or “portal”) is of paramount importance in determining the optimal experimental search strategies. This proposal aims to focus on the implications of “Hidden Sector” models of new physics where the dominant interaction is from the neutrino-mixing portal. Such scenarios give rise to new neutrino self-interactions and DM-neutrino interactions with a variety of consequences. By undertaking the thorough study of neutrino portal (NP) theories of DM this proposal seeks to achieve the following specific objectives: (1) Expand and deepen the theoretical foundations of self-interacting neutrinos in the early universe with application to sterile neutrino DM, (2) establish theoretical and phenomenological bases for DM-modified neutrino oscillations, (3) establish sensitivity projections for neutrino self-interactions from future IceCube data. The proposed research will intellectually contribute to the physics literature by thoroughly characterizing the consequences of “hidden sector” models incorporating DM and neutrinos. Theoretical consistency and a complete description of their impact on the early universe will constitute one component of this work. As a result of the additional structure in the neutrino sector the detectability of various aspects of neutrino physics will be impacted. This includes modifications to oscillation experiments, ultra-high-energy neutrino telescopes (e.g. IceCube, KM3NET), and Cosmic Microwave Background data. The success of this research proposal will depend on the successful integration of these complementary probes of new physics in the proper context of theoretically motivated scenarios determined from their impact on the early universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Improving ProtoDUNE pion cross-section measurements with NuGraph Michel-electron tagging

Understanding hadron-argon interactions is essential for precise neutrino energy reconstruction and final-state interaction modeling in liquid-argon time projection chamber (LArTPC) experiments such as DUNE. In particular, pion absorption and charge-exchange processes constitute significant sources of systematic uncertainty in neutrino oscillation measurements. ProtoDUNE-SP, a large-scale LArTPC prototype operated at the CERN Neutrino Platform and exposed to charged-particle test beams in the few-GeV range, enables direct measurements of these processes. This work focuses on the measurement of differential cross sections for pion absorption and charge exchange using the 2 GeV/c pion beam data from the ProtoDUNE-SP run. A key component of this analysis is the identification of Michel electrons from $\pi \rightarrow \mu \rightarrow e$ decay chains, which helps separate different interaction topologies and improves background rejection. Michel electron identification will also assist in reliably calibrating the electromagnetic response in ProtoDUNE-SP data and for the future DUNE detectors. In this analysis, we apply NuGraph to identify Michel electrons. NuGraph is a graph neural network that models detector hits as nodes connected by spatial and temporal edges for particle and topology classification in LArTPC detectors. We first benchmark NuGraph’s Michel electron classification performance using ICEBERG data, a small-scale LArTPC prototype used for DUNE electronics and reconstruction development, and then transfer the approach to ProtoDUNE-SP. This poster presents the analysis strategy, NuGraph-based classification studies, and discusses how these developments are expected to improve the pion cross-section measurement.

Razafinime, Soamasina Herilala [Cincinnati U.] (OR↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multi‐scale simulation of rock compaction through breakage models with microstructure evolution

Abstract. Regional subsidence due to fluid depletion includes the interaction among multiple physical processes. Specifically, rock compaction is governed by coupled hydro-mechanical feedbacks involving fluid flow, effective stress change and pore collapse. Although poroelastic models are often used to explain the delay between depletion and subsidence, recent evidence indicates that inelastic effects could alter the rock microstructure, thus exacerbating coupling effects. Here, a constitutive law built within the framework of Breakage Mechanics is proposed to account for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. Furthermore, it is embedded into a 1-D hydromechanical coupled finite element analysis (FEA) to explore the effects of micro-structure rearrangement on the development of reservoir compaction. Numerical examples with the proposed model are compared with simulations under constant hydraulic conductivity to illustrate the model capability to capture the non-linear processes of reservoir compaction induced by fluid depletion.

Buscarnera, Giuseppe↗
DOE OSTI2021

Experimental modal analysis of SSR1 cryomodule for numerical model tuning and validation

The present activity aims to test the Cryomodule SSR1 to validate the corresponding mixed multibody-finite element numerical model and to determine the range of frequencies in which the first natural modes of the structure lie. The experimental activity consisted of a hammer test of the cryomodule at an early assembly stage at which the internal parts, constituting the cold mass, are accessible and, thus, many accelerometers can be attached to the beam string. The comparison between the experimental results obtained for the sub-assembly and the numerical results obtained for the corresponding sub-model, derived from the previously developed complete model, allowed to assess the validity of the numerical model and to get feedback on the effectiveness of the placement of the sensors in capturing the dynamic response of the system during transportation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Probing Nuclear Effects with Transverse Kinematic Imbalance in Muon-neutrino Induced Charged-Current $π^0$ Production on Argon with the MicroBooNE Detector

Neutrino-nucleus cross-section measurements are needed to improve interaction modeling and to enable precision neutrino oscillation measurements in upcoming experiments such as the Deep Underground Neutrino Experiment (DUNE), Hyper-Kamiokande, and the Short-Baseline Neutrino program. Baryon-resonance neutrino interactions constitute a dominant contribution near the peak of the DUNE neutrino energy spectrum. We present the first measurement of muon neutrino charged-current resonance-like interactions on argon using transverse kinematic imbalance variables with the MicroBooNE detector. These observables are highly sensitive to the modeling of final-state interactions. This measurement probes kinematic imbalances using the reconstructed momenta of the muon, leading proton, and neutral pion. A comprehensive characterization of the $π^0$-proton final state is presented; however, none of the models considered are able to simultaneously reproduce all measured observables.

Abratenko, P. [Tufts U. (main)]↗
DOE OSTI2022

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials

This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.

42 ENGINEERING↗
DOE OSTI2025

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2021

Looking forward to lepton-flavor-violating ALPs

We assess the status of past and future experiments on lepton flavor violating (LFV) muon and tau decays into a light, invisible, axion-like particle (ALP), a . We propose a new experimental setup for MEG II, the MEGII-fwd , with a forward calorimeter placed downstream from the muon stopping target. Searching for μ → ea decays MEGII-fwd is maximally sensitive to LFV ALPs, if these have nonzero couplings to right-handed leptons. The experimental set-up suppresses the (left-handed) Standard Model background in the forward direction by controlling the polarization purity of the muon beam. The reach of MEGII-fwd is compared with the present constraints, the reach of Mu3e and the Belle-II reach from τ → ℓa decays. We show that a dedicated experimental campaign for LFV muon decays into ALPs at MEG II and Mu3e will be able to probe the ALP parameter space in an unexplored region well beyond the existing astrophysical constraints. We study the implications of these searches for representative LFV ALP models, where the presence of a light ALP is motivated by neutrino masses, the strong CP problem and/or the SM flavor puzzle. To this extent we discuss the majoron in low-scale seesaw setups and introduce the LFV QCD axion, the LFV axiflavon and the leptonic familon, paying particular attention to the cases where the LFV ALPs constitute cold dark matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Clarifying the formation of equiaxed grains and microstructural refinement in the additive manufacturing of Ti-Cu

Controlling microstructural evolution in metallic additive manufacturing (AM) is difficult, especially in producing refined as-built grains instead of coarse, directional grains. Traditional solutions involve adding inoculants to AM feedstocks, but titanium (Ti) alloys cannot employ this approach without producing detrimental secondary phases. Ti-Cu (Ti-copper) alloys offer a solution through constitutional supercooling and/or solid state thermal cycling under AM conditions. This work analyzes a compositionally graded directed energy deposition (DED) Ti-Cu build, single-melt laser tracks, and dilatometric heat treatments to evaluate if, when, and by what mechanism(s) microstructural refinement occurs. Refinement by inoculation of unmelted powder particles was also considered. Constitutional supercooling produced no net microstructural refinement as any equiaxed dendrites which form are remelted with new deposition. This finding agreed with solidification modeling of powder bed fusion-laser beam (PBF-LB) and DED builds. Solid state thermal cycling refined microstructures only during ex-situ dilatometric heat treatments, suggesting build parameter optimization is needed to achieve refinement in-situ. Accidental heterogeneous nucleation on unmelted Ti powder, originating from the different thermophysical properties of Ti and Cu, provided the most significant microstructural refinement. This work systematically assesses the microstructural refinement mechanisms of Ti-Cu in AM builds and offers insights into microstructural control in eutectoid alloys.

36 MATERIALS SCIENCE↗