Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid dynamics”

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 19 records

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗

Impact of particle loading and phase coupling on gas–solid flow dynamics: A case study of a two-phase, gas–solid flow in an annular pipe

The present study is devoted to a two-phase, gas–solid flow in an annular pipe (hollow cylinder) at an elevated pressure of 15 bars and moderate Reynolds number of circa 6000. The influence of the particle loading, the interaction between the phases, and turbulence dispersion on the flow dynamics is systematically studied by means of computational fluid dynamics simulations, employing the Ansys FLUENT commercial package. The cases with a particle volumetric fraction of 1.2% are referred to as “high particle loading,” and those with 0.13% are denoted as “low particle loading.” The following cases are investigated: (1) pure gas flow; (2) low particle loading two-phase flow with one-way coupling and with turbulence dispersion; (3) low particle loading two-phase flow with two-way coupling but without turbulence dispersion; (4) low particle loading two-phase flow with two-way coupling and with turbulence dispersion; (5) high particle loading two-phase flow with one-way coupling and with turbulence dispersion; (6) high particle loading two-phase flow with two-way coupling but without turbulence dispersion; and (7) high particle loading two-phase flow with two-way coupling and with turbulence dispersion. Additionally, the boundary layer is found to grow without fluctuations of the turbulent kinetic energy (TKE) for cases 1, 2, and 5. For case 4, the TKE fluctuations have been identified, although they appear to be less substantial than those in cases 6 and 7. The authors attribute the semi-chaotic nature of the TKE fluctuations to the particle loading and two-way coupling. In addition, the onset and development of the flow instability have been observed at a random axial distance in cases 4, 6, and 7. Such instability is also attributed to the two-way coupling with turbulence dispersion in the flow. It is concluded that the particle loading, one-way, or two-way coupling between the phases, and the turbulence dispersion models significantly influence the development of the flow dynamics with the same inlet and boundary conditions. Consequently, it is not a trivial question, which result a user should trust. The present computational results inspire to perform verification as well as experimental validation of the simulations, so the simulation results can subsequently be used with confidence for design analysis.

42 ENGINEERING↗

Imaging solid–electrolyte interphase dynamics using operando reflection interference microscopy

The quality of the solid–electrolyte interphase is crucial for the performance of most battery chemistries, but its formation dynamics during operation are not well understood due to a lack of reliable operando characterization techniques. Herein, we report a dynamic, non-invasive, operando reflection interference microscope to enable the real-time imaging of the solid–electrolyte interphase during its formation and evolution processes with high sensitivity. The stratified structure of the solid–electrolyte interphase formed during four distinct steps includes the emergence of a permanent inner inorganic layer enriched in LiF, a transient assembly of an interfacial electrified double layer and a consequent emergence of a temporary outer organic-rich layer whose presence is reversible with electrochemical cycling. Reflection interference microscope imaging reveals an inverse correlation between the thicknesses of two interphasial subcomponents, implying that the permanent inorganic-rich inner layer dictates the organic-rich outer layer formation and lithium nucleation. In conclusion, the real-time visualization of solid–electrolyte interphase dynamics provides a powerful tool for the rational design of battery interphases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Ab initio ultrafast spin dynamics in solids

Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and are capable of nanosecond to millisecond simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe, and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent electronic spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over nanoseconds at any temperatures and doping levels. We first apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes dominant at low temperatures for spin relaxation in n-type GaAs. We further examine ultrafast dynamics in novel spin-valleytronic materials: monolayer and bilayer WSe 2 with realistic defects. We find that spin relaxation is highly sensitive to local symmetry and chemical bonds around defects. For the bilayer WSe 2 , we identify the scattering pathways in ultrafast dynamics and determine relevant dynamical properties, essential to its utilization of unique spin-valley-layer locking effects. In conclusion, our work provides a predictive computational platform for spin dynamics in solids, which has potential for designing new materials ideal for spintronics and quantum information technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical study of implosion instability mitigation in magnetically driven solid liner dynamic screw pinches

Magnetized liner inertial fusion experiments on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that compromise integrity of the imploding cylindrical liner, limiting achievable fusion fuel conditions and ultimately reducing magneto-inertial fusion target performance. Dynamic screw pinches (DSP) provide a method to reduce MRTI in-flight via application of magnetic field line tension to the imploding liner outer surface. In contrast with z-pinches that drive implosions with an azimuthal magnetic field, dynamic screw pinches enforce an additional axial drive magnetic field component, making the overall drive magnetic field helical. As the liner implodes, cumulative MRTI development is reduced by dynamically shifting the orientation of the fastest growing instability modes. Three-dimensional magnetohydrodynamic simulations show that the DSP mechanism effectively stabilizes initially solid cylindrical liner implosions driven by Z-scale current pulses, indicating that MRTI mitigation increases with the ratio of axial to azimuthal drive magnetic field components (i.e., the drive field ratio). We also performed a spectral analysis of the simulated imploding density distributions, extracting wavelength and pitch angle of the simulated MRTI structures to study their dynamics during the implosion. Simulations of liners initially perturbed with drive-field-aligned sinusoidal structures indicate that MRTI mitigation in DSP implosions decreases with perturbation wavelength, once again suggestive of magnetic field line tension effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

X-ray thermal diffuse scattering as a texture-robust temperature diagnostic for dynamically compressed solids

We present a model of x-ray thermal diffuse scattering (TDS) from a cubic polycrystal with an arbitrary crystallographic texture, based on the classic approach of Warren [B. E. Warren, Acta Crystallogr. 6, 803 (1953)]. We compare the predictions of our model with femtosecond x-ray diffraction patterns gathered from ambient and dynamically compressed rolled copper foils obtained at the High Energy Density instrument of the European X-Ray Free-Electron Laser facility and find that the texture-aware TDS model yields more accurate results than does the conventional powder model owed to Warren. Nevertheless, we further show: with sufficient angular detector coverage, the TDS signal is largely unchanged by sample orientation and in all cases strongly resembles the signal from a perfectly random powder; shot-to-shot fluctuations in the TDS signal resulting from grain-sampling statistics are at the percent level, in stark contrast to the fluctuations in the Bragg-peak intensities (which are over an order of magnitude greater); and TDS is largely unchanged even following texture evolution caused by compression-induced plastic deformation. We conclude that TDS is robust against texture variation, making it a flexible temperature diagnostic applicable just as well to off-the-shelf commercial foils as to ideal powders.

Crystal lattices↗

Off-harmonic optical probing of high intensity laser plasma expansion dynamics in solid density hydrogen jets

Due to the non-linear nature of relativistic laser induced plasma processes, the development of laser-plasma accelerators requires precise numerical modeling. Especially high intensity laser-solid interactions are sensitive to the temporal laser rising edge and the predictive capability of simulations suffers from incomplete information on the plasma state at the onset of the relativistic interaction. Experimental diagnostics utilizing ultra-fast optical backlighters can help to ease this challenge by providing temporally resolved inside into the plasma density evolution. We present the successful implementation of an off-harmonic optical probe laser setup to investigate the interaction of a high-intensity laser at 5.4 x 10 21 W/cm 2 peak intensity with a solid-density cylindrical cryogenic hydrogen jet target of 5μm diameter as a target test bed. The temporal synchronization of pump and probe laser, spectral filtering and spectrally resolved data of the parasitic plasma self-emission are discussed. The probing technique mitigates detector saturation by self-emission and allowed to record a temporal scan of shadowgraphy data revealing details of the target ionization and expansion dynamics that were so far not accessible for the given laser intensity. Plasma expansion speeds of up to (2.3 ± 0.4)x 10 7 m/s followed by full target transparency at 1.4ps after the high intensity laser peak are observed. A three dimensional particle-in-cell simulation initiated with the diagnosed target pre-expansion at -0.2ps and post processed by ray tracing simulations supports the experimental observations and demonstrates the capability of time resolved optical diagnostics to provide quantitative input and feedback to the numerical treatment within the time frame of the relativistic laser-plasma interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement and Simulation of Ultra-Low-Energy Ion–Solid Interaction Dynamics

Ion implantation is a key capability for the semiconductor industry. As devices shrink, novel materials enter the manufacturing line, and quantum technologies transition to being more mainstream. Traditional implantation methods fall short in terms of energy, ion species, and positional precision. Here, we demonstrate 1 keV focused ion beam Au implantation into Si and validate the results via atom probe tomography. We show the Au implant depth at 1 keV is 0.8 nm and that identical results for low-energy ion implants can be achieved by either lowering the column voltage or decelerating ions using bias while maintaining a sub-micron beam focus. We compare our experimental results to static calculations using SRIM and dynamic calculations using binary collision approximation codes TRIDYN and IMSIL. A large discrepancy between the static and dynamic simulation is found, which is due to lattice enrichment with high-stopping-power Au and surface sputtering. Additionally, we demonstrate how model details are particularly important to the simulation of these low-energy heavy-ion implantations. Finally, we discuss how our results pave a way towards much lower implantation energies while maintaining high spatial resolution.

36 MATERIALS SCIENCE↗

Electrolyte Role in SEI Evolution at Si in the Pre-lithiation Stage vs the Post-lithiation Stage

The formation and evolution of the dynamic solid electrolyte interphase (SEI) at the Si anode/electrolyte interface are yet to be completely understood to solve irreversible capacity loss and increase battery cycle life. Herein, the evolution of SEI and its dynamic properties at the Si anode/electrolyte interface are investigated in two electrolyte systems, a 1.2 M LiPF6 in EC: EMC 3:7 (wt%) electrolyte (referred to as Gen2) and a 1.2 M LiTFSI in EC: EMC 3:7 (wt%) electrolyte (referred to as LiTFSI). Two lithiation stages are studied: the pre-lithiation (pre-Li) SEI stage and the post-lithiation (post-Li) stage. Findings reveal at the pre-Li, SEI formation starts at an early potential and contributes to the greater mass gain in the Si/Gen2, and it is dominated by the formation of a non-uniform F- and P-rich layer in Si/Gen2, in contrast to a homogeneous F- and C-containing layer at the Si/LiTFSI interphase. The initially formed SEI in LiTFSI further benefits the charge transfer kinetics. At the post-Li stage, a more substantial SEI evolution is observed at Si/LiTFSI. This paper offers a foundational understanding of the SEI dynamic evolution with electrolyte dependence. Findings from this report offer important insights into solving the complex SEI stability issues on Si.

25 ENERGY STORAGE↗

Experimental validation of the mechanistic scale-up methodology of gas–solid spouted beds using radioactive particle tracking (RPT)

The very high-temperature reactors (VHTRs) are highly ranked among candidates of Generation IV of nuclear reactors due to their high efficiency, safety, the resistance to proliferation, and reliability. The VHTRs are preferentially fueled by Tristructural-isometric (TRISO) coated fuel particles which has fuel kernels of fissile material coated by four coating layers: a porous buffer pyrolysis carbon layer (buffer PyC), an inner dense pyrolysis carbon layer (IPyC), a silicon carbide layer (SiC) and an outer dense pyrolysis carbon layer (OPyC). The heart of the operation and safety of the VHTRs significantly depends on the reliability of the coating layers of TRISO particles to retain metallic and gaseous fission products within the particles. The technique used for coating TRISO particles are gas-solids spouted beds via chemical vapor deposition (CVD). Fabrication of high-quality low-defect TRISO fuel particles fuel at larger scale spouted beds is required to support the commercialization of the VHTRs. In this work, our new developed mechanistic scale-up methodology of gas-solids spouted beds based on matching the radial profile of gas-holdup has been demonstrated and validated using radioactive particle tracking (RPT). Two spouted beds of small and larger scales were used in the study. Three sets of conditions were carried out which include the conditions of the reference case in the large scale, conditions that provide similar gas holdup radial profile to that of the reference case and conditions that provided dissimilar gas holdup radial in the small-scale spouted beds. The results confirm the validation of the scale-up methodology in terms of the dimensionless values of the spout diameter, cumulative probability distribution of the solids particles penetration into the spout, fraction of cycle time in each region of the bed, the radial profiles of the dimensionless values of the root-mean-square particle velocities and solids eddy diffusivity. Finally, the results further advance the knowledge and understanding of the gas-solids spouted beds provide deeper insight into their solids dynamics and presenting important benchmarking data for validating computational fluid dynamics codes and models. At last, procedures are established for the implementation of the new scale-up methodology.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗