【技术实现步骤摘要】
一种基于变光滑长度的SPH流体模拟方法
本专利技术属于计算机图形学领域，具体涉及一种基于变光滑长度的SPH流体模拟方法。
技术介绍
流体遇到冲击会产生很大的变形，针对模拟大变形流体问题，SPH方法不仅可以处理大变形问题而且能够对粒子进行效果很好的追踪，非常适合对大变形流体的模拟。在SPH方法中，使用粒子近似法主要依赖于一个以光滑半径为长度构造出的粒子支持域，计算效率和精度取决于邻近粒子的数量。光滑长度的选取决定局部核函数近似的精度，直接影响流体模拟的计算速度和真实性。若光滑长度过小，则在支持域中没有足够的粒子对给定粒子提供插值信息和施加作用力，导致模拟结果的精度降低。若光滑长度过大，则用于更新给定粒子的邻近粒子部分信息或局部特性将会丢失，也会导致模拟结果的精度降低。传统SPH方法从光滑长度与粒子密度关系入手，将光滑长度看成受恒定临近粒子数约束条件的独立坐标变量，提出全守恒SPH方程组。这种方法存在方程组复杂、数值误差较大以及适用范围窄的问题，不具有通用性和实用性。基于对前人的研究，本专利技术提出一种从对称核近似角度出...
【技术保护点】
1.一种基于变光滑长度的SPH流体模拟方法，其特征在于包括如下步骤：/n步骤一，从对称核近似角度出发修正SPH密度演化方程、动量方程和能量方程，建立一组修正后的变光滑长度SPH方程组；/n步骤二，由建立的SPH方程组计算出粒子所受的压力、重力、黏性力，再根据牛顿第二定律计算出合力得到加速度；/n步骤三，更新下一步粒子的速度和位置；/n步骤四，基于Marching Cubes算法实现流体表面提取，利用GPU纹理缓存的位置、速度和颜色数据进行三维可视化渲染；/nMarching Cubes流体表面提取算法涉及网格节点处的密度计算、交点坐标和法矢量计算。采用CUDA并行计算网格节...
【技术特征摘要】
1.一种基于变光滑长度的SPH流体模拟方法，其特征在于包括如下步骤：
步骤一，从对称核近似角度出发修正SPH密度演化方程、动量方程和能量方程，建立一组修正后的变光滑长度SPH方程组；
步骤二，由建立的SPH方程组计算出粒子所受的压力、重力、黏性力，再根据牛顿第二定律计算出合力得到加速度；
步骤三，更新下一步粒子的速度和位置；
步骤四，基于MarchingCubes算法实现流体表面提取，利用GPU纹理缓存的位置、速度和颜色数据进行三维可视化渲染；
MarchingCubes流体表面提取算法涉及网格节点处的密度计算、交点坐标和法矢量计算。采用CUDA并行计算网格节点处的密度值和粒子网格节点密度值，计算思想由影响域改为支持域，即当前网格节点的密度值由邻域范围内所有粒子计算所得，而不是计算一个粒子点对邻域范围内产生的影响。
2.根据权利要求1所述的一种基于变光滑长度的SPH流体模拟方法，其特征在于：所述步骤一具体修正过程如下：
(1)密度方程
传统SPH法中，某点ri处密度ρ(ri)的核函数近似式为
其中，W(ri-rj,h)称为核函数，h为光滑长度描述核近似的程度，N为i粒子邻近粒子总数。
光滑长度h的选取，应尽量使得粒子的邻近数目在计算中保持不变，以获得全场一致的核近似精度，通常设定粒子i的光滑长度hi半径区域内包含恒定质量Msph，即：
其中，d为维数，σ对应一维到三维分别取2、π、公式(2)的引入使得光滑长度在空间分布不一致，并随同粒子的运动变化，即光滑长度同为空间和时间的函数hi＝h(ri,t)。
对于变光滑长度核函数近似，采用对称核近似的方法可以得到更高精度的近似值，本发明采用对称光滑长度的方法来对对称函数近似，此时密度核函数近似式为：
其中，
公式(3)两边分别对时间求导数得到：
公式(4)即为用对称核近似的变光滑长度SPH密度演化方程，其中νij＝νi-νj，在恒定光滑长度时该方程退化为传统SPH密...
【专利技术属性】
技术研发人员：倪维涛，孔凡玉，沈红威，
申请(专利权)人：中国计量大学，
类型：发明
国别省市：浙江;33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。