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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1984

Two time scale design of output feedback systems

Singular perturbation techniques are used in this paper to develop a two time scale procedure for designing static gain output feedback controllers. It is shown that, if certain control spillover conditions are satisfied, control designs based on reduced-order models will stabilize the corrgsponding full system. Optimal output feedback control theory is used to derive the necessary conditions for the stabilizing gain matrix. The problem of stabilizing a model for a large space structure is used to illustrate the practicality of the approach.

Calise, A. J.↗
DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Influence of Incidence Angle on Sound Generation by Airfoils Interacting with High-Frequency Gusts

A theoretical model is developed for the sound generated when a convected vortical or entropic gust encounters an airfoil at nonzero angle of attack. The theory is based on a linearization of the Euler equations about the steady, subsonic flow past the airfoil. High-frequency gusts, whose wavelengths are short compared to the airfoil chord, but long compared to the displacement of the mean-flow stagnation point from the leading edge, are considered. The analysis utilizes singular perturbation techniques and involves four asymptotic regions. Local regions, which scale on the gust wavelength, are present at the airfoil leading and trailing edges. Behind the airfoil a "transition" region, which is similar to the transition zone between illuminated and shadow zones in optical problems, is present. In the outer region, far away from the airfoil edges and wake, the solution has a geometric-acoustics form. The primary sound generation is found to be concentrated in the local leading-edge region. The trailing edge plays a secondary role as a scatterer of the sound generated in the leading-edge region. Parametric calculations are presented which illustrate that moderate levels of airfoil steady loading can significantly affect the sound field produced by airfoil-gust interactions.

Matthew R. Myers↗
DOE OSTI2023

Machine Learning-Based Model Predictive Control of Two-Time-Scale Systems

In this study, we present a general form of nonlinear two-time-scale systems, where singular perturbation analysis is used to separate the dynamics of the slow and fast subsystems. Machine learning techniques are utilized to approximate the dynamics of both subsystems. Specifically, a recurrent neural network (RNN) and a feedforward neural network (FNN) are used to predict the slow and fast state vectors, respectively. Moreover, we investigate the generalization error bounds for these machine learning models approximating the dynamics of two-time-scale systems. Next, under the assumption that the fast states are asymptotically stable, our focus shifts toward designing a Lyapunov-based model predictive control (LMPC) scheme that exclusively employs the RNN to predict the dynamics of the slow states. Additionally, we derive sufficient conditions to guarantee the closed-loop stability of the system under the sample-and-hold implementation of the controller. A nonlinear chemical process example is used to demonstrate the theory. In particular, two RNN models are constructed: one to model the full two-time-scale system and the other to predict solely the slow state vector. Both models are integrated within the LMPC scheme, and we compare their closed-loop performance while assessing the computational time required to execute the LMPC optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Topological orders of monopoles and hedgehogs: From electronic and magnetic spin-orbit coupling to quarks

Topological states of matter are, generally, quantum liquids of conserved topological defects. We establish this by constructing and analyzing topological field theories which introduce gauge fields to describe the dynamics of singularities in the original field configurations. In this work, homotopy groups are utilized to identify topologically protected singularities, and the conservation of their protected number is captured by a topological action term that unambiguously obtains from the given set of symmetries. Stable phases of these theories include quantum liquids with emergent massless Abelian and non-Abelian gauge fields, as well as topological orders with long-range quantum entanglement, fractional excitations, boundary modes, and unconventional responses to external perturbations. This paper focuses on the derivation of topological field theories and basic phenomenological characterization of topological orders associated with homotopy groups π n (S n ), n ≥ 1 . These homotopies govern monopole and hedgehog topological defects in d = n + 1 dimensions, and enable the generalization of both weakly interacting and fractional quantum Hall liquids of vortices to d > 2 . Hedgehogs have not been in the spotlight so far, but they are particularly important defects of magnetic moments because they can be stimulated in realistic systems with spin-orbit coupling, such as chiral magnets and d = 3 topological materials. We predict topological orders in systems with U(1) × Spin(d) symmetry in which fractional electric charge attaches to hedgehogs. Monopoles, the analogous defects of charge or generic U(1) currents, may bind to hedgehogs via Zeeman effect, or effectively emerge in purely magnetic systems. The latter can lead to spin liquids with different topological orders than that of the resonant valence bond spin liquid. Charge fractionalization of quarks in atomic nuclei is also seen as possibly arising from the charge-hedgehog attachment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Resurgence of the renormalization group equation

We show how the renormalons emerge from the renormalization group equation with a priori no reference to any Feynman diagrams. The proof is rather given by recasting the renormalization group equation as a resurgent equation studied in the mathematical literature, which describes a function with an infinite number of singularities in the positive axis of the Borel plane. Consistency requires a one-to-one correspondence between the existence of such kind of equation and the actual (generalized) Borel resummation of the renormalons through a one-parameter transseries. Our finding suggests how non-perturbative contributions can affect the running couplings. We also discuss these concepts within the context of gauge theories, making use of the large number of flavor expansion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Eigenvalues of singular differential operators by finite difference methods. II.

Note is made of an earlier paper which defined finite difference operators for the Hilbert space L2(m), and gave the eigenvalues for these operators. The present work examines eigenvalues for higher order singular differential operators by using finite difference methods. The two self-adjoint operators investigated are defined by a particular value in the same Hilbert space, L2(m), and are strictly positive with compact inverses. A class of finite difference operators is considered, with the idea of application to the theory of Toeplitz matrices. The approximating operators consist of a good approximation plus a perturbing operator.

Baxley, J. V.↗
NASA NTRS1974

Theoretical prediction of thick wing and pylon-fuselage-fanpod-nacelle aerodynamic characteristics at subcritical speeds. Part 1: Theory and results

The theoretical development and the comparison of results with data of a thick wing and pylon-fuselage-fanpod-nacelle analysis are presented. The analysis utilizes potential flow theory to compute the surface velocities and pressures, section lift and center of pressure, and the total configuration lift, moment, and vortex drag. The skin friction drag is also estimated in the analysis. The perturbation velocities induced by the wing and pylon, fuselage and fanpod, and nacelle are represented by source and vortex lattices, quadrilateral vortices, and source frustums, respectively. The strengths of these singularities are solved for simultaneously including all interference effects. The wing and pylon planforms, twists, cambers, and thickness distributions, and the fuselage and fanpod geometries can be arbitrary in shape, provided the surface gradients are smooth. The flow through nacelle is assumed to be axisymmetric. An axisymmetric center engine hub can also be included. The pylon and nacelle can be attached to the wing, fuselage, or fanpod.

Tulinius, J. R.↗
NASA NTRS2009

Stability of Detached Solidification

Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.

Mazuruk, K.↗
DOE OSTI2023

Integrable field theories and their CCFT duals

We compute the Mellin transforms of various two-dimensional integrable S-matrices, providing the first explicit, non-perturbative realizations of celestial CFT. In two dimensions, the Mellin transform is simply the Fourier transform in rapidity space, and the “celestial correlator” has no position dependence. The simplified setting allows us to study the analytic properties of CCFT correlators exactly as a function of the conformal dimensions. We find that the correlators exist as real distributions of the conformal weights, with asymptotics controlled by the mass spectrum and three-point couplings of the model. Coupling these models to a flat space limit of JT gravity preserves integrability and dresses the amplitudes by a rapidly varying gravitational phase. We find that the coupling to gravity smooths out certain singular aspects of the Mellin-transformed correlators.

2D Gravity↗
NASA NTRS1990

Theory of finite disturbances in a centrifugal compression system with a vaneless radial diffuser

A previous small perturbation analysis of circumferential waves in circumferential compression systems, assuming inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, when boundary layer displacement effect is included. The Moore-Greitzer analysis for finite strength transients containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under certain assumptions, and except for a new second order swirl, the diffuser velocity field, including resonant singularities, can be carried over from the previous inviscid linear analysis. Nonlinear transient equations are derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded by an outer surge response on a longer time scale. Surge may occur purely as result of circumferential wave resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the characteristic slope. A number of circumferential modes considered simultaneously is briefly discussed.

Moore, F. K.↗
DOE OSTI2023

Quasi-Fermi liquid behavior in a one-dimensional system of interacting spinless fermions

We present numerical evidence for a paradigm in one-dimensional interacting fermion systems, whose phenomenology has traits of both Luttinger liquids and Fermi liquids. This state, dubbed a quasi-Fermi liquid, possesses a discontinuity in its fermion occupation number at the Fermi momentum. The excitation spectrum presents particlelike quasiparticles and absence of holelike quasiparticles, giving rise instead to edge singularities. Such a state is realized in a one-dimensional spinless fermion lattice Hamiltonian by fine-tuning the interactions to a regime where they become irrelevant in the renormalization group sense. We show, using uniform infinite matrix products states and finite-entanglement scaling analysis, that the system ground state is characterized by a Luttinger parameter K = 1 and a discontinuous jump in the fermion occupation number. We support the characterization with calculations of the spectral function that show a particle-hole asymmetry reflected in the existence of well-defined Landau quasiparticles above the Fermi level and edge singularities without the associated quasiparticles below. Furthermore, these results indicate that the quasi-Fermi liquid paradigm can be realized beyond the low-energy perturbative realm.

1-dimensional systems↗
NASA NTRS1986

Influence of airfoil camber on convected gust interaction noise

This paper investigates the effect of airfoil steady loading on the sound generated by the interaction of an airfoil with a convected disturbance. A previous theory, which included only the incidence angle contribution to the mean loading, is extended to include camber. The theory is based on a linearization of the Euler equations about a nonuniform, 0(1) Mach number subsonic mean flow. The discussion concentrates on the case of a slightly cambered airfoil at small incidence angle, interacting with a gust whose wavelength is short compared to the airfoil chord. The small parameter representing the amount of camber and incidence, and the large parameter representing the ratio of airfoil chord to disturbance wavelength, are utilized in a singular perturbation solution to the governing equations. Acoustic power calculations reveal that the amount of sound generated increases significantly with increased loading. More importantly, it is shown that the radiated acoustic power correlates very well with the strength of the mean flow around the leading edge.

Myers, M. R.↗
DOE OSTI2024

Renormalization of asymmetric staple-shaped Wilson-line operators in lattice and continuum perturbation theory

In this work, we study the renormalization of nonlocal quark bilinear operators containing an asymmetric staple-shaped Wilson line at the one-loop level in both lattice and continuum perturbation theory. These operators enter the first-principle calculation of transverse momentum-dependent parton distribution functions (TMDPDFs) in lattice QCD using the formulation of large momentum effective theory. We provide appropriate RI ′ -type conditions that address the power and logarithmic divergences, as well as the mixing among staple operators of different Dirac structures, using a number of different possible projectors. A variant of RI ′ , including calculations of rectangular Wilson loops, which cancel the pinch-pole singularities of the staple operators at infinite length and reduce residual power divergences, is also employed. We calculate at one-loop order the conversion matrix, which relates the quasi-TMDPDFs in the RI ′ -type schemes to the reference scheme MS ¯ for arbitrary values of the renormalization momentum scale and of the dimensions of the staple. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Ladder symmetries and Love numbers of Reissner-Nordström black holes

It is well known that asymptotically flat black holes in general relativity have vanishing tidal Love numbers. In the case of Schwarzschild and Kerr black holes, this property has been shown to be a consequence of a hidden structure of ladder symmetries for the perturbations. In this work, we extend the ladder symmetries to non-rotating charged black holes in general relativity. As opposed to previous works in this context, we adopt a more general definition of Love numbers, including quadratic operators that mix gravitational and electromagnetic perturbations in the point-particle effective field theory. We show that the calculation of a subset of those couplings in full general relativity is affected by an ambiguity in the split between source and response, which we resolve through an analytic continuation. As a result, we derive a novel master equation that unifies scalar, electromagnetic and gravitational perturbations around Reissner-Nordström black holes. The equation is hypergeometric and can be obtained from previous formulations via nontrivial field redefinitions, which allow to systematically remove some of the singularities and make the presence of the ladder symmetries more manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Perturbative unitarity and NEC violation in genesis cosmology

Explorations of the violation of null energy condition (NEC) in cosmology could enrich our understanding of the very early universe and the related gravity theories. Although a fully stable NEC violation can be realized in the “beyond Horndeski” theory, it remains an open question whether a violation of the NEC is allowed by some fundamental properties of UV-complete theories or the consistency requirements of effective field theory (EFT). We investigate the tree-level perturbative unitarity for stable NEC violations in the contexts of both Galileon and “beyond Horndeski” genesis cosmology, in which the universe is asymptotically Minkowskian in the past. We find that the constraints of perturbative unitarity imply that we may need some unknown new physics below the cut-off scale of the EFT other than that represented by the “beyond Horndeski” operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗