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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗
NASA NTRS1996

Forced Rolling Oscillation of a 65 deg-Delta Wing in Transonic Vortex-Breakdown Flow

Unsteady, transonic, vortex dominated flow over a 65 deg. sharp-edged, cropped-delta wing of zero thickness undergoing forced rolling oscillations is investigated computationally. The wing angle of attack is 20 deg. and the free stream Mach number and Reynolds number are 0.85 and 3.23 x 10(exp 6), respectively. The initial condition of the flow is characterized by a transverse terminating shock which induces vortex breakdown of the leading edge vortex cores. The computational investigation uses the time accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux difference splitting, finite-volume scheme. While the maximum roll amplitude is kept constant at 4.0 deg., both Reynolds number and roll frequency are varied covering three cases of forced sinusoidal rolling. First, the Reynolds number is held at 3.23 x 10(exp 6) and the wing is forced to oscillate in roll around the axis of geometric symmetry at a reduced frequency of 2(pi). Second, the Reynolds number is reduced to 0.5 x 10(exp 6) to observe the effects of added viscosity on the vortex breakdown. Third, with the Reynolds number held at 0.5 x 10(exp 6), the roll frequency is reduced to 1(pi) to complete the study.

Menzies, Margaret A.↗
NASA NTRS1995

A Solution Adaptive Structured/Unstructured Overset Grid Flow Solver with Applications to Helicopter Rotor Flows

This paper summarizes a method that solves both the three dimensional thin-layer Navier-Stokes equations and the Euler equations using overset structured and solution adaptive unstructured grids with applications to helicopter rotor flowfields. The overset structured grids use an implicit finite-difference method to solve the thin-layer Navier-Stokes/Euler equations while the unstructured grid uses an explicit finite-volume method to solve the Euler equations. Solutions on a helicopter rotor in hover show the ability to accurately convect the rotor wake. However, isotropic subdivision of the tetrahedral mesh rapidly increases the overall problem size.

Duque, Earl P. N.↗
NASA NTRS1996

Assessment of an Unstructured-Grid Method for Predicting 3-D Turbulent Viscous Flows

A method Is presented for solving turbulent flow problems on three-dimensional unstructured grids. Spatial discretization Is accomplished by a cell-centered finite-volume formulation using an accurate lin- ear reconstruction scheme and upwind flux differencing. Time is advanced by an implicit backward- Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the sublayer region of the boundary layer. A systematic assessment of the method is presented to devise guidelines for more strategic application of the technology to complex problems. The assessment includes the accuracy In predictions of skin-friction coefficient, law-of-the-wall behavior, and surface pressure for a flat-plate turbulent boundary layer, and for the ONERA M6 wing under a high Reynolds number, transonic, separated flow condition.

Frink, Neal T.↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗
NASA NTRS1996

Natural Rolling Responses of a Delta Wing in Transonic and Subsonic Flows

The unsteady, three-dimensional, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the natural rolling response of slender delta wings of zero thickness at moderate to high angles of attack, to transonic and subsonic flows. The governing equations of fluid flow and dynamics of the present multi-disciplinary problem are solved using the time-accurate solution of the NS equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The main focus is to analyze the effect of Mach number and angle of attack on the leading edge vortices and their breakdown, the resultant rolling motion, and overall aerodynamic response of the wing. Three cases demonstrate the natural response of a 65 deg swept, cropped delta wing in a transonic flow with breakdown of the leading edge vortices and an 80 deg swept delta wing in a subsonic flow undergoing either damped or self-excited limit-cycle rolling oscillations as a function of angle of attack. Comparisons with an experimental investigation completes this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret A.↗
NASA NTRS1996

Far-Field Turbulent Vortex-Wake/Exhaust Plume Interaction for Subsonic and HSCT Airplanes

Computational study of the far-field turbulent vortex-wake/exhaust plume interaction for subsonic and high speed civil transport (HSCT) airplanes is carried out. The Reynolds-averaged Navier-Stokes (NS) equations are solved using the implicit, upwind, Roe-flux-differencing, finite-volume scheme. The two-equation shear stress transport model of Menter is implemented with the NS solver for turbulent-flow calculation. For the far-field study, the computations of vortex-wake interaction with the exhaust plume of a single engine of a Boeing 727 wing in a holding condition and two engines of an HSCT in a cruise condition are carried out using overlapping zonal method for several miles downstream. These results are obtained using the computer code FTNS3D. The results of the subsonic flow of this code are compared with those of a parabolized NS solver known as the UNIWAKE code.

Kandil, Osama A.↗
NASA NTRS1996

Computation of Radiation Heat Transfer in Aeroengine Combustors

In this report the highlights of the research completed for the NASA are summarized. This research has been completed in the form of two Ph.D. theses by Chai (1994) and Parthasarathy (1996). Readers are referred to these theses for a complete details of the work and lists of references. In the following sections, first objectives of this research are introduced, then the finite-volume method for radiation heat transfer is described, and finally computations of radiative heat transfer in non-gray participating media is presented.

Patankar, S. V.↗
NASA NTRS1994

Hyperbolic Prismatic Grid Generation and Solution of Euler Equations on Prismatic Grids

A hyperbolic grid generation method is used to generate prismatic grids and an approach using prismatic grids to solve the Euler equations is presented. The theory of the stability and feasibility of the hyperbolic grid generation method is presented. The hyperbolic grid generation method of Steger et al for structured grids is applied to a three dimensional triangularized surface definition to generate a grid that is unstructured on each successive layer. The grid, however, retains structure in the body-normal direction and has a computational cell shaped like a triangular prism. In order to take advantage of the structure in the normal direction, a finite-volume scheme that treats the unknowns along the normal direction implicitly is introduced and the flow over a sphere is simulated.

Pandya, S. A.↗
NASA NTRS1994

The Monostatic/Bistatic Approximation

Many Radar Cross Section (RCS) prediction codes are limited to one monostatic return per run. However, such codes can calculate multiple bistatic returns per incident angle for a relatively small amount of additional computer resources. This note describes a method of using bistatic returns to generate multiple monostatic predictions for each incident angle computed. Typical results are presented and show the accuracy is initially good and then degrades as the separation angle between the incident and viewing angles becomes large. Introduction Since 1990, the monostatic/bistatic approximation has been used to reduce the number of runs required by finite-volume time-domain (FVTD) codes for making RCS versus azimuth plots. This approximation was spawned by the observation of a range test where the transmit and receive antennas were separated by a few degrees to prevent cross talk between the antennas. The measurements from this range are presented as monostatic RCS rather than bistatic: RCS. The procedure of reporting experimental bistatic RCS as the monostatic RCS at the angle bisecting the transmit and receive antennas was extended to FVTD codes and produces excellent results.

Schuh, Michael J.↗
NASA NTRS1995

Fractional Step and Pseudocompressibility Methods for the Solution of Incompressible Navier-Stokes Equations

The fractional step and the pseudocompressibility methods for the solution of the incompressible Navier-Stokes equations are outlined. The fractional step method is based on finite-volume formulation and uses the pressure and the volume fluxes across the faces of each cell as dependent variables. The momentum equations are solved implicitly and the Poisson equation for the pressure is solved by using the multigrid method. The pseudocompressibility approach uses an implicit-higher-order-upwind differencing scheme for the convective terms together with the Gauss-Seidel line relaxation method. The dependent variables in the pseudocompressibility approach are the pressure and the cartesian velocity components in unstaggered mesh orientation. The 90-degree square duct flow, the wing-tip vortex wake flow and unsteady turbulent flows over an oscillating NACA 0015 airfoil are computed using both the fractional step and the pseudocompressibility methods. The results obtained from two different schemes are compared against experimental measurements.

Kiris, Cetin↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS1996

Challenges and Accomplishments of Viscous Incompressible Flow Computations in Aerospace and Beyond

Recent progress in incompressible Navier-Stokes solution methods will be presented. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. Both steady- and unsteady-solution algorithms and their salient features are discussed. A summary of our recent progress in flow solver development is given followed by numerical studies of a few example problems of our current interest. Solvers discussed here are based on structured-grid approach using finite-difference or finite-volume frame work. This short course will be delivered in three one-our lectures. The material in the course are collected from the work performed by the Incompressible Navier-Stokes group at NASA Ames Research Center over the past several years, and can be found in our publications widely disseminated in the US and abroad. This short course is sponsored by AGARD Consultant and Exchange Program under Support Project P-110.

Kwak, Dochan↗
NASA NTRS1992

Characteristics of the Shuttle Orbiter Leeside Flow During A Reentry Condition

A study of the leeside flow characteristics of the Shuttle Orbiter is presented for a reentry flight condition. The flow is computed using a point-implicit, finite-volume scheme known as the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). LAURA is a second-order accurate, laminar Navier-Stokes solver, incorporating finite-rate chemistry with a radiative equilibrium wall temperature distribution and finite-rate wall catalysis. The resulting computational solution is analyzed in terms of salient flow features and the surface quantities are compared with flight data.

Kleb, William L.↗
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↗
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↗