Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

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

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Emergence of p+ip superconductivity in two-dimensional doped Dirac systems

Searching for the $p + ip$ superconducting (SC) state has become a fascinating subject in condensed matter physics recently, as a dream application awaiting in topological quantum computation. Unfortunately, so far there is no universal principle for realizing $p + ip$ in generic solid-state systems. Here we report a theoretical discovery of a $p + ip$ SC ground state (coexisting with ferromagnetic order) in the honeycomb lattice Hubbard model in the extremely strong-coupling limit (e.g., infinite U) at low doping $(δ < 0.2)$, by using both the state-of-art Grassmann tensor product state approach and a continuum quantum field theory approach. Our discovery suggests a mechanism for the $p + ip$ SC state in generic strongly correlated systems based on spin-charge separation and the charge current–current coupling scenario, which opens a door towards experimental realization. The $p + ip$ SC state has an instability towards a potential non-Fermi liquid with a large but finite U . Nevertheless, by applying an in-plane Zeeman field, such a $p + ip$ SC state can be stabilized with finite U in a very wide range of doping. Relevant realistic materials are also proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Tomography of Atomic Nuclei

We carry out continuum quantum Monte Carlo calculations of the quantum-mechanical Wigner distribution functions of selected nuclei, up to $^{16}$O. These distributions provide a form of quantum tomography of the spatial and momentum structure of the system. They also help identify the location of high-momentum regions in atomic nuclei and provide insight into the onset of alpha clustering. Besides their intrinsic interest, these distributions will be useful for neutrino event generators, as they correlate the positions and momenta of nucleons in the initial target state. To facilitate their application, we address the need to store them compactly by developing an accurate Gaussian Process emulator that automatically preserves their normalization.

Rocco, Noemi [Valencia U., IFIC; Fermilab]↗
DOE OSTI2020

Impacts of Hurricane Disturbance on Water Quality across the Aquatic Continuum of a Blackwater River to Estuary Complex

Hurricanes cause landscape-scale disturbances that affect biogeochemical cycling and water quality in coastal ecosystems. During Hurricane Irma’s passage through northern Florida, water movements driven by wind velocities up to 105 km h-1 caused a salinity peak in an estuary/blackwater river complex. Water quality was monitored across the 15 km site to detect the magnitude and duration of disturbance. Saline water intruded 15 km inland into a freshwater portion of the river that peaked at a salinity of 2 psu. Due to the volume of precipitation from the hurricane, significant runoff of freshwater and dissolved organic matter (DOM) caused a decrease in salinity, dissolved oxygen (DO), and Chlorophyll-a concentrations while increasing turbidity and fluorescent dissolved organic matter (fDOM). The disturbance caused rapid changes observed by in-situ water quality monitors over a 3-week period, but some effects persisted for longer periods as shown by 3-month weekly water sampling. This disturbance caused shifts in DOM loading, altered salinity dynamics, and reshaped landscapes due to wind and wave surge both in upland marsh and downstream estuary. Hurricane disturbance temporarily and abruptly alters the aquatic continuum, and observations of system response can help us understand the mechanisms associated with ecosystem resilience and recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Progress on continuum modeling of lithium–sulfur batteries

While lithium–sulfur batteries are a promising next-generation chemistry devices due to their high theoretical energy density, commercialization has been slow due to low coulombic efficiency and poor cycle life. This review explores the ways in which continuum modeling contributes to the understanding of lithium–sulfur (LiS) battery mechanisms and cell-level performance through the lens of micro- and macroscale phenomena. We examine different approaches to modeling important physical phenomena such as reaction mechanisms, cathode microstructure, shuttling, nucleation and precipitation, and transport limitations. This paper also emphasizes the significance and challenge of connecting typical modeling parameters and assumptions to systems-level metrics of a standard state-of-art high performing lithium–sulfur cell. Particularly important, the considerations for high energy density cells and the areas where continuum models can facilitate better collaboration are discussed. Here, we also summarize a few selected works to highlight experimentally-driven modeling, use of electroanalytical techniques, and parameter identification approaches to enable model-based design and advanced battery management systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Twisted moiré photonic crystal enabled optical vortex generation through bound states in the continuum

The twisted stacking of two layered crystals has led to the emerging moiré physics as well as intriguing chiral phenomena such as chiral phonon and photon generation. In this work, we identified and theoretically formulated a non-trivial twist-enabled coupling mechanism in twisted bilayer photonic crystal (TBPC), which connects the bound state in the continuum (BIC) mode to the free space through the twist-enabled channel. Moreover, the radiation from TBPC hosts an optical vortex in the far field with both odd and even topological orders. We quantitatively analyzed the twist-enabled coupling between the BIC mode and other non-local modes in the photonic crystals, giving rise to radiation carrying orbital angular momentum. The optical vortex generation is robust against geometric disturbance, making TBPC a promising platform for well-defined vortex generation. As a result, TBPCs not only provide a new approach to manipulating the angular momentum of photons, but may also enable novel applications in integrated optical information processing and optical tweezers. Our work broadens the field of moiré photonics and paves the way toward the novel application of moiré physics.

42 ENGINEERING↗
DOE OSTI2023

A materials data framework and dataset for elastomeric foam impact mitigating materials

The availability of materials data for impact-mitigating materials has lagged behind applications-based data. For example, data describing on-field helmeted impacts are available, whereas material behaviors for the constituent impact-mitigating materials used in helmet designs lack open datasets. Here, we describe a new FAIR (findable, accessible, interoperable, reusable) data framework with structural and mechanical response data for one example elastic impact protection foam. The continuum-scale behavior of foams emerges from the interplay of polymer properties, internal gas, and geometric structure. This behavior is rate and temperature sensitive, therefore, describing structure-property characteristics requires data collected across several types of instruments. Data included are from structure imaging via micro-computed tomography, finite deformation mechanical measurements from universal test systems with full-field displacement and strain, and visco-thermo-elastic properties from dynamic mechanical analysis. These data facilitate modeling and design efforts in foam mechanics, e.g., homogenization, direct numerical simulation, or phenomenological fitting. The data framework is implemented using data services and software from the Materials Data Facility of the Center for Hierarchical Materials Design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Measurement of the anisotropy power spectrum of the radio synchrotron background

We present the first targeted measurement of the power spectrum of anisotropies of the radio synchrotron background, at 140 MHz, where it is the overwhelmingly dominant photon background. This measurement is important for understanding the background level of radio sky brightness, which is dominated by steep-spectrum synchrotron radiation at frequencies below ν ~ 0.5 GHz and has been measured to be significantly higher than that produced by known classes of extragalactic sources and most models of Galactic halo emission. We determine the anisotropy power spectrum on scales ranging from 2° to 0.2 arcmin with Low-Frequency Array observations of two 18-deg 2 fields – one centred on the Northern hemisphere’s coldest patch of radio sky where the Galactic contribution is smallest and the other offset from that location by 15°. We find that the anisotropy power is higher than that attributable to the distribution of point sources above 100 μJy in flux. Furthermore, this level of radio anisotropy power indicates that if it results from point sources, those sources are likely at low fluxes and incredibly numerous, and likely clustered in a specific manner.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Young but fading radio sources: searching for remnants among compact steep-spectrum radio sources

The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results of the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA (Very Large Array) CSS sample, characterized by a steep optically thin spectrum (α ≥ 1.0) and no core detection in earlier studies. Our deep multifrequency VLA, pc-scale Very Long Baseline Array, and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations allowed us to locate the core component in 10 objects. In three CSS sources, there is no clear evidence of present-time active regions, suggesting that they are likely in a remnant phase. Among sources with core detection, we find three objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in subpopulations of radio sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Diffuse sources, clustering, and the excess anisotropy of the radio synchrotron background

ABSTRACT We present the largest low frequency (120 MHz) arcminute resolution image of the radio synchrotron background (RSB) to date, and its corresponding angular power spectrum of anisotropies (APS) with angular scales ranging from 3° to 0.3 arcmin. We show that the RSB around the north celestial pole has a significant excess anisotropy power at all scales over a model of unclustered point sources based on source counts of known source classes. This anisotropy excess, which does not seem attributable to the diffuse Galactic emission, could be linked to the surface brightness excess of the RSB. To better understand the information contained within the measured APS, we model the RSB varying the brightness distribution, size, and angular clustering of potential sources. We show that the observed APS could be produced by a population of faint clustered point sources only if the clustering is extreme and the size of the Gaussian clusters is ≲1 arcmin. We also show that the observed APS could be produced by a population of faint diffuse sources with sizes ≲1 arcmin, and this is supported by features present in our image. Both of these cases would also cause an associated surface brightness excess. These classes of sources are in a parameter space not well probed by even the deepest radio surveys to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

HPC4Mfg with Samsung: Making semiconductor devices cool through HPC ab initio simulations

For decades, the semiconductor technology has followed the Moore’s law, making newer devices more powerful and energy efficient. Recently, however, it has reached a point where the performance and the energy efficiency of the device do not improve with the shrinking device size. One of the fundamental reasons of this deviation from the past trend is the interconnect resistance, which becomes larger with the shrinking size. The devices size is so small that the quantum mechanical effects can no longer be ignored and the traditional continuum simulation tools such as TCAD become inadequate. In this project, Samsung Semiconductor Inc. and Lawrence Berkeley National Laboratory has collaborated to perform first of kind device-scale ab initio simulations to optimize materials and interconnect morphology to minimize interconnect resistance. We have tested the use of LS3DF method and the PEtot_trans approach on top of the folded spectrum method (FSM) Escan code to calculate the scattering state, and to study various effects influence the interconnect conductivity. We found that, the LS3DF can be used to calculate such metallic system. On the other hand, the use of Escan code to solve the linear equation is not practical due to the slow convergence. We have implemented a Chebyshev filter technique to calculate a few hundred eigen states near the scattering state energy E, then use these eigen states as preconditioner to solve the linear equation. We have used this approach to study the different factors which affect the interconnect conductivity, including the shape, the point defect, the temperature, and the grain boundary.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanical analysis of heterogeneous materials with higher-order parameters

Abstract Even though heterogeneous porous materials are widely used in a variety of engineering and scientific fields, such as aerospace, energy-storage technology, and bio-engineering, the relationship between effective material properties of porous materials and their underlying morphology is still not fully understood. To contribute to this knowledge gap, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we use the second-order scheme that is an extension of the first-order computational homogenization framework, where a generalized continuum enables us to introduce length scale into the material constitutive law and capture both pore size and pore distribution. Through several numerical case studies with different combinations of porosity, pore shapes, and distributions, we systematically studied the relationship between the underlying morphology and effective mechanical properties. The results highlight the necessity of higher-order homogenization in understanding the mechanical properties and reveal that higher-order parameters are required to capture the role of realistic pore morphologies on effective mechanical properties. Furthermore, for specific pore shapes, higher-order parameters exhibit dominant influence over the first-order continuum.

42 ENGINEERING↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗