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

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗
DOE OSTI2024

5-2428: Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior [Slides]

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameter uncertainty and better integrate friction constitutive laws for faults to improve our ability to predict and mitigate seismic hazard at Utah FORGE and future EGS sites. We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, using lab measurements of friction parameters, herein referred to as “THMC+E” models. These simulations will enable exploration of various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗
DOE OSTI2022

What Is the Right Level of Activation of a High-Spin {FeNO} 7 Complex to Enable Direct N–N Coupling? Mechanistic Insight into Flavodiiron NO Reductases

Flavodiiron nitric oxide reductases (FNORs), found in pathogenic bacteria, are capable of reducing nitric oxide (NO) to nitrous oxide (N 2 O) to detoxify NO released by the human immune system. Previously, we reported the first FNOR model system that mediates direct NO reduction (Dong, H. T.; et al. J. Am. Chem. Soc. 2018, 140, 13429-13440), but no intermediate of the reaction could be characterized. Here, we present a new set of model complexes that, depending on the ligand substitution, can either mediate direct NO reduction or stabilize a highly activated high-spin (hs) {FeNO} 7 complex, the first intermediate of the reaction. The precursors, [{Fe II (MPA-(RPhO) 2 )} 2 ] (1, R = H and 2, R = t Bu, Me), were prepared first and fully characterized. Complex 1 (without steric protection) directly reduces NO to N 2 O almost quantitatively, which constitutes only the second example of this reaction in model systems. Contrarily, the reaction of sterically protected 2 with NO forms the stable mononitrosyl complex 3, which shows one of the lowest N-O stretching frequencies (1689 cm -1 ) observed so far for a mononuclear hs-{FeNO} 7 complex. Here this study confirms that an N-O stretch & LE;1700 cm -1 represents the appropriate level of activation of the FeNO unit to enable direct NO reduction. The higher activation level of these hs-{FeNO} 7 complexes required for NO reduction compared to those formed in FNORs emphasizes the importance of hydrogen bonding residues in the active sites of FNORs to activate the bound NO ligands for direct N-N coupling and N2O formation. The implications of these results for FNORs are further discussed.

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

Micro-structural features and material properties impact on adhesive metal joints via computational modeling and machine learning

The quality of structural bonding in practical applications depends on various factors arising from materials, pre-processing conditions, and manufacturing. Understanding how these factors influence bonding performance and determining their relative importance are of significant interest. Thus, this study evaluates the effects of microstructural features and material properties on the structural strength of adhesively-bonded metal joints at the submillimeter scale, utilizing a combination of Finite Element Modeling (FEM) and Machine Learning (ML) with Gradient Boosting Regression (GBR). The microstructural features include adhesive thickness, internal voids within the adhesive, adherend-adhesive interfacial voids, void size and volume fraction, and surface roughness. The material properties include the constitutive behavior of the adhesive, as well as the adherend-adhesive interfacial strength and fracture energy. The changes in structural strength and morphologies of the bonded metal structures with respect to different microstructural features and material properties were clarified by FEM. By further leveraging ML-GBR, the sequence of importance of these factors affecting bonding performance across various scenarios was summarized. This work provides valuable insights into the development of improved structural bonding for adhesive joints in industries such as automotive , aerospace, and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2025

Arctic Black Carbon Aerosol Deposition Study North Slope of Alaska 2020- UHSAS Measurements

Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. Within this data set, we measured fluxes of size resolved particles with an ultrahigh sensitivity aerosol spectrometer (UHSAS) at the ARM facility’s North Slope Alaska site. Measurements were made at the meteorological tower between 9 Sept and 25 Oct 2021. This site is a coastal tundra location and received snow during the project.

Aerosol concentration↗
DOE OSTI2025

Arctic Black Carbon Aerosol Deposition Study North Slope of Alaska 2020- SP2 Measurements

Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. Within this data set, we measured refractory black carbon with a single particle soot photometer (SP2) at the ARM facility’s North Slope of Alaska site. Measurements were made at the meteorological tower between 9 Sept and 27 Oct 2021. This site is a coastal tundra location, and received snow during the project.

Black carbon mass concentration↗
DOE OSTI2025

Arctic Black Carbon Aerosol Deposition Study North Slope of Alaska 2020- Sonic Anemometer Measurements

Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. Within this data set, we measured wind properties using a sonic anemometer. Measurements were made at the meteorological tower between 9 Sept and 25 Oct 2021. This site is a coastal tundra location, and received snow during the project.

Gill Windmaster 3D Ultrasonic Anemometer↗
DOE OSTI2026

Arctic Black Carbon Aerosol Deposition Study North Slope of Alaska 2020- ACSM Measurements

Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. In support of black carbon (SP2) and particles (UHSAS) flux measurements at the meteorological tower, we measured non-refractory chemically-speciated aerosol composition with an Aerosol Chemical Speciation Monitor at the the NSA Central Facility (C1). ACSM measurements were made between 9 Sept and 26 Oct 2021. This site is a coastal tundra location, and received snow during the project.

Chl_ACSM↗
DOE OSTI2026

Arctic Black Carbon Aerosol Deposition Study North Slope of Alaska 2020- ACSM Measurements

Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. In support of black carbon (SP2) and particles (UHSAS) flux measurements at the meteorological tower, we measured non-refractory chemically-speciated aerosol composition with an Aerosol Chemical Speciation Monitor at the the NSA Central Facility (C1). ACSM measurements were made between 9 Sept and 26 Oct 2021. This site is a coastal tundra location, and received snow during the project.

Aerosol Chemical Speciation Monitor↗
DOE OSTI2025

Online and Offline Data Quality Monitoring for the Mu2e Calorimeter

This thesis presents the design, implementation, and validation of a calorimeter Data Quality Monitoring (DQM) toolchain for the Mu2e experiment at Fermilab. Mu2e searches for charged lepton flavor violation via coherent muon-to-electron conversion in the field of an aluminum nucleus, $\mu^- Al \rightarrow e^-Al$, a process whose observation would constitute clear evidence of physics beyond the Standard Model. Achieving target sensitivity requires stringent control of detector performance and data integrity during acquisition, as subtle issues in readout configuration, data formatting, or electronics behavior can compromise reconstruction and bias downstream analyzes. To address these challenges, this work develops a multi-layer DQM approach spanning both raw data validation and reconstructed digi-level diagnostics. At the low level, a fragment analysis component performs word- and bit-field decoding of calorimeter readout blocks, enabling sanity checks of the expected structure and producing detailed error and integrity statistics useful for commissioning and troubleshooting. At the digi level, the CaloDigiDQM analyzer is implemented within the art framework and transforms each CaloDigiCollection into a structured hierarchy of ROOT histograms designed for fast drill-down diagnostics. The module generates coherent monitoring views at global, disk, board, and channel granularity, including occupancy, waveform-derived features (baseline, RMS, peak amplitude and position), and left-right sensor consistency metrics. Detector-aware channel-to-electronics mapping is performed through the conditions system (CaloDAQMap), ensuring that diagnostics remain aligned with hardware identifiers used in operations. For end-to-end testing without reliance on live DAQ data, a synthetic CaloDigi producer is developed to generate realistic waveforms with controlled noise and pulse shapes. The resulting system supports both offline ROOT-file production and online operation, including optional histogram streaming through otsdaq via ots::HistoSender. This toolchain provides a practical and scalable foundation for calorimeter commissioning and stable data collection, enabling early detection of anomalies and reducing operational risk for Mu2e.

Vakulenko, Mark [Drew U.] (ORCID:0009000276197818)↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2025

Noble dark matter: Surprising elusiveness of dark baryons

Dark matter could be a baryonic composite of stronglyc upled constituents transforming under SU⁡(2) 𝐿 . We classify the SU⁡(2) 𝐿 representations of baryons in a class of simple confining dark sectors and find that the lightest state can be a pure singlet or a singlet that mixes with other neutral components of SU⁡(2) 𝐿 representations, which strongly suppresses the dark matter candidate’s interactions with the Standard Model. We focus on models with a confining SU⁡(𝑁 𝑐 ) and heavy dark quarks constituting vectorlike 𝑁 𝑓 -plet of SU⁢(2) 𝐿 . For benchmark 𝑁 𝑐 and 𝑁 𝑓 , we calculate baryon mass spectra, incorporating electroweak gauge boson exchange in the nonrelativistic quark model, and demonstrate that above TeV mass scales, dark matter is dominantly a singlet state. The combination of this singlet nature with the recently discovered ℋ-parity results in an inert state analogous to noble gases, hence we coin the term noble dark matter. Our results can be understood in the nonrelativistic effective theory that treats the dark baryons as elementary states, where we find singlets accompanying triplets, 5-plets, or more exotic representations. This generalization of WIMP-like theories is more difficult to find or rule out than dark matter models that include only a single SU⁡(2) 𝐿 multiplet (such as a wino), motivating new searches in colliders and a reanalysis of direct and indirect detection prospects in astrophysical observations.

dark matter direct detection↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
DOE OSTI2020

Multiscale Thermomechanical Damage Model with Internal State Variables for Ceramic Matrix Composites

A length scale-dependent matrix damage model has been developed to simulate the response of ceramic matrix composites (CMCs). This model captures damage initiation and propagation in the brittle matrix material by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress-strain constitutive relationships at each length scale. Material degradation due to matrix cracking is captured by a damage variable determined using fracture mechanics and the self-consistent scheme. An additional state variable is introduced to capture the nucleation and growth of matrix porosity, a key mechanism contributing to matrix nonlinearity. Porosity occurs as a result of material diffusion around grain boundaries and the growth of the porosity state variable is related to the material entropy dissipation and the volumetric strain. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide, unidirectional (UD) carbon fiber/silicon carbide matrix (C/SiC) CMC, and five-harness satin (5HS) woven C/SiC CMC.

Skinner, Travis D.↗
DOE OSTI2023

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗