SEARCH · Engineering Papers

Results for “Finite Volume”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Edge-Based Viscous Method for Mixed-Element Node-Centered Finite-Volume Solvers

A novel, efficient, edge-based viscous (EBV) discretization method has been recently developed, implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver, and applied to viscous-kernel computations that include evaluations of meanflow viscous fluxes, turbulence-model and chemistry-model diffusion terms, and the corresponding Jacobian contributions. Initially, the EBV method had been implemented for tetrahedral grids and demonstrated multifold acceleration of all viscous-kernel computations. This paper presents an extension of the EBV method for mixed-element grids. In addition to the primal edges of a given mixed-element grid, virtual edges are introduced to connect cell nodes that are not connected by a primal edge. The EBV method uses an efficient loop over all (primal and virtual) edges and features a compact discretization stencil based on the nearest neighbors. This study verifies the EBV method and assesses its efficiency on mixed-element grids by comparing the EBV solution accuracy and iterative convergence with those of well-established solutions obtained using a cell-based viscous (CBV) discretization method. The EBV solver’s memory footprint is optimized and often smaller than the memory footprint of the CBV solver. A multifold speedup is demonstrated for all viscous-kernel computations resulting in significant reduction of the time to solutions for several benchmark mixed-element-grid computations, including simulations of a flow around NASA’s juncture-flow model and a hypersonic, chemically reacting flow around a blunt body.

Edge-based viscous method↗

NASA NTRS1995

Direct Numerical Simulation of Acoustic Waves Interacting with a Shock Wave in a Quasi-1D Convergent-Divergent Nozzle Using an Unstructured Finite Volume Algorithm

Numerical simulation of a very small amplitude acoustic wave interacting with a shock wave in a quasi-1D convergent-divergent nozzle is performed using an unstructured finite volume algorithm with a piece-wise linear, least square reconstruction, Roe flux difference splitting, and second-order MacCormack time marching. First, the spatial accuracy of the algorithm is evaluated for steady flows with and without the normal shock by running the simulation with a sequence of successively finer meshes. Then the accuracy of the Roe flux difference splitting near the sonic transition point is examined for different reconstruction schemes. Finally, the unsteady numerical solutions with the acoustic perturbation are presented and compared with linear theory results.

Bui, Trong T.↗

NASA NTRSDate not supplied

A Simplified FANG Cell-Centered Finite-Volume Method and Comparison with Other Methods for Trouble-Prone Grids

We propose a simplication of the face-averaged nodal-gradient (FANG) method for a cell-centered finite-volume Euler/Navier-Stokes solver on arbitrary grids, and compare it with other gradient methods for trouble-prone grids in two dimensions. The implementation of the FANG method is simplified by adding the face-neighbor cells of the cells around a node to a least-squares gradient stencil. The resulting method is stable for both triangular and quadrilateral grids. Although it increases the residual stencil for triangular grids, it allows the solver to work seamlessly for mixed grids and greatly simplifies the implementation, especially in three dimensions. For comparison, only explicit weighted/unweighted least-squares cell-centered and nodal gradient methods are considered. These gradients are used in both inviscid and viscous schemes, and we investigate their impact on the iterative convergence of an implicit defect-correction solver on difficult grids such as highly-curved-and-thin grids and highly distorted anisotropic grids. Finally, we will also consider a face-stencil-based limiter and compare it with a conventional cell-stencil-based limiter.

Hiroaki Nishikawa↗

NASA NTRS1989

A nonoscillatory, characteristically convected, finite volume scheme for multidimensional convection problems

A new, nonoscillatory upwind scheme is developed for the multidimensional convection equation. The scheme consists of an upwind, nonoscillatory interpolation of data to the surfaces of an intermediate finite volume; a characteristic convection of surface data to a midpoint time level; and a conservative time integration based on the midpoint rule. This procedure results in a convection scheme capable of resolving discontinuities neither aligned with, nor convected along, grid lines.

Yokota, Jeffrey W.↗

NASA NTRS1990

A nonoscillatory, characteristically convected, finite volume scheme for multidimensional convection problems

Yokota, Jeffrey W.↗

NASA NTRS1992

A finite-volume numerical method to calculate fluid forces and rotordynamic coefficients in seals

A numerical method to calculate rotordynamic coefficients of seals is presented. The flow in a seal is solved by using a finite-volume formulation of the full Navier-Stokes equations with appropriate turbulence models. The seal rotor is perturbed along a diameter such that the position of the rotor is a sinusoidal function of time. The resulting flow domain changes with time, and the time-dependent flow in the seal is solved using a space conserving moving grid formulation. The time-varying fluid pressure reaction forces are then linked with the rotor center displacement, velocity and acceleration to yield the rotordynamic coefficients. Results for an annular seal are presented, and compared with experimental data and other more simplified numerical methods.

Athavale, M. M.↗

NASA NTRS2011

Notes on Accuracy of Finite-Volume Discretization Schemes on Irregular Grids

Truncation-error analysis is a reliable tool in predicting convergence rates of discretization errors on regular smooth grids. However, it is often misleading in application to finite-volume discretization schemes on irregular (e.g., unstructured) grids. Convergence of truncation errors severely degrades on general irregular grids; a design-order convergence can be achieved only on grids with a certain degree of geometric regularity. Such degradation of truncation-error convergence does not necessarily imply a lower-order convergence of discretization errors. In these notes, irregular-grid computations demonstrate that the design-order discretization-error convergence can be achieved even when truncation errors exhibit a lower-order convergence or, in some cases, do not converge at all.

Diskin, Boris↗

NASA NTRS1996

A Full Multi-Grid Method for the Solution of the Cell Vertex Finite Volume Cauchy-Riemann Equations

The system of inhomogeneous Cauchy-Riemann equations defined on a square domain and subject to Dirichlet boundary conditions is considered. This problem is discretised by using the cell vertex finite volume method on quadrilateral meshes. The resulting algebraic problem is overdetermined and the solution is defined in a least squares sense. By this approach a consistent algebraic problem is obtained which differs from the original one by O(h(exp 2)) perturbations of the right-hand side. A suitable cell-based convergent smoothing iteration is presented which is naturally linked to the least squares formulation. Hence, a standard multi-grid algorithm is reported which combines the given smoother and cell-based transfer operators. Some remarkable reduction properties of these operators are shown. A full multi-grid method is discussed which solves the discrete problem to the level of truncation error by employing one multi-grid cycle at each current level of discretisation. Experiments and applications of the full multi-grid scheme are presented.

Borzi, A.↗

NASA NTRSDate not supplied

Towards Robust and Accurate Implicit Gradient Methods for Second- and Third-Order Nodal-Gradient Cell-Centered Finite-Volume Discretizations on Tetrahedral Grids

In this paper, we introduce implicit gradient methods as alternatives to conventional least-squares gradient methods for second- and third-order nodal-gradient cell-centered finite-volume discretizations, where solutions are stored at cells but gradients are stored at nodes. Because of the unique configuration of solutions and gradients, implicit gradient systems developed for the node-centered edge-based discretization method can be directly applied once the numerical solutions are interpolated from cells to nodes with sufficient accuracy. The resulting defect-correction solver can be loosely coupled with a flow-equation solver, and at convergence, solutions and gradients that satisfy the corresponding residual equations are obtained. Each iteration is relatively cheap compared with least-squares methods involving hundreds of neighbors. Numerical results are presented for accuracy verification studies and some simple but realistic flow problems.

Computational Fluid Dynamics↗

NASA NTRS1987

A finite volume Euler calculation of the aerodynamics of transonic airfoil-vortex interaction

Unsteady inviscid transonic airfoil-vortex interaction is numerically analyzed by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables is used to form the numerical dissipation terms to stabilize the scheme. A velocity decomposition technique is applied to alleviate the problem of vortex diffusion by the numerical dissipation terms and to treat the interaction of a Rankine vortex with an airfoil accurately. Results obtained are compared with available numerical data.

Damodaran, Murali↗

NASA NTRS1995

Simulating the interaction of the heliosphere with the local interstellar medium: MHD results from a finite volume approach, first bidimensional results

We have designed a full compressible MHD code working on unstructured meshes in order to be able to compute accurately sharp structures embedded in large scale simulations. The code is based on a finite volume method making use of a kinetic flux splitting. A bidimensional version of the code has been used to simulate the interaction of a moving interstellar medium, magnetized or unmagnetized with a rotating and magnetized heliopspheric plasma source. Being aware that these computations are not realistic due to the restriction to two dimensions, we present it to demonstrate the ability of this new code to handle this problem. An axisymetric version, now under development, will be operational in a few months. Ultimately we plan to run a full 3d version.

Chanteur, G.↗

NASA NTRS1992

A finite-volume Euler solver for computing rotary-wing aerodynamics on unstructured meshes

An unstructured-grid solver for the unsteady Euler equations has been developed for predicting the aerodynamics of helicopter rotor blades. This flow solver is a finite-volume scheme that computes flow quantities at the vertices of the mesh. Special treatments are used for the flux differencing and boundary conditions in order to compute rotary-wing flowfields, and these are detailed in the paper. The unstructured-grid solver permits adaptive grid refinement in order to improve the resolution of flow features such as shocks, rotor wakes and acoustic waves. These capabilities are demonstrated in the paper. Example calculations are presented for two hovering rotors. In both cases, adaptive-grid refinement is used to resolve high gradients near the rotor surface and also to capture the vortical regions in the rotor wake. The computed results show good agreement with experimental results for surface airloads and wake geometry.

Strawn, Roger C.↗

NASA NTRS1973

Time-split finite-volume method for three-dimensional blunt-body flow

An efficient numerical method for calculating plane, axisymmetric, and fully three-dimensional blunt-body flow is presented. It is a second-order-accurate, time-dependent finite-volume procedure that solves the Euler equations in integral conservation-law form. These equations are written with respect to a Cartesian coordinate system in which an embedded mesh adjusts in time to the motion of the bow shock that is automatically captured as part of the weak solution. With such an adjusting mesh, oscillations in flow properties near the shock are shown to be virtually eliminated. The scheme uses a time-splitting concept that accelerates the convergence appreciably. Comparisons are made between computed and experimental results.

Rizzi, A. W.↗

NASA NTRS1990

Micro Blowing Simulations Using a Coupled Finite-Volume Lattice-Boltzman n L ES Approach

Three dimensional large-eddy simulations (LES) of single and multiple jet-in-cross-flow (JICF) are conducted using the 19-bit Lattice Boltzmann Equation (LBE) method coupled with a conventional finite-volume (FV) scheme. In this coupled LBE-FV approach, the LBE-LES is employed to simulate the flow inside the jet nozzles while the FV-LES is used to simulate the crossflow. The key application area is the use of this technique is to study the micro blowing technique (MBT) for drag control similar to the recent experiments at NASA/GRC. It is necessary to resolve the flow inside the micro-blowing and suction holes with high resolution without being restricted by the FV time-step restriction. The coupled LBE-FV-LES approach achieves this objectives in a computationally efficient manner. A single jet in crossflow case is used for validation purpose and the results are compared with experimental data and full LBE-LES simulation. Good agreement with data is obtained. Subsequently, MBT over a flat plate with porosity of 25% is simulated using 9 jets in a compressible cross flow at a Mach number of 0.4. It is shown that MBT suppresses the near-wall vortices and reduces the skin friction by up to 50 percent. This is in good agreement with experimental data.

Menon, S.↗

NASA NTRS1988

Finite volume computation of unsteady inviscid rotational transonic flows past airfoils in rigid body motion

Unsteady inviscid transonic flow over airfoils in arbitrary rigid body motion is analyzed numerically by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables are used to form the numerical dissipation terms to stabilize the scheme. Unsteady calculations are started from converged steady-state solutions as initial conditions. Nonreflective boundary conditions are imposed on the far-field boundaries. Results are presented and, where possible, validated against available numerical and experimental data for airfoils subjected to a step change in angle of attack, airfoils oscillating and plunging in transonic flow, and airfoils immersed in a time-varying free stream.

Damodaran, Murali↗

NASA NTRS1987

Free wake analysis of helicopter rotor blades in hover using a finite volume technique

A fully compressible method for determining helicopter rotor wake effects is described which computes the wake without requiring external specification of the wake, or separate computations for the wake and blade region. The method is a modification of a compressible finite volume Potential Flow technique, and it has been implemented in a program, HELIX I, for computing compressible rotor flow fields in hover with free wakes. Wake positions in substantial agreement with experiment have been calculated for cases including subsonic and transonic flows, high and low aspect ratios, and two- and four-bladed rotors.

Ramachandran, K.↗

NASA NTRS1988

Effects of numerical dissipation on finite-volume solutions of compressible flow problems

The effects of numerical dissipation upon solutions to the Euler equations are considered, and results for transonic flows past airfoils are presented to demonstrate the effects of the dissipative terms. The equations are approximated using a finite-volume spatial approximation with added dissipation provided by an adaptive mixture of second and fourth differences. The resulting difference equations are solved using either an explicit multistage Runge-Kutta method or a diagonalized implicit method. It is found that errors in surface values can be introduced by the averaging required to calculate derived quantities of interest.

Caughey, David A.↗

NASA NTRS2001

Recent Developments in DAO's Finite-Volume Data Assimilation System

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. The data assimilation method implemented in CODAS incorporates a simplified version of the model bias estimation and correction algorithm, as described by Dee and da Silva (1998). In this talk we will briefly describe the general system formulation, and focus on the impact of 3 data types recently introduced, namely: 1) cloud tracks winds from the Multi-angle Imaging Spectrometer by the US Air Force, and 3) temperature and moisture information derived from GPS refractivity occultation measurements. The impact of these data types on observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills will be discussed. In addition we will assess the impact of cloud assimilation on top of the atmosphere radiation fields estimated from CERES measurements.

daSilva, Arlindo↗