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

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

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

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

Computational Fluid Dynamics: Principles and Applications, Third Edition presents students, engineers, and scientists with all they need to gain a solid understanding of the numerical methods and principles underlying modern computation techniques in fluid dynamics By providing complete coverage of the essential knowledge required in order to write codes or understand commercial codes, the book gives the reader an overview of fundamentals and solution strategies in the early chapters before moving on to cover the details of different solution techniques

	This updated edition includes new worked programming examples, expanded coverage and recent literature regarding incompressible flows, the Discontinuous Galerkin Method, the Lattice Boltzmann Method, higher-order spatial schemes, implicit Runge-Kutta methods and parallelization

	An accompanying companion website contains the sources of 1-D and 2-D Euler and Navier-Stokes flow solvers (structured and unstructured) and grid generators, along with tools for Von Neumann stability analysis of 1-D model equations and examples of various parallelization techniques



	
		Will provide you with the knowledge required to develop and understand modern flow simulation codes
	
		Features new worked programming examples and expanded coverage of incompressible flows, implicit Runge-Kutta methods and code parallelization, among other topics
	
		Includes accompanying companion website that contains the sources of 1-D and 2-D flow solvers as well as grid generators and examples of parallelization techniques

Computational Fluid Dynamics: Principles and Applications Ed. 3

An elliptic relaxation model is proposed for the strongly inhomogeneous region near the wall in wall-bounded turbulent shear flow. This model enables the correct kinematic boundary condition to be imposed on the normal component of turbulent intensity. Hence, wall blocking is represented. Means for enforcing the correct boundary conditions on the other components of intensity and on the k — ɛ equations are discussed. The present model agrees quite well with direct numerical simulation (DNS) data. The virtue of the present approach is that arbitrary “damping functions” are not required.

Virendra C. Patel

Papers

Turbulence models for near-wall and low Reynolds number flows - A review