SEARCH · Engineering Papers

Results for “Finite Volume Scheme”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Three-dimensional unsteady Euler equations solution using flux vector splitting

A method for numerically solving the three-dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of 2 and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far-field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three-dimensional Euler code and experimental data.

Whitfield, D. L.↗

NASA NTRS1992

Combined forced and free convection in a curved duct

The purpose of this study is to investigate the flow and heat transfer characteristics of a combined forced and free convection flow in a curved duct. Solutions are obtained by solving the low Mach number model of the Navier-Stokes equation using a control volume method. The finite-volume method was developed with the use of a predictor-corrector numerical scheme and some new variations of the classical projection method. Solutions indicated that the existence of a buoyancy force has changed the entire flow structure inside a curved duct. Reversed flow at both inner and outer bend is observed. For moderate Reynolds number, the upstream section of the duct was significantly influenced by the free convection processes. In general, heat transfer is strong at the inner bend of the beginning of the heated section and at the outer bend on the last half of the heated section. The maximum velocity location is strongly influenced by the combined effects of buoyancy and centrifugal forces. A strong buoyancy force can reduce the strength of the secondary flow where it plays an important role in mixing.

Yam, Clement G.↗

NASA NTRS1993

Three-dimensional unstructured grid Euler method applied to turbine blades

Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.

Kwon, Oh J.↗

NASA NTRS1993

A comparison of two incompressible Navier-Stokes algorithms for unsteady internal flow

A comparative study of two different incompressible Navier-Stokes algorithms for solving an unsteady, incompressible, internal flow problem is performed. The first algorithm uses an artificial compressibility method coupled with upwind differencing and a line relaxation scheme. The second algorithm uses a fractional step method with a staggered grid, finite volume approach. Unsteady, viscous, incompressible, internal flow through a channel with a constriction is computed using the first algorithm. A grid resolution study and parameter studies on the artificial compressibility coefficient and the maximum allowable residual of the continuity equation are performed. The periodicity of the solution is examined and several periodic data sets are generated using the first algorithm. These computational results are compared with previously published results computed using the second algorithm and experimental data.

Wiltberger, N. Lyn↗

NASA NTRS2005

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing During Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on simulations using the NASA Global Modeling and Assimilation Office (GMAO) finite-volume general circulation model (fvGCM) with Microphyics of clouds in Relaxed Arakawa Schubert Scheme (McRAS). The aerosol loading are prescribed from three-dimensional monthly distribution of tropospheric aerosols viz., sulfate, black carbon, organic carbon, soil dust, and sea salt from output of the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The aerosol extinction coefficient, single scattering albedo, and asymmetric factor are computed as wavelength-dependent radiative forcing in the radiative transfer scheme of the fvGCM, and as a function of the aerosol loading and ambient relative humidity. We find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excites a planetary scale teleconnection pattern in sea level pressure, temperature and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, i.e., South Asia, East Asia, and northern and western Africa. Additionally, atmospheric heating is found in regions with large loading of dust (over Northern Africa, and Middle East), and black carbon (over South-East Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east-west dipole anomaly with strong cooling over the Caspian Sea, and warming over central and northeastern Asia, where aerosol concentration are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations.

Kim, Maeng-Ki↗

NASA NTRS1986

An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry

High-order (p-type) finite element methods combine the geometric flexibility of standard low-order finite element (or finite volume) techniques with the rapid convergence properties of spectral methods. Various p-type schemes have been proposed for elliptic problems. However, an application of p-type methods to more complex equations, in particular the passive scalar or Navier-Stokes equations, is complicated by the nature of the equations and by efficiency considerations as regards the solution procedure. Patera (1984) has proposed a p-type method for the Navier-Stokes equations, taking into account the arising additional complications. The present study is concerned with an isoparametric spectral element discretization for the solution of flow problems in fully general, curvy geometries. Attention is given to the isoparametric formulation for elliptic equations and a time-splitting scheme for the Navier-Stokes equations.

Korczak, K. Z.↗

NASA NTRS1987

Unsteady Euler solution of transonic helicopter rotor flow

A new time-accurate Euler method is presented for calculating the transonic flow over a lifting or nonlifting rotor blade in both hover and forward flight. The approach is to solve the conservative Euler equations in a rotor-fixed frame of reference using a finite-volume method. The discretized equations are solved by using a Runge-Kutta multistage scheme with a new higher-order implicit residual smoothing procedure. Preliminary results were compared with wind-tunnel data. In all the cases considered, good agreement was found with experimental data.

Chang, I-Chung↗

NASA NTRS1989

Numerical modeling of turbulent supersonic reacting coaxial jets

The paper considers the mixing and subsequent combustion within turbulent reacting shear layers. A computer program was developed to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method integrates the Reynolds averaged Navier-Stokes equations using a finite volume approach while advancing the solution forward in time using a Runge-Kutta scheme. Three separate flowfields are investigated and it is found that no single turbulence model considered could accurately predict the degree of mixing for all three cases.

Eklund, Dean R.↗

NASA NTRS2006

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing during Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite a planetary-scale teleconnection pattern in sea level pressure, temperature, and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east west dipole anomaly with strong cooling over the Caspian Sea and warming over central and northeastern Asia, where aerosol concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Kim, Maeng-Ki↗

NASA NTRS1989

A Simple Algebraic Grid Adaptation Scheme with Applications to Two- and Three-dimensional Flow Problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented. The scheme locally adjusts the grid density based on gradients of selected flow variables from either finite difference or finite volume calculations. A user-prescribed grid stretching can be specified such that control of the grid spacing can be maintained in areas of known flowfield behavior. For example, the grid can be clustered near a wall for boundary layer resolution and made coarse near the outer boundary of an external flow. A grid smoothing technique is incorporated into the adaptive grid routine, which is found to be more robust and efficient than the weight function filtering technique employed by other researchers. Since the present algebraic scheme requires no iteration or solution of differential equations, the computer time needed for grid adaptation is trivial, making the scheme useful for three-dimensional flow problems. Applications to two- and three-dimensional flow problems show that a considerable improvement in flowfield resolution can be achieved by using the proposed adaptive grid scheme. Although the scheme was developed with steady flow in mind, it is a good candidate for unsteady flow computations because of its efficiency.

Hsu, Andrew T.↗

NASA NTRS2004

Agglomeration Multigrid for an Unstructured-Grid Flow Solver

An agglomeration multigrid scheme has been implemented into the sequential version of the NASA code USM3Dns, tetrahedral cell-centered finite volume Euler/Navier-Stokes flow solver. Efficiency and robustness of the multigrid-enhanced flow solver have been assessed for three configurations assuming an inviscid flow and one configuration assuming a viscous fully turbulent flow. The inviscid studies include a transonic flow over the ONERA M6 wing and a generic business jet with flow-through nacelles and a low subsonic flow over a high-lift trapezoidal wing. The viscous case includes a fully turbulent flow over the RAE 2822 rectangular wing. The multigrid solutions converged with 12%-33% of the Central Processing Unit (CPU) time required by the solutions obtained without multigrid. For all of the inviscid cases, multigrid in conjunction with an explicit time-stepping scheme performed the best with regard to the run time memory and CPU time requirements. However, for the viscous case multigrid had to be used with an implicit backward Euler time-stepping scheme that increased the run time memory requirement by 22% as compared to the run made without multigrid.

Frink, Neal↗

NASA NTRS1982

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗

NASA NTRS1988

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗

NASA NTRS1992

CFD Simulations Of Sonic Booms In Near And Mid Fields

Report discusses computational fluid dynamics (CFD) to simulate generation and propagation of sonic booms in near- and mid-field regions of supersonic flows about simplified bodies representative of advanced airplanes. Parabolized Navier-Stokes equations integrated by implicit, approximate-factorization, finite-volume algorithm in which crossflow inviscid fluxes evaluated by Roe's flux-difference-splitting scheme. Near-field solutions obtained by applying algorithm to flows immediately surrounding bodies. Solutions transferred to computer codes based on Whitham"s F-function theory for extrapolation to far-field.

Cheung, Samson H.↗

NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗

NASA NTRS1986

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.↗

NASA NTRS1986

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.↗

NASA NTRS1989

A simple algebraic grid adaptation scheme with applications to two- and three-dimensional flow problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented which provides high flowfield resolution when applied to two-dimensional and three-dimensional flow problems. The present scheme locally adjusts the grid density using gradients of selected flow variables from either finite difference or finite volume calculations, and it can specify user-prescribed grid stretching such that control of the grid spacing can be maintained in areas of known flowfield behavior. A robust and efficient grid smoothing technique has been incorporated into the adaptive grid routine.

Hsu, Andrew T.↗