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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Explicit Weierstrass traveling wave solutions and bifurcation analysis for dissipative Zakharov–Kuznetsov modified equal width equation

We consider the dissipative Zakharov–Kuznetsov modified equal width (ZK-MEW) equation and discuss the effect of the dissipation on the existence and nature of traveling wave solutions of the equation. We use Lyapunov function and dynamical system theory in order to show that when viscosity term is added to the ZK-MEW equation, yet there exists still bounded traveling wave solutions in certain regions. Subsequently, we obtain general solution of the ZK-MEW equation in the presence and absence of viscosity in terms of Weirstrass ℘ functions and Jacobi elliptic functions. We also derive a new form of kink-type solution by exploring a factorization method based on functional transformation and the Abel’s first order nonlinear equation. Finally, we use the phase plane analysis and examine the stability of the viscous waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Dust-acoustic Gardner solitons and double layers in dusty plasmas with nonthermally distributed ions of two distinct temperatures

A rigorous theoretical investigation has been performed on dust-acoustic (DA) solitary structures in an unmagnetized dusty plasma, consisting of negatively charged mobile dust grains, Boltzmann distributed electrons, and nonthermally distributed ions of two distinct temperatures. The Korteweg-de Vries (K-dV), modified K-dV (mK-dV) and Gardner equations, and their solitary waves (SWs) and double layer (DL) (in case of Gardner equation) solutions are derived by using the reductive perturbation technique. The basic features of the DA Gardner solitons (GSs) and DLs are studied analytically as well as numerically. It has been observed that the GSs significantly differ from K-dV and mK-dV solitons, and only positive potential DLs exist in the system. It is also studied that two-temperature nonthermal ions significantly modify the nature and basic properties of the DA SWs. The present investigation can be very effective for understanding and studying the nonlinear characteristics of the DA waves in laboratory and space dusty plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Coherent distributions for the rigid rotator

Coherent solutions of the classical Liouville equation for the rigid rotator are presented as positive phase-space distributions localized on the Lagrangian submanifolds of Hamilton-Jacobi theory. These solutions become Wigner-type quasiprobability distributions by a formal discretization of the left-invariant vector fields from their Fourier transform in angular momentum. The results are consistent with the usual quantization of the anisotropic rotator, but the expected value of the Hamiltonian contains a finite “zero point” energy term. It is shown that during the time when a quasiprobability distribution evolves according to the Liouville equation, the related quantum wave function should satisfy the time-dependent Schrödinger equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

High Order Discontinuous Galerkin Method and Applications (Final Technical Report)

This project aims at the developments and improvements of high order accurate discontinuous Galerkin finite element methods for solving partial differential equations arising from DOE applications. Algorithm development, improvement, analysis, implementation and applications aspects are addressed. The partial differential equations under study involve solutions with shocks or other discontinuities, rapidly changing regions, or multiscale solutions. These include the so-called convection dominated cases, such as the hyperbolic equations (e.g. Euler equations of gas dynamics), parabolic equations with small viscosities (e.g. high Reynolds number Navier-Stokes equations), KdV type equations with small dispersion describing water waves (small third derivative terms), time-dependent bi-harmonic equations and equations with hyper-viscosity (small fourth or higher even derivative terms), problems with highly oscillatory coefficients, or problems for which the macro-scale PDE models may not be correct in localized regions and a micro-scale model must be used instead. Furthermore, we are interested in the long time solutions of such equations with waves traveling a long distance, such as the electromagnetic waves. The requirement on numerical methods is thus high order accuracy, low dissipation/dispersion errors over long time, and a clean and sharp discontinuity resolution. It would also be highly desirable to have the numerical methods to be flexible for complicated geometry and boundary conditions, and to be easy in an adaptive and massively parallel environment. The discontinuous Galerkin finite element method has this potential, however the design and analysis of the discontinuous Galerkin method for each specific class of partial differential equations would require detailed research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Quantum ring solitons and nonlocal effects in plasma wake field excitations

A theoretical investigation of the quantum transverse beam motion for a cold relativistic charged particle beam travelling in a cold, collisionless, strongly magnetized plasma is carried out. This is done by taking into account both the individual quantum nature of the beam particles (single-particle uncertainty relations and spin) and the self consistent interaction generated by the plasma wake field excitation. By adopting a fluid model of a strongly magnetized plasma, the analysis is carried out in the overdense regime (dilute beams) and in the long beam limit. It is shown that the quantum description of the collective transverse beam dynamics is provided by a pair of coupled nonlinear governing equations. It comprises a Poisson-like equation for the plasma wake potential (driven by the beam density) and a 2D spinorial Schroedinger equation for the wave function, whose squared modulus is proportional to the beam density, that is obtained in the Hartree's mean field approximation, after disregarding the exchange interactions. The analysis of this pair of equations, which in general exhibits a strong nonlocal character, is carried out analytically as well as numerically in both the linear and the nonlinear regimes, showing the existence of the quantum beam vortices in the form of Laguerre-Gauss modes and ring envelope solitons, respectively. In particular, when the relation between the plasma wake field response and the beam probability density is strictly local, the pair of the governing equations is reduced to the 2D Gross-Pitaevskii equation that allows one to establish the conditions for the self focusing and collapse. These conditions include the quantum nature of the beam particles. Finally, when the relation between the plasma wake field response and the beam probability density is moderately nonlocal, the above pair of equations permits to follow the spatio-temporal evolution of a quantum ring envelope soliton. Such a structure exhibits small or violent breathing, but it remains very stable for long time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Dynamical model for longitudinal wave functions in light-front holographic QCD

We construct a Schrödinger-like equation for the longitudinal wave function of a meson in the valence qq{sup -bar} sector, based on the ’t Hooft model for large-N two-dimensional QCD, and combine this with the usual transverse equation from light-front holographic QCD, to obtain a model for mesons with massive quarks. The computed wave functions are compared with the wave function ansatz of Brodsky and de Téramond and used to compute decay constants and parton distribution functions. The basis functions used to solve the longitudinal equation may be useful for more general calculations of meson states in QCD. -- Highlights: •Provide relativistic quark model based on light-front holographic QCD. •Incorporate dependence on quark mass. •Consistent with the Brodsky–de Téramond quark-wave-function ansatz. •Compute meson decay constants and parton distribution functions. •Illustrate use of basis functions that could be convenient for more general numerical calculations in light-front QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

On the existence and stability of electrostatic structures in non-Maxwellian electron-positron-ion plasmas

Electrostatic solitary waves in plasmas are the focus of many current studies of localized electrostatic disturbances in both laboratory and astrophysical plasmas. Here, an investigation of the nonlinear dynamics of plasma evolving in two dimensions, in the presence of excess superthermal background electrons and positrons, is undertaken. We investigate the effect of a magnetic field on weakly nonlinear ion acoustic waves. Deviation from the Maxwellian distribution is effectively modelled by the kappa model. A linear dispersion relation is derived, and a decrease in frequency and phase speed in both parallel and perpendicular modes can be seen, when the proportion of positrons to electrons increases. We show that ion acoustic solitary waves can be generated during the nonlinear evolution of a plasma fluid, and their nonlinear propagation is governed by a Zakharov-Kuznetsov (ZK) type equation. A multiple scales perturbation technique is used to derive the ZK equation. The solitary wave structures are dependent on the relation between the system parameters, specifically the superthermality of the system, the proportion of positron content, magnetic field strength, and the difference between electron and positron temperature. The parametric effect of these on electrostatic shock structures is investigated. In particular, we find that stronger superthermality leads to narrower excitations with smaller potential amplitudes. Increased positron concentration also suppresses both the amplitude and the width of solitary wave structures. However, the structures are only weakly affected by temperature differentials between electrons and positrons in our model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Linear analysis on the growth of non-spherical perturbations in supersonic accretion flows

We analyzed the growth of non-spherical perturbations in supersonic accretion flows. We have in mind an application to the post-bounce phase of core-collapse supernovae (CCSNe). Such non-spherical perturbations have been suggested by a series of papers by Arnett, who has numerically investigated violent convections in the outer layers of pre-collapse stars. Moreover, Couch and Ott demonstrated in their numerical simulations that such perturbations may lead to a successful supernova even for a progenitor that fails to explode without fluctuations. This study investigated the linear growth of perturbations during the infall onto a stalled shock wave. The linearized equations are solved as an initial and boundary value problem with the use of a Laplace transform. The background is a Bondi accretion flow whose parameters are chosen to mimic the 15 M {sub ☉} progenitor model by Woosley and Heger, which is supposed to be a typical progenitor of CCSNe. We found that the perturbations that are given at a large radius grow as they flow down to the shock radius; the density perturbations can be amplified by a factor of 30, for example. We analytically show that the growth rate is proportional to l, the index of the spherical harmonics. We also found that the perturbations oscillate in time with frequencies that are similar to those of the standing accretion shock instability. This may have an implication for shock revival in CCSNe, which will be investigated in our forthcoming paper in more detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Electromagnetic soliton-particle with spin and magnetic moment

Electromagnetic soliton-particle with both quasi-static and quick-oscillating wave parts is considered. Its mass, spin, charge, and magnetic moment appear naturally when the interaction with distant solitons is considered. The substantiation of Dirac equation for the wave part of the interacting soliton-particle is given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Multi-level coupled cluster theory

We present a general formalism where different levels of coupled cluster theory can be applied to different parts of the molecular system. The system is partitioned into subsystems by Cholesky decomposition of the one-electron Hartree-Fock density matrix. In this way the system can be divided across chemical bonds without discontinuities arising. The coupled cluster wave function is defined in terms of cluster operators for each part and these are determined from a set of coupled equations. The total wave function fulfills the Pauli-principle across all borders and levels of electron correlation. We develop the associated response theory for this multi-level coupled cluster theory and present proof of principle applications. The formalism is an essential tool in order to obtain size-intensive complexity in the calculation of local molecular properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Solution of the Burnett equations for hypersonic flows near the continuum limit

The INCA code, a three-dimensional Navier-Stokes code for analysis of hypersonic flowfields, was modified to analyze the lower reaches of the continuum transition regime, where the Navier-Stokes equations become inaccurate and Monte Carlo methods become too computationally expensive. The two-dimensional Burnett equations and the three-dimensional rotational energy transport equation were added to the code and one- and two-dimensional calculations were performed. For the structure of normal shock waves, the Burnett equations give consistently better results than Navier-Stokes equations and compare reasonably well with Monte Carlo methods. For two-dimensional flow of Nitrogen past a circular cylinder the Burnett equations predict the total drag reasonably well. Care must be taken, however, not to exceed the range of validity of the Burnett equations.

Imlay, Scott T.↗
NASA NTRS2002

Backscattering and Nonparaxiality Arrest Collapse of Damped Nonlinear Waves

The critical nonlinear Schrodinger equation (NLS) models the propagation of intense laser light in Kerr media. This equation is derived from the more comprehensive nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. It is known that if the input power of the laser beam (i.e., L(sub 2) norm of the initial solution) is sufficiently high, then the NLS model predicts that the beam will self-focus to a point (i.e.. collapse) at a finite propagation distance. Mathematically, this behavior corresponds to the formation of a singularity in the solution of the NLS. A key question which has been open for many years is whether the solution to the NLH, i.e., the 'parent' equation, may nonetheless exist and remain regular everywhere, in particular for those initial conditions (input powers) that lead to blowup in the NLS. In the current study, we address this question by introducing linear damping into both models and subsequently comparing the numerical solutions of the damped NLH (boundary-value problem) with the corresponding solutions of the damped NLS (initial-value problem). Linear damping is introduced in much the same way as done when analyzing the classical constant-coefficient Helmholtz equation using the limiting absorption principle. Numerically, we have found that it provides a very efficient tool for controlling the solutions of both the NLH and NHS. In particular, we have been able to identify initial conditions for which the NLS solution does become singular. whereas the NLH solution still remains regular everywhere. We believe that our finding of a larger domain of existence for the NLH than that for the NLS is accounted for by precisely those mechanisms, that have been neglected when deriving the NLS from the NLH, i.e., nonparaxiality and backscattering.

Fibich, G.↗
NASA NTRS1980

Unsteady transonic small disturbance theory with strong shock waves

A theory to correct the transonic small disturbance (TSD) equation to treat strong shock waves in unsteady flow is developed. The technique involves the addition of higher order terms, which are formally of negligible magnitude, to the low frequency TSD equation. These terms are then chosen such that any shock waves in the flow have strengths approximately equal to the appropriate Rankine-Hugoniot shock strength. Two correcting approaches are investigated. The first is to derive a correction for the mean steady flow and then simply use this corrected form for oscillatory flows. The second is to derive a correction for both steady and oscillatory parts of the flow. This second development is the most satisfactory and comparisons of the present results with Euler equation results are generally favorable, particularly regarding shock location, although there are some discrepancies in the pressure distribution in the leading edge region.

Kerlick, G. D.↗
DOE OSTI2026

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗
DOE OSTI2018

On the Momentum Budget of the Quasi-Biennial Oscillation in the Whole Atmosphere Community Climate Model

Abstract This study documents the contribution of equatorial waves and mesoscale gravity waves to the momentum budget of the quasi-biennial oscillation (QBO) in a 110-level version of the Whole Atmosphere Community Climate Model. The model has high vertical resolution, 500 m, above the boundary layer and through the lower and middle stratosphere, decreasing gradually to about 1.5 km near the stratopause. Parameterized mesoscale gravity waves and resolved equatorial waves contribute comparable easterly and westerly accelerations near the equator. Westerly acceleration by resolved waves is due mainly to Kelvin waves of zonal wavenumber in the range k = 1–15 and is broadly distributed about the equator. Easterly acceleration near the equator is due mainly to Rossby–gravity (RG) waves with zonal wavenumbers in the range k = 4–12. These RG waves appear to be generated in situ during both the easterly and westerly phases of the QBO, wherever the meridional curvature of the equatorial westerly jet is large enough to produce reversals of the zonal-mean barotropic vorticity gradient, suggesting that they are excited by the instability of the jet. The RG waves produce a characteristic pattern of Eliassen–Palm flux divergence that includes strong easterly acceleration close to the equator and westerly acceleration farther from the equator, suggesting that the role of the RG waves is to redistribute zonal-mean vorticity such as to neutralize the instability of the westerly jet. Insofar as unstable RG waves might be present in the real atmosphere, mixing due to these waves could have important implications for transport in the tropical stratosphere.

Garcia, Rolando R.↗
DOE OSTI2019

Propagation of surface waves along a dielectric layer in a photorefractive crystal with a diffusion mechanism for the nonlinearity formation

It is shown that TM-polarised surface waves, which exhibit a different character of their attenuation, can propagate along a dielectric layer separating photorefractive crystals with a diffusion mechanism of the nonlinearity formation. The wave profiles can be symmetric and antisymmetric with respect to the centre of the three-layer structure. For all considered wave types, dispersion equations are derived. The dependences of the propagation constant on the characteristics of the layered structure are found for long-wavelength surface waves in an explicit analytical form, and the conditions for their existence are revealed. The influence of the temperature and thickness of the dielectric layer on the propagation regimes of surface waves and their characteristics is analysed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2005

Experimental Study of Large-Amplitude Faraday Waves in Rectangular Cylinders

Experiment on single-mode Faraday waves having two, thee, and four wavelengths across a rectangular cylinder of high aspect ratio is the subject of discussion. Previous experiments recently done by Henderson & Miles (1989) and by Lei Jiang et. a1 (1996) focused on Faraday waves with one and two wavelengths across rectangular cylinders. In this experimental study the waves steepness ranges from small at threshold levels to a large amplitude which according to Penny & Price theory (1952) approaches the maximum sustainable amplitude for a standing wave. The waves characteristics for small amplitudes are evaluated against an existing well known linear theory by Benjamin & Ursell (l954) and against a weakly nonlinear theory by J. Miles (1984) which includes the effect of viscous damping. The evaluation includes the wave neutral stability and damping rate. In addition, a wave amplitude differential equation of a linear theory including viscous effect by Cerda & Tirapegui (1998) is solved numerically to yield prediction of temporal profiles of both wave damping and wave formation at the threshold. An interesting finding from this exercise is that the fluid kinematic viscosity needs to increase ten times in order to obtain good agreement between the theoretical prediction and the experimental data for both wave damping and wave starting. For large amplitude waves, the experimental data are evaluated against the theory of Penny & Price which predicts wave characteristics of any amplitude up to the point at which the wave reaches its maximum amplitude attainable for a standing wave. The theory yields two criteria to show the maximum wave steepness, the vertical acceleration at the wave crest of half the earth gravity field acceleration and the including angle at the crest of 90 degrees. Comparison with experimental data shows close agreement for the wave crest acceleration but a large discrepancy for the including angle. Additional information is included in the original extended abstract.

Iek, Chanthy↗