Robot motion using Delaunay triangulation.
Robust a posteriori error estimates for finite element discretizations of the heat equation with discontinuous coefficients
In this work we derive a posteriori error estimates based on equations residuals for the heat equation with discontinuous diffusivity coefficients. The estimates are based on a fully discrete scheme based on conforming finite elements in each time slab and on the A-stable -scheme with . Following remarks of [Picasso, Comput. Methods Appl. Mech. Engrg. 167 (1998) 223–237; Verfürth, Calcolo 40 (2003) 195–212] it is easy to identify a time-discretization error-estimator and a space-discretization...
Robust a posteriori error estimates for finite element discretizations of the heat equation with discontinuous coefficients
In this work we derive a posteriori error estimates based on equations residuals for the heat equation with discontinuous diffusivity coefficients. The estimates are based on a fully discrete scheme based on conforming finite elements in each time slab and on the A-stable θ-scheme with 1/2 ≤ θ ≤ 1. Following remarks of [Picasso, Comput. Methods Appl. Mech. Engrg. 167 (1998) 223–237; Verfürth, Calcolo40 (2003) 195–212] it is easy to identify a time-discretization error-estimator and a space-discretization...
Robust a priori error analysis for the approximation of degree-one Ginzburg-Landau vortices
This article discusses the numerical approximation of time dependent Ginzburg-Landau equations. Optimal error estimates which are robust with respect to a large Ginzburg-Landau parameter are established for a semi-discrete in time and a fully discrete approximation scheme. The proofs rely on an asymptotic expansion of the exact solution and a stability result for degree-one Ginzburg-Landau vortices. The error bounds prove that degree-one vortices can be approximated robustly while unstable higher...
Robust a priori error analysis for the approximation of degree-one Ginzburg-Landau vortices
This article discusses the numerical approximation of time dependent Ginzburg-Landau equations. Optimal error estimates which are robust with respect to a large Ginzburg-Landau parameter are established for a semi-discrete in time and a fully discrete approximation scheme. The proofs rely on an asymptotic expansion of the exact solution and a stability result for degree-one Ginzburg-Landau vortices. The error bounds prove that degree-one vortices can be approximated robustly while unstable higher...
Robust domain decomposition preconditioners for abstract symmetric positive definite bilinear forms
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations. The constant in the corresponding abstract...
Robust domain decomposition preconditioners for abstract symmetric positive definite bilinear forms
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations....
Robust local problem error estimation for a singularly perturbed problem on anisotropic finite element meshes
Singularly perturbed problems often yield solutions with strong directional features, e.g. with boundary layers. Such anisotropic solutions lend themselves to adapted, anisotropic discretizations. The quality of the corresponding numerical solution is a key issue in any computational simulation. To this end we present a new robust error estimator for a singularly perturbed reaction–diffusion problem. In contrast to conventional estimators, our proposal is suitable for anisotropic finite element...
Robust local problem error estimation for a singularly perturbed problem on anisotropic finite element meshes
Singularly perturbed problems often yield solutions with strong directional features, e.g. with boundary layers. Such anisotropic solutions lend themselves to adapted, anisotropic discretizations. The quality of the corresponding numerical solution is a key issue in any computational simulation. To this end we present a new robust error estimator for a singularly perturbed reaction-diffusion problem. In contrast to conventional estimators, our proposal is suitable for anisotropic finite element...
Robust monotone iterates for nonlinear singularly perturbed boundary value problems.
Robust numerical approximation of coupled Stokes' and Darcy's flows applied to vascular hemodynamics and biochemical transport*
The fully coupled description of blood flow and mass transport in blood vessels requires extremely robust numerical methods. In order to handle the heterogeneous coupling between blood flow and plasma filtration, addressed by means of Navier-Stokes and Darcy's equations, we need to develop a numerical scheme capable to deal with extremely variable parameters, such as the blood viscosity and Darcy's permeability of the arterial walls. In this paper, we describe a finite element method for...
Robust numerical approximation of coupled Stokes' and Darcy's flows applied to vascular hemodynamics and biochemical transport*
The fully coupled description of blood flow and mass transport in blood vessels requires extremely robust numerical methods. In order to handle the heterogeneous coupling between blood flow and plasma filtration, addressed by means of Navier-Stokes and Darcy's equations, we need to develop a numerical scheme capable to deal with extremely variable parameters, such as the blood viscosity and Darcy's permeability of the arterial walls. In this paper, we describe a finite element method for...
Robust operator estimates and the application to substructuring methods for first-order systems
We discuss a family of discontinuous Petrov–Galerkin (DPG) schemes for quite general partial differential operators. The starting point of our analysis is the DPG method introduced by [Demkowicz et al., SIAM J. Numer. Anal. 49 (2011) 1788–1809; Zitelli et al., J. Comput. Phys. 230 (2011) 2406–2432]. This discretization results in a sparse positive definite linear algebraic system which can be obtained from a saddle point problem by an element-wise Schur complement reduction applied to the test space....
Robust preconditioners for the matrix free truncated Newton method
New positive definite preconditioners for the matrix free truncated Newton method are given. Corresponding algorithms are described in detail. Results of numerical experiments that confirm the efficiency and robustness of the preconditioned truncated Newton method are reported.
Robust quasi-linear system identification
Robust semi-coarsening multilevel preconditioning of biquadratic FEM systems
While a large amount of papers are dealing with robust multilevel methods and algorithms for linear FEM elliptic systems, the related higher order FEM problems are much less studied. Moreover, we know that the standard hierarchical basis two-level splittings deteriorate for strongly anisotropic problems. A first robust multilevel preconditioner for higher order FEM systems obtained after discretizations of elliptic problems with an anisotropic diffusion tensor is presented in this paper. We study...
Robustness of multi-grid applied to anisotropic equations on convex domains with re-entrant corners.
Roe Scheme for Two-layer Shallow Water Equations: Application to the Strait of Gibraltar
The flow trough the Strait of Gibraltar could be analyzed as a problem of two-layer hydraulic exchange between the Atlantic ocean and the Mediterranean sea. The shallow water equations in both layers coupled together are an important tool to simulate this phenomenon. In this paper we perform an upwind schemes for hyperbolic equations based on the Roe approximate Riemann solver, to study the resulting model. The main goal assigned was to predict the location of the interface between the two layers....
Rolling-ball method for estimating the boundary of the support of a point-process intensity
Root Finding by Divided Differences.