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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Reduction of Ducted Fan Engine Noise Via A Boundary Integral Equation Method

The development of a Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is discussed. The method is motivated by the need for an efficient and versatile computational tool to assist in parametric noise reduction studies. In this research, the work in reference 1 was extended to include passive noise control treatment on the duct interior. The BEM considers the scattering of incident sound generated by spinning point thrust dipoles in a uniform flow field by a thin cylindrical duct. The acoustic field is written as a superposition of spinning modes. Modal coefficients of acoustic pressure are calculated term by term. The BEM theoretical framework is based on Helmholtz potential theory. A boundary value problem is converted to a boundary integral equation formulation with unknown single and double layer densities on the duct wall. After solving for the unknown densities, the acoustic field is easily calculated. The main feature of the BIEM is the ability to compute any portion of the sound field without the need to compute the entire field. Other noise prediction methods such as CFD and Finite Element methods lack this property. Additional BIEM attributes include versatility, ease of use, rapid noise predictions, coupling of propagation and radiation both forward and aft, implementable on midrange personal computers, and valid over a wide range of frequencies.

Tweed, J.↗
NASA NTRS2002

Analytical Modal Analysis for Thin-Film Flat Lenses

Due to strong potential applications and more demanding requirements imposed upon thin-film structures for space deployable, there has been increasing research and development activities during recent years in the field of vibration analysis of these types of structures. Moreover, interests in employing these structural components have received renewed emphasis in recent years within NASA and the Air Force. This is due to their inherent lightweight, low packaging and launch volume, and relative simplicity of deployment. Among the potential mission concepts for which these structural elements are included, one can mention solar sails, space solar power generation systems, solar thermal propulsion vehicles, large space telescopes, and inflatable communication antennas. This paper presents analytical procedures to determine vibration and physical characteristics of thin film lenses with circular and elliptical shapes membranes considered in design of a solar concentrator. In general, three methods are used to obtain approximate solutions of Helmholtz boundary value problems. One method requires that solution satisfy the differential equation exactly and the boundary condition approximately. Another method demands a solution that satisfies the boundary conditions exactly and the governing equations approximately. The third method sees a solution that satisfies both the governing equation and boundary conditions approximately. Extensive reviews of vibrations of membrane and plates are provided by Leissa and Mazumdar.

Hamid R. Hamidzadeh↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
NASA NTRS1950

The Boundary Layers in Fluids with Little Friction

The vortices forming in flowing water behind solid bodies are not represented correctly by the solution of the potential theory nor by Helmholtz's jets. Potential theory is unable to satisfy the condition that the water adheres at the wetted bodies, and its solutions of the fundamental hydrodynamic equations are at variance with the observation that the flow separates from the body at a certain point and sends forth a highly turbulent boundary layer into the free flow. Helmholtz's theory attempts to imitate the latter effect in such a way that it joins two potential flows, jet and still water, nonanalytical along a stream curve. The admissibility of this method is based on the fact that, at zero pressure, which is to prevail at the cited stream curve, the connection of the fluid, and with it the effect of adjacent parts on each other, is canceled. In reality, however, the pressure at these boundaries is definitely not zero, but can even be varied arbitrarily. Besides, Helmholtz's theory with its potential flows does not satisfy the condition of adherence nor explain the origin of the vortices, for in all of these problems, the friction must be taken into account on principle, according to the vortex theorem.

Blasius, H.↗
NASA NTRS1977

On the origin of the magnetic field in type-1 comet tails - A reply to Ershkovich

Ershkovich (1976) employed the observed geometry of the type 1 tail of comet Arend-Roland 1957 III to show that any magnetic field present there cannot be of 'internal' origin but must have been 'captured' from the magnetized solar wind. Such an 'external' origin for the magnetic field in type 1 comet tails is not disputed, but several other points raised by Ershkovich are discussed. These concern the generation of substantial magnetic fields within the coma during periods when the comet's heliocentric distance is typically not greater than 1 AU, Ershkovich's use of the pressure balance across the tail 'tangential discontinuity' to estimate that the 'captured' magnetic field is of the order of the solar wind-field, his interpretation of certain observations of motions in the tail of comet Arend-Roland, and his assertion that the distant comet tail will be transformed into a turbulent wake due to the Kelvin-Helmholtz instability. It is argued that the 'captured' magnetic field should be substantially amplified relative to the solar-wind field, that the pressure-balance equation is of limited use, and that the real configuration of a type 1 tail must be treated in detail before any definite conclusions can be reached about the Kelvin-Helmholtz instability.

Mendis, D. A.↗
NASA NTRS1972

Design and testing of a magnetic suspension and damping system for a space telescope

The basic equations of motion are derived for a two dimensional, three degree of freedom simulation of a space telescope coupled to a spacecraft by means of a magnetic suspension and isolation system. The system consists of paramagnetic or ferromagnetic discs confined to the magnetic field between two Helmholtz coils. Damping is introduced by varying the magnetic field in proportion to a velocity signal derived from the telescope. The equations of motion are nonlinear, similar in behavior to the one-dimensional Van der Pol equation. The computer simulation was verified by testing a 264-kilogram air bearing platform which simulates the telescope in a frictionless environment. The simulation demonstrated effective isolation capabilities for disturbance frequencies above resonance. Damping in the system improved the response near resonance and prevented the build-up of large oscillatory amplitudes.

Ockman, N. J.↗
NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗
NASA NTRS1988

Exact solutions and low-frequency instability of the adiabatic auroral arc model

The adiabatic auroral arc model couples a kinetic theory parallel current driven by mirror forces to horizontal ionospheric currents; the resulting equations are nonlinear. Some exact stationary solutions to these equations, some of them based on the Liouville equation, are developed, with both latitudinal and longitudinal spatial variations. These Liouville equation exact solutions are related to stability boundaries of low-frequency instabilities such as Kelvin-Helmholtz, as shown by a study of a simplified model.

Cornwall, John M.↗
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
NASA NTRS1988

The nonlinear breakup of a magnetic layer - Instability to interchange modes

Motivated by considerations of the solar toroidal magnetic field, the behavior of a layer of uniform magnetic field embedded in a convectively stable atmosphere is studied. Since the field can support extra mass, such a configuration is top-heavy and thus instabilities of the Rayleigh-Taylor type can occur. For both static and rotating basic states, the evolution of the interchange modes (no bending of the field lines) is followed by integrating numerically the nonlinear compressible MHD equations. The initial Rayleigh-Taylor instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The subsequent motions are determined primarily by the strong interactions between vortices which are responsible for the rapid disruption of the magnetic layer.

Cattaneo, F.↗
NASA NTRS1974

Effect of the cluster integrals on three particles on the calculated electron density of a hydrogen plasma

The effect of the calculation of the cluster integrals on three particles is analyzed and evaluated for a hydrogen plasma where a pairwise-additive hard sphere-Coulomb potential is assumed. The Mayer cluster integral method was used to calculate the Helmholtz free energy which was then applied to the calculation of the electron number density through an iterative technique using a corrected Saha equation. It is seen that the three particle integrals provide a substantial correction to the calculations in the low energy-high density region of the hydrogen plasma.

Mcintyre, R. G.↗
NASA NTRS1991

Kelvin-Helmholtz instability in an Alfven resonant layer of a solar coronal loop

A Kelvin-Helmholtz instability has been identified numerically on an azimuthally symmetric Alfven resonant layer in an axially bounded, straight cylindrical coronal loop. The set of equations is solved numerically as an initial value problem. The linear growth rate of this instability is shown to be approximately proportional to the Alfven driving amplitude and inversely proportional to the width of the Alfven resonant layer. It is also shown that the linear growth rate increases linearly with m - 1 up to a certain m, reaches its maximum value for the mode whose half wavelength is comparable to the Alfven resonant layer width, and decreases at higher azimuthal mode number.

Uchimoto, E.↗
NASA NTRS1972

Rocket Combustion Instability Suppression With Cavity Resonators

Properly tuned Helmholtz resonator damping devices are an effective method for eliminating high frequency combustion instability in rocket engine thrust chambers. A program was implemented with the objective of extending Helmholtz resonator theory to thrust chamber conditions to establish a standard design procedure. By using systematic cold-flow testing and hot firings, empirical design equations were formulated. An effort is currently underway to extend the design theory to other types of acoustical damping devices such as quarter-wave tubes and slots.

R H Counts↗
NASA NTRS1992

Drift Kelvin-Helmholtz instabilities in space plasmas

Drift Kelvin-Helmholtz instabilities of a finite-beta plasma in equilibrium electric and magnetic fields which are perpendicular to each other are studied using two fluid equations. Three types of these instabilities are considered including the magnetosonic instability of a finite beta-homogeneous plasma, the electrostatic drift instability of an inhomogeneous low-beta plasma, and the magneto-acoustic instability of a high-beta inhomogeneous isothermal plasma. It is shown that the electric field has either stabilizing or destabilizing effect depending on conditions under consideration.

Sharma, Avadhesh C.↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
NASA NTRS1985

Pulsar magnetospheres in binary systems

The criterion for stability of a tangential discontinuity interface in a magnetized, perfectly conducting inviscid plasma is investigated by deriving the dispersion equation including the effects of both gravitational and centrifugal acceleration. The results are applied to neutron star magnetospheres in X-ray binaries. The Kelvin-Helmholtz instability appears to be important in determining whether MHD waves of large amplitude generated by instability may intermix the plasma effectively, resulting in accretion onto the whole star as suggested by Arons and Lea and leading to no X-ray pulsar behavior.

Ershkovich, A. I.↗
NASA NTRS2008

Application of Fast Multipole Methods to the NASA Fast Scattering Code

The NASA Fast Scattering Code (FSC) is a versatile noise prediction program designed to conduct aeroacoustic noise reduction studies. The equivalent source method is used to solve an exterior Helmholtz boundary value problem with an impedance type boundary condition. The solution process in FSC v2.0 requires direct manipulation of a large, dense system of linear equations, limiting the applicability of the code to small scales and/or moderate excitation frequencies. Recent advances in the use of Fast Multipole Methods (FMM) for solving scattering problems, coupled with sparse linear algebra techniques, suggest that a substantial reduction in computer resource utilization over conventional solution approaches can be obtained. Implementation of the single level FMM (SLFMM) and a variant of the Conjugate Gradient Method (CGM) into the FSC is discussed in this paper. The culmination of this effort, FSC v3.0, was used to generate solutions for three configurations of interest. Benchmarking against previously obtained simulations indicate that a twenty-fold reduction in computational memory and up to a four-fold reduction in computer time have been achieved on a single processor.

Dunn, Mark H.↗
NASA NTRS1977

Increase in effectiveness of low frequency acoustic liners by use of coupled Helmholtz resonators

Coupling of Helmholtz resonators in a low-frequency absorber array was studied as a means for increasing the effectiveness for absorbing low-frequency core engine noise. The equations for the impedance of the coupled-resonator systems were developed in terms of uncoupled-resonator parameters, and the predicted impedance for a parallel-coupled scheme is shown to compare favorably with measurements from a test model. In addition, attenuation measurements made in a flow duct on test coupled-resonator panels are shown to compare favorably with predicted values. Finally, the parallel-coupled concept is shown to give significantly more attenuation than that of a typical uncoupled resonator array of the same total volume.

Dean, L. W.↗