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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Analysis of opposed-jet hydrogen-air counter flow diffusion flame

An opposed-jet counterflow diffusion-flame configuration is considered for the analysis of a nitrogen-diluted hydrogen-air diffusion flame. A boundary-layer similarity solution is employed in order to reduce the governing equations to a set of equations in one independent variable. The equation set is written in the time-dependent form and solved by the finite-volume time-marching technique. This model uses detailed chemistry and accounts for the variations of Prandtl number and Lewis number as well as the effect of thermal diffusion on the flame. It is noted that a one-step model can predict several features of the flame, while the detailed-chemistry model can be used for fine-tuning the results. The present results indicate that thermal diffusion has negligible effect on the characteristics of the flame.

Ho, Y. H.↗

NASA NTRS1991

Unstructured Euler flow solutions using hexahedral cell refinement

An attempt is made to extend grid refinement into three dimensions by using unstructured hexahedral grids. The flow solver is developed using the TIGER (topologically Independent Grid, Euler Refinement) as the starting point. The program uses an unstructured hexahedral mesh and a modified version of the Jameson four-stage, finite-volume Runge-Kutta algorithm for integration of the Euler equations. The unstructured mesh allows for local refinement appropriate for each freestream condition, thereby concentrating mesh cells in the regions of greatest interest. This increases the computational efficiency because the refinement is not required to extend throughout the entire flow field.

Melton, John E.↗

NASA NTRS1991

Computation of steady and unsteady compressible quasi-axisymmetric vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux-difference splitting finite-volume scheme. The developed three-dimensional solver has been verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, isolated quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown have been captured. The problem has also been calculated using the Euler solver of the same code and the results are compared with those of the Navier-Stokes solver. The effect of the initial swirl has been tentatively studied.

Kandil, Osama A.↗

NASA NTRS1990

A transonic/supersonic CFD analysis of a generic fighter

Flow fields about a generic fighter model have been computed using FLO57, a three-dimensional, finite-volume Euler code. Computed pressure coefficients, forces, and moments at several Mach numbers - 0.6, 0.8, 1.2, 1.4, and 1.6 - are compared with wind tunnel data over a wide range of angles of attack in order to determine the applicability of the code for the analysis of fighter configurations. Two configurations were studied, a wing/body and a wing/body/chine. FLO57 predicted pressure distributions, forces, and moments well at low angles of attack, at which the flow was fully attached. The FLO57 predictions were also accurate for some test conditions once the leading-edge vortex became well established. At the subsonic speeds, FLO57 predicted vortex breakdown earlier than that seen in the experimental results. Placing the chine on the forebody delayed the onset of bursting and improved the correlation between numerical and experimental data at the subsonic conditions.

Goodsell, Aga M.↗

NASA NTRS1991

Numerical simulation of Shuttle ascent transonic flow using an unstructured-grid approach

An unstructured-grid, finite-volume method has been developed for simulating the inviscid flow over spacecrafts of realistic configuration. The grid generation is accomplished by a new technique on the basis of the advancing-front concept. This simple technique is shown to be equally as powerful for a complex multibody as for a single vehicle. Second- or third-order accuracy is obtained via an innovative interpolation procedure similar to the conventional MUSCL approach. This method has been applied to the Shuttle orbiter and a representative Shuttle launch vehicle consisting of the orbiter, the external tank, and the solid rocket boosters. A comparison is discussed between the present results and other results obtained from structured- and unstructured-grid methods.

Wey, T. C.↗

NASA NTRS1990

A numerical study of hypersonic stagnation heat transfer predictions at a coordinate singularity

The problem of grid induced errors associated with a coordinate singularity on heating predictions in the stagnation region of a three-dimensional body in hypersonic flow is examined. The test problem is for Mach 10 flow over an Aeroassist Flight Experiment configuration. This configuration is composed of an elliptic nose, a raked elliptic cone, and a circular shoulder. Irregularities in the heating predictions in the vicinity of the coordinate singularity, located at the axis of the elliptic nose near the stagnation point, are examined with respect to grid refinement and grid restructuring. The algorithm is derived using a finite-volume formulation. An upwind-biased total-variation diminishing scheme is employed for the inviscid flux contribution, and central differences are used for the viscous terms.

Grasso, Francesco↗

NASA NTRS1990

Numerical simulation of steady and unsteady asymmetric vortical flow

The unsteady, compressible, thin-layer, Navier-Stokes equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme. Locally conical flows are assumed and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round and sharp-edged delta wing flows.

Kandil, O. A.↗

NASA NTRS1991

Generalized conjugate-gradient methods for the Navier-Stokes equations

A generalized conjugate-gradient method is used to solve the two-dimensional, compressible Navier-Stokes equations of fluid flow. The equations are discretized with an implicit, upwind finite-volume formulation. Preconditioning techniques are incorporated into the new solver to accelerate convergence of the overall iterative method. The superiority of the new solver is demonstrated by comparisons with a conventional line Gauss-Siedel Relaxation solver. Computational test results for transonic flow (trailing edge flow in a transonic turbine cascade) and hypersonic flow (M = 6.0 shock-on-shock phenoena on a cylindrical leading edge) are presented. When applied to the transonic cascade case, the new solver is 4.4 times faster in terms of number of iterations and 3.1 times faster in terms of CPU time than the Relaxation solver. For the hypersonic shock case, the new solver is 3.0 times faster in terms of number of iterations and 2.2 times faster in terms of CPU time than the Relaxation solver.

Ajmani, Kumud↗

NASA NTRS1991

Computations of the three-dimensional flow and heat transfer within a coolant passage of a radial turbine blade

A numerical code is developed for computing three-dimensional, turbulent, compressible flow within coolant passages of turbine blades. The code is based on a formulation of the compressible Navier-Stokes equations in a rotating frame of reference in which the velocity dependent variable is specified with respect to the rotating frame instead of the inertial frame. The algorithm employed to obtain solutions to the governing equation is a finite-volume LU algorithm that allows convection, source, as well as diffusion terms to be treated implicitly. In this study, all convection terms are upwind differenced by using flux-vector splitting, and all diffusion terms are centrally differenced. This paper describes the formulation and algorithm employed in the code. Some computed solutions for the flow within a coolant passage of a radial turbine are also presented.

Shih, T. I.-P.↗

NASA NTRS1991

Supersonic quasi-axisymmetric vortex breakdown

An extensive computational study of supersonic quasi-axisymmetric vortex breakdown in a configured circular duct is presented. The unsteady, compressible, full Navier-Stokes (NS) equations are used for the present computational study. The NS equations are solved for quasi-axisymmetric flow using an implicit, upwind, flux-difference splitting, finite-volume scheme. The quasi-axisymmetric solutions are time accurate and are obtained by forcing the components of the flowfield vector to be equal on two axial planes, which are in close proximity of each other. The computational study addresses the effect of the Reynolds number, for laminar flows, on the evolution and persistence of vortex breakdown. The effect of boundary conditions at the duct exit on vortex breakdown is also studied. Finally, the effect of swirl ratio at the duct inlet is investigated.

Kandil, Osama A.↗

NASA NTRS1990

Navier-Stokes computations for aerodynamic configurations at high angles of attack

Navier-Stokes computations at high angles of attack over aerodynamic configurations are presented using an implicit finite-volume algorithm. The spatial differencing is upwind-biased for the convective and pressure terms and central for the shear stress and heat transfer terms. The equations are relaxed to steady state with a spatially factored implicit algorithm. In order to treat general geometries, a multiblock patched-grid framework is implemented. Applications and detailed comparisons with experimental data are made for two simple but representative geometric shapes: (1) a highly swept delta wing and (2) a prolate spheroid of 6:1 length-to-diameter. Recent extensions of the algorithm to compute the flow over an F-18 forebody-strake configuration are shown, including comparisons with wind tunnel and flight test results. Comparisons across the range of Reynolds number for all cases indicate that either low Reynolds number or high Reynolds number flows are simulated well, but that flows at intermediate Reynolds number cannot be simulated accurately without a detailed knowledge of the transitional zone between laminar and turbulent flow.

Thomas, J. L.↗

NASA NTRS1991

Time-dependent solutions of viscous incompressible flows in moving co-ordinates

A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

NASA NTRS1991

Simulation of 3-D shear flows around a nozzle-afterbody at high speeds

3D, compressible, unsteady, Reynolds-averaged Navier-Stokes equations are presently solved by a finite-volume and alternating-direction-implicit method in order to simulate supersonic and hypersonic turbulent shear flows. The effect of turbulence is incorporated via a modified Baldwin-Lomax eddy viscosity model which reflects the influence of high-speed compressibility, multiple walls, near-wall vortices, and turbulent memory effects, as well as local equilibrium effects. Attention is given to the simulation of the flow around the nozzle-afterbody of a generic, scramjet-propelled hypersonic vehicle; computed pressure distributions are consonant with experimental surface and off-surface flow surveys.

Baysal, Oktay↗

NASA NTRS1991

Turbulence and stall in plane diffusers - Computational study

The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.

Demuren, A. O.↗

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 NTRS1991

Application of program LAURA to thermochemical nonequilibrium flow through a nozzle

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Gnoffo, Peter A.↗

NASA NTRS1989

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

NASA NTRSDate not supplied

Aerothermal Shape Optimization of Actively-Cooled Battery Packs using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

thermal management↗