A Combined Finite Element and Marker and Cell Method for Solving Navier-Stokes Equations.
The aim of this paper is to demonstrate how the variational equations from can be formulated and solved in some abstract Banach spaces without any a priori construction of special linearization schemes. This should be useful e.g. in the analysis of heat conduction problems and modelling of flow in porous media.
We analyze a nonlinear discrete scheme depending on second-order finite differences. This is the two-dimensional analog of a scheme which in one dimension approximates a free-discontinuity energy proposed by Blake and Zisserman as a higher-order correction of the Mumford and Shah functional. In two dimension we give a compactness result showing that the continuous problem approximating this difference scheme is still defined on special functions with bounded hessian, and we give an upper and a lower...
We analyze a nonlinear discrete scheme depending on second-order finite differences. This is the two-dimensional analog of a scheme which in one dimension approximates a free-discontinuity energy proposed by Blake and Zisserman as a higher-order correction of the Mumford and Shah functional. In two dimension we give a compactness result showing that the continuous problem approximating this difference scheme is still defined on special functions...
In this paper, we study a model for the magnetization in thin ferromagnetic films. It comes as a variational problem for -valued maps (the magnetization) of two variables : . We are interested in the behavior of minimizers as . They are expected to be -valued maps of vanishing distributional divergence , so that appropriate boundary conditions enforce line discontinuities. For finite , these line discontinuities are approximated by smooth transition layers, the so-called Néel walls. Néel...
We examine different approaches to an efficient solution of the stochastic Galerkin (SG) matrix equations coming from the Darcy flow problem with different, uncertain coefficients in apriori known subdomains. The solution of the SG system of equations is usually a very challenging task. A relatively new approach to the solution of the SG matrix equations is the reduced basis (RB) solver, which looks for a low-rank representation of the solution. The construction of the RB is usually done iteratively...
Even for a well-trained statistician the construction of a histogram for a given real-valued data set is a difficult problem. It is even more difficult to construct a fully automatic procedure which specifies the number and widths of the bins in a satisfactory manner for a wide range of data sets. In this paper we compare several histogram construction procedures by means of a simulation study. The study includes plug-in methods, cross-validation, penalized maximum likelihood and the taut string...
In this paper some of the cointegration tests applied to a single equation are compared. Many of the existent cointegration tests are simply extensions of the unit root tests applied to the residuals of the cointegrating regression and the habitual is no cointegration. However, some non residual-based tests and some tests of the opposite null hypothesis have recently appeared in literature. Monte Carlo simulations have been used for the power comparison of the nine selected tests (, , , ,...
The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is...
The paper deals with formulation and numerical solution of problems of identification of material parameters for continuum mechanics problems in domains with heterogeneous microstructure. Due to a restricted number of measurements of quantities related to physical processes, we assume additional information about the microstructure geometry provided by CT scan or similar analysis. The inverse problems use output least squares cost functionals with values obtained from averages of state problem quantities...
Domain decomposition techniques provide a powerful tool for the numerical approximation of partial differential equations. We focus on mortar finite element methods on non-matching triangulations. In particular, we discuss and analyze dual Lagrange multiplier spaces for lowest order finite elements. These non standard Lagrange multiplier spaces yield optimal discretization schemes and a locally supported basis for the associated constrained mortar spaces. As a consequence, standard efficient iterative...
Domain decomposition techniques provide a powerful tool for the numerical approximation of partial differential equations. We focus on mortar finite element methods on non-matching triangulations. In particular, we discuss and analyze dual Lagrange multiplier spaces for lowest order finite elements. These non standard Lagrange multiplier spaces yield optimal discretization schemes and a locally supported basis for the associated constrained mortar spaces. As a consequence, standard efficient iterative...
A lot of papers and books analyze analytical a posteriori error estimates from the point of view of robustness, guaranteed upper bounds, global efficiency, etc. At the same time, adaptive finite element methods have acquired the principal position among algorithms for solving differential problems in many physical and technical applications. In this survey contribution, we present and compare, from the viewpoint of adaptive computation, several recently published error estimation procedures for...