Three Mach-uniform algorithms

We focus on the construction of Mach-uniform algorithms. The basic idea is to remove the severe time step restrictions for low speed flows, by treating the acoustic terms implicitly. The way to solve the obtained semi-implicit system can be chosen. Three different solution techniques are presented, varying between a fully coupled algorithm and a fully segregated pressure correction algorithm. We show that the number of time steps to reach steady state is comparable for the fully coupled as well as the fully segregated method. Therefore, the more segregation is introduced, the more efficient the calculation can be done.     

Convergence to steady state of solutions of the euler equations,I

In this paper we analyze the convergence to steady state of solutions of the compressible and the incompressible isentropic Euler equations in two space dimensions. In the compressible case, the original equations do not converge. We replace the equation of continuity with an elliptic equation for the density, obtaining a new set of equations, which have the same steady solution. In the incompressible case, the equation of continuity is replaced by a Poisson equation for the pressure. In both cases, we linearize the equations around a steady solution and show that the unsteady solution of the linearized equations converges to the steady solution, if the steady solution is sufficiently smooth. In the proof we consider how the energy of the time dependent part developes with time, and find that it decrease exponentially.     

Global existence of weak solutions for a viscous two-phase model

The purpose of this paper is to explore a viscous two-phase liquid-gas model relevant for well and pipe flow. Our approach relies on applying suitable modifications of techniques previously used for studying the single-phase isothermal Navier-Stokes equations. A main issue is the introduction of a novel two-phase variant of the potential energy function needed for obtaining fundamental a priori estimates. We derive an existence result for weak solutions in a setting where transition to single-phase flow is guaranteed not to occur when the initial state is a true mixture of both phases. Some numerical examples are also included in order to demonstrate characteristic behavior of solutions. In particular, we illustrate how two-phase flow is genuinely different compared to single-phase flow concerning the behavior of an initial mass discontinuity.     

A Galerkin boundary node method and its convergence analysis

The boundary node method (BNM) exploits the dimensionality of the boundary integral equation (BIE) and the meshless attribute of the moving least-square (MLS) approximations. However, since MLS shape functions lack the property of a delta function, it is difficult to exactly satisfy boundary conditions in BNM. Besides, the system matrices of BNM are non-symmetric.     

On global and local mesh refinements by a generalized conforming bisection algorithm

We examine a generalized conforming bisection (GCB-)algorithm which allows both global and local nested refinements of the triangulations without generating hanging nodes. It is based on the notion of a mesh density function which prescribes where and how much to refine the mesh. Some regularity properties of generated sequences of refined triangulations are proved. Several numerical tests demonstrate the efficiency of the proposed bisection algorithm. It is also shown how to modify the GCB-algorithm in order to generate anisotropic meshes with high aspect ratios.     

A multiwave approximate Riemann solver for ideal MHD based on relaxation. I: theoretical framework

We present a relaxation system for ideal magnetohydrodynamics (MHD) that is an extension of the Suliciu relaxation system for the Euler equations of gas dynamics. From it one can derive approximate Riemann solvers with three or seven waves, that generalize the HLLC solver for gas dynamics. Under some subcharacteristic conditions, the solvers satisfy discrete entropy inequalities, and preserve positivity of density and internal energy. The subcharacteristic conditions are nonlinear constraints on the relaxation parameters relating them to the initial states and the intermediate states of the approximate Riemann solver itself. The 7-wave version of the solver is able to resolve exactly all material and Alfven isolated contact discontinuities. Practical considerations and numerical results will be provided in another paper.     

On the boundary element method for the Signorini problem of the Laplacian

Summary For the Laplace equation with Signorini boundary conditions two equivalent boundary variational inequality formulations are deduced. We investigate the discretization by a boundary element Galerkin method and obtain quasi-optimal asymptotic error estimates in the underlying Sobolev spaces. An algorithm based on the decomposition-coordination method is used to solve the discretized problems. Numerical examples confirm the predicted rate of convergence.     

Analysis and finite element approximation of an optimal control problem in electrochemistry with current density controls

Summary. An optimal control problem for impressed cathodic systems in electrochemistry is studied. The control in this problem is the current density on the anode. A matching objective functional is considered. We first demonstrate the existence and uniqueness of solutions for the governing partial differential equation with a nonlinear boundary condition. We then prove the existence of an optimal solution. Next, we derive a necessary condition of optimality and establish an optimality system of equations. Finally, we define a finite element algorithm and derive optimal error estimates. Received March 10, 1993 / Revised version received July 4, 1994     

Global weak solutions for a compressible gas-liquid model with well-formation interaction

The objective of this work is to explore a compressible gas-liquid model designed for modeling of well flow processes. We build into the model well-reservoir interaction by allowing flow of gas between well and formation (surrounding reservoir). Inflow of gas and subsequent expansion of gas as it ascends towards the top of the well (a so-called gas kick) represents a major concern for various well operations in the context of petroleum engineering. We obtain a global existence result under suitable assumptions on the regularity of initial data and the rate function that controls the flow of gas between well and formation. Uniqueness is also obtained by imposing more regularity on the initial data. The key estimates are to obtain appropriate lower and upper bounds on the gas and liquid masses. For that purpose we introduce a transformed version of the original model that is highly convenient for analysis of the original model. In particular, in the analysis of the transformed model additional terms, representing well-formation interaction, can be treated by natural extensions of arguments that previously have been employed for the single-phase Navier-Stokes model. The analysis ensures that transition to single-phase regions do not appear when the initial state is a true gas-liquid mixture.     

On the modeling and simulation of boundary flow through partially open pipe ends

The determination of boundary conditions for the Euler equations of gas dynamics in a pipe with partially open pipe ends is considered. The boundary problem is formulated in terms of the exact solution of the Riemann problem and of the St. Venant equation for quasi-steady flow so that a pressure-driven calculation of boundary conditions is defined. The resulting set of equations is solved by a Newton scheme. The proposed algorithm is able to solve for all inflow and outflow situations including choked and supersonic flow.Received: August 7, 2002; revised: November 11, 2002     

Infinite reload options: Pricing and analysis

Infinite reload options allow the user to exercise his reload right as often as he chooses during the lifetime of the contract. Each time a reload occurs, the owner receives new options where the strike price is set to the current stock price. We consider a modified version of the infinite reload option contract where the strike price of the new options received by the owner is increased by a certain percentage; we refer to this new contract as an increased reload option. The pricing problem for this modified contract is characterized as an impulse control problem resulting in a Hamilton–Jacobi–Bellman equation. We use fully implicit timestepping and prove that the discretized equations are monotone, stable and consistent, implying convergence to the viscosity solution. We also derive a globally convergent iterative method for solving the non-linear discrete equations. Numerical examples show that both the exercise policy and the option value are very sensitive to the percentage increase in the reload strike.     

A penalty method for a mixed nonlinear complementarity problem

We propose a power penalty method for a mixed nonlinear complementarity problem (MNCP) and show that the solution to the penalty equation converges to that of the MNCP exponentially as the penalty parameter approaches infinity, provided that the mapping involved in the MNCP is both continuous and ξ-monotone. Furthermore, a convergence theorem is established when the monotonicity assumption on the mapping is removed. To demonstrate the usefulness and the convergence rates of this method, we design a non-trivial test MNCP problem arising in shape-preserving bi-harmonic interpolation and apply our method to this test problem. The numerical results confirm our theoretical findings.     

High-order discretization and multigrid solution of elliptic nonlinear constrained optimal control problems

High(-mixed)-order finite difference discretization of optimality systems arising from elliptic nonlinear constrained optimal control problems are discussed. For the solution of these systems, an efficient and robust multigrid algorithm is presented. Theoretical and experimental results show the advantages of higher-order discretization and demonstrate that the proposed multigrid scheme is able to solve efficiently constrained optimal control problems also in the limit case of bang-bang control.     

An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method

An approximation scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using two methods. Standard mixed finite element is used for the Darcy velocity equation. A characteristics-mixed finite element method is presented for the concentration equation. Characteristic approximation is applied to handle the convection part of the concentration equation, and a lowest-order mixed finite element spatial approximation is adopted to deal with the diffusion part. Thus, the scalar unknown concentration and the diffusive flux can be approximated simultaneously. In order to derive the optimal L²$L^2$-norm error estimates, a post-processing step is included in the approximation to the scalar unknown concentration. This scheme conserves mass globally; in fact, on the discrete level, fluid is transported along the approximate characteristics. Numerical experiments are presented finally to validate the theoretical analysis.     

Convergence result for the constraint preserving mid-point scheme for micromagnetism

An important progress was recently done in numerical approximation of weak solutions to a micromagnetic model equation. The problem with the nonconvex side-constraint of preserving the length of the magnetization was tackled by using reduced integration. Several schemes were proposed and their convergence to weak solutions was proved. All schemes were derived from the Landau–Lifshitz–Gilbert form of the micromagnetic equation. However, when the precessional term in the original Landau–Lifshitz (LL) form of the micromagnetic equation tends to zero, the above schemes become unusable.     

Error estimate for the upwind finite volume method for the nonlinear scalar conservation law

In this paper we estimate the error of upwind first order finite volume schemes applied to scalar conservation laws. As a first step, we consider standard upwind and flux finite volume scheme discretization of a linear equation with space variable coefficients in conservation form. We prove that, in spite of their lack of consistency, both schemes lead to a first order error estimate. As a final step, we prove a similar estimate for the nonlinear case. Our proofs rely on the notion of geometric corrector, introduced in our previous paper by Bouche et al. (2005) [24] in the context of constant coefficient linear advection equations.     

A posteriori error estimators for the Stokes equations II non-conforming discretizations

Summary We present an a posteriori error estimator for the non-conforming Crouzeix-Raviart discretization of the Stokes equations which is based on the local evaluation of residuals with respect to the strong form of the differential equation. The error estimator yields global upper and local lower bounds for the error of the finite element solution. It can easily be generalized to the stationary, incompressible Navier-Stokes equations and to other non-conforming finite element methods. Numerical examples show the efficiency of the proposed error estimator.     

On variable stepsize Runge-Kutta approximations of a Cauchy problem for the evolution equation

We consider a Cauchy problem for the sectorial evolution equation with generally variable operator in a Banach space. Variable stepsize discretizations of this problem by means of a strongly A(φ)-stable Runge-Kutta method are studied. The stability and error estimates of the discrete solutions are derived for wider families of nonuniform grids than quasiuniform ones (in particular, if the operator in question is constant or Lipschitz-continuous, for arbitrary grids).     

Monotonicity of control volume methods

Robustness of numerical methods for multiphase flow problems in porous media is important for development of methods to be used in a wide range of applications. Here, we discuss monotonicity for a simplified problem of single-phase flow, but where the simulation grids and media are allowed to be general, posing challenges to control-volume methods. We discuss discrete formulations of the maximum principle and derive sufficient criteria for discrete monotonicity for arbitrary nine-point control-volume discretizations for conforming quadrilateral grids in 2D. These criteria are less restrictive than the M-matrix property. It is shown that it is impossible to construct nine-point methods which unconditionally satisfy the monotonicity criteria when the discretization satisfies local conservation and exact reproduction of linear potential fields. Numerical examples are presented which show the validity of the criteria for monotonicity. Further, the impact of nonmonotonicity is studied. Different behavior for different discretization methods is illuminated, and simple ideas are presented for improvement in terms of monotonicity.     

An optimal order convergence for a weak formulation of the compressible Stokes system with inflow boundary condition

Summary. We formulate the compressible Stokes system given in (1.1) into a (new) weak formulation (2.1). A finite element method for this is presented. Existence and uniqueness of the finite element method is shown. An optimal error estimate for the numerical approximation is obtained. Numerical examples are given, showing its efficiency and rates of convergence of the approximate solutions that results from the discrete problem (3.1). Received October 20, 1996 / Revised version received January 21, 1999 / Published online: April 20, 2000