Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space phase model”

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 937 records · Page 52

Measurement and feedback-driven adaptive dynamics in the classical and quantum kicked top

In classical dynamical systems, stochastic feedback can stabilize otherwise unstable periodic orbits, giving rise to distinct controlled and uncontrolled phases as the rate of control application is varied. In this work, we apply these control protocols in classical, semiclassical, and quantum regimes to the kicked top, a paradigmatic model of quantum chaos. The quantum kicked top, modeled as the dynamics of a spin-S object, naturally interpolates between these regimes with the spin size S acting as an effective Planck constant. We show that the dynamics of the kicked top in classical, semiclassical, and fully quantum limits can all be controlled using stochastic feedback protocols. Comparing the full quantum dynamics to a truncated Wigner approximation that captures quantum noise but neglects interference beyond the Ehrenfest time, we find that low-moment observables are largely accounted for semiclassically, while the remaining discrepancy in higher moments is consistent with contributions from interference and possibly nonlinearities in rare trajectories that explore the compact phase space. We also find rapid purification in the numerics studied for all rates of control considered, suggesting that control quenches the top's ability to encode a qubit of quantum information even in the uncontrolled phase.

Prasad, Mahaveer [Tata Inst.; Singapore U. Tech. D↗

Fractionalized Excitations Revealed by Entanglement Entropy

Fractionalized excitations develop in many unusual many-body states such as quantum spin liquids, disordered phases that cannot be described using any local order parameter. Also, because these exotic excitations correspond to emergent degrees of freedom, how to probe them and establish their existence is a long-standing challenge. We present a general procedure to reveal the fractionalized excitations using real-space entanglement entropy in critical spin liquids that are particularly relevant to experiments. Furthermore, we show how to use the entanglement entropy to construct the corresponding spinon Fermi surface. Our work defines a new pathway to establish and characterize exotic excitations in novel quantum phases of matter.

2-dimensional systems↗

Measurement of nuclear effects in neutrino-argon interactions using generalized kinematic imbalance variables with the MicroBooNE detector

We present a set of new generalized kinematic imbalance variables that can be measured in neutrino scattering. These variables extend previous measurements of kinematic imbalance on the transverse plane and are more sensitive to modeling of nuclear effects. We demonstrate the enhanced power of these variables using simulation and then use the MicroBooNE detector to measure them for the first time. We report flux-integrated single- and double-differential measurements of charged-current muon neutrino scattering on argon using a topology with one muon and one proton in the final state as a function of these novel kinematic imbalance variables. These measurements allow us to demonstrate that the treatment of charged current quasielastic interactions in genie version 2 is inadequate to describe data. As a result, they reveal tensions with more modern generator predictions particularly in regions of phase space where final state interactions are important.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Synthesis and magnetic properties of the Shastry-Sutherland family R 2 Be 2 SiO 7 ( R = Nd , Sm , Gd-Yb )

Compounds forming the quasi-two-dimensional Shastry-Sutherland lattice (SSL) have attracted significant experimental and theoretical attention in the field of frustrated magnetism. This is primarily due to their realization of an exactly soluble J 1 –J 2 orthogonal dimer model capable of hosting magnetic order, dimer singlet, and plaquette singlet phases in zero applied field and their complex magnetic phase diagrams with fractional magnetization plateaus and possible superfluid and supersolid phases found between the plateau states. The discovery and characterization of SSL compounds based on rare-earth magnetic ions provide a direct route to study the stability and properties of these exotic magnetic phases in systems with a variety of different magnetic anisotropies. In this paper, we discuss the synthesis and magnetic characterization of polycrystalline samples of the R 2 ⁢Be 2 ⁢SiO 7 family, where R=Nd, Sm, and Gd-Yb. All family members crystallize in the space group $P\bar{⁢4}⁢2_1⁢m$ (113) and show no signs of long-range magnetic order above 2 K, except for R=Tb which orders antiferromagnetically at 2.6 K.

36 MATERIALS SCIENCE↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Disorder Dynamics in Battery Nanoparticles During Phase Transitions Revealed by Operando Single-Particle Diffraction

Structural and ion-ordering phase transitions limit the viability of sodium-ion intercalation materials in grid scale battery storage by reducing their lifetime. However, the combination of phenomena in nanoparticulate electrodes creates complex behavior that is difficult to investigate, especially on the single-nanoparticle scale under operating conditions. In this work, operando single-particle X-ray diffraction (oSP-XRD) is used to observe single-particle rotation, interlayer spacing, and layer misorientation in a functional sodium-ion battery. oSP-XRD is applied to Na 2/3 [Ni 1/3 Mn 2/3 ]O 2 , an archetypal P2-type sodium-ion-positive electrode material with the notorious P2-O2 phase transition induced by sodium (de)intercalation. Here, it is found that during sodium extraction, the misorientation of crystalline layers inside individual particles increases before the layers suddenly align just prior to the P2-O2 transition. The increase in the long-range order coincides with an additional voltage plateau signifying a phase transition prior to the P2-O2 transition. To explain the layer alignment, a model for the phase evolution is proposed that includes a transition from localized to correlated Jahn-Teller distortions. The model is anticipated to guide further characterization and engineering of sodium-ion intercalation materials with P2-O2 type transitions. oSP-XRD, therefore, opens a powerful avenue for revealing complex phase behavior in heterogeneous nanoparticulate systems.

25 ENERGY STORAGE↗

Evidence for WW/WZ vector boson scattering in the decay channel ℓνqq produced in association with two jets in proton-proton collisions at s = 13 TeV

Evidence is reported for electroweak (EW) vector boson scattering in the decay channel $\ell\nu$qq of two weak vector bosons WV (V = W or Z), produced in association with two parton jets. The search uses a data set of proton-proton collisions at 13 TeV collected with the CMS detector during 2016-2018 with an integrated luminosity of 138 fb$^{-1}$. Events are selected requiring one lepton (electron or muon), moderate missing transverse momentum, two jets with a large pseudorapidity separation and a large dijet invariant mass, and a signature consistent with the hadronic decay of a W/Z boson. The cross section is computed in a fiducial phase space defined at parton level requiring all parton transverse momenta $p_\mathrm{T}$$\gt$ 10 GeV and at least one pair of outgoing partons with invariant mass $m_\mathrm{qq}$$\gt$ 100 GeV. The measured and expected EW WV production cross sections are 1.90 $^{+0.53}_{-0.46}$ pb and 2.23 $^{+0.08}_{-0.11}$ (scale) $\pm$ 0.05 (PDF) pb, respectively, where PDF is the parton distribution function. The observed EW signal strength is $m_\mathrm{EW}$ = 0.85 $\pm$ 0.12 (stat) $^{+0.19}_{-0.17}$ (syst), corresponding to a signal significance of 4.4 standard deviations with 5.1 expected, and it is measured keeping the quantum chromodynamics (QCD) associated diboson production fixed to the standard model prediction. This is the first evidence of vector boson scattering in the $\ell\nu$qq decay channel at LHC. The simultaneous measurement of the EW and QCD associated diboson production agrees with the standard model prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fingerprinting shock-induced deformations via diffraction

Abstract During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

36 MATERIALS SCIENCE↗

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗

A combined machine learning and density functional theory study of binary Ti-Nb and Ti-Zr alloys: Stability and Young’s modulus

The multicomponent Ti alloys, specifically the -phase, have experienced a strong growth over the last decades, due to their outstanding properties of ultra-high strength and low Young’s modulus. These properties play a significant role in many aerospace and biomedical applications. Selection and optimization of multicomponent alloys is challenging due to the vast chemical and compositional space. Here we investigate the use of machine learning techniques informed by density functional calculations to guide the selection of Nb- and Zr-based Ti binary alloys. From the cubic structures obtained from high throughput calculations and literature, we identify several structures with Young’s moduli below 40 GPa. The multivariant decision tree methods provide efficient surrogate models to identify structure variables have high influences on the energetic stability and Young’s modulus. Furthermore, we implement a workflow of incorporating DFT provided results and machine learning method to explore the chemical and composition space of other binary and multicomponent alloys, to eventually accelerate the material design via taking advantages of identified key variables.

36 MATERIALS SCIENCE↗

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗

Slowly decaying zero mode in a weakly nonintegrable boundary impurity model

The transverse field Ising model (TFIM) on the half-infinite chain possesses an edge zero mode. This work considers an impurity model: TFIM perturbed by a boundary integrability-breaking interaction. For sufficiently large transverse field, but in the ordered phase of the TFIM, the zero mode is observed to decay. The decay is qualitatively different from zero modes where the integrability-breaking interactions are nonzero all along the chain. It is shown that for the impurity model, the zero mode decays by relaxing to a nonlocal quasiconserved operator, the latter being exactly conserved when the opposite edge of the chain has no noncommuting perturbations so as to ensure perfect degeneracy of the spectrum. In the thermodynamic limit, the quasiconserved operator vanishes, and a regime is identified where the decay of the zero mode obeys Fermi's golden rule. Furthermore, a toy model for the decay is constructed in Krylov space and it is highlighted how Fermi's golden rule may be recovered from this toy model.

1-dimensional systems↗

Electronic Structure, Spectroscopy and Correlation Effects in Novel Materials

This project was targeted on electronic structure, spectroscopic studies, and correlation effects in a variety of novel materials of great current interest. Spectroscopies resolved highly in momentum, energy or spatial dimensions are playing a key role in unraveling the nature of the ground state and excitation properties in wide classes of novel materials. The seminal insights thus obtained are of critical importance not only for answering some of the fundamental questions facing condensed matter physics and materials science today, but also for understanding and thus helping to design and develop new materials with desirable properties, which will continue to be the key to the survival of mankind and its energy needs as a technological society long into the future. However, spectroscopies do not provide a direct map of electronic states, but act as a very complex ‘filter’ or ‘mapping’ of the underlying spectrum. This connection between the electronic states and measured spectra—called the ‘matrix element effect’—is in general an extremely complex function of the phase space of the experiment (e.g. energy/polarization of photons in photoemission), presenting both a challenge and an opportunity. So motivated, this project pursued techniques for realistic treatment of electronic spectra of a wide variety of materials, which served as a prelude to formulating and implementing methodologies for making direct connection with various spectroscopies such as ARPES, STS/STM, and inelastic light scattering. Specific systems are topological materials, including 2D ultrathin films beyond graphene, novel superconductors, and nanocrystals, among others. Our goal was to exploit the strengths of various spectroscopies to piece together the most complete picture of electronic states in systems of current interest, enabling direct and sharpened confrontations with theoretical models, and also to help advance the reach of various spectroscopies.

36 MATERIALS SCIENCE↗

The phase of the gravitational path integral

The gravitational path integral on S 2 × S 2 can be interpreted either as evaluating a contribution to the norm of the Hartle-Hawking wavefunction conditional on spatial S 1 × S 2 topology, or the pair creation rate of black holes in de Sitter. Both interpretations are distinguished at the quantum level. The former requires the path integral to be real and the latter to be imaginary. We develop a formalism to efficiently compute the phase of the gravitational path integral on Einstein spaces. We apply it to a broad class of spacetimes and in particular S 2 × S D−2 , finding it to be real and positive. We generalize some of the analysis to cases with charge and rotation.

AdS-CFT correspondence↗

Inverse design of hypoeutectoid pearlite steel microstructures using a deep learning and genetic algorithm optimization framework

Goal-oriented microstructure design in metallic materials is a challenging task due to complex structure-property relationships. Traditional experimental and computational approaches are time-intensive and economically inefficient, limiting their applicability for large-scale design space exploration. Here, in this work, we propose an end-to-end framework that integrates deep learning models with genetic optimization to design microstructures with targeted mechanical properties. Deep learning models enable accurate forward design, while their integration with genetic optimization enables efficient inverse design within a few hours, compared to days or weeks using conventional finite element simulations. The framework combines experimental characterization and finite element modeling to analyze the influence of microstructural features on the mechanical behavior of hypoeutectoid steels. Data from both experiments and simulations are used to train the deep learning models. To demonstrate its effectiveness, we apply the framework to 0.63% carbon steel with proeutectoid ferrite and pearlite phases, commonly used in industrial applications. In this study, 2D microstructures were used for modeling, selected primarily for computational efficiency and to establish proof of concept. The framework successfully optimizes microstructures for targeted yield strength, ultimate strength, and stress concentration factors while significantly reducing computational time. Beyond hypoeutectoid steels, this scalable framework can be extended to other material systems and integrated with additive manufacturing, offering an efficient approach for accelerating microstructure design for specific engineering applications.

ConvLSTM↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here, in this work, we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗