Source author record

François Dubois

François Dubois appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

On hyperbolic systems with entropy velocity covariant under the action of a group

For hyperbolic systems of conservation laws in one space dimension with a mathematical entropy, we define the notion of entropy velocity. Then we give sufficient conditions for such a system to be covariant under the action of a group of space-time transformations. These conditions naturally introduce a representation of the group in the space of states. We construct such hyperbolic system from the knowledge of data on the manifold of null velocity. We apply these ideas for Galileo, Lorentz and circular groups. We focus on particular non trivial examples for two by two systems of conservation laws.

preprint2022arXiv

Pontryagin calculus in Riemannian geometry

In this contribution, we study systems with a finite number of degrees of freedom as in robotics. A key idea is to consider the mass tensor associated to the kinetic energy as a metric in a Riemannian configuration space. We apply Pontryagin's framework to derive an optimal evolution of the control forces and torques applied to the mechanical system. This equation under covariant form uses explicitly the Riemann curvature tensor. This contribution is dedicated to the memory of Claude Vall{é}e (1945-2014).

preprint2020arXiv

On anti bounce back boundary condition for lattice Boltzmann schemes

In this contribution, we recall the derivation of the anti bounce back boundary condition for the D2Q9 lattice Boltzmann scheme. We recall various elements of the state of the art for anti bounce back applied to linear heat and acoustics equations and in particular the possibility to take into account curved boundaries. We present an asymptotic analysis that allows an expansion of all the fields in the boundary cells. This analysis based on the Taylor expansion method confirms the well known behaviour of anti bounce back boundary for the heat equation. The analysis puts also in evidence a hidden differential boundary condition in the case of linear acoustics. Indeed, we observe discrepancies in the first layers near the boundary. To reduce these discrepancies, we propose a new boundary condition mixing bounce back for the oblique links and anti bounce back for the normal link. This boundary condition is able to enforce both pressure and tangential velocity on the boundary. Numerical tests for the Poiseuille flow illustrate our theoretical analysis and show improvements in the quality of the flow.

preprint2020arXiv

Unexpected convergence of lattice Boltzmann schemes

In this work, we study numerically the convergence of the scalar D2Q9 lattice Boltzmann scheme with multiple relaxation times when the time step is proportional to the space step and tends to zero. We do this by a combination of theory and numerical experiment. The classical formal analysis when all the relaxation parameters are fixed and the time step tends to zero shows that the numerical solution converges to solutions of the heat equation, with a constraint connecting the diffusivity, the space step and the coefficient of relaxation of the momentum. If the diffusivity is fixed and the space step tends to zero, the relaxation parameter for the momentum is very small, causing a discrepency between the previous analysis and the numerical results. We propose a new analysis of the method for this specific situation of evanescent relaxation, based on the dispersion equation of the lattice Boltzmann scheme. A new asymptotic partial differential equation, the damped acoustic system, is emergent as a result of this formal analysis. Complementary numerical experiments establish the convergence of the scalar D2Q9 lattice Boltzmann scheme with multiple relaxation times and acoustic scaling in this specific case of evanescent relaxation towards the numerical solution of the damped acoustic system.

preprint2016arXiv

Coupling Linear Sloshing with Six Degrees of Freedom Rigid Body Dynamics

Fluid motion in tanks is usually described in space industry with the so-called Lomen hypothesis which assumes the vorticity is null in the moving frame. We establish in this contribution that this hypothesis is valid only for uniform rotational motions. We give a more general formulation of this coupling problem, with a compact formulation. We consider the mechanical modeling of a rigid body with a motion of small amplitude, containing an incompressible fluid in the linearized regime. We first establish that the fluid motion remains irrotational in a Galilean referential if it is true at the initial time. When continuity of normal velocity and pressure are prescribed on the free surface, we establish that the global coupled problem conserves an energy functional composed by three terms. We introduce the Stokes - Zhukovsky vector fields, solving Neumann problems for the Laplace operator in the fluid in order to represent the rotational rigid motion with irrotational vector fields. Then we have a good framework to consider the coupled problem between the fluid and the rigid motion. The coupling between the free surface and the ad hoc component of the velocity potential introduces a "Neumann to Dirichlet" operator that allows to write the coupled system in a very compact form. The final expression of a Lagrangian for the coupled system is derived and the Euler-Lagrange equations of the coupled motion are presented.

preprint2015arXiv

Lattice Boltzmann schemes with relative velocities

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.

preprint2015arXiv

On the stability of a relative velocity lattice Boltzmann scheme for compressible Navier-Stokes equations

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.

preprint2015arXiv

Stability of a bidimensional relative velocity lattice Boltzmann scheme

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.

preprint2015arXiv

Third order equivalent equation for the relative velocity lattice Boltzmann schemes with one conservation law

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.

preprint2014arXiv

Boundary Element and Finite Element Coupling for Aeroacoustics Simulations

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.

preprint2014arXiv

Finite volumes and mixed Petrov-Galerkin finite elements : the unidimensional problem

For Laplace operator in one space dimension, we propose to formulate the heuristic finite volume method with the help of mixed Petrov-Galerkin finite elements. Weighting functions for gradient discretization are parameterized by some universal function. We propose for this function a compatibility interpolation condition and we prove that such a condition is equivalent to the inf-sup property when studying stability of the numerical scheme. In the case of stable scheme and under two distinct hypotheses concerning the regularity of the solution, we demonstrate convergence of the finite volume method in appropriate Hilbert spaces and with optimal order of accuracy.

preprint2014arXiv

Kron's method and cell complexes for magnetomotive and electromotive forces

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.

preprint2014arXiv

On a superconvergent lattice Boltzmann boundary scheme

In a seminal paper Ginzburg and Adler analyzed the bounce-back boundary conditions for the lattice Boltzmann scheme and showed that it could be made exact to second order for the Poiseuille flow if some expressions depending upon the parameters of the method were satisfied, thus defining so-called "magic parameters". Using the Taylor expansion method that one of us developed, we analyze a series of simple situations (1D and 2D) for diffusion and for linear fluid problems using bounce-back and "anti bounce-back" numerical boundary conditions. The result is that "magic parameters" depend upon the detailed choice of the moments and of their equilibrium values. They may also depend upon the way the flow is driven.

preprint2014arXiv

On rotational invariance of lattice Boltzmann schemes

We propose the derivation of acoustic-type isotropic partial differential equations that are equivalent to linear lattice Boltzmann schemes with a density scalar field and a momentum vector field as conserved moments. The corresponding linear equivalent partial differential equations are generated with a new "Berliner version" of the Taylor expansion method. The details of the implementation are presented. These ideas are applied for the D2Q9, D2Q13, D3Q19 and D3Q27 lattice Boltzmann schemes. Some limitations associated with necessary stability conditions are also presented.

preprint2014arXiv

On triangular lattice Boltzmann schemes for scalar problems

We propose to extend the d'Humiéres version of the lattice Boltzmann scheme to triangular meshes. We use Bravais lattices or more general lattices with the property that the degree of each internal vertex is supposed to be constant. On such meshes, it is possible to define the lattice Boltzmann scheme as a discrete particle method, without need of finite volume formulation or Delaunay-Voronoi hypothesis for the lattice. We test this idea for the heat equation and perform an asymptotic analysis with the Taylor expansion method for two schemes named D2T4 and D2T7. The results show a convergence up to second order accuracy and set new questions concerning a possible super-convergence.

preprint2014arXiv

Simulation of strong nonlinear waves with vectorial lattice Boltzmann schemes

We show that an hyperbolic system with a mathematical entropy can be discretized with vectorial lattice Boltzmann schemes with the methodology of kinetic representation of the dual entropy. We test this approach for the shallow water equations in one and two space dimensions. We obtain interesting results for a shock tube, reflection of a shock wave and unstationary two-dimensional propagation. This contribution shows the ability of vectorial lattice Boltzmann schemes to simulate strong nonlinear waves in unstationary situations.

preprint2013arXiv

Some results on energy-conserving lattice Boltzmann models

We consider the problem of ''energy conserving'' lattice Boltzmann models. A major difficulty observed in previous studies is the coupling between the viscous and thermal waves even at moderate wave numbers. We propose a theoretical framework based founded on the knowledge of the partial equivalent equations of the lattice Boltzmann scheme at several orders of precision. With the help of linearized models (inviscid and dissipative advective acoustics and classical acoustics), we suggest natural sets of relations for the parameters of lattice Boltzmann schemes. The application is proposed for three two-dimensional schemes. Numerical test cases for simple linear and nonlinear waves establish that the main difficulty in the previous contributions can now be overcome.

preprint2013arXiv

Stable lattice Boltzmann schemes with a dual entropy approach for monodimensional nonlinear waves

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.

preprint2011arXiv

An introduction to finite volumes for gas dynamics

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.

preprint2011arXiv

Isotropy conditions for lattice Boltzmann schemes. Application to D2Q9

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.

preprint2011arXiv

Linear Lattice Boltzmann Schemes for Acoustic: parameters choices and isotropy properties

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.

preprint2011arXiv

Lorentz Transform and Staggered Finite Differences for Advective Acoustics

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.

preprint2011arXiv

Mathematical modeling of antigenicity for HIV dynamics

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.

preprint2011arXiv

Partial Riemann problem, boundary conditions, and gas dynamics

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.

preprint2011arXiv

Quartic Parameters for Acoustic Applications of Lattice Boltzmann Scheme

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.

preprint2010arXiv

Mathematical model for coupling a quasi-unidimensional perfect flow

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.

preprint2010arXiv

Nonlinear Interpolation and Total Variation Diminishing Schemes

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.

preprint2009arXiv

Towards higher order lattice Boltzmann schemes

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.