SEARCH · Engineering Papers

Results for “Finite Volume”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1984

An implict LU scheme for the Euler equations applied to arbitrary cascades

An implicit scheme for solving the Euler equations is derived and demonstrated. The alternating-direction implicit (ADI) technique is modified, using two implicit-operator factors corresponding to lower-block-diagonal (L) or upper-block-diagonal (U) algebraic systems which can be easily inverted. The resulting LU scheme is implemented in finite-volume mode and applied to 2D subsonic and transonic cascade flows with differing degrees of geometric complexity. The results are presented graphically and found to be in good agreement with those of other numerical and analytical approaches. The LU method is also 2.0-3.4 times faster than ADI, suggesting its value in calculating 3D problems.

Buratynski, E. K.↗

NASA NTRS1984

Three-dimensional unsteady Euler equations solution using flux vector splitting

A method for numerically solving the three-dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of 2 and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far-field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three-dimensional Euler code and experimental data.

Whitfield, D. L.↗

NASA NTRS1984

Computational problems on composite grids

Most currently used algorithms for the numerical solution of the partial differential equations encountered in fluid flow problems can be implemented on composite grid systems. Finite volume formulations are easier to derive on composite grids, and may in principle be derived for partial differential equations of all types. Except when using Alternating Difference Implicit-type schemes, the overlapping of grids is an alternative to the more common grid construction procedure where grid lines continue smoothly from one subregion to the next. In the solution of model problems, the correct choice of an interpolation formula has been found able to reduce errors by a factor of two.

Mastin, C. W.↗

NASA NTRS1984

Conservative streamtube solution of steady-state Euler equations

This paper presents a new method for solving the steady state Euler equations. The method is similar to streamline curvature methods but has a conserative finite volume formulation which ensures correct shock capturing. Either wall position or wall pressure may be prescribed as boundary conditions, permitting both direct and inverse calculations. In supersonic applications the solution is obtained by space-marching while in subsonic and transonic applications iterative relaxation methods are used. Numerical results are given for: (1) supersonic diffuser with oblique shocks (direct calculation); (2) supersonic jet entering still reservoir (inverse calculation); (3) subsonic bump in a channel with 25 percent blockage (direct and inverse); (4) subsonic high-work turbine cascade (direct); and (5) transonic bump in a channel with 12 percent blockage (direct calculation).

Drela, M.↗

NASA NTRS1984

Investigation of the three-dimensional flow field within a transonic fan rotor - Experiment and analysis

An experimental investigation of the three-dimensional flow field through a low aspect ratio, transonic, axial flow fan rotor has been conducted, using an advanced laser anemometer (LA) system. Laser velocimeter measurements of the rotor flow field at the design operating speed and over a range of throughflow conditions are compared to analytical solutions. The numerical technique used herein yields the solution to the full, three-dimensional, unsteady Euler equations using an explicit time-marching, finite volume approach. The numerical analysis, when coupled with a simplified boundary layer calculation, generally yields good agreement with the experimental data. The test rotor has an aspect ratio of 1.56, a design total pressure ratio of 1.629 and a tip relative Mach number of 1.38. The high spatial resolution of the LA data matrix (9 radial x 30 axial x 50 blade-to-blade) permits details of the transonic flow field such as shock location, turning distribution, and blade loading levels to be investigated and compared to analytical results.

Pierzga, M. J.↗

NASA NTRS1985

A multistage time-stepping scheme for the Navier-Stokes equations

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Swanson, R. C.↗

NASA NTRS1985

Unsteady Flow in Multistage Turbines

The development of an efficient model for the response of a multistage turbine to either a total pressure or total temperature distortion is discussed. Each turbine blade row is modeled as an actuator disk and is often used to describe the flow through turbine and compressor blade rows. The flow approaching the blade row is assumed to be subsonic and inviscid. The distortion occurs at upstream infinity and is in the form of total pressure or total temperature. A solution was found to unsteady Euler equations by using a finite volume method. With prescribed inlet and boundary conditions flow variables such as density, pressure, and velocities can be calculated at cell centers in the entire flow region. The inlet and exit of the blade row coincide with lines of constant.

Surampudi, S. P.↗

NASA NTRS1985

A multistage time-stepping scheme for the Navier-Stokes equations

Swanson, R. C.↗

NASA NTRS1985

Numerical solution of the Euler equations for high-speed, blended wing-body configurations

Simulation of high-speed three-dimensional flow about blended wing-body combinations is investigated. A finite-volume explicit scheme with Runge-Kutta time integration is used to solve the compressible Euler equations in order to simulate the flow. The method, augmented by carefully chosen dissipative terms and convergence accelerators such as enthalpy damping and maximum local time-stepping, has been found to be very efficient in solving high-speed flows involving strong shocks. An analytic method for generating the body geometry and an algebraic method for quasi-three-dimensional grid generation are described. Results are presented for various blended wing-body configurations at different Mach numbers and angles of attack.

Moitra, A.↗

NASA NTRS1985

Euler calculations for multielement airfoils using Cartesian grids

A finite volume formulation for the Euler equations using Cartesian grids is presented and used to study complex two-dimensional configurations. The formulation extends methods developed for the potential equation to the Euler equations. Results using this approach for single element airfoils are shown to be competitive with and as accurate as other methods that employ mapped grids. Further, it is demonstrated that this method provides a simple and accurate procedure for solving flow problems involving multielement airfoils.

Clarke, D. K.↗

NASA NTRS1985

Solution method for a hovering helicopter rotor using the Euler equations

A method for the calculation of the flow field of a helicopter rotor blade in hover is presented. The approach uses a finite volume solution of the three dimensional Euler equations for the blade near field. In the Euler solver the velocity field is decomposed into two parts. One is the induced velocity of the vortex wake extending below the blade, found from a free wake calculation procedure. The other part is the unknown additional velocity field of the rotor blade. This approach eliminates numerical diffusion of the rolled up wake vorticity due to truncation error and artificial viscosity. Also, the effects of the far wake are included in the limited computational domain. Solutions are presented for an isolated wing and a model helicopter rotor and compared to experiment.

Roberts, T. W.↗

NASA NTRS1984

A kinetic energy analysis of the meso beta-scale severe storm environment

Analyses are performed of the meso beta-scale (20-200 km wavelengths and several hours to one-day periods) severe storm kinetic energy balance on the fifth day of the AVE SESAME campaign of May 1979. A 24-hr interval covering the antecedent, active and post-convective outbreak activity over Oklahoma are considered. Use is made of the kinetic energy budget equation (KEBE) for a finite volume in an isobaric coordinate system. Rawindsonde data with 75 km resolution were treated. The KEBE model covered changes in kinetic energy due to the cross contour flows, horizontal and vertical components of flux divergence, and volumic mass changes on synoptic and subsynoptic scales. The greatest variability was concentrated above 400 mb height and over the most intense storm activity. Energy was generated at the highest rates in divergence and decreased the most in convection. The meso beta-scale lacked sufficient resolution for analyzing mesoscale activity.

Fuelberg, H. E.↗

NASA NTRS1985

Application of Runge Kutta time marching scheme for the computation of transonic flows in turbomachines

Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The present study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions, and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in the solution procedure.

Subramanian, S. V.↗

NASA NTRS1985

Three-dimensional unsteady Euler equations solutions on dynamic grids

A method is presented for solving the three-dimensional unsteady Euler equations on dynamic grids based on flux vector splitting. The equations are cast in curvilinear coordinates and a finite volume discretization is used for handling arbitrary geometries. The discretized equations are solved using an explicit upwind second-order predictor corrector scheme that is stable for a CFL of 2. Characteristic variable boundary conditions are developed and used for unsteady impermeable surfaces and for the far-field boundary. Dynamic-grid results are presented for an oscillating air-foil and for a store separating from a reflection plate. For the cases considered of stores separating from a reflection plate, the unsteady aerodynamic forces on the store are significantly different from forces obtained by steady-state aerodynamics with the body inclination angle changed to account for plunge velocity.

Belk, D. M.↗

NASA NTRS1985

Exterior dimension of fat fractals

Geometric scaling properties of fat fractal sets (fractals with finite volume) are discussed and characterized via the introduction of a new dimension-like quantity which is called the exterior dimension. In addition, it is shown that the exterior dimension is related to the 'uncertainty exponent' previously used in studies of fractal basin boundaries, and it is shown how this connection can be exploited to determine the exterior dimension. Three illustrative applications are described, two in nonlinear dynamics and one dealing with blood flow in the body. Possible relevance to porous materials and ballistic driven aggregation is also noted.

Grebogi, C.↗

NASA NTRS1985

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

NASA NTRS1985

Assessment of wind tunnel corrections for multielement airfoils at transonic speeds

A finite volume formulation of the Euler equations using Cartesian grids is used to calculate the transonic flow over multielement airfoils and to use the resulting solutions to assess wall interference effects in wind tunnels. Available methods and recommendations for evaluating such effects, which are based on shifts in Mach number and angle of attack, are examined and the results are compared with measurements using the flapped supercritical SKF 1.1 airfoil. Based on the calculations, it is concluded that shifts in Mach number and angle of attack cannot by themselves account for viscous and wall effects on multielement airfoils at transonic speeds.

Gaffney, R. L., Jr.↗

NASA NTRS1985

Theoretical Problems in Materials Science

Interactions between theoretical physics and material sciences to identify problems of common interest in which some of the powerful theoretical approaches developed for other branches of physics may be applied to problems in materials science are presented. A unique structure was identified in rapidly quenched Al-14% Mn. The material has long-range directed bonds with icosahedral symmetry which does not form a regular structure but instead forms an amorphous-like quasiperiodic structure. Finite volume fractions of second phase material is advanced and is coupled with nucleation theory to describe the formation and structure of precipitating phases in alloys. Application of the theory of pattern formation to the problem of dendrite formation is studied.

Langer, J. S.↗