CN105224726B - 结构网格动网格技术用于非结构网格流场求解器的方法 - Google Patents
结构网格动网格技术用于非结构网格流场求解器的方法 Download PDFInfo
- Publication number
- CN105224726B CN105224726B CN201510577643.7A CN201510577643A CN105224726B CN 105224726 B CN105224726 B CN 105224726B CN 201510577643 A CN201510577643 A CN 201510577643A CN 105224726 B CN105224726 B CN 105224726B
- Authority
- CN
- China
- Prior art keywords
- grid
- unstrctured
- structured
- structured grid
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 title claims abstract description 22
- 238000013459 approach Methods 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims abstract description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 230000001052 transient effect Effects 0.000 claims description 3
- 238000005212 lattice dynamic Methods 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 210000004709 eyebrow Anatomy 0.000 description 1
- 210000000720 eyelash Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种结构网格动网格技术应用于非结构网格流场求解器的方法,包括以下步骤:在合理拓扑结构设计基础上,采用代数法与椭圆型方程优化相结合生成初始外形结构网格;分别输出结构网格形式文件和非结构网格形式文件;建立非结构网格点和结构网格点的对应关系,确定每个非结构网格点对应的结构网格点序号;对于边界变形情况,运用结构网格动网格技术更新网格,根据对应关系得到新的非结构网格用于非结构网格流场求解器计算。该方法通过将高效鲁棒的结构网格动网格方法应用于非结构网格流场求解器,提高了非结构网格流场求解器在处理边界运动问题时的效率和能力。
Description
技术领域
本发明涉及一种计算流体力学技术领域的方法,具体讲涉及一种结构网格动网格技术用于非结构网格流场求解器的方法。
背景技术
计算流体动力学(Computational Fluid Dynamics)是通过计算机数值计算和图象显示,对包含有流体流动和热传导等相关物理现象的系统所做的分析。它可以看作是流动基本方程控制下对流动的数值模拟,通过这种数值模拟可以得到极其复杂问题的流场内各个位置上的基本物理量,以及这些物理量随时间的变化情况,并可以进行优化设计。采用CFD求解控制方程时,都是想办法将控制方程在空间区域上进行离散,然后求解得到离散方程组。要想在空间域上离散控制方程,必须使用网格,现在已经发展出多种对各种区域进行离散以生成网格的方法,统称网格生成技术。
网格主要分为结构网格和非结构网格两大类,它们的区别在于节点位置能否用固定的法则有序地命名。结构网格比非结构网格捕捉附面层流动方面更精确,尤其是在高雷诺数的情况下。CFD流场求解器也可分为结构网格流场求解器和非结构网格求解器,但结构网格可以输出为非结构网格形式用于非结构网格流场求解器。含边界运动的非定常流动问题在航空航天领域普遍存在,如飞机机翼的静态变形、颤振、阵风响应等,这种情况下,需要运用动网格技术来实现求解域的运动。
动网格技术可分为非结构网格动网格技术和结构网格动网格技术,非结构网格动网格技术如弹簧法需要迭代计算,变形效率低,变形能力差,适用范围小,而结构网格动网格技术如基于扰动衰减规律的弹性变形技术变形效率高,变形能力强,适用范围广,在工程中得到了广泛使用;提高非结构网格求解器在处理边界运动问题时的效率和能力迫在眉睫。
为解决上述技术问题,本发明提供一种结构网格动网格技术用于非结构网格流场求解器的方法。
发明内容
为克服上述现有技术的不足,本发明提供了一种结构网格动网格技术用于非结构网格流场求解器的方法。
实现上述目的所采用的解决方案为:
一种结构网格动网格技术用于非结构网格流场求解器的方法,所述方法包括以下步骤:
I、在合理拓扑结构设计基础上,生成初始外形结构网格;
II、根据所述初始外形结构网格,分别输出非结构网格形式文件和结构网格形式文件;
III、建立非结构网格点和结构网格点的对应关系,确定每个非结构网格点对应的结构网格点序号;
IV、针对网格边界变形情况,运用结构网格动网格技术更新网格,根据对应关系得到新的非结构网格用于非结构网格流场求解器计算。
优选地,所述步骤I中,采用代数法和椭圆形方程相结合的方法生成初始外形结构网格。
优选地,所述步骤III包括以下步骤:
S301、找出所述非结构网格形式文件中每个非结构网格点对应的所述结构网格形式文件中结构网格点的序号(b,i,j,k);
其中,b表示结构网格点所在的网格块号,i,j,k表示网格点在网格块三个方向的序号;
S302、判断所述非结构网格点和所述结构网格点的对应关系是否成功建立;若所述非结构网格点和所述结构网格点坐标之差小于判断阈值,则认为建立对应关系成功;否则认为建立失败。
优选地,所述非结构网格形式文件为UCD格式,所述结构网格形式文件为plot3d格式。
优选地,所述判断阈值为10-6。
优选地,所述步骤IV中,通过基于扰动衰减规律的结构网格动网格技术更新网格,新的网格点坐标表示如下:
其中,表示初始网格点坐标,表示静态网格点随物面边界作刚性运动的瞬态坐标值,g表示网格点序号的函数;
其中,iw,jw,kw表示物面网格点序号,if,jf,kf表示相应的远边界网格点序号。
与最接近的现有技术相比,本发明具有以下有益效果:
通过将高效鲁棒的结构网格动网格方法应用于非结构网格流场求解器,解决运用传统非结构网格动网格技术存在的变形效率低、变形能力差的技术问题,提高了非结构网格流场求解器在处理边界运动问题时的效率和能力。
附图说明
图1为本发明结构网格动网格技术用于非结构网格流场求解器的方法的流程图;
图2为本实施例中NACA0012翼型计算网格;
图3为本实施例中NACA0012翼型升力系数计算值和实验值的比较。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制,其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他改变。实施例仅代表可能的变化。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。本文中,本发明的这些实施方案可以被单独地或总地用术语“发明”来表示,这仅仅是为了方便,并且如果事实上公开了超过一个的发明,不是要自动的限制该应用的范围为任何单个发明或者发明构思。
图1为本发明结构网格动网格技术用于非结构网格流场求解器的方法的流程图;
本发明提供的一种结构网格动网格技术用于非结构网格流场求解器的方法包括以下步骤:
一、在合理的基础拓扑结构设计基础上,采用代数法与椭圆型方程优化相结合生成初始外形结构网格;
二、根据所述初始外形结构网格,分别输出非结构网格形式文件和结构网格形式文件;
三、建立非结构网格点和结构网格点的对应关系,确定每个非结构网格点对应的结构网格点序号;
四、针对网格边界变形情况,运用结构网格动网格技术更新网格,根据对应关系得到新的非结构网格用于非结构网格流场求解器计算。
步骤一中,采用代数法和椭圆形方程相结合的方法生成初始外形结构网格。
步骤二中,非结构网格形式文件为UCD格式,结构网格形式文件为plot3d格式。
步骤三中,具体包括以下步骤:
S301、找出非结构网格形式文件中每个非结构网格点对应的结构网格形式文件中结构网格点的序号(b,i,j,k);
其中,b表示结构网格点所在的网格块号,i,j,k表示网格点在网格块三个方向的序号;
S302、判断所述非结构网格点和所述结构网格点的对应关系是否成功建立;若所述非结构网格点和所述结构网格点坐标之差小于判断阈值,则认为建立对应关系成功;否则认为建立失败。
判断阈值建议取值为10-6。
步骤四中,通过基于扰动衰减规律的结构网格动网格技术更新网格,新的网格点坐标表示如下:
其中,表示初始网格点坐标,表示静态网格点随物面边界作刚性运动的瞬态坐标值,g表示网格点序号的函数;
其中,iw,jw,kw表示物面网格点序号,if,jf,kf表示相应的远边界网格点序号。
图2为本实施例中NACA0012翼型计算网格、图3为本实施例中NACA0012翼型升力系数计算值和实验值的比较;结合附图2、3对本发明的技术方案做进一步的详细说明,方法流程如图1所示:
以NACA0012翼型俯仰运动为例,本实施例中,提供的方法包括如下:
步骤一、网格拓扑结构设计为O型,采用代数法和椭圆形方程相结合的方法生成初始外形结构网格,计算网格如图2所示,网格数为17712;
步骤二、将生成好的结构网格分别输出成非结构网格UCD文件和结构网格plot3d文件;
步骤三、建立非结构网格点和结构网格点的对应关系,找到UCD文件中每个非结构网格点对应的plot3d文件中结构网格点序号(b,i,j,k),其中,b表示结构网格点所在的网格块号,i,j,k表示网格点在网格块三个方向的序号。
判断非结构网格点和结构网格点的对应关系是否成功建立,方法为判断坐标之差是否小于一个小量;
若非结构网格点和结构网格点坐标之差小于10-6,则认为建立对应关系成功;否则认为建立失败。
步骤四、运动高效鲁棒的基于扰动衰减规律的结构网格动网格技术更新非结构网格,新的非结构网格用于非定常计算。
本实施例中,取计算状态为马赫数0.755,α=2.89°+2.41°sin 2kt,其中,α是翼型迎角,t为时间,k为减缩频率取为0.0808。
图3给出了本实施例中方法计算得到的翼型升力系数和实验值的比较,二者吻合良好,说明本发明方法是可行的。
最后应当说明的是:以上实施例仅用于说明本申请的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解:本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换,但这些变更、修改或者等同替换,均在申请待批的权利要求保护范围之内。
Claims (5)
1.一种结构网格动网格技术用于非结构网格流场求解器的方法,其特征在于:所述方法包括以下步骤:
I、在合理拓扑结构设计基础上,生成初始外形结构网格;
II、根据所述初始外形结构网格,分别输出非结构网格形式文件和结构网格形式文件;
III、建立非结构网格点和结构网格点的对应关系,确定每个非结构网格点对应的结构网格点序号;
IV、针对网格边界变形情况,运用结构网格动网格技术更新网格,根据对应关系得到新的非结构网格用于非结构网格流场求解器计算;
其中,所述步骤III包括以下步骤:
S301、找出所述非结构网格形式文件中每个非结构网格点对应的所述结构网格形式文件中结构网格点的序号(b,i,j,k);
其中,b表示结构网格点所在的网格块号,i,j,k表示网格点在网格块三个方向的序号;
S302、判断所述非结构网格点和所述结构网格点的对应关系是否成功建立;若所述非结构网格点和所述结构网格点坐标之差小于判断阈值,则认为建立对应关系成功;否则认为建立失败。
2.如权利要求1所述的方法,其特征在于:所述步骤I中,采用代数法和椭圆形方程相结合的方法生成所述初始外形结构网格。
3.如权利要求1所述的方法,其特征在于:所述非结构网格形式文件为UCD格式,所述结构网格形式文件为plot3d格式。
4.如权利要求1所述的方法,其特征在于:所述判断阈值为10-6。
5.如权利要求1所述的方法,其特征在于:所述步骤IV中,通过基于扰动衰减规律的结构网格动网格技术更新网格,新的网格点坐标表示如下:
其中,表示初始网格点坐标,表示静态网格点随物面边界作刚性运动的瞬态坐标值,g表示网格点序号的函数;
其中,iw,jw,kw表示物面网格点序号,if,jf,kf表示相应的远边界网格点序号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510577643.7A CN105224726B (zh) | 2015-09-11 | 2015-09-11 | 结构网格动网格技术用于非结构网格流场求解器的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510577643.7A CN105224726B (zh) | 2015-09-11 | 2015-09-11 | 结构网格动网格技术用于非结构网格流场求解器的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105224726A CN105224726A (zh) | 2016-01-06 |
CN105224726B true CN105224726B (zh) | 2018-08-31 |
Family
ID=54993692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510577643.7A Expired - Fee Related CN105224726B (zh) | 2015-09-11 | 2015-09-11 | 结构网格动网格技术用于非结构网格流场求解器的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105224726B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096190A (zh) * | 2016-06-28 | 2016-11-09 | 赫新 | 一种结构/非结构耦合求解cfd的方法 |
CN112329321B (zh) * | 2020-12-17 | 2022-03-01 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种空间离散网格点重归整化方法 |
CN114036872B (zh) * | 2021-11-26 | 2024-05-10 | 天津大学 | 一种基于动网格模拟隧道列车运动的结构网格划分方法 |
CN116227043B (zh) * | 2023-05-10 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞行器数值模拟方法、系统、设备及计算机存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102490909A (zh) * | 2011-11-25 | 2012-06-13 | 中国航天空气动力技术研究院 | 一种飞行器多体分离模拟方法 |
CN103226634A (zh) * | 2013-04-19 | 2013-07-31 | 华南理工大学 | 基于三维动网格的旋喷泵非定常流场的计算方法 |
CN104331584A (zh) * | 2014-11-28 | 2015-02-04 | 北京航空航天大学 | 一种面向附面层计算的二维混合网格自动化生成方法 |
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
-
2015
- 2015-09-11 CN CN201510577643.7A patent/CN105224726B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102490909A (zh) * | 2011-11-25 | 2012-06-13 | 中国航天空气动力技术研究院 | 一种飞行器多体分离模拟方法 |
CN103226634A (zh) * | 2013-04-19 | 2013-07-31 | 华南理工大学 | 基于三维动网格的旋喷泵非定常流场的计算方法 |
CN104331584A (zh) * | 2014-11-28 | 2015-02-04 | 北京航空航天大学 | 一种面向附面层计算的二维混合网格自动化生成方法 |
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
Non-Patent Citations (6)
Title |
---|
A HybridMethod for Dynamic Mesh Generation Based on Radial Basis Functions and Delaunay GraphMapping;Ding L, Guo T, Lu Z.;《Advances in Applied Mathematics and Mechanics》;20150630;第338-356页 * |
An Efficient DynamicMesh GenerationMethod for Complex Multi-Block Structured Grid;Ding L, Lu Z, Guo T.;《Advances in AppliedMathematics and Mechanics》;20131213;第120-134页 * |
Generation of Dynamic Grids and Computation of Unsteady Transonic Flows around Assemblies;LU Zhi-liang;《CHINESE JOURNAL OF AERONAUTICS》;20010228;第1-5页 * |
基于重心坐标插值的网格变形方法研究;孙树立;邱鑫;隋杰;《北京力学会第17届学术年会论文集》;20110108;第439-440页 * |
大飞机跨音速气动弹性问题数值模拟;郭同庆,陆志良;《第十一届全国空气弹性学术交流会》;20091231;第623-627页 * |
气动弹性计算中网格变形方法研究进展;张伟伟,高传强,叶正寅;《航空学报》;20140228;第303-319页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105224726A (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Turbulence modeling of deep dynamic stall at relatively low Reynolds number | |
CN105224726B (zh) | 结构网格动网格技术用于非结构网格流场求解器的方法 | |
Gao et al. | Numerical study on closed-loop control of transonic buffet suppression by trailing edge flap | |
CN107391885A (zh) | 基于有限体积法的剪切滑移动网格方法 | |
Zhu et al. | The time asymmetric pitching effects on the energy extraction performance of a semi-active flapping wing power generator | |
CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
Medida et al. | Numerical investigation of 3-D dynamic stall using delayed detached eddy simulation | |
CN117217108A (zh) | 一种基于cfd的通航飞机飞行仿真的气动力分析方法 | |
CN106326514A (zh) | 模拟小展弦比飞行器不可压缩旋流场的数值方法 | |
US9390205B2 (en) | Vorticity-refinement based numerical method for simulating aircraft wing-tip vortex flows | |
CN109918744A (zh) | 一种基于半拉格朗日和径向基函数的无网格格子波尔兹曼方法 | |
Drofelnik et al. | Rapid calculation of unsteady aircraft loads | |
CN104866712B (zh) | 一种新的投影浸入边界法的速度校正方法 | |
Cottet et al. | A vortex penalization method for flows with moving immersed obstacles | |
Yan et al. | Aerodynamic force reduction of rectangular cylinder using deep reinforcement learning-controlled multiple jets | |
Kurniawan | Numerical study of flutter of a two-dimensional aeroelastic system | |
CN103177162A (zh) | 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法 | |
Zhou et al. | Gridless method for unsteady flows involving moving discrete points and its applications | |
Jahangirian et al. | An implicit solution of the unsteady navier-stokes equations on unstructured moving grids | |
Namura et al. | Multipoint design of vortex generators on a swept infinite-wing under cruise and critical condition | |
Yu et al. | The effects of wing planforms on the aerodynamic performance of thin finite-span flapping wings | |
You et al. | A blended model for simulating massive flow separation and laminar-turbulence transition | |
Yaoyuan et al. | Research on computational grid generation method | |
Demasi et al. | High-Fidelity Simulations of a Flexible Flapping Wing in Forward Flight | |
Somasekhar et al. | Adaptive cloud refinement (ACR)-adaptation in meshless framework |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180831 Termination date: 20210911 |
|
CF01 | Termination of patent right due to non-payment of annual fee |