Contents:
The methods considered are those at the foundation of real-world atmospheric or ocean models, with the focus being on the essential mathematical properties of each method. The fundamental character of each scheme is examined in prototypical fluid-dynamical problems like tracer transport, chemically reacting flow, shallow-water waves, and waves in an internally stratified fluid.
The book includes exercises and is well illustrated with figures linking theoretical analyses to results from actual computations.
Changes from the first edition include new chapters, discussions and updates throughout. It is a contemporary and worthy addition to the still-sparse list of quality graduate-level references on the numerical solution of PDEs. Durran 1 1. Buy options. As a first section, the first Stoke problem is considered numerically by introducing the finite difference method.
In the second section, natural convection heat transfer heated from a vertical plate with uniform heat flux is introduced together with the method how to obtain the system of ordinary differential equations. In the third example, linear stability analysis for the onset of secondary flow during the Taylor-Couette flow is numerically treated using the HSMAC method. The governing equation for the fluid flow is known as Navier-Stokes equation, which is however difficult to solve analytically; and therefore, a lot of numerical techniques have been proposed and developed.
Nevertheless various complex flow phenomena such as turbulent flow, multi-phase flow, compressible flow, combustion, and phase change encountered in the fields of engineering would have still difficulties to circumvent even using both present computational resources and numerical techniques. The present chapter devotes not to elucidate such complex phenomena, but to introduce rather simplified fluid flow by using the finite difference method.
One focuses on incompressible flows, in which physical properties such as the viscosity, the thermal conductivity, the specific heat are constant and even the fluid density is not a thermodynamic variable. The former situation is known the low Mach number approximation , while the latter one the Boussinesq approximation. Another simplification on the incompressible flows is the reduction of dimension due to the characteristic of similarity and periodicity. At the stage of onset of instability, the non-linear term is negligible and therefore the function of flow field is separated into the amplitude part and periodic part, respectively.
This makes the effort on numerical analysis to reduce significantly and also to contribute the augmentation of accuracy of the results. This chapter consists of three main bodies. First, a numerical technique for solving the boundary value problem called the first Stokes problem or the Rayleigh problem [ 1 ] is introduced.
The differential equation is transferred into an ordinary equation and it is solved by a finite difference method using the Jacobi method. Second, similar solution of natural convection heat transfer heated from a vertical plate with uniform heat flux is introduced together with the method how to obtain the system of ordinary differential equations. The obtained Nusselt numbers are compared with some previous studies.
Third, for example, of the linear stability analysis, one shows that the HSMAC method can be applied to obtain the critical values for the onset of secondary flow such as the Taylor-Couette flow. The Eigen functions of flow and pressure fields are visualized. An infinite length plate is set in a stationary fluid as an initial condition. Let us consider the situation that the infinite length plate suddenly moves along its parallel direction at a constant speed uw.
Read the latest chapters of Handbook of Numerical Analysis at ScienceDirect. com, Elsevier's leading platform of peer-reviewed scholarly literature. Handbook of Numerical Analysis · Latest volumeAll Search in this handbook. Numerical Methods for Solids (Part 3) Numerical Methods for Fluids (Part 1).
This problem was first solved by Stokes [ 2 ] in his famous treatment of the pendulum. Since Lord Rayleigh [ 3 ] also treated this flow, it is often called the Rayleigh problem in the literature. One takes that x is the plate movement direction and y is distance from the plate.
Since the velocity component perpendicular to the plate v is zero, the momentum equation is simplified and is shown as a diffusion equation.
The boundary conditions for this partial differential equation are as follows:. As a consequence, one needs to solve this boundary value problem. The theoretical solution can be easily obtained and expressed using the error function. Hence, the boundary condition shown below is used instead of Eq.
The approximated velocity profile is expressed by connecting these values smoothly. When the second-order central difference method is used, Eq. Here, N is total number of grids and in this chapter, the first grid point starts from 1 as its definition. The above equation becomes.
The boundary condition 7 is modified. This kind of tridiagonal matrix is often seen and can be solved by a direct numerical method, such as Tomas method. However, the rank of the matrix is usually extremely large and one introduces an iterative method for solving the king-size matrix. In general, the rank of the matrix appearing in computational fluid dynamics CFD is large and iterative methods such as Jacobi, Gauss-Seidel, or successive over relaxation SOR methodare employed.
In this subsection, the Jacobi method is explained. Objectives: To realise an ambitious project close to the real engineer work.
To use a wide scope of knowlege for an concrete project. To be part of a team work. Program: The choice and the definition of the project is left to the initiave groups of two persons gathered in team of about ten, in the framework of well defined constraints concerning scientific domains and organisational aspects. For instance, the two-people group projects must find integration in the ten-people team subject.
Each team has to choose and received an invited speaker for a seminar and an expertise of the team project.
Jomaa and C. Advanced Transport Phenomena Heitmann, Stefan However, only four equations among the eight equations are necessary to solve in this problem because of the symmetricity and anti-symmetricity of the complex variables. This makes the effort on numerical analysis to reduce significantly and also to contribute the augmentation of accuracy of the results. The material is meant to serve as a prerequisite for students who might go on to take additional courses in computational mechanics, computational fluid dynamics, or computational electromagnetics.
The presentation of the two-people groups and team work is made through a multimedia support with hypertext file Intranet in order to enhance the group interactions during the course. The final presentation of the course projects is made in front of a large public. Bibliography: [1] Dynamique des fluides. Mecanique des fluides: Elements d'un premier parcours [2] P. Physical fluid dynamics. Objectives: To group and give physical and mathematical knowledge needed for mastering the statistical approach of incompressible turbulence.
Program: Physical properties of the solutions of Navier-Stokes equations: Advection, diffusion, dissipation, instability and transition. Diffusion by continuous motion: turbulent mixing, vortex stretching, wall transfer. Examples of turbulent motions. Stochastic tools :correlation, spectrum. One point characteristic function transport equation: application to mean and second order moments. Keywords: Turbulent structure, coherent motion, energy transfer, homogeneous turbulence, correlation, spectrum. Bibliography: [1] A first course in Turbulence. Objectives: To get a phenomenological insight on eddy structures to achieve the understanding of energy transfer in incompressible turbulence.
Coherent structures: origin, pairing. Homogeneous and isotropic turbulence: correlations and spectra dynamics. Bibliography: [1] The structure of turbulent shear flow. Press, RODI, J. Eyrolles, Objectives: Knowledge of single point closures to get the skills to decide on the suitable model to predict high Reynolds number turbulent flows.
Program: Classification of turbulent flows. Homogeneous turbulence, turbulent shear flows. Basic mecanisms: isotropic relaxation, mean motion effects contraction, shear,pure deformation,rotation ,pressure- redistribution, diffusion, dissipation. Overview of currents methods Direct numerical simulation, Large eddy simulation, spectral and single point closures.