Search Results

Now showing 1 - 5 of 5
Loading...
Thumbnail Image
Item

Anisotropic Finite Element Mesh Adaptation via Higher Dimensional Embedding

2015, Dassi, Franco, Si, Hang, Perotto, Simona, Streckenbach, Timo

In this paper we provide a novel anisotropic mesh adaptation technique for adaptive finite element analysis. It is based on the concept of higher dimensional embedding, which was exploited in [1], [2], [3], [4] to obtain an anisotropic curvature adapted mesh that fits a complex surface in R3. In the context of adaptive finite element simulation, the solution (which is an unknown function f : Ω ⊂ d → ) is sought by iteratively modifying a finite element mesh according to a mesh sizing field described via a (discrete) metric tensor field that is typically obtained through an error estimator. We proposed to use a higher dimensional embedding, Φf (x):= (x1, …, xd, s f (x1, …, xd), s ▿ f (x1, …, xd))t, instead of the mesh sizing field for the mesh adaption. This embedding contains both informations of the function f itself and its gradient. An isotropic mesh in this embedded space will correspond to an anisotropic mesh in the actual space, where the mesh elements are stretched and aligned according to the features of the function f. To better capture the anisotropy and gradation of the mesh, it is necessary to balance the contribution of the components in this embedding. We have properly adjusted Φf (x) for adaptive finite element analysis. To better understand and validate the proposed mesh adaptation strategy, we first provide a series of experimental tests for piecewise linear interpolation of known functions. We then applied this approach in an adaptive finite element solution of partial differential equations. Both tests are performed on two-dimensional domains in which adaptive triangular meshes are generated. We compared these results with the ones obtained by the software BAMG – a metric-based adaptive mesh generator. The errors measured in the L2 norm are comparable. Moreover, our meshes captured the anisotropy more accurately than the meshes of BAMG.

Loading...
Thumbnail Image
Item

The full Keller-Segel model is well-posed on fairly general domains

2016, Horstmann, Dirk, Rehberg, Joachim, Meinlschmidt, Hannes

In this paper we prove the well-posedness of the full Keller-Segel system, a quasilinear strongly coupled reaction-crossdiffusion system, in the spirit that it always admits a unique local-in-time solution in an adequate function space, provided that the initial values are suitably regular. Apparently, there exists no comparable existence result for the full Keller-Segel system up to now. The proof is carried out for general source terms and is based on recent nontrivial elliptic and parabolic regularity results which hold true even on fairly general spatial domains, combined with an abstract solution theorem for nonlocal quasilinear equations by Amann. Nous considèrons le système de Keller et Segel dans son intégralité, un système quasilinéaire à réaction-diffusion fortement couplé. Le résultat principal montre que ce syst`eme est bien posé, cest-à-dire il admet une solution unique existant localement en temps à valeurs dans un espace fonctionnel approprié, pourvu que les valeurs initiales sont réguliers. Apparemment, il nexiste pas encore des résultats comparables. Pour la demonstration, nous utilisons des résultats récents de régularité elliptique et parabolique applicable à des domaines assez générals, combiné avec un théorème abstrait dAmann concernant les équations quasilinéaires non locales.

Loading...
Thumbnail Image
Item

Optimal control for the thermistor problem

2008, Hömberg, Dietmar, Meyer, Christian, Rehberg, Joachim, Ring, Wolfgang

This paper is concerned with the state-constrained optimal control of the two-dimensional thermistor problem, a quasi-linear coupled system of a parabolic and elliptic PDE with mixed boundary conditions. This system models the heating of a conducting material by means of direct current. Existence, uniqueness and continuity for the state system are derived by employing maximal elliptic and parabolic regularity. By similar arguments the linearized state system is discussed, while the adjoint system involving measures is investigated using a duality argument. These results allow to derive first-order necessary conditions for the optimal control problem.

Loading...
Thumbnail Image
Item

Optimal control of the thermistor problem in three spatial dimensions

2016, Meinlschmidt, Hannes, Meyer, Christian, Rehberg, Joachim

This paper is concerned with the state-constrained optimal control of the three-dimensional thermistor problem, a fully quasilinear coupled system of a parabolic and elliptic PDE with mixed boundary conditions. This system models the heating of a conducting material by means of direct current. Local existence, uniqueness and continuity for the state system are derived by employing maximal parabolic regularity in the fundamental theorem of Pr¨uss. Global solutions are addressed, which includes analysis of the linearized state system via maximal parabolic regularity, and existence of optimal controls is shown if the temperature gradient is under control. The adjoint system involving measures is investigated using a duality argument. These results allow to derive first-order necessary conditions for the optimal control problem in form of a qualified optimality system. The theoretical findings are illustrated by numerical results

Loading...
Thumbnail Image
Item

Optimal regularity for elliptic transmission problems including C1 interfaces

2006, Elschner, Johannes, Rehberg, Joachim, Schmidt, Gunther

We prove an optimal regularity result for elliptic operators $-nabla cdot mu nabla:W^1,q_0 rightarrow W^-1,q$ for a $q>3$ in the case when the coefficient function $mu$ has a jump across a $C^1$ interface and is continuous elsewhere. A counterexample shows that the $C^1$ condition cannot be relaxed in general. Finally, we draw some conclusions for corresponding parabolic operators.