Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical description and modelling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

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

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Onset of fast magnetic reconnection and particle energization in laboratory and space plasmas

The onset of magnetic reconnection in space, astrophysical and laboratory plasmas is reviewed discussing results from theory, numerical simulations and observations. After a brief introduction on magnetic reconnection and approach to the question of onset, in this paper we first discuss recent theoretical models and numerical simulations, followed by observations of reconnection and its effects in space and astrophysical plasmas from satellites and ground-based detectors, as well as measurements of reconnection in laboratory plasma experiments. Mechanisms allowing reconnection spanning from collisional resistivity to kinetic effects as well as partial ionization are described, providing a description valid over a wide range of plasma parameters, and therefore applicable in principle to many different astrophysical and laboratory environments. Finally, we summarize the implications of reconnection onset physics for plasma dynamics throughout the Universe and illustrate how capturing the dynamics correctly is important to understanding particle acceleration. The goal of this review is to give a view on the present status of this topic and future interesting investigations, offering a unified approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Influence of Hyperfine Coupling Strain on Two-Dimensional ESEEM Spectra from I = 1/2 Nuclei

HYSCORE, a variant of two-dimensional ESEEM spectroscopy, is currently one of the most effective tools in high-resolution solid-state EPR spectroscopy. Our systematic studies of hydrogen bonds with semiquinones in 15 N-labeled proteins have previously provided us with an extensive collection of 15 N HYSCORE spectra produced by N–H…O nitrogen donors. Examination of these spectra has indicated in several cases lineshape distortions not described by available theoretical models. We were able to explain all observed artificial phenomena in the 15 N spectra of semiquinones considering a single mechanism—strain of the isotropic hyperfine interaction. In this review article, we introduce the HYSCORE experiment and the theoretical background for spectra of I = 1/2 nuclei in orientation-disordered samples. We then present a description of the observed lineshape distortions and an explanation of the hyperfine strain’s influence on the experimental spectra. In the final part of the article, we discuss (i) similar lineshape distortions found in previously published spectra from 15 N, 29 Si, 31 P, and 57 Fe nuclei that can be explained using the approach applied for the analysis of the 15 N spectra of semiquinones, and (ii) a recently published analysis of 31P spectra with a different type of lineshape distortion based on a more sophisticated method, where the strength of anisotropic coupling is correlated with the sign and value of the isotropic constant in different conformations of a Cu(II) complex.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Model-free estimation of completeness, uncertainties, and outliers in atomistic machine learning using information theory

Abstract An accurate description of information is relevant for a range of problems in atomistic machine learning (ML), such as crafting training sets, performing uncertainty quantification (UQ), or extracting physical insights from large datasets. However, atomistic ML often relies on unsupervised learning or model predictions to analyze information contents from simulation or training data. Here, we introduce a theoretical framework that provides a rigorous, model-free tool to quantify information contents in atomistic simulations. We demonstrate that the information entropy of a distribution of atom-centered environments explains known heuristics in ML potential developments, from training set sizes to dataset optimality. Using this tool, we propose a model-free UQ method that reliably predicts epistemic uncertainty and detects out-of-distribution samples, including rare events in systems such as nucleation. This method provides a general tool for data-driven atomistic modeling and combines efforts in ML, simulations, and physical explainability.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structures in the heavy-ion fusion excitation function at and above the Coulomb barrier

Contrary to descriptions from coupled-channels or other model calculations, heavy-ion fusion excitation functions are not smooth near and above the Coulomb barrier. There appear to be weak but noticeable oscillations or structures within the excitation functions that can be observed clearly in representation $d (σ E)/d E$ and in the comparison with theoretical calculations $σ (E) – σ_{th}(E)$. Importantly, a rather similar phenomenon has been studied before and can be explained in light symmetric systems as the influence of the centrifugal barrier penetration only on even angular momentum, but it cannot be extended to heavier, more asymmetric fusion systems. A more thorough investigation may be required to investigate this newly indicated behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Toward a quantitative description of solvation structure: a framework for differential solution scattering measurements

Appreciating that the role of the solute–solvent and other outer-sphere interactions is essential for understanding chemistry and chemical dynamics in solution, experimental approaches are needed to address the structural consequences of these interactions, complementing condensed-matter simulations and coarse-grained theories. High-energy X-ray scattering (HEXS) combined with pair distribution function analysis presents the opportunity to probe these structures directly and to develop quantitative, atomistic models of molecular systems in situ in the solution phase. However, at concentrations relevant to solution-phase chemistry, the total scattering signal is dominated by the bulk solvent, prompting researchers to adopt a differential approach to eliminate this unwanted background. Though similar approaches are well established in quantitative structural studies of macromolecules in solution by small- and wide-angle X-ray scattering (SAXS/WAXS), analogous studies in the HEXS regime—where sub-ångström spatial resolution is achieved—remain underdeveloped, in part due to the lack of a rigorous theoretical description of the experiment. To address this, herein we develop a framework for differential solution scattering experiments conducted at high energies, which includes concepts of the solvent-excluded volume introduced to describe SAXS/WAXS data, as well as concepts from the time-resolved X-ray scattering community. Our theory is supported by numerical simulations and experiment and paves the way for establishing quantitative methods to determine the atomic structures of small molecules in solution with resolution approaching that of crystallography.

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

Experimental and theoretical description of the process of contact laser surgery with a titanium-doped optothermal fibre converter

In an in vitro experiment simulating a surgeon’s actions in the process of contact laser surgery of soft biological tissue, the dependences of the temperature of a titanium-doped optothermal fibre converter (TOTFC) and the depths of coagulation and ablation of biological tissue on the average radiation power of a diode laser with a wavelength of 980 nm and on the speed of the converter movement along biological tissue are obtained. The structural, optical, and thermophysical models of TOTFC are discussed, as well as the thermophysical model of the interaction of a laser-heated converter with biological tissue, which takes into account the temperature dependences of the basic thermophysical parameters of the converter and biological tissue, as well as the contribution of the thickness h{sub int} of the water vapour layer between the converter and biotissue. It is shown that the proposed model allows describing the result of contact laser surgery of soft biotissue with TOTFC adequately to the experiment. (laser biophotonics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy

The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Applied-field magnetic structure and spectroscopy shifts of the effective spin-$\frac{1}{2}$ 𝑋⁢𝑌-like magnet Li 2 ⁢CoCl 4

Insulators containing chains of magnetic transition metal cations provide platforms for probing spin-$\frac{1}{2}$ dynamics and quantum critical behavior. Li 2 ⁢ CoCl 4 contains edge-sharing CoCl 6 octahedra that form chains along the crystallographic 𝑐 axis and orders antiferromagnetically at zero field, but questions remain about its applied-field magnetic structure and the Co 2+ spin state. Here, in this study, we show with neutron diffraction on a polycrystalline sample how the antialigned chains of cobalt moments begin to transition to a ferromagnetic state above 1.6 T. Further, using magnetic resonance absorption measurements and noninteracting spin models, we reveal the strongly anisotropic nature of the Co 2+ ion's 𝑋⁢𝑌-like magnetic behavior (𝑔 ∥ = 2.77 and 𝑔 ⊥ = 5.23) and its 𝐽 = $\frac{1}{2}$ ground state. We therefore supply the magnetic structures and anisotropic description needed to explore the dynamics of the field-driven magnetic phases, laying the foundation for further experimental and theoretical studies.

Riedel, Zachary W. [Univ. of Illinois at Urbana-Ch↗
DOE OSTI2022

Dielectric strength of noble and quenched gases for high pressure time projection chambers

Abstract Dielectric breakdown strength is one of the critical performance metrics for pure gases and gas mixtures used in large, high pressure gas time projection chambers. In this paper we experimentally study dielectric breakdown strengths of several important time projection chamber working gases and gas-phase insulators over the pressure range 100 mbar to 10 bar, and gap sizes ranging from 0.1 to 10 mm. Gases characterized include argon, xenon, CO $$_2$$ 2 , CF $$_4$$ 4 , and mixtures 90-10 argon-CH $$_4$$ 4 , 90-10 argon-CO $$_2$$ 2 and 99-1 argon-CF $$_4$$ 4 . We develop a theoretical model for high voltage breakdown based on microphysical simulations that use electron swarm Monte Carlo results as input to Townsend- and Meek-like discharge criteria. This model is shown to be highly predictive at high pressure, out-performing traditional Paschen–Townsend and Meek–Raether models significantly. At lower pressure-times-distance, the Townsend-like model is an excellent description for noble gases whereas the Meek-like model provides a highly accurate prediction for insulating gases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Exploring nonequilibrium phases of photo-doped Mott insulators with generalized Gibbs ensembles

Many experiments show that strong excitations of correlated quantum materials can cause non-thermal phases without equilibrium analogues. Understanding the origin and properties of these nonequilibrium states has been challenging due to the limitations of theoretical methods for nonequilibrium strongly correlated systems. In this work, we introduce a generalized Gibbs ensemble description that enables a systematic analysis of the long-time behavior of photo-doped states in Mott insulators based on equilibrium methods. We demonstrate the power of the method by mapping out the nonequilibrium phase diagram of the one-dimensional extended Hubbard model, which features η -pairing and charge density wave phases in a wide photo-doping range. We furthermore clarify that the peculiar kinematics of photo-doped carriers, and the interaction between them, play an essential role in the formation of these non-thermal phases. Our results establish a new path for the systematic analysis of nonequilibrium strongly correlated systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A numerically exact description of ultrafast vibrational decoherence in vibration-coupled electron transfer

Broadband pump–probe spectroscopy has been widely used to measure vibrational decoherence associated with the reaction coordinate in photoinduced ultrafast vibration-coupled electron transfer (VCET) reactions. These experiments provide insight into the interplay of intramolecular coordinates along the reaction coordinate. However, a general theoretical foundation for analyzing, and even for explaining rigorously, these data is lacking. In this work, we study vibrational decoherence in a model VCET reaction using the nearly exact time-dependent density matrix renormalization group simulation method. We explore how analyzing the density matrix with quantum information measures can help elucidate the evolution of vibrational coherence in simulations of dynamics. We examine how vibrational coherence is affected by electron transfer on the timescale of approximately 100 femtoseconds. Our results suggest that electron transfer, in the nonadiabatic model, changes the vibrational equilibrium position abruptly—an example of a “quantum quench” event. This explains the concomitant vibrational decoherence. We find that abrupt vibrational decoherence can be mitigated by wavepacket motion occurring on the timescale of the electron transfer.

Science & Technology - Other Topics↗
DOE OSTI2026

Extending near-axis equilibria in $\mathtt{DESC}$

The near-axis description of optimised stellarator fields has proven to be a powerful tool both for the design and understanding of this magnetic confinement concept. The description consists of an asymptotic model of the equilibrium in the distance from its centremost axis, and is thus only approximate. Any practical application therefore requires the eventual construction of a global equilibrium. This paper presents a novel way of constructing global equilibria using the DESC code that guarantees the correct asymptotic behaviour imposed by a given near-axis construction. The theoretical underpinnings of this construction are carefully presented, and benchmarking examples provided. This opens the door to an efficient coupling of the near-axis framework and that of global equilibria for future optimisation efforts.

fusion plasma↗
DOE OSTI2023

An alternative viewpoint on the nuclear structure towards 100 Sn: Lifetime measurements in 105 Sn

This work aims at presenting an alternative approach to the long standing problem of the B(E2) values in Sn isotopes in the vicinity of the N=Z double-magic nucleus 100 Sn, until now predominantly measured with relativistic and intermediate-energy Coulomb excitation reactions. The direct measurement of the lifetime of low-lying excited states in odd-even Sn isotopes provides a new and precise guidance for the theoretical description of the nuclear structure in this region. Lifetime measurements have been performed in 105 Sn for the first time with the coincidence Recoil Distance Doppler Shift technique. The lifetime results for the $7/2^+_1$ first excited state and the $11/2^+_1$ state, $2^+ (^{104}$Sn) Ⓧ $v1g_{7/2}$ multiplet member, are discussed in comparison with state-of-the-art shell model and mean field calculations, highlighting the crucial contribution of proton excitation across the core of 100 Sn. The reduced transition probability B(E2) of the $11/2^+_1$ core-coupled state points out an enhanced staggering with respect to the B(E2; $2^+_1$ → $0^+_1$) in the even-mass 104 Sn and 106 Sn isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗