SEARCH · Engineering Papers

Results for “finite volume discretization”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1999

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

Barth, Timothy↗

NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗

NASA NTRS1984

Dynamic characteristics of a hydrostatic gas bearing driven by oscillating exhaust pressure

Vibration of a statically loaded, inherently compensated hydrostatic journal bearing due to oscillating exhaust pressure is investigated. Both angular and radial vibration modes are analyzed. The time-dependent Reynolds equation governing the pressure distribution between the oscillating journal and sleeve is solved together with the journal equation of motion to obtain the response characteristics of the bearing. The Reynolds equation and the equation of motion are simplified by applying regular perturbation theory for small displacements. The numerical solutions of the perturbation equations are obtained by discretizing the pressure field using finite-difference aproximations with a discrete, nonuniform line-source model which excludes effects due to feeding hole volume. An iterative scheme is used to simultaneously satisfy the equations of motion for the journal. The results presented include Bode plots of bearing-oscillation gain and phase for a particular bearing configuration for various combinations of parameters over a range of frequencies, including the resonant frequency.

Watkins, C. B.↗

NASA NTRS2018

Unstructured Grid Adaptation and Solver Technology for Turbulent Flows

Unstructured grid adaptation is a tool to control Computational Fluid Dynamics (CFD) discretization error. However, adaptive grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Issues that prevent the use of adaptive grid methods are identified by applying unstructured grid adaptation methods to a series of benchmark cases. Once identified, these challenges to existing adaptive workflows can be addressed. Unstructured grid adaptation is evaluated for test cases described on the Turbulence Modeling Resource (TMR) web site, which documents uniform grid refinement of multiple schemes. The cases are turbulent flow over a Hemisphere Cylinder and an ONERA M6Wing. Adaptive grid force and moment trajectories are shown for three integrated grid adaptation processes with Mach interpolation control and output error based metrics. The integrated grid adaptation process with a finite element (FE) discretization produced results consistent with uniform grid refinement of fixed grids. The integrated grid adaptation processes with finite volume schemes were slower to converge to the reference solution than the FE method. Metric conformity is documented on grid/metric snapshots for five grid adaptation mechanics implementations. These tools produce anisotropic boundary conforming grids requested by the adaptation process.

Park, Michael A.↗

NASA NTRS2019

Numerical Modeling of Boiling in a Heated Tube

This paper presents numerical models of boiling in a heated tube using the Generalized Fluid System Simulation Program (GFSSP), a finite-volume-based general-purpose flow network code developed at NASA/Marshall Space Flight Center. The heated tube is discretized into a one-dimensional array of nodes and branches to represent the flow of liquid and vapor in a tube with a prescribed pressure differential. The solid wall is also discretized into solid nodes and conductors to allow for heat transfer between the wall and the fluid. The conservation equations of mass, momentum, and energy of the fluid are solved simultaneously with the energy conservation equation for the solid wall. Two experimental configurations of fluid flowing in a vertical tube have been simulated, one with water and the other with liquid hydrogen. This paper compares experimental data with numerical predictions based on four different published correlations for boiling heat transfer coefficients. Three of these correlations are applicable to the saturated vertical flow conditions of the experiments. One of them is applicable to film boiling and has been used for the liquid hydrogen experiment, which was in film boiling regime. For the case of boiling water, the predictions of wall temperatures using the boiling heat transfer correlations agreed well with the experimental results. However, in the case of boiling hydrogen larger discrepancies were observed between the experimental data and numerical predictions.

Majumdar, Alok↗

NASA NTRS2019

Accuracy, Scalability, and Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities including parallelization, line generation for general unstructured grids, improved discretization scheme, and optimized iterative solver. The paper reports on the new developments to the flow solver and assesses the accuracy, scalability, and efficiency. The USM3D assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Two benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere-cylinder configuration and a transonic flow around the ONERA M6 wing are considered.

Pandya, Mohagna J.↗

NASA NTRS1991

Development of iterative techniques for the solution of unsteady compressible viscous flows

Efficient iterative solution methods are being developed for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. Thus, the extra work required by iterative schemes can also be designed to perform efficiently on current and future generation scalable, missively parallel machines. An obvious candidate for iteratively solving the system of coupled nonlinear algebraic equations arising in CFD applications is the Newton method. Newton's method was implemented in existing finite difference and finite volume methods. Depending on the complexity of the problem, the number of Newton iterations needed per step to solve the discretized system of equations can, however, vary dramatically from a few to several hundred. Another popular approach based on the classical conjugate gradient method, known as the GMRES (Generalized Minimum Residual) algorithm is investigated. The GMRES algorithm was used in the past by a number of researchers for solving steady viscous and inviscid flow problems with considerable success. Here, the suitability of this algorithm is investigated for solving the system of nonlinear equations that arise in unsteady Navier-Stokes solvers at each time step. Unlike the Newton method which attempts to drive the error in the solution at each and every node down to zero, the GMRES algorithm only seeks to minimize the L2 norm of the error. In the GMRES algorithm the changes in the flow properties from one time step to the next are assumed to be the sum of a set of orthogonal vectors. By choosing the number of vectors to a reasonably small value N (between 5 and 20) the work required for advancing the solution from one time step to the next may be kept to (N+1) times that of a noniterative scheme. Many of the operations required by the GMRES algorithm such as matrix-vector multiplies, matrix additions and subtractions can all be vectorized and parallelized efficiently.

Sankar, Lakshmi N.↗

NASA NTRS1991

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗

NASA NTRS1992

Numerical simulation of conservation laws

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

Chang, Sin-Chung↗

NASA NTRS1988

Computations of vortical supersonic flows over a blunt-nose-cylinder at high angles of attack

Three-dimensional, viscous, and separated flows over a blunt-nose-cylinder at 20, 32, and 44-deg attack angles were computed. The approaching freestream was at a Mach number of 1.6 and a unit Reynolds number of 2 x 10 to the 6th/ft with a total temperature of 585 R. The cylinder used for the computations had a length-to-diameter ratio of 6.67 with a base diameter of 3 in. The flowfield was dominated by large-scale and multiple vortices generated by crossflow separation. The effect of turbulence on the flow structure of one case was modeled algebraically with modifications to correct the length and velocity scales in the regions of separation. The mass averaged Navier-Stokes equations were solved by an approximately factored, upwind-biased, implicit, finite volume scheme. The initialization of the flows was enhanced by a mesh sequencing strategy applied to the diagonalized form of the discretized equations. The convergence to steady-state was accelerated by a multigrid algorithm and using the block inversions for the discretized equations. Calculations were compared with experimental results.

Baysal, O.↗

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, Dimitri J.↗

NASA NTRS1990

Multigrid solution of the Navier-Stokes equations on triangular meshes

Mavriplis, D. J.↗

NASA NTRS1997

Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow

It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.

Barth, Timothy J.↗

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A new Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite-element Galerkin approximation, which can be shown to be equivalent to a finite-volume approximation for regular equilateral triangular meshes. Integration to steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. The convergence efficiency of the present method is shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, D. J.↗

NASA NTRS1989

An investigation of cell centered and cell vertex multigrid schemes for the Navier-Stokes equations

Two efficient and robust finite-volume multigrid schemes for solving the Navier-Stokes equations are investigated. These schemes employ either a cell centered or a cell vertex discretization technique. An explicit Runge-Kutta algorithm is used to advance the solution in time. Acceleration techniques are applied to obtain faster steady-state convergence. Accuracy and convergence of the schemes are examined. Computational results for transonic airfoil flows are essentially the same, even for a coarse mesh. Both schemes exhibit good convergence rates for a broad range of artificial dissipation coefficients.

Radespiel, R.↗

NASA NTRS2003

Particle Size Distributions in Atmospheric Clouds

In this note, we derive a transport equation for a spatially integrated distribution function of particles size that is suitable for sparse particle systems, such as in atmospheric clouds. This is done by integrating a Boltzmann equation for a (local) distribution function over an arbitrary but finite volume. A methodology for evolving the moments of the integrated distribution is presented. These moments can be either tracked for a finite number of discrete populations ('clusters') or treated as continuum variables.

Paoli, Roberto↗

NASA NTRS2014

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity↗

NASA NTRS1992

Thermal residual stresses in silicon-carbide/titanium (0/90) laminate

The current work formulated a micromechanical analysis of a cross-ply laminate and calculated the thermal residual stress in a very thick (0/90)(sub 2n) silicon-carbide/titanium laminate. Results were also shown for a unidirectional laminate of the same material. Discrete fiber-matrix models assuming a rectangular array of fibers with a fiber volume fraction of 32.5 percent and a three-dimensional, finite-element analysis were used. Significant differences in the trends and magnitudes for the fiber, matrix, and interface stresses were calculated for unidirectional and (0/90) models. Larger hoop stresses calculated for the (0/90) model indicate that it may be more susceptible to radial cracking when subjected to mechanical loading than the unidirectional model. The axial stresses in the matrix were calculated to be slightly larger for the (0/90) model. The compressive axial stresses in the fiber were significantly larger in the (0/90) model. The presence of the cross-ply in the (0/90) model reduced the constraint on the fiber, producing radial interface stresses that were less compressive, which could lead to earlier failure of the fiber-matrix interface.

Bigelow, C. A.↗