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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Combining 3-Momentum and Kinetic Energy on Galilei/Newton Spacetime

Without the mass-energy equivalence available on Minkowski spacetime M, it is not possible on 4-dimensional non-relativistic Galilei/Newton spacetime G to combine 3-momentum and total mass-energy in a single tensor object. However, given a fiducial frame, it is possible to combine 3-momentum and kinetic energy into a linear form (particle) or (1,1) tensor (continuum) in a manner that exhibits increased unity of classical mechanics on flat relativistic and non-relativistic spacetimes M and G. As on M, for a material continuum on G, the first law of thermodynamics can be considered a consequence of a unified dynamical law for energy-momentum rather than an independent postulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On the origins of the continuum radiation of an underwater nanosecond pulsed discharge: an absolute-intensity optical emission spectroscopy study

The plasma properties of an underwater nanosecond pulsed discharge remain not fully understood despite being extensively studied for several decades. In this work, we focus on the continuum radiation generated in such discharges. The discharge is characterized by power measurements as well as by absolute emission spectroscopy. When observed, Stark broadenings of H α , H β and O (777 nm) are employed for electron number densities measurements. The discharge was generated by a 10 ns main voltage pulse followed by multiple secondary pulses, which last up to 4 μs after the primary pulse. It is shown that a peak power of 3.5 MW and energy of 35 mJ is coupled during the main voltage pulse. A quantitative estimation of the different possible continuum radiation sources is performed through analytical calculations. This includes emission (blackbody, free–bound and free–free bremsstrahlung radiations) and absorption (electron–ion and electron–neutral free–free inverse bremsstrahlung) mechanisms. Our results suggest that electron–neutral free–free bremsstrahlung is the principal mechanism responsible for the strong continuum radiation observed experimentally during the primary pulse. We also show that self-absorption through electron–neutral (and electron–ion) inverse bremsstrahlung plays an important role in the main discharge pulse. Further, our results indicate the non-negligible additional contribution of the H 2 continuum during the first reflected pulse which is likely ignited in bubbles generated by the first discharge pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2023

Composite Lithium Metal Structure to Mitigate Pulverization and Enable Long‐Life Batteries

In lithium metal batteries, non‐uniform stripping of lithium results in pit formation, which promotes subsequent non‐uniform, dendritic deposition. This viscous cycle leads to pulverization of lithium which promotes cell shorting or capacity degradation, symptoms further exaggerated by high electrode areal loading and lean electrolytes. Here, to address this challenge, a composite lithium metal anode is engineered that contains uniformly distributed, nanometer‐sized carbon particles. This composite lithium is shown to strip more uniformly since the growth of non‐uniform pits is intercepted by the carbon particles. This mechanism is corroborated by a continuum electrochemical model. Subsequent lithium deposition on carbon particles is also found to be more uniform than on the surface with irregular pits. Notably, the pulverization rate of composite lithium is 26 times slower than that of commercial lithium. Moreover, in a Li‐S battery with sulfurized polyacrylonitrile cathode, the use of the composite anode extends the cycle life by three times when the areal capacity is 8 mAh cm −2 . The approach of using an engineered lithium composite structure to address challenges during both stripping and plating can inform future designs of lithium metal anodes for high areal capacity operations.

high areal capacity↗
DOE OSTI2021

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

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

Peridynamic Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. Here, the computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

42 ENGINEERING↗
DOE OSTI2023

Lyman continuum leaker candidates among highly ionised, low-redshift dwarf galaxies selected from He II

Contemporary research suggests that the reionisation of the intergalactic medium (IGM) in the early Universe was predominantly realised by star-forming (proto-)galaxies (SFGs). Due to observational constraints, our knowledge on the origins of sufficient amounts of ionising Lyman continuum (LyC) photons and the mechanisms facilitating their transport into the IGM remains sparse. Recent efforts have thus focussed on the study of local analogues to these high-redshift objects. We aim to acquire a set of very low-redshift SFGs that exhibit signs of a hard radiation field being present. A subsequent analysis of their emission line properties is intended to shed light on how the conditions prevalent in these objects compare to those predicted to be present in early SFGs that are thought to be LyC emitters (LCEs). We used archival spectroscopic SDSS DR12 data to select a sample of low-redshift He II 4686 emitters and restricted it to a set of SFGs with an emission line diagnostic sensitive to the presence of an active galactic nucleus, which serves as our only selection criterion. We performed a population spectral synthesis with FADO to reconstruct these galaxies’ star-formation histories (SFHs). Utilising the spectroscopic information at hand, we constrained the predominant ionisation mechanisms in these galaxies and inferred information on interstellar medium (ISM) conditions relevant for the escape of LyC radiation. Our final sample consists of eighteen ionised, metal-poor galaxies (IMPs). These low-mass (6.2 ≤ log(M * /M ⊙ ) ≤ 8.8), low-metallicity (7.54 ≤log(O/H) + 12 ≤ 8.13) dwarf galaxies appear to be predominantly ionised by stellar sources. We find large [OIII] 5007/[OII] 3727 ratios and [SII] 6717,6731/Hα deficiencies, which provide strong indications for these galaxies to be LCEs. At least 40% of these objects are candidates for featuring cosmologically significant LyC escape fractions ≳10%. The IMPs’ SFHs exhibit strong similarities and almost all galaxies appear to contain an old (> 1 Gyr) stellar component, while also harbouring a young, two-stage (~10 Myr and < 1 Myr) starburst, which we speculate might be related to LyC escape. The properties of the compact emission line galaxies presented here align well with those observed in many local LCEs. In fact, our sample may prove as an extension to the rather small catalogue of local LCEs, as the extreme ISM conditions we find are assumed to facilitate LyC leakage. Notably, all of our eighteen candidates are significantly closer (z < 0.1) than most established LCEs. If the inferred LyC photon loss is genuine, this demonstrates that selecting SFGs from He II 4686 is a powerful selection criterion in the search for LCEs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗
DOE OSTI2025

Manufacturing of Al-Ce/Metal-Matrix-Composites (MMC) for Harsh Environments

An integrated modeling/simulation and experimental approach will be used to design and synthesize Al-Ce/MMC for use in key technologies where improved energy efficiency and performance is harsh environments is required. Modeling and simulation work including thermodynamic modeling, DFT calculations, and continuum modeling of casting processes and mechanical behavior, will be combined with lab-scale synthesis to identify potential matrix and reinforcing phase compositions for large-scale synthesis. The experimental synthesis work will rely on lab-scale casting and gas atomization to identify the optimal matrix compositions for pilot/commercial-scale synthesis is Al-Ce/MMC. Selected Al-Ce/MMC will be synthesized by commercial scale melt processing, powder metallurgy and thermomechanical processing. The corrosion resistance, mechanical properties as a function of temperature, and thermal conductivity of the composites will be characterized to validate target properties are met or exceeded. This work will focus on two selected applications; however, the design framework will have broad applicability for composites for numerous harsh service condition applications. We will have synthesized Al-Ce/MMC that exhibit better performance in harsh environments than currently used materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗
DOE OSTI2021

Revealing the structure of light pseudoscalar mesons at the electron–ion collider

The questions of how the bulk of the Universe's visible mass emerges and how it is manifest in the existence and properties of hadrons are profound, and probe the heart of strongly interacting matter. Paradoxically, the lightest pseudoscalar mesons appear to be key to a further understanding of the emergent mass and structure mechanisms. These mesons, namely, the pion and kaon, are the Nambu-Goldstone boson modes of quantum chromodynamics (QCD). Unravelling their partonic structure and the interplay between emergent and Higgs-boson mass mechanisms is a common goal of three interdependent approaches - continuum QCD phenomenology, lattice-regularised QCD, and the global analysis of parton distributions - linked to experimental measurements of hadron structure. Experimentally, the anticipated electron-ion collider will enable a revolution in our ability to study pion and kaon structures, accessed by scattering from the 'meson cloud' of the proton through the Sullivan process. With the goal of enabling a suite of measurements that can address these questions, we examine key reactions that identify the critical detector-system requirements needed to map tagged pion and kaon cross-sections over a wide range of kinematics. The excellent prospects for extracting pion structural, functional, and form-factor data are outlined, and similar prospects for kaon structures are discussed in the context of a worldwide programme. Finally, the successful completion of the programme outlined herein will deliver deep, far-reaching insights into the emergence of pions and kaons, their properties, and their role as QCD's Goldstone boson modes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2024

ACS Spotlight: Bipolar Membranes for Electrochemical Energy Conversion, Chemical Manufacturing, and Separations

Sustainable energy conversion, chemical manufacturing, and separations are central to addressing the world’s energy and environmental challenges. Electrochemical platforms stand as a cornerstone in addressing these challenges because they are low exergy and can be powered on renewable electrons. In electrochemical systems, bipolar membranes (BPMs) are emerging as a unique class of ion exchange membranes poised to revolutionize various electrochemical processes via pH control of anode and cathode chambers and in situ pH adjustment. In this Spotlight Review, we provide a comprehensive review of electrochemical platforms utilizing BPMs for energy conversion (water electrolyzers for hydrogen production, fuel cells, and flow batteries), chemical manufacturing (electrolyzers that convert carbon dioxide into value-added chemicals and nitrate into ammonia), and separations. The motivation for using BPMs, as well as their performance and durability, in electrochemical platforms are disseminated. We also discuss current challenges that impede BPM electrochemical systems from competing with state-of-the-art electrochemical systems using monopolar ion-exchange membranes (e.g., anion/hydroxide exchange membranes and cation/proton exchange membranes). Here, the review also covers molecular modeling and continuum modeling efforts to understand the basic mechanisms that govern BPM performance.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

In Situ Raman Mapping of Si Island Electrodes and Stress Modeling as a Function of Lithiation and Size

Si is known for cracking and delamination during electrochemical cycling of a battery due to the large volume change associated with Li insertion and extraction. However, it has been found experimentally that patterned Si island electrodes that are 200 nm thick and less than 7 μm wide can deform in a purely elastic manner. Inspired by this, we performed in situ Raman stress characterization of model poly-crystalline Si island electrodes using an electrochemical cell coupled with an immersion objective lens and designed for a short working distance. A 5 μm wide Si island electrode showed a parabolic stress profile during lithiation, while for a 15 μm Si island electrode, a stress plateau in the center of the electrode was observed. A continuum model with coupled electro-chemo-mechanical (ECM) physics was established to understand the stress measurement. A qualitative agreement was reached between modeling and experimental data, and the critical size effect could be explained by the Li diffusive flux as governed by competition between the Li concentration and hydrostatic stress gradients. Below the critical size, the stress gradient drives Li toward the edges, where the electrode volume is free to expand, while above the critical size, the stress plateau inhibits Li diffusion to the edge and forces destructive stress relief by cracking. Furthermore, this work represents a promising methodology for in situ characterization of ECM coupling in battery electrodes, with suggestions provided for further improvement.

25 ENERGY STORAGE↗