Search Results

Now showing 1 - 10 of 1371
  • Item
    Classification and clustering: models, software and applications
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Mucha, Hans-Joachim; Ritter, Gunter
    We are pleased to present the report on the 30th Fall Meeting of the working group ``Data Analysis and Numerical Classification'' (AG-DANK) of the German Classification Society. The meeting took place at the Weierstrass Institute for Applied Analysis and Stochastics (WIAS), Berlin, from Friday Nov. 14 till Saturday Nov. 15, 2008. Already 12 years ago, WIAS had hosted a traditional Fall Meeting with special focus on classification and multivariate graphics (Mucha and Bock, 1996). This time, the special topics were stability of clustering and classification, mixture decomposition, visualization, and statistical software.
  • Item
    Optimal and robust a posteriori error estimates in L∞(L2) for the approximation of Allen-Cahn equations past singularities
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Bartels, Sören; Müller, Rüdiger
    Optimal a posteriori error estimates in L∞(L2) are derived for the finite element approximation of Allen-Cahn equations. The estimates depend on the inverse of a small parameter only in a low order polynomial and are valid past topological changes of the evolving interface. The error analysis employs an elliptic reconstruction of the approximate solution and applies to a large class of conforming, nonconforming, mixed, and discontinuous Galerkin methods. Numerical experiments illustrate the theoretical results.
  • Item
    Impact of slippage on the morphology and stability of a dewetting rim
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Münch, Andreas; Wagner, Barbara
    In this study lubrication theory is used to describe the stability and morphology of the rim that forms as a thin polymer film dewets from a hydrophobized silicon wafer. Thin film equations are derived from the governing hydrodynamic equations for the polymer to enable the systematic mathematical and numerical analysis of the properties of the solutions for different regimes of slippage and for a range of time scales. Dewetting rates and the cross sectional profiles of the evolving rims are derived for these models and compared to experimental results. Experiments also show that the rim is typically unstable in the spanwise direction and develops thicker and thinner parts that may grow into ``fingers''. Linear stability analysis as well as nonlinear numerical solutions are presented to investigate shape and growth rate of the rim instability. It is demonstrated that the difference in morphology and the rate at which the instability develops can be directly attributed to the magnitude of slippage. Finally, a derivation is given for the dominant wavelength of the bulges along the unstable rim.
  • 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
    Error control for the approximation of Allen-Cahn and Cahn-Hilliard equations with a logarithmic potential
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Bartels, Sören; Müller, Rüdiger
    A fully computable upper bound for the finite element approximation error of Allen-Cahn and Cahn-Hilliard equations with logarithmic potentials is derived. Numerical experiments show that for the sharp interface limit this bound is robust past topological changes. Modifications of the abstract results to derive quasi-optimal error estimates in different norms for lowest order finite element methods are discussed and lead to weaker conditions on the residuals under which the conditional error estimates hold.
  • Item
    Anisotropic growth of random surfaces in 2 + 1 dimensions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2008) Borodin, Alexei; Ferrari, Patrik L.
    We construct a family of stochastic growth models in $2+1$ dimensions, that belong to the anisotropic KPZ class. Appropriate projections of these models yield $1+1$ dimensional growth models in the KPZ class and random tiling models. We show that correlation functions associated to our models have determinantal structure, and we study large time asymptotics for one of the models. The main asymptotic results are: (1) The growing surface has a limit shape that consists of facets interpolated by a curved piece. (2) The one-point fluctuations of the height function in the curved part are asymptotically normal with variance of order $ln(t)$ for time $tgg 1$. (3) There is a map of the $(2+1)$-dimensional space-time to the upper half-plane $H$ such that on space-like submanifolds the multi-point fluctuations of the height function are asymptotically equal to those of the pullback of the Gaussian free (massless) field on $H$.
  • Item
    On the evolution of the empirical measure for the hard-sphere dynamics
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Pulvirenti, Mario; Simonella, Sergio
    We prove that the evolution of marginals associated to the empirical measure of a finite system of hard spheres is driven by the BBGKY hierarchical expansion. The usual hierarchy of equations for L1 measures is obtained as a corollary. We discuss the ambiguities arising in the corresponding notion of microscopic series solution to the Boltzmann-Enskog equation.
  • Item
    Modeling and simulations of beam stabilization in edge-emitting broad area semiconductor devices
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Radziunas, Mindaugas; Cˇ iegis, Raimondas
    A 2+1 dimensional PDE traveling wave model describing spatial-lateral dynamics of edge-emitting broad area semiconductor devices is considered. A numerical scheme based on a split-step Fourier method is presented and implemented on a parallel compute cluster. Simulations of the model equations are used for optimizing of existing devices with respect to the emitted beam quality, as well as for creating and testing of novel device design concepts
  • Item
    Moment asymptotics for branching random walks in random environment
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Gün, Onur; König, Wolfgang; Sekulov´c, Ozren
    We consider the long-time behaviour of a branching random walk in random environment on the lattice Zd. The migration of particles proceeds according to simple random walk in continuous time, while the medium is given as a random potential of spatially dependent killing/branching rates. The main objects of our interest are the annealed moments m_np , i.e., the p-th moments over the medium of the n-th moment over the migration and killing/branching, of the local and global population sizes. For n = 1, this is well-understood citeGM98, as m_1 is closely connected with the parabolic Anderson model. For some special distributions, citeA00 extended this to ngeq2, but only as to the first term of the asymptotics, using (a recursive version of) a Feynman-Kac formula for m_n. In this work we derive also the second term of the asymptotics, for a much larger class of distributions. In particular, we show that m_n^p m_1^np are asymptotically equal, up to an error e^o(t). The cornerstone of our method is a direct Feynman-Kac-type formula for mn, which we establish using the spine techniques developed in citeHR1.1
  • Item
    From nonlinear to linear elasticity in a coupled rate-dependent/independent system for brittle delamination
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Rossi, Riccarda; Thomas, Marita
    We revisit the weak, energetic-type existence results obtained in [RT15] for a system for rateindependent, brittle delamination between two visco-elastic, physically nonlinear bulk materials and explain how to rigorously extend such results to the case of visco-elastic, linearly elastic bulk materials. Our approximation result is essentially based on deducing the MOSCO-convergence of the functionals involved in the energetic formulation of the system. We apply this approximation result in two different situations at small strains: Firstly, to pass from a nonlinearly elastic to a linearly elastic, brittle model on the time-continuous level, and secondly, to pass from a time-discrete to a time-continuous model using an adhesive contact approximation of the brittle model, in combination with a vanishing, super-quadratic regularization of the bulk energy. The latter approach is beneficial if the model also accounts for the evolution of temperature.