A Variational r-Adaption and Shape-Optimization Method for Finite-Deformation Elasticity

This paper is concerned with the formulation of a variational r-adaption method for finite-deformation elastostatic problems. The distinguishing characteristic of the method is that the variational principle simultaneously supplies the solution, the optimal mesh and, in problems of shape optimization, the equilibrium shapes of the system. This is accomplished by minimizing the energy functional with respect to the nodal field values as well as with respect to the triangulation of the domain of analysis. Energy minimization with respect to the referential nodal positions has the effect of equilibrating the energetic or configurational forces acting on the nodes. We derive general expressions for the configuration forces for isoparametric elements and nonlinear, possibly anisotropic, materials under general loading. We illustrate the versatility and convergence characteristics of the method by way of selected numerical tests and applications, including the problem of a semi-infinite crack in linear and nonlinear elastic bodies; and the optimization of the shape of elastic inclusions

Numerical approximation of phase field based shape and topology optimization for fluids

We consider the problem of finding optimal shapes of fluid domains. The fluid obeys the Navier--Stokes equations. Inside a holdall container we use a phase field approach using diffuse interfaces to describe the domain of free flow. We formulate a corresponding optimization problem where flow outside the fluid domain is penalized. The resulting formulation of the shape optimization problem is shown to be well-posed, hence there exists a minimizer, and first order optimality conditions are derived. For the numerical realization we introduce a mass conserving gradient flow and obtain a Cahn--Hilliard type system, which is integrated numerically using the finite element method. An adaptive concept using reliable, residual based error estimation is exploited for the resolution of the spatial mesh. The overall concept is numerically investigated and comparison values are provided

A Phase Field Model for Continuous Clustering on Vector Fields

A new method for the simplification of flow fields is presented. It is based on continuous clustering. A well-known physical clustering model, the Cahn Hilliard model, which describes phase separation, is modified to reflect the properties of the data to be visualized. Clusters are defined implicitly as connected components of the positivity set of a density function. An evolution equation for this function is obtained as a suitable gradient flow of an underlying anisotropic energy functional. Here, time serves as the scale parameter. The evolution is characterized by a successive coarsening of patterns-the actual clustering-during which the underlying simulation data specifies preferable pattern boundaries. We introduce specific physical quantities in the simulation to control the shape, orientation and distribution of the clusters as a function of the underlying flow field. In addition, the model is expanded, involving elastic effects. In the early stages of the evolution shear layer type representation of the flow field can thereby be generated, whereas, for later stages, the distribution of clusters can be influenced. Furthermore, we incorporate upwind ideas to give the clusters an oriented drop-shaped appearance. Here, we discuss the applicability of this new type of approach mainly for flow fields, where the cluster energy penalizes cross streamline boundaries. However, the method also carries provisions for other fields as well. The clusters can be displayed directly as a flow texture. Alternatively, the clusters can be visualized by iconic representations, which are positioned by using a skeletonization algorithm.