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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Relaxation of the loss-cone by quasi-linear diffusion of the electron-cyclotron maser instability in the solar corona

A self-consistent numeric two-dimensional code using kinetic-wave particle equations has been developed and applied to the maser dynamics of the solar corona. Time histories generated by the code reveal details of the evolution of the linear phase, the saturation and subsequent relaxation of the loss cone. Quantitative values are obtained for the saturation time, the amount of the converted free energy, the critical wave energy level for onset of maser diffusion, the velocity and pitch angle range of the resonant particles, and the relative importance of the different magnetoionic modes and harmonics. A wide range of initial conditions is covered by varying the loss-cone distribution and the ambient cold plasma parameters.

Aschwanden, Markus J.↗
DOE OSTI2013

Schematic way to find solution of the outcoupled matter wave with a source term

We propose a schematic way to obtain solution of the outcoupled atom laser beam wave function in the presence of a source term where the beam is influenced by gravity. In this case, we only focus on the external potentials inside the region of Bose-Einstein condensate that are generated by electromagnetic source and gravity. Since the evolution of the atom laser beam can be portrayed through the ordinary Schrödinger equation with a source, we are allowed to express the general solution as the superposition of the homogeneous solution and particular solution. With the given external potentials and ansatz solutions, we attain that the obtained energy depends on the parameter constituting to the ratio between the longitudinal frequency and transverse frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A comparative study of diffraction of shallow-water waves by high-level IGN and GN equations

This work is on the nonlinear diffraction analysis of shallow-water waves, impinging on submerged obstacles, by two related theories, namely the classical Green–Naghdi (GN) equations and the Irrotational Green–Naghdi (IGN) equations, both sets of equations being at high levels and derived for incompressible and inviscid flows. Recently, the high-level Green–Naghdi equations have been applied to some wave transformation problems. The high-level IGN equations have also been used in the last decade to study certain wave propagation problems. However, past works on these theories used different numerical methods to solve these nonlinear and unsteady sets of differential equations and at different levels. Moreover, different physical problems have been solved in the past. Therefore, it has not been possible to understand the differences produced by these two sets of theories and their range of applicability so far. We are thus motivated to make a direct comparison of the results produced by these theories by use of the same numerical method to solve physically the same wave diffraction problems. We focus on comparing these two theories by using similar codes; only the equations used are different but other parts of the codes, such as the wave-maker, damping zone, discretion method, matrix solver, etc., are exactly the same. This way, we eliminate many potential sources of differences that could be produced by the solution of different equations. The physical problems include the presence of various submerged obstacles that can be used for example as breakwaters or to represent the continental shelf. A numerical wave tank is created by placing a wavemaker on one end and a wave absorbing beach on the other. The nonlinear and unsteady sets of differential equations are solved by the finite-difference method. The results are compared with different equations as well as with the available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

On the fractional Harry Dym equation

We present and discuss the general case of a fractional nonlinear partial differential equation using similarity reductions and recover results associated with Harry Dym-type equations. Particular cases are discussed and some plots illustrate the results. The fractional derivative is considered in the Caputo sense.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

The limitation and applicability of Musher-Sturman equation to two dimensional lower hybrid wave collapse

The existence of localized regions of intense lower hybrid waves in the auroral ionosphere recently observed by rocket and satellite experiments can be understood by the study of a non-linear two-timescale coupling process. In this Letter, we demonstrate that the leading non-linear term in the standard Musher-Sturman equation vanishes identically in strict two-dimensions (normal to the magnetic field). Instead, the new two-dimensional equation is characterized by a much weaker non-linear term which arises from the ponderomotive force perpendicular to the magnetic field, particularly that due to the ions. The old and new equations are compared by means of time-evolution calculations of wave fields. The results exhibit a remarkable difference in the evolution of the waves as governed by the two equations. Such dissimilar outcomes motivate our investigation of the limitation of Musher-Sturman equation in quasi-two-dimensions. Only within all these limits can Musher-Sturman equation adequately describe the collapse of lower hybrid waves.

Tam, Sunny W. Y.↗
NASA NTRS1976

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

Current equation of state formulations, used for finite difference cratering flow calculations, are cast into a framework permitting comparison of peak pressures attained upon impact of a sphere, with a half-space, along the impact symmetry axis, to one dimensional impedance match solutions. On the basis of this formulation and application of thermochemical data, the regimes of melting and vaporization are examined. For the purpose of identifying material which will, upon isentropic release from the impact-induced shock state, result in a solid just brought to its melting point, i.e., incipiently melted (IM); completely melted (CM); just brought to its boiling point, i.e., incipiently vaporized (IV); and completely vaporized (CV) state, the pressures at which the critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite (GA) are examined in detail. The latter rock type is assumed to be representative of the lunar highlands. The Hugoniot pressures, for which IM, CM, IV, and CV will occur upon isentropic expansion, are calculated to be 2.2, 2.6, 4.2, and 16.8 Mbar, respectively.

Ahrens, T. J.↗
DOE OSTI2020

A Scoping Study to Determine the Location-Specific WEC Threshold Size for Wave-Powered AUV Recharging

The aim of this study is to determine the threshold wave energy converter (WEC) type and size to charge a fleet of U.S. Navy autonomous underwater vehicles (AUVs) in various geographic locations of interest. The U.S. Navy deploys AUVs in locations around the world that must be charged manually, decreasing their operational endurance and creating operational limitations. Ocean waves are a potential power source that can be converted into electricity using a WEC and stored using a battery. It would be beneficial to develop a WEC that could autonomously charge AUVs offshore. Numerous locations were analyzed to determine the minimum size of a WEC capable of providing sufficient charging power and offering a strategic advantage. By predicting the WEC efficiency (based on empirical equations) and wave resource (based on available data), electrical power generation across numerous WEC types and locations was compared in MATLAB. The generalized process developed here could be used to determine the required size and type of WECs to charge a fleet of AUVs in different locations around the world.

16 TIDAL AND WAVE POWER↗
NASA NTRS1994

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass and momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.

Lin, Yuh-Lang↗
DOE OSTI2026

Exact solution of the DeWitt-Brehme-Hobbs equation in copropagating electromagetic and gravitational waves

An accelerated charge interacts with its own electromagnetic field, a phenomenon known as electromagnetic radiation reaction. The DeWitt-Brehme-Hobbs (DWBH) equation describes the motion of a charged mass in the presence of combined electromagnetic and gravitational fields, taking into account electromagnetic radiation-reaction effects. Here, we find the first exact analytical solution of the DWBH equation in the case of a charged mass in the presence of copropagating and otherwise arbitrary electromagnetic and gravitational plane waves. As a consequence of the Penrose limit, the scenario considered here can be seen as a local limit around ultrarelativistic trajectories in a general curved spacetime. Finally, the paradigmatic example of an electromagnetic wave in the presence of a constant-amplitude gravitational wave is worked out explicitly and it is shown how the presence of the gravitational wave can qualitatively change electromagnetic radiation-reaction effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Accurate adiabatic correction in the hydrogen molecule

A new formalism for the accurate treatment of adiabatic effects in the hydrogen molecule is presented, in which the electronic wave function is expanded in the James-Coolidge basis functions. Systematic increase in the size of the basis set permits estimation of the accuracy. Numerical results for the adiabatic correction to the Born-Oppenheimer interaction energy reveal a relative precision of 10{sup −12} at an arbitrary internuclear distance. Such calculations have been performed for 88 internuclear distances in the range of 0 < R ⩽ 12 bohrs to construct the adiabatic correction potential and to solve the nuclear Schrödinger equation. Finally, the adiabatic correction to the dissociation energies of all rovibrational levels in H{sub 2}, HD, HT, D{sub 2}, DT, and T{sub 2} has been determined. For the ground state of H{sub 2} the estimated precision is 3 × 10{sup −7} cm{sup −1}, which is almost three orders of magnitude higher than that of the best previous result. The achieved accuracy removes the adiabatic contribution from the overall error budget of the present day theoretical predictions for the rovibrational levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

On the Rayleigh–Kuo criterion for the tertiary instability of zonal flows

This study reports the stability conditions for intense zonal flows (ZFs) and the growth rate γ TI of the corresponding “tertiary” instability (TI) within the generalized Hasegawa–Mima plasma model. The analytical calculation extends and revises Kuo's analysis of the mathematically similar barotropic vorticity equation for incompressible neutral fluids on a rotating sphere [H.-L. Kuo, J. Meteor. 6, 105 (1949)]; then, the results are applied to the plasma case. An error in Kuo's original result is pointed out. An explicit analytical formula for γ TI is derived and compared with numerical calculations. It is shown that, within the generalized Hasegawa–Mima model, a sinusoidal ZF is TI-unstable if and only if it satisfies the Rayleigh–Kuo criterion (known from geophysics) and that the ZF wave number exceeds the inverse ion sound radius. For non-sinusoidal ZFs, the results are qualitatively similar. As a corollary, there is no TI in the geometrical-optics limit, i.e., when the perturbation wavelength is small compared to the ZF scale. Finally, this also means that the traditional wave kinetic equation, which is derived under the geometrical-optics assumption, cannot adequately describe the ZF stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Non-Hermitian wave packet approximation of Bloch optical equations

We introduce a non-Hermitian approximation of Bloch optical equations. This approximation provides a complete description of the excitation, relaxation, and decoherence dynamics of ensembles of coupled quantum systems in weak laser fields, taking into account collective effects and dephasing. In the proposed method, one propagates the wave function of the system instead of a complete density matrix. Relaxation and dephasing are taken into account via automatically adjusted time-dependent gain and decay rates. As an application, we compute the numerical wave packet solution of a time-dependent non-Hermitian Schroedinger equation describing the interaction of electromagnetic radiation with a quantum nano-structure, and compare the calculated transmission, reflection, and absorption spectra with those obtained from the numerical solution of the Liouville-von Neumann equation. It is shown that the proposed wave packet scheme is significantly faster than the propagation of the full density matrix while maintaining small error. We provide the key ingredients for easy-to-use implementation of the proposed scheme and identify the limits and error scaling of this approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Amplitude various angles (AVA) phenomena in thin layer reservoir: Case study of various reservoirs

Amplitude various offset is widely used in petroleum exploration as well as in petroleum development field. Generally, phenomenon of amplitude in various angles assumes reservoir’s layer is quite thick. It also means that the wave is assumed as a very high frequency. But, in natural condition, the seismic wave is band limited and has quite low frequency. Therefore, topic about amplitude various angles in thin layer reservoir as well as low frequency assumption is important to be considered. Thin layer reservoir means the thickness of reservoir is about or less than quarter of wavelength. In this paper, I studied about the reflection phenomena in elastic wave which considering interference from thin layer reservoir and transmission wave. I applied Zoeppritz equation for modeling reflected wave of top reservoir, reflected wave of bottom reservoir, and also transmission elastic wave of reservoir. Results show that the phenomena of AVA in thin layer reservoir are frequency dependent. Thin layer reservoir causes interference between reflected wave of top reservoir and reflected wave of bottom reservoir. These phenomena are frequently neglected, however, in real practices. Even though, the impact of inattention in interference phenomena caused by thin layer in AVA may cause inaccurate reservoir characterization. The relation between classes of AVA reservoir and reservoir’s character are different when effect of ones in thin reservoir and ones in thick reservoir are compared. In this paper, I present some AVA phenomena including its cross plot in various thin reservoir types based on some rock physics data of Indonesia.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

An examination of a group-velocity criterion for the breakdown of an idealized vortex flow

The phenomenon of vortex breakdown is believed to be associated with a finite amplitude wave that has become trapped at the critical or breakdown location. The conditions at which the propagating waves become trapped at a certain axial location were examined by use of a group-velocity criterion implied by Landahl's general theory of wave trapping. An ideal vortex having constant vorticity and uniform axial velocity at the inlet of a slowly diverging duct was studied. The linear wave propagation analysis is applied to the base flow at several axial stations for several values of the ratio of swirl velocity to axial velocity at the inlet of the divergent duct, assuming a locally parallel flow. The dipsersion relations and hence the group velocities of both the symmetric (n = 0) and asymmetric modes (n = + or - 1) were investigated. The existence of a critical state in the flow (at which the group velocity vanishes), and its relationship to the stagnation point on the axis of the duct and to the occurrence of an irregular singularity in the equations governing wave propagation in the flow field are discussed.

Tsai, C. Y.↗
DOE OSTI2017

Coulomb Problem for Z > Z{sub cr} in Doped Graphene

The dynamics of charge carriers in doped graphene, i.e., graphene with a gap in the energy spectrum depending on the substrate, in the presence of a Coulomb impurity with charge Z is considered within the effective two-dimensional Dirac equation. The wave functions of carriers with conserved angular momentum J = M + 1/2 are determined for a Coulomb potential modified at small distances. This case, just as any two-dimensional physical system, admits both integer and half-integer quantization of the orbital angular momentum in plane, M = 0, ±1, ±2, …. For J = 0, ±1/2, ±1, critical values of the effective charge Z{sub cr}(J, n) are calculated for which a level with angular momentum J and radial quantum numbers n = 0 and n = 1 reaches the upper boundary of the valence band. For Z < Z{sub cr} (J, n = 0), the energy of a level is presented as a function of charge Z for the lowest values of orbital angular momentum M, the level with J = 0 being the first to descend to the band edge. For Z>Z{sub cr} (J, n = 0), scattering phases are calculated as a function of hole energy for several values of supercriticality, as well as the positions ε{sub 0} and widths γ of quasistationary states as a function of supercriticality. The values of ε{sub 0}* and width γ* are pointed out for which quasidiscrete levels may show up as Breit–Wigner resonances in the scattering of holes by a supercritical impurity. Since the phases are real, the partial scattering matrix is unitary, so that the radial Dirac equation is consistent even for Z > Z{sub cr}. In this single-particle approximation, there is no spontaneous creation of electron–hole pairs, and the impurity charge cannot be screened by this mechanism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗
NASA NTRS1985

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗