François Dubois

François Dubois

We show that when we formulate the lattice Boltzmann equation with a small time step $Δ$t and an associated space scale $Δ$x, a Taylor expansion joined with the so-called equivalent equation methodology leads to establish macroscopic fluid equations as a formal limit. We recover the Euler equations of gas dynamics at the first order and the compressible Navier-Stokes equations at the second order. 1) Discrete geometry $\bullet$ We denote by d the dimension of space and by L a regular d-dimensional lattice. Such a lattice is composed by a set L 0 of nodes or vertices and a set L 1 of links or edges between two vertices. From a practical point of view, given a vertex x, there exists a set V (x) of neighbouring nodes, including the node x itself. We consider here that the lattice L is parametrized by a space step $Δ$x > 0. For the fundamental example called D2Q9 (see e.g. Lallemand and Luo, 2000), the set V (x) is given with the help of the family of vectors (e j) 0$\le$j$\le$J defined by J = 8, (1.1) e j = 0 0 , 1 0 , 0 1 , --1 0 , 0 --1 , 1 1 , --1 1 , --1 --1 , 1 --1 and the vicinity (1.2) V (x) = { x + $Δ$x e j , 0 $\le$ j $\le$ J } .

Pierre Lallemand, François Dubois

We show that a single particle distribution for the D2Q13 lattice Boltzmann scheme can simulate coupled effects involving advection and diffusion of velocity and temperature. We consider various test cases: non-linear waves with periodic boundary conditions, a test case with buoyancy, propagation of transverse waves, Couette and Poiseuille flows. We test various boundary conditions and propose to mix bounce-back and anti-bounce-back numerical boundary conditions to take into account velocity and temperature Dirichlet conditions. We present also first results for the de Vahl Davis heated cavity. Our results are compared with the coupled D2Q9-D2Q5 lattice Boltzmann approach for the Boussinesq system and with an elementary finite differences solver for the compressible Navier-Stokes equations.

François Dubois, Pierre Lallemand, Christian Obrecht et al.

We show that the asymptotic properties of the link-wise artificial compressibility method are not compatible with a correct approximation of fluid properties. We propose to adapt the previous method through a framework suggested by the Taylor expansion method and to replace first order terms in the expansion by appropriate three or five points finite differences and to add non linear terms. The "FD-LBM" scheme obtained by this method is tested in two dimensions for shear wave, Stokes modes and Poiseuille flow. The results are compared with the usual lattice Boltzmann method in the framework of multiple relaxation times.

Bruce Boghosian, François Dubois, Benjamin Graille et al.

We consider the D1Q3 lattice Boltzmann scheme with an acoustic scale for the simulation of diffusive processes. When the mesh is refined while holding the diffusivity constant, we first obtain asymptotic convergence. When the mesh size tends to zero, however, this convergence breaks down in a curious fashion, and we observe qualitative discrepancies from analytical solutions of the heat equation. In this work, a new asymptotic analysis is derived to explain this phenomenon using the Taylor expansion method, and a partial differential equation of acoustic type is obtained in the asymptotic limit. We show that the error between the D1Q3 numerical solution and a finite-difference approximation of this acoustic-type partial differential equation tends to zero in the asymptotic limit. In addition, a wave vector analysis of this asymptotic regime demonstrates that the dispersion equation has nontrivial complex eigenvalues, a sign of underlying propagation phenomena, and a portent of the unusual convergence properties mentioned above.

François Dubois, Benjamin Graille, Pierre Lallemand

We consider multi relaxation times lattice Boltzmann scheme with two particle distributions for the thermal Navier Stokes equations formulated with conservation of mass and momentum and dissipation of volumic entropy.Linear stability is taken into consideration to determine a coupling between two coefficients of dissipation.We present interesting numerical results for one-dimensional strong nonlinear acoustic waves with shocks.

François Dubois

We follow the mathematical framework proposed by Bouchut and present in this contribution a dual entropy approach for determining equilibrium states of a lattice Boltzmann scheme. This method is expressed in terms of the dual of the mathematical entropy relative to the underlying conservation law. It appears as a good mathematical framework for establishing a "H-theorem" for the system of equations with discrete velocities. The dual entropy approach is used with D1Q3 lattice Boltzmann schemes for the Burgers equation. It conducts to the explicitation of three different equilibrium distributions of particles and induces naturally a nonlinear stability condition. Satisfactory numerical results for strong nonlinear shocks and rarefactions are presented. We prove also that the dual entropy approach can be applied with a D1Q3 lattice Boltzmann scheme for systems of linear and nonlinear acoustics and we present a numerical result with strong nonlinear waves for nonlinear acoustics. We establish also a negative result: with the present framework, the dual entropy approach cannot be used for the shallow water equations.

François Dubois, Tony Fevrier, Benjamin Graille

This paper studies the stability properties of a two dimensional relative velocity scheme for the Navier-Stokes equations. This scheme inspired by the cascaded scheme has the particularity to relax in a frame moving with a velocity field function of space and time. Its stability is studied first in a linear context then on the non linear test case of the Kelvin-Helmholtz instability. The link with the choice of the moments is put in evidence. The set of moments of the cascaded scheme improves the stability of the d'Humières scheme for small viscosities. On the contrary, a relative velocity scheme with the usual set of moments deteriorates the stability.

François Dubois, Tony Fevrier, Benjamin Graille

In this contribution, a new class of lattice Boltzmann schemes is introduced and studied. These schemes are presented in a framework that generalizes the multiple relaxation times method of d'Humières. They extend also the Geier's cascaded method. The relaxation phase takes place in a moving frame involving a set of moments depending on a given relative velocity field. We establish with the Taylor expansion method that the equivalent partial differential equations are identical to the ones obtained with the multiple relaxation times method up to the second order accuracy. The method is then performed to derive the equivalent equations up to third order accuracy.

François Dubois, Abdelkader Saïdi

We present a numerical scheme for the resolution of matrix Riccati equation used in control problems. The scheme is unconditionnally stable and the solution is definite positive at each time step of the resolution. We prove the convergence in the scalar case and present several numerical experiments for classical test cases.

François Dubois

We introduce in this contribution the notion of partial Riemann problem. Recall that the Riemann problem describes a shock tube interaction between two given states ; the partial Riemann problem is a generalization of the previous concept and introduces the notion of boundary manifold. In what follows, we first recall very classical notions concerning gas dynamics and the associated Riemann problem. In a second part, we introduce the partial Riemann problem for general systems of conservation laws and proves that this problem admits a solution in some class of appropriate nonlinear waves. In section 3, we recall the linearized analysis with the method of characteristics, introduce the weak formulation of the Dirichlet boundary condition for nonlinear situations in terms of the partial Riemann problem and show that lot of physically relevant situations are described with this theoretical framework. In the last paragraph, we propose a practical implementation of the previous onsiderations with the finite volume method.

François Dubois

The Van Leer approach for the approximation of nonlinear scalar conservation laws is studied in one space dimension. The problem can be reduced to a nonlinear interpolation and we propose a convexity property for the interpolated values. We prove that under general hypotheses the method of lines in well posed in $\ ell^{\infty} \cap {\rm BV} $ and we give precise sufficient conditions to establish that the total variation is diminishing. We observe that the second order accuracy can be maintained even at non sonic extrema. We establish also that both the TVD property and second order accuracy can be maintained after discretization in time with the second order accurate Heun scheme. Numerical illustration for the advection equation is presented.

Adeline Augier, François Dubois, Benjamin Graille

In this paper, we investigate the numerous parameters choices for linear lattice Boltzmann schemes according to the definition of the isotropic order given in \cite{ADG11}. This property---written in a general framework including all of the \ddqq schemes---can be read through a group operation. It implies some relations on the parameters of the scheme (equilibrium states and relaxation times) that give rigorous methodology to select them according to the desired order of isotropy. For acoustic applications in two spaces dimensions (namely \ddqn and \ddqt schemes) this methodology is used to propose a full description of the sets of parameters that involve isotropy of order $m$ ($m\in\{1,2,3,5\}$ for \ddqn and $m\in\{1,2\}$ for \ddqt). We then propose numerical illustrations for the \ddqn scheme.

François Dubois, Pierre Lallemand

In this contribution we extend the Taylor expansion method proposed previously by one of us and establish equivalent partial differential equations of DDH lattice Boltzmann scheme at an arbitrary order of accuracy. We derive formally the associated dynamical equations for classical thermal and linear fluid models in one to three space dimensions. We use this approach to adjust relaxation parameters in order to enforce fourth order accuracy for thermal model and diffusive relaxation modes of the Stokes problem. We apply the resulting scheme for numerical computation of associated eigenmodes and compare our results with analytical references.

François Dubois, Hervé Le Meur, Claude Reiss

This contribution is devoted to a new model of HIV multiplication motivated by the patent of one of the authors. We take into account the antigenic diversity through what we define "antigenicity", whether of the virus or of the adapted lymphocytes. We model the interaction of the immune system and the viral strains by two processes. On the one hand, the presence of a given viral quasi-species generates antigenically adapted lymphocytes. On the other hand, the lymphocytes kill only viruses for which they have been designed. We consider also the mutation and multiplication of the virus. An original infection term is derived. So as to compare our system of differential equations with well-known models, we study some of them and compare their predictions to ours in the reduced case of only one antigenicity. In this particular case, our model does not yield any major qualitative difference. We prove mathematically that, in this case, our model is biologically consistent (positive fields) and has a unique continuous solution for long time evolution. In conclusion, this model improves the ability to simulate more advanced phases of the disease.

Adeline Augier, François Dubois, Benjamin Graille

In this paper, we recall the linear version of the lattice Boltzmann schemes in the framework proposed by d'Humiéres. According to the equivalent equations we introduce a definition for a scheme to be isotropic at some order. This definition is chosen such that the equivalent equations are preserved by orthogonal transformations of the frame. The property of isotropy can be read through a group operation and then implies a sequence of relations on relaxation times and equilibrium states that characterizes a lattice Boltzmann scheme. We propose a method to select the parameters of the scheme according to the desired order of isotropy. Applying it to the D2Q9 scheme yields the classical constraints for the first and second orders and some non classical for the third and fourth orders.

François Dubois, Isabelle Greff, Charles Pierre

The mixed finite element method for the Poisson problem with the Raviart-Thomas elements of low-level can be interpreted as a finite volume method with a non-local gradient. In this contribution, we propose a variant of Petrov-Galerkin type for this problem to ensure a local computation of the gradient at the interfaces of the elements. The shape functions are the Raviart-Thomas finite elements. Our goal is to define test functions that are in duality with these shape functions: Precisely, the shape and test functions will be asked to satisfy a L2-orthogonality property. The general theory of Babuška brings necessary and sufficient stability conditions for a Petrov-Galerkin mixed problem to be convergent. We propose specific constraints for the dual test functions in order to ensure stability. With this choice, we prove that the mixed Petrov-Galerkin scheme is identical to the four point finite volumes scheme of Herbin, and to the mass lumping approach developed by Baranger, Maitre and Oudin. Finally, we construct a family of dual test functions that satisfy the stability conditions. Convergence is proven with the usual techniques of mixed finite elements.

François Dubois, Eric Duceau, Frédéric Marechal et al.

We study acoustic wave propagation in a uniform stationary flow. We develop a method founded on the Lorentz transform and a hypothesis of irrotationality of the acoustic perturbation. After a transformation of the space-time and of the unknown fields, we derive a system of partial differential equations that eliminates the external flow and deals with the classical case of non advective acoustics. A sequel of the analysis is a new set of perfectly matched layers equations in the spirit of the work of Berenger and Collino. The numerical implementation of the previous ideas is presented with the finite differences method HaWAY on cartesian staggered grids. Relevant numerical tests are proposed.

François Dubois

For an elliptic problem with two space dimensions, we propose to formulate the finite volume method with the help of Petrov-Galerkin mixed finite elementsthat are based on the building of a dual Raviart-Thomas basis.

François Dubois

We propose a notion of hyperbolic system of conservation laws invariant for the Galileo group of transformations. We show that with natural physical and mathematical hypotheses, such a system conducts to the gas dynamics equations or to exotic systems that are detailed in this contribution to Cemracs 99.

François Dubois, Abdelkader Saïdi

In this paper we present a numerical scheme for the resolution of matrix Riccati equation, usualy used in control problems. The scheme is unconditionnaly stable and the solution is definite positive at each time step of the resolution. We prove the convergence in the scalar case and present several numerical experiments for classical test cases.

François Dubois

For an elliptic problem with two space dimensions, we propose to formulate the finite volume method with the help of Petrov-Galerkin mixed finite elementsthat are based on the building of a dual Raviart-Thomas basis.

Benjamin Graille, François Dubois, Tony Fevrier

We study the formal precision of the relative velocity lattice Boltzmann schemes. They differ from the d'Humières schemes by their relaxation phase: it occurs for a set of moments parametrized by a velocity field function of space and time. We deal with the asymptotics of the relative velocity schemes for one conservation law: the third order equivalent equation is exposed for an arbitrary number of dimensions and velocities.

François Dubois, Isabelle Greff, Charles Pierre

The general theory of Babuška ensures necessary and sufficient conditions for a mixed problem in classical or Petrov-Galerkin form to be well posed in the sense of Hadamard. Moreover, the mixed method of Raviart-Thomas with low-level elements can be interpreted as a finite volume method with a non-local gradient. In this contribution, we propose a variant of type Petrov-Galerkin to ensure a local computation of the gradient at the interfaces of the elements. The in-depth study of stability leads to a specific choice of the test functions. With this choice, we show on the one hand that the mixed Petrov-Galerkin obtained is identical to the finite volumes scheme "volumes finis à 4 points" ("VF4") of Faille, Galloüet and Herbin and to the condensation of mass approach developed by Baranger, Maitre and Oudin. On the other hand, we show the stability via an inf-sup condition and finally the convergence with the usual methods of mixed finite elements.

François Dubois

We propose an elementary introduction to the finite volume method in the context of gas dynamics conservation laws. Our approach is founded on the advection equation, the exact integration of the associated Cauchy problem, and the so-called upwind scheme in one space dimension. It is then extended in three directions: hyperbolic linear systems and particularily the system of acoustics, gas dynamics with the help of the Roe matrix and two space dimensions by following the approach proposed by Van Leer. A special emphasis on boundary conditions is proposed all along the text.

Régis Msallam, François Dubois

Nonlinear acoustics of wind instruments conducts to study unidimensional fluid flows. From physically relevant approximations that are modelized with the thin layer Navier Stokes equations, we propose a coupled model where perfect fluid flow is described by the Euler equations of gas dynamics and viscous and thermal boundary layer is modelized by a linear equation. We describe numerical discretization, validate the associated software by comparison with analytical solutions and consider musical application of strongly nonlinear waves in the trombone.

Zeno Toffano, François Dubois

This work proposes a formulation of propositional logic, named Eigenlogic, using quantum observables as propositions. The eigenvalues of these operators are the truth-values and the associated eigenvectors the interpretations of the propositional system. Fuzzy logic arises naturally when considering vectors outside the eigensystem, the fuzzy membership function is obtained by the Born rule of the logical observable.This approach is then applied in the context of quantum robots using simple behavioral agents represented by Braitenberg vehicles. Processing with non-classical logic such as multivalued logic, fuzzy logic and the quantum Eigenlogic permits to enlarge the behavior possibilities and the associated decisions of these simple agents.

François Dubois, Tony Février, Benjamin Graille

In this contribution, we study the theoretical and numerical stability of a bidimensional relative velocity lattice Boltzmann scheme. These relative velocity schemes introduce a velocity field parameter called "relative velocity" function of space and time. They generalize the d'Humières multiple relaxation times scheme and the cascaded automaton. This contribution studies the stability of a four velocities scheme applied to a single linear advection equation according to the value of this relative velocity. We especially compare when it is equal to 0 (multiple relaxation times scheme) or to the advection velocity ("cascaded like" scheme). The comparison is made in terms of L1 and L2 stability. The L1 stability area is fully described in terms of relaxation parameters and advection velocity for the two choices of relative velocity. These results establish that no hierarchy of these two choices exists for the L1 notion. Instead, choosing the parameter equal to the advection velocity improves the numerical L2 stability of the scheme. This choice cancels some dispersive terms and improve the numerical stability on a representative test case. We theoretically strengthen these results with a weighted L2 notion of stability.

Olivier Maurice, Alain Reineix, Philippe Durand et al.

Starting from topological principles we first recall the elementary ones giving Kirchhoff's laws for current conservation. Using in a second step the properties of spaning tree, we show that currents are under one hypothesis intrinsically boundaries of surfaces flux. Naturally flux appears as the object from which the edge comes from. The current becomes the magnetomotive force (mmf) that creates the flux in the magnetostatic representation. Using a metric and an Hodge's operator, this flux creates an electromotive force (emf). This emf is finally linked with the current to give the fundamental tensor - or "metric" - of the Kron's tensorial analysis of networks. As it results in a link between currents of cycles (surface boundaries) and energy sources in the network, we propose to symbolize this cross talk using chords between cycles in the graph structure on which the topology is based. Starting then from energies relations we show that this metric is the Lagrange's operator of the circuit. But introducing moment space, the previous results can be extended to non local interactions as far field one. And to conclude, we use the same principle to create general relation of information exchange between networks as functors between categories.

Nolwenn Balin, Fabien Casenave, François Dubois et al.

We consider the scattering of acoustic perturbations in a presence of a flow. We suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential. In the first zone, we apply a Prandtl-Glauert transformation to recover the Helmholtz equation. The well-known setting of boundary element method for the Helmholtz equation is available. In the second zone, the flow quantities are space dependent, we have to consider a local resolution, namely the finite element method. Herein, we carry out the coupling of these two methods and present various applications and validation test cases. The source term is given through the decomposition of an incident acoustic field on a section of the computational domain's boundary.

François Dubois, Pierre Lallemand

Using the Taylor expansion method, we show that it is possible to improve the lattice Boltzmann method for acoustic applications. We derive a formal expansion of the eigenvalues of the discrete approximation and fit the parameters of the scheme to enforce fourth order accuracy. The corresponding discrete equations are solved with the help of symbolic manipulation. The solutions are explicited in the case of D3Q27 lattice Boltzmann scheme. Various numerical tests support the coherence of this approach.