SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Vibrations of three-dimensional pipe systems with acoustic coupling

A general algorithm is developed to calculate the beam-type dynamic response of three dimensional multiplane finite length pipe systems, consisting of elbow and straight ducts with continuous interfaces. Emphasis is on secondary acoustic wave effects giving rise to coupling mechanisms; and the simulation accounts for one-dimensional elastoacoustic coupling from a plane acoustic wave and secondary loads resulting from wave asymmetries. The transfer matrix approach is adopted in modeling the elastodynamics of each duct, with allowance for distribution loads. Secondary loads from plane wave distortion are considered with a solution of the Helmholtz equation in an equivalent rigid waveguide, and effects of path imperfection are introduced as a perturbation from the hypothetical perfectly straight pipe. Computations indicate that the one-dimensional acoustic assumption is valid for frequencies below one-half the first cut-off frequency, and the three-dimensional acoustic effects produce an increase in response levels near and above cut-off.

El-Raheb, M.↗

NASA NTRS1981

Inverse scattering for an exterior Dirichlet program

Scattering due to a metallic cylinder which is in the field of a wire carrying a periodic current is considered. The location and shape of the cylinder is obtained with a far field measurement in between the wire and the cylinder. The same analysis is applicable in acoustics in the situation that the cylinder is a soft wall body and the wire is a line source. The associated direct problem in this situation is an exterior Dirichlet problem for the Helmholtz equation in two dimensions. An improved low frequency estimate for the solution of this problem using integral equation methods is presented. The far field measurements are related to the solutions of boundary integral equations in the low frequency situation. These solutions are expressed in terms of mapping function which maps the exterior of the unknown curve onto the exterior of a unit disk. The coefficients of the Laurent expansion of the conformal transformations are related to the far field coefficients. The first far field coefficient leads to the calculation of the distance between the source and the cylinder.

Hariharan, S. I.↗

NASA NTRS1982

Prediction of sound radiation from different practical jet engine inlets

Computer codes which were based upon a special integral representation of the external solutions of the Helmholtz equation were upgraded so that they would yield accurate results for the acoustic radiation patterns in the field surrounding an axisymmetric body for nondimensional wave numbers, based on duct radius, of up to twenty. The accuracy of these computer programs was checked by the use of the point source method for the generation of exact solutions and then by comparison with the results of other experimental and theoretical studies. These computer codes were used in a parametric study of the dependence of the radiated sound field on input modal distribution, wave number, and inlet lip shape. The results of this study show that: (1) as the wave number (i.e., cut off ratio) is increased for a given input modal distribution that the acoustic radiation peak moves towards the inlet centerline and becomes more compact (i.e., narrows); (2) as the input mode number is increased for a given cut off ratio the acoustic radiation peak moves away from the inlet centerline and becomes more compact.

Meyer, W. L.↗

NASA NTRS1982

Inverse scattering for an exterior Dirichlet problem

Scattering caused by a metallic cylinder in the field of a wire carrying a periodic current is studied, with a view to determining the location and shape of the cylinder in light of far field measurements between the cylinder and the wire. The associated direct problem is the exterior Dirichlet problem for the Helmholtz equation in two dimensions, and an improved low frequency estimate for its solution by integral equation methods is shown by inverse scattering calculations to be accurate to this estimate. The far field measurements are related to low frequency boundary integral equations whose solutions may be expressed in terms of a mapping function for the exterior of the unknown curve onto the exterior of a unit disk. The conformal transformation's Laurent expansion coefficients can be related to those of the far field, the first of which leads to the calculation of the distance between the source and the cylinder, while the other coefficients are determined by placing the source in a different location.

Hariharan, S. I.↗

NASA NTRS1983

A theoretical investigation of the sound radiation fields associated with a Bellmouth inlet

Analytical results are obtained by numerical integration of a cylindrically symmetric integral representation of the external solutions of the Helmholtz equation. The accuracy of this method is checked by comparisons of computed results with 'exact' solutions generated by the point source method. In all cases, the average error for the amplitude and phase of the points calculated in the field is found to be less than ten percent. Theoretical studies which model experiments run for the NASA Langely Bellmouth inlet configuration are presented and comparisons are made with the experimental results. In all cases very good agreement is obtained between the experimental and theoretically calculated values.

Meyer, W. L.↗

NASA NTRS1983

Spinning mode acoustic radiation from the flight inlet

A mathematical model was developed for spinning mode acoustic radiation from a thick wall duct without flow. This model is based on a series of experiments (with and without flow). A nearly pure azimuthal spinning mode was isolated and then reflection coefficients and far field pressure (amplitude and phase) were measured. In our model the governing boundary value problem for the Helmholtz equation is first converted into an integral equation for the unknown acoustic pressure over a disk, S1, near the mouth of the duct and over the exterior surface, S2, of the duct. Assuming a pure azimuthal mode excitation, the azimuthal dependence is integrated out which yields an integral equation over the generator C1 of S1 and the generator C2 of S2. The sound pressure on C1 was approximated by a truncated modal expansion of the interior acoustic pressure. Piecewise linear spline approximation on C2 was used.

Moss, W. F.↗

NASA NTRS1984

Conjugate gradient coupled with multigrid for an indefinite problem

An iterative algorithm for the Helmholtz equation is presented. This scheme was based on the preconditioned conjugate gradient method for the normal equations. The preconditioning is one cycle of a multigrid method for the discrete Laplacian. The smoothing algorithm is red-black Gauss-Seidel and is constructed so it is a symmetric operator. The total number of iterations needed by the algorithm is independent of h. By varying the number of grids, the number of iterations depends only weakly on k when k(3)h(2) is constant. Comparisons with a SSOR preconditioner are presented.

Gozani, J.↗

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 NTRS1984

A parametric study of the effect of inlet lip shape upon the radiated sound field

Far field sound radiation predictions for four different axisymmetric inlet lips excited by different tangential acoustic modes at several cut-off ratios are presented. These results were obtained by numerical integration of a special cylindrically symmetric integral representation of the external solutions of the Helmholtz equation which yields unique solutions at all wave numbers. The paper presents plots which detail the dependence of the relative SPL (dB) in the field upon the engine inlet lip shape, the modal input, and the cut-off ratio. Examination of these data indicate that: (1) as the inlet lip becomes larger the predominant acoustic radiation peak in the field becomes narrower and moves towards the centerline of the inlet; (2) as the order of the tangential acoustic mode of the driver increases the radiated sound peak again becomes narrower but moves away from the inlet centerline; and (3) as the cut-off ratio is increased for a specific tangential acoustic mode the predominant radiation peak becomes narrower and moves towards the centerline of the inlet.

Meyer, W. L.↗

NASA NTRS1985

On accuracy conditions for the numerical computation of waves

Bayliss, A.↗

NASA NTRS1988

Light scattering of large rough particles application to cometary grains

While the electromagnetic field scattered by a spherical particle is classically obtained by the Helmholtz equation, the general case of an arbitrary particle may be investigated in the general framework of the interaction of a wave with a scattering potential. A wave function then satisfies the Schroedinger equation. The general solution of the Schroedinger equation is given. The main disadvantage of this approach are its restriction to large particles and its scalar nature preventing the calculation of the polarization. However, Perrin and Lamy have shown how to avoid the second limitation and retrieve a vectorial description. They proved that in the case of large spheres when the ad hoc assumptions are satisfied, the expression of the scattering amplitude may be approximated by an expansion series in partial waves, i.e., on a discrete basis. The analogy may be generalized, and the ratio of the two components for a rough particle obtained by taking the ratio of the reflectivities for the two directions of polarization. These reflectivities involve the simple and double reflections calculated following the method developed by Wolff for rough surfaces. The theory is further detailed.

Lamy, Philippe L.↗

NASA NTRS1988

On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices

Conjugate gradient type methods are considered for the solution of large linear systems Ax = b with complex coefficient matrices of the type A = T + i(sigma)I where T is Hermitian and sigma, a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidian error minimization, respectively, are investigated. In particular, numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices are proposed. Error bounds for all three methods are derived. It is shown how the special shift structure of A can be preserved by using polynomial preconditioning. Results on the optimal choice of the polynomial preconditioner are given. Also, some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation are reported.

Freund, Roland↗

NASA NTRS1989

Adaptation of boundary elements to the dynamics of elastic solids

The method of boundary elements is adapted to the dynamics of elastic axisymmetric solids. Three scalar potentials are used, each satisfying the Helmholtz equation and each represented by its own source distribution. The sources are distributed over a surface enclosing the solid boundary. Stress or displacement conditions on the boundary are realized by third-order derivatives of the potentials. Discrete ring elements on the source enclosure combine their influences at discrete ring elements on the solid boundary.

El-Raheb, Michael↗

NASA NTRS1992

The Pekeris problem, and the interior of the pulsating sphere, first success of exact ray theory, and what rays?

This paper presents what is probably the first application of exact ray theory to provide nontrivial exact solutions of the Helmholtz equation. First, the Pekeris problem is considered. This is followed by consideration of the problem of the field maintained within a uniformly pulsating sphere. In the Pekeris problem, an iterative solution provides the exact solution on the second iteration. In the second problem, it is shown that the exact ray theory eliminates the focus of the geometric theory.

Mcaninch, Gerry L.↗

NASA NTRS1998

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗

NASA NTRS1993

Mode Profiles in Waveguide-Coupled Resonators

Surface acoustic wave (SAW) waveguide-coupled resonators are of considerable interest for narrow-band filter applications, though to date there has been very little published on the acoustic details of their operation. As in any resonator, one must fully understand its mode structure and herein we study the SAW mode profiles in these devices. Transverse mode profiles in the resonant cavity of the device were measured at various frequencies of interest using a knife-edge laser probe. In addition we predict the mode profiles for the device structure by two independent methods. One is a stack-matrix approach adapted from integrated optics and the other is a conventional analytical eigenmode analysis of the Helmholtz equation. Both modeling techniques are in good agreement with the measured results.

Hunt, William D.↗

NASA NTRS2001

Radiation Losses Due to Tapering of a Double-Core Optical Waveguide

The theoretical model we designed parameterizes the power losses as a function of .the profile shape for a tapered, single mode, optical dielectric coupler. The focus of this project is to produce a working model that determines the power losses experienced by the fibers when light crosses a taper region. This phenomenon can be examined using coupled mode theory. The optical directional coupler consists of a parallel, dual-channel, waveguide with minimal spacing between the channels to permit energy exchange. Thus, power transfer is essentially a function of the taper profile. To find the fields in the fibers, the approach used was that of solving the Helmholtz equation in cylindrical coordinates involving Bessel and modified Bessel functions depending on the location.

Lyons, Donald R.↗

NASA NTRS2001

Turbofan Duct Propagation Model

The CDUCT code utilizes a parabolic approximation to the convected Helmholtz equation in order to efficiently model acoustic propagation in acoustically treated, complex shaped ducts. The parabolic approximation solves one-way wave propagation with a marching method which neglects backwards reflected waves. The derivation of the parabolic approximation is presented. Several code validation cases are given. An acoustic lining design process for an example aft fan duct is discussed. It is noted that the method can efficiently model realistic three-dimension effects, acoustic lining, and flow within the computational capabilities of a typical computer workstation.

Lan, Justin H.↗