Non-Local Scattering Kernel and the Hydrodynamic Limit

In this paper we study the interaction of a fluid with a wall in the framework of the kinetic theory. We consider the possibility that the fluid molecules can penetrate the wall to be reflected by the inner layers of the wall. This results in a scattering kernel which is a non-local generalization of the classical Maxwell scattering kernel. The proposed scattering kernel satisfies a global mass conservation law and a generalized reciprocity relation. We study the hydrodynamic limit performing a Knudsen layer analysis, and derive a new class of (weakly) nonlocal boundary conditions to be imposed to the Navier–Stokes equations.     

Navier–Stokes Transport Coefficients of <Emphasis Type="Italic">d</Emphasis>-Dimensional Granular Binary Mixtures at Low Density

The Navier–Stokes transport coefficients for binary mixtures of smooth inelastic hard disks or spheres under gravity are determined from the Boltzmann kinetic theory by application of the Chapman–Enskog method for states near the local homogeneous cooling state. It is shown that the Navier–Stokes transport coefficients are not affected by the presence of gravity. As in the elastic case, the transport coefficients of the mixture verify a set of coupled linear integral equations that are approximately solved by using the leading terms in a Sonine polynomial expansion. The results reported here extend previous calculations (Garzó, V., Dufty, J.W. in Phys. Fluids 14:1476–1490, 2002) to an arbitrary number of dimensions and provide explicit expressions for the seven Navier–Stokes transport coefficients in terms of the coefficients of restitution and the masses, composition, and sizes of the constituents of the mixture. In addition, to check the accuracy of our theory, the inelastic Boltzmann equation is also numerically solved by means of the direct simulation Monte Carlo method to evaluate the diffusion and shear viscosity coefficients for hard disks. The comparison shows a good agreement over a wide range of values of the coefficients of restitution and the parameters of the mixture (masses and sizes).     

Global Existence of Classical Solutions to the Vlasov-Poisson-Boltzmann System

The time evolution of the distribution function for the charged particles in a dilute gas is governed by the Vlasov–Poisson–Boltzmann system when the force is self-induced and its potential function satisfies the Poisson equation. In this paper, we give a satisfactory global existence theory of classical solutions to this system when the initial data is a small perturbation of a global Maxwellian. Moreover, the convergence rate in time to the global Maxwellian is also obtained through the energy method. The proof is based on the theory of compressible Navier–Stokes equations with forcing and the decomposition of the solutions to the Boltzmann equation with respect to the local Maxwellian introduced in [23] and elaborated in [31].     

Gas Flow in Microchannels – A Lattice Boltzmann Method Approach

Gas flow in microchannels can often encounter tangential slip motion at the solid surface even under creeping flow conditions. To simulate low speed gas flows with Knudsen numbers extending into the transition regime, alternative methods to both the Navier–Stokes and direct simulation Monte Carlo approaches are needed that balance computational efficiency and simulation accuracy. The lattice Boltzmann method offers an approach that is particularly suitable for mesoscopic simulation where details of the molecular motion are not required. In this paper, the lattice Boltzmann method has been applied to gas flows with finite Knudsen number and the tangential momentum accommodation coefficient has been implemented to describe the gas-surface interactions. For fully-developed channel flows, the results of the present method are in excellent agreement with the analytical slip-flow solution of the Navier–Stokes equations, which are valid for Knudsen numbers less than 0.1. The present paper demonstrates that the lattice Boltzmann approach is a promising alternative simulation tool for the design of microfluidic devices.     

Asymptotic Analysis of Lattice Boltzmann Boundary Conditions

In this article, we use a general method for the analysis of finite difference schemes to investigate lattice Boltzmann algorithms for Navier–Stokes problems with Dirichlet boundary conditions. Several link based boundary conditions for commonly used lattice Boltzmann BGK models are considered. With our method, the accuracy of the algorithms can be exactly predicted. Moreover, the analytical results can be used to construct new algorithms which is demonstrated with a corrected bounce back rule that requires only local evaluations but still yields second order accuracy for the velocity. The analysis is applicable to general geometries and instationary flows     

Transition to Longitudinal Instability of Detonation Waves is Generically Associated with Hopf Bifurcation to Time-Periodic Galloping Solutions

We show that transition to longitudinal instability of strong detonation solutions of reactive compressible Navier–Stokes equations is generically associated with Hopf bifurcation to nearby time-periodic “galloping”, or “pulsating”, solutions, in agreement with physical and numerical observation. In the process, we determine readily numerically verifiable stability and bifurcation conditions in terms of an associated Evans function, and obtain the first complete nonlinear stability result for strong detonations of the reacting Navier–Stokes equations, in the limit as amplitude (hence also heat release) goes to zero. The analysis is by pointwise semigroup techniques introduced by the authors and collaborators in previous works.     

A particle–particle hybrid method for kinetic and continuum equations

We present a coupling procedure for two different types of particle methods for the Boltzmann and the Navier–Stokes equations. A variant of the DSMC method is applied to simulate the Boltzmann equation, whereas a meshfree Lagrangian particle method, similar to the SPH method, is used for simulations of the Navier–Stokes equations. An automatic domain decomposition approach is used with the help of a continuum breakdown criterion. We apply adaptive spatial and time meshes. The classical Sod’s 1D shock tube problem is solved for a large range of Knudsen numbers. Results from Boltzmann, Navier–Stokes and hybrid solvers are compared. The CPU time for the hybrid solver is 3–4 times faster than for the Boltzmann solver.     

Boundary Layers for the Navier–Stokes Equations¶of Compressible Fluids

The global unique solvability is proved for the Navier–Stokes equations of compressible fluids for the one-dimensional spiral flows between two circular cylinders. The zero shear viscosity limit μ→ 0 is justified. The value O(μ^α), 0 < α < 1/2, is established for the boundary layer thickness. Received: 11 December 1998 / Accepted: 9 June 1999     

Long-Time Stability of Multi-Dimensional Noncharacteristic Viscous Boundary Layers

We establish long-time stability of multi-dimensional noncharacteristic boundary layers of a class of hyperbolic–parabolic systems including the compressible Navier–Stokes equations with inflow [outflow] boundary conditions, under the assumption of strong spectral, or uniform Evans, stability. Evans stability has been verified for small-amplitude layers by Guès, Métivier, Williams, and Zumbrun. For large-amplitude layers, it may be efficiently checked numerically, as done in the one-dimensional case by Costanzino, Humpherys, Nguyen, and Zumbrun.     

Stationary Flow Past a Semi-Infinite Flat Plate: Analytical and Numerical Evidence for a Symmetry-Breaking Solution

We consider the question of the existence of stationary solutions for the Navier Stokes equations describing the flow of a incompressible fluid past a semi-infinite flat plate at zero incidence angle. By using ideas from the theory of dynamical systems we analyze the vorticity equation for this problem and show that a symmetry-breaking term fits naturally into the downstream asymptotic expansion of a solution. Finally, in order to check that our asymptotic expressions can be completed to a symmetry-breaking solution of the Navier–Stokes equations we solve the problem numerically by using our asymptotic results to prescribe artificial boundary conditions for a sequence of truncated domains. The results of these numerical computations a clearly compatible with the existence of such a solution. Mathematics Subject Classification (2000). 76D05, 76D25, 76M10, 41A60, 35Q35 Supported in part by the Fonds National Suisse.     

Simulation of Shock Wave Diffraction over 90° Sharp Corner in Gases of Arbitrary Statistics

The unsteady shock wave diffraction over a 90° sharp corner in gases of arbitrary particle statistics is simulated using an accurate and direct algorithm for solving the semiclassical Boltzmann equation with relaxation time approximation in phase space. The numerical method is based on the usage of discrete ordinate method for discretizing the velocity space of the distribution function; whereas a second order accurate TVD scheme (Harten in J. Comput. Phys. 49(3):357–393, 1983) with Van Leer’s limiter (J. Comput. Phys. 32(1):101–136, 1979) is used for evolving the solution in physical space and time. The specular reflection surface boundary condition is assumed. The complete diffraction patterns are recorded using various flow property contours. Different range of relaxation times approximately corresponding to continuum, slip and transitional regimes are considered and the equilibrium Euler limit solution is also computed for comparison. The effects of gas particles that obey the Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac statistics are examined and depicted.     

Comparing the first- and second-order theories of relativistic dissipative fluid dynamics using the 1+1 dimensional relativistic flux corrected transport algorithm

Focusing on the numerical aspects and accuracy we study a class of bulk viscosity driven expansion scenarios using the relativistic Navier–Stokes and truncated Israel–Stewart form of the equations of relativistic dissipative fluids in 1+1 dimensions. The numerical calculations of conservation and transport equations are performed using the numerical framework of flux corrected transport. We show that the results of the Israel–Stewart causal fluid dynamics are numerically much more stable and smoother than the results of the standard relativistic Navier–Stokes equations.     

Asymptotic analysis of the lattice Boltzmann equation

In this article we analyze the lattice Boltzmann equation (LBE) by using the asymptotic expansion technique. We first relate the LBE to the finite discrete-velocity model (FDVM) of the Boltzmann equation with the diffusive scaling. The analysis of this model directly leads to the incompressible Navier–Stokes equations, as opposed to the compressible Navier–Stokes equations obtained by the Chapman–Enskog analysis with convective scaling. We also apply the asymptotic analysis directly to the fully discrete LBE, as opposed to the usual practice of analyzing a continuous equation obtained through the Taylor-expansion of the LBE. This leads to a consistency analysis which provides order-by-order information about the numerical solution of the LBE. The asymptotic technique enables us to analyze the structure of the leading order errors and the accuracy of numerically derived quantities, such as vorticity. It also justifies the use of Richardson’s extrapolation method. As an example, a two-dimensional Taylor-vortex flow is used to validate our analysis. The numerical results agree very well with our analytic predictions.     

Integral Representation of the Linear Boltzmann Operator for Granular Gas Dynamics with Applications

We investigate the properties of the collision operator Q associated to the linear Boltzmann equation for dissipative hard-spheres arising in granular gas dynamics. We establish that, as in the case of non-dissipative interactions, the gain collision operator is an integral operator whose kernel is made explicit. One deduces from this result a complete picture of the spectrum of Q in an Hilbert space setting, generalizing results from T. Carleman (Publications Scientifiques de l’Institut Mittag-Leffler, vol. 2, 1957) to granular gases. In the same way, we obtain from this integral representation of Q ⁺ that the semigroup in L ¹(ℝ³×ℝ³,dx ⊗ dv) associated to the linear Boltzmann equation for dissipative hard spheres is honest generalizing known results from Arlotti (Acta Appl. Math. 23:129–144, 1991).     

Symmetry of the Unsteady Linearized Boltzmann Equation in a Fixed Bounded Domain

A symmetric relation between time-dependent problems described by the linearized Boltzmann equation is obtained for a gas in a fixed bounded domain. General representations of the total mass, momentum, and energy in the domain, as well as their fluxes through the boundary, in terms of an appropriate Green function are derived from that relation. Several application examples are presented. Similarities to the fluctuation–dissipation theorem in the linear response theory and its generalization to gas systems of arbitrary Knudsen numbers are also discussed. The present paper is an extension of the previous work of the author (Takata in J. Stat. Phys. 136: 751–784, 2009) to time-dependent problems.     

Entire Solutions of Hydrodynamical Equations with Exponential Dissipation

We consider a modification of the three-dimensional Navier–Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as e^|k|/k_d{{{\rm e}^{|k|/k_{\rm d}}}} at high wavenumbers |k|. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than e^{-C(k/k_d) ln(|k|/k_d)}{{{\rm e}^{-C(k/k_{\rm d})\,{\rm ln}(|k|/k_{\rm d})}}} for any C < 1/(2 ln 2). The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with C = C _* = 1/ ln 2. The same behavior with a universal constant C _* is conjectured for the Navier–Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier–Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.     

A high-order gas-kinetic Navier–Stokes flow solver

The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge–Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge–Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to its spatial and temporal decoupling. Many recently developed high-order methods require a Navier–Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier–Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier–Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such a high-order method. The necessity to couple spatial and temporal evolution nonlinearly in the flux evaluation can be clearly observed through the numerical performance of the scheme for the viscous flow computations.     

One-Dimensional Hard-Rod Caricature of Hydrodynamics: “Navier–Stokes Correction” for Local Equilibrium Initial States

A one-dimensional system, consisting of identical hard-rod particles of length $a$ is studied in the hydrodynamical limit. A “Navier–Stokes correction” to the Euler equation is found for an initial local equilibrium family of states , of constant density. The correction is given, at , by a non-linear second order differential operator acting on f(q,v), the hydrodynamical density at a point of the “species” of fluid with velocity . Received: 14 October 1996 / Accepted: 13 February 1997     

Open and traction boundary conditions for the incompressible Navier–Stokes equations

We present numerical schemes for the incompressible Navier–Stokes equations (NSE) with open and traction boundary conditions. We use pressure Poisson equation (PPE) formulation and propose new boundary conditions for the pressure on the open or traction boundaries. After replacing the divergence free constraint by this pressure Poisson equation, we obtain an unconstrained NSE. For Stokes equation with open boundary condition on a simple domain, we prove unconditional stability of a first order semi-implicit scheme where the pressure is treated explicitly and hence is decoupled from the computation of velocity. Using either boundary condition, the schemes for the full NSE that treat both convection and pressure terms explicitly work well with various spatial discretizations including spectral collocation and C⁰ finite elements. Moreover, when Reynolds number is of O(1) and when the first order semi-implicit time stepping is used, time step size of O(1) is allowed in benchmark computations for the full NSE. Besides standard stability and accuracy check, various numerical results including flow over a backward facing step, flow past a cylinder and flow in a bifurcated tube are reported. Numerically we have observed that using PPE formulation enables us to use the velocity/pressure pairs that do not satisfy the standard inf–sup compatibility condition. Our results extend that of Johnston and Liu [H. Johnston, J.-G. Liu, Accurate, stable and efficient Navier–Stokes solvers based on explicit treatment of the pressure term. J. Comp. Phys. 199 (1) (2004) 221–259] which deals with no-slip boundary conditions only.