Serendipity Virtual Elements for General Elliptic Equations in Three Dimensions

The authors study the use of the virtual element method(VEM for short) of order k for general second order elliptic problems with variable coefficients in three space dimensions. Moreover, they investigate numerically also the serendipity version of the VEM and the associated computational gain in terms of degrees of freedom.

Virtual Element Discretization of Optimal Control Problem Governed by Brinkman Equations

In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L2 and H1 norm are derived. The theoretical findings are illustrated by the numerical experiments.

A Priori Error Analysis for NCVEM Discretization of Elliptic Optimal Control Problem

In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H1 and L2 norms. Finally, some numerical experiments are carried out to support the theoretical results.

Shape optimization of plate with static and dynamic constraints via virtual laminated element

The virtual laminated element method (VLEM) can resolve structural shap e optimization problems with a new method. According to the characteristics of V LEM , only some characterized layer thickness values need be defined as design v ariables instead of boundary node coordinates or some other parameters determini ng the system boundary. One of the important features of this method is that it is not necessary to regenerate the FE(finite element) grid during the optimizati on process so as to avoid optimization failures resulting from some distortion grid elements. Th e thickness distribution in thin plate optimization problems in other studies be fore is of stepped shape. However, in this paper, a continuous thickness distrib ution can be obtained after optimization using VLEM, and is more reasonable. Fur thermore, an approximate reanalysis method named ″behavior model technique″ ca n be used to reduce the amount of structural reanalysis. Some typical examples are offered to prove the effectiveness and practicality of the proposed method.

Virtual Element Formulation for Finite Strain Elastodynamics Dedicated to Professor Karl Stark Pister for his 95th birthday

The virtual element method(VEM)can be seen as an extension of the classical finite element method(FEM)based on Galerkin projection.It allows meshes with highly irregular shaped elements,including concave shapes.So far the virtual element method has been applied to various engineering problems such as elasto-plasticity,multiphysics,damage and fracture mechanics.This work focuses on the extension of the virtual element method to efficient modeling of nonlinear elasto-dynamics undergoing large deformations.Within this framework,we employ low-order ansatz functions in two and three dimensions for elements that can have arbitrary polygonal shape.The formulations considered in this contribution are based on minimization of potential function for both the static and the dynamic behavior.Generally the construction of a virtual element is based on a projection part and a stabilization part.While the stiffness matrix needs a suitable stabilization,the mass matrix can be calculated using only the projection part.For the implicit time integration scheme,Newmark-Method is used.To show the performance of the method,various two-and three-dimensional numerical examples in are presented.

Allen-Cahn方程的虚拟元方法

In this article,the virtual element method of the Allen-Cahn equation on a polygon grid is discussed in the fully discrete formulation.With the help of the energy projection operator,we give the corresponding error estimates in the L2 norm and H1 norm.

ERROR ANALYSIS OF THE SECOND-ORDER SERENDIPITY VIRTUAL ELEMENT METHOD FOR SEMILINEAR PSEUDO-PARABOLIC EQUATIONS ON CURVED DOMAINS

The second-order serendipity virtual element method is studied for the semilinear pseudo-parabolic equations on curved domains in this paper.Nonhomogeneous Dirichlet boundary conditions are taken into account,the existence and uniqueness are investigated for the weak solution of the nonhomogeneous initial-boundary value problem.The Nitschebased projection method is adopted to impose the boundary conditions in a weak way.The interpolation operator is used to deal with the nonlinear term.The Crank-Nicolson scheme is employed to discretize the temporal variable.There are two main features of the proposed scheme:(i)the internal degrees of freedom are avoided no matter what type of mesh is utilized,and(ii)the Jacobian is simple to calculate when Newton’s iteration method is applied to solve the fully discrete scheme.The error estimates are established for the discrete schemes and the theoretical results are illustrated through some numerical examples.

Divergence-Free Virtual Element Method for the Stokes Equations with Damping on Polygonal Meshes

In this paper,we construct,analyze,and numerically validate a family of divergence-free virtual elements for Stokes equations with nonlinear damping on polygonal meshes.The virtual element method is H1-conforming and exact divergence-free.By virtue of these properties and the topological degree argument,we rigorously prove the well-posedness of the proposed discrete scheme.The con-vergence analysis is carried out,which imply that the error estimate for the velocity in energy norm does not explicitly depend on the pressure.Numerical experiments on various polygonal meshes validate the accuracy of the theoretical analysis and the asymptotic pressure robustness of the proposed scheme.

NON-UNIFORM LINEAR ARRAY CONFIGURATION FOR MIMO RADAR

Array configuration of multiple-input multiple-output （MIMO） radar with non-uniform linear array （NLA） is proposed. Unlike a standard phased-array radar where NLA is used to generate thinner beam patterns, in MIMO radar the property of NLA is exploited to get more distinct virtual array elements so as to improve pa- rameter identifiability, which means the maximum number of targets that can be uniquely identified by the radar. A class of NLA called minimum redundancy linear array （MRLA） is employed and a new method to construct large MRLAs is descrihed. The numerical results verify that compared to uniform linear array （ULA） MIMO radars, NLA MIMO radars can retain the same parameter identifiability with fewer physical antennas and achieve larger aperture length and lower Cramer-Rao bound with the same number of the physical antennas.

Three-dimensional nonlinear analysis of creep in concrete filled steel tube columns

This paper proposes a based on 3D-VLE （three-dimensional nonlinear viscoelastic theory） three-parameters viscoelastic model for studying the time-dependent behaviour of concrete filled steel tube （CFT） columns. The method of 3D-VLE was developed to analyze the effects of concrete creep behavior on CFT structures. After the evaluation of the parameters in the proposed creep model, experimental measurements of two prestressed reinforced concrete beams were used to investigate the creep phenomenon of three CFT columns under long-term axial and eccentric load was investigated. The experimentally obtained time-dependent creep behaviour accorded well with the cu~＇es obtained from the proposed method. Many factors （such as ratio of long-term load to strength, slenderness ratio, steel ratio, and eccentricity ratio） were considered to obtain the regularity of influence of concrete creep on CFT structures. The analytical results can be consulted in the engineering practice and design.

The Bulk-Surface Virtual Element Method for Reaction-Diffusion PDEs:Analysis and Applications

Bulk-surface partial differential equations(BS-PDEs)are prevalent in manyapplications such as cellular,developmental and plant biology as well as in engineeringand material sciences.Novel numerical methods for BS-PDEs in three space dimensions(3D)are sparse.In this work,we present a bulk-surface virtual elementmethod(BS-VEM)for bulk-surface reaction-diffusion systems,a form of semilinearparabolic BS-PDEs in 3D.Unlike previous studies in two space dimensions(2D),the3D bulk is approximated with general polyhedra,whose outer faces constitute a flatpolygonal approximation of the surface.For this reason,the method is restricted tothe lowest order case where the geometric error is not dominant.The BS-VEM guaranteesall the advantages of polyhedral methods such as easy mesh generation andfast matrix assembly on general geometries.Such advantages are much more relevantthan in 2D.Despite allowing for general polyhedra,general nonlinear reaction kineticsand general surface curvature,the method only relies on nodal values without needingadditional evaluations usually associated with the quadrature of general reactionkinetics.This latter is particularly costly in 3D.The BS-VEM as implemented in thisstudy retains optimal convergence of second order in space.

UNIFIED WAY FOR DEALING WITH THREE-DIMENSIONAL PROBLEMS OF SOLID ELASTICITY

Unified way for dealing with the problems of three dimensional solid, each type of plates and shells etc. was presented with the virtual boundary element least, squares method(VBEM). It proceeded from the differential equations of three-dimensional theory of elasticity and employs the Kelvin solution and the least squares method, It is advantageous to the establishment of the models of a software for general application to calculate each type of three-dimensional problems of elasticity. Owing to directly employing the Kelvin solution and not citing any hypothesis, the numerical results of the method should be better than any others. The merits of the method are highlighted in comparison with the direct formulation of boundary element method (BEM). It is shown that coefficient matrix is symmetric and the treatment of singular integration is rendered unnecessary in the presented method. The examples prove the efficiency and calculating precision of the method.

A New Locking-Free Virtual Element Method for Linear Elasticity Problems

This paper devises a new lowest-order conforming virtual element method(VEM)for planar linear elasticity with the pure displacement/traction boundary condition.The main trick is to view a generic polygon K as a new one K with additional vertices consisting of interior points on edges of K,so that the discrete admissible space is taken as the V1 type virtual element space related to the partition{K}instead of{K}.The method is proved to converge with optimal convergence order both in H^(1)and L^(2)norms and uniformly with respect to the Lam´e constantλ.Numerical tests are presented to illustrate the good performance of the proposed VEM and confirm the theoretical results.

CONSERVATIVE CONFORMING AND NONCONFORMING VEMS FOR FOURTH ORDER NONLINEAR SCHRODINGER EQUATIONS WITH TRAPPED TERM

This paper aims to construct and analyze the conforming and nonconforming virtual element methods for a class of fourth order nonlinear Schrodinger equations with trapped term.We mainly consider three types of virtual elements,including H^(2) conforming virtual element,C^(0) nonconforming virtual element and Morley-type nonconforming virtual element.The fully discrete schemes are constructed by virtue of virtual element methods in space and modified Crank-Nicolson method in time.We prove the mass and energy conservation,the boundedness and the unique solvability of the fully discrete schemes.After introducing a new type of the Ritz projection,the optimal and unconditional error estimates for the fully discrete schemes are presented and proved.Finally,two numerical examples are investigated to confirm our theoretical analysis.

UNCONDITIONAL ERROR ANALYSIS OF VEMS FOR A GENERALIZED NONLINEAR SCHRODINGER EQUATION

In this work,we focus on the conforming and nonconforming leap-frog virtual element methods for the generalized nonlinear Schrodinger equation,and establish their unconditional stability and optimal error estimates.By constructing a time-discrete system,the error between the solutions of the continuous model and the numerical scheme is separated into the temporal error and the spatial error,which makes the spatial error τ-independent.The inverse inequalities in the existing conforming and new constructed nonconforming virtual element spaces are utilized to derive the L^(∞)-norm uniform boundedness of numerical solutions without any restrictions on time-space step ratio,and then unconditionally optimal error estimates of the numerical schemes are obtained naturally.What needs to be emphasized is that if we use the pre-existing nonconforming virtual elements,there is no way to derive the L^(∞)-norm uniform boundedness of the functions in the nonconforming virtual element spaces so as to be hard to get the corresponding inverse inequalities.Finally,several numerical examples are reported to confirm our theoretical results.

A MIXED VIRTUAL ELEMENT METHOD FOR THE BOUSSINESQ PROBLEM ON POLYGONAL MESHES

In this work we introduce and analyze a mixed virtual element method(mixed-VEM)for the two-dimensional stationary Boussinesq problem.The continuous formulation is based on the introduction of a pseudostress tensor depending nonlinearly on the velocity,which allows to obtain an equivalent model in which the main unknowns are given by the aforementioned pseudostress tensor,the velocity and the temperature,whereas the pressure is computed via a postprocessing formula.In addition,an augmented approach together with a fixed point strategy is used to analyze the well-posedness of the resulting continuous formulation.Regarding the discrete problem,we follow the approach employed in a previous work dealing with the Navier-Stokes equations,and couple it with a VEM for the convection-diffusion equation modelling the temperature.More precisely,we use a mixed-VEM for the scheme associated with the fluid equations in such a way that the pseudostress and the velocity are approximated on virtual element subspaces of H(div)and H^(1),respectively,whereas a VEM is proposed to approximate the temperature on a virtual element subspace of H^(1).In this way,we make use of the L^(2)-orthogonal projectors onto suitable polynomial spaces,which allows the explicit integration of the terms that appear in the bilinear and trilinear forms involved in the scheme for the fluid equations.On the other hand,in order to manipulate the bilinear form associated to the heat equations,we define a suitable projector onto a space of polynomials to deal with the fact that the diffusion tensor,which represents the thermal conductivity,is variable.Next,the corresponding solvability analysis is performed using again appropriate fixed-point arguments.Further,Strang-type estimates are applied to derive the a priori error estimates for the components of the virtual element solution as well as for the fully computable projections of them and the postprocessed pressure.The corresponding rates of convergence are also established.Finally,several numerical examples illustrating the performance of the mixed-VEM scheme and confirming these theoretical rates are presented.