Hydrodynamic Limit of a B.G.K. Like Model on Domains with Boundaries and Analysis of Kinetic Boundary Conditions for Scalar Multidimensional Conservation Laws

In this paper we study the hydrodynamic limit of a B.G.K. like kinetic model on domains with boundaries via BV _loc theory. We obtain as a consequence existence results for scalar multidimensional conservation laws with kinetic boundary conditions. We require that the initial and boundary data satisfy the optimal assumptions that they all belong to L ¹∩L ^∞ with the additional regularity assumptions that the initial data are in BV _loc . We also extend our hydrodynamic limit analysis to the case of a generalized kinetic model to account for forces effects and we obtain as a consequence the existence theory for conservation laws with source terms and kinetic boundary conditions.     

Time-dependent aerodynamic non-LTE radiative transfer

A method is developed for solving simultaneously in one dimension the equation of transfer for non-LTE spectral line radiation and the time-dependent equations specifying conservation of mass, energy and linear momentum. In particular, we illustrate the method on a ‘simple’ time-dependent problem in which a pulsating disturbance at some point in a model homogeneous atmosphere propagates towards the surface and steepens into a shock. The resulting emergent intensities show rather dramatic changes over very small time intervals due to the effect of the velocity, density and temperature distributions on the radiative absorption properties of the gas, and thus emphasises the need to solve the above-mentioned four basic equations if one is to obtain physically realistic model atmospheres experiencing initial disturbances.     

A scalar model for MHD turbulence

A model for homogeneous MHD turbulence is proposed. Nonlinear interactions acts between nearest neighbours in a discretized wavenumbers space and conservation properties (total energy and v-b correlation) are verified. The model can be truncated at will. With three modes, a bifurcation analysis is given. In the turbulent case (dissipation and kinetic forcing are present) one obtains time fluctuations at all scales and time-averaged power law spectra, the small scales exhibit intermittency effects. Typical MHD phenomena such as the dynamo effect or the increase of v-b correlation in the decaying cases are also observed.     

Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

Uncertainty quantification for systems of conservation laws

Uncertainty quantification through stochastic spectral methods has been recently applied to several kinds of non-linear stochastic PDEs. In this paper, we introduce a formalism based on kinetic theory to tackle uncertain hyperbolic systems of conservation laws with Polynomial Chaos (PC) methods. The idea is to introduce a new variable, the entropic variable, in bijection with our vector of unknowns, which we develop on the polynomial basis: by performing a Galerkin projection, we obtain a deterministic system of conservation laws. We state several properties of this deterministic system in the case of a general uncertain system of conservation laws. We then apply the method to the case of the inviscid Burgers’ equation with random initial conditions and we present some preliminary results for the Euler system. We systematically compare results from our new approach to results from the stochastic Galerkin method. In the vicinity of discontinuities, the new method bounds the oscillations due to Gibbs phenomenon to a certain range through the entropy of the system without the use of any adaptative random space discretizations. It is found to be more precise than the stochastic Galerkin method for smooth cases but above all for discontinuous cases.     

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. These observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.     

A statistical limit in the solution of the nonlinear Schrödinger equation under deterministic initial conditions

We investigate the semiclassical limit for the nonlinear Schrödinger equation in the case of a defocusing medium under oscillating nonperiodic initial conditions specified on the entire x axis. We formulate a system of integral conservation laws in terms of an infinite number of spatially averaged densities explicitly calculated from the initial conditions. We study the direct scattering problem and show that the scattering phase is a uniformly distributed random variable. The evolution of this system leads to the development of nonlinear oscillations, which become statistical in nature on long time scales. A modified inverse scattering method based on constructing a maximizer of the N-soliton solution in the continuum limit for N → is used to obtain an asymptotic solution. Using the maximizer, we found an infinite set of conserved averaged densities in the statistical state. This allowed us to couple the initial state with the limiting statistical steady (for t → ∞) state and, thus, to unambiguously determine the level spectrum in the statistical limit.     

Solution of the BGK model kinetic equation for very hard particle interaction

We study the Bhatnagar-Gross-Krook model kinetic equation with a velocity-dependent collision frequency. We derive the conditions that must be verified in order to keep the main physical properties of the Boltzmann equation, i.e.,H-theorem and conservation laws. The particular case of the so-called VHP interaction is considered, and the resulting kinetic equation is solved for a homogeneous and isotropic gas. Overpopulation phenomena are observed and analyzed for some kinds of initial conditions. The results are compared, where possible, with the exact solution of the Boltzmann equation.     

Mass properties of active and sterile neutrinos in a phenomenological (3 + 1 + 2) model

A generalized model involving three active neutrinos and three sterile neutrinos of different mass, one being relatively heavy [(3 + 1 + 2) model], is considered on the basis of experimental data, which admit the existence of anomalies beyond the minimally extended standard model featuring three active neutrinos of different mass. Basic properties used to describe massive active and sterile neutrinos are studied along with methods for determining the absolute scale of neutrino masses and for estimating neutrino masses on the basis of available experimental data. In the approximation of CP conservation, admissible values of the elements of the neutrino mass matrix are found from numerical calculations versus the possible values of the mass of one of the sterile neutrinos. The dependences of the mass properties of the neutrinos on the sterile-neutrino mass are constructed with allowance for possible sterile-neutrino contributions. The respective results can be used to interpret and predict results of various neutrino experiments.     

Properties of hydrodynamic nuclei for quantum field models

The transition to the hydrodynamic limit for a many-particle quantum system withN local conservation laws is made in a specified class of external effects. It is shown that the hydrodynamic equations are nonlocal in time and space and the hydrodynamic model is equivalent to the initial quantum statistical model. The nuclei appearing in the material relations are expressed in terms of the Green functions for the currents. It follows from the properties of the Green functions that the hydrodynamic model satisfies the dissipation conditions. When the quantum field model isT-invariant, the nuclei are related by reciprocal relations (analogous to the Onsager relations).Institute of Earth Physics, Russian Academy of Sciences. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 7, pp. 34–38, July, 1995.     

Condensation de l'acide sulfurique dans un conduit d'evacuation des gaz de combustion

A method for computing the condensation of water vapour and sulphuric acid in a removal gas conduct is proposed. It utilizes a theoretical approach to determine the condensate production during the phase change from the numerically computed thermal and dynamic properties of the steam in the conduct. The temperature and velocity distributions are given from a k — ε model in one-phase incompressible flow, taking into account the initial temperature and flow rate in the conduct. The quality of the model is tested on a full-scale experimentation pilot equipped with thermocouples and collectors of liquid condensate. The deposit production is obtained for different conditions of temperature and concentration, and the model proves to be satisfactory in domestic boiler conditions.     

Diffusion for a Quantum Particle Coupled to Phonons in d ≥ 3

We prove diffusion for a quantum particle coupled to a field of bosons (phonons or photons). The importance of this result lies in the fact that our model is fully Hamiltonian and randomness enters only via the initial (thermal) state of the bosons. This model is closely related to the one considered in (De Roeck, Fröhlich, Commun Math Phys 303:613–707, 2011) but various restrictive assumptions of the latter have been eliminated. In particular, depending on the dispersion relation of the bosons, the present result holds in dimension d ≥ 3 and no severe infrared conditions on the coupling are necessary.     

Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation

We present a class of augmented approximate Riemann solvers for the shallow water equations in the presence of a variable bottom surface. These belong to the class of simple approximate solvers that use a set of propagating jump discontinuities, or waves, to approximate the true Riemann solution. Typically, a simple solver for a system of m conservation laws uses m such discontinuities. We present a four wave solver for use with the the shallow water equations—a system of two equations in one dimension. The solver is based on a decomposition of an augmented solution vector—the depth, momentum as well as momentum flux and bottom surface. By decomposing these four variables into four waves the solver is endowed with several desirable properties simultaneously. This solver is well-balanced: it maintains a large class of steady states by the use of a properly defined steady state wave—a stationary jump discontinuity in the Riemann solution that acts as a source term. The form of this wave is introduced and described in detail. The solver also maintains depth non-negativity and extends naturally to Riemann problems with an initial dry state. These are important properties for applications with steady states and inundation, such as tsunami and flood modeling. Implementing the solver with LeVeque’s wave propagation algorithm [R.J. LeVeque, Wave propagation algorithms for multi-dimensional hyperbolic systems, J. Comput. Phys. 131 (1997) 327–335] is also described. Several numerical simulations are shown, including a test problem for tsunami modeling.     

Left-flat spaces and their associated space-times

It is already known that for an asymptotically flat space-time the metric coefficients and the other Newman-Penrose variables (in a suitable frame) can be constructed, in principle, by specifying certain initial data at conformal null infinity (and one further function on another null hypersurface), integrating the Newman-Penrose equations in the conformally rescaled “unphysical” space, and then transforming the results back to the physical space-time. If this is done approximately near ?⁺, for vacuum, the well-known Newman-Unti expansion is obtained. In this paper, after complexifying null infinity ?⁺ we generate, in a similar fashion, a left-flat spaceH using as much of the initial data of a given asymptotically flat space-timeM as possible, and show that the left-flat spaceH thus constructed is, in fact, the H-space corresponding toM. The advantage of our method is that it allows a reversal of procedure. Under suitable conditions we can generate from a given left-flat spaceH a class of physical space-times whose H-space is precisely the given left-flat spaceH. We shall see that the formal procedure requires only the local but not the global properties of ?⁺.     

Discretized kinetic theory on scale-free networks

The network of interpersonal connections is one of the possible heterogeneous factors which affect the income distribution emerging from micro-to-macro economic models. In this paper we equip our model discussed in [1, 2] with a network structure. The model is based on a system of n differential equations of the kinetic discretized-Boltzmann kind. The network structure is incorporated in a probabilistic way, through the introduction of a link density P(α) and of correlation coefficients P(β|α), which give the conditioned probability that an individual with α links is connected to one with β links. We study the properties of the equations and give analytical results concerning the existence, normalization and positivity of the solutions. For a fixed network with P(α) = c/α^q, we investigate numerically the dependence of the detailed and marginal equilibrium distributions on the initial conditions and on the exponent q. Our results are compatible with those obtained from the Bouchaud-Mezard model and from agent-based simulations, and provide additional information about the dependence of the individual income on the level of connectivity.     

A deterministic sandpile automaton revisited

The Bak-Tang-Wiesenfeld (BTW) sandpile model is a cellular automaton which has been intensively studied during the last years as a paradigm for self-organized criticality. In this paper, we reconsider a deterministic version of the BTW model introduced by Wiesenfeld, Theiler and McNamara, where sand grains are added always to one fixed site on the square lattice. Using the Abelian sandpile formalism we discuss the static properties of the system. We present numerical evidence that the deterministic model is only in the BTW universality class if the initial conditions and the geometric form of the boundaries do not respect the full symmetry of the square lattice. Received 19 August 1999     

Sensitivity to initial conditions, entropy production, and escape rate at the onset of chaos

We analytically link three properties of nonlinear dynamical systems, namely sensitivity to initial conditions, entropy production, and escape rate, in z-logistic maps for both positive and zero Lyapunov exponents. We unify these relations at chaos, where the Lyapunov exponent is positive, and at its onset, where it vanishes. Our result unifies, in particular, two already known cases, namely (i) the standard entropy rate in the presence of escape, valid for exponential functionality rates with strong chaos, and (ii) the Pesin-like identity with no escape, valid for the power-law behavior present at points such as the Feigenbaum one.     

Gauge-invariant lattice gas for the microcanonical Ising model

We introduce a lattice gas for particles with discrete momenta (1, 0, –1) and local deterministic microdynamics, which exactly reproduces Creutz's microcanonical algorithm for the ferromagnetic Ising model. However, because of the manifest gauge invariance of our variables, both the Ising ferromagnetic and spin-glass systems share precisely the same dynamics with different initial conditions. Additional conservation laws in the 1D Ising case result in a completely integrable system in the limit of zero or unbounded demon energy cutoff. Numerical investigations of ergodicity are presented for the pure Ising lattice gas in one and two dimensions.     

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.