Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adiabatic approximation”

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 235 records · Page 13

Energy and pitch angle distributions for auroral ions using the current sheet acceleration model

Using a dipole-plus-tail magnetic field model, H(+), He(+2), and O(+6) ions are followed numerically, backward in time, from an output plane perpendicular to the axis of the geomagnetic tail, to their point of entrance to the magnetosphere as solar wind particles in the magnetosheath. An adiabatic or guiding center approximation is used in regions where the particles do not interact directly with the current sheet. A Maxwellian distribution with bulk flow is assumed for solar wind particles in the magnetosheath. Bulk velocity, density, and temperature along the magnetopause are taken from the fluid calculations of Spreiter. Using Liouville's theorem, and varying initial conditions at the output plane, the distribution function is found as a function of energy and pitch angle at the output plane. These results are then mapped to the auroral ionosphere using guiding center theory.

Jaeger, E. F.↗

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Guiding-center simulations of stormtime transport of ring-current and radiation-belt ions having first adiabatic invariants mu is approximately greater than 15 MeV/G (E is approximately greater than 165 keV at L is approximately 3) are surprisingly well described (typically within a factor of approximately less than 4) by the quasilinear theory of radial diffusion. This holds even for the case of an individual model storm characterized by substorm-associated impulses in the convection electric field, provided that the actual spectrum of the electric field is incorporated in the quasilinear theory. Correction of the quasilinear diffusion coefficient D(sub LL)(sup ql) for drift-resonance broadening (so as to define D(sub LL)(sup ql)) reduced the typical discrepancy with the diffusion coefficients D(sub LL)(sup sim) deduced from guiding-center simulations of representative-particle trajectories to a factor of approximately 3. The typical discrepancy was reduced to a factor of approximately 1.4 by averaging D(sub LL)(sup sim), D(sub LL)(sup ql), and D(sub LL)(sup rb) over an ensemble of model storms characterized by different (but statistically equivalent) sets of substorm-onset times.

Chen, Margaret W.↗

Asymptotic theory of intermediate- and high-degree solar acoustic oscillations

A second-order asymptotic approximation is developed for adiabatic nonradial p-modes of a spherically symmetric star. The exact solutions of adiabatic oscillations are assumed in the outermost layers, where the asymptotic description becomes invalid, which results in a eigenfrequency equation with model-dependent surface phase shift. For lower degree modes, the phase shift is a function of frequency alone; for high-degree modes, its dependence on the degree is explicitly taken into account.

Brodsky, M.↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

Transition Delay via Vortex Generators in a Hypersonic Boundary Layer at Flight Conditions

The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.

Incompressible boundary layer↗

EFFECT OF ADIABATIC FOCUSING ON THE DIFFUSION OF ENERGETIC CHARGED PARTICLES

An approximate analytic description of a diffusion coefficient, including the effect of adiabatic focusing, has been developed. This description is formulated with the aid of stochastic differential equations and the steady perturbation solution of the Fokker–Plank transport equation. The analytical formula is based on three important assumptions. First, the pitch-angle diffusion coefficient is set to be separable from the spatial coordinate and the pitch-angle cosine. Second, the spatial dependence of the ratio between the mean-free path and focusing length is assumed to be weak. Third, the pitch-angle distribution relaxes quickly to a steady state. The new analytic formula could be applied to calculate the spatial diffusion coefficient in the interplanetary and interstellar space.

79 ASTRONOMY AND ASTROPHYSICS↗

Extended adiabatic blast waves and a model of the soft X-ray background

An analytical approximation is generated which follows the development of an adiabatic spherical blast wave in a homogeneous ambient medium of finite pressure. An analytical approximation is also presented for the electron temperature distribution resulting from coulomb collisional heating. The dynamical, thermal, ionization, and spectral structures are calculated for blast waves of energy E sub 0 = 5 x 10 to the 50th power ergs in a hot low-density interstellar environment. A formula is presented for estimating the luminosity evolution of such explosions. The B and C bands of the soft X-ray background, it is shown, are reproduced by such a model explosion if the ambient density is about .000004 cm, the blast radius is roughly 100 pc, and the solar system is located inside the shocked region. Evolution in a pre-existing cavity with a strong density gradient may, it is suggested, remove both the M band and OVI discrepancies.

Cox, D. P.↗

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy–Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn–Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree–Fock and Bethe ansatz local density approximation) to show their successes and limitations. Here, we also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Full-dimensional quantum stereodynamics of the non-adiabatic quenching of OH(A2Σ+) by H2

Abstract The Born–Oppenheimer approximation, assuming separable nuclear and electronic motion, is widely adopted for characterizing chemical reactions in a single electronic state. However, the breakdown of the Born–Oppenheimer approximation is omnipresent in chemistry, and a detailed understanding of the non-adiabatic dynamics is still incomplete. Here we investigate the non-adiabatic quenching of electronically excited OH(A 2 Σ + ) molecules by H 2 molecules using full-dimensional quantum dynamics calculations for zero total nuclear angular momentum using a high-quality diabatic-potential-energy matrix. Good agreement with experimental observations is found for the OH(X 2 Π) ro-vibrational distribution, and the non-adiabatic dynamics are shown to be controlled by stereodynamics, namely the relative orientation of the two reactants. The uncovering of a major (in)elastic channel, neglected in a previous analysis but confirmed by a recent experiment, resolves a long-standing experiment–theory disagreement concerning the branching ratio of the two electronic quenching channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of a qubit in a high-impedance transmission line from a bath perspective

We investigate the quantum dynamics of a generic model of light-matter interaction in the context of high-impedance waveguides, focusing on the behavior of the photonic states generated in the waveguide. The model treated consists simply of a two-level system coupled to a bosonic bath (the Ohmic spin-boson model). Quantum quenches as well as scattering of an incident coherent pulse are studied using two complementary methods. First, we develop an approximate ansatz for the electromagnetic waves based on a single multimode coherent state wave function; formally, this approach combines in a single framework ideas from adiabatic renormalization, the Born-Markov approximation, and input-output theory. Second, we present numerically exact results for scattering of a weak intensity pulse by using numerical renormalization group (NRG) calculations. NRG provides a benchmark for any linear response property throughout the ultrastrong-coupling regime. We find that in a sudden quantum quench, the coherent state approach produces physical artifacts, such as improper relaxation to the steady state. These previously unnoticed problems are related to the simplified form of the ansatz that generates spurious correlations within the bath. In the scattering problem, NRG is used to find the transmission and reflection of a single photon, as well as the inelastic scattering of that single photon. Simple analytical formulas are established and tested against the NRG data that predict quantitatively the transport coefficients for up to moderate environmental impedance. These formulas resolve pending issues regarding the presence of inelastic losses in the spin-boson model near absorption resonances, and could be used for comparison to experiments in Josephson waveguide quantum electrodynamics. Finally, the scattering results using the coherent state wave-function approach are compared favorably to the NRG results for very weak incident intensity. Finally, we end our study by presenting results at higher power where the response of the system is nonlinear.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

INTO THE LAIR: GRAVITATIONAL-WAVE SIGNATURES OF DARK MATTER

The nature and properties of dark matter (DM) are both outstanding issues in physics. Besides clustering in halos, the universal character of gravity implies that self-gravitating compact DM configurations-predicted by various models-might be spread throughout the universe. Their astrophysical signature can be used to probe fundamental particle physics, or to test alternative descriptions of compact objects in active galactic nuclei. Here, we discuss the most promising dissection tool of such configurations: the inspiral of a compact stellar-size object and consequent gravitational-wave (GW) emission. The inward motion of this ''test probe'' encodes unique information about the nature of the supermassive configuration. When the probe travels through some compact region we show, within a Newtonian approximation, that the quasi-adiabatic inspiral is mainly driven by DM accretion and by dynamical friction, rather than by radiation reaction. When accretion dominates, the frequency and amplitude of the GW signal produced during the latest stages of the inspiral are nearly constant. In the exterior region we study a model in which the inspiral is driven by GW and scalar-wave emission, described at a fully relativistic level. Resonances in the energy flux appear whenever the orbital frequency matches the effective mass of the DM particle, corresponding to the excitation of the central object's quasinormal frequencies. Unexpectedly, these resonances can lead to large dephasing with respect to standard inspiral templates, to such an extent as to prevent detection with matched filtering techniques. We discuss some observational consequences of these effects for GW detection.

79 ASTRONOMY AND ASTROPHYSICS↗

PHYSICAL-CONSTRAINT-PRESERVING CENTRAL DISCONTINUOUS GALERKIN METHODS FOR SPECIAL RELATIVISTIC HYDRODYNAMICS WITH A GENERAL EQUATION OF STATE

The ideal gas equation of state (EOS) with a constant adiabatic index is a poor approximation for most relativistic astrophysical flows, although it is commonly used in relativistic hydrodynamics (RHD). This paper develops high-order accurate, physical-constraints-preserving (PCP), central, discontinuous Galerkin (DG) methods for the one- and two-dimensional special RHD equations with a general EOS. It is built on our theoretical analysis of the admissible states for RHD and the PCP limiting procedure that enforce the admissibility of central DG solutions. The convexity, scaling invariance, orthogonal invariance, and Lax–Friedrichs splitting property of the admissible state set are first proved with the aid of its equivalent form. Then, the high-order central DG methods with the PCP limiting procedure and strong stability-preserving time discretization are proved, to preserve the positivity of the density, pressure, specific internal energy, and the bound of the fluid velocity, maintain high-order accuracy, and be L {sup 1}-stable. The accuracy, robustness, and effectiveness of the proposed methods are demonstrated by several 1D and 2D numerical examples involving large Lorentz factor, strong discontinuities, or low density/pressure, etc.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure and circulation of the Venus atmosphere

The Pioneer Venus data relevant to the dynamics and thermodynamics of the atmosphere is summarized and interpreted. On the day side there is a thermosphere in which temperatures increase with height to an exospheric temperature of about 300 K. On the night side there is a cryosphere in which temperatures decrease with height to an exospheric temperature of about 100 K. The atmosphere is stratified stably from the highest altitudes down to about 28 km except for a layer in the clouds between about 50 and 55 km which is nearly adiabatic. Horizontal thermal contrasts are approximately 1 to 2% in the deep atmosphere and 100% in the upper atmosphere. The temperatures generally decrease with latitude at and below the clouds on constant pressure surfaces. Above the clouds there is a reversed zonally averaged latitudinal temperature gradient. The dominant circulation of the atmosphere above the lowest one or two scale heights is a zonal retrograde motion with 100 m/s winds at 60 km altitude. There is also a superrotation at altitudes of 150 km and above.

Schubert, G.↗

Seasonal changes in the structure of the stratospheric circulation

The seasonal evolution of the stratospheric circulation in both hemispheres is studied using a coordinate independent diagnostic tied to the vortex, wherever it lies, rather than to latitude circles. Zonally averaged quantities are highly variable in the winter stratosphere, as the vortex moves on and off the pole, masking systematic changes in the structure of the vortex. The quantities shown here are the areas A(Q), where Ertel's potential vorticity is greater than a specified value Q, for various Q. In these calculations, Q is evaluated on the 850K isentropic surface which lies near 10 mb in the middle of the stratosphere. The advantage of using contours of Q (rather than contours of geopotential height, say) is that an area enclosed by such contours is approximately constant for quasi nondivergent, adiabatic, inviscid motion. Changes in area will occur if any of the above restrictions on the motion are invalid, and also if scales of motions are generated below those that can be adequately resolved. The evolution of area shows the systematic change in the average structure of the vortex. At a given time, the rate of change of Q with A(Q) gives a measure of the average gradient of Q. The gradient of Q is relevant for the propagation of large-scale disturbances into the stratosphere. Experience has also shown that winds tend to be strong (weak) when the gradient of Q is strong (weak).

Oneill, A.↗

The Effect of Ultrapolish on a Transonic Axial Rotor

Back-to-back testing has been done using NASA fan rotor 67 in the Glenn Research Center W8 Axial Compressor Test Facility. The rotor was baseline tested with a normal industrial RMS surface finish of 0.5-0.6 m (20-24 microinches) at 60, 80 and 100% of design speed. At design speed the tip relative Mach number was 1.38. The blades were then removed from the facility and ultrapolished to a surface finish of 0.125 m (5 microinch) or less and retested. At 100% speed near the design point, the ultrapolished blades showed approximately 0.3 - 0.5% increase in adiabatic efficiency. The difference was greater near maximum flow. Due to increased relative measurement error at 60 and 80% speed, the performance difference between the normal and ultrapolished blades was indeterminate at these speeds.

Roberts, William B.↗