CN113158531A - 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 - Google Patents
一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 Download PDFInfo
- Publication number
- CN113158531A CN113158531A CN202110174788.8A CN202110174788A CN113158531A CN 113158531 A CN113158531 A CN 113158531A CN 202110174788 A CN202110174788 A CN 202110174788A CN 113158531 A CN113158531 A CN 113158531A
- Authority
- CN
- China
- Prior art keywords
- component
- fluid
- simulation
- particle
- equation
- 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.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 139
- 238000004088 simulation Methods 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 69
- 239000002245 particle Substances 0.000 claims abstract description 78
- 230000008859 change Effects 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 239000007787 solid Substances 0.000 claims description 9
- 230000004048 modification Effects 0.000 claims 1
- 238000012986 modification Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 22
- 239000000203 mixture Substances 0.000 description 10
- 238000012733 comparative method Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000006835 compression Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 244000304337 Cuminum cyminum Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 229910052738 indium Inorganic materials 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种利用形变梯度的单组分与多组分不可压缩流体仿真方法,利用形变梯度,计算单组分或多组分粒子的体积变化量,并基于形变梯度建立流体不可压缩条件的方程并用松弛的雅各比迭代方法求解。本发明方法实现了对多组分流体的不可压缩求解,实现了真实的多组分流体仿真效果,并提供了艺术化控制的方案,对单组分流体的仿真效率和效果与当前业界先进方法的效果可比。
Description
技术领域
本发明涉及计算机图形学物理模拟领域,尤其涉及一种单组分与多组分不可压缩流体 的动力学仿真方法与***。
背景技术
不可压缩的流体仿真是流体动力学的研究重点,它对于流体仿真的视觉真实性,物理 可靠性都有重要影响。多组分的流体仿真,是流体仿真中重要的研究领域。
对于不可压缩性的流体仿真,前人有着长时间的研究。借鉴于基于网格的不可压缩流 体求解器,早期方法,如S.J.Cummins and M.Rudman在“An sph projection
Method”一文中,通过直接将光滑粒子映射到离散的网格上,并求解不可压缩条件, 以维持流体的不可压缩性。但是这种方法对采样方式敏感,而且尺度伸缩性差。近年来, M.Ihmsen等人在“Implicit incompressible sph”提出的隐式不可压缩的sph方法,改进了离散化方式,使得整个方案在并行图形处理单元上可以运行。
在多组分的流体仿真领域,前人也有出色的研究成果。如B.Ren等人在“Multiple-fluid sph simulation using a mixture model”提出的使用混合的SPH模 型对多组分流体进行仿真,实现了多组分流体的分层效果。但是由于混合流体的密度不是 常量,速度场也不是无散的,该方法难以保证流体的可压缩性。
发明内容
本发明提出了一种新颖的保证流体仿真中流体不可压缩性的算法,借助了形变梯度 (deformation gradient)来衡量物理量的体积变化。提出的该算法框架不仅可以实现传 统的单组分不可压缩流体仿真任务,还可以实现新的多组分不可压缩流体仿真任务。该算 法框架在实现了高真实感的流体仿真的同时,保证了流体运动的高效求解。
本发明的技术方案:
一种利用形变梯度的单组分与多组分不可压缩流体仿真方法,本发明方法可以进行下 述流体仿真任务:
1)对单组分流体进行基于光滑粒子动力学方法的不可压缩流体仿真;
2)对多组分流体进行基于光滑粒子动力学方法的不可压缩流体仿真;对于单组分或 多组分流体的仿真步骤如下:
第1步:初始化流体场景的粒子参数,包括位置,速度,各组分体积分数,密度;
第2步:在每个时间步开始,根据当前场景中的各粒子参数,预测计算下一时刻的粒 子位置、速度,并计算出与预测对应的形变梯度及其雅可比行列式;对于多组分流体仿真 任务,额外提前根据混合流模型的质量方程计算粒子的组分与密度变化及每组分的偏移速 度;
第3步:利用第2步的计算结果,通过使用松弛雅可比迭代方式求解一个线性方程组 的形式,求解获得一个当前时间步内的粒子压强修正值;所述求解线性方程组的步骤如下:
第3.1步:对于每个粒子,利用第2步中的预测计算结果各项,推导出该粒子压强修正值与邻域粒子压强修正值需满足的方程等式,联立所有粒子需满足的方程等式构成一个线性方程组;构建线性方程组的方法如下:
第3.1.1步:借助形变梯度及雅可比行列式计算各流体粒子在仿真过程中的体积变化;
第3.1.2步:对于固体边界处的仿真,在流体粒子邻域存在边界固体粒子时,对所述 的线性方程组进行修正以包含边界固体粒子的影响;该修正不改变方程的整体数学形式和 求解方法;
第3.2步:对于单组分或多组分流体模拟任务,所求解的线性方程组及其求解方法具 有统一的数学形式,二者的区别仅在于少量的具体项的计算方法不同;
第3.3步:从一个初始估计压强修正值开始,根据松弛雅可比方法不断迭代获得更好 的压强修正值,直到仿真中所有粒子的方程左右两边的误差的最大值小于给定最小量,或 迭代次数超过给定上限;
第4步:利用第3步获得的压强修正值修正预测计算中的粒子位置,速度,获得不可压缩流体的仿真结果;
第5步:重复第2-4步,获得一个时间序列中每时间步对应的不可压缩流体的仿真结 果。
本发明的优点和有益效果:
本发明的创新点包括:
1)基于形变梯度的衡量流体形变的技术
物质物理上的压缩性直接与物质的体积相关,这个体积变化可以用形变梯度张量来衡 量,它的Jacobian式衡量了体积的变化率。形变梯度的概念起源于固体动力学仿真,用 来描述运动过程中的不可压缩固体。本发明引入了形变梯度的概念,并提出了新颖的仿真 方案,来保证基于SPH方法的多流体仿真中的不可压缩性。
2)多流体混合物仿真的不可压缩增强技术
本发明提出了一种新颖的保证多流体混合物不可压缩性的技术。借助形变梯度,直接 测定流体混合物中,由于物理性和非物理性导致的流体体积变化,从而推导出关于使得流 体不可压缩的压强修正值的线性方程组。这项新技术可以被直接应用在现有单流体和多流 体混合物的仿真框架中,并实现真实的视觉效果以及高效的计算。与现有的多流体混合物 仿真框架相比,本发明能够获得不可压缩的仿真结果;对单流体仿真来说,本发明与现有 先进框架的计算性能相当。
综合而言,本发明从理论上,将流体运动中的流体密度,速度散度与非物理性的压缩 进行解耦。本发明大大地增强了多流体仿真中的视觉真实性,避免了前人的技术方案中流 体压缩的问题。本发明应用在多流体仿真中,在避免了原有求解器的复杂的同时,仿真出 更加逼真的流体场景;应用在单流体仿真过程中,本发明可以集成到现有的流体仿真算法 框架中,而无需引入过多的计算量。
附图说明
图1为本发明利用相场理论进行二维与三维晶体生长模拟的***方框图。
图2是二维5角对称情形不同对称破缺的模拟效果。
图3是三维结晶生长结果及其对应的对称函数的图示。
图4是艺术控制模拟效果。
图5是艺术控制模拟效果。
图6是使用本发明方法恢复的二维晶体三角网格输出示意图。
图7是本发明方法在曲面流型上的效果。
图8展示了本发明方法相较对比方法能更好地实现不可压缩性。
图9展示了本发明方法相较对比方法能降低粒子杂乱度,提高分层效果。
具体实施方式
1.初始化
我们的方法可以使用现有标准的粒子方法仿真方案的通用初始化方法。对每个粒子,存储它的位置,速度,各组分体积分数,密度。
2.单/多组分流体仿真方法中的预测计算
根据标准的单组分流体纳维-斯托克斯动量方程(参见下公式2)和多组分流体混合 流模型纳维-斯托克斯方程(参见下公式15)的显示求解方法,根据当前时间步开始时所有粒子上存储的参量,进行显式地计算,获得一个预测的流体粒子加速度,并进一步根据该加速度,使用标准流体仿真的平流算法获得预测的下一时间步该粒子的位置、速度,并计算出与预测对应的形变梯度及其雅可比行列式;对于多组分流体仿真任务,额外提前根据混合流模型的质量方程(公式14)计算粒子的组分与密度变化及每组分的偏移速度;
3.单/多组分流体仿真方法中所需求解的线性方程组的推导过程、形式及求解方法
先以单流体为例,推导并解释本发明方案的基本思想,然后再扩展到多流体仿真应用 上。
3.1.借助形变梯度衡量体积变化
以单流体成分的仿真为例,流体运动由质量方程和动量方程来描述:
其中,ρ,p,u,M分别表示密度,压力,速度和流体受到的非压力外力。通常情况 下,单流体的不可压缩性是通过恒定的密度来实现的,通过公式1可以推出流体无散速 度的条件,也就是基于网格的解法里,传统的求解无散条件的方法通过求解带 条件的动量方程。但是在基于SPH的方法里,由于粒子的位置在不停发生变化,所以一种 高效的方式是使用当前时间步已知的周围粒子的数值来求解不可压缩条件。
一般来说,在流体仿真中,视觉上的可压缩效果通常是由两个原因导致的,一个是由 于弱可压缩条件或者求解器中的数值误差;另一个是由于流体物理属性发生变化导致的体 积变化,例如多流体混合物的反应,溶解等。因此,多流体仿真模型的设计应该允许由于 物理现象导致的流体密度变化,速度散度的同时,消除非物理性原因导致的压缩效果。
本发明方案设计了一个新颖的保证不可压缩的SPH仿真技术,通过当前时间步的粒子 拓扑计算形变梯度,并用形变梯度来衡量流体体积上的变化。形变梯度是通过形变位置与 初始时间步位置的梯度计算的,公式如下:
形变梯度的雅各比式衡量了体积的变化,其中V和V0分别表示形变后和 原始的元素体积。|F|表示F矩阵的行列式。因此,对于恒定密度的流体来说,要保证不可压缩性就相当于在仿真过程中,设置:J=1。对于多流体的仿真来说,J指示了流体粒 子局部体积的变化量,直接衡量了局部粒子体积变化是否与物理现象导致的体积变化一致。
3.2.单组分流体线性方程组
用xadv,vadv和Fadv分别表示预测的粒子位置,速度及形变梯度,这三个物理量是通过已知的当前帧的物理量,向前平流一个时间步得到的预测量。然后我们计算一个使J满足不可压缩条件的压力校正项Δp,来更正预测量。根据公式2和SPH速度更新方案,Δp会带来的速度改变量up为:
其中,
在后续推导中,任一项都默认是使用t时刻的数值进行计算的,为了简便,下文中的 下标t将省略,下标t+1会明确写出。公式7右侧可以通过下式计算:
J=|Fadv|+d|Fadv|=|Fadv|+tr(Fadv*dFadv) 公式8
基于上述式子,我们可以进一步替换up中的Δp。首先,为了方便起见,可以重写下式:
将公式11代入到公式10中,可以得到如下关于Δp的线性方程:
使用松弛雅各比方法,如下式迭代地求解Δp:
其中ω是松弛因子。
3.3.多组分流体线性方程组
对于多组分流体,求解混合物的SPH模型需要求解下列公式:
其中,ρm是粒子总密度,αk是流体中第k组分的体积分数,um,umk分别是粒子速度 和第k组分的漂移速度,Mm表示导致粒子速度改变的力。组分速度为:uk=um+umk。
多组分流体中,每个粒子的组成处于不断变化之中,要为一个流体粒子定义J值并不 容易,取而代之的是,我们为每个组分定义Jk和粒子总密度ρm=∑kαk,t+1ρk,来推导不可压缩求解方法。首先,混合粒子的密度可以表示为:
其中Jk是第k组分的形变梯度张量的雅各比式。公式16的第二个等式是通过假设各 个组分都是均匀地分布在SPH粒子上得到的。粒子上的一个元体积中,每个组分的体积分 数的升高使得粒子中该组分的比例升高,反之亦然。这就相当于流体局部区域的组分体积 的升高或收缩。因此,
上式与单组分仿真的对应公式10相比具有完全一致的形式。进一步地推导得到公式 12变为:
公式13变为:
比较公式12与公式19,公式13与公式20可见,单组份流体与多组分流体的线性方程组与其求解方法具有统一的数学形式。
3.4.边界处理
本发明采用了与Versatile rigid-fluid coupling for incompressible sph类似的 流固耦合方案。边界粒子对流体粒子i施加的压力为:
与单组分流体模拟相比,通过用ψ中的ρmi=∑kαkiρki替换ρ0i,得到:
对于单流体求解器,将公式22代入到公式10可以得到:
4.修正预测计算,保证流体不可压缩性
在第3步中已经计算出Δp后,重新进行第2步的显式计算过程,但在此时用到粒子压强p的时候以p+Δp作为取值,即利用Δp修正原先的p值。这样获得的下一时间步 的位置、速度即可保证流体的不可压缩性。
图1第一行展示了本发明方案的效果,第二行展示的是Multiple-fluid sphsimulation using a mixture model中方法的结果。展示的场景是三种不可互溶的液体的溃堤场景。可以从黄线的指示看出,本发明的方法更好地保证了流体体积上的不可压缩性。
图2展示仿真过程中测量到的平均密度误差。本发明方案显著地减少了密度误差,并 实现了更好的视觉效果。本发明方案的仿真效果,使得流体体积得到更好的维持,流体运 动过程中的外观更有序。
图3设置了一个立方体形状的流体掉入一个由两种流体构成的水池中,四种组分的液 体可互溶。第一行展示了本发明方案的效果,本发明方案实现了具有真实感的不可压缩流 体仿真。第二行展示了对比方法的结果,可以观察到,流体压缩性导致了粒子边界的剧烈 动荡,此外,可压缩性和剧烈的粒子动荡导致了混合后不自然的颜色。第三行展示了使用 对比方法但使用了1/100的时间步大小和100倍的stiffness。可以观察到,通过调小时间 步,流体整体体积得到了维持,但是更大的stiffness导致了边界动荡问题更加严重,甚至 会使仿真过程崩溃,粒子边界动荡同样地导致了不同流体组分的不自然的混合。在实验统 计中,第二行的仿真过程中,有48.7%-41.2%的体积损失,第三行的仿真过程中,通过使 用更小的时间步,体积损失减少到21.9%-9.1%,但是第三行的结果仍然存在着明显的不 自然的痕迹:剧烈的粒子间动荡导致了不同组分间流体的不自然的混合。本发明方案提供 了自然的不可压缩仿真方案。
图4展示了一个离心力场景。图中的圆柱形容器是完全密封并且完全注满流体的,也 就是说,如果流体的不可压缩性能被完美维持,那么容器中不会有任何额外空间存在。图 中第一行展示的是本发明方案效果,第二行展示的是对比方案的效果。可见,本发明方案 的仿真实现了高度的不可压缩性。而对比方案的仿真结果中,由于流体中出现了体积压缩, 导致容器中心出现了明显空缺。
图5的例子中,通过改变形变梯度来进行艺术化效果的控制。也就是说,在不改变组 分密度的前提下,通过改变J,粒子体积会相应发生改变。并将下一个时间步的粒子体积设置为期望值,通过重写公式16。替 换Jm,t为JcontJm,t,并将公式18中的ρmi,t+1替换为在仿真过程中,将Jcont设置为小 于1的数值会使相应位置的压力值更大,反之亦然。所以通过在场景中S型区域,在仿真 过程中途,人为设置Jcont=2,密度更大的深色组分向上运动,形成一个S型的喷泉,而其 它流体区域仍维持平稳。
图6,展示了本发明方法以及IISPH方法对单组分流体仿真的效果。为了比较两种方 法的效果,我们设置了一样的仿真参数,并设置了固定的迭代步数。在图6场景中,使用了52000个流体粒子,在一个Intel Xeon 2.6GHz CPU的单核上,本发明方案每个时间步 使用1.73s,IISPH方法耗时1.53s。仿真过程中的平均密度误差和最大密度误差展示在图 7中。本发明方案与IISPH的收敛速度相近。
图8展示了与对比方案的效果。第一行,设置了更高的stiffness系数,使用了更小的时间步,对比方案的流体仿真仍然开始有出现崩溃的倾向。第二行,使用对比方法并调小时间步为原来1/4,并调大stiffness系数为原来14倍,但是容器中间由于流体压缩 仍出现明显空缺。
图9中左侧展示了本发明方案的效果,左侧展示了对比方法的效果。从粒子排布上可 以观察到,本发明方案的粒子排布更加规整有序。而右侧中,相邻的粒子之间的颜色差异 很大,流体静止后的分层间也存在着无序粒子。本发明方法在仿真过程中,粒子处于有序 状态,流体静止后,分层效果更好。
Claims (3)
1.一种利用形变梯度的单组分与多组分不可压缩流体仿真方法,其特征在于,该仿真方法能够进行下述流体仿真任务:1):对单组分流体进行基于光滑粒子动力学方法的不可压缩流体仿真;2):对多组分流体进行基于光滑粒子动力学方法的不可压缩流体仿真;对于单组分或多组分流体的仿真步骤如下:
第1步:初始化流体场景的粒子参数,包括位置,速度,各组分体积分数,密度;
第2步:在每个时间步开始,根据当前场景中的各粒子参数,预测计算下一时刻的粒子位置、速度,并计算出与预测对应的形变梯度及其雅可比行列式;对于多组分流体仿真任务,额外提前根据混合流模型的质量方程计算粒子的组分与密度变化及每组分的偏移速度;
第3步:利用第2步的计算结果,通过使用松弛雅可比迭代方式求解一个线性方程组的形式,求解获得一个当前时间步内的粒子压强修正值;
第4步:利用第3步获得的压强修正值修正预测计算中的粒子位置,速度,获得不可压缩流体的仿真结果;
第5步:重复第2-4步,获得一个时间序列中每时间步对应的不可压缩流体的仿真结果。
2.根据权利要求1所述利用形变梯度的单组分与多组分不可压缩流体仿真方法,其特征在于,第3步所述求解线性方程组的步骤如下:
第3.1:对于每个粒子,利用第2步中的预测计算结果各项,推导出该粒子压强修正值与邻域粒子压强修正值需满足的方程等式,联立所有粒子需满足的方程等式构成一个线性方程组;
第3.2:对于单组分或多组分流体模拟任务,所求解的线性方程组及其求解方法具有统一的数学形式;
第3.3:从一个初始估计压强修正值开始,根据松弛雅可比方法不断迭代获得更好的压强修正值,直到仿真中所有粒子的方程左右两边的误差的最大值小于给定最小量,
或迭代次数超过给定上限。
3.根据权利2所述利用形变梯度的单组分与多组分不可压缩流体仿真方法中,其特征在于:第3.1步的构建线性方程组的方法如下,
第3.1.1步:借助形变梯度及雅可比行列式计算各流体粒子在仿真过程中的体积变化;
第3.1.2步:对于固体边界处的仿真,在流体粒子邻域存在边界固体粒子时,对权利要求2中的线性方程组进行修正以包含边界固体粒子的影响;该修正不改变方程的整体数学形式和求解方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110174788.8A CN113158531B (zh) | 2021-02-07 | 2021-02-07 | 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110174788.8A CN113158531B (zh) | 2021-02-07 | 2021-02-07 | 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113158531A true CN113158531A (zh) | 2021-07-23 |
CN113158531B CN113158531B (zh) | 2022-06-21 |
Family
ID=76883048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110174788.8A Active CN113158531B (zh) | 2021-02-07 | 2021-02-07 | 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113158531B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116797673A (zh) * | 2023-08-22 | 2023-09-22 | 四川沃耐稀新材料科技有限公司 | 一种硝酸锆冷却结晶状态监测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120123754A1 (en) * | 2009-06-03 | 2012-05-17 | Algoryx Simulation Ab | Method, an apparatus and computer program product for simulating dynamic fluids |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真*** |
CN104200015A (zh) * | 2014-08-20 | 2014-12-10 | 清华大学 | 一种流体模拟方法及装置 |
CN105389839A (zh) * | 2015-11-06 | 2016-03-09 | 北京航空航天大学 | 基于流体分析的流体参数估计方法 |
CN107045571A (zh) * | 2017-03-27 | 2017-08-15 | 杭州远算科技有限公司 | 一种基于云计算的压铸仿真控制方法 |
CN108763841A (zh) * | 2018-07-24 | 2018-11-06 | 北京航空航天大学青岛研究院 | 一种基于对偶边界元和应变能优化分析的弹性断裂仿真方法 |
CN109408836A (zh) * | 2017-08-16 | 2019-03-01 | 上海索辰信息科技有限公司 | 利用Boltzmann方程进行流体仿真的方法 |
CN111104753A (zh) * | 2019-12-30 | 2020-05-05 | 中国科学院软件研究所 | 一种基于sph的粘性不可压缩流体仿真方法 |
-
2021
- 2021-02-07 CN CN202110174788.8A patent/CN113158531B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120123754A1 (en) * | 2009-06-03 | 2012-05-17 | Algoryx Simulation Ab | Method, an apparatus and computer program product for simulating dynamic fluids |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真*** |
CN104200015A (zh) * | 2014-08-20 | 2014-12-10 | 清华大学 | 一种流体模拟方法及装置 |
CN105389839A (zh) * | 2015-11-06 | 2016-03-09 | 北京航空航天大学 | 基于流体分析的流体参数估计方法 |
CN107045571A (zh) * | 2017-03-27 | 2017-08-15 | 杭州远算科技有限公司 | 一种基于云计算的压铸仿真控制方法 |
CN109408836A (zh) * | 2017-08-16 | 2019-03-01 | 上海索辰信息科技有限公司 | 利用Boltzmann方程进行流体仿真的方法 |
CN108763841A (zh) * | 2018-07-24 | 2018-11-06 | 北京航空航天大学青岛研究院 | 一种基于对偶边界元和应变能优化分析的弹性断裂仿真方法 |
CN111104753A (zh) * | 2019-12-30 | 2020-05-05 | 中国科学院软件研究所 | 一种基于sph的粘性不可压缩流体仿真方法 |
Non-Patent Citations (2)
Title |
---|
BO REN等: "Multiple-Fluid SPH Simulation Using a Mixture Model", 《ACM TRANSACTIONS ON GRAPHICS》 * |
吴德阳等: "基于物理的不可压缩流体模拟技术综述", 《高技术通讯》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116797673A (zh) * | 2023-08-22 | 2023-09-22 | 四川沃耐稀新材料科技有限公司 | 一种硝酸锆冷却结晶状态监测方法 |
CN116797673B (zh) * | 2023-08-22 | 2023-10-31 | 四川沃耐稀新材料科技有限公司 | 一种硝酸锆冷却结晶状态监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113158531B (zh) | 2022-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eisenschmidt et al. | Direct numerical simulations for multiphase flows: An overview of the multiphase code FS3D | |
Müller et al. | Point based animation of elastic, plastic and melting objects | |
Yan et al. | Multiphase SPH simulation for interactive fluids and solids | |
Bender et al. | Position-based simulation of continuous materials | |
Liu et al. | Three-dimensional lattice Boltzmann model for immiscible two-phase flow simulations | |
Solenthaler et al. | A unified particle model for fluid–solid interactions | |
Bleyer | Advances in the simulation of viscoplastic fluid flows using interior-point methods | |
CN110298105A (zh) | 饱和多孔介质大变形分析的ccpdi-impm方法 | |
US10210287B2 (en) | Augmented material point method for simulating phase changes and varied materials | |
de Jesus et al. | A 3D front-tracking approach for simulation of a two-phase fluid with insoluble surfactant | |
Hua et al. | Level set, phase-field, and immersed boundary methods for two-phase fluid flows | |
CN108269299A (zh) | 一种基于sph方法近似求解的粘性流体建模方法 | |
CN108491619A (zh) | 基于物理与非物理混合的复杂场景流固耦合高效模拟方法 | |
Richefeu et al. | Lattice Boltzmann modelling of liquid distribution in unsaturated granular media | |
KR101113301B1 (ko) | 부피비를 이용한 다종 유체 해석 방법 및 기록 매체 | |
Shao et al. | Stable and fast fluid–solid coupling for incompressible SPH | |
Zhang et al. | A sharp interface method for SPH | |
CN113158531B (zh) | 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 | |
CN107423511A (zh) | 满足无滑边界条件和连续性条件浸入边界隐式迭代求解法 | |
Chen et al. | A multi-resolution SPH-FEM method for fluid–structure interactions | |
Fai et al. | Lubricated immersed boundary method in two dimensions | |
Zhang et al. | On large deformation granular strain measures for generating stress–strain relations based upon three-dimensional discrete element simulations | |
Caseiro et al. | A systematic development of EAS three-dimensional finite elements for the alleviation of locking phenomena | |
Joneidi et al. | Isogeometric boundary integral analysis of drops and inextensible membranes in isoviscous flow | |
CN113051842A (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 |