Search NASA⌕ Search

SEARCH · Search NASA

Results for “singular perturbation theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Control Design of Series-connected PV-powered Grid-forming Converters via Singular Perturbation

We report modular architectures that consist of several series-connected dc-ac converters have been a focal point of recent innovations in transformerless medium-voltage applications. In this paper, we consider an architecture consisting of dc-ac modules containing a quadruple active bridge dc-dc converter, which generates three floating dc links that feed grid-side dc-ac inverters. Practical implementation of such a converter module in PV systems requires a variety of controllers that collectively achieve maximum power point tracking, dc-link regulation, and ac-side power control. Design of such multi-loop systems is generally quite challenging due to the potential for destabilizing interactions among loops. Here we propose a design approach where singular perturbation theory is used to decompose the timescales at which each control loop operates and provides a systematic framework for parametric selection. Our approach also ensures system stability of multiple modules with identical controls connected in series across a grid. The paper concludes with experimental results of three 1000 W series-connected converter modules across a stiff grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Singular Perturbation-Based Large-Signal Order Reduction of Microgrids for Stability and Accuracy Synthesis With Control

The increasing penetration of distributed energy resources (DERs) highlights the growing importance of microgrids (MGs) in enhancing power system reliability. Employing electromagnetic transient (EMT) analysis in MGs becomes crucial for controlling the rapid transients. However, this requires an accurate but high-order model of power electronics and their underlying control loops, complexifying the stability analysis from the viewpoint of a higher control level. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging our derived reduced-order model (ROM) and boundary layer model (BLM). Furthermore, this theorem furnishes a set of conditions that determine the accuracy of the developed ROM. Lastly, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. Case studies are conducted on MG systems to show the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Interacting Dirac magnons in the van der Waals ferromagnet CrBr 3

We study the effects of magnon-magnon interactions in the two-dimensional van der Waals ferromagnet CrBr 3 focusing on its honeycomb lattice structure. Motivated by earlier theoretical predictions of temperature-induced spectral shifts and van Hove singularities in the magnon dispersion [S. S. Pershoguba et al., Phys. Rev. X 8, 011010 (2018)], we go beyond the commonly used thermal magnon approximation by applying second-order perturbation theory in a fully numerical framework. Our analysis uncovers significant deviations from previous analysis: in particular, the predicted singularities are absent, consistent with recent inelastic neutron scattering measurements [S. E. Nikitin et al ., Phys. Rev. Lett. 129, 127201 (2022)]. Moreover, we find that the temperature dependence of the renormalized magnon spectrum exhibits a distinct 𝑇 3 behavior for the optical magnon branch, while retaining 𝑇 2 behavior for the acoustic or down magnon band. This feature sheds light on the collective dynamics of Dirac magnons and their interactions. We further compare the honeycomb case with a triangular Bravais lattice, relevant for ferromagnetic monolayer MnBi 2 ⁢Te 4 , and show that both systems lack singular features while displaying quite distinct thermal trends.

Holstein-Primakoff method↗
DOE OSTI2020

Classical double copy of nonsingular black holes

We apply the classical double copy procedure to a class of regular, non-singular black hole solutions. We give several examples, paying particular attention to a string-theory-corrected black hole solution emerging from T-duality. Non-perturbative stringy corrections introduce an ultraviolet (UV) zero-point length cutoff which results in non-singular black hole spacetimes. Apart from the UV regulator, the solution is equivalent to the Bardeen black hole spacetime. We extend this solution to include an asymptotic de Sitter background. All Yang-Mills field theory quantities associated with the double copy are well-behaved and finite for all values of parameters. Furthermore, we present a thorough analysis of the black hole horizon structure, additionally uncovering a simple yet new connection between horizons on the gravity side and electric fields on the gauge theory side of the double copy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Nonperturbative resolution of strong coupling singularities in 4D 𝒩 = 1 heterotic M-theory

We investigate the interior of the moduli space of four-dimensional 𝒩 = 1 theories of gravity arising from compactifications of the 𝐸 8 × 𝐸 8 heterotic string on Calabi-Yau threefolds. By studying the threshold corrections to the coupling of the heterotic gauge groups, we infer the existence of a strong coupling singularity for one of the perturbative heterotic gauge groups, which effectively yields an additional finite distance boundary of the classical scalar field space. In heterotic M-theory, this boundary maps to a domain wall solution for which the gauge coupling and the warp factor on one of the Horava-Witten 9-branes diverge, thus highlighting the gravitational origin of the classical strong coupling singularity. The divergence of the warp factor is, however, regulated once nonperturbative effects are taken into account, as we demonstrate by studying the instanton corrections to the 5D BPS domain wall equations. This regularization implies that the classical strong coupling boundary of the scalar field 4D 𝒩 = 1 heterotic M-theory is resolved, indicating that, at the quantum level, the field space can be extended beyond this classical boundary.

Cvetič, Mirjam↗
DOE OSTI2021

Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: $\phi$ 3 QFT in 6 Dimensions

We analyze the asymptotically free massless scalar $\phi$ 3 quantum field theory in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer–Connes Hopf-algebraic Dyson–Schwinger equations for the anomalous dimension. This scalar conformal field theory is asymptotically free and has a real Lipatov instanton. In the Hopf-algebraic approach we find a trans-series having an intricate Borel singularity structure, with three distinct but resonant non-perturbative terms, each repeated in an infinite series. These expansions are in terms of the renormalized coupling. The resonant structure leads to powers of logarithmic terms at higher levels of the trans-series, analogous to logarithmic terms arising from interactions between instantons and anti-instantons, but arising from a purely perturbative formalism rather than from a semi-classical analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Finite-size effects in periodic coupled cluster calculations

Here, we provide the first rigorous study of the finite-size error in the simplest and representative coupled cluster theory, namely the coupled cluster doubles (CCD) theory, for gapped periodic systems. Given exact Hartree-Fock orbitals and their corresponding orbital energies, we demonstrate that the correlation energy obtained from the approximate CCD method, after a finite number of fixed-point iterations over the amplitude equation, exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Here $N_k$ is the number of discretization points in the Brillouin zone and characterizes the system size. Under additional assumptions ensuring the convergence of the fixed-point iterations, we demonstrate that the CCD correlation energy also exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Our analysis shows that the dominant error lies in the coupled cluster amplitude calculation, and the convergence of the finite-size error in energy calculation can be boosted to $\mathcal{O}(N^{-1}_k)$ with accurate amplitudes. This also provides the first proof of the scaling of the finite-size error in the third order Møller-Plesset perturbation theory (MP3) for periodic systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Cosmological perturbations from a new approach to inflation

In a previous paper, we proposed a new approach to the beginning of inflation—a lingering universe [Melcher et al. Phys. Rev. D 110, 063517 (2024).]—that is a quasistatic phase prior to inflation. The Universe begins in a lingering state with a nearly vanishing Hubble parameter. This calls into question the absolute age of the Universe, as the Hubble time can be nearly infinite. It also provides promise for addressing the initial singularity of inflation and issues with quantum field theory in de Sitter spacetime. Such models arise in classical cosmologies with nonvanishing spatial curvature (inspired by PLANCK 2018 data) and independently by models that arise in string cosmology. Here, in this paper, we consider the importance of cosmological perturbations for the stability of the lingering phase and how this influences cosmological observations. Our goal is to establish observables in this new paradigm for the origin of inflation, which is in contrast to eternal inflation and cyclic cosmologies. We also address questions of stability and the transition to inflation.

cosmology↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2023

Gluon transverse-momentum-dependent distributions from large-momentum effective theory

We demonstrate that gluon transverse-momentum-dependent parton distribution functions (TMDPDFs) can be extracted from lattice calculations of appropriate Euclidean correlations in large-momentum effective theory (LaMET). Based on perturbative calculations of gluon unpolarized and helicity TMDPDFs, we present a matching formula connecting them and their LaMET counterparts, where the latter are renormalized in a scheme facilitating lattice calculations and converted to the $\overline{MS}$ scheme. The hard matching kernel is given up to one-loop level. We also show that the perturbative result is independent of the prescription used for the pinch-pole singularity in the relevant correlations. Our results offer a guidance for the extraction of gluon TMDPDFs from lattice simulations, and have the potential to greatly facilitate perturbative calculations of the hard matching kernel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗