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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Umbral oscillations as resonant modes of magneto-atmospheric waves

Umbral oscillations in sunspots are identified as a resonant response of the umbral atmosphere to forcing by oscillatory convection in the subphotosphere. The full, linearized equations for magnetoatmospheric waves are solved numerically for a detailed model of the umbral atmosphere, for both forced and free oscillations. Resonant 'fast' modes are found, the lowest mode having a period of 153 s, typical of umbral oscillations. A comparison is made with a similar analysis by Uchida and Sakurai (1975), who calculated resonant modes using an approximate ('quasi-Alfven') form of the wave equations. Whereas both analyses give an appropriate value for the period of oscillation, several new features of the motion follow from the full equations. The resonant modes are due to upward reflection in the subphotosphere (due to increasing sound speed) and downward reflection in the photosphere and low chromosphere (due to increasing Alfven speed); downward reflection at the chromosphere-corona transition is unimportant for these modes.

Scheuer, M. A.↗
NASA NTRS1972

Shock-wave structure using nonlinear model Boltzmann equations.

The structure of strong plane shock waves in a perfect monatomic gas was studied using four nonlinear models of the Boltzmann equation. The models involved the use of a simplified collision operator with velocity-independent collision frequency, in place of the complicated Boltzmann collision operator. The models employed were the BGK and ellipsoidal models developed by earlier authors, and the polynomial and trimodal gain function models developed during the work. An exact set of moment equations was derived for the density, velocity, temperature, viscous stress, and heat flux within the shock. This set was reduced to a pair of coupled nonlinear integral equations and solved using specially adapted numerical techniques. A new and simple Gauss-Seidel iteration was developed during the work and found to be as efficient as the best earlier iteration methods.

Segal, B. M.↗
NASA NTRS1984

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗
NASA NTRS1984

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗
NASA NTRS1989

Nonlinear evolution of the auroral electron beam

The nonlinear spatial evolution, from the source to the atmosphere, of the auroral electron beam and the beam-generated electrostatic whistler noise was studied, calculating changes in beam parameters from equations for the conservation of total particle and wave energy and momentum flux density. Wave power fluxes were calculated by numerically integrating the wave kinetic equations, and the levels of beam-generated noise were determined by using thermal levels of Cerenkov radiation as a source. It was found that beam parameters evolve on ionospheric scale lengths, and their positive slope feature in velocity space is maintained over altitudes measured in thousands of kilometers of altitude, even though they can generate wave energy density fluxes sufficient to modify the ionospheric density profile.

Maggs, James E.↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Stochastic partial differential equations in turbulence related problems

The theory of stochastic partial differential equations (PDEs) and problems relating to turbulence are discussed by employing the theories of Brownian motion and diffusion in infinite dimensions, functional differential equations, and functional integration. Relevant results in probablistic analysis, especially Gaussian measures in function spaces and the theory of stochastic PDEs of Ito type, are taken into account. Linear stochastic PDEs are analyzed through linearized Navier-Stokes equations with a random forcing. Stochastic equations for waves in random media as well as model equations in turbulent transport theory are considered. Markovian models in fully developed turbulence are discussed from a stochastic equation viewpoint.

Chow, P.-L.↗
NASA NTRS1989

The sensitivity of stationary waves to variations in the basic state zonal flow

A linear, primitive equation stationary wave model having high vertical and meridional resolution is used to examine the sensitivity of orographically forced (primarily by Himalayas) stationary waves at middle and high latitudes to variations in the basic state zonal wind distribution. We find relatively little sensitivity to the winds in high latitudes, but remarkable sensitivity to small variations in the subtropical jet. Fluctuations well within the range of observed variability in the jet can lead to large variations in the stationary waves of the high latitude stratosphere, and to large changes even in tropospheric stationary waves. Implications for both sudden warmings and large-scale weather are discussed.

Nigam, Sumant↗
NASA NTRS2001

A Self-Consistent Model of the Interacting Ring Current Ions with Electromagnetic ICWs

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of ring current ions and ion cyclotron waves in a quasilinear approach. These two equations were solved on a global scale under non steady-state conditions during the May 2-5, 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the wave active zones at three time cuts around initial, main, and late recovery phases of the May 4, 1998 storm phase are presented and discussed in detail. Comparisons of the model wave-ion data with the Polar/HYDRA and Polar/MFE instruments results are presented..

Khazanov, G. V.↗
DOE OSTI2014

Effect of strongly coupled plasma on photoionization cross section

The effect of strongly coupled plasma on the ground state photoionization cross section is studied. In the non relativistic dipole approximation, cross section is evaluated from bound-free transition matrix element. The bound and free state wave functions are obtained by solving the radial Schrodinger equation with appropriate plasma potential. We have used ion sphere potential (ISP) to incorporate the plasma effects in atomic structure calculation. This potential includes the effect of static plasma screening on nuclear charge as well as the effect of confinement due to the neighbouring ions. With ISP, the radial equation is solved using Shooting method approach for hydrogen like ions (Li{sup +2}, C{sup +5}, Al{sup +12}) and lithium like ions (C{sup +3}, O{sup +5}). The effect of strong screening and confinement is manifested as confinement resonances near the ionization threshold for both kinds of ions. The confinement resonances are very much dependent on the edge of the confining potential and die out as the plasma density is increased. Plasma effect also results in appearance of Cooper minimum in lithium like ions, which was not present in case of free lithium like ions. With increasing density the position of Cooper minimum shifts towards higher photoelectron energy. The same behaviour is also true for weakly coupled plasma where plasma effect is modelled by Debye-Huckel potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Waves guided by a thin viscoelastic layer between elastic solids

The propagation of ultrasonic waves guided along a viscoelastic layer which separates two dissimilar elastic solid half spaces is described. For the limiting cases of rigid and soft bonding by a layer thin compared to acoustic wavelength, Stoneley and independent Rayleigh wave characteristic equations, respectively, result. For combinations of layer rigidity and wave length that correspond to guided wave propagation, external cyclic loading of the layer produces mechanical hysteresis and a resulting dynamic change in bond properties. Measurements of velocity hysteresis in an aluminum polymer adhesive aluminum system are described.

Claus, R. O.↗
NASA NTRS1986

Vector solution for the mean electromagnetic fields in a layer of random particles

The mean electromagnetic fields are found in a layer of randomly oriented particles lying over a half space. A matrix-dyadic formulation of Maxwell's equations is employed in conjunction with the Foldy-Lax approximation to obtain equations for the mean fields. A two variable perturbation procedure, valid in the limit of small fractional volume, is then used to derive uncoupled equations for the slowly varying amplitudes of the mean wave. These equations are solved to obtain explicit expressions for the mean electromagnetic fields in the slab region in the general case of arbitrarily oriented particles and arbitrary polarization of the incident radiation. Numerical examples are given for the application to remote sensing of vegetation.

Lang, R. H.↗
DOE OSTI2019

Matter wave soliton solutions for driven Gross–Pitaevskii equation with distributed coefficients

Highlights: •Bright/dark matter wave soliton solutions are found for distributed coefficient driven GP equation. •Dynamics of solitons has been revealed for different type of experimentally realizable nonlinear interactions. •Effect of source term on the obtained soliton solutions has been shown. -- Abstract: We have obtained the matter wave soliton solution for the driven nonautonomous GP equation with quadratic–cubic nonlinearities. We have revealed their interesting dynamical features for the experimentally realizable form of nonlinear interactions. Additionally, we have described that by suitably choosing some system parameters, we can control their characteristics such as amplitude and width etc. We have found that our results/solutions are stable and can be applied for future applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SUPRATHERMAL SOLAR WIND ELECTRONS AND LANGMUIR TURBULENCE

The steady-state model recently put forth for the solar wind electron velocity distribution function during quiet time conditions, was originally composed of three population electrons (core, halo, and superhalo) with the core remaining nonresonant with any plasma waves while the halo and superhalo separately maintained steady-state resonance with whistler- and Langmuir-frequency range fluctuations, respectively. However, a recent paper demonstrates that whistler-range fluctuations in fact have no significant contribution. The present paper represents a consummation of the model in that a self-consistent model of the suprathermal electron population, which encompasses both the halo and the superhalo, is constructed solely on the basis of the Langmuir fluctuation spectrum. Numerical solutions to steady-state particle and wave kinetic equations are obtained on the basis of an initial trial electron distribution and Langmuir wave spectrum. Such a finding offers a self-consistent explanation for the observed steady-state electron distribution in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Shock wave structure using nonlinear model Boltzmann equations

The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.

Segal, Ben Maurice↗
DOE OSTI2013

Standing waves for coupled nonlinear Schrödinger equations with decaying potentials

We study the following singularly perturbed problem for a coupled nonlinear Schrödinger system which arises in Bose-Einstein condensate: −ε{sup 2}Δu + a(x)u = μ{sub 1}u{sup 3} + βuv{sup 2} and −ε{sup 2}Δv + b(x)v = μ{sub 2}v{sup 3} + βu{sup 2}v in R{sup 3} with u, v > 0 and u(x), v(x) → 0 as |x| → ∞. Here, a, b are non-negative continuous potentials, and μ{sub 1}, μ{sub 2} > 0. We consider the case where the coupling constant β > 0 is relatively large. Then for sufficiently small ε > 0, we obtain positive solutions of this system which concentrate around local minima of the potentials as ε → 0. The novelty is that the potentials a and b may vanish at someplace and decay to 0 at infinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗