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 37 records · Page 2

Accelerated finite-volume calculation of transonic potential flows

A fully conservative finite-volume algorithm is used to calculate transonic potential flows past isolated airfoils and through two-dimensional channels. The difference equations are solved by a multigrid technique which uses an alternating direction implicit method as a smoothing algorithm. The finite-volume formulation provides a framework within which to treat flows past complicated geometries, while the multigrid/alternating direction scheme provides rapid convergence of the solution to very small residuals.

Jameson, A.↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume numerical techniques with finite thickness models for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the one-dimensional semi -infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody were investigated. An array of streamwise orientated heating striations were generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients due to the striation patterns two-dimensional heat transfer techniques were necessary to obtain accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations--potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomeclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations: potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomenclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗

Split Space-Marching Finite-Volume Method for Chemically Reacting Supersonic Flow

A space-marching finite-volume method employing a nonorthogonal coordinate system and using a split differencing scheme for calculating steady supersonic flow over aerodynamic shapes is presented. It is a second-order-accurate mixed explicit-implicit procedure that solves the inviscid adiabatic and nondiffusive equations for chemically reacting flow in integral conservation-law form. The relationship between the finite-volume and differential forms of the equations is examined and the relative merits of each discussed. The method admits initial Cauchy data situated on any arbitrary surface and integrates them forward along a general curvilinear coordinate, distorting and deforming the surface as it advances. The chemical kinetics term is split from the convective terms which are themselves dimensionally split, thereby freeing the fluid operators from the restricted step size imposed by the chemical reactions and increasing the computational efficiency. The accuracy of this splitting technique is analyzed, a sufficient stability criterion is established, a representative flow computation is discussed, and some comparisons are made with another method.

Rizzi, Arthur W.↗

Finite volume solultion of the compressible boundary layer equations

A box-type finite volume discretization is applied to the integral form of the compressible boundary layer equations. Boundary layer scaling is introduced through the grid construction: streamwise grid lines follow eta = y/h = const., where y is the normal coordinate and h(x) is a scale factor proportional to the boundary layer thickness. With this grid, similarity can be applied explicitly to calculate initial conditions. The finite volume method preserves the physical transparency of the integral equations in the discrete approximation. The resulting scheme is accurate, efficient, and very simple. Computations for similar and non-similar flows show excellent agreement with tabulated results, solutions computed with Keller's Box scheme, and experimental data.

Loyd, B.↗

A vectorized, finite-volume, adaptive grid algorithm applied to planetary entry problems

An adaptive grid, finite-volume method has been applied to problems in planetary entry for computing complete flowfields. The adaption algorithm is implicit in nature and is keyed to resolve user specified gradients. The finite-volume algorithm is explicit, utilizing a maximum time step advancement at each grid point to accelerate convergence to the steady state. The present version of the code is for the laminar flow of a perfect gas. The role of the adaption algorithm in resolving various features of blunt body/wake flow for planetary entry conditions is emphasized.

Gnoffo, P. A.↗

A time-split finite-volume algorithm for three-dimensional flow-field simulation

A general finite-volume algorithm is developed for solving three-dimensional, time-dependent, compressible Navier-Stokes equations for high Reynolds number flows over an arbitrary geometry. This algorithm adapts MacCormack's (1982) explicit-implicit scheme to a time-split, three-dimensional finite-volume concept in a general coordinate system. It is shown that the thin-layer approximation in all three spatial directions significantly reduces the evaluation of viscous terms and allows the algorithm to solve more complicated geometries with all boundaries in two or all three directions. The calculated results using this method are found to be in good agreement with the experimental measurements of a blunt-fin induced shock wave and boundary-layer interaction problems. Observations of the existence of peak pressure, primary horseshoe and secondary vortices, and reversed supersonic zones show that computational fluid dynamics can effectively supplement the wind tunnel tests for aerodynamic design as well as for understanding basic fluid dynamics.

Hung, C. M.↗

Finite volume solution of the compressible boundary-layer equations

A box-type finite volume discretization is applied to the integral form of the compressible boundary layer equations. Boundary layer scaling is introduced through the grid construction: streamwise grid lines follow eta = y/h = const., where y is the normal coordinate and h(x) is a scale factor proportional to the boundary layer thickness. With this grid, similarity can be applied explicity to calculate initial conditions. The finite volume method preserves the physical transparency of the integral equations in the discrete approximation. The resulting scheme is accurate, efficient, and conceptually simple. Computations for similar and non-similar flows show excellent agreement with tabulated results, solutions computed with Keller's Box scheme, and experimental data.

Loyd, B.↗

Development of a finite volume time marching method

Progress made in the development of a time-marching finite volume method from September 1984 to December 1984 is reported. The objective of the work is to develop and demonstrate a Navier-Stokes approach for transonic flow which includes viscous terms in the finite volume method. The accuracy of the computational method will be verified using a transonic diffuser as a test case. The computational goal is to calculate the flow in sufficient detail and with sufficient accuracy that the loss generating mechanisms can be studied to assess the sources of inefficiency in the transonic diffuser.

Nicholson, Stephen↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume heat conduction techniques for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the standard one-dimensional semi-infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the NASA Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody was investigated. An array of streamwise-orientated heating striations was generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients caused by striation patterns multi-dimensional heat transfer techniques were necessary to obtain more accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates compared to 2-D analysis because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

Application of Local Discretization Methods in the NASA Finite-Volume General Circulation Model

We present the basic ideas of the dynamics system of the finite-volume General Circulation Model developed at NASA Goddard Space Flight Center for climate simulations and other applications in meteorology. The dynamics of this model is designed with emphases on conservative and monotonic transport, where the property of Lagrangian conservation is used to maintain the physical consistency of the computational fluid for long-term simulations. As the model benefits from the noise-free solutions of monotonic finite-volume transport schemes, the property of Lagrangian conservation also partly compensates the accuracy of transport for the diffusion effects due to the treatment of monotonicity. By faithfully maintaining the fundamental laws of physics during the computation, this model is able to achieve sufficient accuracy for the global consistency of climate processes. Because the computing algorithms are based on local memory, this model has the advantage of efficiency in parallel computation with distributed memory. Further research is yet desirable to reduce the diffusion effects of monotonic transport for better accuracy, and to mitigate the limitation due to fast-moving gravity waves for better efficiency.

Yeh, Kao-San↗

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit scheme can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit schemes can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗

Viscous analysis of high speed flows using an upwind finite volume technique

A numerical study using an upwind finite volume technique to analyze high speed viscous flows is described. Two test cases were considered: a Mach 4.0 shock wave/laminar-boundary-layer interaction, and Mach 14.1 compression corners (15-deg and 24-deg ramp angle). PNS, TLNS, and NS equation results were obtained on each test case. As expected, PNS solutions are not adequate to resolve the flow physics. For the grids used in this study, there is no noticeable difference between the TLNS and full NS solutions. Both the TLNS and full NS solutions compared favorably with experimental data.

Ng, W. F.↗

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

Lin, Shian-Jiann↗