一种三维流场模拟方法技术

本发明专利技术提供了一种三维流场模拟方法

【技术实现步骤摘要】
一种三维流场模拟方法、系统及电子设备


[0001]本专利技术涉及环境模拟领域，特别是涉及一种基于地理环境与光滑粒子耦合的三维流场模拟方法
、
系统及电子设备
。

技术介绍

[0002]流场可视分析有利于地理环境专业知识科学表达与共享，帮助研究者对地理环境现状分布和时空变化规律有较为直观的理解，有效提高人们对时空规律的认知效能
。
近年来，粒子系统在复杂流场三维可视化过程中得到广泛研究，利用粒子离散化表达现实世界，在水动力矢量场可视化方面，国内外学者围绕流场数据动态模拟和仿真方面开展了大量研究，尤其在海洋水环境三维流场可视化方面开展了诸多研究和应用
。
[0003]但是，湖泊中小尺度流场可视分析缺乏研究
。
一方面，湖泊流场极易受到外部水文
、
气象环境等条件的影响，导致大气或水环境流场稳定性差
、
规律特征不明显，三维可视分析较为困难
。
另一方面，在中小尺度下，湖泊水环境流场可视分析对精度要求更高，如果仅仅考虑宏观欧拉网格速度场，将导致分析结果与实际情况产生偏差
。
因此，需要从微观角度进行三维流场物理过程仿真模拟和试验验证
。

技术实现思路

[0004]本专利技术的目的是提供一种三维流场模拟方法
、
系统及电子设备，可提高湖泊三维流场的模拟精度
。
[0005]为实现上述目的，本专利技术提供了如下方案：
[0006]一种三维流场模拟方法，包括：
[0007]获取待模拟湖泊的流域范围及地形；
[0008]基于所述待模拟湖泊的流域范围及地形，初始化三维欧拉网格，并在所述三维欧拉网格中初始化粒子数量
、
各粒子的位置及质量；所述三维欧拉网格包括多个子网格；
[0009]针对任一欧拉流场更新周期，获取所述欧拉流场更新周期内初始时刻的欧拉网格流场；
[0010]根据所述欧拉流场更新周期内初始时刻的欧拉网格流场，确定各子网格在初始时刻的速度及各粒子在初始时刻的速度；
[0011]将所述欧拉流场更新周期划分为多个子时段，并确定多个微观更新时刻；
[0012]针对任一微观更新时刻，基于所述三维欧拉网格
、
各粒子的质量
、
各粒子的当前位置及各粒子在前一微观更新时刻的速度，采用光滑粒子流体动力学方法，计算各粒子在所述微观更新时刻的受力情况；
[0013]根据各粒子在所述微观更新时刻的受力情况及各粒子在前一微观更新时刻的速度，更新各粒子的速度及位置，直至到达所述欧拉流场更新周期内的最后一个微观更新时刻；
[0014]获取所述欧拉流场更新周期内结束时刻的欧拉网格流场，并根据所述欧拉流场更
新周期内结束时刻的欧拉网格流场，更新各子网格的速度，以模拟粒子在水环境中的运动情况
。
[0015]为实现上述目的，本专利技术还提供了如下方案：
[0016]一种三维流场模拟系统，包括：
[0017]湖泊数据获取模块，用于获取待模拟湖泊的流域范围及地形；
[0018]网格粒子初始化模块，用于基于所述待模拟湖泊的流域范围及地形，初始化三维欧拉网格，并在所述三维欧拉网格中初始化粒子数量
、
各粒子的位置及质量；所述三维欧拉网格包括多个子网格；
[0019]流场获取模块，用于针对任一欧拉流场更新周期，获取所述欧拉流场更新周期内初始时刻的欧拉网格流场；
[0020]速度初始化模块，用于根据所述欧拉流场更新周期内初始时刻的欧拉网格流场，确定各子网格在初始时刻的速度及各粒子在初始时刻的速度；
[0021]时段划分模块，用于将所述欧拉流场更新周期划分为多个子时段，并确定多个微观更新时刻；
[0022]受力计算模块，用于针对任一微观更新时刻，基于所述三维欧拉网格
、
各粒子的质量
、
各粒子的当前位置及各粒子在前一微观更新时刻的速度，采用光滑粒子流体动力学方法，计算各粒子在所述微观更新时刻的受力情况；
[0023]微观更新模块，用于根据各粒子在所述微观更新时刻的受力情况及各粒子在前一微观更新时刻的速度，更新各粒子的速度及位置，直至到达所述欧拉流场更新周期内的最后一个微观更新时刻；
[0024]宏观更新模块，用于获取所述欧拉流场更新周期内结束时刻的欧拉网格流场，并根据所述欧拉流场更新周期内结束时刻的欧拉网格流场，更新各子网格的速度，以模拟粒子在水环境中的运动情况
。
[0025]为实现上述目的，本专利技术还提供了如下方案：
[0026]一种电子设备，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行上述的三维流场模拟方法
。
[0027]根据本专利技术提供的具体实施例，本专利技术公开了以下技术效果：
[0028]本专利技术结合欧拉与拉格朗日方法进行流场模拟，宏观上，通过三维欧拉流场更新各子网格的速度，并确定各粒子在初始时刻的速度，微观上，在欧拉流场更新周期内，利用光滑粒子流体动力学方法计算各粒子的受力情况，进而根据粒子的受力情况更新各粒子的速度及位置，实现水环境污染物扩散和沉积物集聚过程的模拟，提高了湖泊三维流场的模拟精度
。
附图说明
[0029]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图
。
[0030]图1为本专利技术提供的三维流场模拟方法的流程图；
[0031]图2为欧拉网格和光滑粒子耦合的流场三维分析过程示意图；
[0032]图3为本专利技术提供的三维流场模拟系统的示意图
。
[0033]符号说明：
[0034]1‑
湖泊数据获取模块，2‑
网格粒子初始化模块，3‑
流场获取模块，4‑
速度初始化模块，5‑
时段划分模块，6‑
受力计算模块，7‑
微观更新模块，8‑
宏观更新模块
。
具体实施方式
[0035]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚
、
完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例
。
基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围
。
[0036]本专利技术的目的是提供一种三维流场模拟方法
、
系统及电子设备，将地理环境与粒子系统耦合，通过三维网格数值模拟方式创建大型时空粒子系统，结合水动力模型
Delft3D...

【技术保护点】

【技术特征摘要】
1.
一种三维流场模拟方法，其特征在于，所述三维流场模拟方法包括：获取待模拟湖泊的流域范围及地形；基于所述待模拟湖泊的流域范围及地形，初始化三维欧拉网格，并在所述三维欧拉网格中初始化粒子数量
、
各粒子的位置及质量；所述三维欧拉网格包括多个子网格；针对任一欧拉流场更新周期，获取所述欧拉流场更新周期内初始时刻的欧拉网格流场；根据所述欧拉流场更新周期内初始时刻的欧拉网格流场，确定各子网格在初始时刻的速度及各粒子在初始时刻的速度；将所述欧拉流场更新周期划分为多个子时段，并确定多个微观更新时刻；针对任一微观更新时刻，基于所述三维欧拉网格
、
各粒子的质量
、
各粒子的当前位置及各粒子在前一微观更新时刻的速度，采用光滑粒子流体动力学方法，计算各粒子在所述微观更新时刻的受力情况；根据各粒子在所述微观更新时刻的受力情况及各粒子在前一微观更新时刻的速度，更新各粒子的速度及位置，直至到达所述欧拉流场更新周期内的最后一个微观更新时刻；获取所述欧拉流场更新周期内结束时刻的欧拉网格流场，并根据所述欧拉流场更新周期内结束时刻的欧拉网格流场，更新各子网格的速度，以模拟粒子在水环境中的运动情况
。2.
根据权利要求1所述的三维流场模拟方法，其特征在于，所述受力情况包括粒子的密度
、
粒子所受的粘性力及粒子所受的压力
。3.
根据权利要求2所述的三维流场模拟方法，其特征在于，采用以下公式计算粒子
i
的密度：其中，
ρ
i
为粒子
i
的密度，
m
i
为粒子
i
的质量，
h
为平滑半径，表示粒子
i
与粒子
j
之间的距离向量，根据粒子
i
与粒子
j
的当前位置确定
。4.
根据权利要求3所述的三维流场模拟方法，其特征在于，采用以下公式计算粒子
i
所受的粘性力：其中，
F
i,v
为粒子
i
所受的粘性力，
μ
为动力黏度，
m
j
为粒子
j
的质量，
v
i,t
‑1为粒子
i
在前一微观更新时刻的速度，
v
j,t
‑1为粒子
j
在前一微观更新时刻的速度，粒子
j
为粒子
i
的邻域粒子，
ρ
j
为粒子
i
的密度，
W()
为速度平滑核函数，为
Nabla
算子
。5.
根据权利要求4所述的三维流场模拟方法，其特征在于，采用以下公式计算粒子
i
所受的压力：其中，
F
i,p
为粒子
i
所受的压力
。6.
根据权利要求2所述的三维流场模拟方法，其特征在于，根据各粒子在所述微观更新
时刻的受力情况及各粒子在前一微观更新时刻的速度，更新各粒子的速度及位置，直至到达所述欧拉流场更新周期内的最后一个微观更新时刻，具体包括：针对任一粒子，根据所述粒子在所述微观更新时刻的受力情况，计算所述粒子的加速度；根据所述粒子的加速度及所述粒子在前一微观更新时刻的速度，计算所述粒子在所...

【专利技术属性】
技术研发人员：何望君，石丽红，张福浩，蔡鹏，廖振良，
申请(专利权)人：中国测绘科学研究院，
类型：发明
国别省市：