A new and very general technique for simulating solid–fluid suspensions is described; its most important feature is that the computational cost scales linearly with the number of particles. The method combines Newtonian dynamics of the solid particles with a discretized Boltzmann equation for the fluid phase; the many-body hydrodynamic interactions are fully accounted for, both in the creeping-flow regime and at higher Reynolds numbers. Brownian motion of the solid particles arises spontaneously from stochastic fluctuations in the fluid stress tensor, rather than from random forces or displacements applied directly to the particles. In this paper, the theoretical foundations of the technique are laid out, illustrated by simple analytical and numerical examples; in a companion paper (Part 2), extensive numerical tests of the method, for stationary flows, time-dependent flows, and finite-Reynolds-number flows, are reported.

Numerical simulations of particulate suspensions via a discretized boltzmann equation: part 1. theoretical foundation

Some Preliminary Thoughts * Part I: Inviscid Hypersonic Flow * Hypersonic Shock and Expansion-Wave Relations * Local Surface Inclination Methods * Hypersonic Inviscid Flowfields: Approximate Methods * Hypersonic Inviscid Flowfields: Exact Methods * Part II: Viscous Hypersonic Flow * Viscous Flow: Basic Aspects, Boundary Layer Results, and Aerodynamic Heating * Hypersonic Viscous Interactions * Computational Fluid Dynamic Solutions of Hypersonic Viscous Flows * Part III: High-Temperature Gas Dynamics * High-Temperature Gas Dynamics: Some Introductory Considerations * Some Aspects of the Thermodynamics of Chemically Reacting Gases (Classical Physical Chemistry) * Elements of Statistical Thermodynamics * Elements of Kinetic Theory * Chemical Vibrational Nonequilibrium * Inviscid High-Temperature Equilibrium Flows * Inviscid High-Temperature Nonequilibrium Flows * Kinetic Theory Revisited: Transport Properties in High-Temperature Gases * Viscous High-Temperature Flows * Introduction to Radiative Gas Dynamics.

Hypersonic and high temperature gas dynamics

Monte Carlo is one of the most versatile and widely used numerical methods. Its convergence rate, O(N−1/2), is independent of dimension, which shows Monte Carlo to be very robust but also slow. This article presents an introduction to Monte Carlo methods for integration problems, including convergence theory, sampling methods and variance reduction techniques. Accelerated convergence for Monte Carlo quadrature is attained using quasi-random (also called low-discrepancy) sequences, which are a deterministic alternative to random or pseudo-random sequences. The points in a quasi-random sequence are correlated to provide greater uniformity. The resulting quadrature method, called quasi-Monte Carlo, has a convergence rate of approximately O((logN)kN−1). For quasi-Monte Carlo, both theoretical error estimates and practical limitations are presented. Although the emphasis in this article is on integration, Monte Carlo simulation of rarefied gas dynamics is also discussed. In the limit of small mean free path (that is, the fluid dynamic limit), Monte Carlo loses its effectiveness because the collisional distance is much less than the fluid dynamic length scale. Computational examples are presented throughout the text to illustrate the theory. A number of open problems are described.

http://dsec.pku.edu.cn/~tieli/notes/numer_anal/MCQMC_Caflisch.pdf

Monte Carlo and quasi-Monte Carlo methods

A moving-particle semi-implicit (MPS) method for simulating fragmentation of incompressible fluids is presented. The motion of each particle is calculated through interactions with neighboring particles covered with the kernel function. Deterministic particle interaction models representing gradient, Laplacian, and free surfaces are proposed. Fluid density is implicitly required to be constant as the incompressibility condition, while the other terms are explicitly calculated. The Poisson equation of pressure is solved by the incomplete Cholesky conjugate gradient method. Collapse of a water column is calculated using MPS. The effect of parameters in the models is investigated in test calculations. Good agreement with an experiment is obtained even if fragmentation and coalescence of the fluid take place.

Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

The direct simulation Monte Carlo (or DSMC) method has, in recent years, become widely used in engineering and scientific studies of gas flows that involve low densities or very small physical dimensions. This method is a direct physical simulation of the motion of representative molecules, rather than a numerical solution of the equations that provide a mathematical model of the flow. These computations are no longer expensive and the period since the 1976 publication of the original Molecular Gas Dynamics has seen enormous improvements in the molecular models, the procedures, and the implementation strategies for the DSMC method. The molecular theory of gas flows is developed from first principles and is extended to cover the new models and procedures. Note: The disk that originally came with this book is no longer available. However, the same information is available from the author's website (http://gab.com.au/)

Molecular Gas Dynamics and the Direct Simulation of Gas Flows

Molecular gas dynamics

The computational and physical requirements for the successful application of the direct simulation Monte-Carlo method to typical engineering flows, as opposed to earlier applications to artificial test cases, are discussed. An outline is given of the changes that have occurred in the program structure as a result of both altered requirements and advances in computer technology. The accumulated experience from applications of the classical molecular models is reviewed. It is shown that requirements for the accurate modeling of engineering flows are best met by a new molecular model called the variable hard sphere, or VHS, model, The necessary kinetic theory for the application of this model is presented.

Monte-carlo simulation in an engineering context

Approach to Translational Equilibrium in a Rigid Sphere Gas

The direct simulation Monte Carlo method for the numerical solution of problems in rarefied gas dynamics is described and discussed. It is shown that the procedures adopted in this method can be directly related to the Boltzmann equation and that the two are entirely consistent. It is concluded that the results obtained from the method constitute a solution of the Boltzmann equation.

Graeme A. Bird

Papers

Molecular Gas Dynamics and the Direct Simulation of Gas Flows

Molecular gas dynamics

Monte-carlo simulation in an engineering context

Approach to Translational Equilibrium in a Rigid Sphere Gas

Direct Simulation and the Boltzmann Equation