In this article, we show the convergence of a class of numerical schemes for certain maximal monotone evolution systems; a by-product of this resultsis the existence of solutions in cases which had not been previouslytreated. The order of these schemes is 1/2 in general and 1 when the only non Lipschitz continuous term is the subdifferential of the indicatrix of a closed convex set. In the case of Prandtl'srheological model, our estimates in maximum norm do not dependon spatial dimension.

For a nonconforming finite element approximation of an elliptic model problem, we propose a posteriori error estimates in the energy norm which use as an additive term the “post-processing error” between the original nonconforming finite element solution and an easy computable conforming approximation of that solution.Thus, for the error analysis, the existing theory from the conformingcase can be used together with some simple additional arguments.As an essential point, the property is exploited that the nonconformingfinite element space contains as a subspace a conforming finite element space of first order. This property is fulfilled for many knownnonconforming spaces. We prove local lower and global upper a posteriori error estimates for an enhanced error measure which is the discretization error in the discrete energy norm plus the error of the best representation of the exact solution by a function in the conforming space used for thepost-processing. We demonstrate that the idea to use a computed conforming approximation ofthe nonconforming solution can be applied also to derive an a posteriorierror estimate for a linear functional of the solution which representssome quantity of physical interest.

Let Φ : H → R be a C 2 function on a real Hilbert space and ∑ ⊂ H x R the manifold defined by ∑ := Graph (Φ). We study the motion of a material point with unit mass, subjected to stay on Σ and which moves under the action of the gravity force(characterized by g>0), the reaction force and the friction force ( $\gamma>0$ is the friction parameter). For any initial conditions at time t=0, we prove the existence of a trajectory x(.) defined on R +. We are then interested in the asymptotic behaviour of the trajectories when t → +∞. More precisely, we prove the weak convergence of the trajectories when Φ is convex. When Φ admits a strong minimum, we show moreover that the mechanical energy exponentially decreases to its minimum.

In this paper we consider the Maxwell resolvent operator and its finite elementapproximation. In this framework it is natural the use of the edge elementspaces and to impose the divergence constraint in a weaksense with the introduction of a Lagrange multiplier, followingan idea by Kikuchi [14].We shall review some of the known properties for edge elementapproximations and prove some new result. In particular we shall prove auniform convergence in the L 2 norm for the sequence of discrete operators.These results, together with a general theory introduced by Brezzi, Rappaz andRaviart [8], allow an immediate proof of convergence for thefinite element approximation of the time-harmonicMaxwell system.

Different effective boundary conditions or wall laws for unsteady incompressible Navier-Stokes equations over rough domains are derived in the laminar setting. First and second order unsteady wall laws are proposed using two scale asymptotic expansion techniques. The roughness elements are supposed to be periodic and the influence of the rough boundary is incorporated through constitutive constants. These constants are obtained by solving steady Stokes problems and so they are calculated only once. Numerical tests are presented to validate and compare the proposed boundary conditions.

In this paper, we are interested in finding the optimal shapeof a magnet. The criterion to maximize is the jump of theelectromagnetic field between two different configurations.We prove existence of an optimal shape into a natural classof domains. We introduce a quasi-Newton type algorithm whichmoves the boundary. This method is very efficient to improvean initial shape. We give some numerical results.

This article is devoted to the numerical study of a flame ball model, derived by Joulin, which obeys to a singular integro-differential equation. The numerical scheme that we analyze here, is based upon a one step method, and we are interested in its long-time behaviour. We recover the same dynamics as in the continuous case: quenching, or stabilization of the flame, depending on heat losses, and an energy input parameter.

This paper is devoted to the study of a turbulentcirculation model. Equations are derived from the “Navier-Stokes turbulentkinetic energy” system. Some simplifications are performed but attentionis focused on non linearities linked to turbulent eddy viscosity  $\nu _{t}$ . The mixing length $\ell $ acts as a parameter which controls theturbulent part in $\nu _{t}$ . The main theoretical results that we haveobtained concern the uniqueness of the solution for bounded eddy viscositiesand small values of $\ell $ and its asymptotic decreasing as $\ell\rightarrow \infty $ in more general cases. Numerical experimentsillustrate but also allow to extend these theoretical results: uniqueness isproved only for $\ell $ small enough while regular solutions are numericallyobtained for any values of $\ell $ . A convergence theorem is proved forturbulent kinetic energy: $k_{\ell }\rightarrow 0$ as $\ell \rightarrow\infty ,$ but for velocity $u_{\ell }$ we obtain only weaker results.Numerical results allow to conjecture that $k_{\ell }\rightarrow 0,$ $\nu_{t}\rightarrow \infty $ and $u_{\ell }\rightarrow 0$ as $\ell \rightarrow\infty .$ So we can conjecture that this classical turbulent model obtainedwith one degree of closure regularizes the solution.

The numerical solution of the flow of a liquid crystal governedby a particular instance of the Ericksen–Leslie equations is considered.Convergence results for this system rely crucially upon energyestimates which involve H2 (Ω) norms of the director field. Weshow how a mixed method may be used to eliminate the need forHermite finite elements and establish convergence of the method.

This paper addresses the recovery of piecewise smooth functions from their discrete data.Reconstruction methods using both pseudo-spectral coefficients andphysical space interpolants have been discussed extensively in theliterature, and it is clear that an a priori knowledge of the jumpdiscontinuity location is essential for any reconstruction techniqueto yield spectrally accurate results with high resolution near thediscontinuities. Hence detection of the jump discontinuities iscritical for all methods. Here we formulate a new localized reconstruction method adapted from themethod developed in Gottlieb and Tadmor (1985) and recently revisited in Tadmor and Tanner (in press). Our procedure incorporates the detection of edges into the reconstruction technique. The methodis robust and highly accurate, yielding spectral accuracy up to a smallneighborhood of the jump discontinuities. Results are shown inone and two dimensions.

In this paper we combine the dual-mixed finite element method with a Dirichlet-to-Neumann mapping(given in terms of a boundary integral operator) to solve linear exterior transmission problems inthe plane. As a model we consider a second order elliptic equation in divergence form coupled withthe Laplace equation in the exterior unbounded region. We show that the resulting mixed variationalformulation and an associated discrete scheme using Raviart-Thomas spaces are well posed, and derivethe usual Cea error estimate and the corresponding rate of convergence. In addition, we develop twodifferent a-posteriori error analyses yielding explicit residual and implicit Bank-Weiser typereliable estimates, respectively. Several numerical results illustrate the suitability of theseestimators for the adaptive computation of the discrete solutions.

The aim of this work is to deduce the existence of solutionof a coupled problem arising in elastohydrodynamiclubrication. The lubricant pressure and concentration aremodelled by Reynolds equation, jointly with the free-boundaryElrod-Adams model in order to take into account cavitationphenomena. The bearing deformation is solution of Koitermodel for thin shells. The existence of solution to thevariational problem presents some difficulties: the coupledcharacter of the equations, the nonlinear multivaluedoperator associated to cavitation and the fact of writing theelastic and hydrodynamic equations on two different domains.In a first step, we regularize the Heaviside operator.Additional difficulty related to the differentdomains is circumvented by means of prolongation andrestriction operators, arriving to a regularized coupledproblem. This one is decoupled into elastic and hydrodynamicparts, and we prove the existence of a fixed point for theglobal operator. Estimations obtained for theregularized problem allow us to prove the existence ofsolution to the original one. Finally, a numerical method is proposed in orderto simulate a real journal-bearing device and illustrate the qualitative andquantitative properties of the solution.

We consider a new formulation for finite volume element methods, which is satisfied byknown finite volume methods and itcan be used to introduce new ones. This framework results by approximating the test function in theformulation of finite element method.We analyze piecewise linear conforming or nonconforming approximations on nonuniform triangulations andprove optimal order H 1-norm and L 2-norm errorestimates.

This paper is concerned with the coupling of two models for the propagation of particles in scattering media. The first model is a linear transport equation of Boltzmann type posed in the phase space (position and velocity). It accurately describes the physics but is very expensive to solve. The second model is a diffusion equation posed in the physical space. It is only valid in areas of high scattering, weak absorption, and smooth physical coefficients, but its numerical solution is much cheaper than that of transport. We are interested in the case when the domain is diffusive everywhere except in some small areas, for instance non-scattering or oscillatory inclusions. We present a natural coupling of the two models that accounts for both the diffusive and non-diffusive regions. The interface separating the models is chosen so that the diffusive regime holds in its vicinity to avoid the calculation of boundary or interface layers. The coupled problem is analyzed theoretically and numerically. To simplify the presentation, the transport equation is written in the even parity form. Applications include, for instance, the treatment of clear or spatially inhomogeneous regions in near-infra-red spectroscopy, which is increasingly being used in medical imaging for monitoring certain properties of human tissues.

Using systematically a tricky idea of N.V. Krylov, we obtain general results on the rate of convergence of a certain class of monotone approximation schemes for stationary Hamilton-Jacobi-Bellman equations with variable coefficients. This result applies in particular to control schemes based on the dynamic programming principle and to finite difference schemes despite, here, we are not able to treat the most general case. General results have been obtained earlier by Krylov for finite difference schemes in the stationary case with constant coefficients and in the time-dependent case with variable coefficients by using control theory and probabilistic methods. In this paper we are able to handle variable coefficients by a purely analytical method. In our opinion this way is far simpler and, for the cases we can treat, it yields a better rate of convergence than Krylov obtains in the variable coefficients case.

Estimates for the combined effect of boundaryapproximation and numerical integration on the approximation of(simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficientsin convex domains with curved boundary by an isoparametric mixed finite element method, which,in the particular case of bending problems ofaniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending momenttensor field $\Psi= (\psi_{ij})_{1 \le i,j \le 2}$ anddisplacement field `u', have been developed.

In this work, we investigate the PerfectlyMatched Layers (PML) introduced by Bérenger [3] for designing efficient numerical absorbing layers in electromagnetism.We make a mathematical analysis of this model, first via a modalanalysis with standard Fourier techniques, then via energytechniques. We obtain uniform in time stability results (that makeprecise some results known in the literature) and state some energydecay results that illustrate the absorbing properties of themodel. This last technique allows us to prove the stability of theYee's scheme for discretizing PML's.

We study the asymptotic behavior of a semi-discrete numerical approximation for a pair of heat equations ut = Δu, vt = Δv in Ω x (0,T); fully coupled by the boundary conditions $\frac{\partial u}{\partial\eta} = u^{p_{11}}v^{p_{12}}$ , $\frac{\partial v}{\partial\eta} = u^{p_{21}}v^{p_{22}}$ on ∂Ω x (0,T), where Ω is a bounded smooth domain in ${\mathbb{R}}^d$ . We focus in the existence or not of non-simultaneous blow-up for a semi-discrete approximation (U,V). We prove that if U blows up in finite time then V can fail to blow up if and only if p 11 > 1 and p 21 < 2(p 11 - 1), which is the same condition as the one for non-simultaneous blow-up in the continuous problem. Moreover, we find that if the continuous problem has non-simultaneous blow-up then the same is true for the discrete one. We also prove some results about the convergence of the scheme and the convergence of the blow-up times.

A modal synthesis method to solve the elastoacoustic vibration problemis analyzed. A two-dimensional coupled fluid-solid system is considered;the solid is described by displacement variables, whereas displacement potential is used for the fluid. A particular modal synthesis leading to a symmetric eigenvalue problem is introduced. Finite element discretizations with Lagrangian elements are considered for solving the uncoupled problems.Convergence for eigenvalues and eigenfunctions is proved, error estimates are given, and numerical experiments exhibiting the good performance of the method are reported.

This paper deals with the numerical approximation of mild solutions of elliptic-parabolic equations, relying on the existence results of Bénilan and Wittbold (1996). We introduce a new and simple algorithm based on Halpern's iteration for nonexpansive operators (Bauschke, 1996; Halpern, 1967; Lions, 1977), which is shown to be convergent in the degenerate case, and compare it with existing schemes (Jäger and Kačur, 1995; Kačur, 1999).