[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"doc-detail-84765-en":3,"doc-seo-84765-105":29,"detail-sidebar-cat-0-en-105":91},{"code":4,"msg":5,"data":6},0,"success",{"doc_id":7,"user_id":8,"nickname":9,"user_avatar":10,"doc_module":4,"category_id":11,"category_name":12,"doc_title":13,"doc_description":14,"doc_content":15,"file_id":16,"file_url":17,"file_type":18,"file_size":19,"view_count":20,"is_deleted":4,"is_public":20,"is_downloadable":20,"audit_status":20,"page_count":21,"language":22,"language_code":23,"site_id":24,"html_lang":23,"table_of_contents":25,"faqs":26,"seo_title":13,"seo_description":14,"update_tm":27,"read_time":28},84765,4398048950312,"Violet","https://ap-avatar.wpscdn.com/avatar/400002538284de19e3c?_k=1778320343897328908",8,"Research & Report","Hydro-mechanical Model for Slope Stability Assessment A Polygonal Stabilization-free Discretization","Rainfall-induced landslides result from the coupled action of subsurface water flow and soil mechanics, requiring reliable numerical simulation of variably saturated porous media. The study formulates a semi-coupled hydro-mechanical model using Richards’ equation and linear elasticity, then discretizes it with a stabilization-free Virtual Element Method on general polygonal meshes. The framework avoids problem-dependent stabilization terms, uses mass lumping for the storage term, and applies Nitsche’s method to weakly impose infiltration and seepage-face boundary conditions with automatic Neumann/Dirichlet switching. Backward Euler time stepping and Picard iterations handle nonlinearities, and experiments confirm stability for rainfall infiltration and slope stability via the Local Factor of Safety.","arXiv :2607 .04778v1 [math .NA] 6 Jul 2026  \nHYDRO-MECHANICAL MODEL FOR SLOPE STABILITY ASSESSMENT: A POLYGONAL STABILIZATION-FREE  \nDISCRETIZATION  \nA PREPRINT  \n Stefano Berrone  \nDipartimento di Scienze Matematiche“G. L. Lagrange”  \nPolitecnico di Torino, Italy stefano .berrone@polito .it  \n Francesca Marcon  \nDipartimento di Scienze Matematiche“G. L. Lagrange”  \nPolitecnico di Torino, Italy francesca .marcon@polito .it  \n Gioana Teora  \nDipartimento di Scienze Matematiche  \n“G. L. Lagrange”  \nPolitecnico di Torino, Italy  \ngioana.teora@polito.it  \nJuly 7, 2026  \nABSTRACT  \nRainfall-induced landslides are governed by the interaction between subsurface water flow and soil mechanics, requiring robust numerical methods for the simulation of variably saturated porous media. In this work, we consider a semi-coupled hydro-mechanical model based on Richards’equation and linear elasticity and propose a numerical framework based on a stabilization-free Virtual Element Method for its spatial discretization. The proposed approach naturally accommodates general polygonal meshes while avoiding problem-dependent stabilization terms, whose design may become challenging when heterogeneous and strongly non-linear coefficients are involved. The approach is combined with a mass-lumping technique to improve stability in the treatment of the storage term and with Nitsche’s method to weakly impose seepage-face and infiltration boundary conditions, allowing for the automatic switching between Neumann and Dirichlet conditions. Time integration is performed using the backward Euler scheme, while non-linearities are handled through a Picard iteration. Numerical experiments demonstrate the stability and robustness of the proposed methodology and show its effectiveness in simulating rainfall infiltration and evaluating slope stability through the Local Factor of Safety.  \nKeywords Polygonal mesh; stabilization-free; infiltration; seepage; soil-stability  \n1 Introduction  \nLandslides are geological processes involving the downslope movement of soil and rock masses, often triggered by intense or prolonged rainfall events that increase soil saturation and reduce slope stability [1] . Rainwater infiltration modifies both the unit weight of the soil and the pore water pressure, thereby altering the stress distribution within hill-slopes and reducing the available shear strength [2] . In particular, an increase in pore water pressure decreases matric suction and effective stress, leading to a reduction in soil cohesion and, consequently, to a higher susceptibility to slope failure [3] .  \nTo investigate these processes, hydro-mechanical multi-physics models are commonly employed, coupling subsurface flow and soil mechanics [4, 3] . Since soils are generally unsaturated or partially saturated, water flow is modelled  \nusing Richards’ equation, which combines Darcy’s law, mass conservation, and constitutive relationships linking saturation and permeability to the pressure-head [1, 3] . The mechanical behaviour is described by the linear momentum equilibrium equation together with a linear elasticity constitutive relation. The resulting mathematical model describes the quasi-static consolidation of variably saturated porous media, commonly referred to as unsaturated poro-elasticity [5] .  \nSlope stability is commonly assessed through the Local Factor of Safety (LFS), a Coulomb stress-field-based indicator defined as the ratio between the Coulomb stress for the potential failure state and the Coulomb stress for the current state of stress under the Mohr-Coulomb criterion [4] . The LFS is evaluated pointwise throughout the computational domain and depends on the effective stress tensor, soil cohesion, and internal friction angle. Unlike traditional limit-equilibrium approaches, the LFS does not require prior assumptions regarding the geometry or location of the failure surface. Consequently, it can be naturally computed on unstructured meshes, providing detailed","cbCaitO9EXDra5RS","https://ap.wps.com/l/cbCaitO9EXDra5RS","pdf",2432728,1,26,"English","en",105,"# Introduction\n## Hydro-mechanical processes in rainfall-triggered landslides\n## Slope stability via Local Factor of Safety (LFS)\n## Numerical methods: FEM and its limitations\n## Virtual Element Method (VEM) and stabilization challenges\n## Stabilization-free VEM framework for the proposed model","[{\"question\":\"What physical processes does the proposed model capture for rainfall-triggered landslides?\",\"answer\":\"It couples subsurface water flow with soil mechanics, modeling variably saturated porous media under rainfall infiltration using Richards’ equation and linear elasticity.\"},{\"question\":\"Why does the method emphasize a stabilization-free Virtual Element Method (SFVEM)?\",\"answer\":\"Standard VEM formulations often need problem-dependent stabilization to regain coercivity, which becomes difficult for strongly non-linear problems such as Richards’ equation; SFVEM avoids these stabilization terms.\"},{\"question\":\"How are infiltration and seepage-face boundary conditions handled?\",\"answer\":\"Nitsche’s method is used to weakly impose infiltration and seepage-face boundary conditions, enabling automatic switching between Neumann and Dirichlet conditions.\"}]",1784198108,66,{"code":4,"msg":30,"data":31},"ok",{"site_id":24,"language":23,"slug":32,"title":13,"keywords":33,"description":14,"schema_data":34,"social_meta":86,"head_meta":88,"extra_data":90,"updated_unix":27},"hydro-mechanical-model-for-slope-stability-assessment-a-polygonal-stabilization-free-discretization","",{"@graph":35,"@context":85},[36,53,68],{"@type":37,"itemListElement":38},"BreadcrumbList",[39,43,47,50],{"item":40,"name":41,"@type":42,"position":20},"https://docshare.wps.com","Home","ListItem",{"item":44,"name":45,"@type":42,"position":46},"https://docshare.wps.com/document/","Document",2,{"item":48,"name":12,"@type":42,"position":49},"https://docshare.wps.com/document/research-report/",3,{"item":51,"name":13,"@type":42,"position":52},"https://docshare.wps.com/document/hydro-mechanical-model-for-slope-stability-assessment-a-polygonal-stabilization-free-discretization/84765/",4,{"url":51,"name":13,"@type":54,"author":55,"headline":13,"publisher":57,"fileFormat":60,"inLanguage":23,"description":14,"dateModified":61,"datePublished":62,"encodingFormat":60,"isAccessibleForFree":63,"interactionStatistic":64},"DigitalDocument",{"name":9,"@type":56},"Person",{"url":40,"name":58,"@type":59},"DocShare","Organization","application/pdf","2026-07-17","2026-07-16",true,{"@type":65,"interactionType":66,"userInteractionCount":20},"InteractionCounter",{"@type":67},"ViewAction",{"@type":69,"mainEntity":70},"FAQPage",[71,77,81],{"name":72,"@type":73,"acceptedAnswer":74},"What physical processes does the proposed model capture for rainfall-triggered landslides?","Question",{"text":75,"@type":76},"It couples subsurface water flow with soil mechanics, modeling variably saturated porous media under rainfall infiltration using Richards’ equation and linear elasticity.","Answer",{"name":78,"@type":73,"acceptedAnswer":79},"Why does the method emphasize a stabilization-free Virtual Element Method (SFVEM)?",{"text":80,"@type":76},"Standard VEM formulations often need problem-dependent stabilization to regain coercivity, which becomes difficult for strongly non-linear problems such as Richards’ equation; SFVEM avoids these stabilization terms.",{"name":82,"@type":73,"acceptedAnswer":83},"How are infiltration and seepage-face boundary conditions handled?",{"text":84,"@type":76},"Nitsche’s method is used to weakly impose infiltration and seepage-face boundary conditions, enabling automatic switching between Neumann and Dirichlet conditions.","https://schema.org",{"og:url":51,"og:type":87,"og:title":13,"og:site_name":58,"og:description":14},"article",{"robots":89,"canonical":51},"index,follow",{"doc_id":7,"site_id":24},{"code":4,"msg":5,"data":92},[93,97,101,105,110,115,120,123,128,131,135],{"id":20,"doc_module":4,"doc_module_name":45,"category_name":94,"show_sort_weight":95,"slug":96},"Story & Novel",90,"story-novel",{"id":46,"doc_module":4,"doc_module_name":45,"category_name":98,"show_sort_weight":99,"slug":100},"Literature",80,"literature",{"id":52,"doc_module":4,"doc_module_name":45,"category_name":102,"show_sort_weight":103,"slug":104},"Exam",70,"exam",{"id":106,"doc_module":4,"doc_module_name":45,"category_name":107,"show_sort_weight":108,"slug":109},5,"Comic",60,"comic",{"id":111,"doc_module":4,"doc_module_name":45,"category_name":112,"show_sort_weight":113,"slug":114},6,"Technology",50,"technology",{"id":116,"doc_module":4,"doc_module_name":45,"category_name":117,"show_sort_weight":118,"slug":119},7,"Healthcare",40,"healthcare",{"id":11,"doc_module":4,"doc_module_name":45,"category_name":12,"show_sort_weight":121,"slug":122},30,"research-report",{"id":124,"doc_module":4,"doc_module_name":45,"category_name":125,"show_sort_weight":126,"slug":127},9,"Religion & Spirituality",20,"religion-spirituality",{"id":126,"doc_module":4,"doc_module_name":45,"category_name":129,"show_sort_weight":126,"slug":130},"World Cup","world-cup",{"id":132,"doc_module":4,"doc_module_name":45,"category_name":133,"show_sort_weight":132,"slug":134},10,"Lifestyle","lifestyle",{"id":136,"doc_module":4,"doc_module_name":45,"category_name":137,"show_sort_weight":106,"slug":138},19,"General","general"]