计算机科学

计算机科学 ›› 2023, Vol. 50 ›› Issue (11A): 230200019-8.doi: 10.11896/jsjkx.230200019

• 计算机软件&体系架构 • 上一篇    下一篇

一种基于流场物理信息的CFD网格密度优化方法

刘江1,2, 曾志勇1   

  1. 1 重庆邮电大学计算机科学与技术学院 重庆 400065
    2 中国科学院重庆绿色智能技术研究院 重庆 400722
  • 发布日期:2023-11-09
  • 通讯作者: 曾志勇(mailzengzhiyong@gmail.com)
  • 作者简介:(liujiang@cigit.ac.cn)
  • 基金资助:
    国家重点研发计划(2018YFC0116704);中国科学院科技服务网络计划区域重点项目(KFJ-STS-QYZD-2021-01-001);国家自然科学基金(61672488);雷达资料同化关键技术及数值预报客观订正技术研究(E190600801)

CFD Mesh Density Optimization Method Based on Characteristic Flow Distributions

LIU Jiang1,2, ZENG Zhiyong1   

  1. 1 School of Computer Science and Technology,Chongqing University of Posts and Telecommunications,Chongqing 400065,China
    2 Chongqing Institute of Green and Intelligent Technology,Chinese Academy of Sciences,Chongqing 400722,China
  • Published:2023-11-09
  • About author:ZENG Zhiyong,born in 1998,postgra-duate.His main research interests include computational fluid dynamics and CFD mesh optimization.
  • Supported by:
    National Key R & D Program of China(2018YFC0116704),Science and Technology Service Network Initiative(KFJ-STS-QYZD-2021-01-001),National Natural Science Foundation of China(61672488) and Research on Key Techniques of Radar Data Assimilation and Objective Correction of Numerical Forecast(E190600801).

PDF (PC)

160

摘要： CFD计算网格的生成及优化技术是计算流体力学数值计算的关键技术之一,很大程度上决定了数值计算的最终精度及计算效率。CFD模拟实际工程问题所采用的网格规模可达千万级别,在给定的计算时间内获得更高的数值计算精度,是当前CFD急需发展的关键技术。已有研究表明,不同类型的流场误差产生的主要区域与其特征物理量的梯度正相关。有鉴于此,提出了一种基于流场物理信息的网格密度优化方法,使用该方法以及OpenFOAM与cfMesh工具对不可压缩流体案例、燃烧流案例和多相流案例分别进行数值计算。实验结果表明,对于具有不同特征的案例,基于所提方法优化的网格均能够在计算时长差别不大的情况下,比较明显地提升计算精度。

关键词: CFD, 网格优化, 特征物理量, OpenFOAM

Abstract: The generation and optimization of CFD mesh is a key technology in numerical computation of CFD,which significantly determines the final accuracy and efficiency of numerical simulation.The scale of the mesh used in CFD simulation of practical engineering problems can be tens of millions.How to obtain higher numerical accuracy within a given calculation time is a key technology urgently needed to develop CFD.It is shown that the errors of numerical simulation are positively correlated with the gradients of their characteristic physical quantities.In this paper,a mesh density optimization method based on the gradients of chara-cteristic flow distributions is proposed.The mesh density optimization method,the OpenFOAM and cfMesh tools are used to simulate an incompressible fluid,a combustion flow,and a multiphase flow.Simulation results show that for various cases with differentcharacteristics,the optimized mesh by the proposed method can significantly improve the calculation accuracy within almost same calculation time.

Key words: CFD, Mesh optimization, Characteristic physical quantity, OpenFOAM

中图分类号: 

  • TP311
[1]DOMANIN M,VERGARA C,BISSACCO D,et al.Computa-tional fluid dynamics to assess hemodynamic forces in abdominal aortic aneurysm[J].Vascular and Endovascular Surgery,2022,56(3):349-349.
[2]LIAO J W,ZHANG L T,ZHOU X,et al.Improvement of dynamic mesh relaxation method based on double mesh[J].Chinese Journal of Computational Mechanics,2022,39(1):37-41.
[3]WANG B J,LI H,YANG X Y,et al.Application of CFD numerical simulation technology in droplet microfluidic multiphase flow characteristics[J].Chemical Industry Progress,2021,40(4):1715-1735.
[4]YANG D,HAN J,PAN L.A method of heat exchange structure optimization of the cricoid plastic parts[J].Sensors & Transducers,2014,163(1):301-307.
[5]FIELD D A.Laplacian smoothing and delaunay triangulations[J].Communications in Applied Numerical Methods,1988,4(6):709-712.
[6]TSAI H M,WONG A S F,CAI J,et al.Unsteady flow calculations with a parallel multiblock moving mesh algorithm[J].AIAA Journal,2001,39(6):1021-1029.
[7]ZHAO Y,FORHAD A.A general method for simulation of fluid flows with moving and compliant boundaries on unstructured grids[J].Computer Methods in Applied Mechanics and Engineering,2003,192(39/40):4439-4466.
[8]SELIM M M,KOOMULLIL R.Incremental matrix inversionapproach for radial basis function mesh deformation[C]//Proceedings of the Fifteenth Annual Early Career Technical Confe-rence.2015.
[9]LIU X,QIN N,XIA H.Fast dynamic grid deformation based on delaunay graph mapping[J].Journal of Computational Physics,2006,211(2):405-423.
[10]WANG Y,QIN N,ZHAO N.Delaunay graph and radial basis function for fast quality mesh deformation[J].Journal of Computational Physics,2015,294:149-172.
[11]TANG J.Research on adaptive technology of unstructured mesh for engineering[D].Xi’An:Northwestern Polytechnical University,2019.
[12]SEO I H,LEE H J,WI S H,et al.Validation of an air temperature gradient using computational fluid dynamics in a semi-open type greenhouse and determination of kimchi cabbage physiological responses to temperature differences[J].Horticulture,Environment,and Biotechnology,2021,62(5):737-750.
[13]HUANG M,HÖHNE T.Numerical simulation of multicomponent flows with the presence of density gradients for the upgrading of advanced turbulence models[J].Nuclear Engineering and Design,2019,344:28-37.
[14]WANG N H,ZHANG L P,MA R,et al.Influence of unstructured mesh quality on gradient reconstruction and invisco-free flow simulation accuracy[J].Chinese Journal of Computational Mechanics,2017,34(5):555-563.
[15]DE SCHRYVER R,EL CHEIKH K,LESAGE K,et al.Numerical reliability study based on rheological input for bingham paste pumping using a finite volume approach in OpenFOAM[J].Materials,2021,14(17):5011.
[16]LI Y F,WANG Y W,DU C.Performance simulation of plate heat exchanger based on ANSYS ICEM[J].IOP Conference Series:Earth and Environmental Science,2020,546(5):052046.
[17]WILHELMSEN N C A,ANFINSEN H,AAMO O M.Minimum time observer designs for n+m linear hyperbolic systems with unilateral,bilateral or pointwise in-domain sensing[J].European Journal of Control,2021,61:50-67.
[18]LI J,ZHU S.On distributed H1 shape gradient flows in optimalshape design of stokes flows:convergence analysis and numerical applications[J].Journal of Computational Mathematics,2022,40(2):232-258.
[19]YAN C,YU J,XU J L,et al.Progress and prospect of CFD simulation methods[J].Advances in Mechanics,2011,41(5):562-589.
[20]EMILIA F,MARIASOLE B,SANDRO T.Automation techno-logy and computer technology[J].China Radio Electronics Abstract,2010,26(1):152-234.
[21]GHIASI P,NAJAFI G,GHOBADIAN B,et al.Analytical and numerical solution for H-type darrieus wind turbine performance at the tip speed ratio of below one[J].International Journal of Renewable Energy Development,2021,10(2):269-281.
[22]LÖHNER R,BAUM J D.Adaptive h-refinement on 3D unstructured grids for transient problems[J].International Journal for Numerical Methods in Fluids,1992,14(12):1407-1419.
[23]DWYER H A.Grid adaptation for problems in fuid dynamics[J].AIAA Journal,1984,22:1705-1712.
[24]LI J,ZHU S.Shape optimization of the Stokes eigenvalue problem[J].SIAM Journal on Scientific Computing,2022,45(2):A798-A828.
[25]GRASSO G,WU H,ORESTANO S,et al.CFD-based prediction of wall-pressure spectra under a turbulent boundary layer with adverse pressure gradient[J].CEAS Aeronautical Journal,2021,12(1):125-133.
[26]KHAYRULLINA A,VAN HOOFF T,ALANIS RUIZ C,et al.Minimum momentum flux ratio required to prevent air curtain breakthrough in case of cross-curtain pressure gradients:CFD versus analytical equation[J].Building Simulation,2020,13(4):943-960.
[27]LIU M L,ZHANG S Q,LI S N.Analysis of influence of grid on CFD simulation results[J].Hydropower and Pumped Storage,2016,2(4):41-47.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
渝ICP备16013121号-4
地址：重庆市渝北区洪湖西路18号    邮编：401121  
电话：023-63500828（编辑部） 023-67039612（会议/专题） 023-67039623（财务/发票）
E-mail：jsjkx12@163.com
本系统由北京玛格泰克科技发展有限公司设计开发