In this paper, we design and analyze a Hybrid High-Order discretization method for the steady motion of non-Newtonian, incompressible fluids in the Stokes approximation of small velocities. The proposed method has several appealing features including the support of general meshes and high-order, unconditional inf-sup stability, and orders of convergence that match those obtained for scalar Leray-Lions problems. A complete well-posedness and convergence analysis of the method is carried out under new, general assumptions on the strain rate-shear stress law, which encompass several common examples such as the power-law and Carreau-Yasuda models. Numerical examples complete the exposition.
In this paper we present a rigorous derivation of the reduced MHD models with and without parallel velocity that are implemented in the non-linear MHD code JOREK. The model we obtain contains some terms that have been neglected in the implementation but might be relevant in the nonlinear phase. These are necessary to guarantee exact conservation with respect to the full MHD energy. For the second part of this work, we have replaced the linearized time stepping of JOREK by a nonlinear solver based on the Inexact Newton method including adaptive time stepping. We demonstrate that this approach is more robust especially with respect to numerical errors in the saturation phase of an instability and allows to use larger time steps in the non-linear phase.
We study here the water waves problem for uneven bottoms in a highly nonlinear regime where the small amplitude assumption of the Korteweg-de Vries (KdV) equation is enforced. It is known that, for such regimes, a generalization of the KdV equation (somehow linked to the Camassa-Holm equation) can be derived and justified [Constantin and Lannes, Arch. Ration. Mech. Anal. 192 (2009) 165–186] when the bottom is flat. We generalize here this result with a new class of equations taking into account variable bottom topographies. Of course, many variable depth KdV equations existing in the literature are recovered as particular cases. Various regimes for the topography regimes are investigated and we prove consistency of these models, as well as a full justification for some of them. We also study the problem of wave breaking for our new variable depth and highly nonlinear generalizations of the KdV equations.
We focus here on the water waves problem for uneven bottoms in the long-wave regime, on an unbounded two or three-dimensional domain. In order to derive asymptotic models for this problem, we consider two different regimes of bottom topography, one for small variations in amplitude, and one for strong variations. Starting from the Zakharov formulation of this problem, we rigorously compute the asymptotics expansion of the involved Dirichlet-Neumann operator. then, following the global strategy introduced by Bona, Colin and Lannes, new symetric asymptotic models are derived for each regime of bottom topography. Solutions of these systems are proved to give good approximations of solutions of the water waves problem. These results hold for solutions that evanesce at infinity as well as for spatially periodic ones.
In order to handle the flow of a viscous incompressible fluid in a porous medium with cracks, the thickness of which cannot be neglected, we consider a model which couples the Darcy equations in the medium with the Stokes equations in the cracks by a new boundary condition at the interface, namely the continuity of the pressure. We prove that this model admits a unique solution and propose a mixed formulation of it. Relying on this formulation, we describe a finite element discretization and derive a priori and a posteriori error estimates. We present some numerical experiments that are in good agreement with the analysis.
The purpose of this work is to study an example of low Mach (Froude) number limit of compressible flows when the initial density (height) is almost equal to a function depending on x. This allows us to connect the viscous shallow water equation and the viscous lake equations. More precisely, we study this asymptotic with well prepared data in a periodic domain looking at the influence of the variability of the depth. The result concerns weak solutions. In a second part, we discuss the general low Mach number limit for standard compressible flows given in P.–L. Lions' book that means with constant viscosity coefficients.
In order to describe a solid which deforms smoothly in some region, but non smoothly in some other region, many multiscale methods have recently been proposed.They aim at coupling an atomistic model (discrete mechanics) with a macroscopic model (continuum mechanics).We provide here a theoretical ground for such a coupling in a one-dimensional setting.We briefly study the general case of a convex energy, and next concentrate on a specific example of a nonconvex energy, the Lennard-Jones case.In the latter situation, we prove that the discretization needs to account in an adequate way for the coexistence of a discrete model and a continuous one.Otherwise, spurious discretization effects may appear.We provide a numerical analysis of the approach.
In this paper we investigate the motion of a rigid ball in an incompressible perfect fluid occupying R-2. We prove the global in time existence and the uniqueness of the classical solution for this fluid-structure problem. The proof relies mainly on weighted estimates for the vorticity associated with the strong solution of a fluid-structure problem obtained by incorporating some dissipation.
Many problems in quantum chemistry deal with the computation of fundamental or excited states of molecules and lead to the resolution of eigenvalue problems. One of the major difficulties in these computations lies in the very large dimension of the systems to be solved. Indeed these eigenfunctions depend on 3n variables where n stands for the number of particles (electrons and/or nucleari) in the molecule. In order to diminish the size of the systems to be solved, the chemists have proposed many interesting ideas. Among those stands the adiabatic variable method; we present in this paper a mathematical analysis of this approximation and propose, in particular, an a posteriori estimate that might allow for verifying the adiabaticity hypothesis that is done on some variables; numerical simulations that support the a posteriori estimators obtained theoretically are also presented.
We present a finite element method to compute guided modes in a stratified medium. The major difficulty to overcome is related to the unboundedness of the stratified medium. Our method is an alternative to the use of artificial boundary conditions and to the use of integral representation formulae. The domain is bounded in such a way we can write the solution on its lateral boundaries in terms of Fourier series. The series is then truncated for the computations over the bounded domain. The problem is scalar and 2-dimensional.
Lorsque le domaine de calcul est un polygone non convexe, c'est-a-dire avec un ou plusieurs coins rentrantes, nous donnons une decomposition de la solution des equations de Maxwell en une partie reguliere et une partie singuliere. Nous prouvons que l'espace des parties singulieres est engendre par les solutions d'un probleme stationnaire simple. La regularite exacte de la solution est determinee en fonction de l'angle aux coins rentrants. Cette decomposition mathematique permet alors de construire un algorithme de resolution numerique des equations de Maxwell dans un polygone non convexe. Cet article est la suite de la note [3]. Cette methodologie peut egalement s'appliquer a l'equation de Helmholtz ou au systeme de Lame.
accord avec les conditions générales d'utilisation (http:/
