A robust monolithic nonlinear Newton method for the compressible Reynolds averaged Navier–Stokes Equations

IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Computers & Fluids Pub Date : 2025-03-15 Epub Date: 2025-01-09 DOI:10.1016/j.compfluid.2025.106549
Hulya Sukas, Mehmet Sahin
{"title":"A robust monolithic nonlinear Newton method for the compressible Reynolds averaged Navier–Stokes Equations","authors":"Hulya Sukas,&nbsp;Mehmet Sahin","doi":"10.1016/j.compfluid.2025.106549","DOIUrl":null,"url":null,"abstract":"<div><div>The goal of this work is the development of a fully monolithic nonlinear Newton algorithm for the unstructured vertex-based finite volume algorithm presented in [Akkurt and Sahin, An efficient edge based data structure for the compressible RANS equations on hybrid unstructured meshes. International Journal for Numerical Methods in Fluids, 94:13-31, (2022)]. A special attention is paid for the highly accurate construction of the first- and second-order Jacobian matrices for the Navier–Stokes equations and the one-equation negative Spalart–Allmaras turbulence model. The inviscid flux Jacobian matrices are evaluated exactly using the source code transformations provided by the INRIA Tapenade library. The implementation of the nonlinear Newton method is carried out using the PETSc Scalable Nonlinear Equations Solvers (SNES) with a line search technique. The method can utilize both the Jacobian-free finite difference and direct approaches for the Jacobian vector product. The INRIA pyAMG anisotropic mesh adaptation library has been integrated in order to further improve its numerical accuracy at a lower computational cost. The algorithm has been applied to the two- and three-dimensional mesh convergence test cases from the 4th AIAA CFD High Lift Prediction Workshop in order to demonstrate its robustness in achieving machine precision for a realistic high-lift system. The highly accurate Jacobian evaluation process and high CFL numbers are found to be very critical for achieving machine precision. The mesh adaptation is also proven to be very robust in refining regions with high gradients, which is especially important for the simulations at high angles of attack, where it is rather difficult to determine refinement regions <em>a priori</em>. The anisotropic mesh adaptation studies are carried out for up to 92,993,470 vertices in three-dimensions, and the numerical results indicate very good agreement with the committee’s experimental data.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"289 ","pages":"Article 106549"},"PeriodicalIF":3.0000,"publicationDate":"2025-03-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers & Fluids","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0045793025000106","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2025/1/9 0:00:00","PubModel":"Epub","JCR":"Q3","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0

Abstract

The goal of this work is the development of a fully monolithic nonlinear Newton algorithm for the unstructured vertex-based finite volume algorithm presented in [Akkurt and Sahin, An efficient edge based data structure for the compressible RANS equations on hybrid unstructured meshes. International Journal for Numerical Methods in Fluids, 94:13-31, (2022)]. A special attention is paid for the highly accurate construction of the first- and second-order Jacobian matrices for the Navier–Stokes equations and the one-equation negative Spalart–Allmaras turbulence model. The inviscid flux Jacobian matrices are evaluated exactly using the source code transformations provided by the INRIA Tapenade library. The implementation of the nonlinear Newton method is carried out using the PETSc Scalable Nonlinear Equations Solvers (SNES) with a line search technique. The method can utilize both the Jacobian-free finite difference and direct approaches for the Jacobian vector product. The INRIA pyAMG anisotropic mesh adaptation library has been integrated in order to further improve its numerical accuracy at a lower computational cost. The algorithm has been applied to the two- and three-dimensional mesh convergence test cases from the 4th AIAA CFD High Lift Prediction Workshop in order to demonstrate its robustness in achieving machine precision for a realistic high-lift system. The highly accurate Jacobian evaluation process and high CFL numbers are found to be very critical for achieving machine precision. The mesh adaptation is also proven to be very robust in refining regions with high gradients, which is especially important for the simulations at high angles of attack, where it is rather difficult to determine refinement regions a priori. The anisotropic mesh adaptation studies are carried out for up to 92,993,470 vertices in three-dimensions, and the numerical results indicate very good agreement with the committee’s experimental data.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
可压缩Reynolds平均Navier-Stokes方程的鲁棒整体非线性牛顿法
这项工作的目标是为[Akkurt和Sahin]中提出的基于非结构化顶点的有限体积算法开发一种完全单一的非线性牛顿算法,一种用于混合非结构化网格上可压缩RANS方程的有效的基于边缘的数据结构。流体数值计算方法[j].流体力学学报,2016,34(4):555 - 557。对Navier-Stokes方程和单方程负Spalart-Allmaras湍流模型的一阶和二阶雅可比矩阵的高精度构造给予了特别的关注。使用INRIA Tapenade库提供的源代码转换精确地求出了无粘通量雅可比矩阵。利用PETSc可扩展非线性方程求解器(SNES)的线性搜索技术实现了非线性牛顿法。该方法既可以利用无雅可比矩阵有限差分法,也可以利用直接法求解雅可比矩阵向量积。为了在较低的计算成本下进一步提高数值精度,集成了INRIA pyAMG各向异性网格自适应库。将该算法应用于第四届AIAA CFD高扬程预测研讨会的二维和三维网格收敛测试案例,以验证其在实现实际高扬程系统的机器精度方面的鲁棒性。高精度的雅可比矩阵计算过程和高CFL数是实现机器精度的关键。网格自适应在高梯度的细化区域也被证明是非常鲁棒的,这对于大攻角的模拟尤其重要,因为在大攻角的模拟中,很难先验地确定细化区域。在三维空间中对92,993,470个顶点进行了各向异性网格自适应研究,数值结果与委员会的实验数据吻合良好。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
Computers & Fluids
Computers & Fluids 物理-计算机:跨学科应用
CiteScore
5.30
自引率
7.10%
发文量
242
审稿时长
10.8 months
期刊介绍: Computers & Fluids is multidisciplinary. The term ''fluid'' is interpreted in the broadest sense. Hydro- and aerodynamics, high-speed and physical gas dynamics, turbulence and flow stability, multiphase flow, rheology, tribology and fluid-structure interaction are all of interest, provided that computer technique plays a significant role in the associated studies or design methodology.
期刊最新文献
Variational cell-centered ALE scheme for compressible flows: Improved consistency in the Lagrangian limit Development of a nonlinear reduced-order model for spatio-temporal bluff body airwakes An explicit time-integration framework for accelerated direct numerical simulations of stiff reacting flows Analysis of the impact of numerically induced high-order filters and regularization on LES modelling of an incompressible 3D Taylor-Green vortex flow undergoing laminar-turbulent transition An open-source finite-volume LES code for indoor ventilation and aerosol transport simulations
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1