A network programming approach in solving Darcy's equations by mixed finite-element methods.
In this work we design a new domain decomposition method for the Euler equations in 2 dimensions. The starting point is the equivalence with a third order scalar equation to whom we can apply an algorithm inspired from the Robin-Robin preconditioner for the convection-diffusion equation [Achdou and Nataf, C. R. Acad. Sci. Paris Sér. I325 (1997) 1211–1216]. Afterwards we translate it into an algorithm for the initial system and prove that at the continuous level and for a decomposition into 2 sub-domains,...
A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations is presented. Applying the orthogonal projection technique, we introduce two local Gauss integrations as a stabilizing term in the error correction method, and derive a new error correction method. In both the coarse solution computation step and the error computation step, a locally stabilizing term based on two local Gauss integrations is introduced. The stability and convergence of the...
We analyze a new formulation of the Stokes equations in three-dimensional axisymmetric geometries, relying on Fourier expansion with respect to the angular variable: the problem for each Fourier coefficient is two-dimensional and has six scalar unknowns, corresponding to the vector potential and the vorticity. A spectral discretization is built on this formulation, which leads to an exactly divergence-free discrete velocity. We prove optimal error estimates.
We analyze a new formulation of the Stokes equations in three-dimensional axisymmetric geometries, relying on Fourier expansion with respect to the angular variable: the problem for each Fourier coefficient is two-dimensional and has six scalar unknowns, corresponding to the vector potential and the vorticity. A spectral discretization is built on this formulation, which leads to an exactly divergence-free discrete velocity. We prove optimal error estimates.
We introduce a new stable MINI-element pair for incompressible Stokes equations on quadrilateral meshes, which uses the smallest number of bubbles for the velocity. The pressure is discretized with the P1-midpoint-edge-continuous elements and each component of the velocity field is done with the standard Q1-conforming elements enriched by one bubble a quadrilateral. The superconvergence in the pressure of the proposed pair is analyzed on uniform rectangular meshes, and tested numerically on uniform...
We study a model of motion of a passive tracer particle in a turbulent flow that is strongly mixing in time variable. In [8] we have shown that there exists a probability measure equivalent to the underlying physical probability under which the quasi-Lagrangian velocity process, i.e. the velocity of the flow observed from the vintage point of the moving particle, is stationary and ergodic. As a consequence, we proved the existence of the mean of the quasi-Lagrangian velocity, the so-called Stokes...
There are many problems of groundwater flow in a disrupted rock massifs that should be modelled using numerical models. It can be done via “standard approaches” such as increase of the permeability of the porous medium to account the fracture system (or double-porosity models), or discrete stochastic fracture network models. Both of these approaches appear to have their constraints and limitations, which make them unsuitable for the large- scale long-time hydrogeological calculations. In the article,...
A high resolution finite volume method for the computation of unsteady solutions of the Euler equations in two space dimensions is presented and validated. The scheme is of Godunov-type. The first order part of the flux function uses the approximate Riemann problem solver of Pandolfi and here a new derivation of this solver is presented. This construction paves the way to understand the conditions under which the scheme satisfies an entropy condition. The extension to higher order is done by applying...
The results of a workshop concerning the numerical simulation of the liquid flow around a hydrofoil in non-cavitating and cavitating conditions are presented. This workshop was part of the conference “Mathematical and Numerical aspects of Low Mach Number Flows” (2004) and was aimed to investigate the capabilities of different compressible flow solvers for the low Mach number regime and for flows in which incompressible and supersonic regions are simultaneously present. Different physical models...
The results of a workshop concerning the numerical simulation of the liquid flow around a hydrofoil in non-cavitating and cavitating conditions are presented. This workshop was part of the conference “Mathematical and Numerical aspects of Low Mach Number Flows” (2004) and was aimed to investigate the capabilities of different compressible flow solvers for the low Mach number regime and for flows in which incompressible and supersonic regions are simultaneously present. Different physical models...
The penalty method when applied to the Stokes problem provides a very efficient algorithm for solving any discretization of this problem since it gives rise to a system of two equations where the unknowns are uncoupled. For a spectral or spectral element discretization of the Stokes problem, we prove a posteriori estimates that allow us to optimize the penalty parameter as a function of the discretization parameter. Numerical experiments confirm the interest of this technique.
The penalty method when applied to the Stokes problem provides a very efficient algorithm for solving any discretization of this problem since it gives rise to a system of two equations where the unknowns are uncoupled. For a spectral or spectral element discretization of the Stokes problem, we prove a posteriori estimates that allow us to optimize the penalty parameter as a function of the discretization parameter. Numerical experiments confirm the interest of this technique.
A general construction of test functions in the Petrov-Galerkin method is described. Using this construction; algorithms for an approximate solution of the Dirichlet problem for the differential equation are presented and analyzed theoretically. The positive number is supposed to be much less than the discretization step and the values of . An algorithm for the corresponding two-dimensional problem is also suggested and results of numerical tests are introduced.