Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

CFD spinoff - Computational electromagnetics for radar cross section (RCS) studies

A finite-volume discretization procedure derived from proven CFD methods is used to solve the conservation form of the time-domain Maxwell's equations, in order to compute EM scattering from layered objects. This time-domain approach handles both single-frequency/continuous wave and broadband-frequency/pulse incident excitation. Arbitrarily shaped objects are modeled by means of a body-fitted coordinate transformation; complex internal/external structures with many material layers are treated through the implementation of a multizone framework capable of handling any type of zonal boundary condition. Results are presented for various two- and three-dimensional problems.

Shankar, Vijaya↗
NASA NTRS1990

An upwind approach to unsteady flowfield simulation

A numerical method to determine unsteady solutions of the laminar, perfect gas Navier-Stokes equations has been developed. The structured finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the viscous terms are computed by central differencing. Comparisons of numerical, analytical, and experimental results are made in two and three dimensions. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing.

Atwood, Christopher A.↗
NASA NTRS1990

Computation of vortex-dominated flow for a delta wing undergoing pitching oscillation

The conservative, unsteady Euler equations for the flow relative to a moving frame of reference are used to solve for the three-dimensional steady and unsteady flows around a sharp-edged delta wing. The resulting equations are solved by using an implicit, approximately factored, finite-volume scheme. Implicit second-order and explicit second- and fourth-order dissipations are added to the scheme. The boundary conditions are explicitly satisfied. The grid is generated by locally using a modified Joukowski transformation in crossflow planes at the grid-chord stations. The computational applications cover a steady flow around a delta wing, whose results serve as the initial conditions for the unsteady flow around a pitching delta wing at a large mean angle of attack. The steady results are compared with the experimental data, and the unsteady results are compared with results of a flux-difference splitting scheme.

Kandil, Osama A.↗
NASA NTRS1990

Numerical computation of viscous flow about unconventional airfoil shapes

A new two-dimensional computer code was developed to analyze the viscous flow around unconventional airfoils at various Mach numbers and angles of attack. The Navier-Stokes equations are solved using an implicit, upwind, finite-volume scheme. Both laminar and turbulent flows can be computed. A new nonequilibrium turbulence closure model was developed for computing turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present code was evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the present code in conjunction with the new nonequilibrium turbulence model gives excellent results.

Ahmed, S.↗
NASA NTRS1990

Interactive real time flow simulations

An interactive real time flow simulation technique is developed for an unsteady channel flow. A finite-volume algorithm in conjunction with a Runge-Kutta time stepping scheme was developed for two-dimensional Euler equations. A global time step was used to accelerate convergence of steady-state calculations. A raster image generation routine was developed for high speed image transmission which allows the user to have direct interaction with the solution development. In addition to theory and results, the hardware and software requirements are discussed.

Sadrehaghighi, I.↗
NASA NTRS1989

Application of unstructured grid methods to steady and unsteady aerodynamic problems

The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Batina, John T.↗
NASA NTRS1989

Navier-Stokes solutions for flows related to store separation

The objective is developing CFD capabilities to obtain solutions for viscous flows about generic configurations of internally and externally carried stores. The emphasis is placed on the supersonic flow regime with extensions being made to the transonic regime. The project is broken into four steps: (1) Cavity flows for internal carriage configurations; (2) High angle of attack flows, which may be experienced during the separation of the stores: (3) Flows about a body near a flat plate for external carriage configurations; and (4) Flows about a body inside or in the proximity of a cavity. Three-dimensional unsteady cavity flow solutions are obtained by an explicit, MacCormack algorithm, EMCAV3, for open, close, and transitional cavities. High angle of attack flows past cylinders are solved by an implicit, upwind algorithm. All the results compare favorably with the experimental data. For flows about multiple body configurations, the Chimera embedding scheme is modified for finite-volume and multigrid algorithms, MaGGiE. Then a finite volume, implicit, upwind, multigrid Navier-Stokes solver which uses on overlapped/embedded and zonal grids, VUMXZ3, is developed from the CFL3D code. Supersonic flows past a cylinder near a flat plate are computed using this code. The results are compared with the experimental data. Currently the VUMXZ3 code is being modified to accomplish step 4 of this project. Wind tunnel experiments are also being conducted for validation purposes.

Baysal, Oktay↗
NASA NTRS1989

A comparative study of Navier-Stokes codes for high-speed flows

A comparative study was made with four different codes for solving the compressible Navier-Stokes equations using three different test problems. The first of these cases was hypersonic flow through the P8 inlet, which represents inlet configurations typical of a hypersonic airbreathing vehicle. The free-stream Mach number in this case was 7.4. This 2-D inlet was designed to provide an internal compression ratio of 8. Initial calculations were made using two state-of-the-art finite-volume upwind codes, CFL3D and USA-PG2, as well as NASCRIN, a code which uses the unsplit finite-difference technique of MacCormack. All of these codes used the same algebraic eddy-viscosity turbulence model. In the experiment, the cowl lip was slightly blunted; however, for the computations, a sharp cowl leading edge was used to simplify the construction of the grid. The second test problem was the supersonic (Mach 3.0) flow in a three-dimensional corner formed by the intersection of two wedges with equal wedge angles of 9.48 degrees. The flow in such a corner is representative of the flow in the corners of a scramjet inlet. Calculations were made for both laminar and turbulent flow and compared with experimental data. The three-dimensional versions of the three codes used for the inlet study (CFL3D, USA-PG3, and SCRAMIN, respectively) were used for this case. For the laminar corner flow, a fourth code, LAURA, which also uses recently-developed upwind technology, was also utilized. The final test case is the two-dimensional hypersonic flow over a compression ramp. The flow is laminar with a free-stream Mach number of 14.1. In the experiment, the ramp angle was varied to change the strength of the ramp shock and the extent of the viscous-inviscid interaction. Calculations were made for the 24-degree ramp configuration which produces a large separated-flow region that extends upstream of the corner.

Rudy, David H.↗
NASA NTRS1990

Development of a time-dependent incompressible Navier-Stokes solver based on a fractional-step method

The development, validation and application of a fractional step solution method of the time-dependent incompressible Navier-Stokes equations in generalized coordinate systems are discussed. A solution method that combines a finite-volume discretization with a novel choice of the dependent variables and a fractional step splitting to obtain accurate solutions in arbitrary geometries was previously developed for fixed-grids. In the present research effort, this solution method is extended to include more general situations, including cases with moving grids. The numerical techniques are enhanced to gain efficiency and generality.

Rosenfeld, Moshe↗
NASA NTRS1990

Asymmetric supersonic flow around cones with noncircular sections

The unsteady, compressible, single and double thin-layer, Navier-Stokes equations are used to solve for steady and unsteady, asymmetric, supersonic flow around pointed bodies with noncircular sections at high incidence and zero side slip. The equations are solved by using an implicit, upwind, flux-difference splitting finite-volume scheme. Since the flow is locally conical, the solutions are presented on a cross-flow plane at the axial station of unit. The grid is generated by using a modified Joukowski transformation. The computational applications cover noncircular sections with elliptic and diamond shapes. Unsteady asymmetric vortex shedding has been captured at large angles of attack. It is shown that for the same flow conditions and same cross-section fineness ratio, the diamond-section cones develop less flow asymmetry than the elliptic-section cones. Passive control of flow asymmetry was also demonstrated computationally.

Kandil, Osama A.↗
NASA NTRS1990

A multiblock multigrid three-dimensional Euler equation solver

Current aerodynamic designs are often quite complex (geometrically). Flexible computational tools are needed for the analysis of a wide range of configurations with both internal and external flows. In the past, geometrically dissimilar configurations required different analysis codes with different grid topologies in each. The duplicity of codes can be avoided with the use of a general multiblock formulation which can handle any grid topology. Rather than hard wiring the grid topology into the program, it is instead dictated by input to the program. In this work, the compressible Euler equations, written in a body-fitted finite-volume formulation, are solved using a pseudo-time-marching approach. Two upwind methods (van Leer's flux-vector-splitting and Roe's flux-differencing) were investigated. Two types of explicit solvers (a two-step predictor-corrector and a modified multistage Runge-Kutta) were used with multigrid acceleration to enhance convergence. A multiblock strategy is used to allow greater geometric flexibility. A report on simple explicit upwind schemes for solving compressible flows is included.

Cannizzaro, Frank E.↗
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 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

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 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↗
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↗
DOE OSTI2025

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗