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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Efficient and accurate computation of the incomplete Airy functions

The incomplete Airy integrals serve as canonical functions for the uniform ray optical solutions to several high-frequency scattering and diffraction problems that involve a class of integrals characterized by two stationary points that are arbitrarily close to one another or to an integration endpoint. Integrals with such analytical properties describe transition region phenomena associated with composite shadow boundaries. An efficient and accurate method for computing the incomplete Airy functions would make the solutions to such problems useful for engineering purposes. In this paper a convergent series solution for the incomplete Airy functions is derived. Asymptotic expansions involving several terms are also developed and serve as large argument approximations. The combination of the series solution with the asymptotic formulae provides for an efficient and accurate computation of the incomplete Airy functions. Validation of accuracy is accomplished using direct numerical integration data.

Constantinides, E. D.↗

NASA NTRS1985

Propagation and stability of wavelike solutions of finite difference equations with variable coefficients

The propagation and dissipation of wavelike solutions to finite difference equations is analyzed on the basis of an asymptotic approach in which a wave solution is expressed as a product of a complex amplitude and an oscillatory phase function whose frequency and wavenumber may also be complex. An asymptotic expansion leads to a local dispersion relation for wavenumber and frequency; the first-order terms produce an equation for the amplitude in which the local group velocity appears as the convection velocity of the amplitude. Equations for the motion of wavepackets and their interaction at boundaries are derived, and a global stability analysis is carried out.

Giles, M. B.↗

NASA NTRS1982

Optimal two-impulse transfer between specified terminal states of Keplerian orbits

The integrals of motion for an inverse-square force field and a necessary condition for optimality are used in a transfer orbit co-ordinate system to formulate the solution of the optimal two-impulse transfer between fixed position and velocity vectors on Keplerian orbits. In this co-ordinate system, the equations reveal two asymptotes that are useful in analysing the solution. It is shown that there are only two real extremals (both minimums) which are separated by an asymptote. An approximate analytic solution is also obtained by ignoring a quadratic term whose coefficient is approximately zero for a large class of orbit transfer problems.

Jezewski, D. J.↗

DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Space–time structure of weak magnetohydrodynamic turbulence

The two-time energy spectrum of weak magnetohydrodynamic turbulence is found by applying a wave-turbulence closure to the cumulant hierarchy constructed from the dynamical equations. Solutions are facilitated via asymptotic expansions in terms of the small parameter $\varepsilon$ , describing the ratio of time scales corresponding to Alfvénic propagation and nonlinear interactions between counter-propagating Alfvén waves. The strength of nonlinearity at a given spatial scale is further quantified by an integration over all possible delta-correlated modes compliant in a given set of three-wave interactions that are associated with energy flux through the said scale. The wave-turbulence closure for the two-time spectrum uncovers a secularity occurring on a time scale of order $\varepsilon ^{-2}$ , and the asymptotic expansion for the spectrum is reordered in a manner comparable to the one-time case. It is shown that for the regime of stationary turbulence, the two-time energy spectrum exponentially decays on a lagged time scale $(\varepsilon ^2 \gamma _k^s)^{-1}$ in proportion to the strength of the associated three-wave interactions, characterized by nonlinear decorrelation frequency $\gamma _k^s$ . The scaling of the form $k_{\perp } v_0 \chi _0$ exhibited by this frequency is reminiscent of random sweeping by the outer scale with characteristic fluctuation velocity $v_0$ that is modified due to competition with Alfvénic propagation (characterized by $\chi _0$ ) at the said scale. A brief calculation of frequency broadening of the power spectrum due to nonlinear interactions is also presented.

Physics↗

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 OSTI2021

Energy-entropy relation for asymptotically Lifshitz spacetimes with universal horizons

We numerically solve for 2 + 1 asymptotically Lifshitz universal horizon solutions in Horava-Lifshitz gravity for dynamical exponents z = 2 through z = 16 . We find that for all z there is a thermodynamical first law. Furthermore, we find that the energy-entropy relation or Smarr formula expected for a thermal state in a two-dimensional Lifshitz field theory, E = 2 z + 2 T S , is also satisfied for universal horizons, including the correct z scaling.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Separability of the Dirac equation in the STU black hole spacetime: Pairwise-equal charge case study

We study the separability of the Dirac equation in the background of four dimensional charged rotating asymptotically flat black hole solution of supergravity and the low energy string theory known as STU black hole. In particular, we analyze in detail the separability conditions in the pairwise equal charge STU black hole space-time. While in the latter case the minimally coupled Dirac equation is not separable, the introduction of a specific torsion term ensures the separability. The source of the torsion is the Kalb-Ramond field, which is an integral part of string theory, but further aspects of its properties and coupling to fermionic fields remain to be studied. To derive the torsion, two different approaches are used in conformally related frames, showing that the torsion is not unique. The correspondingly modified Dirac equations in the Einstein and string frames are shown to be separable. Furthermore, the massless radial and angular wave equations are examined; they show close similarity with corresponding equations for the standard Kerr background. A generalization of the Teukolsky equation for the pairwise equal case is conjectured. We also briefly analyze a technically sophisticated radial equation in the massive case.

Dirac equation↗

NASA NTRS1973

Numerical experiments in climate stability.

Two semiempirical climate models are formulated as time-dependent problems in order to study the stability of their asymptotic steady state equilibrium solutions to perturbations in internal (in this case, initial) conditions. For fixed external conditions the time-dependent versions of these models exhibit fully transitive behavior to positive perturbations in initial conditions, slight intransivity to negative perturbations up to -18 K, and an ice-covered earth regime is obtained for extremely large negative perturbations in initial temperatures (below -18 K). The parametrization found to be most critical in these models is the albedo-temperature coupling, especially in tropical regions. The temperature albedo parametrization from an intransitive climate model used in a time-dependent form of one of the preceding models is found to be highly sensitive to negative perturbations in both internal and external conditions. Numerical experimentation with these semiempirical models shows the important role of the tropics in maintaining the equilibrium climate and suggests that the radiation balance in equatorial latitudes might have a controlling influence on the equatorward extent of the polar ice cap.

Schneider, S. H.↗

NASA NTRS1976

Multiple scattering in planetary atmospheres

Certain simple procedures for solving radiative transfer problems in planetary atmospheres are reviewed: the similarity relations relating isotropic to anisotropic scattering, asymptotic results relating known solutions for semi-infinite layers to desired solutions for optical depths, and expansions relating known solutions for conservative scattering to desired solutions for the nearly conservative case. The complications introduced by atmospheric inhomogeneity, surface reflection, and spectral features are also discussed.

Irvine, W. M.↗

NASA NTRS1977

On the hydrodynamic expansion of a relativistic gas

An exact solution to the one-dimensional relativistic Euler equations is presented which is not of a self-similar nature and which satisfies certain initial conditions found to be appropriate in describing high-energy cosmic-ray collisions and models of radio galaxies. A full numerical analysis is performed to demonstrate the temporal behavior of all relevant thermodynamic quantities at any given point in the expansion. The results indicate that the motion is self-similar at the very beginning of free expansion, that this description ceases to be valid when the entire bulk of fluid is in motion, and that with the passage of time, the nonsimilar solution prevails over more and more of the fluid except for a small region near the leading edge where the fluid borders on the vacuum. The very bulk of the expansion is shown to become nonsimilar because the leading edge of the self-similar solution moves into the vacuum at the speed of light while that of the nonsimilar solution approaches c only asymptotically. The relativistic expansion of a gaseous disk is analyzed as an example.

Canuto, V.↗

NASA NTRS1977

On the gravitational potential and field anomalies due to thin mass layers

The gravitational potential and field anomalies for thin mass layers are derived using the technique of matched asymptotic expansions. An inner solution is obtained using an expansion in powers of the thickness and it is shown that the outer solution is given by a surface distribution of mass sources and dipoles. Coefficients are evaluated by matching the inner expansion of the outer solution with the outer expansion of the inner solution. The leading term in the inner expansion for the normal gravitational field gives the Bouguer formula. The leading term in the expansion for the gravitational potential gives an expression for the perturbation to the geoid. The predictions given by this term are compared with measurements by satellite altimetry. The second-order terms in the expansion for the gravitational field are required to predict the gravity anomaly at a continental margin. The results are compared with observations.

Ockendon, J. R.↗

NASA NTRS1979

Numerical method for predicting flow characteristics and performance of nonaxisymmetric nozzles, theory

The theoretical foundation and formulation of a numerical method for predicting the viscous flowfield in and about isolated three dimensional nozzles of geometrically complex configuration are presented. High Reynolds number turbulent flows are of primary interest for any combination of subsonic, transonic, and supersonic flow conditions inside or outside the nozzle. An alternating-direction implicit (ADI) numerical technique is employed to integrate the unsteady Navier-Stokes equations until an asymptotic steady-state solution is reached. Boundary conditions are computed with an implicit technique compatible with the ADI technique employed at interior points of the flow region. The equations are formulated and solved in a boundary-conforming curvilinear coordinate system. The curvilinear coordinate system and computational grid is generated numerically as the solution to an elliptic boundary value problem. A method is developed that automatically adjusts the elliptic system so that the interior grid spacing is controlled directly by the a priori selection of the grid spacing on the boundaries of the flow region.

Thomas, P. D.↗

NASA NTRS1980

Evaluation of approximate methods for the prediction of noise shielding by airframe components

An evaluation of some approximate methods for the prediction of shielding of monochromatic sound and broadband noise by aircraft components is reported. Anechoic-chamber measurements of the shielding of a point source by various simple geometric shapes were made and the measured values compared with those calculated by the superposition of asymptotic closed-form solutions for the shielding by a semi-infinite plane barrier. The shields used in the measurements consisted of rectangular plates, a circular cylinder, and a rectangular plate attached to the cylinder to simulate a wing-body combination. The normalized frequency, defined as a product of the acoustic wave number and either the plate width or cylinder diameter, ranged from 4.6 to 114. Microphone traverses in front of the rectangular plates and cylinders generally showed a series of diffraction bands that matched those predicted by the approximate methods, except for differences in the magnitudes of the attenuation minima which can be attributed to experimental inaccuracies. The shielding of wing-body combinations was predicted by modifications of the approximations used for rectangular and cylindrical shielding. Although the approximations failed to predict diffraction patterns in certain regions, they did predict the average level of wing-body shielding with an average deviation of less than 3 dB.

Ahtye, W. F.↗

NASA NTRS1981

Fully implicit shock tracking

A tracking method for solving hyperbolic free boundary problems arising in fluid dynamics is presented. The principal feature of the new approach is that the implicit form is used to simultaneously treat both interior points and boundary conditions. In particular, the location of the free boundary (shock) surface is treated implicitly and coupled with all other unknowns. The method is presented in the context of unsteady one-dimensional flow in a variable area duct with an internal shock wave. The fully implicit method and other strategies for advancing the shock are compared for computing a steady solution via a time asymptotic approach. Issues regarding extension of the method to multiple dimensions are also discussed.

Bell, J. B.↗

NASA NTRS1984

Implementation of uniform perturbation method for potential flow past axisymmetric and two-dimensional bodies

The aerodynamic characteristics of potential flow past an axisymmetric slender body and a thin airfoil are calculated using a uniform perturbation analysis method. The method is based on the superposition of potentials of point singularities distributed inside the body. The strength distribution satisfies a linear integral equation by enforcing the flow tangency condition on the surface of the body. The complete uniform asymptotic expansion of its solution is obtained with respect to the slenderness ratio by modifying and adapting an existing technique. Results calculated by the perturbation analysis method are compared with the existing surface singularity panel method and some available analytical solutions for a number of cases under identical conditions. From these comparisons, it is found that the perturbation analysis method can provide quite accurate results for bodies with small slenderness ratio. The present method is much simpler and requires less memory and computation time than existing surface singularity panel methods of comparable accuracy.

Wong, T. C.↗

NASA NTRS1987

A hybrid perturbation Galerkin technique with applications to slender body theory

A two step hybrid perturbation-Galerkin method to solve a variety of applied mathematics problems which involve a small parameter is presented. The method consists of: (1) the use of a regular or singular perturbation method to determine the asymptotic expansion of the solution in terms of the small parameter; (2) construction of an approximate solution in the form of a sum of the perturbation coefficient functions multiplied by (unknown) amplitudes (gauge functions); and (3) the use of the classical Bubnov-Galerkin method to determine these amplitudes. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is applied to some singular perturbation problems in slender body theory. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the degree of applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

NASA NTRS1986

Enclosed Gas and Liquid with Nonuniform Heating from Above

Buoyancy-driven flows of gases above liquids in a common enclosure with nonuniform heating from above are studied via finite-difference solutions of the governing equations. Unsteady solutions are calculated, and steady-state solutions are sought as asymptotes. Grashof numbers between 10 to the 3rd and 10 to the 8th are examined, and multicellular circulatory flow structure is found at the higher Grashof numbers. Convective transport dominates for higher Grashof numbers, while conductive transport is the primary mechanism at the lower Grashof numbers. Surface tension has a major effect upon the gas flow field only at lower Grashof numbers but, since conduction dominates there, it does not significantly affect transport.

Aggarwal, S. K.↗