Research database

PYGEOM - Polyhedral Galerkin methods for engineering applications to improve disaster risk forecast and management: stabilization-free operator-preserving methods and optimal stabilization methods

Duration:
24 months (2025)
Principal investigator(s):
Project type:
Nationally funded research - PRIN
Funding body:
MINISTERO (MUR)
Project identification number:
P2022BH5CB
PoliTo role:
Coordinator

Abstract

The proposal concerns the development of numerical methods suitable for simulations for safety assessments in several fields of civil and environmental engineering. Civil structures are subject to aging and the evolution of the surrounding environment increases the possibility of dramatic events. Moreover, most natural catastrophic events are hardly predictable and local or global human activities can strongly impact the stability of the environment surrounding civil structures, for example by changing the hydrogeological properties of the soil. This project aims to develop new numerical tools suitable to tackle all the engineering models and the physical phenomena needed to describe such catastrophic events overcoming usual difficulties such as, for example, mesh generation problems, mesh adaptation problems, mesh quality problems, and incompressibility conditions. The Virtual Element Method (VEM) is a polygonal/polyhedral Finite Element Method that was developed for this aim. Although its great success in managing many of the overmentioned problems, some issues still affect its applicability in fields like non-linear elasticity, elastoplasticity, viscoplasticity and in many engineering fields that deal with strongly non-linear problems, strongly anisotropic problems or problems governed by uncommon differential operators. The main issue is related to the non-consistent stabilization term needed by these methods to get well-posed discrete problems for the simulation. This project will focus on stabilization-free VEMs with the aim of retaining the key features of VEM but without the side-effect of the pollution introduced by the stabilizing terms. Stabilization-free VEMs can be developed for general discretizations, but for some specific engineering applications, “ad hoc” methods can be more effective and will be developed. Stabilization-free methods are preferable for low-order methods, whereas for high-order methods can be less effective. For high-order methods, optimal stabilizations requiring a low computational cost will be investigated in the project. Stabilization-free VEMs will be derived and applied to elasticity, elastoplasticity and Richards equations for modelling unsaturated soils. Moreover, optimal cheap stabilization methods, mainly for high-order VEMs and optimal preconditioner and solvers will be investigated from a mathematical viewpoint. These mathematical tools will be applied to understand slope and soil stability and to forecast natural catastrophic events like landslides and to water resources conservation and preservation simulating the roots-soil interaction in unsaturated and fractured porous media. The developed methods will be useful for simulations in the field of geothermal energy exploitation and storage that can help reducing the human impact on the earth's environment.

Structures

Partners

  • C.N.R. - CONSIGLIO NAZIONALE DELLE RICERCHE
  • POLITECNICO DI MILANO
  • POLITECNICO DI TORINO - Coordinator
  • UNIVERSITA' DI MILANO BICOCCA

Keywords

ERC sectors

PE1_17 - Numerical analysis
PE1_21 - Application of mathematics in industry and society
PE8_3 - Civil engineering, architecture, maritime/hydraulic engineering, geotechnics, waste treatment

Budget

Total cost: € 230,449.00
Total contribution: € 230,449.00
PoliTo total cost: € 91,501.00
PoliTo contribution: € 91,501.00