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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Supersonic propeller noise in a uniform flow

The reported investigation is concerned with an extension of a theory presented by Hawkings and Lowson (1974). The extension considers also forward flight. Attention is given to the retarded potential solution of the convected wave equation, the Fourier components of acoustic pressure, the zones of relative silence and Doppler amplification, and the asymptotic evaluation of a noise field for a large wave number. The discussed analysis provides the Fourier coefficients of the acoustic signature of a supersonic rotor in a uniform flow. The results are represented by an integral over the blade planform which must be evaluated numerically. The fast Fourier transform can be employed for the Fourier inversion.

Jou, W.-H.↗

NASA NTRS1981

Two-dimensional analytical model of twin jet shielding

Estimation of the impact of jet noise requires analytical modeling not only of the noise source mechanism, but also of the propagation of noise from the source to receiver. One factor which influences the jet noise path is shielding of one jet by another in a twin jet configuration. An analytical model is developed to investigate the shielding phenomenon. The two-dimensional wave equation is solved for a stationary line source impinging upon a cylinder of heated flow. The solution estimates the diffraction and scattering of the incident sound wave by the shielding jet in a plane normal to the jet axis. The frequency dependence of the normalized sound pressure estimated by the model is found to agree in form with empirical data. The azimuthal redistribution of the incident sound shows that, as the frequency of maximum shielding is approached, the scattered sound becomes more highly concentrated into lobes adjacent to the shielding zone.

Gerhold, C. H.↗

DOE OSTI2023

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Heating of the solar corona by the resonant absorption of Alfven waves

An improved method for calculating the resonance absorption heating rate is discussed and the results are compared with observations in the solar corona. To accomplish this, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients. For parameters representative of the solar corona, it is found that a two-scale description of the wave motion is appropriate. The large-scale motion, which can be approximated as nearly ideal, has a scale which is on the order of the width of the loop. The small-scale wave, however, has a transverse scale much smaller than the width of the loop, with a width of about 0.3-250 km, and is highly dissipative. These two wave motions are coupled in a narrow resonance region in the loop where the global wave frequency equals the local Alfven wave frequency. Formally, this coupling comes about from using the method of matched asymptotic expansions to match the inner and outer (small and large scale) solutions. The resultant heating rate can be calculated from either of these solutions. A formula derived using the outer (ideal) solution is presented, and shown to be consistent with observations of heating and line broadening in the solar corona.

Davila, Joseph M.↗

NASA NTRS1989

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗

NASA NTRS1976

Acoustic radiation from sources convected at transonic speeds

Each supersonic jet contains a long segment of transonic flow, and it is known through acoustic experiments that the strongest noise production region is embedded in this segment. The linearized convected wave equation breaks down in the transonic regime and a nonlinear equation must be used. Amplitude, directivity, and the transonic similarity rule for noise radiation can be established through proper matching of the displacement amplitude between the source and the wave field. This paper is a first step towards establishing a general theory of transonic aeroacoustics.

Pao, S. P.↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

NASA NTRS1980

Experimental determination of the admittances of aluminized propellants by the impedance tube method

The adaptation of the impedance tube concept for the determination of the pressure coupled admittances and response functions of burning solid propellants is discussed. The results obtained in experiments with UTP-3001 and UTP-19360 aluminized propellants are presented. It is shown that the admittance Y remains constant during the quasi-steady burn period of a test, indicating constant driving of the gas phase disturbance by the burning solid propellant. The measured real part of the admittance is positive, indicating that the burning aluminized propellant is driving the gas phase oscillations. In addition, the measured high gas phase damping, provided by the aluminum oxide particles in the gas phase, suggests that the latter can significantly increase the damping in unstable solid rockets over the investigated frequency range. Finally, it is shown that the wave structure obtained by numerically solving the impedance tube wave equations which utilize the determined propellant admittance as an initial condition and the determined value of G to describe the gas phase losses is in excellent agreement with the measured wave structure.

Baum, J. D.↗

NASA NTRS1986

MHD bending waves in a current sheet

Transverse MHD bending waves are considered in an isothermal and compressible two-dimensional current sheet of finite thickness in which the magnetic field changes direction and strength. The general form of the wave equation is obtained. It is shown that rotation of the magnetic field across the current sheet prevents the existence of singular points so that continuous spectrum solutions and the concomitant wave decay disappear. Instead, normal modes exist and closed integral solution for arbitrary current sheet structure are found. The results are discussed in terms of small-scale waves on the heliospheric current sheet.

Musielak, Z. E.↗

NASA NTRS1990

Development of a linearized unsteady aerodynamic analysis for cascade gust response predictions

A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Verdon, Joseph M.↗

DOE OSTI2018

Mode-conversion of the extraordinary wave at the upper hybrid resonance in the presence of small-amplitude density fluctuations

We report in spherical tokamaks, the electron plasma frequency is greater than the electron cyclotron (EC) frequency. Electromagnetic waves in the EC range of frequencies are unsuitable for directly heating such plasmas due to their reduced accessibility. However, mode-conversion of the extraordinary wave to the electron Bernstein wave (X–B mode-conversion) at the upper hybrid resonance makes it possible to efficiently couple externally-launched electromagnetic wave energy into an overdense plasma core. Traditional mode-conversion models describe an X-mode wave propagating in a potential containing two cutoffs that bracket a single wave resonance. Often, however, the mode-conversion region is in the edge, where turbulent fluctuations and blobs can generate abrupt cutoffs and scattering of the incident X-mode wave. We present a new framework for studying the X–B mode-conversion which makes the inclusion of these fluctuations analytically tractable. In the new approach, the high-field cutoff is modeled as an infinite barrier, which manifests as a boundary condition applied to a wave equation involving only one cutoff adjacent to the resonance on the low-field side. The new model reproduces the main features of the previous approach, yet is more suitable for analyzing experimental observations and extrapolating to higher dimensions. We then develop an analytical estimate for the effect of small-amplitude, quasi-monochromatic density fluctuations on the X–B mode-conversion efficiency using perturbation theory. We find that Bragg backscattering of the launched X-mode wave reduces the mode-conversion efficiency significantly when the fluctuation wavenumber is resonant with the wavenumber of the incident X-mode wave. These analytical results are corroborated by numerically integrating the mode-conversion equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Mode-conversion of the extraordinary wave at the upper hybrid resonance in the presence of small-amplitude density fluctuations

In spherical tokamaks, the electron plasma frequency is greater than the electron cyclotron frequency. Electromagnetic waves in the electron cyclotron range of frequencies are unsuitable for directly heating such plasmas due to their reduced accessibility. However, mode-conversion of the extraordinary wave to the electron Bernstein wave (X-B mode-conversion) at the upper hybrid resonance makes it possible to efficiently couple externally-launched electromagnetic wave energy into an overdense plasma core. Traditional mode-conversion models describe an X-mode wave propagating in a potential containing two cutoffs that bracket a single wave resonance. Often, however, the mode-conversion region is in the edge, where turbulent fluctuations and blobs can generate abrupt cutoffs and scattering of the incident X-mode wave. We present a new framework for studying the X-B mode-conversion which makes the inclusion of these fluctuations analytically tractable. In the new approach, the high- field cutoff is modelled as an infinite barrier, which manifests as a boundary condition applied to a wave equation involving only one cutoff adjacent to the resonance on the low-field side. The new model reproduces the main features of the previous approach, yet is more suitable for analyzing experimental observations and extrapolating to higher dimensions. We then develop an analytical estimate for the effect of small-amplitude, quasi-monochromatic density fluctuations on the X-B mode-conversion efficiency using perturbation theory. We find that Bragg backscattering of the launched X-mode wave reduces the mode-conversion efficiency significantly when the fluctuation wavenumber is resonant with the wavenumber of the incident X-mode wave. These analytical results are corroborated by numerically integrating the mode-conversion equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1981

Equations of motion for a free-electron laser with an electromagnetic pump field and an axial electrostatic field

The equations of motion for a free-electron laser (FEL) with an electromagnetic pump field and a static axial electric field are derived using a Hamiltonian formalism. Equations governing the energy transfer between the electron beam and each of the electromagnetic fields are given, and the phase shift for each of the electromagnetic fields is derived from a linearized Maxwell wave equation. The relation between the static axial electric field and the resonant phase is given. Laser gain and the fraction of the electron energy converted to photon energy are determined using a simplified resonant particle model. These results are compared to those of a more exact particle simulation code.

Hiddleston, H. R.↗

NASA NTRS1974

The influence of the computational mesh on accuracy for initial value problems with discontinuous or nonunique solutions

Discontinuous, or weak, solutions of the wave equation, the inviscid form of Burgers equation, and the time-dependent, two-dimensional Euler equations are studied. A numerical method of second-order accuracy in two forms, differential and integral, is used to calculate the weak solutions of these equations for several initial value problems, including supersonic flow past a wedge, a double symmetric wedge, and a sphere. The effect of the computational mesh on the accuracy of computed weak solutions including shock waves and expansion phenomena is studied. Modifications to the finite-difference method are presented which aid in obtaining desired solutions for initial value problems in which the solutions are nonunique.

Maccormack, R. W.↗

NASA NTRS1990

Long-time behaviour of absorbing boundary conditions

A new class of computational far-field boundary conditions for hyperbolic partial differential equations was recently introduced by the authors. These boundary conditions combine properties of absorbing conditions for transient solutions and properties of far-field conditions for steady states. This paper analyses the properties of the wave equation coupled with these new boundary conditions: well-posedness, dissipativity and convergence in time.

Engquist, B.↗

NASA NTRS1981

Analytical and experimental investigation of the propagation and attenuation of sound in extended reaction lined ducts

Results are presented of an analytical and experimental study of the attenuation and propagation of harmonically excited sound waves in an extended reaction lined cylindrical duct. The duct geometry considered consisted of an annular outer region of bulk material surrounding an inner cylinder of air. The coupled wave equations governing the motion of the sound in both the inner and annular regions were solved numerically. The numerically predicted attenuation and propagation constants were in excellent agreement with measured values using Kevlar as the liner material for plane-wave mode (0,0) excitation over the frequency from 100 to 7,000 Hz. Although the numerical model was verified using Kevlar, it can be used with any fibrous constructed bulk liner. The results of this study demonstrate that a good start has been made on the numerical modelling of the acoustic performance of extended reaction liners.

Hersh, A. S.↗

NASA NTRS1988

Far field boundary conditions for computation over long time

A new class of computational far field boundary conditions for hyperbolic partial differential equations is developed. These boundary conditions combine properties of absorbing boundary conditions for transient solutions and properties of far field boundary conditions for steady-state problems. The conditions can be used to limit the computational domain when both traveling waves and evanescent waves are present. Boundary conditions for scalar wave equations are derived and analyzed. Extensions to systems of equations are discussed and results from numerical experiments are presented.

Engquist, Bjorn↗

NASA NTRS1990

Matrix-based, finite-difference algorithms for computational acoustics

A compact numerical algorithm is introduced for simulating multidimensional acoustic waves. The algorithm is expressed in terms of a set of matrix coefficients on a three-point spatial grid that approximates the acoustic wave equation with a discretization error of O(h exp 5). The method is based on tracking a local phase variable and its implementation suggests a convenient coordinate splitting along with natural intermediate boundary conditions. Results are presented for oblique plane waves and compared with other procedures. Preliminary computations of acoustic diffraction are also considered.

Davis, Sanford↗