SEARCH · Search NASA

Results for “volume integral 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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Preconditioned conjugate gradient methods for the Navier-Stokes equations

A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Ajmani, Kumud↗

NASA NTRS1995

Three-dimensional unstructured method for flows past bodies in 6-DOF relative motion

A three dimensional, unstructured-mesh methodology was developed to simulate unsteady flows past bodies in relative motion, where the trajectory was determined from the instantaneous aerodynamics. The method coupled the equations of fluid flow and those of rigid-body dynamics, and captured the time-dependent interference between stationary and moving boundaries. The unsteady, compressible Euler equations were solved on dynamic, unstructured meshes by an explicit, finite-volume, upwind method. The grid adaptation was performed within a window placed around the moving body. The Euler equations of dynamics were solved by a Runge-Kutta integration scheme. The flow solver and the adaptation scheme were validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic pitching motion, then comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom (DOF) motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, K. P.↗

NASA NTRS1998

Refraction of Sound Emitted Near Solid Boundaries from a Sheared Jet

A mathematical model is developed to describe the sound emitted from an arbitrary point within a turbulent flow near solid boundaries. A unidirectional, transversely sheared mean flow is assumed, and the cross-section of the cold jet is of arbitrary shape. The analysis begins with Lilley's formulation of aerodynamic noise and, depending upon the specific model of turbulence used, leads via Fourier analysis to an expression for the spectral density of the intensity of the far-field sound emitted from a unit volume of turbulence. The expressions require solution of a reduced Green's function of Lilley's equation as well as certain moving axis velocity correlations of the turbulence. Integration over the entire flow field is required in order to predict the sound emitted by the complete flow. Calculations are presented for sound emitted from a plugflow jet exiting a semi-infinite flat duct. Polar plots of the far-field directivity show the dependence upon frequency and source position within the duct. Certain model problems are suggested to investigate the effect of duct termination, duct geometry, and mean flow shear upon the far-field sound.

Dill, Loren H.↗

NASA NTRS1986

The evolution of methods for noise prediction of high speed rotors and propellers in the time domain

Linear wave equation models which have been used over the years at NASA Langley for describing noise emissions from high speed rotating blades are summarized. The noise sources are assumed to lie on a moving surface, and analysis of the situation has been based on the Ffowcs Williams-Hawkings (FW-H) equation. Although the equation accounts for two surface and one volume source, the NASA analyses have considered only the surface terms. Several variations on the FW-H model are delineated for various types of applications, noting the computational benefits of removing the frequency dependence of the calculations. Formulations are also provided for compact and noncompact sources, and features of Long's subsonic integral equation and Farassat's high speed integral equation are discussed. The selection of subsonic or high speed models is dependent on the Mach number of the blade surface where the source is located.

Farassat, F.↗

NASA NTRS1983

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity is studied in detail via the eigenstrain approach. A complete formulation and a treatment of both types of eigenstrains for equivalence between the inhomogeneity problem and the inclusion problem are given. This approach is shown to be mathematically identical to other approaches such as the direct volume integral formulation. Expanding the eigenstrains and applied strains in the polynomial form in the position vector and satisfying the equivalence conditions at every point, the governing simultaneous algebraic equations for the unknown coefficients in the eigenstrain expansion are derived. The elastodynamic field outside an ellipsoidal inhomogeneity in a linear elastic isotropic medium is given as an example. The angular and frequency dependence of the induced displacement field, as well as the differential and total cross sections are formally given in series expansion form for the case of uniformly distributed eigenstrains.

Fu, L. S.↗

NASA NTRS1990

Conical Euler solution for a highly-swept delta wing undergoing wing-rock motion

Modifications to an unsteady conical Euler code for the free-to-roll analysis of highly-swept delta wings are described. The modifications involve the addition of the rolling rigid-body equation of motion for its simultaneous time-integration with the governing flow equations. The flow solver utilized in the Euler code includes a multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization on an unstructured mesh made up of triangles. Steady and unsteady results are presented for a 75 deg swept delta wing at a freestream Mach number of 1.2 and an angle of attack of 30 deg. The unsteady results consist of forced harmonic and free-to-roll calculations. The free-to-roll case exhibits a wing rock response produced by unsteady aerodynamics consistent with the aerodynamics of the forced harmonic results. Similarities are shown with a wing-rock time history from a low-speed wind tunnel test.

Lee, Elizabeth M.↗

NASA NTRS1993

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗

NASA NTRS1979

A calculation procedure for viscous flow in turbomachines, volume 1

A method for analyzing the nonadiabatic viscous flow through turbomachine rotors is presented. The field analysis is based upon the numerical integration of the full incompressible stream function vorticity form of the Navier-Stokes equations, together with the energy equation, over the rotor blade-to-blade stream channels. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system that suits the most complicated blade geometries. A numerical scheme is used to carry out the necessary integration of the elliptic governing equations. The flow characteristics within the rotor of a radial inflow turbine are investigated over a wide range of operating conditions. The calculated results are compared to existing experimental data. The flow in a radial compressor is analyzed in order to study the behavior of viscous flow in diffusing cascades. The results are compared qualitatively to known experimental trends. The solution obtained provides insight into the flow phenomena in this type of turbomachine. It is concluded that the method of analysis is quite general and gives a good representation of the actual flow behavior within turbomachine passages.

Khalil, I.↗

NASA NTRS1983

A generalized computer code for developing dynamic gas turbine engine models (DIGTEM)

This paper describes DIGTEM (digital turbofan engine model), a computer program that simulates two spool, two stream (turbofan) engines. DIGTEM was developed to support the development of a real time multiprocessor based engine simulator being designed at the Lewis Research Center. The turbofan engine model in DIGTEM contains steady state performance maps for all the components and has control volumes where continuity and energy balances are maintained. Rotor dynamics and duct momentum dynamics are also included. DIGTEM features an implicit integration scheme for integrating stiff systems and trims the model equations to match a prescribed design point by calculating correction coefficients that balance out the dynamic equations. It uses the same coefficients at off design points and iterates to a balanced engine condition. Transients are generated by defining the engine inputs as functions of time in a user written subroutine (TMRSP). Closed loop controls can also be simulated. DIGTEM is generalized in the aerothermodynamic treatment of components. This feature, along with DIGTEM's trimming at a design point, make it a very useful tool for developing a model of a specific turbofan engine.

Daniele, C. J.↗

NASA NTRS1984

A generalized computer code for developing dynamic gas turbine engine models (DIGTEM)

Daniele, C. J.↗

NASA NTRS1996

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗

NASA NTRS1987

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

NASA NTRS1991

Evaluation of Rotordynamic Coefficients of Look-through Labyrinths by Means of a Three Volume Bulk Flow Model

To describe the compressible, turbulent flow in a labyrinth seal, a three volume bulk flow model is presented. The conservation equations for mass, momentum, and energy are established in every control volume. A perturbation analysis is performed, yielding zeroth order equations for centric rotor position and first order equations describing the flow field for small rotor motions around the seal center. The equations are integrated numerically. From perturbation pressure, the forces on the shaft and the dynamic coefficients are calculated.

Nordmann, R.↗

NASA NTRS1991

Hypersonic cone flow predictions using an implicit upwind space-marching code

Solutions have been computed and results are presented for Problem 1, the case of Mach 9 transitional flow past a 7 deg half-angle cone at zero incidence. The solutions were computed using a code developed for the integration of the parabolized Navier-Stokes equations. The algorithm employed in the code is based on a Roe-type flux-difference-splitting scheme applied following a finite-volume approach. The basic algorithm has been modified to make it implicit and second-order accurate in the crossflow directions. Results are presented in terms of surface pressure and heat transfer as well as boundary layer profiles of pitot pressure, Mach number, and tangential velocity. The case was recalculated several times in an effort to determine sensitivities to such parameters as grid density, wall temperature, turbulence model parameters, as well as freestream expansion. Comparisons with the experimental data are presented and discussed.

Lawrence, Scott L.↗

NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗

NASA NTRS1983

Digital computer program for generating dynamic turbofan engine models (DIGTEM)

This report describes DIGTEM, a digital computer program that simulates two spool, two-stream turbofan engines. The turbofan engine model in DIGTEM contains steady-state performance maps for all of the components and has control volumes where continuity and energy balances are maintained. Rotor dynamics and duct momentum dynamics are also included. Altogether there are 16 state variables and state equations. DIGTEM features a backward-differnce integration scheme for integrating stiff systems. It trims the model equations to match a prescribed design point by calculating correction coefficients that balance out the dynamic equations. It uses the same coefficients at off-design points and iterates to a balanced engine condition. Transients can also be run. They are generated by defining controls as a function of time (open-loop control) in a user-written subroutine (TMRSP). DIGTEM has run on the IBM 370/3033 computer using implicit integration with time steps ranging from 1.0 msec to 1.0 sec. DIGTEM is generalized in the aerothermodynamic treatment of components.

Daniele, C. J.↗

NASA NTRS1986

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗

NASA NTRS2005

Permeable Surface Corrections for Ffowcs Williams and Hawkings Integrals

The acoustic prediction methodology discussed herein applies an acoustic analogy to calculate the sound generated by sources in an aerodynamic simulation. Sound is propagated from the computed flow field by integrating the Ffowcs Williams and Hawkings equation on a suitable control surface. Previous research suggests that, for some applications, the integration surface must be placed away from the solid surface to incorporate source contributions from within the flow volume. As such, the fluid mechanisms in the input flow field that contribute to the far-field noise are accounted for by their mathematical projection as a distribution of source terms on a permeable surface. The passage of nonacoustic disturbances through such an integration surface can result in significant error in an acoustic calculation. A correction for the error is derived in the frequency domain using a frozen gust assumption. The correction is found to work reasonably well in several test cases where the error is a small fraction of the actual radiated noise. However, satisfactory agreement has not been obtained between noise predictions using the solution from a three-dimensional, detached-eddy simulation of flow over a cylinder.

Lockard, David P.↗