CN114880890B - 一种远场边界条件多窗口处理方法、设备及介质 - Google Patents

一种远场边界条件多窗口处理方法、设备及介质 Download PDF

Info

Publication number
CN114880890B
CN114880890B CN202210809742.3A CN202210809742A CN114880890B CN 114880890 B CN114880890 B CN 114880890B CN 202210809742 A CN202210809742 A CN 202210809742A CN 114880890 B CN114880890 B CN 114880890B
Authority
CN
China
Prior art keywords
boundary
grid
flow
field
far
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.)
Active
Application number
CN202210809742.3A
Other languages
English (en)
Other versions
CN114880890A (zh
Inventor
赵凡
齐琛
张健
王显焯
武文军
马洪林
王应宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202210809742.3A priority Critical patent/CN114880890B/zh
Publication of CN114880890A publication Critical patent/CN114880890A/zh
Application granted granted Critical
Publication of CN114880890B publication Critical patent/CN114880890B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种远场边界条件多窗口处理方法、设备及介质,属于计算流体力学领域,包括步骤:在求解器输入参数,通过输入的参数对每个网格分别进行边界条件的赋值实现边界条件的重定义,再使用不可压流动求解算法对流场进行计算。本发明解决了不可压缩流动数值模拟中难以使用“远场”边界的难题,提升了求解器的通用性和用户友好性。

Description

一种远场边界条件多窗口处理方法、设备及介质
技术领域
本发明涉及计算流体力学领域,更为具体的,涉及一种远场边界条件多窗口处理方法、设备及介质。
背景技术
CFD仿真模拟时经常设有“远场”边界条件,该类边界条件主要用于可压缩流动问题,即整个计算区域的外边界是一整块网格。此类边界条件在网格生成、改变来流方向时非常方便,但需要很大的计算域。在不可压缩流动问题的研究中,也时常会遇到这种圆形/球形的远场网格,但由于计算域通常比较小,远场类边界条件并不适合不可压缩流动模拟,通常需要在外边界网格中明确区分出入口、出口边界。为了在不可压流动模拟中可以处理此类网格,需要一种技术,可以在模拟不可压缩流动问题时针对性设置“远场”边界条件。
现有的边界条件处理技术基本都是单窗口技术,即,相同区域的同种类型的网格整体定义一个边界条件,例如Fluent等商业软件均是如此操作。
单窗口边界条件处理技术满足大部分边界条件处理需求,然而,在求解不可压缩流动时,“远场”边界就无法使用。此时需要将“远场”边界转化为“入口”、“出口”、“对称”等边界条件才可以,然而单窗口技术导致一个区域网格的边界条件是整体定义的,故而无法将边界条件进行转化。
因此,需要一种新的边界条件处理技术来改进现有的缺陷,使得边界条件可以转化。
发明内容
本发明的目的在于克服现有技术的不足,提供一种远场边界条件多窗口处理方法、设备及介质,解决了不可压缩流动数值模拟中难以使用“远场”边界的难题,提升了求解器的通用性和用户友好性。
本发明的目的是通过以下方案实现的:
一种远场边界条件多窗口处理方法,包括步骤:
在求解器输入参数,通过输入的参数对每个网格分别进行边界条件的赋值实现边界条件的重定义,再使用不可压流动求解算法对流场进行计算。
进一步地,所述对每个网格分别进行边界条件的赋值实现边界条件的重定义,包括子步骤:利用每个网格的网格单元流量的正负,在求解器内部定义入口、出口、对称类边界条件。
进一步地,所述对每个网格分别进行边界条件的赋值实现边界条件的重定义,包括子步骤:
步骤一,将边界条件分别定义到每个网格的网格单元,即多窗口化处理边界条件;
步骤二,判断各个单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”。
进一步地,所述网格为任意形状的边界网格。
进一步地,所述不可压流动求解算法包SIMPLE算法、SIMPLEC算法、SIMPLEX算法、SIMPLER算法、PISO算法中的任一种。
进一步地,包括子步骤:
S1,创建出一套流场计算的网格并定义其边界条件;
S2,将网格导入到求解器中;
S3,对网格进行前处理;
S4,输入参数,通过输入的参数对每个网格进行边界条件的赋值实现边界条件的重定义;
S5,使用不可压流动求解算法对流场进行计算。
进一步地,在步骤S5之后,包括步骤:S6,对收敛的流场进行后处理显示。
进一步地,在步骤S4中,包括子步骤:
S41,将边界条件分别定义到每个网格单元,即多窗口化处理边界条件;
S42,判断各个网格单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;
S43,若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”。
一种计算机设备,所述计算机设备包括处理器和存储器,所述存储器中存储有计算机程序,当所述计算机程序被所述处理器加载并执行如上任一项所述的方法。
一种计算机可读存储介质,在可读存储介质中存储有计算机程序,所述计算机程序被处理器加载并执行如上任一项所述的方法。
本发明的有益效果包括:
本发明通过分别给每个边界网格赋了边界条件,使得“远场”边界可以转化为传统的“入口”、“出口”、“对称”等边界条件,成功解决了不可压缩流动数值模拟中难以使用“远场”边界的难题,提升了求解器的通用性和用户友好性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种适用于不可压流动求解的远场边界条件多窗口处理方法计算流程示意图;
图2为带有“远场”边界的0012翼型网格划分示意图;
图3为0012翼型模型计算结果压力场示意图,图3中cp代表压力系数,图3中cp下面的数值代表不同压力系数数值。
具体实施方式
本说明书中所有实施例公开的所有特征,或隐含公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合和/或扩展、替换。
本发明实施例提供一种适用于不可压流动求解的远场边界条件多窗口处理方法,包括以下步骤:
1)创建出一套流场计算的网格并定义其边界条件;
2)将网格导入到求解器中;
3)对网格进行前处理;
4)输入计算参数,通过输入的参数对边界条件进行重定义,具体如下:
① 将边界条件分别定义到每个网格的网格单元,即多窗口化处理边界条件;
② 判断各个单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;
③ 若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”;
5)使用不可压流动求解算法对流场进行计算;
6)对收敛的流场进行后处理显示。
本发明实施例提供的远场边界条件多窗口处理办法,可以将“远场”边界处理为传统的“入口”、“出口”、“对称”边界,成功解决了不可压缩流动数值模拟中难以使用“远场”边界的难题,提升了求解器的通用性和用户友好性。
实施例1:一种远场边界条件多窗口处理方法,在具体实施过程中,包括步骤:
在求解器输入参数,不对一个区域的网格整体进行边界条件的赋值,而是通过输入的参数对每个网格分别进行边界条件的赋值实现边界条件的重定义,再使用不可压流动求解算法对流场进行计算。
实施例2:在实施例1的基础上,所述对每个网格分别进行边界条件的赋值实现边界条件的重定义,包括子步骤:利用每个网格的网格单元流量的正负,在求解器内部定义入口、出口、对称类边界条件。
实施例3:在实施例1的基础上,所述对每个网格分别进行边界条件的赋值实现边界条件的重定义,包括子步骤:
步骤一,将边界条件分别定义到每个网格的网格单元,即多窗口化处理边界条件;
步骤二,判断各个单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”。
实施例4:在实施例2~实施例3的基础上,所述边界网格为任意形状的边界网格。
实施例5:在实施例1的基础上,所述不可压流动求解算法包SIMPLE算法、SIMPLEC算法、SIMPLEX算法、SIMPLER算法、PISO算法中的任一种。
实施例6:在实施例1的基础上,包括子步骤:
S1,创建出一套流场计算的网格并定义其边界条件;
S2,将网格导入到求解器中;
S3,对网格进行前处理;
S4,输入参数,通过输入的参数对每个网格进行边界条件的赋值实现边界条件的重定义;
S5,使用不可压流动求解算法对流场进行计算。
实施例7:在实施例6的基础上,在步骤S5之后,包括步骤:S6,对收敛的流场进行后处理显示。
实施例8:在实施例6的基础上,在步骤S4中,包括子步骤:
S41,将边界条件分别定义到每个网格单元,即多窗口化处理边界条件;
S42,判断各个网格单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;
S43,若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”。
实施例9:一种计算机设备,所述计算机设备包括处理器和存储器,所述存储器中存储有计算机程序,当所述计算机程序被所述处理器加载并执行如实施例1~实施例8任一项所述的方法。
实施例10:一种计算机可读存储介质,在可读存储介质中存储有计算机程序,所述计算机程序被处理器加载并执行如实施例1~实施例8任一项所述的方法。
实施例11:下面将结合图1~图3,对本实施例的具体实施方式作进一步的说明。本实施例提供了一种适用于不可压流动求解的远场边界条件多窗口处理方法。操作步骤如图1所示,具体包括如下步骤:
步骤1:建立带有“远场”边界的0012翼型几何模型,并划分网格,网格总量19065。
步骤2:将网格导入到SIMPLE求解器中。
步骤3:对网格进行前处理。
步骤4:输入计算参数,通过输入的参数对边界条件进行重定义,具体如下:
① 将边界条件分别定义到每个网格的网格单元,即多窗口化处理边界条件;
② 判断各个单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;
③ 若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”;
步骤5:使用SIMPLE算法对流场进行计算;
步骤6:对收敛的流场进行后处理显示。
图2为带有“远场”边界的0012翼型网格划分示意图。
图3为0012翼型模型计算结果压力场示意图。
从上述技术方案可以看出,本实施例提供的一种适用于不可压流动求解的远场边界条件多窗口处理方法,根据输入参数,将“远场”边界条件转化为传统的“入口”、“出口”、“对称”等边界条件,成功解决了不可压缩流动数值模拟中难以使用“远场”边界的难题,并准确计算出了流场结果。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围内。因此,本发明的保护范围应该以权力要求的保护范围为准。
本说明书中未作详细描述的内容属于本领域普通技术人员公知的现有技术。
除以上实例以外,本领域技术人员根据上述公开内容获得启示或利用相关领域的知识或技术进行改动获得其他实施例,各个实施例的特征可以互换或替换,本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。

Claims (6)

1.一种远场边界条件多窗口处理方法,其特征在于,包括步骤:
在求解器输入参数,通过输入的参数对每个网格分别进行边界条件的赋值实现边界条件的重定义,具体包括如下步骤:
S1,创建出一套流场计算的网格并定义其边界条件;
S2,将网格导入到求解器中;
S3,对网格进行前处理;
S4,输入参数,通过输入的参数对每个网格进行边界条件的赋值实现边界条件的重定义;
S5,使用不可压流动求解算法对流场进行计算;
其中,在步骤S4中,包括子步骤:
S41,将边界条件分别定义到每个网格单元,即多窗口化处理边界条件;
S42,判断各个网格单元的边界条件类型,若为“远场”边界条件,则对该网格边界处的流量进行计算;
S43,若流量为正,重新定义该网格处边界为“入口”;若流量为负,重新定义该网格处边界为“出口”;若流量为零,重新定义该网格处边界为“对称”;
再使用不可压流动求解算法对流场进行计算。
2.根据权利要求1所述的远场边界条件多窗口处理方法,其特征在于,所述网格为任意形状的边界网格。
3.根据权利要求1所述的远场边界条件多窗口处理方法,其特征在于,所述不可压流动求解算法包括SIMPLE算法、SIMPLEC算法、SIMPLEX算法、SIMPLER算法、PISO算法中的任一种。
4.根据权利要求1所述的远场边界条件多窗口处理方法,其特征在于,在步骤S5之后,包括步骤:S6,对收敛的流场进行后处理显示。
5.一种计算机设备,其特征在于,所述计算机设备包括处理器和存储器,所述存储器中存储有计算机程序,当所述计算机程序被所述处理器加载并执行如权利要求1~4任一项所述的方法。
6.一种计算机可读存储介质,其特征在于,在可读存储介质中存储有计算机程序,所述计算机程序被处理器加载并执行如权利要求1~4任一项所述的方法。
CN202210809742.3A 2022-07-11 2022-07-11 一种远场边界条件多窗口处理方法、设备及介质 Active CN114880890B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210809742.3A CN114880890B (zh) 2022-07-11 2022-07-11 一种远场边界条件多窗口处理方法、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210809742.3A CN114880890B (zh) 2022-07-11 2022-07-11 一种远场边界条件多窗口处理方法、设备及介质

Publications (2)

Publication Number Publication Date
CN114880890A CN114880890A (zh) 2022-08-09
CN114880890B true CN114880890B (zh) 2022-10-14

Family

ID=82682856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210809742.3A Active CN114880890B (zh) 2022-07-11 2022-07-11 一种远场边界条件多窗口处理方法、设备及介质

Country Status (1)

Country Link
CN (1) CN114880890B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115859480B (zh) * 2023-02-08 2023-05-02 中国空气动力研究与发展中心计算空气动力研究所 基于确定发动机入口边界条件的气动分析方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105183975A (zh) * 2015-09-01 2015-12-23 中国空气动力研究与发展中心高速空气动力研究所 一种基于尾涡流场传递的多飞行器飞行编队数值模拟方法
CN105631068A (zh) * 2014-10-31 2016-06-01 北京临近空间飞行器系统工程研究所 一种非结构网格cfd计算的网格边界条件处理方法
CN112182771A (zh) * 2020-10-11 2021-01-05 中国运载火箭技术研究院 基于数值模拟的数据处理方法、存储介质、电子装置
CN112989680A (zh) * 2021-05-14 2021-06-18 中国空气动力研究与发展中心计算空气动力研究所 减少网格使用量的fvfd远场积分边界条件计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210357555A1 (en) * 2018-09-14 2021-11-18 Northwestern University Data-driven representation and clustering discretization method and system for design optimization and/or performance prediction of material systems and applications of same
US11580282B2 (en) * 2020-03-23 2023-02-14 Hasl Co., Ltd. Multilayer fluid analysis program, and multilayer fluid analysis system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631068A (zh) * 2014-10-31 2016-06-01 北京临近空间飞行器系统工程研究所 一种非结构网格cfd计算的网格边界条件处理方法
CN105183975A (zh) * 2015-09-01 2015-12-23 中国空气动力研究与发展中心高速空气动力研究所 一种基于尾涡流场传递的多飞行器飞行编队数值模拟方法
CN112182771A (zh) * 2020-10-11 2021-01-05 中国运载火箭技术研究院 基于数值模拟的数据处理方法、存储介质、电子装置
CN112989680A (zh) * 2021-05-14 2021-06-18 中国空气动力研究与发展中心计算空气动力研究所 减少网格使用量的fvfd远场积分边界条件计算方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FVM与LBM分区耦合模拟远场边界条件下的圆柱绕流;栾辉宝 等;《中国计算力学大会"2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集》;20100820;第1-7页 *
Time to redefine PD? Introductory statement of the MDS Task Force on the definition of Parkinson"s disease. Movement Disorders;Daniela Berg 等;《Featured Article》;20140311;第29卷(第4期);第454-462页 *
一种新的定压边界处理方法;王辉 等;《中国科学院研究生院学报》;20111115;第28卷(第6期);第740-745页 *
有限基本解法中不可压源的远场近似;范家义;《南京航空航天大学学报》;19780430;第15-22页 *
计算流体力学中变形运动边界处理方法及应用;敬军 等;《中国舰船研究》;20070415;第2卷;第57-62页 *

Also Published As

Publication number Publication date
CN114880890A (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
CN108038040B (zh) 计算机集群性能指标检测方法、电子设备及存储介质
Haelterman et al. Improving the performance of the partitioned QN-ILS procedure for fluid–structure interaction problems: Filtering
CN114880890B (zh) 一种远场边界条件多窗口处理方法、设备及介质
ur Rehman et al. A quadrature method for numerical solutions of fractional differential equations
CN112965903A (zh) 测试方法、装置、电子设备及计算机可读存储介质
Shete et al. Area of scalar isosurfaces in homogeneous isotropic turbulence as a function of Reynolds and Schmidt numbers
Tseluiko et al. A global attracting set for nonlocal Kuramoto–Sivashinsky equations arising in interfacial electrohydrodynamics
Cheng et al. IB-WENO method for incompressible flow with elastic boundaries
CN115658021A (zh) 动态模型的确定方法、装置、存储介质及电子设备
CN114595646A (zh) 一种基于计算流体动力学的分析方法
US20150356132A1 (en) Methods And Systems For Processing Data
CN112732685A (zh) 金融数据处理方法、装置、计算机设备及存储介质
Peng et al. Estimating distribution sensitivity using generalized likelihood ratio method
CN113553407B (zh) 事件追溯方法、装置、电子设备及存储介质
Konovalov et al. Developing a Technique for Constructing Surrogate Models of Equipment on the Example of a Hydraulic Model of a Gas-Removal Device
Hernández-Maldonado et al. Fast procedure to compute empirical and Bernstein copulas
CN116187034B (zh) 基于不确定度量化的压气机仿真可信度评估方法
CN114676656B (zh) 多响应cfd模型的一致性度量方法、装置、设备及存储介质
US20230142773A1 (en) Method and system for real-time simulations using convergence stopping criterion
CN116150438B (zh) 基于交易图谱的数据处理方法及装置
CN117375577B (zh) 声传播问题的数值滤波方法、装置、电子设备及存储介质
CN114595231B (zh) 一种数据库表生成方法、装置、电子设备及存储介质
CN115495854B (zh) 计算机辅助工程模型的参数标定方法、装置、设备及介质
JP2015049724A (ja) 連立一次方程式の計算方法及び数値計算装置
CN115455863A (zh) 一种流体模型的三维稳定性分析方法及装置、电子设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant