Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cascading from $\mathscr{N}$ = 2 supersymmetric Yang–Mills theory to confinement and chiral symmetry breaking in adjoint QCD

We argue that adjoint QCD in 3 + 1 dimensions, with any SU(N) gauge group and two Weyl fermion flavors (i.e. one adjoint Dirac fermion), confines and spontaneously breaks its chiral symmetries via the condensation of a fermion bilinear. We flow to this theory from pure $\mathscr{N}$ = 2 SUSY Yang–Mills theory with the same gauge group, by giving a SUSY-breaking mass M to the scalars in the $\mathscr{N}$ = 2 vector multiplet. This flow can be analyzed rigorously at small M, where it leads to a deconfined vacuum at the origin of the $\mathscr{N}$ = 2 Coulomb branch. The analysis can be extended to all M using an Abelian dual description that arises from the N multi-monopole points of the $\mathscr{N}$ = 2 theory. At each such point, there are N −1 hypermultiplet Higgs fields h$^{i=1,2}_m$, which are SU(2) R doublets. We provide a detailed study of the phase diagram as a function of M, by analyzing the semi-classical phases of the dual using a combination of analytic and numerical techniques. The result is a cascade of first-order phase transitions, along which the Higgs fields h i m successively turn on, and which interpolates between the Coulomb branch at small M, where all h$^{i}_m$ = 0, and a maximal Higgs branch, where all h$^{i}_m$ ≠ 0, at sufficiently large M. We show that this maximal Higgs branch precisely matches the confining and chiral symmetry breaking phase of two-flavor adjoint QCD, including its broken and unbroken symmetries, its massless spectrum, and the expected large-N scaling of various observables. The spontaneous breaking pattern SU(2) R → U(1) R , consistent with the Vafa–Witten theorem, is ensured by an intricate alignment mechanism for the h$^{i}_m$ in the dual, and leads to a CP 1 sigma model of increasing radius along the cascade.

D’Hoker, Eric [Univ. of California, Los Angeles, C↗
DOE OSTI2022

Simulating self-powered neutron detector responses to infer burnup-induced power distribution perturbations in next-generation light water reactors

Understanding how 3D power distribution will be monitored throughout reactor core volumetric space in next-generation nuclear power reactors is crucial to the design, deployment, and licensing of these reactors. Although numerous techniques exist for 3D power distribution monitoring based on the response of both in situ and ex situ sensors currently implemented or proposed for use in the US reactor fleet, crucial details about these techniques are often unclear. The publicly available documentation does not include information such as how well these techniques are characterized and optimized in their implementations and the levels of uncertainty in the inferred 3D power distribution. The work described herein investigated a recently developed 3D power distribution inferencing method as applied to two next-generation reactor simulations: (1) the NuScale small modular reactor design and (2) the Westinghouse AP1000 design, both of which contain in-core strings of vanadium self-powered neutron detectors (SPNDs). This investigation considered a range of SPND string sensor densities, as well as a range of 3D power distribution axial segment sizes. In this work, SPND response simulation is informed by neutron flux calculations in representative homogenized cores. For the different sensor densities and power distribution axial segment sizes in these simulations, the average solution error, solver iterations, and run time were tracked to parameterize the sensor-core configuration.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Optimizing the location and configuration of disaster resilience hubs under transportation and electric power network failures

Natural disasters often result in failures of transportation network components and blackouts that imperil the wellbeing of vulnerable populations. In response to these events, resilience hubs have been proposed as a pre-disaster planning strategy to improve access to critical services. This paper introduces an optimization-based approach to locate and configure electric power-generating resilience hubs considering the possibility of failures in transportation and electric power systems. The model's objective is to identify hub locations and configurations that maximize transportation accessibility to the hubs and maximize the satisfaction of basic energy needs through hub-generated electric power. Besides a budget constraint, the model accounts for limits on the levels of hub energy generation vis-à-vis community energy demands, and on the transportation network distance of communities to hubs. Three heuristics are presented for the proposed planning problem. The first heuristic is a genetic algorithm (GA) with problem-specific solution generation procedures. The other two heuristics implement greedy search techniques. Numerical experiments were conducted, using data from rural Puerto Rico, to illustrate the application of the proposed model and heuristics, and examine their performance. In the numerical experiments, the GA heuristic found better solutions than the greedy heuristics. Additionally, design solutions consisting of spatially dispersed hubs with low energy generation capacity were better than solutions with spatially concentrated high-capacity hubs. Lastly, across a wide range of hub demand scenarios, only a small number of candidate hub locations consistently ranked among the best locations for establishing a hub.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Real-Time Time-Dependent Electronic Structure Theory

Real-time electronic structure methods provide an unprecedented view of electron dynamics and ultrafast spectroscopy on the atto- and femtosecond time scale with vast potential to yield new insights into the electronic behavior of molecules and materials. Here, we discuss the fundamental theory underlying various real-time electronic structure methods as well as advantages and disadvantages of each. We give an overview of the numerical techniques that are widely used for real-time propagation of the quantum electron dynamics with an emphasis on Gaussian basis set methods. We also showcase many of the chemical applications and scientific advances made by using real-time electronic structure calculations and provide an outlook of possible new directions.

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

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2021

Structure-Activity Relationships of Hydrothermally Aged Titania-Supported Vanadium-Tungsten Oxide Catalysts for SCR of NOx Emissions with NH3

Supported V2O5-WO3/TiO2 materials employed as selective catalytic reduction (SCR) catalysts for NOx emission control from power plants. Fresh SCR catalysts usually receive exposure to harsh treatments in industry to accelerate catalyst activation (calcination in air at 650oC) and catalyst aging (hydrothermal aging at 650oC) in a way that represents various points in the catalyst/product lifetime. The present study investigates the catalyst structural and chemical changes occurring during such harsh treatments. Three series of supported V2O5-WO3/TiO2 catalysts were prepared by incipient-wetness impregnation of aqueous ammonium metavanadate and metatungstate precursors. The catalysts were subsequently dried and calcined at 550oC in O2, 650oC in O2, and hydrothermal conditions (10% O2, 8% H2O, 7% CO2 and 75% N2) at 650oC. The resulting catalysts were physically characterized by numerous techniques (in situ Raman, in situ IR, in situ High Field-High Spinning solid-state 51V MAS NMR, in situ EPR, XRD, BET surface area and ICP) and chemically probed with adsorbed ammonia, SCR-TPSR, and the SCR reaction. The surface WOx sites on the TiO2 support behave as a textural promoter that stabilizes the TiO2 (anatase) phase from sintering and transforming to the undesirable crystalline TiO2 (rutile) phase that can lead to formation of a Ti1-xVxO2 (rutile) solid solution with reduced V4+ cations (~7-15%). The surface VOx sites are mostly oligomerized as surface V+5Ox sites (~50-85% oligomers) and the extent of oligomerization tends to increase with surface WOx coverage and calcination temperature. A major difference between the calcined and hydrothermally treated catalysts was the low concentration of surface NH3* species on Lewis acid sites for the hydrothermally treated catalysts, yet the SCR activity was almost comparable for both catalysts. This finding suggests that surface NH4+*, primarily associated with the surface VOx sites, are able to efficiently perform the SCR reaction. Given that multiple catalyst parameters were simultaneously varying during these treatments, it was difficult to correlate the SCR activity with any single catalyst parameter. A correlation, however, was found between the SCR TOF/activity and the sum of the surface NH3* and NH4+* species, which is dominated by the surface NH4+* species.

SCR, NO, NH3, catalyst, V2O5, WO3, TiO2, promotion↗
DOE OSTI2021

On the Reliability of Parameter Inferences in a Multiscale Model for Transport in Stream Corridors

Nonreacting tracer tests capture information about physical processes in transient storage zones including the hyporheic zone (HZ). However, reliably extracting this information from breakthrough curves (BTCs) and distinguishing the effects of in-channel dispersion and transient storage are well-known challenges. Using BTCs from a nonreacting tracer test monitored at multiple locations, we explore ways for reliable parameter estimations. The identifiability of parameters is greatly influenced by the choice of forward and inverse modeling frameworks in addition to the quality of the data. Our forward model is a recently proposed multiscale model that uses subgrid transport models written in the Lagrangian form to represent transport along a diverse set of HZ pathways with a shape-free distribution of travel times. Joint distributions of HZ and channel parameters are estimated using the Markov Chain Monte Carlo technique. Numerical experiments show ambiguity between channel dispersion and HZ transport when the reach length is too short to allow significant solute-HZ interaction, the observation period is too brief to observe the tailing behavior, or the solute source is spread in time. In contrast, we obtained reliable parameter estimates by simultaneously fitting BTCs observed at different locations in the test reach using a single set of HZ parameters and section-specific channel areas and dispersion coefficients. Furthermore, this study demonstrates the estimation of travel time distributions, HZ exchange rates, and channel parameters in a new multiscale approach and offers guidance for extracting reliable parameter estimates from multiple BTCs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Phase diagram of a distorted kagome antiferromagnet and application to Y-kapellasite

Abstract We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings J ⬡ , J , and $$J^{\prime}$$ J ′ , and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear $$\overrightarrow{Q}=0$$ Q → = 0 magnetic phase, two unusual non-collinear coplanar $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y 3 Cu 9 (OH) 19 Cl 8 is a realization of this model and predict its ground state to lie in the region of $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) order, which remains stable even after the inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spectral localizer for line-gapped non-Hermitian systems

Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.

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

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamics of E6 chiral gauge theories

We present exact nonperturbative vacuum solutions to chiral gauge theories based on the gauge group and several matter fermions in the fundamental -dimensional representation. They are obtained when supersymmetric versions are perturbed by small supersymmetry breaking by anomaly mediation. The universality classes obtained are very different from what can be conjectured by the tumbling hypothesis. In particular, the case with three may have an unbroken SU(3) symmetry with massless composite fermions in of SU(3). For this case, we employed numerical techniques to obtain the exact ground state.

Goh, Andrew↗
DOE OSTI2021

Various facets of magnetic charge correlation: Micromagnetic and distorted-wave Born approximation simulations study

The emergent concept of the magnetic charge quasiparticle provides a new realm to study the evolution of magnetic properties in two-dimensional artificially frustrated magnets. Here we report on the exploration of magnetic phases due to various magnetic charge correlation using the complementary numerical techniques of micromagnetic and distorted-wave Born approximation simulations in artificial permalloy honeycomb lattice. The honeycomb element length varies between 10 nm and 100 nm, while the width and thickness are kept within the single domain limit. In addition to the charge ordered loop state, we observe disordered charge arrangement, characterized by the random distribution of ±Q charges, in single domain size honeycomb lattice. As the length of the honeycomb element increases, low multiplicity magnetic charges tend to form contiguous bands in thinner lattice. Thin honeycomb lattice with 100 nm element length exhibits a perfect spin ice pattern, which remains unaffected by the modest increase in the width of element size. We simulate scattering profiles under the pretext of distorted-wave Born approximation formalism for the micromagnetic phases. The results are expected to provide useful guidance in the experimental investigation of magnetic phases in an artificial honeycomb magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Frustrated mixed-spin ladders: Evidence for a bond-order wave phase between rung-singlet and Haldane phases

In frustrated spin ladders, the interplay of frustration and correlations leads to the familiar Haldane (H) and rung-singlet (RS) phases. The nature of the transition between these two phases is still under debate. In this paper we tackle this issue using tools of quantum information theory. We consider frustrated mixed-spin-(1, 1/2) ladders with antiferromagnetic leg, rung and diagonal couplings, and calculate various quantities, such as the entanglement entropy (EE), the Schmidt gap, and the level degeneracy of the entanglement spectrum (ES). We use two numerical techniques, the infinite time-evolving block decimation (iTEBD) and the density matrix renormalization group (DMRG). We demonstrate that there exists an intermediate phase in which the ES levels do not exhibit the characteristic degeneracies of the H and RS phases. To understand the underlying physics in this phase, here we investigate short-range spin correlations along legs, rungs and diagonals and show that in this intermediate phase long-wavelength modulations occur, akin to bond-order waves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Heavy-ion fusion cross section formula and barrier height distribution

Methods for obtaining fusion cross section formulas are discussed, especially for those that take on an empirical form. A new expression starting with σ( E) = $\frac{1}{E}$ $\int ^{E}_{E_0}$ $(\int^{E'}_{E_0} B (E^\shortparallel) dE^\shortparallel)dE'$ has been explored. Here, B(E) is a reasonable, assumed function of the barrier height distribution, $\frac{d^2 (\sigma E)}{dE^2}$. Further, the resulting analytic cross section formula reproduces very well the excitation functions for many light and heavy fusion systems across wide energy ranges, when B(E) is assumed to be a multi-Gaussian function. This study offers an improved determination of the fusion barrier height distribution over other numerical techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

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↗