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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Flowfield computations and comparison with Shuttle aerodynamic data

An in-house developed flow solver, E3D, has been applied to investigate the flow field around the Shuttle Orbiter over an angle-of-attack range of 0 deg to 60 deg at Mach 3.5 and 10.0. The 3D Euler equations are integrated by means of a time-marching finite-volume shock-capturing method, based on cell-centered and upwind evaluation of the cell face fluxes with van Leer or Roe's flux-splitting techniques. A comparison of CFD, wind tunnel, and aeronautical data base is presented. Solutions are in very good agreement with available experimental data. Through the comparisons, it is demonstrated that the E3D code is capable of evaluating total aerodynamic performance.

Ma, E. C.↗

NASA NTRS1991

Effects of leading-edge flap oscillation on unsteady delta wing flow and rock control

The isolated and interdisciplinary problems of unsteady fluid dynamics and rigid-body dynamics and control of delta wings with and without leading-edge flap oscillation are considered. For the fluid dynamics problem, the unsteady, compressible, thin-layer Navier-Stokes (NS) equations, which are written relative to a moving frame of reference, are solved along with the unsteady, linearized, Navier-displacement (ND) equations. The NS equations are solved for the flowfield using an implicit finite-volume scheme. The ND equations are solved for the grid deformation, if the leading-edge flaps oscillate, using an ADI scheme. For the dynamics and control problem, the Euler equation of rigid-body rolling motion for a wing and its flaps are solved interactively with the fluid dynamics equations for the wing-rock motion and subsequently for its control. A four-stage Runge-Kutta scheme is used to explicitly integrate the dynamics equation.

Kandil, Osama A.↗

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

Io's volcanic and sublimation atmospheres

Fully 3D axisymmetric gasdynamic equations simulating SO2 and H2S frost sublimation and SO2 dayside and nightside volcanic atmospheres on Io are numerically solved, using a time-explicit finite-volume formulation. Both the sublimation and volcanic atmospheres generate horizontal supersonic winds away from the subsolar point or the volcanic vent. While the sublimation atmosphere is primarily driven by horizontal pressure gradients determined by surface temperatures, the volcanic atmosphere is driven by pressure gradients that are determined by the source rate. Sublimation and condensation produce patterns of surface deposits which are characteristic of the two types of atmospheres. The volcanic model is quantitatively consistent with Voyager observations of ring deposits.

Moreno, Miguel A.↗

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 NTRS1991

Time-dependent aerodynamic analysis of ducted and unducted propfans at angle of attack

A three-dimensional unsteady aerodynamic analysis is described for predicting the time-dependent flow about ducted and unducted propfans operating at angle of attack. Although the freestream is assumed to be uniform, the flow relative to the rotating blades varies with circumferential position, resulting in an inherent unsteadiness due to the nonaxial inflow. The time-dependent Euler equations are solved utilizing a Runge-Kutta time-stepping scheme. The analysis is based on a finite-volume discretization employing a multiple-block grid network. To permit the use of large calculation time steps, an implicit residual smoothing scheme previously tested for unsteady flow calculations in two dimensions is extended to three spatial dimensions. For unducted propfans, a single H-type grid block is used for each blade passage to determine the time-periodic flowfield. For ducted propfans (ultra-high bypass fans) a body-centered C-type grid is wrapped about the cowl to improve the accuracy of the analysis in the high gradient flow region near the cowl leading edge. Numerical results are compared with available data for both ducted and unducted propfans operating at angle of attack.

Hall, Edward J.↗

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

Time accurate Navier-Stokes computations of solid rocket motor internal field joint cavity flows

The low-speed flow of combustion gases in the Space Shuttle Solid Rocket Motor (SRM) following ignition are investigated with a time-accurate CFD technique. The flow of the combustion gases is examined in the forward joint section of the SRM at 80 s after ignition with attention given to the amplification of existing 'organ-pipe' modes of 15, 30, and 45 Hz resulting from frequency lock-in with vortex-shedding modes. The CFD code employed is a multizone compressible Navier-Stokes flow solver using a finite-volume, upwind-biased, total-variation-diminishing scheme. The vortex-generated oscillations generated by the inhibitor and the cavity at the forward field joint are described in terms of amplitude and frequency and related to the flow dynamics in the SRM. Oscillations are shown to occur spontaneously when the trailing edge of the cavity at the joint is sharp and can be amplified if the unburned portion of the inhibitor protrudes into the flow.

Farr, Rebecca A.↗

NASA NTRS1991

Aerobrake aerothermodynamic environment predictions in support of the Aeroassist Flight Experiment

A computational fluid dynamic (CFD) simulation of the flow past an Aeroassist Flight Experiment (AFE) aerobrake under wind tunnel test conditions has been made. The Navier-Stokes equations in three dimensions were solved numerically using a finite-volume, implicit approach. Comparisons with experimental data include surface pressures and heat transfer rates and aerodynamic coefficients. Agreement with experiment is shown to be excellent. The dependence of the solution on the computational grid is explored. The present work is preliminary to simulation of the flow past the AFE under flight conditions, which at present cannot be duplicated with ground-based experimental facilities.

Fay, John F.↗

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 NTRS1991

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗

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↗

NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗