Computer Methods in Applied Mechanics and Engineering最新文献

筛选
英文 中文
The generalized shifted boundary method for geometry-parametric PDEs and time-dependent domains 几何参数偏微分方程和时变域的广义位移边界法
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-15 DOI: 10.1016/j.cma.2026.118748
Oriol Colomés , Jan Modderman , Guglielmo Scovazzi
{"title":"The generalized shifted boundary method for geometry-parametric PDEs and time-dependent domains","authors":"Oriol Colomés ,&nbsp;Jan Modderman ,&nbsp;Guglielmo Scovazzi","doi":"10.1016/j.cma.2026.118748","DOIUrl":"10.1016/j.cma.2026.118748","url":null,"abstract":"<div><div>Many engineering and scientific problems require the solution of partial differential equations in complex geometries. Often, these problems involve parametrized geometries, e.g. design optimization, or moving domains, e.g. fluid-structure interaction problems. For such cases, traditional methods based on body-fitted grids require time-consuming mesh generation or re-meshing techniques. Unfitted finite element methods, e.g. CutFEM of AgFEM, are appealing techniques that address these challenges. However, they require ad-hoc integration methods and stabilization techniques to prevent instabilities for small cut cells. Recently, the Shifted Boundary Method (SBM), was introduced to prevent integration over cut cells and small cut-cell instabilities. An extension of the SBM was recently introduced, the Weighted Shifted Boundary Method (WSBM), where the variational form is weighted by the elemental active volume fraction, improving discrete mass/momentum conservation properties in simulations with moving domains. In this work we introduce the Generalized Shifted Boundary Method (GSBM), a geometry-agnostic generalization of the SBM and WSBM formulations that avoids the need of redefinition of integration domains and finite element spaces. The GSBM enables a unified formulation for problems with evolving geometries, supports gradient-based optimization of problems with varying geometries including topological changes, and unifies SBM, WSBM, and optimal-surrogate variants within a single framework. In this work we describe the formulation, and corresponding tests, for three model problems, namely: the Poisson problem, linear elasticity and transient Stokes flow.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118748"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980092","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}
引用次数: 0
Discovering neural cohesive zone laws from displacement fields 从位移场发现神经内聚带规律
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-17 DOI: 10.1016/j.cma.2026.118733
Georgios Barkoulis Gavris, WaiChing Sun
{"title":"Discovering neural cohesive zone laws from displacement fields","authors":"Georgios Barkoulis Gavris,&nbsp;WaiChing Sun","doi":"10.1016/j.cma.2026.118733","DOIUrl":"10.1016/j.cma.2026.118733","url":null,"abstract":"<div><div>Conventional inverse problems for cohesive zones often utilize homogenized responses of the effective media to identify a fixed set of material parameters prescribed <em>a priori</em>. However, the mixed-mode loading conditions of composites or natural materials may exhibit interfacial relations that are difficult to anticipate. This article presents a model-discovery framework for directly identifying cohesive zone models inferred from displacement fields across the interface, without fixing on a specific form of equations. We develop a differentiable version of the Material Point Method (MPM) with interface elements formulated to capture the traction-separation law at a pre-existing crack or bonded interface. Ensuring the differentiability of the MPM solver enables us to backpropagate the mismatch between simulated and measured (e.g., DIC/DVC) displacement fields through the time integrator and interface physics. Using only kinematics and equilibrium as constraints, numerical experiments suggest that the method may recover (i) a Mode-I traction-separation curve in a double-cantilever-beam test and (ii) a mixed-mode law for a circular interface shear test. These numerical results demonstrate that displacement-only experiments, combined with a differentiable solver, offer a promising pathway for identifying rich and potentially nonparametric cohesive laws.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118733"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980090","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}
引用次数: 0
An efficient three-step subgrid stabilized method for the steady natural convection equations 稳定自然对流方程的一种有效的三步子网格稳定方法
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-07 DOI: 10.1016/j.cma.2025.118710
Bo Zheng , Chuanqin Zheng , Yueqiang Shang , Yinnian He
{"title":"An efficient three-step subgrid stabilized method for the steady natural convection equations","authors":"Bo Zheng ,&nbsp;Chuanqin Zheng ,&nbsp;Yueqiang Shang ,&nbsp;Yinnian He","doi":"10.1016/j.cma.2025.118710","DOIUrl":"10.1016/j.cma.2025.118710","url":null,"abstract":"<div><div>This article is concerned with an efficient three-step subgrid stabilized method for the steady natural convection equations with moderate to high Rayleigh numbers in which two levels of finite element meshes are required. Within this method, we first solve one small, nonlinear coarse mesh natural convection problem with subgrid stabilizations, and then solve two large, subgrid stabilized and Newton-linearized fine mesh problems which have identical coefficient matrices with only different right-hand sides. This indicates that we avoid reassembly of the linearized problem in the third step, making the presented method easy-to-implement. Among them, both the coarse mesh and fine mesh problems are stabilized by the subgrid-scale models defined by two elliptic projections into lower-order finite element spaces of the velocity and temperature to improve the stability and convergence and thus, our present method has broad potential applications in simulating moderate to high Rayleigh number flows. Under the weak uniqueness condition, error bounds of the approximate solutions from the proposed method are strictly established. Scalings of the algorithmic parameters concerning the mesh sizes and stabilization parameters are also derived. In the end, manufactured solution examples are carried out to conform the theoretical analysis, showing a higher precision of the stabilized solutions calculated by our present three-step subgrid stabilized method than that of its counterpart method only staying at the second step. Besides, some numerical simulations and results are presented for the buoyancy-driven square cavity flow, sinusoidal hot cylinder flow and the isolated island problem with the Rayleigh numbers up to <span><math><mrow><mi>R</mi><mi>a</mi><mo>=</mo><msup><mn>10</mn><mn>7</mn></msup></mrow></math></span> representative of the convection-dominated regime to illustrate the high efficiency and effectiveness of our present method.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118710"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145940982","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}
引用次数: 0
Structural optimization of a stack of elastic rings under gravity 重力作用下一叠弹性环的结构优化
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-14 DOI: 10.1016/j.cma.2025.118698
Luis Andres Mollericon Titirico , Sylvain Lefebvre , Ole Sigmund , Jonàs Martínez
{"title":"Structural optimization of a stack of elastic rings under gravity","authors":"Luis Andres Mollericon Titirico ,&nbsp;Sylvain Lefebvre ,&nbsp;Ole Sigmund ,&nbsp;Jonàs Martínez","doi":"10.1016/j.cma.2025.118698","DOIUrl":"10.1016/j.cma.2025.118698","url":null,"abstract":"<div><div>We study the structural optimization of a stack of elastic rings subjected to gravity. The aim is to optimize for minimum or maximum volume enclosing structures while preventing their collapse under their own weight. We formulate the problem using a parameterization of the cross-section geometry of the rings and a tailored optimization scheme that considers axisymmetric finite elements. As demonstrated through numerical examples, our method produces structurally sound shapes that can then be fabricated by stacking ring components or through extrusion-based manufacturing.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118698"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145962659","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}
引用次数: 0
Physics informed surface autoencoders for thin shell analysis 物理通知表面自编码器薄壳分析
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-24 DOI: 10.1016/j.cma.2026.118764
Aswanth Thani , Adrian Buganza Tepole
{"title":"Physics informed surface autoencoders for thin shell analysis","authors":"Aswanth Thani ,&nbsp;Adrian Buganza Tepole","doi":"10.1016/j.cma.2026.118764","DOIUrl":"10.1016/j.cma.2026.118764","url":null,"abstract":"<div><div>We present a physics-informed surface autoencoder (PISA) framework for Kirchhoff-Love thin shell analysis. The method constructs global <em>C</em><sup>1</sup> surface parameterizations directly from unstructured point clouds for both single-patch surfaces homeomorphic to disks, and multi-patch parameterizations for closed genus-zero surfaces. In the multi-patch case, a classification network assigns probabilistic labels to points, and the autoencoder learns overlapping charts with smooth transitions, ensuring global <em>C</em><sup>1</sup> continuity. With the learned parameterizations, we introduce a decoder for the displacement field and compute differential geometric quantities such as the metric and second fundamental form in the reference and deformed surfaces. Then, we enforce equilibrium by minimizing the total potential energy. The approach is validated on classical shell benchmarks, including the Scordelis-Lo roof, pinched cylinder, and hemisphere under pressure. We showcase the flexibility of the framework with complex geometries such as the Stanford Bunny and dura mater. Compared with traditional spline-based parameterizations and existing machine learning approaches, PISA offers a pipeline for generating smooth surface maps for complex geometries and integrates the surface representation into the physics-informed solver. Importantly, the thin shell analysis pipeline proposed works directly with unstructured point cloud data. Thus, this PISA framework’s potential applications range from engineering structures to biological membranes such as heart valves, skin, and dura mater.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118764"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023859","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}
引用次数: 0
An interface-aware, conservative δR-Plus-SPH for granular slide–water interaction across subaerial–subaqueous regimes 一个界面感知的,保守的δ r - + sph,用于颗粒滑块-水在水下-水下的相互作用
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-29 DOI: 10.1016/j.cma.2026.118756
Yangfan MA , Mitsuteru ASAI , Zheng HAN , Bin SU , Guangqi CHEN
{"title":"An interface-aware, conservative δR-Plus-SPH for granular slide–water interaction across subaerial–subaqueous regimes","authors":"Yangfan MA ,&nbsp;Mitsuteru ASAI ,&nbsp;Zheng HAN ,&nbsp;Bin SU ,&nbsp;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-04-15","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}
引用次数: 0
Generalized Eigenvalue stabilization for immersed explicit dynamics 浸入式显式动力学的广义特征值镇定
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-17 DOI: 10.1016/j.cma.2026.118727
Tim Bürchner , Lars Radtke , Sascha Eisenträger , Alexander Düster , Ernst Rank , Stefan Kollmannsberger , Philipp Kopp
{"title":"Generalized Eigenvalue stabilization for immersed explicit dynamics","authors":"Tim Bürchner ,&nbsp;Lars Radtke ,&nbsp;Sascha Eisenträger ,&nbsp;Alexander Düster ,&nbsp;Ernst Rank ,&nbsp;Stefan Kollmannsberger ,&nbsp;Philipp Kopp","doi":"10.1016/j.cma.2026.118727","DOIUrl":"10.1016/j.cma.2026.118727","url":null,"abstract":"<div><div>Explicit time integration for immersed finite element discretizations severely suffers from the influence of poorly cut elements. In this contribution, we propose a generalized eigenvalue stabilization (GEVS) strategy for the element mass matrices of cut elements to cure their adverse impact on the critical time step size of the global system. We use spectral basis functions, specifically <em>C</em><sup>0</sup> continuous Lagrangian interpolation polynomials defined on Gauss-Lobatto-Legendre (GLL) points, which, in combination with its associated GLL quadrature rule, yield high-order convergent diagonal mass matrices for uncut elements. Moreover, considering cut elements, we combine the proposed GEVS approach with the finite cell method to guarantee definiteness of the system matrices. However, the proposed GEVS stabilization can directly be applied to other immersed boundary finite element methods. Numerical experiments demonstrate that the stabilization strategy achieves optimal convergence rates and recovers critical time step sizes of equivalent boundary-conforming discretizations. This also holds in the presence of weakly enforced Dirichlet boundary conditions using either Nitsche’s method or penalty formulations.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118727"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145995462","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}
引用次数: 0
Reduced-order initialization model for accelerating topology optimization in thermal-fluid system 加速热流体系统拓扑优化的降阶初始化模型
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-02-06 DOI: 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,&nbsp;Melvan Tan Kian Hao,&nbsp;Dong Chen,&nbsp;Lailai Zhu,&nbsp;Poh Seng Lee","doi":"10.1016/j.cma.2026.118794","DOIUrl":"10.1016/j.cma.2026.118794","url":null,"abstract":"<div><div>Topology Optimization (TO) in thermal-fluid systems inevitably requires repeated numerical simulations, resulting in prohibitive computational costs. To overcome this challenge, Reduced-Order Initialization Model (ROIM) is proposed to capture latent-space relationships between previous and current TO iterations, thereby predicting high-quality initial fields that accelerate numerical convergence and reduce the overall TO runtime. The prediction from the solid-fluid configuration to the physical fields is achieved by combining Weighted Proper Orthogonal Decomposition (WPOD) with reduced-order mapping. As a demonstration, ROIM is applied to four representative TO cases, including 2D heat transfer optimization (22.55% runtime reduction), 2D multi-inlet flow optimization (35.28%), 3D flow optimization (33.68%), and 3D heat transfer optimization (32.56%), consistently achieving substantial computational time savings. In terms of prediction accuracy, ROIM attains a low mean error of 2.09% in 2D heat transfer optimization and exhibits robust performance under boundary-condition fluctuations. Furthermore, in the 3D comparative cases, ROIM achieves computational time reductions ranging from 28.89% to 54.32% as the model fidelity increases, demonstrating clear scalability.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118794"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","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}
引用次数: 0
A variational phase-field framework for multiphysics modeling of delayed hydride cracking in zirconium alloys 锆合金延迟氢化物裂纹多物理场模型的变分相场框架
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-01-20 DOI: 10.1016/j.cma.2026.118765
Wei Li , Han Zhao , Minghua Chi , Vincent Beng Chye Tan
{"title":"A variational phase-field framework for multiphysics modeling of delayed hydride cracking in zirconium alloys","authors":"Wei Li ,&nbsp;Han Zhao ,&nbsp;Minghua Chi ,&nbsp;Vincent Beng Chye Tan","doi":"10.1016/j.cma.2026.118765","DOIUrl":"10.1016/j.cma.2026.118765","url":null,"abstract":"<div><div>Delayed hydride cracking (DHC) poses a significant integrity threat to zirconium nuclear cladding, arising from the cyclic interplay of hydrogen diffusion and hydride precipitation/dissolution, which leads to characteristic intermittent subcritical crack advance. This work develops a thermodynamically consistent variational phase-field framework that, in a unified formulation, couples hydrogen diffusion, stress- and temperature-regulated hydride evolution, thermo-elastoplasticity, and a ductile-to-brittle fracture transition. The model reproduces key experimental observations across temperatures—intermittent crack advance, DHC velocity, striation spacing, and incubation time—and quantifies how stress and temperature jointly govern hydride nucleation and cracking. Mechanistically, hydrostatic tension reduces the terminal solid solubility for precipitation/dissolution and concentrates hydrogen; ensuing precipitation relaxes and redistributes stresses, establishing a self-sustaining feedback loop that triggers hydride precipitation–fracture–dissolution–reprecipitation cycles. Parametric studies of pre-cracked cladding reveal a stress-dependent transition between diffusion-controlled and precipitation-controlled initiation. High applied stress induces pronounced thermal sensitivity via stress-assisted hydrogen accumulating, whereas low stress exhibits a weaker temperature response that manifests only above a critical hydride fraction. Beyond DHC, the framework is readily extensible to simulations of oxide–hydride synergistic delayed cracking under reactor-relevant conditions. It thus provides a physics-based foundation for mechanism identification, threshold assessment, and life prediction of zirconium cladding, and establishes a platform for future extensions to oxide–hydride interaction.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118765"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146014810","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}
引用次数: 0
Speeding up an unsteady flow simulation by adaptive BDDC and Krylov subspace recycling 基于自适应BDDC和Krylov子空间循环加速非定常流场模拟
IF 7.3 1区 工程技术
Computer Methods in Applied Mechanics and Engineering Pub Date : 2026-04-15 Epub Date: 2026-02-02 DOI: 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 ,&nbsp;Jan Papež ,&nbsp;Jakub Šístek","doi":"10.1016/j.cma.2026.118788","DOIUrl":"10.1016/j.cma.2026.118788","url":null,"abstract":"<div><div>We deal with accelerating the solution of a sequence of large linear systems solved by preconditioned conjugate gradient method (PCG). The sequence originates from time-stepping within a simulation of an unsteady incompressible flow. We apply a pressure correction scheme and focus on the solution of the Poisson problem for the pressure corrector. Its scalable solution presents the main computational challenge in many applications. The right-hand side of the problem changes in each time step, while the system matrix is constant and symmetric positive definite. The acceleration techniques are studied on a representative problem of flow around a unit sphere. Our baseline approach is based on a parallel solution of each problem in the sequence by nonoverlapping domain decomposition method. The interface problem is solved by PCG with the three-level BDDC preconditioner. As a preliminary step, an appropriate stopping criterion for the PCG iterations is chosen. Next, two techniques for accelerating the solution are gradually added to the baseline approach. Deflation is used within PCG with several approaches to Krylov subspace recycling. Finally, we add the adaptive selection of the coarse space within the three-level BDDC method. The paper is rich in experiments with careful measurements of computational times on a parallel supercomputer. The combination of the acceleration techniques eventually leads to saving more than 40 % of the computational time.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118788"},"PeriodicalIF":7.3,"publicationDate":"2026-04-15","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}
引用次数: 0
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
相关产品
×
本文献相关产品
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信
小红书