Refine
Language
- English (19)
Keywords
- (visco-)elasticity (1)
- 3D mesh compression (1)
- A posteriori, error analysis, finite element method, composite material, multiscale, high contrast, generalized Delaunay, network (1)
- C-stationarity (1)
- Cahn-Hilliard system (1)
- Cahn-Hilliard/Navier-Stokes system (1)
- Context-Based Coding (1)
- Cosserat rod (1)
- Dynamical contact problems (1)
- Existence result (1)
This paper is intended to be a first step towards the continuous dependence of dynamical contact problems on the initial data as well as the uniqueness of a solution. Moreover, it provides the basis for a proof of the convergence of popular time integration schemes as the Newmark method.
We study a frictionless dynamical contact problem between both linearly elastic and viscoelastic bodies which is formulated via the Signorini contact conditions. For viscoelastic materials fulfilling the Kelvin-Voigt constitutive law, we find a characterization of the class of problems which satisfy a perturbation result in a non-trivial mix of norms in function space. This characterization is given in the form of a stability condition on the contact stresses at the contact boundaries.
Furthermore, we present perturbation results for two well-established approximations of the classical Signorini condition: The Signorini condition formulated in velocities and the model of normal compliance, both satisfying even a sharper version of our stability condition.
{We suggest hierarchical a posteriori error estimators for time-discretized
Allen-Cahn and Cahn-Hilliard equations with logarithmic potential and investigate
their robustness numerically.
We observe that the associated effectivity ratios seem to saturate for decreasing mesh size
and are almost independent of the temperature.
Motivated by an obstacle problem for a membrane
subject to cohesion forces, constrained minimization problems involving
a non-convex and non-differentiable objective functional
representing the total potential energy are considered. The associated
first order optimality system leads to a hemi-variational inequality,
which can also be interpreted as a special complementarity
problem in function space. Besides an analytical investigation
of first-order optimality, a primal-dual active set solver is introduced.
It is associated to a limit case of a semi-smooth Newton
method for a regularized version of the underlying problem class.
For the numerical algorithms studied in this paper, global as well as
local convergence properties are derived and verified numerically.
The aim of this paper is to devise an adaptive timestep control in the contact--stabilized Newmark method (CONTACX) for dynamical contact problems between two viscoelastic bodies in the framework of Signorini's condition. In order to construct a comparative scheme of higher order accuracy, we extend extrapolation techniques. This approach demands a subtle theoretical investigation of an asymptotic error expansion of the contact--stabilized Newmark scheme. On the basis of theoretical insight and numerical observations, we suggest an error estimator and a timestep selection which also cover the presence of contact. Finally, we give a numerical example.
This note constructs a local generalized finite element basis for elliptic problems with heterogeneous and highly varying diffusion tensor. The basis functions are solutions of local problems on vertex patches. The error of the corresponding generalized finite element method decays exponentially w.r.t. the number of element layers in the patches. Hence, on a uniform mesh of size H, patches of diameter log(1/H) are sufficient to preserve the convergence rates of the classical P1-FEM for the Poisson problem.
The analysis does not rely on regularity of the solution or scale separation in the coefficient.
The result justifies the use of the class of variational multiscale methods, introduced in [Comput. Methods Appl. Mech. Engrg., 196:2313--2324, 2007].
An adaptive a posteriori error estimator based finite element method for the
numerical solution of a coupled Cahn-Hilliard/Navier-Stokes system with a
double-obstacle homogenous free (interfacial) energy density is proposed. A
semi-implicit Euler scheme for the time-integration is applied which results in
a system coupling a quasi-Stokes or Oseen-type problem for the fluid flow to
a variational inequality for the concentration and the chemical potential according
to the Cahn-Hilliard model [13]. A Moreau-Yosida regularization is
employed which relaxes the constraints contained in the variational inequality
and, thus, enables semi-smooth Newton solvers with locally superlinear convergence
in function space. Moreover, upon discretization this yields a mesh
independent method for a fixed relaxation parameter. For the finite dimensional
approximation of the concentration and the chemical potential piecewise
linear and globally continuous finite elements are used, and for the numerical
approximation of the fluid velocity Taylor-Hood finite elements are employed.
The paper ends by a report on numerical examples showing the efficiency of the
new method
Multiresolution meshes provide an efficient and structured representation of geometric objects. To increase the
mesh resolution only at vital parts of the object, adaptive refinement is widely used. We propose a lossless compression
scheme for these adaptive structures that exploits the parent-child relationships inherent to the mesh
hierarchy. We use the rules that correspond to the adaptive refinement scheme and store bits only where some
freedom of choice is left, leading to compact codes that are free of redundancy. Moreover, we extend the coder to
sequences of meshes with varying refinement. The connectivity compression ratio of our method exceeds that of
state-of-the-art coders by a factor of 2 to 7.
For efficient compression of vertex positions we adapt popular wavelet-based coding schemes to the adaptive
triangular and quadrangular cases to demonstrate the compatibility with our method. Akin to state-of-the-art
coders, we use a zerotree to encode the resulting coefficients. Using improved context modeling we enhanced the
zerotree compression, cutting the overall geometry data rate by 7% below those of the successful Progressive
Geometry Compression. More importantly, by exploiting the existing refinement structure we achieve compression
factors that are 4 times greater than those of coders which can handle irregular meshes.
We introduce hexagonal global parameterizations, a new
type of surface parameterizations in which parameter lines respect six-fold rotational symmetries (6-RoSy). Such parameterizations enable the tiling of surfaces with nearly regular hexagonal or triangular patterns, and can be used for triangular remeshing.
To construct a hexagonal parameterization on a surface, we provide an automatic algorithm to generate a 6-RoSy field that respects directional and singularity features of the surface. This field is then used to direct a hexagonal global parameterization. The framework, called HexCover, extends the QuadCover algorithm and formulates necessary conditions for hexagonal parameterization.
We demonstrate the usefulness of our geometry-aware global parameterization with applications such as surface tiling with nearly regular textures and geometry patterns, as well as triangular and hexagonal remeshing.
Despite the success of quad-based 2D surface parameterization methods, effective parameterization algorithms for 3D volumes with cubes, i.e. hexahedral elements, are still missing. CubeCover is a first approach which provides both, a consistent theoretical framework for volume parameterization plus a full pipeline for generating
a hexahedral tessellation of a given volume with boundary aligned cubes which are guided by a frame field.
The input of CubeCover is a tetrahedral volume mesh. First, a frame field is designed with manual input from
the designer. It guides the interior and boundary layout of the parameterization. Then, the parameterization and
the hexahedral mesh are computed so as to align with the given frame field.
CubeCover has similarities to the QuadCover algorithm and extends it from 2D surfaces to 3D volumes. The
paper also provides theoretical results for 3D hexahedral parameterizations and analyses topological properties
of the appropriate function space.
In this paper we study the distributed optimal control for the Cahn-Hilliard system. A general class of free energy potentials is allowed which, in particular, includes the double-obstacle potential. The latter potential yields an optimal control problem of a parabolic variational inequality which is of fourth-order in space. We show the existence of optimal controls to approximating problems where the potential is replaced by a mollified version of its Moreau-Yosida approximation. Corresponding first order optimality conditions for the mollified problems are given. For this purpose a new result on the continuous Fr\'echet differentiability of superposition operators with values in Sobolev spaces is established. Besides the convergence of optimal controls of the mollified problems to an optimal control of the original problem, we also derive first order optimality conditions for the original problem by a limit process. The newly derived stationarity system corresponds to a function space version of C-stationarity.