Let $\mathcal{Q}$ be a partition of a polygonal domain of the plan into convexe quadrilaterals. Given a regular function f , we construct a function πƒ ∈ C2(Ω), interpolating position values and derivatives of f up of order 2 at vertices of $\mathcal{Q}.$ On each quadrilateral $Q\in\mathcal{Q},$ πƒ|Q is a finite element obtained from a polynomial scheme of FVS type by adding some rational functions.

The aim of this work is to introduce and to analyze new algorithms for solving the transport neutronique equation in 2D geometry. These algorithms present the duplicate favors to be, on the one hand faster than some classic algorithms and easily to be implemented and naturally deviced for parallelisation on the other hand. They are based on a splitting of the collision operator holding amount of caracteristics of the transport operator. Some numerical results are given at the end of this work.

We consider a model problem (with constant coefficients and simplified geometry) for the boundary layer phenomena which appear in thin shell theory as the relative thickness ε of the shell tends to zero. For ε = 0 our problem is parabolic, then it is a model of developpable surfaces. Boundary layers along and across the characteristic have very different structure. It also appears internal layers associated with propagations of singularities along the characteristics. The special structure of the limit problem often implies solutions which exhibit distributional singularities along the characteristics. The corresponding layers for small ε have a very large intensity. Layers along the characteristics have a special structure involving subspaces; the corresponding Lagrange multipliers are exhibited. Numerical experiments show the advantage of adaptive meshes in these problems.

The hydrostatic approximation of the incompressible 3D stationary Navier-Stokes equations is widely used in oceanography and other applied sciences. It appears through a limit process due to the anisotropy of the domain in use, an ocean, and it is usually studied as such. We consider in this paper an equivalent formulation to this hydrostatic approximation that includes Coriolis force and an additional pressure term that comes from taking into account the pressure in the state equation for the density. It therefore models a slight dependence of the density upon compression terms. We study this model as an independent mathematical object and prove an existence theorem by means of a mixed variational formulation. The proof uses a family of finite element spaces to discretize the problem coupled with a limit process that yields the solution. We finish this paper with an existence and uniqueness result for the evolutionary linear problem associated to this model. This problem includes the same additional pressure term and Coriolis force.

We consider solutions to the time-harmonic Maxwell's Equations of a TE (transverse electric) nature. For such solutions we provide a rigorous derivation of the leading order boundary perturbations resulting from the presence of a finite number of interior inhomogeneities of small diameter. We expect that these formulas will form the basis for very effective computational identification algorithms, aimed at determining information about the inhomogeneities from electromagnetic boundary measurements.

We derive in this article some models of Cahn-Hilliard equations in nonisotropic media. These models, based on constitutive equations introduced by Gurtin in [19], take the work of internal microforces and also the deformations of the material into account. We then study the existence and uniqueness of solutions and obtain the existence of finite dimensional attractors.

Fast singular oscillating limits of the three-dimensional "primitive" equations of geophysical fluid flows are analyzed. We prove existence on infinite time intervals of regular solutions to the 3D "primitive" Navier-Stokes equations for strong stratification (large stratification parameter N). This uniform existence is proven for periodic or stress-free boundary conditions for all domain aspect ratios, including the case of three wave resonances which yield nonlinear "$2\frac{1}{2}$ dimensional" limit equations for N → +∞; smoothness assumptions are the same as for local existence theorems, that is initial data in Hα, α ≥ 3/4. The global existence is proven using techniques of the Littlewood-Paley dyadic decomposition. Infinite time regularity for solutions of the 3D "primitive" Navier-Stokes equations is obtained by bootstrapping from global regularity of the limit resonant equations and convergence theorems.

We study a finite volume method, used to approximate the solution of the linear two dimensional convection diffusion equation, with mixed Dirichlet and Neumann boundary conditions, on Cartesian meshes refined by an automatic technique (which leads to meshes with hanging nodes). We propose an analysis through a discrete variational approach, in a discrete H1 finite volume space. We actually prove the convergence of the scheme in a discrete H1 norm, with an error estimate of order O(h) (on meshes of size h).

The present work is a mathematical analysis of two algorithms, namely the Roothaan and the level-shifting algorithms, commonly used in practice to solve the Hartree-Fock equations. The level-shifting algorithm is proved to be well-posed and to converge provided the shift parameter is large enough. On the contrary, cases when the Roothaan algorithm is not well defined or fails in converging are exhibited. These mathematical results are confronted to numerical experiments performed by chemists.

We are interested in a barotropic motion of the non-Newtonian bipolar fluids . We consider a special case where the stress tensor is expressed in the form of potentials depending on eii and $(\frac{\partial e_{ij}}{\partial x_{k}})$. We prove the asymptotic stability of the rest state under the assumption of the regularity of the potential forces.

A new numerical method based on fictitious domain methods for shape optimization problems governed by the Poisson equation is proposed. The basic idea is to combine the boundary variation technique, in which the mesh is moving during the optimization, and efficient fictitious domain preconditioning in the solution of the (adjoint) state equations. Neumann boundary value problems are solved using an algebraic fictitious domain method. A mixed formulation based on boundary Lagrange multipliers is used for Dirichlet boundary problems and the resulting saddle-point problems are preconditioned with block diagonal fictitious domain preconditioners. Under given assumptions on the meshes, these preconditioners are shown to be optimal with respect to the condition number. The numerical experiments demonstrate the efficiency of the proposed approaches.

In this article we consider local solutions for stochastic Navier Stokes equations, based on the approach of Von Wahl, for the deterministic case. We present several approaches of the concept, depending on the smoothness available. When smoothness is available, we can in someway reduce the stochastic equation to a deterministic one with a random parameter. In the general case, we mimic the concept of local solution for stochastic differential equations.

The paper is devoted to analysis of an elliptic-algebraic system of equations describing heat explosion in a two phase medium filling a star-shaped domain. Three types of solutions are found: classical, critical and multivalued. Regularity of solutions is studied as well as their behavior depending on the size of the domain and on the coefficient of heat exchange between the two phases. Critical conditions of existence of solutions are found for arbitrary positive source function.

Consider the domain $Z_\epsilon=\{x\in\mathbb{R}^n ; {dist}(x,\epsilon\mathbb{Z}^n)> \epsilon^\gamma\}$ and let the free path length be defined as $\tau_\epsilon(x,v)=\inf\{t> 0 ; x-tv\in Z_\epsilon\}.$ In the Boltzmann-Grad scaling corresponding to $\gamma=\frac{n}{n-1}$, it is shown that the limiting distribution $\phi_\epsilon$ of $\tau_\epsilon$ is bounded from below by an expression of the form C/t, for some C> 0. A numerical study seems to indicate that asymptotically for large t, $\phi_\epsilon\sim C/t$. This is an extension of a previous work [J. Bourgain et al., Comm. Math. Phys.190 (1998) 491-508]. As a consequence, it is proved that the linear Boltzmann type transport equation is inappropriate to describe the Boltzmann-Grad limit of the periodic Lorentz gas, at variance with the usual case of a Poisson distribution of scatterers treated in [G. Gallavotti (1972)].

We consider a domain decomposition method for some unsteady heat conduction problem in composite structures. This linear model problem is obtained by homogenization of thin layers of fibres embedded into some standard material. For ease of presentation we consider the case of two space dimensions only. The set of finite element equations obtained by the backward Euler scheme is parallelized in a problem-oriented fashion by some noniterative overlapping domain splitting method, eventually enhanced by inexpensive local iterations to reduce the overlap. We present a detailed convergence analysis of this algorithm which is particularly well appropriate to handle fibre layers of nonlinear material. Special emphasis is to take into account the specific regularity properties of the present mathematical model. Numerical experiments show the reliability of the theoretical predictions.

The LBB condition is well-known to guarantee the stability of a finite element (FE) velocity - pressure pair in incompressible flow calculations. To ensure the condition to be satisfied a certain constant should be positive and mesh-independent. The paper studies the dependence of the LBB condition on the domain geometry. For model domains such as strips and rings the substantial dependence of this constant on geometry aspect ratios is observed. In domains with highly anisotropic substructures this may require special care with numerics to avoid failures similar to those when the LBB condition is violated. In the core of the paper we prove that for any FE velocity-pressure pair satisfying usual approximation hypotheses the mesh-independent limit in the LBB condition is not greater than its continuous counterpart, the constant from the Nečas inequality. For the latter the explicit and asymptotically accurate estimates are proved. The analytic results are illustrated by several numerical experiments.

The phase relaxation model is a diffuse interface model with small parameter ε which consists of a parabolic PDE for temperature θ and an ODE with double obstacles for phase variable χ. To decouple the system a semi-explicit Euler method with variable step-size τ is used for time discretization, which requires the stability constraint τ ≤ ε. Conforming piecewise linear finite elements over highly graded simplicial meshes with parameter h are further employed for space discretization. A posteriori error estimates are derived for both unknowns θ and χ, which exhibit the correct asymptotic order in terms of ε, h and τ. This result circumvents the use of duality, which does not even apply in this context. Several numerical experiments illustrate the reliability of the estimators and document the excellent performance of the ensuing adaptive method.

In dimension one it is proved that the solution to a total variation-regularized least-squares problem is always a function which is "constant almost everywhere" , provided that the data are in a certain sense outside the range of the operator to be inverted. A similar, but weaker result is derived in dimension two.

We propose a new formulation of the 3D Boltzmann non linear operator, without assuming Grad's angular cutoff hypothesis, and for intermolecular laws behaving as 1/rs, with s> 2. It involves natural pseudo differential operators, under a form which is analogous to the Landau operator. It may be used in the study of the associated equations, and more precisely in the non homogeneous framework.