Search Results

Now showing 1 - 6 of 6
  • Item
    Anisotropic finite element mesh adaptation via higher dimensional embedding
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 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-4] to obtain an anisotropic curvature adapted mesh that fits a complex surface in ℝ3. In the context of adaptive finite element simulation, the solution (which is an unknown function ƒ: Ω ⊂; ℝ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 ƒ 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 ƒ. 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.
  • Item
    WPM package manager version 1.0 : software documentation
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Streckenbach, Timo
    WPM is a command-line tool designed to support build and installation facilities. It is implemented as a collection of script files, written in Bourne shell syntax. For the sake of portability the code takes care of the common pitfalls of shell programming.
  • Item
    Inverse modeling of thin layer flow cells for detection of solubility, transport and reaction coefficients from experimental data
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Fuhrmann, Jürgen; Linke, Alexander; Merdon, Christian; Neumann, Felix; Streckenbach, Timo; Baltruschat, Helmut; Khodayari, Mehdi
    Thin layer flow cells are used in electrochemical research as experimental devices which allow to perform investigations of electrocatalytic surface reactions under controlled conditions using reasonably small electrolyte volumes. The paper introduces a general approach to simulate the complete cell using accurate numerical simulation of the coupled flow, transport and reaction processes in a flow cell. The approach is based on a mass conservative coupling of a divergence-free finite element method for fluid flow and a stable finite volume method for mass transport. It allows to perform stable and efficient forward simulations that comply with the physical bounds namely mass conservation and maximum principles for the involved species. In this context, several recent approaches to obtain divergence-free velocities from finite element simulations are discussed. In order to perform parameter identification, the forward simulation method is coupled to standard optimization tools. After an assessment of the inverse modeling approach using known real-istic data, first results of the identification of solubility and transport data for O2 dissolved in organic electrolytes are presented. A plausibility study for a more complex situation with surface reactions concludes the paper and shows possible extensions of the scope of the presented numerical tools.
  • Item
    Parameter determination of an evolution model for phase transformations in steel
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Streckenbach, Timo
    The paper is concerned with a general ansatz of a phenomenological evolution model for solid-solid phase transformation kinetics in steel. To model the phase transition of austenite-ferrite, -pearlite or -bainite, a first order nonlinear ordinary differential equation (ODE) is considered. The main goal of this paper is to derive certain conditions for parameters which based on data obtained from transformation diagrams. This leads to a set of independent parameters for which the inverse problem has an unique solution
  • Item
    Numerical simulation of TEM images for In(Ga)As/GaAs quantum dots with various shapes
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Maltsi, Anieza; Niermann, Tore; Streckenbach, Timo; Tabelow, Karsten; Koprucki, Thomas
    We present a mathematical model and a tool chain for the numerical simulation of TEM images of semiconductor quantum dots (QDs). This includes elasticity theory to obtain the strain profile coupled with the Darwin-Howie-Whelan equations, describing the propagation of the electron wave through the sample. We perform a simulation study on indium gallium arsenide QDs with different shapes and compare the resulting TEM images to experimental ones. This tool chain can be applied to generate a database of simulated TEM images, which is a key element of a novel concept for model-based geometry reconstruction of semiconductor QDs, involving machine learning techniques.
  • Item
    Locking free and gradient robust H(div)-conforming HDG methods for linear elasticity
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Fu, Guosheng; Lehrenfeld, Christoph; Linke, Alexander; Streckenbach, Timo
    Robust discretization methods for (nearly-incompressible) linear elasticity are free of volume-locking and gradient-robust. While volume-locking is a well-known problem that can be dealt with in many different discretization approaches, the concept of gradient-robustness for linear elasticity is new. We discuss both aspects and propose novel Hybrid Discontinuous Galerkin (HDG) methods for linear elasticity. The starting point for these methods is a divergence-conforming discretization. As a consequence of its well-behaved Stokes limit the method is gradient-robust and free of volume-locking. To improve computational efficiency, we additionally consider discretizations with relaxed divergence-conformity and a modification which re-enables gradient-robustness, yielding a robust and quasi-optimal discretization also in the sense of HDG superconvergence.