Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave equation”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Ring Current Ion Coupling with Electromagnetic Ion Cyclotron Waves

A new ring current global model has been developed for the first time that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes wave evolution of electromagnetic ion cyclotron waves (EMIC). The coupled model is able to simulate, for the first time self-consistently calculated RC ion kinetic and evolution of EMIC waves that propagate along geomagnetic field lines and reflect from the ionosphere. Ionospheric properties affect the reflection index through the integral Pedersen and Hall coductivities. The structure and dynamics of the ring current proton precipitating flux regions, intensities of EMIC, global RC energy balance, and some other parameters will be studied in detail for the selected geomagnetic storms. The space whether aspects of RC modelling and comparison with the data will also be discussed.

Khazanov, George V.↗
NASA NTRS1995

Symmetries of the TDNLS equations for weakly nonlinear dispersive MHD waves

In this paper we consider the symmetries and conservation laws for the TDNLS equations derived by Hada (1993) and Brio, Hunter and Johnson, to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a(g)(exp 2) = V(A)(exp 2) where a(g) is the gas sound speed and V(A) is the Alfven speed. We discuss Lagrangian and Hamiltonian formulations, and similarity solutions for the equations.

Webb, G. M.↗
NASA NTRS1981

A simple model of the Lagrangian-mean flow produced by dissipating planetary waves

A simple equation for the Lagrangian-mean flow induced by damped planetary waves is derived. The flow computed for stationary planetary waves of a beta-plane is found to be generally poleward and downward during winter and appears to be about twice as strong as the diabatic circulation in the lower stratosphere. An important factor in determining the high-latitude Lagrangian-mean flow field is the subtropical jet stream which blocks planetary wave propagation toward the equatorial regions. Computations using two types of Lagrangian-mean boundary conditions at the surface show that incorrect orographic forcing distorts the Lagrangian-mean flow up to three scale heights or more above ground.

Schoeberl, M. R.↗
DOE OSTI2016

Multiscale real-space quantum-mechanical tight-binding calculations of electronic structure in crystals with defects using perfectly matched layers

Graphical abstract: - Abstract: We consider the scattering of incident plane-wave electrons from a defect in a crystal modeled by the time-harmonic Schrödinger equation. While the defect potential is localized, the far-field potential is periodic, unlike standard free-space scattering problems. Previous work on the Schrödinger equation has been almost entirely in free-space conditions; a few works on crystals have been in one-dimension. We construct absorbing boundary conditions for this problem using perfectly matched layers in a tight-binding formulation. Using the example of a point defect in graphene, we examine the efficiency and convergence of the proposed absorbing boundary condition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

A planar oblique wave model for the Euler equations

A numerical method in which any transition between neighboring states in a multidimensional flowfield is reconstructed by a pattern of oblique waves. The wave pattern is obtained from a superposition of model waves which arise in a local linearization based on Rankine-Hugoniot averaging. The strengths and orientations of these waves are derived from a minimum path length rule applied in state space. Two test cases with comparisons to other methods in planar flow problems are presented to demonstrate the strengths of the method, and to point out its weaknesses.

Parpia, Ijaz H.↗
NASA NTRS1977

Equations of state and impact-induced shock-wave attenuation on the moon

Equation of state formulations are considered in a framework that permits comparison with one-dimensional impedance match solutions. The problem considered is the peak pressures attained along the impact symmetry axis when a sphere impacts with a half-space. The regimes of melting and vaporization - in particular, incipiently melted, completely melted, incipiently vaporized, and completely vaporized states - are examined, and the pressures at which critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite are considered. A means of representing the spatial attenuation of shock pressure along the impact axis by two regimes is introduced, and results for the near-field and far-field regime are presented. It is thought that the treatment can be used to obtain quantitative bounds on the impact velocity of the meteorite.

Ahrens, T. J.↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: Recent advances

The basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of the wave-kinetic theory. Recently, this framework has been reexamined from first principles, which has led to more accurate yet still tractable “improved” wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow “curvature” (the second radial derivative of the flow velocity) on the dynamics and stability of drift waves and zonal flows. Here, we overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Wave normal direction and spectral properties of whistler mode hiss observed on the DE 1 satellite

Hiss is represented by a field distribution function in order to investigate magnetospheric hiss as a spatially and temporally enduring phenomenon. The study takes into account the whistler mode relationships and the linear and spin motion of the satellite. Hiss signals received on September 23, 1983 by the DE-1 electric and magnetic field antennas are analyzed. A wave normal angle of 60 + or - 5 deg with respect to the local geomagnetic field is found near the geomagnetic equator, and wave normal directions from 30-80 deg with respect to the local geomagnetic field are found away from the equator.

Sonwalkar, Vikas S.↗
DOE OSTI2024

Superluminal matter waves

The Dirac equation has resided among the greatest successes of modern physics since its emergence as the first quantum mechanical theory fully compatible with special relativity. This compatibility ensures that the expectation value of the velocity is less than the vacuum speed of light. Here, we show that the Dirac equation admits free-particle solutions where the peak amplitude of the wave function can travel at any velocity, including those exceeding the vacuum speed of light, despite having a subluminal velocity expectation value. The solutions are constructed by superposing basis functions with correlations in momentum space. These arbitrary velocity wave functions feature a near-constant profile and may impact quantum mechanical processes that are sensitive to the local value of the probability density as opposed to expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Simultaneous measurements of waves and precipitating electrons near the equator in the outer radiation belt

An investigation of wave-particle interactions is made using several simultaneous electron and wave measurements performed at near-equatorial positions from the Combined Release and Radiation Effects Satellite (CRRES) satellite. Bursts of electron precipitation were observed, most frequently at local times near dawn. Examples of bursts are presented in which the fluxes of the precipitating electrons and the wave intensities are correlated with coefficients as high as 0.7. During bursts the frequencies of the enhanced waves spanned a wide range from 311 Hz to 3.11 kHz, and the energies of the enhanced electrons were in the range 1.7 keV to 288 keV. The changes of the precipitating fluxes were generally less pronounced at the lowest energies. On the basis of electron-cyclotron resonant calculations using the cold plasma densities and ambient magnetic fields taken from the CRRES measurements it was found that the wave frequencies and precipitating electron energies were generally consistent with those expected from electron resonance with parallel propagating whistler waves. The electron data of principal concern here were acquired in and about the loss cone with narrow angular resolution spectrometers covering the energy range 340 eV to 5 MeV. The wave data included electric field measurements spanning frequencies from 5 Hz to 400 kHz and magnetic field measurements from 5 Hz to 10 kHz.

Imhof, W. L.↗
DOE OSTI2015

Kinetic description of linear wave propagation in inhomogeneous, nonstationary, anisotropic, weakly magnetized, and collisional plasma

This paper addresses the linear propagation of an electron wave in a plasma whose distribution function, at zero order in the wave amplitude, may be chosen arbitrarily, provided that it is not strongly peaked at the wave phase velocity, and that it varies very little over one wave period and one wavelength. Then, from first principles is derived an equation for the wave action density that allows for Landau damping, whose rate is calculated at first order in the variations of the wave number and frequency. Moreover, the effect of collisions is accounted for in a way that adapts to any choice for the collision operator in Boltzmann equation. The wave may also be externally driven, so that the results presented here apply to stimulated Raman scattering.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗
DOE OSTI2016

Solitary waves in the nonlinear Dirac equation in the presence of external driving forces

In this paper, we consider the nonlinear Dirac (NLD) equation in (1 + 1) dimensions with scalar–scalar self interaction g 2 /κ + 1 (Ψ¯Ψ) κ + 1 in the presence of external forces as well as damping of the form f(x) - iμγ 0 Ψ, where both f and Ψ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. For κ = 1 we solve explicitly the CC equations of the variational approximation for slow moving solitary waves in a constant external force without damping and find reasonable agreement with solving numerically the CC equations. Finally, we then compare the results of the variational approximation with no damping with numerical simulations of the NLD equation for κ = 1, when the components of the external force are of the form f j = r j exp(–iΚx) and again find agreement if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

The effect of latent heat release on synoptic-to-planetary wave interactions and its implication for satellite observations: Theoretical modeling

The project objectives are to develop process models to investigate the interaction of planetary and synoptic-scale waves including the effects of latent heat release (precipitation), nonlinear dynamics, physical and boundary-layer processes, and large-scale topography; to determine the importance of latent heat release for temporal variability and time-mean behavior of planetary and synoptic-scale waves; to compare the model results with available observations of planetary and synoptic wave variability; and to assess the implications of the results for monitoring precipitation in oceanic-storm tracks by satellite observing systems. Researchers have utilized two different models for this project: a two-level quasi-geostrophic model to study intraseasonal variability, anomalous circulations and the seasonal cycle, and a 10-level, multi-wave primitive equation model to validate the two-level Q-G model and examine effects of convection, surface processes, and spherical geometry. It explicitly resolves several planetary and synoptic waves and includes specific humidity (as a predicted variable), moist convection, and large-scale precipitation. In the past year researchers have concentrated on experiments with the multi-level primitive equation model. The dynamical part of that model is similar to the spectral model used by the National Meteorological Center for medium-range forecasts. The model includes parameterizations of large-scale condensation and moist convection. To test the validity of results regarding the influence of convective precipitation, researchers can use either one of two different convective schemes in the model, a Kuo convective scheme or a modified Arakawa-Schubert scheme which includes downdrafts. By choosing one or the other scheme, they can evaluate the impact of the convective parameterization on the circulation. In the past year researchers performed a variety of initial-value experiments with the primitive-equation model. Using initial conditions typical of climatological winter conditions, they examined the behavior of synoptic and planetary waves growing in moist and dry environments. Surface conditions were representative of a zonally averaged ocean. They found that moist convection associated with baroclinic wave development was confined to the subtropics.

Branscome, Lee E.↗
DOE OSTI2014

Numerical solution of the time dependent Ginzburg-Landau equations for mixed (d + s)-wave superconductors

The time-dependent Ginzburg-Landau formalism for (d + s)-wave superconductors and their representation using auxiliary fields is investigated. By using the link variable method, we then develop suitable discretization of these equations. Numerical simulations are carried out for a mesoscopic superconductor in a homogeneous perpendicular magnetic field which revealed peculiar vortex states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

The thermal-wave model: A Schroedinger-like equation for charged particle beam dynamics

We review some results on longitudinal beam dynamics obtained in the framework of the Thermal Wave Model (TWM). In this model, which has recently shown the capability to describe both longitudinal and transverse dynamics of charged particle beams, the beam dynamics is ruled by Schroedinger-like equations for the beam wave functions, whose squared modulus is proportional to the beam density profile. Remarkably, the role of the Planck constant is played by a diffractive constant epsilon, the emittance, which has a thermal nature.

Fedele, Renato↗