2016, Reichelt, Sina

This note is devoted to the derivation of quantitative estimates for linear elliptic equations with coefficients that are not exactly ε-periodic and the ellipticity constant may degenerate for vanishing ε. Here ε>0 denotes the ratio between the microscopic and the macroscopic length scale. It is shown that for degenerating and non-degenerating coefficients the error between the original solution and the effective solution is of order √ε. Therefore suitable test functions are constructed via the periodic unfolding method and a gradient folding operator making only minimal additional assumptions on the given data and the effective solution with respect to the macroscopic scale.

2016, Heltai, Luca, Rotundo, Nella

When solving elliptic partial differential equations in a region containing immersed interfaces (possibly evolving in time), it is often desirable to approximate the problem using a uniform background discretisation, not aligned with the interface itself. Optimal convergence rates are possible if the discretisation scheme is enriched by allowing the discrete solution to have jumps aligned with the surface, at the cost of a higher complexity in the implementation. A much simpler way to reformulate immersed interface problems consists in replacing the interface by a singular force field that produces the desired interface conditions, as done in immersed boundary methods. These methods are known to have inferior convergence properties, depending on the global regularity of the solution across the interface, when compared to enriched methods. In this work we prove that this detrimental effect on the convergence properties of the approximate solution is only a local phenomenon, restricted to a small neighbourhood of the interface. In particular we show that optimal approximations can be constructed in a natural and inexpensive way, simply by reformulating the problem in a distributionally consistent way, and by resorting to weighted norms when computing the global error of the approximation.

2015, Colli, Pierluigi, Gilardi, Gianni, Rocca, Elisabetta, Sprekels, Jürgen

This paper is concerned with a phase field system of Cahn-Hilliard type that is related to a tumor growth model and consists of three equations in gianni terms of the variables order parameter, chemical potential and nutrient concentration. This system has been investigated in the recent papers citeCGH and citeCGRS gianni from the viewpoint of well-posedness, long time bhv and asymptotic convergence as two positive viscosity coefficients tend to zero at the same time. Here, we continue the analysis performed in citeCGRS by showing two independent sets of results as just one of the coefficents tends to zero, the other remaining fixed. We prove convergence results, uniqueness of solutions to the two resulting limit problems, and suitable error estimates

2020, Koprucki, Thomas, Maltsi, Anieza, Mielke, Alexander

The Darwin-Howie-Whelan equations are commonly used to describe and simulate the scattering of fast electrons in transmission electron microscopy. They are a system of infinitely many envelope functions, derived from the Schrödinger equation. However, for the simulation of images only a finite set of envelope functions is used, leading to a system of ordinary differential equations in thickness direction of the specimen. We study the mathematical structure of this system and provide error estimates to evaluate the accuracy of special approximations, like the two-beam and the systematic-row approximation.

2016, Reichelt, Sina

Based on previous homogenization results for imperfect transmission problems in two-component domains with periodic microstructure, we derive quantitative estimates for the difference between the microscopic and macroscopic solution. This difference is of order , where > 0 describes the periodicity of the microstructure and 2 (0; 1/2 ] depends on the transmission condition at the interface between the two components. The corrector estimates are proved without assuming additional regularity for the local correctors using the periodic unfolding method.

2012, Jenkins, Eleanor W., John, Volker, Linke, Alexander, Rebholz, Leo G.

Grad-div stabilization has been proved to be a very useful tool in discretizations of incompressible flow problems. Standard error analysis for inf-sup stable conforming pairs of finite element spaces predicts that the stabilization parameter should be optimally chosen to be O(1). This paper revisits this choice for the Stokes equations on the basis of minimizing the H1( ) error of the velocity and the L2( ) error of the pressure. It turns out, by applying a refined error analysis, that the optimal parameter choice is more subtle than known so far in the literature. It depends on the used norm, the solution, the family of finite element spaces, and the type of mesh. Depending on the situation, the optimal stabilization parameter might range from being very small to very large. The analytic results are supported by numerical examples.

2015, Colli, Pierluigi, Gilardi, Gianni, Rocca, Elisabetta, Sprekels, Jürgen

In this paper we perform an asymptotic analysis for two different vanishing viscosity coefficients occurring in a phase field system of Cahn--Hilliard type that was recently introduced in order to approximate a tumor growth model. In particular, we extend some recent results obtained in [Colli-Gilardi-Hilhorst 2015], letting the two positive viscosity parameters tend to zero independently from each other and weakening the conditions on the initial data in such a way as to maintain the nonlinearities of the PDE system as general as possible. Finally, under proper growth conditions on the interaction potential, we prove an error estimate leading also to the uniqueness result for the limit system.

2012, Barrenechea, Gabriel R., John, Volker, Knobloch, Petr

An extension of the local projection stabilization (LPS) finite element method for convection-diffusion-reaction equations is presented and analyzed, both in the steady-state and the transient setting. In addition to the standard LPS method, a nonlinear crosswind diffusion term is introduced that accounts for the reduction of spurious oscillations. The existence of a solution can be proved and, depending on the choice of the stabilization parameter, also its uniqueness. Error estimates are derived which are supported by numerical studies. These studies demonstrate also the reduction of the spurious oscillations.

2015, Ahmed, Naveed

This paper studies the parameter choice in the grad-div stabilization applied to the generalized problems of Oseen type. Stabilization parameters based on minimizing the H1 (Omega) error of the velocity are derived which do not depend on the viscosity parameter. For the proposed parameter choices, the H1 (Omega) error of the velocity is derived that shows a direct dependence on the viscosity parameter. Differences and common features to the situation for the Stokes equations are discussed. Numerical studies are presented which confirm the theoretical results. Moreover, for the Navier-Stokes equations, numerical simulations were performed on a two-dimensional flow past a circular cylinder. It turns out, for the MINI element, that the best results can be obtained without grad-div stabilization.

2014, Reichelt, Sina

We derive quantitative error estimates for coupled reaction-diffusion systems, whose coefficient functions are quasi-periodically oscillating modeling microstructure of the underlying macroscopic domain. The coupling arises via nonlinear reaction terms, and we allow for different diffusion length scales, i.e. whereas some species have characteristic diffusion length of order 1, other species may diffuse much slower, namely, with order of the characteristic microstructure-length scale. We consider an effective system, which is rigorously obtained via two-scale convergence, and we prove that the error of its solution to the original solution is of order 1/2.