Pub Date : 2026-02-07DOI: 10.1016/j.cma.2026.118789
Krishna Prasath Logakannan, Mauricio Aristizabal, Geoffrey Bomarito, Zhitong Xu, Shandian Zhe, Robert M. Kirby, Harry Millwater, Jacob Hochhalter
{"title":"Efficiently training SciML models with derivative-informed training data using order truncated imaginary numbers","authors":"Krishna Prasath Logakannan, Mauricio Aristizabal, Geoffrey Bomarito, Zhitong Xu, Shandian Zhe, Robert M. Kirby, Harry Millwater, Jacob Hochhalter","doi":"10.1016/j.cma.2026.118789","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118789","url":null,"abstract":"","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"9 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146138868","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-06DOI: 10.1016/j.cma.2026.118794
Yuguo Fu, Melvan Tan Kian Hao, Dong Chen, Lailai Zhu, Poh Seng Lee
{"title":"Reduced-order initialization model for accelerating topology optimization in thermal-fluid system","authors":"Yuguo Fu, Melvan Tan Kian Hao, Dong Chen, Lailai Zhu, Poh Seng Lee","doi":"10.1016/j.cma.2026.118794","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118794","url":null,"abstract":"","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"77 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146134449","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-06DOI: 10.1016/j.cma.2026.118791
Zhengfei Ren, Junxiang Yang
{"title":"Energy-stable decoupled numerical approximation with practical correction technique for the binary phase field Darcy fluid system","authors":"Zhengfei Ren, Junxiang Yang","doi":"10.1016/j.cma.2026.118791","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118791","url":null,"abstract":"","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"126 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146134451","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-04DOI: 10.1016/j.cma.2026.118787
Juan C. Àlvarez-Hostos, Fernando R. Urgorri, Javier Principe
{"title":"Topology optimisation-based design of duct cross-sections for fully developed magnetohydrodynamic flows","authors":"Juan C. Àlvarez-Hostos, Fernando R. Urgorri, Javier Principe","doi":"10.1016/j.cma.2026.118787","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118787","url":null,"abstract":"","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"90 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146134459","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-02DOI: 10.1016/j.cma.2026.118788
Martin Hanek, Jan Papež, Jakub Šístek
{"title":"Speeding up an unsteady flow simulation by adaptive BDDC and Krylov subspace recycling","authors":"Martin Hanek, Jan Papež, Jakub Šístek","doi":"10.1016/j.cma.2026.118788","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118788","url":null,"abstract":"","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"1 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146110150","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-02DOI: 10.1016/j.cma.2026.118754
Andreas Neofytou, Lucy T. Zhang, H. Alicia Kim
This work presents level set topology optimization method (LSTO) for fluid-structure interaction (FSI) problems. For the solution of the FSI problem the modified immersed finite element method (mIFEM) is used. The proposed formulation provides several advantages compared to the existing methods that rely on unified formulations or remeshing approaches. First it results in separate solid and fluid domains, thus allowing any discretization to be used for each physics. This allows for well established independent fluid and solid solvers to be utilized. Further, this modularity is possible without requiring re-meshing, which maintains efficiency especially for the fluid solver. The finite element method (FEM) is used to solve the flow equations efficiently on an Eulerian grid. The Lagrangian solid is analyzed with the reproducing kernel particle method (RKPM) which is a Galerkin-based meshfree method. The combination of LSTO and RKPM provides a well-defined solid interface which can be maintained on the computational domain by laying particles in the solid and on the level set boundary. The solid and fluid sensitivities are then computed to optimize the fully coupled problem, which is the fundamental challenge in this modular formulation. To identify and remove solid free-floating volumes that emerge during optimization within the flow field, an algorithm based on neighbor information is also introduced. For verification of the approach, benchmarking examples are solved and analyzed based on the assumption of steady state conditions. Beyond the linear elastic solid case considered by the majority of FSI works, we also test our approach with a nonlinear solid with large deformation.
{"title":"Level set topology optimization for fluid-structure interaction using the modified immersed finite element method","authors":"Andreas Neofytou, Lucy T. Zhang, H. Alicia Kim","doi":"10.1016/j.cma.2026.118754","DOIUrl":"https://doi.org/10.1016/j.cma.2026.118754","url":null,"abstract":"This work presents level set topology optimization method (LSTO) for fluid-structure interaction (FSI) problems. For the solution of the FSI problem the modified immersed finite element method (mIFEM) is used. The proposed formulation provides several advantages compared to the existing methods that rely on unified formulations or remeshing approaches. First it results in separate solid and fluid domains, thus allowing any discretization to be used for each physics. This allows for well established independent fluid and solid solvers to be utilized. Further, this modularity is possible without requiring re-meshing, which maintains efficiency especially for the fluid solver. The finite element method (FEM) is used to solve the flow equations efficiently on an Eulerian grid. The Lagrangian solid is analyzed with the reproducing kernel particle method (RKPM) which is a Galerkin-based meshfree method. The combination of LSTO and RKPM provides a well-defined solid interface which can be maintained on the computational domain by laying particles in the solid and on the level set boundary. The solid and fluid sensitivities are then computed to optimize the fully coupled problem, which is the fundamental challenge in this modular formulation. To identify and remove solid free-floating volumes that emerge during optimization within the flow field, an algorithm based on neighbor information is also introduced. For verification of the approach, benchmarking examples are solved and analyzed based on the assumption of steady state conditions. Beyond the linear elastic solid case considered by the majority of FSI works, we also test our approach with a nonlinear solid with large deformation.","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"97 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146098269","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Corrosion-fatigue interaction is a complex and strongly coupled chemo-mechanical degradation process that critically impacts the durability of reinforced concrete (RC) structures exposed to aggressive environments and cyclic loading. This work presents a holistic multiphysics phase-field modeling framework for simulating the full spectrum of coupled physical mechanisms that govern corrosion-fatigue degradation in RC structures. The proposed framework captures: (i) chloride transport and binding in the concrete matrix leading to corrosion initiation, (ii) reactive transport and precipitation of Fe2+ and Fe3+ ions in concrete pores, (iii) pressure accumulation due to rust formation and precipitation around steel reinforcement, (iv) corrosion diffusion and material degradation in steel representing softening due to film rupture and material dissolution, (v) fatigue degradation in steel reinforcement, (vi) fatigue crack propagation in concrete as well as splitting fracture due to corrosion, and (vii) degradation-dependent diffusivity enabling interaction between mechanical cracking and ionic transport. These processes are fully coupled within a unified chemo-mechanical phase-field formulation. Key components of the corrosion and fatigue submodels are validated against experimental data to ensure physical fidelity. The framework is then used to investigate the bidirectional interaction between corrosion and fatigue in both 2D and 3D settings, demonstrating how corrosion accelerates fatigue failure and, conversely, how early-stage fatigue cracking promotes corrosion progression. This comprehensive approach offers a robust tool for assessing service life and designing more durable RC structures under coupled environmental and mechanical loading. The corresponding source codes are openly available at [https://doi.org/10.25835/3duuzvj4 ], allowing reproducibility by interested researchers.
{"title":"Corrosion-fatigue degradation in reinforced concrete structures: A multiphysics phase-field modeling approach","authors":"Abedulgader Baktheer, Manikandan Gopakumar, Ghandi Kenjo, Fadi Aldakheel","doi":"10.1016/j.cma.2025.118693","DOIUrl":"https://doi.org/10.1016/j.cma.2025.118693","url":null,"abstract":"Corrosion-fatigue interaction is a complex and strongly coupled chemo-mechanical degradation process that critically impacts the durability of reinforced concrete (RC) structures exposed to aggressive environments and cyclic loading. This work presents a holistic multiphysics phase-field modeling framework for simulating the full spectrum of coupled physical mechanisms that govern corrosion-fatigue degradation in RC structures. The proposed framework captures: (i) chloride transport and binding in the concrete matrix leading to corrosion initiation, (ii) reactive transport and precipitation of Fe<mml:math altimg=\"si29.svg\"><mml:msup><mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:math> and Fe<mml:math altimg=\"si30.svg\"><mml:msup><mml:mrow></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:math> ions in concrete pores, (iii) pressure accumulation due to rust formation and precipitation around steel reinforcement, (iv) corrosion diffusion and material degradation in steel representing softening due to film rupture and material dissolution, (v) fatigue degradation in steel reinforcement, (vi) fatigue crack propagation in concrete as well as splitting fracture due to corrosion, and (vii) degradation-dependent diffusivity enabling interaction between mechanical cracking and ionic transport. These processes are fully coupled within a unified chemo-mechanical phase-field formulation. Key components of the corrosion and fatigue submodels are validated against experimental data to ensure physical fidelity. The framework is then used to investigate the bidirectional interaction between corrosion and fatigue in both 2D and 3D settings, demonstrating how corrosion accelerates fatigue failure and, conversely, how early-stage fatigue cracking promotes corrosion progression. This comprehensive approach offers a robust tool for assessing service life and designing more durable RC structures under coupled environmental and mechanical loading. The corresponding source codes are openly available at [<ce:inter-ref xlink:href=\"https://doi.org/10.25835/3duuzvj4\" xlink:type=\"simple\">https://doi.org/10.25835/3duuzvj4</ce:inter-ref> ], allowing reproducibility by interested researchers.","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"90 1","pages":""},"PeriodicalIF":7.2,"publicationDate":"2026-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146098293","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-30DOI: 10.1016/j.cma.2026.118773
Cesare Bracco , Francesco Patrizi , Alessandra Sestini
We introduce a novel quadrature strategy for Isogeometric Analysis (IgA) boundary element discretizations, specifically tailored to collocation methods. Thanks to the dimensionality reduction and the natural handling of unbounded domains, boundary integral formulations are particularly appealing in the IgA framework. However, they require the evaluation of boundary integrals whose kernels exhibit singular or nearly singular behavior. Even when the kernel is not singular, its numerical evaluation becomes challenging whenever the integration region lies close to a collocation point. These integrals of polar and nearly singular functions represent the main computational difficulty of IgA-BEM and motivate the development of efficient and accurate quadrature rules. Unlike traditional methods that classify integrals as singular, nearly singular, or regular, our approach employs a desingularizing change of variables that smoothly adapts to the physical distance from singularities in the boundary integral kernels. The transformation intensifies near the polar point and progressively weakens when integrating over portions of the domain that are farther from it, ultimately leaving the integrand unchanged in the limit of a diametrically opposed region. This automatic calibration enhances accuracy and robustness by eliminating the traditional classification step, to which the approximation quality is often highly sensitive. Moreover, integration is performed directly over B-spline supports rather than over individual elements, reducing computational cost, particularly for higher-degree splines. The proposed method is validated through boundary element benchmarks for the three dimensional Stokes problem, where we achieve excellent convergence rates.
{"title":"A smoothly varying quadrature approach for 3D IgA-BEM discretizations: Application to Stokes flow simulations","authors":"Cesare Bracco , Francesco Patrizi , Alessandra Sestini","doi":"10.1016/j.cma.2026.118773","DOIUrl":"10.1016/j.cma.2026.118773","url":null,"abstract":"<div><div>We introduce a novel quadrature strategy for Isogeometric Analysis (IgA) boundary element discretizations, specifically tailored to collocation methods. Thanks to the dimensionality reduction and the natural handling of unbounded domains, boundary integral formulations are particularly appealing in the IgA framework. However, they require the evaluation of boundary integrals whose kernels exhibit singular or nearly singular behavior. Even when the kernel is not singular, its numerical evaluation becomes challenging whenever the integration region lies close to a collocation point. These integrals of polar and nearly singular functions represent the main computational difficulty of IgA-BEM and motivate the development of efficient and accurate quadrature rules. Unlike traditional methods that classify integrals as singular, nearly singular, or regular, our approach employs a desingularizing change of variables that smoothly adapts to the physical distance from singularities in the boundary integral kernels. The transformation intensifies near the polar point and progressively weakens when integrating over portions of the domain that are farther from it, ultimately leaving the integrand unchanged in the limit of a diametrically opposed region. This automatic calibration enhances accuracy and robustness by eliminating the traditional classification step, to which the approximation quality is often highly sensitive. Moreover, integration is performed directly over B-spline supports rather than over individual elements, reducing computational cost, particularly for higher-degree splines. The proposed method is validated through boundary element benchmarks for the three dimensional Stokes problem, where we achieve excellent convergence rates.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118773"},"PeriodicalIF":7.3,"publicationDate":"2026-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146078657","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-29DOI: 10.1016/j.cma.2026.118771
Hualin ZHANG, Weihong ZHANG
In this work, an integrated shape-topology-layout optimization method is developed for stiffened shell structures to match the designs of the skin, openings and stiffeners. The proposed method takes full advantage of IsoGeometric Analysis (IGA), Finite Cell Method (FCM) and Feature-Driven Optimization (FDO). Specifically, the IGA enables accurate structural analysis and facilitates shape optimization of the skin by directly adjusting the control points of Non-Uniform Rational B-Splines (NURBS) skin surface. The FCM enables the use of the fixed structured mesh during the optimization process. The FDO adaptively inserts openings and stiffener features over the NURBS parametric domain for the optimization of skin topology and stiffener layout as easily as for 2D structures. Both implicit and parametric descriptions are jointly used to model openings and stiffener features. The implicit form concerns the level set function (LSF) in favor of the description of topological changes of the skin and the rapid identification of intersections of stiffeners. The parametric form concerns the use of a small number of design variables to produce optimized results with clear and smooth boundaries. Meanwhile, rigidity contributions of stiffeners are merged into the skin surface without the need of mesh conforming between the stiffeners and the skin surface. Representative examples are presented to demonstrate the effectiveness and advantages of the proposed design approach.
{"title":"Shape-topology-layout optimization of stiffened shell structures via the Feature-Driven Optimization (FDO) method","authors":"Hualin ZHANG, Weihong ZHANG","doi":"10.1016/j.cma.2026.118771","DOIUrl":"10.1016/j.cma.2026.118771","url":null,"abstract":"<div><div>In this work, an integrated shape-topology-layout optimization method is developed for stiffened shell structures to match the designs of the skin, openings and stiffeners. The proposed method takes full advantage of IsoGeometric Analysis (IGA), Finite Cell Method (FCM) and Feature-Driven Optimization (FDO). Specifically, the IGA enables accurate structural analysis and facilitates shape optimization of the skin by directly adjusting the control points of Non-Uniform Rational B-Splines (NURBS) skin surface. The FCM enables the use of the fixed structured mesh during the optimization process. The FDO adaptively inserts openings and stiffener features over the NURBS parametric domain for the optimization of skin topology and stiffener layout as easily as for 2D structures. Both implicit and parametric descriptions are jointly used to model openings and stiffener features. The implicit form concerns the level set function (LSF) in favor of the description of topological changes of the skin and the rapid identification of intersections of stiffeners. The parametric form concerns the use of a small number of design variables to produce optimized results with clear and smooth boundaries. Meanwhile, rigidity contributions of stiffeners are merged into the skin surface without the need of mesh conforming between the stiffeners and the skin surface. Representative examples are presented to demonstrate the effectiveness and advantages of the proposed design approach.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118771"},"PeriodicalIF":7.3,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146072454","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-29DOI: 10.1016/j.cma.2026.118756
Yangfan MA , Mitsuteru ASAI , Zheng HAN , Bin SU , Guangqi CHEN
Modeling landslide-generated tsunamis using Smoothed Particle Hydrodynamics (SPH) is hindered by multiphase interactions, large density ratios, and regime transitions that promote excessive numerical dissipation, interface smearing, and divergence errors. This study develops a conservative, interface-aware δR+-SPH framework that couples incremental density diffusion to suppress oscillations with a low-dissipation Riemann flux for energy preservation. Particle disorder is mitigated through a two-stage Optimized Particle Shifting (OPS) for intra-phase regularity and Volume-Conservation Shifting (VCS) to control long-time volume drift. Incompressibility is enforced by combining Velocity-divergence Error Mitigation (VEM) with Hyperbolic/Parabolic Divergence Cleaning (HPDC), which together control both transient and cumulative divergence. The granular phase employs a regime-consistent visco-inertial rheology that unifies rate-dependent friction with effective-pressure regulation across dry, transitional, and submerged regimes. Verification and validation against two-fluid hydrostatics, a rotating patch, immersed granular collapse, and granular slide–tsunami benchmarks confirm that δR+-SPH achieves sharper interfaces, stronger energy and volume conservation, and reduced divergence compared to existing SPH variants. The numerical campaign also yields practical, phase-aware guidelines for shifting and stabilization that balance numerical accuracy and physical fidelity. Collectively, the framework improves the predictive reliability of landslide–tsunami simulations, supporting robust hazard assessment and engineering design.
{"title":"An interface-aware, conservative δR-Plus-SPH for granular slide–water interaction across subaerial–subaqueous regimes","authors":"Yangfan MA , Mitsuteru ASAI , Zheng HAN , Bin SU , Guangqi CHEN","doi":"10.1016/j.cma.2026.118756","DOIUrl":"10.1016/j.cma.2026.118756","url":null,"abstract":"<div><div>Modeling landslide-generated tsunamis using Smoothed Particle Hydrodynamics (SPH) is hindered by multiphase interactions, large density ratios, and regime transitions that promote excessive numerical dissipation, interface smearing, and divergence errors. This study develops a conservative, interface-aware <em>δ</em>R+-SPH framework that couples incremental density diffusion to suppress oscillations with a low-dissipation Riemann flux for energy preservation. Particle disorder is mitigated through a two-stage Optimized Particle Shifting (OPS) for intra-phase regularity and Volume-Conservation Shifting (VCS) to control long-time volume drift. Incompressibility is enforced by combining Velocity-divergence Error Mitigation (VEM) with Hyperbolic/Parabolic Divergence Cleaning (HPDC), which together control both transient and cumulative divergence. The granular phase employs a regime-consistent visco-inertial rheology that unifies rate-dependent friction with effective-pressure regulation across dry, transitional, and submerged regimes. Verification and validation against two-fluid hydrostatics, a rotating patch, immersed granular collapse, and granular slide–tsunami benchmarks confirm that <em>δ</em>R+-SPH achieves sharper interfaces, stronger energy and volume conservation, and reduced divergence compared to existing SPH variants. The numerical campaign also yields practical, phase-aware guidelines for shifting and stabilization that balance numerical accuracy and physical fidelity. Collectively, the framework improves the predictive reliability of landslide–tsunami simulations, supporting robust hazard assessment and engineering design.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118756"},"PeriodicalIF":7.3,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146072453","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：