In this article the principles of the field operation and manipulation (FOAM) C++ class library for continuum mechanics are outlined. Our intention is to make it as easy as possible to develop reliable and efficient computational continuum-mechanics codes: this is achieved by making the top-level syntax of the code as close as possible to conventional mathematical notation for tensors and partial differential equations. Object-orientation techniques enable the creation of data types that closely mimic those of continuum mechanics, and the operator overloading possible in C++ allows normal mathematical symbols to be used for the basic operations. As an example, the implementation of various types of turbulence modeling in a FOAM computational-fluid-dynamics code is discussed, and calculations performed on a standard test case, that of flow around a square prism, are presented. To demonstrate the flexibility of the FOAM library, codes for solving structures and magnetohydrodynamics are also presented with appropriate test case results given. © 1998 American Institute of Physics.

A tensorial approach to computational continuum mechanics using object-oriented techniques

AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings

The accuracy of numerical simulation algorithms is one of main concerns in modern Computational Fluid Dynamics. Development of new and more accurate mathematical models requires an insight into the problem of numerical errors. In order to construct an estimate of the solution error in Finite Volume calculations, it is first necessary to examine its sources. Discretisation errors can be divided into two groups: errors caused by the discretisation of the solution domain and equation discretisation errors. The first group includes insufficient mesh resolution, mesh skewness and non-orthogonality. In the case of the second order Finite Volume method, equation discretisation errors are represented through numerical diffusion. Numerical diffusion coefficients from the discretisation of the convection term and the temporal derivative are derived. In an attempt to reduce numerical diffusion from the convection term, a new stabilised and bounded second-order differencing scheme is proposed. Three new methods of error estimation are presented. The Direct Taylor Series Error estimate is based on the Taylor series truncation error analysis. It is set up to enable single-mesh single-run error estimation. The Moment Error estimate derives the solution error from the cell imbalance in higher moments of the solution. A suitable normalisation is used to estimate the error magnitude. The Residual Error estimate is based on the local inconsistency between face interpolation and volume integration. Extensions of the method to transient flows and the Local Residual Problem error estimate are also given. Finally, an automatic error-controlled adaptive mesh refinement algorithm is set up in order to automatically produce a solution of pre-determined accuracy. It uses mesh refinement and unrefinement to control the local error magnitude. The method is tested on several characteristic flow situations, ranging from incompressible to supersonic flows, for both steady-state and transient problems.

Error analysis and estimation for the finite volume method with applications to fluid flows

The RNG κ-e turbulence model derived by Yakhot and Orszag (1986) based on the Renormalization Group theory has been modified and applied to variable-density engine flows in the present study. The original RNG-based turbulence transport approximations were developed formally for an incompressible flow. In order to account for flow compressibility the RNG e-equation is modified and closed through an isotropic rapid distortion analysis. Computations were made of engine compressing/expanding flows and the results were compared with available experimental observations in a production diesel engine geometry. The modified RNG κ-e model was also applied to diesel spray combustion computations. It is shown that the use of the RNG model is warranted for spray combustion modeling since the ratio of the turbulent to mean-strain time scales is appreciable due to spray-generated mean flow gradients, and the model introduces a term to account for these effects. Large scale flow structures are predicted which ar...

Turbulence Modeling of Internal Combustion Engines Using RNG κ-ε Models

Turbulent combustion modeling

Turbulence models are developed by supplementing the renormalization group (RNG) approach of Yakhot and Orszag [J. Sci. Comput. 1, 3 (1986)] with scale expansions for the Reynolds stress and production of dissipation terms. The additional expansion parameter (η≡SK/■) is the ratio of the turbulent to mean strain time scale. While low‐order expansions appear to provide an adequate description for the Reynolds stress, no finite truncation of the expansion for the production of dissipation term in powers of η suffices−terms of all orders must be retained. Based on these ideas, a new two‐equation model and Reynolds stress transport model are developed for turbulent shear flows. The models are tested for homogeneous shear flow and flow over a backward facing step. Comparisons between the model predictions and experimental data are excellent.

/pdf/development-of-turbulence-models-for-shear-flows-by-a-double-37jwuzybet.pdf

S. Thangam

Papers

Development of turbulence models for shear flows by a double expansion technique

Development of turbulence models for shear flows by a double expansion technique