The conforming virtual element method for polyharmonic and elastodynamics problems: a review

In this paper, we review recent results on the conforming virtual element approximation of polyharmonic and elastodynamics problems. The structure and the content of this review is motivated by three paradigmatic examples of applications: classical and anisotropic Cahn-Hilliard equation and phase field models for brittle fracture, that are briefly discussed in the first part of the paper. We present and discuss the mathematical details of the conforming virtual element approximation of linear polyharmonic problems, the classical Cahn-Hilliard equation and linear elastodynamics problems.Comment: 30 pages, 7 figures. arXiv admin note: text overlap with arXiv:1912.0712

On the Virtual Element Method for Topology Optimization on polygonal meshes: a numerical study

It is well known that the solution of topology optimization problems may be affected both by the geometric properties of the computational mesh, which can steer the minimization process towards local (and non-physical) minima, and by the accuracy of the method employed to discretize the underlying differential problem, which may not be able to correctly capture the physics of the problem. In light of the above remarks, in this paper we consider polygonal meshes and employ the virtual element method (VEM) to solve two classes of paradigmatic topology optimization problems, one governed by nearly-incompressible and compressible linear elasticity and the other by Stokes equations. Several numerical results show the virtues of our polygonal VEM based approach with respect to more standard methods

On a non-isothermal diffuse interface model for two-phase flows of incompressible fluids

We introduce a diffuse interface model describing the evolution of a mixture of two different viscous incompressible fluids of equal density. The main novelty of the present contribution consists in the fact that the effects of temperature on the flow are taken into account. In the mathematical model, the evolution of the macroscopic velocity is ruled by the Navier-Stokes system with temperature-dependent viscosity, while the order parameter representing the concentration of one of the components of the fluid is assumed to satisfy a convective Cahn-Hilliard equation. The effects of the temperature are prescribed by a suitable form of the heat equation. However, due to quadratic forcing terms, this equation is replaced, in the weak formulation, by an equality representing energy conservation complemented with a differential inequality describing production of entropy. The main advantage of introducing this notion of solution is that, while the thermodynamical consistency is preserved, at the same time the energy-entropy formulation is more tractable mathematically. Indeed, global-in-time existence for the initial-boundary value problem associated to the weak formulation of the model is proved by deriving suitable a-priori estimates and showing weak sequential stability of families of approximating solutions.Comment: 26 page

Virtual Element Methods for hyperbolic problems on polygonal meshes

In the present paper we develop the Virtual Element Method for hyperbolic problems on polygonal meshes, considering the linear wave equations as our model problem. After presenting the semi-discrete scheme, we derive the convergence estimates in H^1 semi-norm and L^2 norm. Moreover we develop a theoretical analysis on the stability for the fully discrete problem by comparing the Newmark method and the Bathe method. Finally we show the practical behaviour of the proposed method through a large array of numerical tests

The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries

We explore diffuse formulations of Nitsche's method for consistently imposing Dirichlet boundary conditions on phase-field approximations of sharp domains. Leveraging the properties of the phase-field gradient, we derive the variational formulation of the diffuse Nitsche method by transferring all integrals associated with the Dirichlet boundary from a geometrically sharp surface format in the standard Nitsche method to a geometrically diffuse volumetric format. We also derive conditions for the stability of the discrete system and formulate a diffuse local eigenvalue problem, from which the stabilization parameter can be estimated automatically in each element. We advertise metastable phase-field solutions of the Allen-Cahn problem for transferring complex imaging data into diffuse geometric models. In particular, we discuss the use of mixed meshes, that is, an adaptively refined mesh for the phase-field in the diffuse boundary region and a uniform mesh for the representation of the physics-based solution fields. We illustrate accuracy and convergence properties of the diffuse Nitsche method and demonstrate its advantages over diffuse penalty-type methods. In the context of imaging based analysis, we show that the diffuse Nitsche method achieves the same accuracy as the standard Nitsche method with sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human vertebral body