Derivation of fluid dynamics from kinetic theory with the 14-moment approximation

We review the traditional derivation of the fluid-dynamical equations from kinetic theory according to Israel and Stewart. We show that their procedure to close the fluid-dynamical equations of motion is not unique. Their approach contains two approximations, the first being the so-called 14-moment approximation to truncate the single-particle distribution function. The second consists in the choice of equations of motion for the dissipative currents. Israel and Stewart used the second moment of the Boltzmann equation, but this is not the only possible choice. In fact, there are infinitely many moments of the Boltzmann equation which can serve as equations of motion for the dissipative currents. All resulting equations of motion have the same form, but the transport coefficients are different in each case.     

The Boltzmann equation from quantum field theory

We show from first principles the emergence of classical Boltzmann equations from relativistic nonequilibrium quantum field theory as described by the Kadanoff–Baym equations. Our method applies to a generic quantum field, coupled to a collection of background fields and sources, in a homogeneous and isotropic spacetime. The analysis is based on analytical solutions to the full Kadanoff–Baym equations, using the WKB approximation. This is in contrast to previous derivations of kinetic equations that rely on similar physical assumptions, but obtain approximate equations of motion from a gradient expansion in momentum space. We show that the system follows a generalized Boltzmann equation whenever the WKB approximation holds. The generalized Boltzmann equation, which includes off-shell transport, is valid far from equilibrium and in a time dependent background, such as the expanding universe.     

Coarsely grained stochastic Boltzmann equation and its moment equations

The stochastic Boltzmann equation is coarsely grained. The coarsely grained stochastic (CGS) Boltzmann equation has fluctuating terms in its collision term. On the basis of the CGS Boltzmann equation, reduced Grad’s 26 moment equations are derived. Coarsely grained moment equations obtained from the CGS Boltzmann equation show that fluctuating terms remain as nonvanishing terms owing to the nonlinearity in the collision term of the CGS Boltzmann equation. The Navier-Stokes-Fourier law obtained using the CGS Boltzmann equation indicates that the pressure deviator and heat flux include fluctuations of their one-order higher moments.     

Transport theory of dissipative heavy-ion collisions

The non-perturbative quantum-statistical theory of dissipative heavy-ion collisions introduced earlier, is generalized by including explicitly the relative motion of the colliding nuclei. We start from the Liouville equation in the Wigner representation which allows for useful and illustrative interpretations of the resulting quantities and equations. Using the randomness of the coupling matrix elements and the semi-classical approximation for the relative motion we derive a general time-dependent transport equation for the macroscopic Wigner functions (phase-space distribution functions). The limits of weak and strong coupling are discussed.     

The Langevin equation on Lie algebras: Maxwell-Boltzmann is not always the equilibrium

We give a geometric formulation of the Fokker-Planck-Kramer equations for a particle moving on a Lie algebra under the influence of a dissipative and a random force. Special cases of interest are fluid mechanics, the Stochastic Loewner equation and the rigid body. We find that the Boltzmann distribution, although a static solution, is not normalizable when the algebra is not unimodular. This is because the invariant measure of integration in momentum space is not the standard one. We solve the special case of the upper half-plane (hyperboloid) explicitly: there is another equilibrium solution to the Fokker-Planck equation, which is integrable. It breaks rotation invariance.     

Comparison between lattice Boltzmann method and Navier–Stokes high order schemes for computational aeroacoustics

Computational aeroacoustic (CAA) simulation requires accurate schemes to capture the dynamics of acoustic fluctuations, which are weak compared with aerodynamic ones. In this paper, two kinds of schemes are studied and compared: the classical approach based on high order schemes for Navier–Stokes-like equations and the lattice Boltzmann method. The reference macroscopic equations are the 3D isothermal and compressible Navier–Stokes equations. A Von Neumann analysis of these linearized equations is carried out to obtain exact plane wave solutions. Three physical modes are recovered and the corresponding theoretical dispersion relations are obtained. Then the same analysis is made on the space and time discretization of the Navier–Stokes equations with the classical high order schemes to quantify the influence of both space and time discretization on the exact solutions. Different orders of discretization are considered, with and without a uniform mean flow. Three different lattice Boltzmann models are then presented and studied with the Von Neumann analysis. The theoretical dispersion relations of these models are obtained and the error terms of the model are identified and studied. It is shown that the dispersion error in the lattice Boltzmann models is only due to the space and time discretization and that the continuous discrete velocity Boltzmann equation yield the same exact dispersion as the Navier–Stokes equations. Finally, dispersion and dissipation errors of the different kind of schemes are quantitatively compared. It is found that the lattice Boltzmann method is less dissipative than high order schemes and less dispersive than a second order scheme in space with a 3-step Runge–Kutta scheme in time. The number of floating point operations at a given error level associated with these two kinds of schemes are then compared.     

Cellular automaton fluids 1: Basic theory

Continuum equations are derived for the large-scale behavior of a class of cellular automaton models for fluids. The cellular automata are discrete analogues of molecular dynamics, in which particles with discrete velocities populate the links of a fixed array of sites. Kinetic equations for microscopic particle distributions are constructed. Hydrodynamic equations are then derived using the Chapman-Enskog expansion. Slightly modified Navier-Stokes equations are obtained in two and three dimensions with certain lattices. Viscosities and other transport coefficients are calculated using the Boltzmann transport equation approximation. Some corrections to the equations of motion for cellular automaton fluids beyond the Navier-Stokes order are given.     

Stochastic theory for classical and quantum mechanical systems

We formulate from first principles a theory of stochastic processes in configuration space. The fundamental equations of the theory are an equation of motion which generalizes Newton's second law and an equation which expresses the condition of conservation of matter. Two types of stochastic motion are possible, both described by the same general equations, but leading in one case to classical Brownian motion behavior and in the other to quantum mechanical behavior. The Schrödinger equation, which is derived here with no further assumption, is thus shown to describe a specific stochastic process. It is explicitly shown that only in the quantum mechanical process does the superposition of probability amplitudes give rise to interference phenomena; moreover, the presence of dissipative forces in the Brownian motion equations invalidates the superposition principle. At no point are any special assumptions made concerning the physical nature of the underlying stochastic medium, although some suggestions are discussed in the last section.     

Dissipative Model of the Universe

In the present study a trial is carried out to generalize the non-relativistic-dynamic model of cosmology put forward by O. Heckmann. The generalization is performed under the framework of the Gyarmati principle of irreversible thermodynamics for anisotropic inhomogeneous and viscous case. The equation of motion of the dissipative Universe will be a Navier-Stokes tensor equation leading to a Riccatian differential equation general in one dimension. The equations of the Heckmann model compatible with standard cosmology can be obtained from this by means of further simplifying assumption.     

Convergence from Boltzmann to Landau Processes with Soft Potential and Particle Approximations

Our aim in this paper is to show how a probabilistic interpretation of the Boltzmann and Landau equations gives a microscopic understanding of these equations. We firstly associate stochastic jump processes with the Boltzmann equations we consider. Then we renormalize these equations following asymptotics which make prevail the grazing collisions, and prove the convergence of the associated Boltzmann jump processes to a diffusion process related to the Landau equation. The convergence is pathwise and also implies a convergence at the level of the partial differential equations. The best feature of this approach is the microscopic understanding of the transition between the Boltzmann and the Landau equations, by an accumulation of very small jumps. We deduce from this interpretation an approximation result for a solution of the Landau equation via colliding stochastic particle systems. This result leads to a Monte-Carlo algorithm for the simulation of solutions by a conservative particle method which enables to observe the transition from Boltzmann to Landau equations. Numerical results are given.     

Deformed quantum harmonic oscillator with diffusion and dissipation

A master equation for the deformed quantum harmonic oscillator interacting with a dissipative environment, in particular with a thermal bath, is derived in the microscopic model by using perturbation theory. The coefficients of the master equation and of equations of motion for observables depend on the deformation function. The steady-state solution of the equation for the density matrix in the number representation is obtained and the equilibrium energy of the deformed harmonic oscillator is calculated in the approximation of small deformation.     

Hilbert Expansion from the Boltzmann Equation to Relativistic Fluids

We study the local-in-time hydrodynamic limit of the relativistic Boltzmann equation using a Hilbert expansion. More specifically, we prove the existence of local solutions to the relativistic Boltzmann equation that are nearby the local relativistic Maxwellians. The Maxwellians are constructed from a class of solutions to the relativistic Euler equations that includes a large subclass of near-constant, non-vacuum fluid states. In particular, for small Knudsen number, these solutions to the relativistic Boltzmann equation have dynamics that are effectively captured by corresponding solutions to the relativistic Euler equations.     

DESCRIPTIONS OF THE EQUILIBRATION PROCESS OF THE INTRINSIC DEGREES OF FREEDOM AND DISSIPATIVE PROCESS OF THE NUCLEAR COLLECTIVE MOTION

In this paper the Hamiltonian model is used for studying the nuclear dynamics by taking both the one-body and two-body interaction mechanisms into account. On the basis of the Von Neuman equation the coupling between the collective motion and the single particle degrees of freedom is discussed. Thus, the equations obtained are physically transparent and easy for numerical computations. They may be useful for describing the dissipative process of the nuclear collective motion as well as the equilibration process of the intrinsic degrees of freedom.     

Quantum noise generated by four-wave mixing process within a fiber ring resonator

We have analyzed the dissipative effect of the entangled photons generated by a nonlinear optical ring resonator from a non-degenerate four-wave mixing (FWM) process. The system and reservoir Hamiltonian are established in the interaction picture. To eliminate the reservoir operators, the Markov approximation is used and result them in a Linblad form in the master equation. Consequently, the positive P representation can recast this equation to the Fokker-Planck equation, and then the stochastic differential equations i.e., the entangled photon state equation of motion for photons propagating in the fiber, are obtained and easy to analyze numerically. Results obtained have shown that the entangled strength measurement depends on three main factors; first the nonlinear susceptibility of the third harmonic generation, second the damping rate that represents loss of energy from the system to the reservoir, and final the diffusion of fluctuations in the reservoir into the entangled photon modes.     

The Fluctuation Theorem

The question of how reversible microscopic equations of motion can lead to irreversible macroscopic behaviour has been one of the central issues in statistical mechanics for more than a century. The basic issues were known to Gibbs. Boltzmann conducted a very public debate with Loschmidt and others without a satisfactory resolution. In recent decades there has been no real change in the situation. In 1993 we discovered a relation, subsequently known as the Fluctuation Theorem (FT), which gives an analytical expression for the probability of observing Second Law violating dynamical fluctuations in thermostatted dissipative non-equilibrium systems. The relation was derived heuristically and applied to the special case of dissipative non-equilibrium systems subject to constant energy 'thermostatting'. These restrictions meant that the full importance of the Theorem was not immediately apparent. Within a few years, derivations of the Theorem were improved but it has only been in the last few of years that the generality of the Theorem has been appreciated. We now know that the Second Law of Thermodynamics can be derived assuming ergodicity at equilibrium, and causality. We take the assumption of causality to be axiomatic. It is causality which ultimately is responsible for breaking time reversal symmetry and which leads to the possibility of irreversible macroscopic behaviour. The Fluctuation Theorem does much more than merely prove that in large systems observed for long periods of time, the Second Law is overwhelmingly likely to be valid. The Fluctuation Theorem quantifies the probability of observing Second Law violations in small systems observed for a short time. Unlike the Boltzmann equation, the FT is completely consistent with Loschmidt's observation that for time reversible dynamics, every dynamical phase space trajectory and its conjugate time reversed 'anti-trajectory', are both solutions of the underlying equations of motion. Indeed the standard proofs of the FT explicitly consider conjugate pairs of phase space trajectories. Quantitative predictions made by the Fluctuation Theorem regarding the probability of Second Law violations have been confirmed experimentally, both using molecular dynamics computer simulation and very recently in laboratory experiments.     

Stable first-order particle-frame relativistic hydrodynamics for dissipative systems

We propose a stable first-order relativistic dissipative hydrodynamic equation in the particle frame (Eckart frame) for the first time. The equation to be proposed was in fact previously derived by the authors and a collaborator from the relativistic Boltzmann equation. We demonstrate that the equilibrium state is stable with respect to the time evolution described by our hydrodynamic equation in the particle frame. Our equation may be a proper starting point for constructing second-order causal relativistic hydrodynamics, to replace Eckart's particle-flow theory.     

A numerical study of an expanding plasma of quarks in a chiral model

We numerically solve the transport equations for a quark gas described by the the Nambu-Jona-Lasinio model. The mean field equations of motion, which consist of the Vlasov equation for the density and the gap equation for the mean field, are discussed, and energy and momentum conservation are proven. Numerical solutions of the partial differential equations are obtained by applying finite difference methods. For an expanding fireball the light quark mass evolves from small values initially to the value of 350 MeV. This leads to a depletion of the high energy part of the quark spectrum and an enhancement at low momenta. When collisions are included one obtains an equation of the Boltzmann type, where the transition amplitudes depend on the properties of the medium. These equations are given for flavor SU(3), i.e. including strangeness. They are solved numerically in the relaxation time approximation and the time evolution of various observables is given. Medium effects in the relaxation times do not significantly influence the shape of the spectra. The mass of the strange quark changes little during the expansion. The strangeness yield and the slope temperatures of the final spectra are studied as a function of the size of the initial fireball.     

A STOCHASTIC OPTIMAL SEMI-ACTIVE CONTROL STRATEGY FOR ER/MR DAMPERS

A stochastic optimal semi-active control strategy for randomly excited systems using electrorheological/magnetorheological (ER/MR) dampers is proposed. A system excited by random loading and controlled by using ER/MR dampers is modelled as a controlled, stochastically excited and dissipated Hamiltonian system with n degrees of freedom. The control forces produced by ER/MR dampers are split into a passive part and an active part. The passive control force is further split into a conservative part and a dissipative part, which are combined with the conservative force and dissipative force of the uncontrolled system, respectively, to form a new Hamiltonian and an overall passive dissipative force. The stochastic averaging method for quasi-Hamiltonian systems is applied to the modified system to obtain partially completed averaged Itô stochastic differential equations. Then, the stochastic dynamical programming principle is applied to the partially averaged Itô equations to establish a dynamical programming equation. The optimal control law is obtained from minimizing the dynamical programming equation subject to the constraints of ER/MR damping forces, and the fully completed averaged Itô equations are obtained from the partially completed averaged Itô equations by replacing the control forces with the optimal control forces and by averaging the terms involving the control forces. Finally, the response of semi-actively controlled system is obtained from solving the final dynamical programming equation and the Fokker-Planck-Kolmogorov equation associated with the fully completed averaged Itô equations of the system. Two examples are given to illustrate the application and effectiveness of the proposed stochastic optimal semi-active control strategy.