54 results on '"Polygonal meshe"'
Search Results
2. Adaptive virtual element methods with equilibrated fluxes
- Author
-
F. Dassi, J. Gedicke, L. Mascotto, Dassi, F, Gedicke, J, and Mascotto, L
- Subjects
Computational Mathematics ,Numerical Analysis ,Hypercircle method ,Applied Mathematics ,FOS: Mathematics ,hp-adaptivity ,Virtual element method ,Mathematics - Numerical Analysis ,Numerical Analysis (math.NA) ,Equilibrated fluxe ,65N12, 65N30, 65N50 ,Polygonal meshe - Abstract
We present an hp-adaptive virtual element method (VEM) based on the hypercircle method of Prager and Synge for the approximation of solutions to diffusion problems. We introduce a reliable and efficient a posteriori error estimator, which is computed by solving an auxiliary global mixed problem. We show that the mixed VEM satisfies a discrete inf-sup condition, with inf-sup constant independent of the discretization parameters. Furthermore, we construct a stabilization for the mixed VEM, with explicit bounds in terms of the local degree of accuracy of the method. The theoretical results are supported by several numerical experiments, including a comparison with the residual a posteriori error estimator. The numerics exhibit the p-robustness of the proposed error estimator. In addition, we provide a first step towards the localized flux reconstruction in the virtual element framework, which leads to an additional reliable a posteriori error estimator that is computed by solving local (cheap-to-solve and parallelizable) mixed problems. We provide theoretical and numerical evidence that the proposed local error estimator suffers from a lack of efficiency.
- Published
- 2022
3. Sharper Error Estimates for Virtual Elements and a Bubble-Enriched Version
- Author
-
Beirao da Veiga, L, Vacca, G, Vacca, G., Beirao da Veiga, L, Vacca, G, and Vacca, G.
- Abstract
In the present contribution we develop a sharper error analysis for the Virtual Element Method, applied to a model elliptic problem, that separates the element boundary and element interior contributions to the error. As a consequence we are able to propose a variant of the scheme that allows one to take advantage of polygons with many edges (such as those composing Voronoi meshes or generated by agglomeration procedures) in order to yield a more accurate discrete solution. The theoretical results are supported by numerical experiments.
- Published
- 2022
4. Adaptive virtual element methods with equilibrated fluxes
- Author
-
Dassi, F, Gedicke, J, Mascotto, L, Dassi, F., Gedicke, J., Mascotto, L., Dassi, F, Gedicke, J, Mascotto, L, Dassi, F., Gedicke, J., and Mascotto, L.
- Abstract
We present an hp-adaptive virtual element method (VEM) based on the hypercircle method of Prager and Synge for the approximation of solutions to diffusion problems. We introduce a reliable and efficient a posteriori error estimator, which is computed by solving an auxiliary global mixed problem. We show that the mixed VEM satisfies a discrete inf-sup condition with inf-sup constant independent of the discretization parameters. Furthermore, we construct a stabilization for the mixed VEM with explicit bounds in terms of the local degree of accuracy of the method. The theoretical results are supported by several numerical experiments, including a comparison with the residual a posteriori error estimator. The numerics exhibit the p-robustness of the proposed error estimator. In addition, we provide a first step towards the localized flux reconstruction in the virtual element framework, which leads to an additional reliable a posteriori error estimator that is computed by solving local (cheap-to-solve and parallelizable) mixed problems. We provide theoretical and numerical evidence that the proposed local error estimator suffers from a lack of efficiency.
- Published
- 2022
5. Extension of the nonconforming Trefftz virtual element method to the Helmholtz problem with piecewise constant wave number
- Author
-
Lorenzo Mascotto, Alexander Pichler, Mascotto, L, and Pichler, A
- Subjects
Degrees of freedom (statistics) ,010103 numerical & computational mathematics ,Space (mathematics) ,01 natural sciences ,Nonconforming virtual element method ,Helmholtz problem ,FOS: Mathematics ,Applied mathematics ,Wavenumber ,Polygon mesh ,Mathematics - Numerical Analysis ,0101 mathematics ,Plane and evanescent wave ,35J05, 65N12, 65N30, 74J20 ,Mathematics ,Numerical Analysis ,Applied Mathematics ,Piecewise constant wave number ,Numerical Analysis (math.NA) ,010101 applied mathematics ,Computational Mathematics ,Range (mathematics) ,Special functions ,Piecewise ,Constant (mathematics) ,Polygonal meshe ,Trefftz methods - Abstract
We extend the nonconforming Trefftz virtual element method introduced in arXiv:1805.05634 to the case of the fluid-fluid interface problem, that is, a Helmholtz problem with piecewise constant wave number. With respect to the original approach, we address two additional issues: firstly, we define the coupling of local approximation spaces with piecewise constant wave numbers, secondly, we enrich such local spaces with special functions capturing the physical behaviour of the solution to the target problem. As these two issues are directly related to an increase of the number of degrees of freedom, we use a reduction strategy inspired by arXiv:1807.11237, which allows to mitigate the growth of the dimension of the approximation space when considering $h$- and $p$-refinements. This renders the new method highly competitive in comparison to other Trefftz and quasi-Trefftz technologies tailored for the Helmholtz problem with piecewise constant wave number. A wide range of numerical experiments, including the $p$-version with quasi-uniform meshes and the $hp$-version with isotropic and anisotropic mesh refinements, is presented., Comment: 23 pages, 15 figures
- Published
- 2020
6. Sharper Error Estimates for Virtual Elements and a Bubble-Enriched Version
- Author
-
L. Beira͂o da Veiga, G. Vacca, Beirao da Veiga, L, and Vacca, G
- Subjects
Numerical Analysis ,Computational Mathematics ,error analysi ,interpolation estimate ,Applied Mathematics ,Virtual Element Method ,polygonal meshe - Abstract
In the present contribution we develop a sharper error analysis for the Virtual Element Method, applied to a model elliptic problem, that separates the element boundary and element interior contributions to the error. As a consequence we are able to propose a variant of the scheme that allows one to take advantage of polygons with many edges (such as those composing Voronoi meshes or generated by agglomeration procedures) in order to yield a more accurate discrete solution. The theoretical results are supported by numerical experiments.
- Published
- 2022
7. A virtual element method for the miscible displacement of incompressible fluids in porous media
- Author
-
Beirao da Veiga, L, Pichler, A, Vacca, G, Beirao da Veiga L., Pichler A., Vacca G., Beirao da Veiga, L, Pichler, A, Vacca, G, Beirao da Veiga L., Pichler A., and Vacca G.
- Abstract
In the present contribution, we construct a virtual element (VE) discretization for the problem of miscible displacement of one incompressible fluid by another, described by a time-dependent coupled system of nonlinear partial differential equations. Our work represents a first study to investigate the premises of virtual element methods (VEM) for complex fluid flow problems. We combine the VEM discretization with a time stepping scheme and develop a complete theoretical analysis of the method under the assumption of a regular solution. The scheme is then tested both on a regular and on a more realistic test case.
- Published
- 2021
8. Vorticity-stabilized virtual elements for the Oseen equation
- Author
-
Beirão da Veiga, L, Dassi, F, Vacca, G, Beirão da Veiga, L., Dassi, F., Vacca, G., Beirão da Veiga, L, Dassi, F, Vacca, G, Beirão da Veiga, L., Dassi, F., and Vacca, G.
- Abstract
In this paper, we extend the divergence-free VEM of [L. Beirão da Veiga, C. Lovadina and G. Vacca, Virtual elements for the Navier-Stokes problem on polygonal meshes, SIAM J. Numer. Anal. 56 (2018) 1210-1242] to the Oseen problem, including a suitable stabilization procedure that guarantees robustness in the convection-dominated case without disrupting the divergence-free property. The stabilization is inspired from [N. Ahmed, G. R. Barrenechea, E. Burman, J. Guzman, A. Linke and C. Merdon, A pressure-robust discretization of Oseen's equation using stabilization in the vorticity equation, SIAM J. Numer. Anal. 59 (2021) 2746-2774] and includes local SUPG-like terms of the vorticity equation, internal jump terms for the velocity gradients, and an additional VEM stabilization. We derive theoretical convergence results that underline the robustness of the scheme in different regimes, including the convection-dominated case. Furthermore, as in the non-stabilized case, the influence of the pressure on the velocity error is moderate, as it appears only through higher-order terms.
- Published
- 2021
9. Equilibrium analysis of an immersed rigid leaflet by the virtual element method
- Author
-
Beirao Da Veiga, L, Canuto, C, Nochetto, R, Vacca, G, Nochetto, RH, Beirao Da Veiga, L, Canuto, C, Nochetto, R, Vacca, G, and Nochetto, RH
- Abstract
We study, both theoretically and numerically, the equilibrium of a hinged rigid leaflet with an attached rotational spring, immersed in a stationary incompressible fluid within a rigid channel. Through a careful investigation of the properties of the domain functional describing the angular momentum exerted by the fluid on the leaflet (which depends on both the leaflet angular position and its thickness), we identify sufficient conditions on the spring stiffness function for the existence (and uniqueness) of equilibrium positions. This study resorts to techniques from shape differential calculus. We propose a numerical technique that exploits the mesh flexibility of the Virtual Element Method (VEM). A (polygonal) computational mesh is generated by cutting a fixed background grid with the leaflet geometry, and the problem is then solved with stable VEM Stokes elements of degrees 1 and 2 combined with a bisection algorithm. We prove quasi-optimal error estimates and present a large array of numerical experiments to document the accuracy and robustness with respect to degenerate geometry of the proposed methodology.
- Published
- 2021
10. SUPG-stabilized virtual elements for diffusion-convection problems: A robustness analysis
- Author
-
Beirao da Veiga, L, Dassi, F, Lovadina, C, Vacca, G, Beirao da Veiga, L, Dassi, F, Lovadina, C, and Vacca, G
- Abstract
The objective of this contribution is to develop a convergence analysis for SUPG-stabilized Virtual Element Methods in diffusion-convection problems that is robust also in the convection dominated regime. For the original method introduced in [Benedetto et al., CMAME 2016] we are able to show an "almost uniform"error bound (in the sense that the unique term that depends in an unfavourable way on the parameters is damped by a higher order mesh-size multiplicative factor). We also introduce a novel discretization of the convection term that allows us to develop error estimates that are fully robust in the convection dominated cases. We finally present some numerical result.
- Published
- 2021
11. p- and hp- virtual elements for the Stokes problem
- Author
-
Chernov, A, Marcati, C, Mascotto, L, Chernov, A, Marcati, C, and Mascotto, L
- Abstract
We analyse the p- and hp-versions of the virtual element method (VEM) for the Stokes problem on polygonal domains. The key tool in the analysis is the existence of a bijection between Poisson-like and Stokes-like VE spaces for the velocities. This allows us to re-interpret the standard VEM for Stokes as a VEM, where the test and trial discrete velocities are sought in Poisson-like VE spaces. The upside of this fact is that we inherit from Beirão da Veiga et al. (Numer. Math. 138(3), 581–613, 2018) an explicit analysis of best interpolation results in VE spaces, as well as stabilization estimates that are explicit in terms of the degree of accuracy p of the method. We prove exponential convergence of the hp-VEM for Stokes problems with regular right-hand sides. We corroborate the theoretical estimates with numerical tests for both the p- and hp-versions of the method.
- Published
- 2021
12. A nonconforming Trefftz virtual element method for the Helmholtz problem: Numerical aspects
- Author
-
Ilaria Perugia, Lorenzo Mascotto, Alexander Pichler, Mascotto, L, Perugia, I, and Pichler, A
- Subjects
Helmholtz equation ,Computer science ,Plane wave ,Computational Mechanics ,General Physics and Astronomy ,Basis function ,010103 numerical & computational mathematics ,Degrees of freedom (mechanics) ,01 natural sciences ,Set (abstract data type) ,Reduction (complexity) ,Nonconforming space ,FOS: Mathematics ,Applied mathematics ,Virtual element method ,Mathematics - Numerical Analysis ,0101 mathematics ,Ill-conditioning ,35J05, 65N12, 65N30, 74J20 ,Mechanical Engineering ,Numerical Analysis (math.NA) ,Computer Science Applications ,010101 applied mathematics ,Mechanics of Materials ,Enhanced Data Rates for GSM Evolution ,Element (category theory) ,Polygonal meshe - Abstract
We discuss the implementation details and the numerical performance of the recently introduced nonconforming Trefftz virtual element method (Mascotto et al., 2018) for the 2D Helmholtz problem. In particular, we present a strategy to significantly reduce the ill-conditioning of the original method; such a recipe is based on an automatic filtering of the basis functions edge by edge, and therefore allows for a notable reduction of the number of degrees of freedom. A widespread set of numerical experiments, including an application to acoustic scattering, the h -, p -, and h p -versions of the method, is presented. Moreover, a comparison with other Trefftz-based methods for the Helmholtz problem shows that this novel approach results in robust and effective performance.
- Published
- 2019
13. p- and hp- virtual elements for the Stokes problem
- Author
-
Carlo Marcati, Lorenzo Mascotto, Alexey Chernov, Chernov, A, Marcati, C, and Mascotto, L
- Subjects
p-and hp-Galerkin methods ,Degree (graph theory) ,Stokes equation ,Virtual element methods ,Polygonal meshes ,p- and hp-Galerkin method ,Exponential convergence ,Applied Mathematics ,Numerical Analysis (math.NA) ,Stokes flow ,Computational Mathematics ,Explicit analysis ,65N12, 65N15, 65N30, 76D07 ,FOS: Mathematics ,Stokes problem ,Bijection ,Applied mathematics ,Virtual element method ,Mathematics - Numerical Analysis ,Element (category theory) ,Polygonal meshe ,Mathematics ,Interpolation - Abstract
We analyse the p- and hp-versions of the virtual element method (VEM) for the Stokes problem on polygonal domains. The key tool in the analysis is the existence of a bijection between Poisson-like and Stokes-like VE spaces for the velocities. This allows us to re-interpret the standard VEM for Stokes as a VEM, where the test and trial discrete velocities are sought in Poisson-like VE spaces. The upside of this fact is that we inherit from Beirão da Veiga et al. (Numer. Math. 138(3), 581–613, 2018) an explicit analysis of best interpolation results in VE spaces, as well as stabilization estimates that are explicit in terms of the degree of accuracy p of the method. We prove exponential convergence of the hp-VEM for Stokes problems with regular right-hand sides. We corroborate the theoretical estimates with numerical tests for both the p- and hp-versions of the method., Advances in Computational Mathematics, 47 (2), ISSN:1019-7168, ISSN:1572-9044
- Published
- 2021
14. The Stokes complex for Virtual Elements in three dimensions
- Author
-
Beirão da Veiga, L, Dassi, F, Vacca, G, Beirão da Veiga, L., Dassi, F., Vacca, G., Beirão da Veiga, L, Dassi, F, Vacca, G, Beirão da Veiga, L., Dassi, F., and Vacca, G.
- Abstract
This paper has two objectives. On one side, we develop and test numerically divergence-free Virtual Elements in three dimensions, for variable "polynomial" order. These are the natural extension of the two-dimensional divergence-free VEM elements, with some modification that allows for a better computational efficiency. We test the element's performance both for the Stokes and (diffusion dominated) Navier-Stokes equation. The second, and perhaps main, motivation is to show that our scheme, also in three dimensions, enjoys an underlying discrete Stokes complex structure. We build a pair of virtual discrete spaces based on general polytopal partitions, the first one being scalar and the second one being vector valued, such that when coupled with our velocity and pressure spaces, yield a discrete Stokes complex.
- Published
- 2020
15. Extension of the nonconforming Trefftz virtual element method to the Helmholtz problem with piecewise constant wave number
- Author
-
Mascotto, L, Pichler, A, Mascotto L., Pichler A., Mascotto, L, Pichler, A, Mascotto L., and Pichler A.
- Abstract
We extend the nonconforming Trefftz virtual element method introduced in [30] to the case of the fluid-fluid interface problem, that is, a Helmholtz problem with piecewise constant wave number. With respect to the original approach, we address two additional issues: firstly, we define the coupling of local approximation spaces with piecewise constant wave numbers; secondly, we enrich such local spaces with special functions capturing the physical behavior of the solution to the target problem. As these two issues are directly related to an increase of the number of degrees of freedom, we use a reduction strategy inspired by [31], which allows to mitigate the growth of the dimension of the approximation space when considering h- and p-refinements. This renders the new method highly competitive in comparison to other Trefftz and quasi-Trefftz technologies tailored for the Helmholtz problem with piecewise constant wave number. A wide range of numerical experiments, including the p-version with quasi-uniform meshes and the hp-version with isotropic and anisotropic mesh refinements, is presented.
- Published
- 2020
16. The p- and hp-versions of the virtual element method for elliptic eigenvalue problems
- Author
-
O., C, Gardini, F, Manzini, G, Mascotto, L, Vacca, G, O. Certik, Gardini F., Manzini G., Mascotto L., Vacca G., O., C, Gardini, F, Manzini, G, Mascotto, L, Vacca, G, O. Certik, Gardini F., Manzini G., Mascotto L., and Vacca G.
- Abstract
We discuss the p- and hp-versions of the virtual element method for the approximation of eigenpairs of elliptic operators with a potential term on polygonal meshes. An application of this model is provided by the Schrödinger equation with a pseudo-potential term. As an interesting byproduct, we present for the first time in literature an explicit construction of the stabilization of the mass matrix. We present in detail the analysis of the p-version of the method, proving exponential convergence in the case of analytic eigenfunctions. The theoretical results are supplied with a wide set of experiments. We also show numerically that, in the case of eigenfunctions with finite Sobolev regularity, an exponential approximation of the eigenvalues in terms of the cubic root of the number of degrees of freedom can be obtained by employing hp-refinements. Importantly, the geometric flexibility of polygonal meshes is exploited in the construction of the hp-spaces.
- Published
- 2020
17. Bricks for the mixed high-order virtual element method: Projectors and differential operators
- Author
-
Dassi, F, Vacca, G, Dassi, F, and Vacca, G
- Abstract
We present the essential tools to deal with virtual element method (VEM) for the approximation of solutions of partial differential equations in mixed form. Functional spaces, degrees of freedom, projectors and differential operators are described emphasizing how to build them in a virtual element framework and for a general approximation order. To achieve this goal, it was necessary to make a deep analysis on polynomial spaces and decompositions. We exploit such “bricks” to construct virtual element approximations of Stokes, Darcy and Navier–Stokes problems and we provide a series of examples to numerically verify the theoretical behaviour of high-order VEM.
- Published
- 2020
18. Vorticity-stabilized virtual elements for the Oseen equation
- Author
-
L. Beirão da Veiga, Franco Dassi, Giuseppe Vacca, Beirão da Veiga, L, Dassi, F, and Vacca, G
- Subjects
Physics ,convection-dominated problem ,Classical mechanics ,Applied Mathematics ,Modeling and Simulation ,Virtual element ,Vorticity ,Oseen equation ,polygonal meshe - Abstract
In this paper, we extend the divergence-free VEM of [L. Beirão da Veiga, C. Lovadina and G. Vacca, Virtual elements for the Navier–Stokes problem on polygonal meshes, SIAM J. Numer. Anal. 56 (2018) 1210–1242] to the Oseen problem, including a suitable stabilization procedure that guarantees robustness in the convection-dominated case without disrupting the divergence-free property. The stabilization is inspired from [N. Ahmed, G. R. Barrenechea, E. Burman, J. Guzman, A. Linke and C. Merdon, A pressure-robust discretization of Oseen’s equation using stabilization in the vorticity equation, SIAM J. Numer. Anal. 59 (2021) 2746–2774] and includes local SUPG-like terms of the vorticity equation, internal jump terms for the velocity gradients, and an additional VEM stabilization. We derive theoretical convergence results that underline the robustness of the scheme in different regimes, including the convection-dominated case. Furthermore, as in the non-stabilized case, the influence of the pressure on the velocity error is moderate, as it appears only through higher-order terms.
- Published
- 2021
19. SUPG-stabilized virtual elements for diffusion-convection problems: A robustness analysis
- Author
-
Carlo Lovadina, Giuseppe Vacca, Franco Dassi, L. Beirão da Veiga, Beirao da Veiga, L, Dassi, F, Lovadina, C, and Vacca, G
- Subjects
Convection ,Discretization ,010103 numerical & computational mathematics ,01 natural sciences ,SUPG stabilization ,Robustness (computer science) ,Convergence (routing) ,FOS: Mathematics ,Applied mathematics ,Mathematics - Numerical Analysis ,Virtual element method ,0101 mathematics ,Mathematics ,Numerical Analysis ,Applied Mathematics ,Multiplicative function ,Diffusion convection ,Numerical Analysis (math.NA) ,Term (time) ,010101 applied mathematics ,Computational Mathematics ,Modeling and Simulation ,Convection dominated problem ,Element (category theory) ,Analysis ,Polygonal meshe - Abstract
The objective of this contribution is to develop a convergence analysis for SUPG-stabilized Virtual Element Methods in diffusion-convection problems that is robust also in the convection dominated regime. For the original method introduced in [Benedetto et al., CMAME 2016] we are able to show an “almost uniform” error bound (in the sense that the unique term that depends in an unfavourable way on the parameters is damped by a higher order mesh-size multiplicative factor). We also introduce a novel discretization of the convection term that allows us to develop error estimates that are fully robust in the convection dominated cases. We finally present some numerical result.
- Published
- 2021
20. A virtual element method for the miscible displacement of incompressible fluids in porous media
- Author
-
Alexander Pichler, L. Beirão da Veiga, Giuseppe Vacca, Beirao da Veiga, L, Pichler, A, and Vacca, G
- Subjects
Discretization ,Computer science ,Computational Mechanics ,Regular solution ,Porous media ,General Physics and Astronomy ,010103 numerical & computational mathematics ,Miscible fluid flow ,01 natural sciences ,Physics::Fluid Dynamics ,Incompressible flow ,FOS: Mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Partial differential equation ,Mechanical Engineering ,Mathematical analysis ,Numerical Analysis (math.NA) ,Computer Science Applications ,010101 applied mathematics ,Nonlinear system ,Virtual element methods ,Flow (mathematics) ,Mechanics of Materials ,Compressibility ,Displacement (fluid) ,Polygonal meshe - Abstract
In the present contribution, we construct a virtual element (VE) discretization for the problem of miscible displacement of one incompressible fluid by another, described by a time-dependent coupled system of nonlinear partial differential equations. Our work represents a first study to investigate the premises of virtual element methods (VEM) for complex fluid flow problems. We combine the VEM discretization with a time stepping scheme and develop a complete theoretical analysis of the method under the assumption of a regular solution. The scheme is then tested both on a regular and on a more realistic test case.
- Published
- 2021
- Full Text
- View/download PDF
21. First-order VEM for Reissner–Mindlin plates
- Author
-
Elio Sacco, Luca Patruno, Antonio Maria D'Altri, S. de Miranda, D'Altri, A. M., Patruno, L., de Miranda, S., Sacco, E., and D’Altri, A. M.
- Subjects
Virtual element method · Shear deformable plates · Locking-free · Polygonal meshes · Reissner–Mindlin plates ,Applied Mathematics ,Mechanical Engineering ,Shear deformable plate ,Mathematical analysis ,Polygonal meshes ,Computational Mechanics ,Boundary (topology) ,Ocean Engineering ,Bending ,Reissner–Mindlin plate ,Displacement (vector) ,Shear (sheet metal) ,Computational Mathematics ,Transverse plane ,Reissner–Mindlin plates ,Computational Theory and Mathematics ,Simple (abstract algebra) ,Degree of a polynomial ,Shear deformable plates ,Virtual element method ,Locking-free ,Mathematics ,Stiffness matrix ,Polygonal meshe - Abstract
In this paper, a first-order virtual element method for Reissner–Mindlin plates is presented. A standard displacement-based variational formulation is employed, assuming transverse displacement and rotations as independent variables. In the framework of the first-order virtual element, a piecewise linear approximation is assumed for both displacement and rotations on the boundary of the element. The consistent term of the stiffness matrix is determined assuming uncoupled polynomial approximations for the generalized strains, with different polynomial degrees for bending and shear parts. In order to mitigate shear locking in the thin-plate limit while keeping the element formulation as simple as possible, a selective scheme for the stabilization term of the stiffness matrix is introduced, to indirectly enrich the approximation of the transverse displacement with respect to that of the rotations. Element performance is tested on various numerical examples involving both thin and thick plates and different polygonal meshes.
- Published
- 2021
22. Equilibrium analysis of an immersed rigid leaflet by the virtual element method
- Author
-
Giuseppe Vacca, Claudio Canuto, L. Beirão da Veiga, Ricardo H. Nochetto, Beirao Da Veiga, L, Canuto, C, Nochetto, R, and Vacca, G
- Subjects
shape calculu ,fluid-structure interaction ,polygonal meshes ,shape calculus ,Virtual elements ,010103 numerical & computational mathematics ,Spring (mathematics) ,01 natural sciences ,Fluid–structure interaction ,FOS: Mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Physics ,Leaflet (botany) ,Applied Mathematics ,Virtual element ,Mechanics ,Numerical Analysis (math.NA) ,polygonal meshe ,010101 applied mathematics ,Modeling and Simulation ,Shape calculus ,Compressibility ,Element (category theory) ,Communication channel - Abstract
We study, both theoretically and numerically, the equilibrium of a hinged rigid leaflet with an attached rotational spring, immersed in a stationary incompressible fluid within a rigid channel. Through a careful investigation of the properties of the domain functional describing the angular momentum exerted by the fluid on the leaflet (which depends on both the leaflet angular position and its thickness), we identify sufficient conditions on the spring stiffness function for the existence (and uniqueness) of equilibrium positions. This study resorts to techniques from shape differential calculus. We propose a numerical technique that exploits the mesh flexibility of the Virtual Element Method (VEM). A (polygonal) computational mesh is generated by cutting a fixed background grid with the leaflet geometry, and the problem is then solved with stable VEM Stokes elements of degrees [Formula: see text] and [Formula: see text] combined with a bisection algorithm. We prove quasi-optimal error estimates and present a large array of numerical experiments to document the accuracy and robustness with respect to degenerate geometry of the proposed methodology.
- Published
- 2020
- Full Text
- View/download PDF
23. The harmonic virtual element method: stabilization and exponential convergence for the Laplace problem on polygonal domains
- Author
-
Alexey Chernov, Lorenzo Mascotto, Chernov, A, and Mascotto, L
- Subjects
General Mathematics ,harmonic polynomial ,Degrees of freedom (statistics) ,Boundary (topology) ,Harmonic (mathematics) ,010103 numerical & computational mathematics ,01 natural sciences ,Trefftz method ,FOS: Mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Mathematics ,Laplace's equation ,Applied Mathematics ,Mathematical analysis ,Order (ring theory) ,Numerical Analysis (math.NA) ,Laplace equation ,polygonal meshe ,Finite element method ,010101 applied mathematics ,Computational Mathematics ,virtual element method ,Rate of convergence ,hp Galerkin method ,Piecewise - Abstract
We introduce the harmonic virtual element method (harmonic VEM), a modification of the virtual element method (VEM) for the approximation of the 2D Laplace equation using polygonal meshes. The main difference between the harmonic VEM and the VEM is that in the former method only boundary degrees of freedom are employed. Such degrees of freedom suffice for the construction of a proper energy projector on (piecewise harmonic) polynomial spaces. The harmonic VEM can also be regarded as an "$H^1$-conformisation" of the Trefftz discontinuous Galerkin-finite element method (TDG-FEM). We address the stabilization of the proposed method and develop an $hp$ version of harmonic VEM for the Laplace equation on polygonal domains. As in Trefftz DG-FEM, the asymptotic convergence rate of harmonic VEM is exponential and reaches order $\mathcal O ( \exp(-b\sqrt[2]{N}))$, where $N$ is the number of degrees of freedom. This result overperformes its counterparts in the framework of $hp$ FEM and $hp$ VEM, where the asymptotic rate of convergence is of order $\mathcal O ( \exp(-b\sqrt[3]{N}) )$., 25 pages, 9 figures
- Published
- 2018
24. An extended B-Rep solid modeling kernel integrating mesh and NURBS faces
- Author
-
Roberto Raffaeli, Giacomo Ferrari, Giulio Casciola, Serena Morigi, Serena Morigi, Giacomo Ferrari, Giulio Casciola, and Roberto Raffaeli
- Subjects
0209 industrial biotechnology ,Computer science ,Computational Mechanics ,020207 software engineering ,02 engineering and technology ,Solid modeling ,polygonal meshes ,polygonal meshe ,B-Rep ,Computer Graphics and Computer-Aided Design ,Computational science ,regularized boolean operations ,Solid modeling kernel ,Computational Mathematics ,020901 industrial engineering & automation ,Kernel (statistics) ,0202 electrical engineering, electronic engineering, information engineering ,Polygon mesh ,ComputingMethodologies_COMPUTERGRAPHICS ,Parametric statistics - Abstract
In several application contexts, virtual solid models require to integrate portions of polygonal meshes with synthetic models, designed by traditional parametric/analytical multipatches systems. The paper reports the research aiming at covering the theoretical and numerical aspects connected with an extended geometric solid modeling system, focusing on the B-Rep models and introducing the new paradigm of Extended B-Rep (EB-Rep), which is able to integrate mesh-faces as part of a B-rep model. This paradigm introduces a notion of continuity between parametric and discrete representations, regularized Boolean Operations, a join operator and an approach to represent a valence semi-regular mesh as an EB-Rep structure. A prototype of the geometric solid modeling kernel has been realized and tested in the OpenCascade library environment.
- Published
- 2018
25. A multigrid algorithm for the p-version of the virtual element method
- Author
-
Paola F. Antonietti, Lorenzo Mascotto, Marco Verani, Antonietti, P, Mascotto, L, and Verani, M
- Subjects
Polynomial ,Discretization ,Uniform convergence ,010103 numerical & computational mathematics ,01 natural sciences ,Multigrid method ,FOS: Mathematics ,Applied mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,multigrid ,Mathematics ,Numerical Analysis ,Sequence ,Preconditioner ,Applied Mathematics ,Numerical analysis ,Polygonal meshes ,Numerical Analysis (math.NA) ,Galerkin methods ,010101 applied mathematics ,Virtual element methods ,Computational Mathematics ,P multigrid ,Modeling and Simulation ,Analysis ,P Galerkin method ,Smoothing ,Polygonal meshe - Abstract
We present a multigrid algorithm for the solution of the linear systems of equations stemming from the $p-$version of the Virtual Element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed decreasing progressively the polynomial approximation degree of the Virtual Element space, as in standard $p$-multigrid schemes. The construction of the interspace operators relies on auxiliary Virtual Element spaces, where it is possible to compute higher order polynomial projectors. We prove that the multigrid scheme is uniformly convergent, provided the number of smoothing steps is chosen sufficiently large. We also demonstrate that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom that can be employed to accelerate the convergence of classical Krylov-based iterative schemes. Numerical experiments validate the theoretical results., 25 pages, 3 figures
- Published
- 2018
26. Veamy: an extensible object-oriented C++ library for the virtual element method
- Author
-
Ortiz-Bernardin, A, Alvarez, C, Hitschfeld-Kahler, N, Russo, A, Silva-Valenzuela, R, Olate-Sanzana, E, Ortiz-Bernardin A., Alvarez C., Hitschfeld-Kahler N., Russo A., Silva-Valenzuela R., Olate-Sanzana E., Ortiz-Bernardin, A, Alvarez, C, Hitschfeld-Kahler, N, Russo, A, Silva-Valenzuela, R, Olate-Sanzana, E, Ortiz-Bernardin A., Alvarez C., Hitschfeld-Kahler N., Russo A., Silva-Valenzuela R., and Olate-Sanzana E.
- Abstract
This paper summarizes the development of Veamy, an object-oriented C++ library for the virtual element method (VEM) on general polygonal meshes, whose modular design is focused on its extensibility. The linear elastostatic and Poisson problems in two dimensions have been chosen as the starting stage for the development of this library. The theory of the VEM, upon which Veamy is built, is presented using a notation and a terminology that resemble the language of the finite element method (FEM) in engineering analysis. Several examples are provided to demonstrate the usage of Veamy, and in particular, one of them features the interaction between Veamy and the polygonal mesh generator PolyMesher. A computational performance comparison between VEM and FEM is also conducted. Veamy is free and open source software.
- Published
- 2019
27. A nonconforming Trefftz virtual element method for the Helmholtz problem
- Author
-
Mascotto, L, Perugia, I, Pichler, A, Mascotto L., Perugia I., Pichler A., Mascotto, L, Perugia, I, Pichler, A, Mascotto L., Perugia I., and Pichler A.
- Abstract
We introduce a novel virtual element method (VEM) for the two-dimensional Helmholtz problem endowed with impedance boundary conditions. Local approximation spaces consist of Trefftz functions, i.e. functions belonging to the kernel of the Helmholtz operator. The global trial and test spaces are not fully discontinuous, but rather interelement continuity is imposed in a nonconforming fashion. Although their functions are only implicitly defined, as typical of the VEM framework, they contain discontinuous subspaces made of functions known in closed form and with good approximation properties (plane-waves, in our case). We carry out an abstract error analysis of the method, and derive h-version error estimates. Moreover, we initiate its numerical investigation by presenting a first test, which demonstrates the theoretical convergence rates.
- Published
- 2019
28. A nonconforming Trefftz virtual element method for the Helmholtz problem: Numerical aspects
- Author
-
Mascotto, L, Perugia, I, Pichler, A, Mascotto L., Perugia I., Pichler A., Mascotto, L, Perugia, I, Pichler, A, Mascotto L., Perugia I., and Pichler A.
- Abstract
We discuss the implementation details and the numerical performance of the recently introduced nonconforming Trefftz virtual element method (Mascotto et al., 2018) for the 2D Helmholtz problem. In particular, we present a strategy to significantly reduce the ill-conditioning of the original method; such a recipe is based on an automatic filtering of the basis functions edge by edge, and therefore allows for a notable reduction of the number of degrees of freedom. A widespread set of numerical experiments, including an application to acoustic scattering, the h-, p-, and hp-versions of the method, is presented. Moreover, a comparison with other Trefftz-based methods for the Helmholtz problem shows that this novel approach results in robust and effective performance.
- Published
- 2019
29. The harmonic virtual element method: Stabilization and exponential convergence for the Laplace problem on polygonal domains
- Author
-
Chernov, A, Mascotto, L, Chernov A., Mascotto L., Chernov, A, Mascotto, L, Chernov A., and Mascotto L.
- Abstract
We introduce the harmonic virtual element method (VEM) (harmonic VEM), a modification of the VEM (Beirão da Veiga et al. (2013) Basic principles of virtual element methods. Math. Models Methods Appl. Sci., 23, 199-214.) for the approximation of the two-dimensional Laplace equation using polygonal meshes. The main difference between the harmonic VEM and the VEM is that in the former method only boundary degrees of freedom are employed. Such degrees of freedom suffice for the construction of a proper energy projector on (piecewise harmonic) polynomial spaces. The harmonic VEM can also be regarded as an '$H^1$-conformisation' of the Trefftz discontinuous Galerkin-finite element method (TDG-FEM) (Hiptmair et al. (2014) Approximation by harmonic polynomials in starshaped domains and exponential convergence of Trefftz hp-DGFEM. ESAIM Math. Model. Numer. Anal., 48, 727-752.). We address the stabilization of the proposed method and develop an hp version of harmonic VEM for the Laplace equation on polygonal domains. As in TDG-FEM, the asymptotic convergence rate of harmonic VEM is exponential and reaches order $mathscr{O}(exp (-bsqrt [2]{N}))$, where $N$ is the number of degrees of freedom. This result overperforms its counterparts in the framework of hp FEM (Schwab, C. (1998)p- and hp-Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics. Clarendon Press Oxford.) and hp VEM (Beirão da Veiga et al. (2018) Exponential convergence of the hp virtual element method with corner singularity. Numer. Math., 138, 581-613.), where the asymptotic rate of convergence is of order $mathscr{O}(exp(-bsqrt [3]{N}))$.
- Published
- 2019
30. Mimetic finite difference methods for Hamiltonian wave equations in 2D
- Author
-
Giuseppe Vacca, Luciano Lopez, L. Beirão da Veiga, Beirao da Veiga, L, Lopez, L, and Vacca, G
- Subjects
Discretization ,010103 numerical & computational mathematics ,01 natural sciences ,Mimetic finite difference method ,Hamiltonian system ,symbols.namesake ,Computational Theory and Mathematic ,FOS: Mathematics ,Covariant Hamiltonian field theory ,Mathematics - Numerical Analysis ,Hamiltonian systems ,0101 mathematics ,Mathematics::Symplectic Geometry ,Mathematics ,Mimetic finite difference methods ,Polygonal meshes ,Mathematical analysis ,Finite difference method ,Finite difference ,Numerical Analysis (math.NA) ,Wave equation ,010101 applied mathematics ,Computational Mathematics ,Computational Theory and Mathematics ,Modeling and Simulation ,symbols ,Hamiltonian (quantum mechanics) ,Polygonal meshe ,Symplectic geometry - Abstract
In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in time to integrate the semi-discrete Hamiltonian system. The main characteristic of MFD methods, when applied to stationary problems, is to mimic important properties of the continuous system. This approach, associated with a symplectic method for the time integration yields a full numerical procedure suitable to integrate Hamiltonian problems. A complete theoretical analysis of the method and some numerical simulations are developed in the paper., Comment: 26 pages, 8 figures
- Published
- 2017
31. The Stokes Complex for Virtual Elements with Application to Navier–Stokes Flows
- Author
-
David Mora, L. Beirão da Veiga, Giuseppe Vacca, Beirao da Veiga, L, Mora, D, and Vacca, G
- Subjects
Numerical Analysis ,Pure mathematics ,Discretization ,Applied Mathematics ,Computability ,General Engineering ,Degrees of freedom (physics and chemistry) ,Structure (category theory) ,010103 numerical & computational mathematics ,Differential operator ,01 natural sciences ,Discrete velocity ,Theoretical Computer Science ,Virtual elements ,010101 applied mathematics ,Computational Mathematics ,Computational Theory and Mathematics ,Navier stokes ,0101 mathematics ,Element (category theory) ,Software ,Mathematics ,Polygonal meshe ,Stokes complex - Abstract
In the present paper, we investigate the underlying Stokes complex structure of the Virtual Element Method for Stokes and Navier–Stokes introduced in previous papers by the same authors, restricting our attention to the two dimensional case. We introduce a Virtual Element space $${\varPhi }_h \subset H^2({\varOmega })$$ Φ h ⊂ H 2 ( Ω ) and prove that the triad $$\{{\varPhi }_h, {\varvec{V}}_h, Q_h\}$$ { Φ h , V h , Q h } (with $${\varvec{V}}_h$$ V h and $$Q_h$$ Q h denoting the discrete velocity and pressure spaces) is an exact Stokes complex. Furthermore, we show the computability of the associated differential operators in terms of the adopted degrees of freedom and explore also a different discretization of the convective trilinear form. The theoretical findings are supported by numerical tests.
- Published
- 2019
32. The Stokes complex for Virtual Elements in three dimensions
- Author
-
Franco Dassi, L. Beirão da Veiga, Giuseppe Vacca, Beirão da Veiga, L, Dassi, F, and Vacca, G
- Subjects
Polynomial ,Computer science ,Applied Mathematics ,Mathematical analysis ,Virtual element ,010103 numerical & computational mathematics ,Extension (predicate logic) ,Numerical Analysis (math.NA) ,polygonal meshe ,01 natural sciences ,Navier-Stokes equation ,010101 applied mathematics ,Modeling and Simulation ,FOS: Mathematics ,Order (group theory) ,Mathematics - Numerical Analysis ,0101 mathematics ,Variable (mathematics) ,Stokes complex - Abstract
This paper has two objectives. On one side, we develop and test numerically divergence-free Virtual Elements in three dimensions, for variable “polynomial” order. These are the natural extension of the two-dimensional divergence-free VEM elements, with some modification that allows for a better computational efficiency. We test the element’s performance both for the Stokes and (diffusion dominated) Navier–Stokes equation. The second, and perhaps main, motivation is to show that our scheme, also in three dimensions, enjoys an underlying discrete Stokes complex structure. We build a pair of virtual discrete spaces based on general polytopal partitions, the first one being scalar and the second one being vector valued, such that when coupled with our velocity and pressure spaces, yield a discrete Stokes complex.
- Published
- 2019
- Full Text
- View/download PDF
33. The $p$- and $hp$-versions of the virtual element method for elliptic eigenvalue problems
- Author
-
O. Certik, Gardini F., Manzini G., Mascotto L., Vacca G., O., C, Gardini, F, Manzini, G, Mascotto, L, and Vacca, G
- Subjects
Eigenvalue problem ,p- and hp-Galerkin method ,FOS: Mathematics ,65N12, 65N30 ,Virtual element method ,Numerical Analysis (math.NA) ,Mathematics - Numerical Analysis ,Polygonal meshe - Abstract
We discuss the $p$- and the $hp$-versions of the virtual element method for the approximation of eigenpairs of elliptic operators with a potential term on polygonal meshes. An application of this model is provided by the Schr\"odinger equation with a pseudo-potential term. We present in details the analysis of the p-version of the method, proving exponential convergence in the case of analytic eigenfunctions. The theoretical results are supplied with a wide set of experiments. We also show numerically that, in the case of eigenfunctions with finite Sobolev regularity, an exponential approximation of the eigenvalues in terms of the cubic root of the number of degrees of freedom can be obtained by employing $hp$-refinements. Importantly, the geometric flexibility of polygonal meshes is exploited in the construction of the $hp$-spaces., Comment: 25 pages, 7 figures
- Published
- 2018
34. Virtual elements for the navier-stokes problem on polygonal meshes
- Author
-
Beirao da Veiga, L, Lovadina, C, Vacca, G, VACCA, GIUSEPPE, Beirao da Veiga, L, Lovadina, C, Vacca, G, and VACCA, GIUSEPPE
- Abstract
A family of virtual element methods for the two-dimensional Navier-Stokes equations is proposed and analyzed. The schemes provide a discrete velocity field which is pointwise divergence-free. A rigorous error analysis is developed, showing that the methods are stable and optimally convergent. Several numerical tests are presented, confirming the theoretical predictions. A comparison with some mixed finite elements is also performed.
- Published
- 2018
35. A multigrid algorithm for the p -version of the virtual element method
- Author
-
Antonietti, P, Mascotto, L, Verani, M, Antonietti P. F., Mascotto L., Verani M., Antonietti, P, Mascotto, L, Verani, M, Antonietti P. F., Mascotto L., and Verani M.
- Abstract
We present a multigrid algorithm for the solution of the linear systems of equations stemming from the p-version of the virtual element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed decreasing progressively the polynomial approximation degree of the virtual element space, as in standard p-multigrid schemes. The construction of the interspace operators relies on auxiliary virtual element spaces, where it is possible to compute higher order polynomial projectors. We prove that the multigrid scheme is uniformly convergent, provided the number of smoothing steps is chosen sufficiently large. We also demonstrate that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom that can be employed to accelerate the convergence of classical Krylov-based iterative schemes. Numerical experiments validate the theoretical results.
- Published
- 2018
36. An H1-conforming virtual element for Darcy and Brinkman equations
- Author
-
Vacca, G, Vacca, G., Vacca, G, and Vacca, G.
- Abstract
The focus of this paper is on developing a virtual element method (VEM) for Darcy and Brinkman equations. In [L. Beirão da Veiga, C. Lovadina and G. Vacca, ESAIM Math. Model. Numer. Anal. 51 (2017)], we presented a family of virtual elements for Stokes equations and we defined a new virtual element space of velocities such that the associated discrete kernel is pointwise divergence-free. We use a slightly different virtual element space having two fundamental properties: The L2-projection onto Pk is exactly computable on the basis of the degrees of freedom, and the associated discrete kernel is still pointwise divergence-free. The resulting numerical scheme for the Darcy equation has optimal order of convergence and H1-conforming velocity solution. We can apply the same approach to develop a robust virtual element method for the Brinkman equation that is stable for both the Stokes and Darcy limit case. We provide a rigorous error analysis of the method and several numerical tests.
- Published
- 2018
37. Bricks for the mixed high-order virtual element method: projectors and differential operators
- Author
-
Franco Dassi, Giuseppe Vacca, Dassi, F, and Vacca, G
- Subjects
Numerical Analysis ,Polynomial ,Partial differential equation ,Series (mathematics) ,Applied Mathematics ,Numerical analysis ,Degrees of freedom (statistics) ,Computational mathematics ,010103 numerical & computational mathematics ,Differential operator ,01 natural sciences ,Mixed problem ,010101 applied mathematics ,Algebra ,Computational Mathematics ,Computational Mathematic ,Virtual element method ,Mathematics - Numerical Analysis ,0101 mathematics ,Element (category theory) ,Projector ,High-order ,Numerical Analysi ,Polygonal meshe ,Mathematics - Abstract
We present the essential tools to deal with virtual element method (VEM) for the approximation of solutions of partial differential equations in mixed form. Functional spaces, degrees of freedom, projectors and differential operators are described emphasizing how to build them in a virtual element framework and for a general approximation order. To achieve this goal, it was necessary to make a deep analysis on polynomial spaces and decompositions. We exploit such “bricks” to construct virtual element approximations of Stokes, Darcy and Navier–Stokes problems and we provide a series of examples to numerically verify the theoretical behaviour of high-order VEM.
- Published
- 2018
38. Virtual Elements for the Navier--Stokes Problem on Polygonal Meshes
- Author
-
Carlo Lovadina, Giuseppe Vacca, L. Beirão da Veiga, Beirao da Veiga, L, Lovadina, C, and Vacca, G
- Subjects
Field (physics) ,Mathematics::Analysis of PDEs ,010103 numerical & computational mathematics ,01 natural sciences ,Navier-Stokes equation ,Physics::Fluid Dynamics ,FOS: Mathematics ,Polygon mesh ,Navier stokes ,Virtual element method ,Mathematics - Numerical Analysis ,0101 mathematics ,Navier–Stokes equations ,Numerical Analysi ,Mathematics ,Pointwise ,Numerical Analysis ,Numerical analysis ,Applied Mathematics ,Mathematical analysis ,Computational mathematics ,Numerical Analysis (math.NA) ,010101 applied mathematics ,Computational Mathematics ,Computational Mathematic ,Element (category theory) ,Polygonal meshe - Abstract
A family of virtual element methods for the two-dimensional Navier-Stokes equations is proposed and analyzed. The schemes provide a discrete velocity field which is pointwise divergence-free. A rigorous error analysis is developed, showing that the methods are stable and optimally convergent. Several numerical tests are presented, confirming the theoretical predictions. A comparison with some mixed finite elements is also performed.
- Published
- 2018
- Full Text
- View/download PDF
39. An H1-conforming virtual element for Darcy and Brinkman equations
- Author
-
Giuseppe Vacca and Vacca, G
- Subjects
Pointwise ,Basis (linear algebra) ,Applied Mathematics ,Mathematical analysis ,Degrees of freedom (physics and chemistry) ,010103 numerical & computational mathematics ,polygonal meshe ,01 natural sciences ,Darcy–Weisbach equation ,Brinkman equation ,010101 applied mathematics ,Kernel (image processing) ,Rate of convergence ,Modeling and Simulation ,Darcy equation ,Limit (mathematics) ,Virtual element method ,0101 mathematics ,Element (category theory) ,Brinkman equations ,Mathematics - Abstract
The focus of this paper is on developing a virtual element method (VEM) for Darcy and Brinkman equations. In [L. Beirão da Veiga, C. Lovadina and G. Vacca, ESAIM Math. Model. Numer. Anal. 51 (2017)], we presented a family of virtual elements for Stokes equations and we defined a new virtual element space of velocities such that the associated discrete kernel is pointwise divergence-free. We use a slightly different virtual element space having two fundamental properties: the [Formula: see text]-projection onto [Formula: see text] is exactly computable on the basis of the degrees of freedom, and the associated discrete kernel is still pointwise divergence-free. The resulting numerical scheme for the Darcy equation has optimal order of convergence and [Formula: see text]-conforming velocity solution. We can apply the same approach to develop a robust virtual element method for the Brinkman equation that is stable for both the Stokes and Darcy limit case. We provide a rigorous error analysis of the method and several numerical tests.
- Published
- 2018
40. A nonconforming Trefftz virtual element method for the Helmholtz problem
- Author
-
Alexander Pichler, Lorenzo Mascotto, Ilaria Perugia, Mascotto, L, Perugia, I, and Pichler, A
- Subjects
Physics ,nonconforming method ,Applied Mathematics ,Mathematical analysis ,Plane wave ,35J05, 65N12, 65N15, 65N30 ,010103 numerical & computational mathematics ,Numerical Analysis (math.NA) ,Impedance boundary condition ,polygonal meshe ,01 natural sciences ,010101 applied mathematics ,Helmholtz problem ,plane-wave ,Trefftz method ,Modeling and Simulation ,FOS: Mathematics ,Mathematics - Numerical Analysis ,Virtual element method ,0101 mathematics ,Element (category theory) - Abstract
We introduce a novel virtual element method (VEM) for the two dimensional Helmholtz problem endowed with impedance boundary conditions. Local approximation spaces consist of Trefftz functions, i.e., functions belonging to the kernel of the Helmholtz operator. The global trial and test spaces are not fully discontinuous, but rather interelement continuity is imposed in a nonconforming fashion. Although their functions are only implicitly defined, as typical of the VEM framework, they contain discontinuous subspaces made of functions known in closed form and with good approximation properties (plane waves, in our case). We carry out an abstract error analysis of the method, and derive $h$-version error estimates. Moreover, we initiate its numerical investigation by presenting a first test, which demonstrates the theoretical convergence rates., Comment: 27 pages, 5 figures
- Published
- 2018
- Full Text
- View/download PDF
41. Mimetic finite difference methods for Hamiltonian wave equations in 2D
- Author
-
Beirao da Veiga, L, Lopez, L, Vacca, G, Beirao da Veiga, L., Lopez, L., Vacca, G., Beirao da Veiga, L, Lopez, L, Vacca, G, Beirao da Veiga, L., Lopez, L., and Vacca, G.
- Abstract
In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimensions. We construct a two step procedure in which we first discretize the space by the Mimetic Finite Difference (MFD) method and then we employ a standard symplectic scheme to integrate the semi-discrete Hamiltonian system derived. The main characteristic of the MFD methods, when applied to stationary problems, is to mimic important properties of the continuous system. This approach yields a full numerical procedure suitable to integrate Hamiltonian problems. A complete theoretical analysis of the method and some numerical simulations are developed in the paper.
- Published
- 2017
42. Arbitrary order 2D virtual elements for polygonal meshes: part I, elastic problem
- Author
-
Artioli, E, BEIRAO DA VEIGA, L, Lovadina, C, Sacco, E, Artioli, E., BEIRAO DA VEIGA, LOURENCO, Lovadina, C., Sacco, E., Artioli, E, BEIRAO DA VEIGA, L, Lovadina, C, Sacco, E, Artioli, E., BEIRAO DA VEIGA, LOURENCO, Lovadina, C., and Sacco, E.
- Abstract
The present work deals with the formulation of a virtual element method for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II (Artioli et al. in Comput Mech, 2017) the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a computer code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.
- Published
- 2017
43. Virtual Element Methods for hyperbolic problems on polygonal meshes
- Author
-
Vacca, G, VACCA, GIUSEPPE, Vacca, G, and VACCA, GIUSEPPE
- Abstract
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 H1 semi-norm and L2 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 set of numerical tests
- Published
- 2017
44. Virtual Elements for a shear-deflection formulation of Reissner-Mindlin plates
- Author
-
Gonzalo Rivera, L. da Veiga, David Mora, Da Veiga, L, Mora, D, and Rivera, G
- Subjects
Error analysi ,Algebra and Number Theory ,business.industry ,Applied Mathematics ,65N30, 65N12, 74K20, 74S05, 65N15 (Primary) ,Computational mathematics ,010103 numerical & computational mathematics ,Structural engineering ,Numerical Analysis (math.NA) ,01 natural sciences ,010101 applied mathematics ,Computational Mathematics ,Deflection (engineering) ,Error analysis ,Computational Mathematic ,FOS: Mathematics ,Virtual element method ,Mathematics - Numerical Analysis ,Reissner-Mindlin plate ,0101 mathematics ,business ,Polygonal meshe ,Mathematics - Abstract
We present a virtual element method for the Reissner-Mindlin plate bending problem which uses shear strain and deflection as discrete variables without the need of any reduction operator. The proposed method is conforming in $[H^{1}(\Omega)]^2 \times H^2(\Omega)$ and has the advantages of using general polygonal meshes and yielding a direct approximation of the shear strains. The rotations are then obtained by a simple postprocess from the shear strain and deflection. We prove convergence estimates with involved constants that are uniform in the thickness $t$ of the plate. Finally, we report numerical experiments which allow us to assess the performance of the method.
- Published
- 2017
45. Arbitrary order 2D virtual elements for polygonal meshes: part I, elastic problem
- Author
-
Artioli, Edoardo, da Veiga, Lourenco Beirao, Lovadina, Carlo, Sacco, Elio, Artioli, E., Beirão da Veiga, L., Lovadina, C., Sacco, E., Artioli, E, BEIRAO DA VEIGA, L, Lovadina, C, and Sacco, E
- Subjects
Virtual element method ,Elasticity ,Static analysis ,Mechanical Engineering ,Applied Mathematics ,Polygonal meshes ,Computational Mechanics ,Ocean Engineering ,Numerical Analysis (math.NA) ,010103 numerical & computational mathematics ,01 natural sciences ,Static analysi ,010101 applied mathematics ,Computational Mathematics ,Computational Theory and Mathematics ,Computational Theory and Mathematic ,Computational Mathematic ,FOS: Mathematics ,Settore ICAR/08 - Scienza delle Costruzioni ,Mathematics - Numerical Analysis ,0101 mathematics ,Polygonal meshe - Abstract
The present work deals with the formulation of a virtual element method for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II (Artioli et al. in Comput Mech, 2017) the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a computer code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.
- Published
- 2017
46. Virtual Element Methods for hyperbolic problems on polygonal meshes
- Author
-
Giuseppe Vacca and Vacca, G
- Subjects
Mathematical optimization ,Wave propagation ,Virtual Element Method ,Stability (learning theory) ,010103 numerical & computational mathematics ,Volume mesh ,Numerical Analysis (math.NA) ,Wave equation ,01 natural sciences ,010101 applied mathematics ,Hyperbolic problem ,Computational Mathematics ,Computational Theory and Mathematics ,Modeling and Simulation ,Convergence (routing) ,FOS: Mathematics ,Applied mathematics ,Newmark-beta method ,Polygon mesh ,Mathematics - Numerical Analysis ,0101 mathematics ,Element (category theory) ,Mathematics ,Polygonal meshe - Abstract
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.
- Published
- 2016
47. A Virtual Element Method for elastic and inelastic problems on polytope meshes
- Author
-
Carlo Lovadina, David Mora, L. Beirão da Veiga, BEIRAO DA VEIGA, L, Lovadina, C, and Mora, D
- Subjects
Virtual Element Method, Nonlinear structural mechanics ,Large deformation ,Computer science ,Mechanical Engineering ,Virtual Element Method ,Constitutive equation ,Mathematical analysis ,Displacement gradient ,Computational Mechanics ,General Physics and Astronomy ,Polytope ,Numerical Analysis (math.NA) ,Elasticity ,Computer Science Applications ,Nonlinear system ,Mechanics of Materials ,Displacement field ,FOS: Mathematics ,Polygon mesh ,Mathematics - Numerical Analysis ,Elasticity (economics) ,Convergence analysi ,Polygonal meshe - Abstract
We present a Virtual Element Method (VEM) for possibly nonlinear elastic and inelastic problems, mainly focusing on a small deformation regime. The numerical scheme is based on a low-order approximation of the displacement field, as well as a suitable treatment of the displacement gradient. The proposed method allows for general polygonal and polyhedral meshes, it is efficient in terms of number of applications of the constitutive law, and it can make use of any standard black-box constitutive law algorithm. Some theoretical results have been developed for the elastic case. Several numerical results within the 2D setting are presented, and a brief discussion on the extension to large deformation problems is included.
- Published
- 2015
48. Virtual element methods for parabolic problems on polygonal meshes
- Author
-
VACCA, GIUSEPPE, BEIRAO DA VEIGA, LOURENCO, Vacca, G, and BEIRAO DA VEIGA, L
- Subjects
virtual elements ,parabolic problem ,polygonal meshe ,virtual element - Abstract
The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time-dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests.
- Published
- 2015
49. Virtual element methods for parabolic problems on polygonal meshes
- Author
-
Vacca, G, BEIRAO DA VEIGA, L, VACCA, GIUSEPPE, BEIRAO DA VEIGA, LOURENCO, Vacca, G, BEIRAO DA VEIGA, L, VACCA, GIUSEPPE, and BEIRAO DA VEIGA, LOURENCO
- Abstract
The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time-dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests.
- Published
- 2015
50. A Virtual Element Method for elastic and inelastic problems on polytope meshes
- Author
-
BEIRAO DA VEIGA, L, Lovadina, C, Mora, D, BEIRAO DA VEIGA, LOURENCO, Mora, D., BEIRAO DA VEIGA, L, Lovadina, C, Mora, D, BEIRAO DA VEIGA, LOURENCO, and Mora, D.
- Abstract
We present a Virtual Element Method (VEM) for possibly nonlinear elastic and inelastic problems, mainly focusing on a small deformation regime. The numerical scheme is based on a low-order approximation of the displacement field, as well as a suitable treatment of the displacement gradient. The proposed method allows for general polygonal and polyhedral meshes, it is efficient in terms of number of applications of the constitutive law, and it can make use of any standard black-box constitutive law algorithm. Some theoretical results have been developed for the elastic case. Several numerical results within the 2D setting are presented, and a brief discussion on the extension to large deformation problems is included.
- Published
- 2015
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.