基于openfoam的高压力低马赫数湍流模型求解器

IF 7.2 2区物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computer Physics Communications Pub Date : 2025-04-02 DOI:10.1016/j.cpc.2025.109600

Danh Nam Nguyen , Chun Sang Yoo

{"title":"基于openfoam的高压力低马赫数湍流模型求解器","authors":"Danh Nam Nguyen , Chun Sang Yoo","doi":"10.1016/j.cpc.2025.109600","DOIUrl":null,"url":null,"abstract":"<div><div>Numerical simulations of non-reacting/reacting flows at supercritical pressure near the critical points with real-fluid models in OpenFOAM often encounter instability and divergence issues unless the solution algorithm incorporates special techniques. In this paper, we develop a novel pressure-based solver, <em>realFluidFoam</em>, tailored for simulations of subsonic turbulent flows at transcritical and supercritical conditions in OpenFOAM. The <em>realFluidFoam</em> solver utilizes unique algorithms to enhance the stability and convergency while taking into account real-fluid effects. Its source code and implementation details are provided to facilitate a comprehensive understanding of integrating real-fluid models into fluid flow simulations in OpenFOAM. The <em>realFluidFoam</em> solver is validated against experimental data by performing large-eddy simulations (LESs) of liquid nitrogen injection and coaxial liquid nitrogen/preheated hydrogen injection under transcritical and supercritical conditions. The LES results show a satisfactory agreement with the experimental data, verifying that the <em>realFluidFoam</em> solver can accurately simulate transcritical and supercritical turbulent fluid flows over the wide range of pressure, especially near the critical points.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> realFluidFoam</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/btzj8b7b8j.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link (OpenFOAM-6 ver.):</em> <span><span>https://github.com/danhnam11/realFluidFoam-6</span><svg><path></path></svg></span></div><div><em>Developer's repository link (OpenFOAM-8 ver.):</em> <span><span>https://github.com/danhnam11/realFluidFoam-8</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> GPLv3</div><div><em>Programming language:</em> C++</div><div><em>Nature of problem:</em> Instability and divergence problems are often encountered in simulations of subsonic fluid flows under high pressure conditions near critical points (i.e., transcritical and supercritical conditions) using real-fluid models in OpenFOAM due to pseudo-boiling effects and high density stratifications.</div><div><em>Solution method:</em> A pressure-based solution method with a modified PIMPLE algorithm is employed to improve the stability while a fast and robust coupling Newton-Bisection algorithm is utilized to guarantee the convergency of fluid flow simulations under transcritical and supercritical conditions in OpenFOAM.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109600"},"PeriodicalIF":7.2000,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"An OpenFOAM-based solver for modeling low Mach number turbulent flows at high pressure with real-fluid effects\",\"authors\":\"Danh Nam Nguyen , Chun Sang Yoo\",\"doi\":\"10.1016/j.cpc.2025.109600\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>Numerical simulations of non-reacting/reacting flows at supercritical pressure near the critical points with real-fluid models in OpenFOAM often encounter instability and divergence issues unless the solution algorithm incorporates special techniques. In this paper, we develop a novel pressure-based solver, <em>realFluidFoam</em>, tailored for simulations of subsonic turbulent flows at transcritical and supercritical conditions in OpenFOAM. The <em>realFluidFoam</em> solver utilizes unique algorithms to enhance the stability and convergency while taking into account real-fluid effects. Its source code and implementation details are provided to facilitate a comprehensive understanding of integrating real-fluid models into fluid flow simulations in OpenFOAM. The <em>realFluidFoam</em> solver is validated against experimental data by performing large-eddy simulations (LESs) of liquid nitrogen injection and coaxial liquid nitrogen/preheated hydrogen injection under transcritical and supercritical conditions. The LES results show a satisfactory agreement with the experimental data, verifying that the <em>realFluidFoam</em> solver can accurately simulate transcritical and supercritical turbulent fluid flows over the wide range of pressure, especially near the critical points.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> realFluidFoam</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/btzj8b7b8j.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link (OpenFOAM-6 ver.):</em> <span><span>https://github.com/danhnam11/realFluidFoam-6</span><svg><path></path></svg></span></div><div><em>Developer's repository link (OpenFOAM-8 ver.):</em> <span><span>https://github.com/danhnam11/realFluidFoam-8</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> GPLv3</div><div><em>Programming language:</em> C++</div><div><em>Nature of problem:</em> Instability and divergence problems are often encountered in simulations of subsonic fluid flows under high pressure conditions near critical points (i.e., transcritical and supercritical conditions) using real-fluid models in OpenFOAM due to pseudo-boiling effects and high density stratifications.</div><div><em>Solution method:</em> A pressure-based solution method with a modified PIMPLE algorithm is employed to improve the stability while a fast and robust coupling Newton-Bisection algorithm is utilized to guarantee the convergency of fluid flow simulations under transcritical and supercritical conditions in OpenFOAM.</div></div>\",\"PeriodicalId\":285,\"journal\":{\"name\":\"Computer Physics Communications\",\"volume\":\"312 \",\"pages\":\"Article 109600\"},\"PeriodicalIF\":7.2000,\"publicationDate\":\"2025-04-02\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Computer Physics Communications\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0010465525001031\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465525001031","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

摘要

在OpenFOAM软件中，对超临界压力下靠近临界点的非反应/反应流动进行实际流体模型数值模拟，往往会遇到不稳定和发散问题，除非在求解算法中加入特殊技术。在本文中，我们开发了一种新的基于压力的求解器realFluidFoam，专门用于在OpenFOAM中模拟跨临界和超临界条件下的亚音速湍流。realFluidFoam求解器采用独特的算法来增强稳定性和收敛性，同时考虑到实际流体效应。提供其源代码和实现细节，以方便全面了解将真实流体模型集成到OpenFOAM中的流体流动模拟中。通过在跨临界和超临界条件下进行液氮注入和同轴液氮/预热氢注入的大涡模拟（LESs），对realFluidFoam求解器进行了实验数据验证。LES计算结果与实验数据吻合较好，验证了realFluidFoam求解器可以准确模拟跨临界和超临界流体在大压力范围内的流动，特别是在临界点附近。程序摘要程序标题：realFluidFoamCPC库链接到程序文件：https://doi.org/10.17632/btzj8b7b8j.1Developer's存储库链接（OpenFOAM-6版本）：https://github.com/danhnam11/realFluidFoam-6Developer's存储库链接（OpenFOAM-8版本）：https://github.com/danhnam11/realFluidFoam-8Licensing规定：gplv3编程语言：c++问题性质：在OpenFOAM中使用实流体模型模拟亚音速流体在靠近临界点（即跨临界和超临界）的高压条件下的流动时，由于伪沸腾效应和高密度分层，经常会遇到不稳定和发散问题。求解方法：采用基于压力的求解方法和改进的PIMPLE算法来提高稳定性，同时采用快速鲁棒的耦合牛顿-对分算法来保证OpenFOAM中跨临界和超临界条件下流体流动模拟的收敛性。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

查看原文本刊更多论文

An OpenFOAM-based solver for modeling low Mach number turbulent flows at high pressure with real-fluid effects

Numerical simulations of non-reacting/reacting flows at supercritical pressure near the critical points with real-fluid models in OpenFOAM often encounter instability and divergence issues unless the solution algorithm incorporates special techniques. In this paper, we develop a novel pressure-based solver, realFluidFoam, tailored for simulations of subsonic turbulent flows at transcritical and supercritical conditions in OpenFOAM. The realFluidFoam solver utilizes unique algorithms to enhance the stability and convergency while taking into account real-fluid effects. Its source code and implementation details are provided to facilitate a comprehensive understanding of integrating real-fluid models into fluid flow simulations in OpenFOAM. The realFluidFoam solver is validated against experimental data by performing large-eddy simulations (LESs) of liquid nitrogen injection and coaxial liquid nitrogen/preheated hydrogen injection under transcritical and supercritical conditions. The LES results show a satisfactory agreement with the experimental data, verifying that the realFluidFoam solver can accurately simulate transcritical and supercritical turbulent fluid flows over the wide range of pressure, especially near the critical points.

Program summary

Program Title: realFluidFoam

CPC Library link to program files: https://doi.org/10.17632/btzj8b7b8j.1

Developer's repository link (OpenFOAM-6 ver.): https://github.com/danhnam11/realFluidFoam-6

Developer's repository link (OpenFOAM-8 ver.): https://github.com/danhnam11/realFluidFoam-8

Licensing provisions: GPLv3

Programming language: C++

Nature of problem: Instability and divergence problems are often encountered in simulations of subsonic fluid flows under high pressure conditions near critical points (i.e., transcritical and supercritical conditions) using real-fluid models in OpenFOAM due to pseudo-boiling effects and high density stratifications.

Solution method: A pressure-based solution method with a modified PIMPLE algorithm is employed to improve the stability while a fast and robust coupling Newton-Bisection algorithm is utilized to guarantee the convergency of fluid flow simulations under transcritical and supercritical conditions in OpenFOAM.

求助全文

通过发布文献求助，成功后即可免费获取论文全文。去求助

来源期刊

Computer Physics Communications 物理-计算机：跨学科应用

CiteScore

12.10

自引率

3.20%

发文量

287

审稿时长

5.3 months

期刊介绍： The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper. Computer Programs in Physics (CPiP) These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged. Computational Physics Papers (CP) These are research papers in, but are not limited to, the following themes across computational physics and related disciplines. mathematical and numerical methods and algorithms; computational models including those associated with the design, control and analysis of experiments; and algebraic computation. Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.