Degenerate two-phase incompressible flow II: regularity, stability and stabilization

In this paper, we analyze a coupled system of highly degenerate elliptic-parabolic partial differential equations for two-phase incompressible flow in porous media. This system involves a saturation and a global pressure (or a total flow velocity). First, we show that the saturation is Hölder continuous both in space and time and the total velocity is Hölder continuous in space (uniformly in time). Applying this regularity result, we then establish the stability of the saturation and pressure with respect to initial and boundary data, from which uniqueness of the solution to the system follows. Finally, we establish a stabilization result on the asymptotic behavior of the saturation and pressure; we prove that the solution to the present system converges (in appropriate norms) to the solution of a stationary system as time goes to infinity. An example is given to show typical regularity of the saturation.     

Spaces of distributions, interpolation by translates of a basis function and error estimates

Summary. Interpolation with translates of a basis function is a common process in approximation theory. The most elementary form of the interpolant consists of a linear combination of all translates by interpolation points of a single basis function. Frequently, low degree polynomials are added to the interpolant. One of the significant features of this type of interpolant is that it is often the solution of a variational problem. In this paper we concentrate on developing a wide variety of spaces for which a variational theory is available. For each of these spaces, we show that there is a natural choice of basis function. We also show how the theory leads to efficient ways of calculating the interpolant and to new error estimates. Received December 10, 1996 / Revised version received August 29, 1997     

A stabilized finite element method for the incompressible magnetohydrodynamic equations

Summary. We propose and analyze a stabilized finite element method for the incompressible magnetohydrodynamic equations. The numerical results that we present show a good behavior of our approximation in experiments which are relevant from an industrial viewpoint. We explain in particular in the proof of our convergence theorem why it may be interesting to stabilize the magnetic equation as soon as the hydrodynamic diffusion is small and even if the magnetic diffusion is large. This observation is confirmed by our numerical tests. Received August 31, 1998 / Revised version received June 16, 1999 / Published online June 21, 2000     

A priori error estimates for approximate solutions to convex conservation laws

Summary. We introduce a new technique for proving a priori error estimates between the entropy weak solution of a scalar conservation law and a finite–difference approximation calculated with the scheme of Engquist-Osher, Lax-Friedrichs, or Godunov. This technique is a discrete counterpart of the duality technique introduced by Tadmor [SIAM J. Numer. Anal. 1991]. The error is related to the consistency error of cell averages of the entropy weak solution. This consistency error can be estimated by exploiting a regularity structure of the entropy weak solution. One ends up with optimal error estimates. Received December 21, 2001 / Revised version received February 18, 2002 / Published online June 17, 2002     

Some remarks on the level set flow by anisotropic curvature

We show that almost all level sets of the unique viscosity solution for general anisotropic mean curvature flow satisfy a weak form of the flow equation. This generalizes the case of isotropic mean curvature flow studied by Evans and Spruck, in which a relation between the viscosity solution and Brakke's varifold mean curvature flow is established. Received June 4, 1998 / Accepted February 26, 1999     

A penalization method to take into account obstacles in incompressible viscous flows

From the Navier-Stokes/Brinkman model, a penalization method has been derived by several authors to compute incompressible Navier-Stokes equations around obstacles. In this paper, convergence theorems and error estimates are derived for two kinds of penalization. The first one corresponds to penalization inducing a Darcy equation in the solid body, the second one corresponds to a penalization and induces a Brinkman equation in the body. Numerical tests are performed to confirm the efficiency and accuracy of the method. Received August 10, 1997     

New anisotropic a priori error estimates

Summary. We prove a priori anisotropic estimates for the and interpolation error on linear finite elements. The full information about the mapping from a reference element is employed to separate the contribution to the elemental error coming from different directions. This new error estimate does not require the “maximal angle condition”. The analysis has been carried out for the 2D case, but may be extended to three dimensions. Numerical experiments have been carried out to test our theoretical results. Received March 3, 2000 / Revised version received June 27, 2000 / Published online April 5, 2001     

Convection–diffusion equation with space–time ergodic random flow

We prove the homogenization of convection-diffusion in a time-dependent, ergodic, incompressible random flow which has a bounded stream matrix and a constant mean drift. We also prove two variational formulas for the effective diffusivity. As a consequence, we obtain both upper and lower bounds on the effective diffusivity. Received: 17 December 1996/Revised revision: 9 February 1998     

Partial regularity of weak solutions for the curve shortening flow

We study the regularity of certain weak solutions for the curve shortening flow in arbitrary codimension. These solutions arise as limits of a regularization process which is related to an approach suggested by Calabi. We prove that the set of times for which such a weak solution is not smooth has Hausdorff dimension at most ?. Received: 23 May 1998 / Revised version: 7 September 1998     

A posteriori error estimate for finite volume approximations to singularly perturbed nonlinear convection-diffusion equations

Summary. This paper is devoted to the study of a posteriori and a priori error estimates for the scalar nonlinear convection diffusion equation . The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the -norm in the situation, where the diffusion parameter is smaller or comparable to the mesh size. Numerical experiments underline the theoretical results. Received February 25, 1999 / Revised version received July 6, 1999 / Published online August 2, 2000     

On positivity preserving finite volume schemes for Euler equations

Summary. We consider the positivity preserving property of first and higher order finite volume schemes for one and two dimensional Euler equations of gas dynamics. A general framework is established which shows the positivity of density and pressure whenever the underlying one dimensional first order building block based on an exact or approximate Riemann solver and the reconstruction are both positivity preserving. Appropriate limitation to achieve a high order positivity preserving reconstruction is described. Received May 20, 1994     

Stability and error estimates for the $\theta$-method for strongly monotone and infinitely stiff evolution equations

Summary.   For evolution equations with a strongly monotone operator we derive unconditional stability and discretization error estimates valid for all . For the -method, with , we prove an error estimate , if , where is the maximal integration step for an arbitrary choice of sequence of steps and with no assumptions about the existence of the Jacobian as well as other derivatives of the operator , and an optimal estimate under some additional relation between neighboring steps. The first result is an improvement over the implicit midpoint method , for which an order reduction to sometimes may occur for infinitely stiff problems. Numerical tests illustrate the results. Received March 10, 1999 / Revised version received April 3, 2000 / Published online February 5, 2001     

Asymptotic formula for the error in cardinal interpolation

Summary. We give the asymptotic formula for the error in cardinal interpolation. We generalize the Mazur Orlicz Theorem for periodic function. Received February 22, 1999 / Revised version received October 15, 1999 / Published online March 20, 2001     

$L^1$-Stability and error estimates for approximate Hamilton-Jacobi solutions

Summary.   We study the -stability and error estimates of general approximate solutions for the Cauchy problem associated with multidimensional Hamilton-Jacobi (H-J) equations. For strictly convex Hamiltonians, we obtain a priori error estimates in terms of the truncation errors and the initial perturbation errors. We then demonstrate this general theory for two types of approximations: approximate solutions constructed by the vanishing viscosity method, and by Godunov-type finite difference methods. If we let denote the `small scale' of such approximations (– the viscosity amplitude , the spatial grad-size , etc.), then our -error estimates are of , and are sharper than the classical -results of order one half, . The main building blocks of our theory are the notions of the semi-concave stability condition and -measure of the truncation error. The whole theory could be viewed as a multidimensional extension of the -stability theory for one-dimensional nonlinear conservation laws developed by Tadmor et. al. [34,24,25]. In addition, we construct new Godunov-type schemes for H-J equations which consist of an exact evolution operator and a global projection operator. Here, we restrict our attention to linear projection operators (first-order schemes). We note, however, that our convergence theory applies equally well to nonlinear projections used in the context of modern high-resolution conservation laws. We prove semi-concave stability and obtain -bounds on their associated truncation errors; -convergence of order one then follows. Second-order (central) Godunov-type schemes are also constructed. Numerical experiments are performed; errors and orders are calculated to confirm our -theory. Received April 20, 1998 / Revised version received November 8, 1999 / Published online August 24, 2000     

A fully discrete numerical scheme for weighted mean curvature flow

Summary. We analyze a fully discrete numerical scheme approximating the evolution of n–dimensional graphs under anisotropic mean curvature. The highly nonlinear problem is discretized by piecewise linear finite elements in space and semi–implicitly in time. The scheme is unconditionally stable und we obtain optimal error estimates in natural norms. We also present numerical examples which confirm our theoretical results. Received October 2, 2000 / Published online July 25, 2001     

On estimation of the logarithmic Sobolev constant and gradient estimates of heat semigroups

Summary. This paper presents some explicit lower bound estimates of logarithmic Sobolev constant for diffusion processes on a compact Riemannian manifold with negative Ricci curvature. Let Ric≧−K for some K>0 and d, D be respectively the dimension and the diameter of the manifold. If the boundary of the manifold is either empty or convex, then the logarithmic Sobolev constant for Brownian motion is not less than max {(d d+2) ^d 1 2(d+1)D ² exp [−1−(3d+2)D ² K],     (d−1 d+1) ^d K exp [−4D√d K]} . Next, the gradient estimates of heat semigroups (including the Neumann heat semigroup and the Dirichlet one) are studied by using coupling method together with a derivative formula modified from [11]. The resulting estimates recover or improve those given in [7, 21] for harmonic functions. Received: 19 September 1995 / In revised form 11 April 1996     

Finite volume schemes for nonhomogeneous scalar conservation laws: error estimate

Summary. In this paper, we study finite volume schemes for the nonhomogeneous scalar conservation law with initial condition . The source term may be either stiff or nonstiff. In both cases, we prove error estimates between the approximate solution given by a finite volume scheme (the scheme is totally explicit in the nonstiff case, semi-implicit in the stiff case) and the entropy solution. The order of these estimates is in space-time -norm (h denotes the size of the mesh). Furthermore, the error estimate does not depend on the stiffness of the source term in the stiff case. Received October 21, 1999 / Published online February 5, 2001     

Improved smoothing-type methods for the solution of linear programs

Summary. We consider a smoothing-type method for the solution of linear programs. Its main idea is to reformulate the primal-dual optimality conditions as a nonlinear and nonsmooth system of equations, and to apply a Newton-type method to a smooth approximation of this nonsmooth system. The method presented here is a predictor-corrector method, and is closely related to some methods recently proposed by Burke and Xu on the one hand, and by the authors on the other hand. However, here we state stronger global and/or local convergence properties. Moreover, we present quite promising numerical results for the whole netlib test problem collection. Received August 9, 2000 / Revised version received September 28, 2000 / Published online June 7, 2001