一种基于径向基函数的DSMC计算结果降噪处理方法技术

本发明专利技术公开了一种基于径向基函数的DSMC计算结果降噪处理方法。本发明专利技术借鉴碰撞网格与抽样网格分开的思想，采用两套不同粗细的两套网格，增大统计样本数来降低统计噪声。利用径向基函数工具，完成原有碰撞网格节点上的物理数据向抽样网格节点上的数据插值，无需增大单元内模拟分子数目或者增大抽样步数，具有不影响主计算程序，计算简单，工作量小和计算精度较高的优点。较高的优点。较高的优点。

【技术实现步骤摘要】
一种基于径向基函数的DSMC计算结果降噪处理方法


[0001]本专利技术涉及稀薄气体动力学数值模拟领域，具体涉及一种基于径向基函数的DSMC计算结果降噪处理方法。

技术介绍

[0002]作为稀薄流动数值计算领域的主要工具，DSMC方法建立在概率统计的基础之上。该方法用大量的模拟分子代替真实的气体分子，通过模拟分子的运动及碰撞，从而实现对物理流动的模拟。宏观的物理信息需要通过对微观的模拟分子信息进行抽样统计而得到。因此，计算结果中存在一定的统计误差，或者称作统计“噪声”，会干扰流动结构并影响对计算结果的读取和分析。而且在DSMC/NS混合算法中，统计误差的存在会给计算的收敛带来困难。统计误差与模拟分子数的平方根成反比。为了降低统计误差，需要单元内的模拟分子数目足够多，或者在一定的时间段内，对单元内的模拟分子进行重复的抽样以增大样本数目。但这两种方式，无疑将增加计算成本，降低计算效率，特别是在Kn数较低的情况下。
[0003]为了保证碰撞计算的准确，碰撞网格的尺度需要与当地的分子平均自由程保持在同一量级。但是宏观量的抽样并没有特别的必要非得限制在碰撞网格内进行。将抽样网格与碰撞网格分开处理，抽样网格可以是碰撞网格的数倍，样本数大大增大，统计噪声从而大大降低。值得指出的是，这样的处理在直角和结构网格中易于操作，但是在非结构网格中却并不方便。

技术实现思路

[0004]针对现有技术中的上述不足，本专利技术提供的一种基于径向基函数的DSMC计算结果降噪处理方法解决了DSMC碰撞计算精度不高的问题。/>[0005]为了达到上述专利技术目的，本专利技术采用的技术方案为：一种基于径向基函数的DSMC计算结果降噪处理方法，包括以下步骤：
[0006]S1、根据计算条件和几何外形生成满足DSMC碰撞算法要求的标准计算网格，进行完整的DSMC计算，得到待处理结果；
[0007]S2、生成网格最小尺度为局部分子平均自由程N倍的壁面抽样网格；
[0008]S3、利用径向基函数，实现壁面物理量在碰撞网格和抽样网格之间的插值，完成待处理结果的降噪处理。
[0009]进一步地：所述步骤S1中的DSMC碰撞算法要求为网格最小尺度小于三分之一的局部分子平均自由程。
[0010]进一步地：所述步骤S2中N的取值为1或2。
[0011]进一步地：所述步骤S3的具体步骤为：
[0012]S31、在d维欧几里得空间给定一个点集X＝{x1,x2,
…
,x
N
}，以及这些点上的标量值g1,g2,
…
g
N
，并定义插值函数g(x)：
[0013][0014]上式中，g(x)为x点处的未知函数值，α
i
为插值系数，x为未知点的坐标，x
i
为已知点的坐标，φ为基函数，||x
‑
x
i
||为两点间的欧几里得距离，p(x)为低阶d维多项式；
[0015]S32、根据插值函数利用定解条件求出未知系数，定解条件为：
[0016]g(x
i
)＝g
i
,i＝1,2,
…
,N
[0017][0018]上式中，q(x)为满足deg(q(x))≤deg(p(x))的所有多项式；
[0019]S33、根据未知系数和抽样网格上的物理量计算传递矩阵；
[0020]S34、通过传递矩阵计算壁面物理量在碰撞网格和抽样网格之间的插值。
[0021]进一步地：所述低阶d维多项式p(x)为：
[0022]p(x)＝γ0+γ1x+γ2y+γ3z
[0023]上式中，γ0、γ1、γ2和γ3分别为多项式系数。
[0024]进一步地：所述传递矩阵的计算公式为：
[0025][0026]上式中，H为传递矩阵，A
sc
、P和M为中间参数；
[0027]其中，A
sc
的计算公式为：
[0028][0029]上式中，为节点N
s
的坐标，为节点N
s
和节点N
c
之间的距离函数，N
s
为抽样网格的节点个数，N
c
为碰撞网格的节点个数；
[0030]P和M的计算公式为：
[0031][0032]上式中，D
cc
为中间参数；
[0033]D
cc
的计算公式为：
[0034][0035]上式中，
[0036]本专利技术的有益效果为：本专利技术借鉴碰撞网格与抽样网格分开的思想，采用两套不同粗细的两套网格，增大统计样本数来降低统计噪声。利用径向基函数工具，完成原有碰撞网格节点上的物理数据向抽样网格节点上的数据插值，无需增大单元内模拟分子数目或者增大抽样步数，具有不影响主计算程序，计算简单，工作量小和计算精度较高的优点。
[0037]本专利技术对DSMC计算过程没有任何影响，通过后处理方式对DSMC壁面结果进行光滑性的改进，独立性强，无需对主程序进行改动，理论上适应任意的非结构四面体DSMC解算结果。
[0038]经过插值处理后，计算结果中的统计误差得到了有效的抑制，等值线变得很光滑。插值前后摩擦力系数最大误差控制在5％以内，压力系数极大值和热流系数极大值的最大误差分别控制在1％和2％以内，在可接受精度范围内。
[0039]本专利技术中具有较强的适应性，插值结果对径向基函数的具体形式和网格大小不敏感。在一定的精度许可范围内，可以用来作为降低DSMC计算结果中统计误差的一种初步处理的方式。
附图说明
[0040]图1为本专利技术流程图；
[0041]图2为标准网格及该网格下的壁面摩擦力等值线分布图；
[0042]图3为插值网格及该网格下基于不同径向基函数的壁面摩擦力等值线分布图。
具体实施方式
[0043]下面对本专利技术的具体实施方式进行描述，以便于本
的技术人员理解本专利技术，但应该清楚，本专利技术不限于具体实施方式的范围，对本
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本专利技术的精神和范围内，这些变化是显而易见的，一切利用本专利技术构思的专利技术创造均在保护之列。
[0044]DSMC计算中，壁面抽样统计的值经常出现较大的统计误差。为了保证计算的精度，碰撞网格必须采用比较小的尺度，而用于抽样的网格应该比碰撞网格的尺度要大，这样才能有效地降低统计的噪声。
[0045]本专利技术针对DSMC方法在固壁面抽样的数据统计散落较大的问题，借鉴碰撞网格与抽样网格分开的思想，采用两套网格，并利用径向基函数插值完成了壁面物理量在两套网格之间的映射。初步的研究结果表明，这种方式可以在一定程度上降低DSMC计算结果中的统计误差，其本质上可以看作是一种滤波或者平均。
[0046]如图1所示，一种基于径向基函数的DSMC计算结果降噪处理方法，包本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于径向基函数的DSMC计算结果降噪处理方法，其特征在于，包括以下步骤：S1、根据计算条件和几何外形生成满足DSMC碰撞算法要求的标准计算网格，进行完整的DSMC计算，得到待处理结果；S2、生成网格最小尺度为局部分子平均自由程N倍的壁面抽样网格；S3、利用径向基函数，实现壁面物理量在碰撞网格和抽样网格之间的插值，完成待处理结果的降噪处理。2.根据权利要求1所述的基于径向基函数的DSMC计算结果降噪处理方法，其特征在于，所述步骤S1中的DSMC碰撞算法要求为网格最小尺度小于三分之一的局部分子平均自由程。3.根据权利要求1所述的基于径向基函数的DSMC计算结果降噪处理方法，其特征在于，所述步骤S2中N的取值为1或2。4.根据权利要求1所述的基于径向基函数的DSMC计算结果降噪处理方法，其特征在于，所述步骤S3的具体步骤为：S31、在d维欧几里得空间给定一个点集X＝{x1,x2,
…
,x
N
}，以及这些点上的标量值g1,g2,
…
g
N
，并定义插值函数g(x)：上式中，g(x)为x点处的未知函数值，α
i
为插值系数，x为未知点的坐标，x
i
为已知点的坐标，φ为基函数，||x
‑
x
i
||为两点间的欧几里得距离，p(...

【专利技术属性】
技术研发人员：李锦，耿湘人，江定武，陈坚强，
申请(专利权)人：中国空气动力研究与发展中心计算空气动力研究所，
类型：发明
国别省市：