Results for “finite volume discretization” — Page 13

NASA NTRS1997

Complex Volume Grid Generation Through the Use of Grid Reusability

This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.

Alter, Stephen J.↗

NASA NTRS1972

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗

NASA NTRS1993

The 3-D numerical study of airflow in the compressor/combustor prediffuser and dump diffuser of an industrial gas turbine

This paper describes the 3D computations of a flow field in the compressor/combustor diffusers of an industrial gas turbine. The geometry considered includes components such as the combustor support strut, the transition piece and the impingement sleeve with discrete cooling air holes on its surface. Because the geometry was complex and 3D, the airflow path was divided into two computational domains sharing an interface region. The body-fitted grid was generated independently in each of the two domains. The governing equations for incompressible Navier-Stokes equations were solved using the finite volume approach. The results show that the flow in the prediffuser is strongly coupled with the flow in the dump diffuser and vice versa. The computations also revealed that the flow in the dump diffuser is highly nonuniform.

Agrawal, Ajay K.↗

NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗

NASA NTRS1976

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗

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 NTRSDate not supplied

Nodal Modeling of Tank Pressurization and Draining using a Multi-Node-Ullage Approach

Pressurized expulsion tests of liquid methane were modeled by a nodal code using a multinode ullage approach. Generalized Fluid System Simulation Program (GFSSP), a finite volume based nodal code was used to model the expulsion of liquid methane from a 1.52-meter (5 ft) diameter spherical tank by pressurizing with helium and gaseous methane. The purpose of the model was to estimate the amount of pressurant required and amount of condensation and evaporation of methane during the operation. The ullage was discretized into multiple nodes and each fluid node was connected to solid nodes. Both the fluid and solid nodes grow as the tank drains. The heat and mass transfer between fluid and solid node in the tank ullage was computed. The heat and mass transfer between the bottom ullage node and the liquid methane was also computed. The model predictions were compared with data from twelve test cases. The predicted pressurant consumption for helium pressurization compares with test data with a Mean Absolute Percent Error (MAPE) of 3 %. For autogenous pressurization, where gaseous methane was used to pressurize liquid methane, the predicted pressurant consumption compares with a MAPE of 6%.

Cryogenic Propellant↗

NASA NTRSDate not supplied

Nodal Modeling of Helium Pressurization and Autogenous Pressurization and Draining using a Multi-Node-Ullage Approach

Pressurized expulsion tests of liquid methane were modeled by a nodal code using a multinode ullage approach. Generalized Fluid System Simulation Program (GFSSP), a finite volume based nodal code was used to model the expulsion of liquid methane from a 1.52-meter (5 ft) diameter spherical tank by pressurizing with helium and gaseous methane. The purpose of the model was to estimate the amount of pressurant required and amount of condensation and evaporation of methane during the operation. The ullage was discretized into multiple nodes and each fluid node was connected to solid nodes. Both the fluid and solid nodes grow as the tank drains. The heat and mass transfer between fluid and solid node in the tank ullage was computed. The heat and mass transfer between the bottom ullage node and the liquid methane was also computed. The model predictions were compared with data from twelve test cases. The predicted pressurant consumption for helium pressurization compares with test data with a Mean Absolute Percent Error (MAPE) of 3%. For autogenous pressurization, where gaseous methane was used to pressurize liquid methane, the predicted pressurant consumption compares with a MAPE of 6%.

Nodal Model↗

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.

heat transfer↗

NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗

NASA NTRS2023

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

NASA NTRS2023

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

NASA NTRS1993

A comparison of high resolution upwind solvers on 3-D inviscid hypersonic flows

A detailed comparison of numerical results obtained by solving the Euler equations for the inviscid flow over delta wings in reentry configurations is presented. The investigation involves a side-by-side comparison between independently developed upwind Euler solvers at VKI (M3D) and NASA Johnson Space Center (E3D) using identical grids. In both solvers the governing equations are integrated by means of time marching finite volume shock capturing methods, based on a cell centered upwind evaluation of the cell face fluxes and nonlinear limiters. High resolution schemes are obtained via MUSCL characteristic variable extrapolation to ensure total variation diminishing (TVD) properties and therefore monotonic discontinuity capturing. Further comparisons are made with several published results by other authors based on both upwind and central discretizations. The present results contribute in asserting the high resolution upwind TVD schemes as the most reliable numerical technique to handle the strong discontinuities typical of high speed flows. The performances of the two upwind solvers are satisfactory and the numerical results in good agreement. However, the important issue of reaching a grid converged solution for the present complex three-dimensional inviscid flow problems was not achieved.

Manna, M.↗

NASA NTRS2013

High-Order Entropy Stable Formulations for Computational Fluid Dynamics

A systematic approach is presented for developing entropy stable (SS) formulations of any order for the Navier-Stokes equations. These SS formulations discretely conserve mass, momentum, energy and satisfy a mathematical entropy inequality. They are valid for smooth as well as discontinuous flows provided sufficient dissipation is added at shocks and discontinuities. Entropy stable formulations exist for all diagonal norm, summation-by-parts (SBP) operators, including all centered finite-difference operators, Legendre collocation finite-element operators, and certain finite-volume operators. Examples are presented using various entropy stable formulations that demonstrate the current state-of-the-art of these schemes.

Carpenter, Mark H.↗

NASA NTRS2014

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

NASA NTRS1992

Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of code, and assessment of performance, as well as demonstration of flexibility.

Elmiligui, Alaa↗

NASA NTRS2001

Visualization of CFD Results in Immersive Virtual Environments

An object-oriented event-driven immersive virtual environment (VE) is described for the visualization of computational fluid dynamics (CFD) results. The VE incorporates the following types of primitive software objects: interface objects, support objects, geometric entities, and finite elements. The fluid domain is discretized using either a multi-block structured grid or an unstructured finite element mesh. The VE allows natural 'fly-through' visualization of the model, the CFD grid, and the model's surroundings. In order to help visualize the flow and its effects on the model, the VE incorporates the following objects: stream objects (lines, surface-restricted lines. ribbons. and volumes); colored surfaces; elevation surfaces; surface arrows; global and local iso-surfaces; vortex cores; and separation/attachment surfaces and lines. Most of these objects can be used for dynamically probing the flow. Particles and arrow animations can be displayed on top of stream objects. Primitive response quantities as well as derived quantities can be used. A recursive tree search algorithm is used for real-time point and value search in the CFD grid.

Wasfy, Tamer M.↗

NASA NTRS1988

Solutions for hypersonic viscous flow over delta wings

A conical Navier-Stokes solver is developed to study hypersonic flow over delta wings. The differencing scheme used is van Leer flux vector splitting with MUSCL differencing. A flux limiter which involves a second difference in pressure was found to provide sharp, nonoscillatory shocks in windward and leeward flow regions. The governing equations are cast in finite-volume form and are solved using an explicit multistage method. A von Neumann stability analysis of a one-dimensional flux split model is used to predict the stability characteristics of the conical solver. The stability behavior of the scheme is found to strongly depend on the stage coefficients of the multistage scheme and on the order of accuracy of the upwind discretization.

Ruffin, Stephen M.↗