一种气动流动数值仿真流场结构的显示方法技术

本发明专利技术涉及一种气动流动数值仿真流场结构的显示方法，包括：S1.获取气动流动仿真的流场结果数据，并基于初始物理判据对所述流场结果数据进行筛选，获取所述流场结果数据中存在激波面的第一位置；S2.基于筛选物理判据对所述第一位置进行筛选处理，获取激波面结构；S3.汇总所述激波面结构形成气动全流场激波波系结构，并显示所述气动全流场激波波系结构。通过采用本发明专利技术的方法，对流场结构的显示无需显卡支持，可直接在计算集群上对流场仿真数据进行流动结构的批量提取，提取结果文件大小比流场仿真数据约小2个数量级，可以方便地在个人电脑上查看流场结构。提取结果可清晰直观地显示出流场的波系结构，对飞行器气动设计和分析提供了极大的便利。

【技术实现步骤摘要】
一种气动流动数值仿真流场结构的显示方法
本专利技术涉及一种气动流动数值仿真流场结构的显示方法，尤其涉及一种基于物理特性的气动流动数值仿真流场结构的显示方法。
技术介绍
基于计算能力的迅速提升，CFD气动仿真技术近年来得到了蓬勃发展，得到了大规模的工程应用。针对复杂外形的跨声速/超声速流动仿真已逐渐替代风洞试验成为气动设计与分析环节中常用的设计分析手段。对于航空航天常见的复杂外形超声速流动中，常常存在复杂的激波波系。激波在流场中发生反射、干扰，形成非常复杂的波系结构，对飞行器的气动力造成影响，同时飞行速度高时，波系干扰会导致严重的气动加热，影响飞行器热防护设计。因此准确识别激波波系有助于理解复杂流场，有助于飞行器防热设计，在此基础上方能有效地开展气动设计与分析工作。目前常见的CFD求解器主要分为基于结构网格的结构求解器和基于非结构网格的非结构求解器两类。对于复杂外形流动，工程实用的计算网格量一般在千万量级，单个流场计算结果文件大小为1～10G量级。在普通电脑上难以直接对完整的流场文件进行后处理和可视化，一般需要采用专业显卡的高配置工作站。
技术实现思路
本专利技术的目的在于提供一种气动流动本文档来自技高网...

【技术保护点】
1.一种气动流动数值仿真流场结构的显示方法，包括：S1.获取气动流动仿真的流场结果数据，并基于初始物理判据对所述流场结果数据进行筛选，获取所述流场结果数据中存在激波面的第一位置；S2.基于筛选物理判据对所述第一位置进行筛选处理，获取激波面结构；S3.汇总所述激波面结构形成气动全流场激波波系结构，并显示所述气动全流场激波波系结构。

【技术特征摘要】
1.一种气动流动数值仿真流场结构的显示方法，包括：S1.获取气动流动仿真的流场结果数据，并基于初始物理判据对所述流场结果数据进行筛选，获取所述流场结果数据中存在激波面的第一位置；S2.基于筛选物理判据对所述第一位置进行筛选处理，获取激波面结构；S3.汇总所述激波面结构形成气动全流场激波波系结构，并显示所述气动全流场激波波系结构。2.根据权利要求1所述的显示方法，其特征在于，步骤S1中包括：S11.读取气动流动仿真获得的流场结果数据；S12.选取流场数据中一个网格单元的单元中心点PE，i为起始，以流线推进方式分别沿上游、下游获取其余所述流场单元的单元中心点PE，n，以及相应的流场数据；S13.根据所述初始物理判据对所述流场数据进行判断，获取所述第一位置。3.根据权利要求2所述的显示方法，其特征在于，步骤S12中，所述流线推进方式满足：Pk+1＝Pk+sdVkλk其中，Pk为当前点，Pk+1为流线上的下一个点，sd取值为±1，代表是向下游搜索还是向上游搜索，Vk是Pk点处的流动速度矢量，λk是Pk点处的网格最小尺度。4.根据权利要求3所述的显示方法，其特征在于，步骤S12中所述流场数据包括：流动速度V、压力p、密度ρ和温度T。5.根据权利要求1或4所述的显示方法，其特征在于，所述初始物理判据包括第一子物理判据和第二子物理判据；其中，所述第一子物理判据为：其中，MaE,i为PE,i点马赫数；所述第二子物理判据为：PE,-N-PE,N≥-λPE,0TE,-N-TE,N≥-λTE,0ρE,-N-ρE,N≥-λρE,0|VE,...

【专利技术属性】
技术研发人员：徐艺哲，左光，张柏楠，杜若凡，周印佳，万千，侯砚泽，
申请(专利权)人：北京空间技术研制试验中心，
类型：发明
国别省市：北京,11