{"title":"A parallel solver framework for fully implicit monolithic fluid-structure interaction","authors":"Yujie Sun \n (, ), Qingshuang Lu \n (, ), Ju Liu \n (, )","doi":"10.1007/s10409-024-24074-x","DOIUrl":null,"url":null,"abstract":"<div><p>We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction (FSI). The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables, which offers a bridge between computational fluid dynamics (CFD) and computational structural dynamics. The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-α method as a means of providing a robust discrete scheme. In particular, the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI. Based on the chosen fully implicit scheme, we systematically develop a combined suite of nonlinear and linear solver strategies. Invoking a block factorization of the Jacobian matrix, the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage. The first is of the elliptic type, indicating that the algebraic multigrid method serves as a well-suited option. The second exhibits a two-by-two block structure that is analogous to the system arising in CFD. Inspired by prior studies, the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem. Since the number of unknowns matches in both subdomains, it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures. We use two representative FSI benchmarks to demonstrate the robustness, efficiency, and scalability of the overall FSI solver framework. In particular, it is found that the developed FSI solver is comparable to the CFD solver in several aspects, including fixed-size and isogranular scalability as well as robustness.</p><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":7109,"journal":{"name":"Acta Mechanica Sinica","volume":"41 2","pages":""},"PeriodicalIF":3.8000,"publicationDate":"2024-09-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Acta Mechanica Sinica","FirstCategoryId":"5","ListUrlMain":"https://link.springer.com/article/10.1007/s10409-024-24074-x","RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, MECHANICAL","Score":null,"Total":0}
引用次数: 0
Abstract
We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction (FSI). The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables, which offers a bridge between computational fluid dynamics (CFD) and computational structural dynamics. The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-α method as a means of providing a robust discrete scheme. In particular, the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI. Based on the chosen fully implicit scheme, we systematically develop a combined suite of nonlinear and linear solver strategies. Invoking a block factorization of the Jacobian matrix, the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage. The first is of the elliptic type, indicating that the algebraic multigrid method serves as a well-suited option. The second exhibits a two-by-two block structure that is analogous to the system arising in CFD. Inspired by prior studies, the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem. Since the number of unknowns matches in both subdomains, it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures. We use two representative FSI benchmarks to demonstrate the robustness, efficiency, and scalability of the overall FSI solver framework. In particular, it is found that the developed FSI solver is comparable to the CFD solver in several aspects, including fixed-size and isogranular scalability as well as robustness.
期刊介绍:
Acta Mechanica Sinica, sponsored by the Chinese Society of Theoretical and Applied Mechanics, promotes scientific exchanges and collaboration among Chinese scientists in China and abroad. It features high quality, original papers in all aspects of mechanics and mechanical sciences.
Not only does the journal explore the classical subdivisions of theoretical and applied mechanics such as solid and fluid mechanics, it also explores recently emerging areas such as biomechanics and nanomechanics. In addition, the journal investigates analytical, computational, and experimental progresses in all areas of mechanics. Lastly, it encourages research in interdisciplinary subjects, serving as a bridge between mechanics and other branches of engineering and the sciences.
In addition to research papers, Acta Mechanica Sinica publishes reviews, notes, experimental techniques, scientific events, and other special topics of interest.
Related subjects » Classical Continuum Physics - Computational Intelligence and Complexity - Mechanics