Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite volume method”

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 361 records · Page 20

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

Viscous simulation method for unsteady flows past multicomponent configurations

The present algorithm for the numerical simulation of flows about complex configurations (whose multiple and nonsimilar components have arbitraty geometries) employs a hybridization of the domain decomposition techniques for grid generation as well as to reduce computer-memory requirements. A fully vectorized, finite-volume, upwind-biased, approximately factored multigrid method is used to solve 3D Reynolds-averaged unsteady and compressible Navier-Stokes equations simulating supersonic flows past an ogive-nose-cylinder near or within a cavity. The time-averaged surface pressures obtained compare favorably with the wind tunnel data.

Fouladi, Kamran↗

Transonic flow simulation of prop-fan configurations

This paper presents the results of a study that used a three-dimensional transonic analysis method to analyze flow around a wing-mounted prop-fan configuration. The transonic analysis method employs a cell-oriented finite volume approach to solve the full potential equation in conservative form over a numerically generated surface-fitted grid. The propeller slipstream effects are simulated using a linearized transpiration boundary condition. The results of analyses of flow over clean wing/body, wing/body/nacelle, and wing/body nacelle with slipstream are compared to test data.

Samant, S. S.↗

Numerical solutions of the complete Navier-Stokes equations

The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.

Hassan, H. A.↗

Navier-Stokes analysis of transonic cascade flow

A new kind of C-type grid is proposed, this grid is non-periodic on the wake and allows minimum skewness for cascades with high turning and large camber. Reynolds-averaged Navier-Stokes equations are solved on this type of grid using a finite volume discretization and a full multigrid method which uses Runge-Kutta stepping as the driving scheme. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A detailed numerical study is proposed for a highly loaded transonic blade. A grid independence analysis is presented in terms of pressure distribution, exit flow angles, and loss coefficient. Comparison with experiments clearly demonstrates the capability of the proposed procedure.

Arnone, A.↗

Navier-Stokes solution of transonic cascade flows using nonperiodic C-type grids

A new kind of C-type grid is proposed, this grid is non-periodic on the wake and allows minimum skewness for cascades with high turning and large camber. Reynolds-averaged Navier-Stokes equations are solved on this type of grid using a finite volume discretization and a full multigrid method which uses Runge-Kutta stepping as the driving scheme. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A detailed numerical study is proposed for a highly loaded transonic blade. A grid independence analysis is presented in terms of pressure distribution, exit flow angles, and loss coefficient. Comparison with experiments clearly demonstrates the capability of the proposed procedure.

Arnone, Andrea↗

Freestream capturing for moving coordinates in three dimensions

A detailed formula is developed that can be used in both finite-volume (FV) and finite-difference (FD) methods for constructing freestream capturing metrics in space and time. It is shown that, considering an FV cell on the FD grid, the freestream capturing metrics in space and time can be constructed from the FD formulation. The approach is costly but guarantees the global conservation for an arbitrary motion of the grid.

Obayashi, Shigeru↗

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↗

A numerical study of the effect of geometry variation, turbulence models, and dissipation on the flow past control surfaces

The steady 3D viscous flow past the ONERA M6 wing and a slender delta wing-body with trailing edge control surfaces has been computed. A cell-centered finite-volume Navier-Stokes patched zonal method has been used for the numerical simulation. Both diagonalized and LUSGS schemes have been implemented. Besides the standard nonplanar zonal interfacing techniques, a new virtual zone capability has been employed. For code validation, the transonic flow past the ONERA M5 wing is calculated for angles-of-attack of 3.06 deg and 5.06 deg and compared with the available experiments. The wing-body computational results are compared with experimental data for both trailing-edge flaps deflected. The experimental flow conditions are M subinfinity = 0.4, a turbulent Reynolds number of 5.41 million based on a mean aerodynamic chord of 25.959 inches, adiabatic wall, and angles-of-attack varying from 0 deg to 23.85 deg. The computational results are presented for the 23.85 deg angle-of-attack case. The effects of the base flow due to a model sting, the varying second and fourth order numerical dissipation, and the turbulence model are all considered.

Chaussee, Denny S.↗

Application of a Third Order Upwind Scheme to Viscous Flow over Clean and Iced Wings

A 3-D compressible Navier-Stokes solver has been developed and applied to 3-D viscous flow over clean and iced wings. This method uses a third order accurate finite volume scheme with flux difference splitting to model the inviscid fluxes, and second order accurate symmetric differences to model the viscous terms. The effects of turbulence are modeled using a Kappa-epsilon model. In the vicinity of the sold walls the kappa and epsilon values are modeled using Gorski's algebraic model. Sampling results are presented for surface pressure distributions, for untapered swept clean and iced wings made of NACA 0012 airfoil sections. The leading edge of these sections is modified using a simulated ice shape. Comparisons with experimental data are given.

Bangalore, A.↗

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗

More About Generating Three-Dimensional Grids About Anything

Report presents additional details of mathematical basis and potential applications of method of generating three-dimensional grids fitted to bodies having arbitrary, complicated shapes. Grids used in finite-difference or finite-volume computations of flows around bodies. Grid-generation method involves numerical solutions of Poisson's equations described in "Generating Three-Dimensional Grids About Anything" (ARC-12620).

Sorenson, Reese L.↗

Aspects of unstructured grids and finite-volume solvers for the Euler and Navier-Stokes equations

Basic algorithms for unstructured mesh generation and fluid flow calculation are discussed. In particular the following are addressed: preliminaries of graphs and meshes; duality and data structures; basic graph operations important in CFD (Computational Fluid Dynamics); triangulation methods, including Varonoi diagrams and Delaunay triangulation; maximum principle analysis; finite volume schemes for scalar conservation law equations; finite volume schemes for the Euler and Navier-Stokes equations; and convergence acceleration for steady state calculations.

Barth, T. J.↗

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗

Numerical Analysis of Cavitation-Induced Two-Phase Fluid Transients in Cryogenic Feedlines Using Finite Volume Procedure

This article presents a numerical investigation of the cavitation induced fluid transient in cryogenic feedline due to sudden closure of the valve A one dimensional numerical model based on finite volume procedure consisting of nodes and branches using Generalized Fluid System Simulation Program ( has been employed for the prediction of pressure variations The unsteady friction formulation accurately determines wall shear stress, and a homogeneous mixture model is used for two phase flow modeling Numerical results obtained with the present formulation are compared with experimental results and Method of Characteristic ( based results It is observed that low residual in flow variable calculation improves the solution providing numerically stable results GFSSP is able to reproduce the fluid transient peak due to vapor collapse Subsequently, it can retain the shape of pressure wave for two phase fluid transients, but it is dampening faster compared to experiments Further investigation is required with set of experiments to validate the present numerical approach for different experimental data.

Fluid Transient↗

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗