Search NASA⌕ Search

SEARCH · Search NASA

Results for “Helmholtz 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗
NASA NTRS1983

Difference equation technique for determining the evolution of radially symmetric ultrasonic fields

A simplified difference equation technique has been developed to calculate the three-dimensional evolution of radially symmetric, complex ultrasonic fields. The parabolic approximation to the Helmholtz equation is obtained by assuming that localized changes in the magnitude of the field with respect to distance along the axis of propagation are small compared to those with respect to radial displacement. The Helmholtz equation in cylindrical coordinates may then be approximated, after separating variables, using truncated series expansions. The resulting difference equation determines complex field values at discrete locations in a plane using corresponding adjacent locations in the preceding plane. Stability criteria which determine acceptable values for various parameters are discussed. The technique has been applied to the propagation of ultrasonic fields generated by two-dimensional focusing and non-focusing, uniform and Gaussian profile, piezoelectric transducers in both profiled and non-profiled media.

Dockery, G. D.↗
NASA NTRS1977

Non-linear duct acoustics and its application to fan noise

Quite often the scalar Helmholtz equation is assumed to be the fundamental equation of duct acoustics and various relations are derived from it. An investigation is, therefore, conducted regarding the underlying assumptions leading to the Helmholtz equation. It is found that serious errors, even well below 186 dB, are possible if the assumption of linearity were to be made. In duct acoustics these errors are highest near the cutoff points and in the case of ducts with transonic flow. It is pointed out that serious errors could arise if a nonlinear problem is analyzed by linear methods. At times, an entire phenomenon remains unpredictable by linear methods. Attention is given to the characteristics of nonlinear systems, the phenomenon of multiple pure tones in the case of noise from the aircraft engine fans, the limitations of the conventional acoustics equation, duct shape and nonlinearity, and the numerical simulation of a fan with 8 blades.

Vaidya, P. G.↗
NASA NTRS1983

Low frequency acoustic and electromagnetic scattering

This paper deals with two classes of problems arising from acoustics and electromagnetics scattering in the low frequency stations. The first class of problem is solving Helmholtz equation with Dirichlet boundary conditions on an arbitrary two dimensional body while the second one is an interior-exterior interface problem with Helmholtz equation in the exterior. Low frequency analysis show that there are two intermediate problems which solve the above problems accurate to 0(k(2) log k) where k is the frequency. These solutions greatly differ from the zero frequency approximations. For the Dirichlet problem numerical examples are shown to verify the theoretical estimates.

Hariharan, S. I.↗
NASA NTRS1986

Low frequency acoustic and electromagnetic scattering

This paper deals with two classes of problems arising from acoustics and electromagnetics scattering in the low frequency stations. The first class of problem is solving Helmholtz equation with Dirichlet boundary conditions on an arbitrary two dimensional body while the second one is an interior-exterior interface problem with Helmholtz equation in the exterior. Low frequency analysis show that there are two intermediate problems which solve the above problems accurate to 0(k/2/ log k) where k is the frequency. These solutions greatly differ from the zero frequency approximations. For the Dirichlet problem numerical examples are shown to verify the theoretical estimates.

Hariharan, S. I.↗
NASA NTRS2016

Effects of Liner Length and Attenuation on NASA Langley Impedance Eduction

This study explores the effects of liner length and attenuation on the CHE (convected Helmholtz equation) impedance eduction method, in which the surface impedance of an acoustic liner is inferred through an iterative process based on repeated solutions to the convected Helmholtz equation. Wire mesh-over-honeycomb and perforate-over-honeycomb acoustic liners are tested in the NASA Langley Grazing Flow Impedance Tube, and the resultant data are processed using two impedance eduction methods. The first is the CHE method, and the second is a direct method (labeled the KT method) that uses the Kumaresan and Tufts algorithm to compute the impedance directly. The CHE method has been extensively used for acoustic liner evaluation, but experiences anomalous behavior under some test conditions. It is postulated that the anomalies are related to the liner length and/or attenuation. Since the KT method only employs data measured over the length of the liner, it is expected to be unaffected by liner length. A comparison of results achieved with the two impedance eduction methods is used to explore the interactive effects of liner length and attenuation on the CHE impedance eduction method.

Jones, M. G.↗
NASA NTRS1989

Conjugate gradient type methods for linear systems with complex symmetric coefficient matrices

We consider conjugate gradient type methods for the solution of large sparse linear system Ax equals b with complex symmetric coefficient matrices A equals A(T). Such linear systems arise in important applications, such as the numerical solution of the complex Helmholtz equation. Furthermore, most complex non-Hermitian linear systems which occur in practice are actually complex symmetric. We investigate conjugate gradient type iterations which are based on a variant of the nonsymmetric Lanczos algorithm for complex symmetric matrices. We propose a new approach with iterates defined by a quasi-minimal residual property. The resulting algorithm presents several advantages over the standard biconjugate gradient method. We also include some remarks on the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

Freund, Roland↗
NASA NTRS1988

A noise source identification technique using an inverse Helmholtz integral equation method

A technique is developed which utilizes numerical models and field pressure information to characterize acoustic fields and identify acoustic sources. The numerical models are based on boundary element numerical procedures. Either pressure, velocity, or passive boundary conditions, in the form of impedance boundary conditions, may be imposed on the numerical model. Alternatively, if no boundary information is known, a boundary condition can be left unspecified. Field pressure data may be specified to overdetermine the numerical problem. The problem is solved numerically for the complete sound field from which the acoustic sources may be determined. The model can then be used to idenfify acoustic intensity paths in the field. The solution can be modified and the model used to evaluate design alternatives. In this investigation the method is tested analytically and verified. In addition, the sensitivity of the method to random and bias error in the input data is demonstrated.

Gardner, B. K.↗
NASA NTRS1994

A fast numerical solution of scattering by a cylinder: Spectral method for the boundary integral equations

It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.

Hu, Fang Q.↗
NASA NTRS1974

Analysis of sound propagation in ducts using the wave envelope concept

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.

Baumeister, K. J.↗
NASA NTRS1975

Generalized wave envelope analysis of sound propagation in ducts with stepped noise source profiles and variable axial impedance

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution tends not to oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude. Example solutions indicate that stepped noise source profiles have much higher attenuation than plane waves in a uniform impedance liner. Also, multiple stepped impedance liners are shown to have higher attenuation than uniform ducts if the impedances are chosen properly. For optimum noise reduction with axial variations in impedance, the numerical analysis indicates that for a plane wave input the resistance should be near zero at the entrance of a suppressor duct, while the reactance should be near the optimum value associated with the least-attenuated mode in a uniform duct.

Baumeister, K. J.↗
NASA NTRS1999

The Prediction of Noise Scattered by a Wing/Ducted Fan Configuration

In this proof of concept research, a computational method was developed for predicting the sound field created by the scattering of ducted fan engine noise by a blended wing-body (BWB). It was assumed that all acoustic processes were linear and time harmonic with excitation frequency co. Inflow effects were neglected and no penetration boundary conditions were applied to the engine nacelle and BWB surfaces. A scattering approach was adopted in which the total acoustic field is written as the sum of known incident (from the engine duct) and unknown scattered parts. We further assume that the incident field is independent of the scattered field. Application of the above conditions to the equations of linearized acoustics yields the Helmholtz equation (reduced wave equation) for the scattered pressure with Neumann boundary conditions.

Tweed, John↗
NASA NTRS1975

Gravity shear waves atop the cirrus layer of intense convective storms

Recent visual satellite photographs of certain intense convective storms have revealed concentric wave patterns. A model for the generation and growth of these waves is proposed. The proposed initial generating mechanism is similar to the effect noticed when a pebble is dropped into a calm pond. The penetration of the tropopause by overshooting convection is analogous to the pebble's penetration of the water's surface. The model for wave growth involves instability due to the wind shear resulting from the cirrus outflow. This model is based on an equation for the waves' phase speed which is similar to the Helmholtz equation. It, however, does not assume an incompressible atmosphere, but rather assumes density is a logarithmic function of height. Finally, the model is evaluated on the two mid-latitude and three tropical cases. The data indicate that shearing instability may be a significant factor in the appearance of these waves.

Stobie, J. G.↗
NASA NTRS1996

A Computer Program for the Computation of Running Gear Temperatures Using Green's Function

A new technique has been developed to study two dimensional heat transfer problems in gears. This technique consists of transforming the heat equation into a line integral equation with the use of Green's theorem. The equation is then expressed in terms of eigenfunctions that satisfy the Helmholtz equation, and their corresponding eigenvalues for an arbitrarily shaped region of interest. The eigenfunction are obtalned by solving an intergral equation. Once the eigenfunctions are found, the temperature is expanded in terms of the eigenfunctions with unknown time dependent coefficients that can be solved by using Runge Kutta methods. The time integration is extremely efficient. Therefore, any changes in the time dependent coefficients or source terms in the boundary conditions do not impose a great computational burden on the user. The method is demonstrated by applying it to a sample gear tooth. Temperature histories at representative surface locatons are given.

GREENS FUNCTIONS↗
NASA NTRS1985

Interaction of finite-amplitude sound with air-filled porous materials

The propagation of high intensity sound waves through an air-filled porus material was studied. The material is assumed: (1) to be rigid, incompressible, and homogeneous, and (2) to be adequately described by two properties: resistivity r and porosity. The resulting wave equation is still nonlinear, however, because of the u sgn(u) term in the resistivity. The equation is solved in the frequency domain as an infinite set of coupled inhomogeneous Helmholtz equations, one for each harmonic. An approximate but analytical solution leads to predictions of excess attenuation, saturation, and phase speed reduction for the fundamental component. A more general numerical solution is used to calculate the propagation curves for the higher harmonics. The u sgn(u) nonlinearity produces a cubic distortion pattern; when the input signal is a pure tone, only odd harmonic distortion products are generated.

Nelson, D. A.↗
NASA NTRS2001

Computation of Nonlinear Backscattering Using a High-Order Numerical Method

The nonlinear Schrodinger equation (NLS) is the standard model for propagation of intense laser beams in Kerr media. The NLS is derived from the nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. In this study we use a fourth-order finite-difference method supplemented by special two-way artificial boundary conditions (ABCs) to solve the NLH as a boundary value problem. Our numerical methodology allows for a direct comparison of the NLH and NLS models and for an accurate quantitative assessment of the backscattered signal.

Fibich, G.↗