A numerical method for solving the problem
A numerical scheme for the two phase Mullins-Sekerka problem.
A numerical study of effects of the axial stress on unsteady liquid pipeline flows.
A numerical study of non-cavitating and cavitating liquid flow around a hydrofoil
The results of a workshop concerning the numerical simulation of the liquid flow around a hydrofoil in non-cavitating and cavitating conditions are presented. This workshop was part of the conference “Mathematical and Numerical aspects of Low Mach Number Flows” (2004) and was aimed to investigate the capabilities of different compressible flow solvers for the low Mach number regime and for flows in which incompressible and supersonic regions are simultaneously present. Different physical models...
A numerical study of non-cavitating and cavitating liquid flow around a hydrofoil
The results of a workshop concerning the numerical simulation of the liquid flow around a hydrofoil in non-cavitating and cavitating conditions are presented. This workshop was part of the conference “Mathematical and Numerical aspects of Low Mach Number Flows” (2004) and was aimed to investigate the capabilities of different compressible flow solvers for the low Mach number regime and for flows in which incompressible and supersonic regions are simultaneously present. Different physical models...
A numerical study of some questions in vortex rings theory
A parallel algorithm for two phase multicomponent contaminant transport
We discuss the formulation of a simulator in three spatial dimensions for a multicomponent, two phase (air, water) system of groundwater flow and transport with biodegradation kinetics and wells with multiple screens. The simulator has been developed for parallel, distributed memory, message passing machines. The numerical procedures employed are a fully implicit expanded mixed finite element method for flow and either a characteristics-mixed method or a Godunov method for transport and reactions...
A parameter choice for Tikhonov regularization for solving nonlinear inverse problems leading to optimal convergence rates
We give a derivation of an a-posteriori strategy for choosing the regularization parameter in Tikhonov regularization for solving nonlinear ill-posed problems, which leads to optimal convergence rates. This strategy requires a special stability estimate for the regularized solutions. A new proof fot this stability estimate is given.
A parameter-free smoothness indicator for high-resolution finite element schemes
This paper presents a postprocessing technique for estimating the local regularity of numerical solutions in high-resolution finite element schemes. A derivative of degree p ≥ 0 is considered to be smooth if a discontinuous linear reconstruction does not create new maxima or minima. The intended use of this criterion is the identification of smooth cells in the context of p-adaptation or selective flux limiting. As a model problem, we consider a 2D convection equation discretized with bilinear finite...
A Particle Method for First-order Symmetric Systems.
A posteriori error analysis for the Crank-Nicolson method for linear Schrödinger equations
We prove a posteriori error estimates of optimal order for linear Schrödinger-type equations in the L∞(L2)- and the L∞(H1)-norm. We discretize only in time by the Crank-Nicolson method. The direct use of the reconstruction technique, as it has been proposed by Akrivis et al. in [Math. Comput. 75 (2006) 511–531], leads to a posteriori upper bounds that are of optimal order in the L∞(L2)-norm, but of suboptimal order in the L∞(H1)-norm. The optimality in the case of L∞(H1)-norm is recovered by using...
A posteriori error analysis for the Crank-Nicolson method for linear Schrödinger equations*
We prove a posteriori error estimates of optimal order for linear Schrödinger-type equations in the L∞(L2)- and the L∞(H1)-norm. We discretize only in time by the Crank-Nicolson method. The direct use of the reconstruction technique, as it has been proposed by Akrivis et al. in [Math. Comput.75 (2006) 511–531], leads to a posteriori upper bounds that are of optimal order in the L∞(L2)-norm, but of suboptimal order in the L∞(H1)-norm. The optimality in the case of L∞(H1)-norm is recovered by using...
A posteriori error analysis of Euler-Galerkin approximations to coupled elliptic-parabolic problems
We analyze Euler-Galerkin approximations (conforming finite elements in space and implicit Euler in time) to coupled PDE systems in which one dependent variable, say , is governed by an elliptic equation and the other, say , by a parabolic-like equation. The underlying application is the poroelasticity system within the quasi-static assumption. Different polynomial orders are used for the - and -components to obtain optimally convergent a priori bounds for all the terms in the error energy norm....
A posteriori error analysis of Euler-Galerkin approximations to coupled elliptic-parabolic problems
We analyze Euler-Galerkin approximations (conforming finite elements in space and implicit Euler in time) to coupled PDE systems in which one dependent variable, say u, is governed by an elliptic equation and the other, say p, by a parabolic-like equation. The underlying application is the poroelasticity system within the quasi-static assumption. Different polynomial orders are used for the u- and p-components to obtain optimally convergent a priori bounds for all the terms in the error energy...
A posteriori error analysis of the fully discretized time-dependent Stokes equations
The time-dependent Stokes equations in two- or three-dimensional bounded domains are discretized by the backward Euler scheme in time and finite elements in space. The error of this discretization is bounded globally from above and locally from below by the sum of two types of computable error indicators, the first one being linked to the time discretization and the second one to the space discretization.
A posteriori error analysis of the fully discretized time-dependent Stokes equations
The time-dependent Stokes equations in two- or three-dimensional bounded domains are discretized by the backward Euler scheme in time and finite elements in space. The error of this discretization is bounded globally from above and locally from below by the sum of two types of computable error indicators, the first one being linked to the time discretization and the second one to the space discretization.
A posteriori error bounds for reduced-basis approximations of parametrized parabolic partial differential equations
In this paper, we extend the reduced-basis methods and associated a posteriori error estimators developed earlier for elliptic partial differential equations to parabolic problems with affine parameter dependence. The essential new ingredient is the presence of time in the formulation and solution of the problem – we shall “simply” treat time as an additional, albeit special, parameter. First, we introduce the reduced-basis recipe – Galerkin projection onto a space spanned by solutions of the...
A posteriori error bounds for reduced-basis approximations of parametrized parabolic partial differential equations
In this paper, we extend the reduced-basis methods and associated a posteriori error estimators developed earlier for elliptic partial differential equations to parabolic problems with affine parameter dependence. The essential new ingredient is the presence of time in the formulation and solution of the problem – we shall “simply” treat time as an additional, albeit special, parameter. First, we introduce the reduced-basis recipe – Galerkin projection onto a space WN spanned by solutions...
A posteriori error control for the Allen–Cahn problem : circumventing Gronwall’s inequality
Phase-field models, the simplest of which is Allen–Cahn’s problem, are characterized by a small parameter that dictates the interface thickness. These models naturally call for mesh adaptation techniques, which rely on a posteriori error control. However, their error analysis usually deals with the underlying non-monotone nonlinearity via a Gronwall argument which leads to an exponential dependence on . Using an energy argument combined with a topological continuation argument and a spectral...
A posteriori error control for the Allen–Cahn problem: circumventing Gronwall's inequality
Phase-field models, the simplest of which is Allen–Cahn's problem, are characterized by a small parameter ε that dictates the interface thickness. These models naturally call for mesh adaptation techniques, which rely on a posteriori error control. However, their error analysis usually deals with the underlying non-monotone nonlinearity via a Gronwall argument which leads to an exponential dependence on ε-2. Using an energy argument combined with a topological continuation argument and...