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 235 records · Page 13

Some implementational issues of convection schemes for finite volume formulations

Two higher-order upwind schemes - second-order upwind and QUICK - are examined in terms of their interpretation, implementation as well as performance for a recirculating flow in a lid-driven cavity, in the context of a control volume formulation using the SIMPLE algorithm. The present formulation of these schemes is based on a unified framework wherein the first-order upwind scheme is chosen as the basis, with the remaining terms being assigned to the source term. The performance of these schemes is contrasted with the first-order upwind and second-order central difference schemes. Also addressed in this study is the issue of boundary treatment associated with these higher-order upwind schemes. Two different boundary treatments - one that uses a two-point scheme consistently within a given control volume at the boundary, and the other that maintains consistency of flux across the interior face between the adjacent control volumes - are formulated and evaluated.

Thakur, Siddharth↗

Some implementational issues of convection schemes for finite-volume formulations

Two higher-order upwind schemes - second-order upwind and QUICK - are examined in terms of their interpretation, implementations, as well as performance for a recirculating flow in a lid-driven cavity, in the context of a control-volume formulation using the SIMPLE algorithm. The present formulation of these schemes is based on a unified framework wherein the first-order upwind scheme is chosen as the basis, with the remaining terms being assigned to the source term. The performance of these schemes is contrasted with the first-order upwind and second-order central difference schemes. Also addressed in this study is the issue of boundary treatment associated with these higher-order upwind schemes. Two different boundary treatments - one that uses a two-point scheme consistently within a given control volume at the boundary, and the other that maintains consistency of flux across the interior face between the adjacent control volumes - are formulated and evaluated.

Thakur, Siddharth↗

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. These algorithms are now used in the practical engineering design of devices such as cars and airplanes as well as theoretical studies of complex phenomena such as fluid turbulence. In past years, limited computer resources have motivated the development of efficient numerical methods in computational fluid dynamics (CFD) utilizing structured meshes. These meshes are comprised of systematic arrays of quadrilateral or hexahedral cells. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Structured meshes also permit the use of highly efficient solution techniques such as alternating direction implicit (ADI) iteration schemes or multigrid. Following the dramatic improvement in computing speed in recent years, emphasis has shifted towards the design of algorithms capable of treating complex geometries. The automatic generation of structured grids about complex geometries is problematic. Unstructured grids offer one promising alternative technique for treating these general geometries. Unstructured meshes have irregular connectivity and usually contain triangles and/or quadrilaterals in two dimensions and tetrahedra and/or hexahedra in three dimensions. The generation and use of unstructured grids poses new challenges in computational fluid dynamics. This is true for both grid generation as well as for the design of algorithms for flow solution. The purpose of these notes is to present recent developments in the unstructured grid generation and flow solution technology.

Barth, T. J.↗

Coupled-Circulation-Chemistry Studies with the Finite-Volume CCM: Trace Gas Transport in the Tropopause Region

A joint project between the Data Assimilation Office at NASA GSFC and NCAR involves linking the physical packages from the Community Climate Model (CCM) with the flux-form semi-Lagrangian dynamical core developed by Lin and Rood in the DAO. A further development of this model includes the implementation of a chemical package developed by Douglass and colleagues in the Atmospheric Chemistry and Dynamics Branch at NASA GSFC. Results from this coupled dynamics-radiation-chemistry model will be presented, focussing on trace gas transport in the tropopause region.

Pawson, Steven↗

Climatology of the Finite-Volume CCM and Sensitivity Studies

A new model has been developed in the Data Assimilation Office at NASA Goddard Space Flight Center. It is intended for data assimilation and general circulation modeling, including representations of radiatively active trace gases such as stratospheric ozone. The model is based on the flux-form semi-Lagrangian' dynamical core, developed by Lin and Rood. In collaboration with the National Center for Atmospheric Research, the physical parametrization package from the Community Climate model, Version 3, has been included. This paper will discuss the climatology of a multi-year integration of the model, with observed boundary conditions (sea-surface temperature, sea ice, and incoming solar radiation). The standard model version has a resolution of 2.5 degrees longitude, 2 degrees latitude and has 55 levels between the surface and 80km. The mean climate and interannual variability of this model will be discussed, with emphasis on the middle atmosphere. The sensitivity to horizontal resolution will be examined, as will the impact of increasing and decreasing the vertical resolution. Changes which differ significantly from the interannual variations of the control run will be highlighted. Impacts of using long-term mean boundary conditions or values which vary from year to year will be discussed. Changes in the tropospheric circulation induced by lowering the upper boundary of the model will be discussed. The emphasis of the paper is on the climatology of the middle atmosphere and its impacts on the troposphere, with particular attention given to the climatologically important region near the tropopause.

Pawson, Steven↗

Spectral (Finite) Volume Method for One Dimensional Euler Equations

Consider a mesh of unstructured triangular cells. Each cell is called a Spectral Volume (SV), denoted by Si, which is further partitioned into subcells named Control Volumes (CVs), indicated by C(sub i,j). To represent the solution as a polynomial of degree m in two dimensions (2D) we need N = (m+1)(m+2)/2 pieces of independent information, or degrees of freedom (DOFs). The DOFs in a SV method are the volume-averaged mean variables at the N CVs. For example, to build a quadratic reconstruction in 2D, we need at least (2+1)(3+1)/2 = 6 DOFs. There are numerous ways of partitioning a SV, and not every partition is admissible in the sense that the partition may not be capable of producing a degree m polynomial. Once N mean solutions in the CVs of a SV are given, a unique polynomial reconstruction can be obtained.

Wang, Z. J.↗

Implementation of the NCAR Community Land Model (CLM) in the NASA/NCAR finite-volume Global Climate Model (fvGCM)

In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).

Radakovich, Jon D.↗

Accuracy Analysis for Finite-Volume Discretization Schemes on Irregular Grids

A new computational analysis tool, downscaling test, is introduced and applied for studying the convergence rates of truncation and discretization errors of nite-volume discretization schemes on general irregular (e.g., unstructured) grids. The study shows that the design-order convergence of discretization errors can be achieved even when truncation errors exhibit a lower-order convergence or, in some cases, do not converge at all. The downscaling test is a general, efficient, accurate, and practical tool, enabling straightforward extension of verification and validation to general unstructured grid formulations. It also allows separate analysis of the interior, boundaries, and singularities that could be useful even in structured-grid settings. There are several new findings arising from the use of the downscaling test analysis. It is shown that the discretization accuracy of a common node-centered nite-volume scheme, known to be second-order accurate for inviscid equations on triangular grids, degenerates to first order for mixed grids. Alternative node-centered schemes are presented and demonstrated to provide second and third order accuracies on general mixed grids. The local accuracy deterioration at intersections of tangency and in flow/outflow boundaries is demonstrated using the DS tests tailored to examining the local behavior of the boundary conditions. The discretization-error order reduction within inviscid stagnation regions is demonstrated. The accuracy deterioration is local, affecting mainly the velocity components, but applies to any order scheme.

Diskin, Boris↗

Improved Finite-Volume Method for Radiative Hydrodynamics

Fully coupled simulations of hydrodynamics and radiative transfer are essential to a number of fields ranging from astrophysics to engineering applications. Of particular interest in this work are hypersonic atmospheric entries and associated experimental apparatus, e.g., shock tubes and high enthalpy testing facilities. The radiative transfer calculations must supply to the CFD a heating term in the energy equation in the form of the divergence of the radiative heat flux and the radiative heat fluxes to bounding surfaces. It is most efficient to solve the radiative transfer equation on the same grid as the CFD solution, and this work presents an algorithm with improved accuracy for such simulations on structured and unstructured grids compared to more conventional approaches. Results will be shown for shock radiation during hypersonic reentry. Issues of parallelization within a radiation sweep will also be discussed.

Wray, Alan↗

Face- and Cell-Averaged Nodal-Gradient Approach to Cell-Centered Finite-Volume Method on Mixed Grids

In this paper, the averaged nodal-gradient approach previously developed for triangular grids is extended to mixed triangular-quadrilateral grids. It is shown that the face- averaged approach leads to deteriorated iterative convergence on quadrilateral grids. To develop a convergent solver, we consider cell-averaging instead of face-averaging for quadri- lateral cells. We show that the cell-averaged approach leads to a convergent solver and can be efficiently combined with the face-averaged approach on mixed grids. The method is demonstrated for various inviscid and viscous problems from low to high Mach numbers on two-dimensional mixed grids.

Nishikawa, Hiroaki↗

A 3-D Nodal-Averaged Gradient Approach for Unstructured-Grid Cell-Centered Finite-Volume Methods for Application to Turbulent Hypersonic Flow

A 2-D nodal weighted least-squares gradient method and a related face-averaged nodal gradient approach that were developed for use with triangular grids are extended to 3-D for use with tetrahedral grids. In addition, a method, developed in 2-D, to stabilize the iterative convergence of these methods on quadrilateral cells is described and extended to 3-D and remedies are investigated to determine the nodal gradient averaging approach most suitable for use with grids made up of hexahedral, prismatic, pyramidal and tetrahedral cells. Moreover, due to an interest in hypersonic flow, a robust multidimensional gradient limiter procedure that is consistent with the stencil used to construct the nodal gradients is described. Finally, we demonstrate that the resulting 3-D methods are sufficiently robust for use in scramjet computations through the solution of three canonical turbulent hypersonic flow problems as well as a physically realistic 3-D scramjet inlet geometry.

Jeffery A White↗

Finite-Volume Diffusion Schemes for Svard's Eulerian Governing Equations

This paper, contains discussion on the implementation of new diffusion schemes for Svard’s Eulerian flow (EF) governing equations. It contains an alpha-damping type scheme and a new hyperbolic Eulerian flow (HEF) discretization paralleling previous work on hyperbolic Navier-Stokes (HNS). It shows these equations have simplified hyperbolic discretizations due to the simplicity of the new stress tensor and compares these schemes to established Navier-Stokes (NS) discretizations to verify the predictive utilities of such diffusion schemes.

Computational Fluid Dynamics↗

Finite-Volume Diffusion Schemes for Svard's Eulerian Governing Equations

This paper, contains discussion on the implementation of new diffusion schemes for Svard’s Eulerian flow (EF) governing equations. It contains an alpha-damping type scheme and a new hyperbolic Eulerian flow (HEF) discretization paralleling previous work on hyperbolic Navier-Stokes (HNS). It shows these equations have simplified hyperbolic discretizations due to the simplicity of the new stress tensor and compares these schemes to established Navier-Stokes (NS) discretizations to verify the predictive utilities of such diffusion schemes.

Computational Fluid Dynamics↗