Engineering Papers⌕ Search

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering

When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.

Fibich, Gadi↗
NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗
NASA NTRS2012

Prediction of Broadband Shock-Associated Noise Including Propagation Effects Originating NASA

An acoustic analogy is developed based on the Euler equations for broadband shock-associated noise (BBSAN) that directly incorporates the vector Green s function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) to describe the mean flow. The vector Green s function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation rate. An adjoint vector Green s function solver is implemented to determine the vector Green s function based on a locally parallel mean flow at different streamwise locations. The newly developed acoustic analogy can be simplified to one that uses the Green s function associated with the Helmholtz equation, which is consistent with a previous formulation by the authors. A large number of predictions are generated using three different nozzles over a wide range of fully-expanded jet Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise experimental facilities. In addition, two models for the so-called fine-scale mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include propagation effects.

Miller, Steven↗
NASA NTRS1993

Mathmatical modeling for diffractive optics

We consider a 'diffractive optic' to be a biperiodic surface separating two half-spaces, each having constant constitutive parameters; within a unit cell of the periodic surface and across the transition zone between the two half-spaces, the constitutive parameters can be a continuous, complex-valued function. Mathematical models for diffractive optics have been developed, and implemented as numerical codes, both for the 'direct' problem and for the 'inverse' problem. In problems of the 'direct' class, the diffractive optic is specified, and the full set of Maxwell's equations is cast in a variational form and solved numerically by a finite element approach. This approach is well-posed in the sense that existence and uniqueness of the solution can be proved and specific convergence conditions can be derived. An example of a metallic grating at a Wood anomaly is presented as a case where other approaches are known to have convergence problems. In problems of the 'inverse' class, some information about the diffracted field (e.g., the far-field intensity) is given, and the problem is to find the periodic structure in some optimal sense. Two approaches are described: phase reconstruction in the far-field approximation; and relaxed optimal design based on the Helmholtz equation. Practical examples are discussed for each approach to the inverse problem, including array generators in the far-field case and antireflective structures for the relaxed optimal design.

Dobson, David↗
NASA NTRS1991

Krylov Subspace Methods for Complex Non-Hermitian Linear Systems

We consider Krylov subspace methods for the solution of large sparse linear systems Ax = b with complex non-Hermitian coefficient matrices. Such linear systems arise in important applications, such as inverse scattering, numerical solution of time-dependent Schrodinger equations, underwater acoustics, eddy current computations, numerical computations in quantum chromodynamics, and numerical conformal mapping. Typically, the resulting coefficient matrices A exhibit special structures, such as complex symmetry, or they are shifted Hermitian matrices. In this paper, we first describe a Krylov subspace approach with iterates defined by a quasi-minimal residual property, the QMR method, for solving general complex non-Hermitian linear systems. Then, we study special Krylov subspace methods designed for the two families of complex symmetric respectively shifted Hermitian linear systems. We also include some results concerning 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 W.↗
NASA NTRS2011

The Prediction of Broadband Shock-Associated Noise Including Propagation Effects

An acoustic analogy is developed based on the Euler equations for broadband shock- associated noise (BBSAN) that directly incorporates the vector Green's function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) as the mean flow. The vector Green's function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation. An adjoint vector Green's function solver is implemented to determine the vector Green's function based on a locally parallel mean flow at streamwise locations of the SRANS solution. However, the developed acoustic analogy could easily be based on any adjoint vector Green's function solver, such as one that makes no assumptions about the mean flow. The newly developed acoustic analogy can be simplified to one that uses the Green's function associated with the Helmholtz equation, which is consistent with the formulation of Morris and Miller (AIAAJ 2010). A large number of predictions are generated using three different nozzles over a wide range of fully expanded Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise labs. In addition, two models for the so-called 'fine-scale' mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include the propagation effects, especially in the upstream direction of the jet.

Miller, Steven↗
NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗
DOE OSTI2023

Phase diagram to demarcate supercritical, transcritical, and continuous phase regimes for binary fluid equilibrium mixing relevant to combustion applications

Here, a robust methodology to develop phase diagrams of binary fluid mixtures at fixed thermo dynamic conditions (pressure, temperature, and mole fraction) as well as of two initially separated fluids undergoing mixing near critical conditions are presented for fluids and con ditions relevant to rockets, gas turbines, and diesel engine applications. Phase equilibria of mixtures is first examined to provide insight into the continuous-phase mixing behavior (including but not limited to supercritical behavior), and to develop a broadly applicable phase-diagram for binary fluid mixtures at fixed conditions. Next, adiabatic mixing theory and reduced Helmholtz equations of state are used to predict the thermodynamic conditions required to attain continuous-phase binary fluid mixing near critical conditions. Then, a 3D surface diagram (P,T fuel ,T amb ) separating single and two-phase regions is constructed by varying the ambient pressure and the initial temperatures of the two fluids. Polynomial fits of the 3D surfaces for 10 different binary mixtures are tabulated for nitrogen-alkane and methane-oxygen blends relevant to air-breathing and propellant based engines, respectively.

42 ENGINEERING↗
DOE OSTI2020

Experimental determination of the solubility constant of kurnakovite, MgB3O3(OH)5·5H2O

In this study, I present experimental results on the equilibrium between boracite [Mg 3 B 7 O 13 C1(cr)] and kurnakovite [chemical formula, Mg 2 B 6 O 11 ·15H 2 O(cr); structural formula, MgB 3 O 3 (OH) 5 ·5H 2 O(cr)] at 22.5 ± 0.5 °C from a long-term experiment up to 1629 days, approaching equilibrium from the direction of supersaturation, Mg 3 B 7 O 13 C1(cr) + H + + 2B(OH)$^-_4$ + 18H 2 O(1) ⇌ 3MgB 3 O 3 (OH) 5 ·5H 2 O(cr) + C1 - . Based on solubility measurements, the 10-based logarithm of the equilibrium constant for the above reaction at 25 °C is determined to be 12.83 ± 0.08 (2σ). Based on the equilibrium constant for dissolution of boracite, Mg 3 B 7 O 13 C1(cr) + 15H 2 O(1) = 3Mg 2+ + 7B(OH)$^-_4$ + C1 - + 2H + at 25 °C measured previously (Xiong et al. 2018) and that for the reaction between boracite and kurnakovite determined here, the equilibrium constant for dissolution of kurnakovite, MgB 3 O 3 (OH) 5 ·5H 2 O(cr) = Mg 2+ + 3B(OH)$^-_4$ + H + + H 2 O(1) is derived as -14.11 ± 0.40 (2σ). Using the equilibrium constant for dissolution of kurnakovite obtained in this study and the experimental enthalpy of formation for kurnakovite from the literature, a set of thermodynamic properties for kurnakovite at 25 °C and 1 bar is recommended as follows: Δ H $^0_f$ = -4813.24 ± 4.92 kJ/mol, Δ G$^0_f$ = -4232.0 ± 2.3 kJ/mol, and S 0 = 414.3 ± 0.9 J/(mol·K). Among them, the Gibbs energy of formation is based on the equilibrium constant for kurnakovite determined in this study; the enthalpy of formation is from the literature (Li et al. 1997), and the standard entropy is calculated internally with the Gibbs-Helmholtz equation in this work. The thermodynamic properties of kurnakovite estimated using the group contribution method for borate minerals based on the sums of contributions from the cations, borate polyanions, and structural water to the thermodynamic properties from the literature (Li et al. 2000) are consistent, within their uncertainties, with the values listed above. Since kurnakovite usually forms in salt lakes rich in sulfate, studying the interactions of borate with sulfate is important to modeling kurnakovite in salt lakes. Finally, for this purpose, I have re-calibrated our previous model (Xiong et al. 2013) describing the interactions of borate with sulfate based on the new solubility data for borax in Na 2 SO 4 solutions presented here.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Semi-Algebraic Two Level Solver

We develop a simple semi-algebraic 2-level solver built on traditional multigrid ideas. It is designed to be easily incorporated into existing simulation software. It exhibits good convergence for many classes of challenging problems including discontinuous diffusion, convection- diffusion, and Helmholtz equations. It has built-in structure that makes it simple to generalize in several interesting directions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

A spectral domain method for remotely probing stratified media

The problem of remotely probing a stratified, lossless, dielectric medium is formulated using the spectral domain method of probing. The response of the medium to a spectrum of plane waves incident at various angles is used to invert the unknown profile. For TE polarization, the electric field satisfies a Helmholtz equation. The inverse problem is solved by means of a new representation for the wave function. The principal step in this inversion is solving a second kind Fredholm equation which is very amenable to numerical computations. Several examples are presented including some which indicate that the method can be used with experimentally obtained data. When the fields exhibit a surface wave behavior, a unique inversion can be obtained only if information about the magnetic field is also available. In this case, the inversion is accomplished by a two-step procedure which employs a formula of Jost and Kohn. Some examples are presented, and an approach which greatly shortens the computations without greatly deteriorating the results is discussed.

Schaubert, D. H.↗
NASA NTRS1978

A finite element method for transmission in non-uniform ducts without flow Comparison with the method of weighted residuals

A finite element method is derived for studying the transmission of sound in nonuniform ducts without flow. The method is based on a weighted residuals technique and eight noded isoparametric elements. The method employs two computational schemes: the Helmholtz equation found through conservation equations, and the conservation equations themselves. Solutions are found for transmission and reflection coefficients for both two-dimensional and circular geometries. Numerical solutions are compared with those for a method of weighted residuals in the form of a modified Galerkin method, showing good agreement.

Astley, R. J.↗
NASA NTRS1978

On the use of co-ordinate stretching in the numeral computation of high frequency scattering

The scattering of the sound of a jet engine by an airplane fuselage is modeled by solving the axially symmetric Helmholtz equation exterior to a long thin ellipsoid. The integral equation method based on the single layer potential formulation is used. A family of coordinate systems on the body is introduced and an algorithm is presented to determine the optimal coordinate system. Numerical results verify that the optimal choice enables the solution to be computed with a grid that is coarse relative to the wavelength.

Bayliss, A.↗
NASA NTRS1979

Numerical spatial marching techniques in duct acoustics

Direct calculation of the internal structure of a ducted noise source from farfield pressure measurements is regarded as an initial value problem, where the pressure and pressure gradient (farfield impedance) are assumed to be known along a line in the farfield. If pressure and impedance are known at the boundary of the farfield, the pressure can be uniquely determined in the vicinity of the inlet and inside the inlet ducting. A marching procedure is developed which, with this information obtained from measurements, enables a description of a ducted noise source. The technique uses a finite difference representation of the homogeneous Helmholtz equation.

Baumeister, K. J.↗
NASA NTRS1980

Acoustic propagation in rigid three-dimensional waveguides

The linear acoustic propagation in finite rigid three-dimensional waveguides is determined analytically using an eigenfunction expansion of the Helmholtz equation. The geometry considered consists of straight and circular bends of rectangular cross section with continuous interfaces (branches and sharp corners are excluded). The phenomena of resonance shift and relocation are explained for a bend-straight duct combination.

El-Raheb, M.↗