Search NASA⌕ Search

SEARCH · Search NASA

Results for “ASYMPTOTIC SOLUTIONS”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Trivializing Flow in 2D-O(3) model

The two dimensional O(3) sigma model, just as quantum chromodynamics, is an asymptotically free theory with a mass gap. Therefore, it is an interesting and simple toy model to investigate algorithms for Markov Chain Monte Carlo simulations of quantum chromodynamics. In this talk, we discuss the construction of a trivializing map, a field transformation from a given theory to a trivial one, through a suitably chosen gradient flow. An analytic solution for the generating functional of this trivializing flow has been obtained by a perturbative expansion in the flow time. Utilizing this solution allows for new approaches to be considered when proposing updates for a Markov Chain Monte Carlo algorithm.

Chamness, Christopher↗
DOE OSTI2022

Fast nonlinear iterative solver for an implicit, energy-conserving, asymptotic-preserving charged-particle orbit integrator

Here recently, an asymptotic-preserving (AP) particle orbit integrator has been proposed with remarkable properties including exact energy conservation, the ability to capture of all first-order drifts (including the ∇B-drift), the ability to capture trapped-passing boundaries with parallel velocity extremely close to the critical velocity, and the ability to transition from strongly to weakly magnetized spatial regions. The new AP orbit integrator is implicit, employing a Crank-Nicolson (CN) temporal discretization to ensure exact energy conservation. This, in turn, requires a local nonlinear iteration involving particle velocities and positions, and the local electromagnetic fields, to obtain the new-time solution. Ref. [1] did not attempt to provide an efficient solver for this system, and employed a brute-force GMRES-driven Jacobian-free Newton-Krylov (JFNK) solver to invert the particle orbit equations at every timestep for expediency. While JFNK is robust and reliable, it is also expensive and very intrusive for practical implementations of the method (it requires having the JFNK machinery available and solving a 6 x 6 Jacobian system iteratively once per iteration per particle).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Stochastic Transport Model for the Cumulative Number of Fissions and Deposited Fission Energy

The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fueling limits in a cylindrical viscosity-limited reactor

Recently, a method to achieve a “natural hot-ion mode” was suggested by utilizing ion viscous heating in a rotating plasma with a fixed boundary. We explore the steady-state solution to the Braginskii equations and find the parameter regime in which a significant temperature difference between ions and electrons can be sustained in a driven steady state. Here, the threshold for this effect occurs at ρ i ≳ 0.1R. An analytic, leading order low flow solution is obtained, and a numerical, moderate Mach number M ≲ 2 is investigated. The limitation is found to be at moderate Mach numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Medium-induced radiative kernel with the Improved Opacity Expansion

We calculate the fully differential medium-induced radiative spectrum at next-to-leading order (NLO) accuracy within the Improved Opacity Expansion (IOE) framework. This scheme allows us to gain analytical control of the radiative spectrum at low and high gluon frequencies simultaneously. The high frequency regime can be obtained in the standard opacity expansion framework in which the resulting power series diverges at the characteristic frequency ω c ~ q^L 2 . In the IOE, all orders in opacity are resumed systematically below ωc yielding an asymptotic series controlled by logarithmically suppressed remainders down to the thermal scale T « ω c , while matching the opacity expansion at high frequency. Furthermore, we demonstrate that the IOE at NLO accuracy reproduces the characteristic Coulomb tail of the single hard scattering contribution as well as the Gaussian distribution resulting from multiple soft momentum exchanges. Finally, we compare our analytic scheme with a recent numerical solution, that includes a full resummation of multiple scatterings, for LHC-inspired medium parameters. Furthermore, we find a very good agreement both at low and high frequencies showcasing the performance of the IOE which provides for the first time accurate analytic formulas for radiative energy loss in the relevant perturbative kinematic regimes for dense media.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

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

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2020

Beurling-Selberg extremization and modular bootstrap at high energies

We consider previously derived upper and lower bounds on the number of operators in a window of scaling dimensions [\Delta - \delta,\Delta + \delta] [ Δ − δ , Δ + δ ] at asymptotically large \Delta Δ in 2d unitary modular invariant CFTs. These bounds depend on a choice of functions that majorize and minorize the characteristic function of the interval [\Delta - \delta,\Delta + \delta] [ Δ − δ , Δ + δ ] and have Fourier transforms of finite support. The optimization of the bounds over this choice turns out to be exactly the Beurling-Selberg extremization problem, widely known in analytic number theory. We review solutions of this problem and present the corresponding bounds on the number of operators for any \delta \geq 0 δ ≥ 0 . When 2\delta \in \mathbb Z_{\geq 0} 2 δ ∈ ℤ ≥ 0 the bounds are saturated by known partition functions with integer-spaced spectra. Similar results apply to operators of fixed spin and Virasoro primaries in c>1 c > 1 theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Novel features of attractors and transseries in nonconformal Bjorken flows

In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the nonconformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the nonconformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. Further, these terms propagate into the perturbative and nonperturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient nonhydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of local attraction at early and late times for the shear viscous tensor and bulk viscous pressure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Exact WKB and the quantum Seiberg-Witten curve for 4d $N$ = 2 pure SU(3) Yang-Mills. Abelianization

We investigate the exact WKB method for the quantum Seiberg-Witten curve of 4d N = 2 pure SU(3) Yang-Mills in the language of abelianization. The relevant differential equation is a third-order equation on $\mathbb{CP}^1$ with two irregular singularities. We employ the exact WKB method to study the solutions to such a third-order equation and the associated Stokes phenomena. We also investigate the exact quantization condition for a certain spectral problem. Moreover, exact WKB analysis leads us to consider new Darboux coordinates on a moduli space of flat SL(3,$\mathbb{C}$)-connections. In particular, in the weak coupling region we encounter coordinates of the higher length-twist type generalizing Fenchel-Nielsen coordinates. The Darboux coordinates are conjectured to admit asymptotic expansions given by the formal quantum periods series and we perform numerical analysis supporting this conjecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ladder symmetries of black holes. Implications for love numbers and no-hair theorems

It is well known that asymptotically flat black holes in general relativity have a vanishing static, conservative tidal response. Here, we show that this is a result of linearly realized symmetries governing static (spin 0,1,2) perturbations around black holes. The symmetries have a geometric origin: in the scalar case, they arise from the (E)AdS isometries of a dimensionally reduced black hole spacetime. Underlying the symmetries is a ladder structure which can be used to construct the full tower of solutions, and derive their general properties: (1) solutions that decay with radius spontaneously break the symmetries, and must diverge at the horizon; (2) solutions regular at the horizon respect the symmetries, and take the form of a finite polynomial that grows with radius. Taken together, these two properties imply that static response coefficients — and in particular Love numbers — vanish. Moreover, property (1) is consistent with the absence of black holes with linear (perturbative) hair. We also discuss the manifestation of these symmetries in the effective point particle description of a black hole, showing explicitly that for scalar probes the worldline couplings associated with a non-trivial tidal response and scalar hair must vanish in order for the symmetries to be preserved.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2025

Asymptotic Relaxation of Moment Equations for a Multi-species, Homogeneous BGK Model

Multi-species BGK models describe the dynamics of rarefied gases with constituent particles of different elements or compounds with potentially nontrivial velocity distributions. Here, in this paper, moment equations for the bulk velocities, energies, and temperatures of a spatially homogeneous multi-species BGK model are examined. A key challenge in analyzing these equations is the fact that the collision frequencies are allowed to depend on the species temperatures, which allows for more realistic simulations of dilute gas flow. Therefore, a positive lower bound is established for the species temperatures. With this lower bound, a global existence and uniqueness of solutions to the coupled velocity-energy ODE system is established. The lower bound also enables a proof of exponential decay to a unique steady-state solution. Numerical results are presented to demonstrate how the bulk velocities and temperatures relax for large times.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Holographic duality for Ising CFT with boundary

We extend the holographic duality between 3d pure gravity and 2d Ising CFT proposed in ref. [1] to CFTs with boundaries. Besides the usual asymptotic boundary, the dual bulk spacetime now has a real cutoff, on which live branes with finite tension, giving Neumann boundary condition on the metric tensor. The strongly coupled bulk theory requires that we dress the well-known semiclassical AdS/BCFT answer with boundary gravitons, turning the partition function into the form of Virasoro characters. Using this duality, we relate the brane tensions to the modular S-matrix elements of the dual BCFT and derive the transformation between gravitational solutions with different brane tensions under modular S action.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Quantum Hamilton-Jacobi theory, spectral path integrals, and exact WKB analysis

We propose a new way to perform path integrals in quantum mechanics by using a quantum version of Hamilton-Jacobi (HJ) theory. In classical mechanics, Hamilton-Jacobi theory is a powerful formalism, however, its utility is not explored in quantum theory beyond approximation schemes. The canonical transformation enables one to set the new Hamiltonian to constant or zero, but keeps the information about solution in Hamilton’s characteristic function. To benefit from this in quantum theory, one must work with a formulation in which classical Hamiltonian is used. This uniquely points to phase space path integral. However, the main variable in HJ formalism is energy, not time. Thus, we are led to consider the Fourier transform of the path integral, the spectral path integral Z ˜ ( E ) . The evaluation of path integrals reduces to determining the quantum Hamilton characteristic functions (which can be achieved via an asymptotic analysis) and a discrete sum over the quantum period lattice, generalizing Gutzwiller’s sum. Published by the American Physical Society 2025

Türe, Mustafa (ORCID:0009000975968618)↗
DOE OSTI2021

Modified Schwarzschild metric from a unitary accelerating mirror analog

We present a modified Schwarzschild solution for a model of evaporation of a black hole with information preservation. By drawing a direct analogy to the quantum pure accelerating mirror (dynamical Casimir effect of a 1D horizon), we derive a Schwarzschild metric with not only the usual Schwarzschild radius but an additional length scale related to the Planck length. The black hole has thermal particle production that leads to complete evaporation of the black hole, resulting in non-divergent entanglement entropy, Page curve turn-over, and an asymptotic quantum pure state with no information loss.

79 ASTRONOMY AND ASTROPHYSICS↗