{"title":"等几何分析的多重网格缩短时间","authors":"R. Tielen, M. Möller, K. Vuik","doi":"10.4995/yic2021.2021.12219","DOIUrl":null,"url":null,"abstract":"Isogeometric Analysis [1] has become increasingly popular as an alternative to the Finite Element Method. Solving the resulting linear systems when adopting higher order B-spline basis functions remains a challenging task, as most (standard) iterative methods have a deteriorating preformance for higher values of the approximation order p.Recently, we succesfully applied p-multigrid methods to discretizations arising in IsogeometricAnalysis [2]. In contrast to h-multigrid methods, where each level of the multigrid hierarchycorresponds to a different mesh width h, the p-multigrid hierarchy is constructed based on different approximation orders. The residual equation is then solved at level p = 1, enabling the use of efficient solution techniques developed for low-order standard FEM. Numerical results show that the number of iterations needed for convergence is independent of both h and p when the p-multigrid method is enhanced with a smoother based on an Incomplete LU factorization with dual treshold (ILUT). However, a slight dependence on the number of patches has been observed for multipatch geometries.Since the resulting system matrix has a block structure in case of a multipatch geometry, weconsider the use of block ILUT as a smoother. Results indicate that the use of block ILUT can be an efficient alternative to ILUT on multipatch geometries within a heterogeneous HPC framework. Prelimenary results for p-multigrid methods adopting a block ILUT smoother will be presented in this talk. Furthermore, we investigate the use of alternative multigrid hierarchies, in particular when considering time-dependent problems.REFERENCES[1] T.J.R. Hughes, J.A. Cottrell and Y. Bazilevs, Isogeometric Analysis: CAD, Finite Elements,NURBS, Exact Geometry and Mesh Refinement, Computer Methods in Applied Mechanicsand Engineering, 194, 4135 - 4195, 2005[2] R.Tielen, M. Möller, D. Göddeke and C. Vuik, p-multigrid methods and their comparison toh-multigrid methods within Isogeometric Analysis, Computer Methods in Applied Mechanicsand Engineering, 372, 2020","PeriodicalId":406819,"journal":{"name":"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference","volume":"1 1","pages":"0"},"PeriodicalIF":0.0000,"publicationDate":"2021-07-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"2","resultStr":"{\"title\":\"Multigrid Reduced in Time for Isogeometric Analysis\",\"authors\":\"R. Tielen, M. Möller, K. Vuik\",\"doi\":\"10.4995/yic2021.2021.12219\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Isogeometric Analysis [1] has become increasingly popular as an alternative to the Finite Element Method. Solving the resulting linear systems when adopting higher order B-spline basis functions remains a challenging task, as most (standard) iterative methods have a deteriorating preformance for higher values of the approximation order p.Recently, we succesfully applied p-multigrid methods to discretizations arising in IsogeometricAnalysis [2]. In contrast to h-multigrid methods, where each level of the multigrid hierarchycorresponds to a different mesh width h, the p-multigrid hierarchy is constructed based on different approximation orders. The residual equation is then solved at level p = 1, enabling the use of efficient solution techniques developed for low-order standard FEM. Numerical results show that the number of iterations needed for convergence is independent of both h and p when the p-multigrid method is enhanced with a smoother based on an Incomplete LU factorization with dual treshold (ILUT). However, a slight dependence on the number of patches has been observed for multipatch geometries.Since the resulting system matrix has a block structure in case of a multipatch geometry, weconsider the use of block ILUT as a smoother. Results indicate that the use of block ILUT can be an efficient alternative to ILUT on multipatch geometries within a heterogeneous HPC framework. Prelimenary results for p-multigrid methods adopting a block ILUT smoother will be presented in this talk. Furthermore, we investigate the use of alternative multigrid hierarchies, in particular when considering time-dependent problems.REFERENCES[1] T.J.R. Hughes, J.A. Cottrell and Y. Bazilevs, Isogeometric Analysis: CAD, Finite Elements,NURBS, Exact Geometry and Mesh Refinement, Computer Methods in Applied Mechanicsand Engineering, 194, 4135 - 4195, 2005[2] R.Tielen, M. Möller, D. Göddeke and C. Vuik, p-multigrid methods and their comparison toh-multigrid methods within Isogeometric Analysis, Computer Methods in Applied Mechanicsand Engineering, 372, 2020\",\"PeriodicalId\":406819,\"journal\":{\"name\":\"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference\",\"volume\":\"1 1\",\"pages\":\"0\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2021-07-07\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"2\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/10.4995/yic2021.2021.12219\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.4995/yic2021.2021.12219","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Multigrid Reduced in Time for Isogeometric Analysis
Isogeometric Analysis [1] has become increasingly popular as an alternative to the Finite Element Method. Solving the resulting linear systems when adopting higher order B-spline basis functions remains a challenging task, as most (standard) iterative methods have a deteriorating preformance for higher values of the approximation order p.Recently, we succesfully applied p-multigrid methods to discretizations arising in IsogeometricAnalysis [2]. In contrast to h-multigrid methods, where each level of the multigrid hierarchycorresponds to a different mesh width h, the p-multigrid hierarchy is constructed based on different approximation orders. The residual equation is then solved at level p = 1, enabling the use of efficient solution techniques developed for low-order standard FEM. Numerical results show that the number of iterations needed for convergence is independent of both h and p when the p-multigrid method is enhanced with a smoother based on an Incomplete LU factorization with dual treshold (ILUT). However, a slight dependence on the number of patches has been observed for multipatch geometries.Since the resulting system matrix has a block structure in case of a multipatch geometry, weconsider the use of block ILUT as a smoother. Results indicate that the use of block ILUT can be an efficient alternative to ILUT on multipatch geometries within a heterogeneous HPC framework. Prelimenary results for p-multigrid methods adopting a block ILUT smoother will be presented in this talk. Furthermore, we investigate the use of alternative multigrid hierarchies, in particular when considering time-dependent problems.REFERENCES[1] T.J.R. Hughes, J.A. Cottrell and Y. Bazilevs, Isogeometric Analysis: CAD, Finite Elements,NURBS, Exact Geometry and Mesh Refinement, Computer Methods in Applied Mechanicsand Engineering, 194, 4135 - 4195, 2005[2] R.Tielen, M. Möller, D. Göddeke and C. Vuik, p-multigrid methods and their comparison toh-multigrid methods within Isogeometric Analysis, Computer Methods in Applied Mechanicsand Engineering, 372, 2020