CN104143027B - 一种基于sph算法的流体热运动仿真系统 - Google Patents
一种基于sph算法的流体热运动仿真系统 Download PDFInfo
- Publication number
- CN104143027B CN104143027B CN201410377151.9A CN201410377151A CN104143027B CN 104143027 B CN104143027 B CN 104143027B CN 201410377151 A CN201410377151 A CN 201410377151A CN 104143027 B CN104143027 B CN 104143027B
- Authority
- CN
- China
- Prior art keywords
- particle
- sph
- fluid
- temperature
- boundary condition
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于SPH算法的流体热运动仿真系统,属于热力学模拟与可视化流动领域。该方法包括如下步骤:首先根据计算实例和计算要求对流体进行建模,然后基于SPH算法求解流体控制方程和修正方程,实现流体拉格朗日热运动的仿真模拟,最后应用可视化处理技术以图形及动画的方式直观的展示流体的热运动过程。与现有网格化CFD仿真系统相比,在模拟某些高度非线性问题(如射流问题、海洋运动问题、溃坝问题、动态边界问题)时,会得到更为合理的仿真结果;与现有无网格CFD仿真系统相比,提升了计算稳定性和计算精度,并且填补了流体热运动仿真和湍流仿真的空白。
Description
技术领域
本发明为一种基于光滑粒子流体动力学算法(SPH)的流体热运动仿真系统,属于热力学模拟与可视化流动领域。
背景技术
流体的热运动仿真在实际工程中有着广泛的应用背景,例如建筑物隔热、核反应堆安全设计、微电子设备的冷却等。
光滑粒子流体动力学方法是一种基于粒子的无网格数值模拟方法.该方法被应用于解决各种流体力学和固体力学问题,包括流固耦合问题、凝固问题、射流问题、多相流问题、热传导问题等.光滑粒子动力学方法作为最经典的无网格方法,其算法被不断完善,多种修正技术被提出并应用,并且发展了相应的前处理与后处理技术,其已经成为较为完善的粒子型无网格方法.
近年来,使用传统的网格化计算流体动力学(CFD)仿真系统在解决实际流体工程问题时遇到了许多难点,其模拟高度非线性问题(如射流问题、海洋运动问题、溃坝问题、动态边界问题)时往往会受到网格变形的制约,无法实现某些高度非线性问题的模拟,这些高度非线性问题往往需要通过采用无网格CFD方法来解决;而现有的无网格CFD仿真系统又无法解决热力学问题。
发明内容
本发明目的是为了解决现有网格化CFD仿真系统在解决某些高度非线性问题时因为网格的变形导致结果不够准确以及现有无网格CFD仿真系统无法对流体的热运动进行仿真的问题,提供一种基于SPH算法的流体热运动仿真系统。
本发明技术方案的思想是:根据计算实例和计算要求,首先对流体进行建模,然后基于SPH算法求解流体控制方程和修正方程,实现流体拉格朗日热运动的仿真模拟,最后应用可视化处理技术以图形及动画的方式直观的展示流体的热运动过程。
本发明是通过以下技术方案实现的:
一种基于SPH算法的流体热运动仿真系统,包括前处理模块、计算模块和后处理模块,整个仿真系统的模拟流程如图1所示。
1.前处理模块
所述前处理模块,用于根据实际模拟目标建立CFD仿真模型,一般CFD仿真模型需具备以下内容:离散方程、离散方法和离散模型,本仿真系统中离散方程为流体的流动控制方程(质量守恒方程、动量守恒方程、能量守恒方程),离散方法为SPH算法,离散模型为根据用户要求构建的粒子模型;因此,系统将前处理模块分为以下两个子模块:SPH算法设置和粒子模型建模;
1)SPH算法设置
前处理模块中的SPH算法设置将决定如何使用SPH算法,其包括以下四部分内容:计算方式设置、SPH算法格式设置、修正方法设置和SPH算法参数设置;
a)计算方式设置决定了该仿真系统的计算方式,以充分利用计算资源,其可选类型包括:中央处理器(CPU)并行、图形处理器(GPU)并行;其中CPU并行通过标准消息传递界面(MPI)技术实现,GPU并行通过通用并行计算架构(CUDA)实现。
b)若要使用SPH算法解算流体控制方程,则需要根据模拟对象来确定如何使用SPH算法,即需要设置SPH算法的格式;SPH算法的格式包括:时间积分格式、核函数格式、粘性解算格式;其中时间积分格式决定了模拟时间推进的方式,可选格式包括:蛙跳法、韦尔莱积分法、辛积分法;核函数为SPH算法的最基本函数,其决定了SPH算法如何离散流动控制方程,可选格式包括:三次样条核函数、五次样条核函数;粘性解算格式为动量守恒方程中流体粘性的计算方式,可选格式包括:层流粘性、人工粘性、亚粒子模型湍流粘性(湍流模型);
c)由于SPH算法易于产生震荡误差,则需要在算法中引入一些修正方法,修正方法包括:密度重构技术,人工位移技术(Xu R,Stansby P,Laurence D.Accuracy andstability in incompressible SPH(ISPH)based on the projection method and a newapproach[J].Journal of Computational Physics,2009,228(18):6703-6725);其中密度重构技术用于修正密度计算的震荡误差,可选类型包括:薛帕德滤波、δ-SPH;人工位移技术用于修正粒子位置的计算误差。
d)SPH算法参数包括:粒子间距离、时间步长、模拟时长、核函数的影响半径(光滑长度)、修正方程参数、人工声速(人工声速为人工状态方程的参数,人工状态方程用于压力求解)。
2)粒子模型建模
粒子模型建模模块用于对粒子模型进行建模,包括粒子坐标、粒子属性(流体或者固体)、边界条件和初始条件;
所述粒子坐标的生成通过输入相应外形参数实现,根据外形参数的几何特征及给定的粒子间距离,便可求得该几何形状内所有粒子的坐标,例如:按顺序输入矩形的四个顶点坐标,即可生成矩形粒子域;
所述粒子属性设置即是将边界粒子属性设置为固体,其他粒子属性设置为流体的过程;
所述边界条件确定了边界固体粒子对内部流体粒子的作用方式,可设置以下两方面:固体壁面边界条件、温度边界条件;
计算域边界部分粒子的属性为固体粒子,边界粒子位置不变,并对靠近的流体粒子产生一定作用力,作用力施加方式可选类型包括:边界力法、镜像粒子法;其中边界力法通过施加排斥力实现固体边界;镜像粒子法通过施加排斥动量实现固体边界;
当存在具有一定温度特性的边界时,该部分边界粒子需要设置温度边界条件,以使其温度特性影响其附近粒子,温度边界条件可选类型包括:恒温边界条件、绝热边界条件;恒温边界条件情况下,该部分固体粒子温度恒定,持续与周围粒子进行热交换;绝热边界条件情况下,该部分固体粒子与其他粒子不发生任何热交换。
所述初始条件为计算域内粒子的初始物理参数,包括:速度、压力、温度、参考温度、密度、运动粘度、重力、瑞利数、普朗特数。
当完成了前处理模块中的SPH算法设置和粒子模型建模,便可将数据导入计算模块,进行流体热运动的解算。
2.计算模块
所述计算模块用于使用SPH算法对流体的流动控制方程和修正方程进行求解,得到各时间步的各物理参数的数值模拟结果。该模块首先读入前处理模块所建立的仿真模型,从而获得SPH算法设置信息和粒子模型信息,进而开始SPH计算循环,在每个计算循环中,首先SPH算法会将流体的流动控制方程(质量守恒方程、动量守恒方程、能量守恒)施加到计算域内所有的粒子上;然后所有粒子的物理参数会按照SPH算法离散后的控制方程产生变化;最后当所有粒子完成一次计算后,求解前处理模块中选取的修正方程,使需要修正的物理参数产生修正变化,由此完成所有粒子物理参数的更新,该时间步完成;当完成粒子状态更新后,按照时间积分格式更新时间步,并同时输出包含所有粒子物理参数的数据文件,然后进入下一循环,直至时间步迭代到模拟时长结束。
3.后处理模块
所述后处理模块,用于对计算模块输出的仿真数据进行动态可视化处理,以直观展现流体热运动过程中流体的运动过程及各个参数的变化过程。该模块的处理过程如下:首先读入计算模块输出的数据文件,进而根据每个时间步仿真数据中的位置、速度、温度、密度等物理参数,将粒子位置分布以坐标点的形式表示其分布,得到整个计算域内粒子分布图,将其他物理参数用连续变化的颜色表示其分布,得到整个计算域内物理参数分布云图,最后利用所有时间步的粒子分布图及物理参数分布云图生成整个模拟过程的粒子运动动画和各物理参数变化动画。
有益效果
本无网格CFD仿真系统由于采用了SPH算法,较之已有的网格化CFD仿真系统,在模拟某些高度非线性问题(如射流问题、海洋运动问题、溃坝问题、动态边界问题)时,会得到更为合理的仿真结果;
本无网格CFD仿真系统较之已有的无网格CFD仿真系统存在如下优点:
1)添加了能量守恒方程的求解,从而使本仿真系统得以解决热力学问题;
2)添加了多种修正方程的求解,从而提升了系统的计算稳定性和计算精度;
3)添加了湍流模型,从而使本仿真系统可以解决流体热运动湍流问题。
附图说明
图1是本发明系统的组成结构以及数据流向示意图;
图2为本发明实施例的封闭方腔自然对流物理模型;
图3为本发明实施例的粒子模型;
图4为本发明实施例某一时间步所生成的流场温度分布云图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清晰明白,下面结合实施例和附图,对本发明实施例做进一步的详细说明。
下面以封闭方腔自然对流模拟为例介绍本发明基于SPH算法的流体热运动仿真系统的实施过程。
按图1所示流程实现该实施例的仿真模拟,具体实现步骤如下:
1.前处理模块包括:
1)SPH算法设置
计算方式:GPU并行。
SPH算法格式设置如下,时间积分格式:蛙跳法;核函数格式:五次样条核函数;粘性结算格式:层流粘性。
修正方法:人工位移技术。
SPH算法参数:粒子间距离dx=0.01、时间步长10-5、模拟时长0.5、核函数的影响半径(光滑长度)h=1.5dx=0.015、人工位移技术强度参数0.04、人工声速400。
2)粒子模型建模
图2为封闭方腔自然对流物理模型,其中L(L=1)表示物理模型边长,其中表示该部分粒子采用了恒温边界条件,TH和TC分别表示恒高温和恒低温,g表示重力的方向。
粒子坐标和粒子属性:以图2中正方形的四个顶点坐标为控制点,建立如图3所示的粒子模型,内部正方形域设置为流体粒子,正方形域外围设置三层边界(固体)粒子。
边界条件:四个壁面均采用固体壁面边界条件,采用边界力法;左右壁面采用恒温温度边界条件,上下壁面采用绝热温度边界条件。
初始条件:速度0、压力0、左侧壁面温度为1,右侧壁面温度为0、参考温度为0.5、密度1.225、运动粘度ν=10-6、重力g=9.81、瑞利数Ra=105、普朗特数Pr=0.71(以上均为无量纲形式)。
2.计算模块进行以下处理:
1)读入前处理模块输入的信息;
2)根据前处理模块所输入的信息,设置SPH算法与粒子模型;
3)开始循环,根据所选计算方法与参数利用SPH算法解算质量守恒方程、动量守恒方程、能量守恒方程和人工位移技术方程;
4)根据解算结果,更新粒子信息,包括:速度、位置、温度、密度;
5)输出该时间步的数据文件;
6)通过时间更新进入下一循环,直至时间步迭代到模拟时长循环结束。
3.后处理模块进行以下处理:
1)读入计算模块生成的数据;
2)生成单个时间步的粒子分布图及物理参数分布云图,如图4所示为根据某一时间步数据生成的流场温度分布云图,可以看出所得温度分布较为合理;
3)用所有时间步的粒子分布图及物理参数分布云图生成整个模拟过程的粒子运动动画和各物理参数变化动画,以动态展示流场的变化情况。
以上所述仅为本发明的具体实施例,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于SPH算法的流体热运动仿真系统,其特征在于,包括前处理模块、计算模块和后处理模块;所述前处理模块用于根据实际模拟目标建立CFD仿真模型,所述计算模块用于使用SPH算法对流体的流动控制方程和修正方程进行求解,得到各时间步的仿真数据,所述后处理模块用于对计算模块输出的仿真数据进行动态可视化处理,以直观展现流体热运动过程中流体的运动过程及各个参数的变化过程;
前处理模块包含SPH算法设置和粒子模型建模两个子模块;
所述SPH算法设置模块决定计算模块如何使用SPH算法,包括以下内容:计算方式设置、SPH算法格式设置、修正方法设置和SPH算法参数设置;
a)计算方式设置决定了该仿真系统的计算方式,包括中央处理器(CPU)并行和图形处理器(GPU)并行两种计算方式;其中CPU并行通过标准消息传递界面技术实现,GPU并行通过通用并行计算架构实现;
b)SPH算法格式设置包括时间积分格式、核函数格式和粘性解算格式三项内容;其中时间积分格式决定了模拟时间推进的方式,可选格式包括:蛙跳法、韦尔莱积分法和辛积分法;核函数为SPH算法的最基本函数,其决定了SPH算法如何离散流动控制方程,可选格式包括:三次样条核函数和五次样条核函数;粘性解算格式为动量守恒方程中流体粘性的计算方式,可选格式包括:层流粘性、人工粘性和亚粒子模型湍流粘性;
c)由于SPH算法易于产生震荡误差,则需要在算法中引入一些修正方法,修正方法包括:密度重构技术和人工位移技术;其中密度重构技术用于修正密度计算的震荡误差,可选类型包括:薛帕德滤波和δ-SPH;人工位移技术用于修正粒子位置的计算误差;
d)SPH算法参数包括:粒子间距离、时间步长、模拟时长、核函数的影响半径、修正方程参数和人工声速;
所述粒子模型建模模块用于对粒子模型进行建模,其包括以下四部分内容:粒子坐标、粒子属性、边界条件和初始条件;
a)粒子坐标的生成通过输入相应外形参数实现,根据外形参数的几何特征及给定的粒子间距离,便可求得该几何形状内所有粒子的坐标;
b)粒子属性设置是将边界粒子属性设置为固体,其他粒子属性设置为流体的过程;
c)边界条件确定了边界固体粒子对内部流体粒子的作用方式,包括两方面内容:固体壁面边界条件和温度边界条件;
固体壁面边界条件是计算域边界部分粒子的属性为固体粒子,边界粒子位置不变,并对靠近的流体粒子产生一定作用力,作用力施加方式包括边界力法和镜像粒子法两种方式;其中边界力法通过施加排斥力实现固体边界;镜像粒子法通过施加排斥动量实现固体边界;
温度边界条件是当存在具有一定温度特性的边界时,该部分边界粒子需要设置温度边界条件,以使其温度特性影响其附近流体粒子,温度边界条件包括恒温边界条件和绝热边界条件两种;恒温边界条件情况下,该部分固体粒子温度恒定,持续与周围粒子进行热交换;绝热边界条件情况下,该部分固体粒子与其他粒子不发生任何热交换;
d)初始条件为计算域内粒子的初始物理参数,包括:速度、压力、温度、参考温度、密度、运动粘度、重力、瑞利数和普朗特数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410377151.9A CN104143027B (zh) | 2014-08-01 | 2014-08-01 | 一种基于sph算法的流体热运动仿真系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410377151.9A CN104143027B (zh) | 2014-08-01 | 2014-08-01 | 一种基于sph算法的流体热运动仿真系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104143027A CN104143027A (zh) | 2014-11-12 |
CN104143027B true CN104143027B (zh) | 2017-03-29 |
Family
ID=51852199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410377151.9A Expired - Fee Related CN104143027B (zh) | 2014-08-01 | 2014-08-01 | 一种基于sph算法的流体热运动仿真系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104143027B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110633528A (zh) * | 2019-09-18 | 2019-12-31 | 河海大学常州校区 | 基于光滑粒子流体动力学的太阳能热水器加热仿真方法 |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104462830A (zh) * | 2014-12-12 | 2015-03-25 | 武汉大学 | 一种基于gpu加速的混合粒子的实时血流血管耦合方法 |
CN105279781B (zh) * | 2015-10-23 | 2018-06-08 | 山东师范大学 | 基于多精度融合的流体动画生成方法 |
KR20170067374A (ko) * | 2015-12-08 | 2017-06-16 | 삼성전자주식회사 | 유체 입자에 기반하여 공기 방울의 움직임을 모델링하는 방법 및 장치 |
CN105631100B (zh) * | 2015-12-23 | 2018-10-26 | 浙江大学 | 水场景目标红外尾迹特性的流体模拟方法 |
CN105701308B (zh) * | 2016-01-20 | 2018-11-20 | 中国民用航空飞行学院 | 一种基于无网格流固耦合的降落伞分析方法 |
CN106528989B (zh) * | 2016-11-03 | 2019-05-03 | 英特工程仿真技术(大连)有限公司 | 一种分布式并行sph仿真方法 |
CN106529063A (zh) * | 2016-11-14 | 2017-03-22 | 宜兴八达流体技术有限公司 | 基于cfd技术的流体系统及其设计方法 |
JP6741252B2 (ja) | 2016-12-06 | 2020-08-19 | 富士通株式会社 | 流脈線可視化装置、流脈線可視化方法、および流脈線可視化プログラム |
CN108171800B (zh) * | 2016-12-07 | 2021-05-25 | 北京乌有园艺术有限公司 | 创造具有无穷造型的制品的方法 |
CN107563064B (zh) * | 2017-09-05 | 2021-03-19 | 河海大学 | 一种海啸波越浪过程的二维数值模拟方法 |
US11042674B2 (en) * | 2017-10-10 | 2021-06-22 | Dassault Systemes Simulia Corp. | Acoustic effects of a mesh on a fluid flow |
CN107908918A (zh) * | 2017-10-19 | 2018-04-13 | 新疆大学 | 一种平坦沙床上沙粒冲击起动的sph数值模拟方法 |
CA3081825A1 (en) * | 2017-11-07 | 2019-05-16 | Aveva Software, Llc | Interactive guidance system for selecting thermodynamics methods in process simulations |
CN108491619B (zh) * | 2018-03-19 | 2020-05-26 | 浙江大学 | 基于物理与非物理混合的复杂场景流固耦合高效模拟方法 |
CN109088776B (zh) * | 2018-09-12 | 2020-07-28 | 西安交通大学 | 一种基于超级计算机的并行cfd计算与通信重叠优化方法 |
CN109408915B (zh) * | 2018-10-11 | 2022-10-14 | 北京动力机械研究所 | 固体火箭超燃冲压发动机燃烧流场仿真方法 |
CN109670212A (zh) * | 2018-11-28 | 2019-04-23 | 北京清华同衡规划设计研究院有限公司 | 天然气烟气余热回收模块化仿真软件 |
CN109783935B (zh) * | 2019-01-15 | 2020-12-11 | 北京航空航天大学 | 一种基于isph提高飞溅流体稳定性的实现方法 |
CN110162821B (zh) * | 2019-03-08 | 2022-09-13 | 南京理工大学 | 一种计算鸟撞高速旋转发动机叶片的方法 |
CN110083866A (zh) * | 2019-03-26 | 2019-08-02 | 中国核电工程有限公司 | 一种建筑物内部水淹的模拟系统及方法 |
CN110096766B (zh) * | 2019-04-15 | 2020-09-18 | 北京航空航天大学 | 一种基于物理的三维云运动演化方法 |
CN110309543B (zh) * | 2019-05-31 | 2023-04-18 | 中国航天空气动力技术研究院 | 一种多介质流体界面运动的仿真流程设计方法 |
CN110378025B (zh) * | 2019-07-22 | 2022-10-25 | 中国空气动力研究与发展中心 | 一种用于计算热化学非平衡绕流的处理系统及方法 |
CN110737935B (zh) * | 2019-09-18 | 2022-07-22 | 同济大学 | 一种基于数字孪生的室内热环境建模方法 |
CN110956005A (zh) * | 2019-11-27 | 2020-04-03 | 北京航空航天大学 | 一种基于粒子聚类和拉格朗日方法的高压消防水枪仿真方法 |
CN110909513B (zh) * | 2019-11-27 | 2021-07-30 | 北京航空航天大学 | 一种基于物理的油水混合现象可视化仿真方法 |
CN111737891B (zh) * | 2020-05-11 | 2024-05-10 | 自然资源部第二海洋研究所 | 一种海底热液流体循环三维仿真方法 |
CN112069745B (zh) * | 2020-09-10 | 2022-02-08 | 西北工业大学 | 一种固体推进剂废料切割处理的数值仿真方法及系统 |
CN112069689B (zh) * | 2020-09-10 | 2022-02-08 | 西北工业大学 | 一种航空发动机燃油雾化特性的仿真方法及系统 |
CN112100939B (zh) * | 2020-09-14 | 2023-06-16 | 福建天晴在线互动科技有限公司 | 一种基于Compute Shader的实时流体仿真方法及其系统 |
CN112685976A (zh) * | 2021-01-05 | 2021-04-20 | 西安交通大学 | 一种运动条件下自然循环系统数值模拟方法 |
CN113158531B (zh) * | 2021-02-07 | 2022-06-21 | 南开大学 | 一种利用形变梯度的单组分与多组分不可压缩流体仿真方法 |
CN115248989A (zh) * | 2021-04-25 | 2022-10-28 | 北京航空航天大学 | 基于屈服准则约束的粘性流体仿真方法 |
CN113378445B (zh) * | 2021-05-10 | 2024-02-02 | 中国科学院过程工程研究所 | 一种基于离散模拟的气液多相系统计算方法及系统 |
CN113486566B (zh) * | 2021-07-13 | 2022-05-27 | 南开大学 | 一种使用虚拟粒子相进行多孔介质和多相流体混合仿真的方法 |
CN113836835B (zh) * | 2021-08-28 | 2022-06-21 | 西安交通大学 | 一种核反应堆燃料熔化熔融物迁徙行为流固耦合分析方法 |
CN113761812B (zh) * | 2021-09-09 | 2024-03-29 | 西安交通大学 | 基于无网格法的单层壁面复杂几何流域的求解方法及系统 |
CN113947003B (zh) * | 2021-10-15 | 2024-05-07 | 西安交通大学 | 一种面向热流耦合场景的粒子型无网格仿真系统 |
CN114580252A (zh) * | 2022-05-09 | 2022-06-03 | 山东捷瑞数字科技股份有限公司 | 一种用于流体仿真的图神经网络仿真方法及系统 |
CN115630559B (zh) * | 2022-12-08 | 2023-03-10 | 中国人民解放军国防科技大学 | 一种基于粒子网格适配算法的流固耦合方法以及装置 |
CN116595884B (zh) * | 2023-05-24 | 2024-05-14 | 上海交通大学 | 微型核反应堆系统超实时计算实现方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101278292A (zh) * | 2004-12-14 | 2008-10-01 | 施卢默格控股有限公司 | 用于油藏模拟器中的耦合应力/液体流动的有限容积方法 |
-
2014
- 2014-08-01 CN CN201410377151.9A patent/CN104143027B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101278292A (zh) * | 2004-12-14 | 2008-10-01 | 施卢默格控股有限公司 | 用于油藏模拟器中的耦合应力/液体流动的有限容积方法 |
Non-Patent Citations (5)
Title |
---|
"一种改进的光滑粒子流体动力学前处理方法";雷娟棉 等;《物理学报》;20140723;第63卷(第14期);第144702-1—144702-11 * |
"光滑粒子动力学方法在复杂流动中的研究进展";周光正 等;《化工学报》;20140430;第65卷(第4期);第1145-1161页 * |
"基于CPU-GPU混合加速的SPH 流体仿真方法";胡鹏飞 等;《计算机工程与科学》;20140731;第36卷(第7期);第1231-1237页 * |
"封闭方腔自然对流的涡结构和传热特性";胡俊 等;《北京理工大学学报》;20140331;第34卷(第3期);第238-303页 * |
"流体运动SPH数值模拟后处理软件的设计与开发";许艳丽 等;《中国水运》;20140331;第14卷(第3期);第109-110页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110633528A (zh) * | 2019-09-18 | 2019-12-31 | 河海大学常州校区 | 基于光滑粒子流体动力学的太阳能热水器加热仿真方法 |
CN110633528B (zh) * | 2019-09-18 | 2022-08-16 | 河海大学常州校区 | 基于光滑粒子流体动力学的太阳能热水器加热仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104143027A (zh) | 2014-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104143027B (zh) | 一种基于sph算法的流体热运动仿真系统 | |
Rasmussen et al. | Directable photorealistic liquids | |
Wang et al. | Numerical investigation of unsteady cloud cavitating flow around the Clark-Y hydrofoil with adaptive mesh refinement using OpenFOAM | |
CN110717269B (zh) | 一种基于网格和粒子耦合的流体表面细节保护方法 | |
CN110298105A (zh) | 饱和多孔介质大变形分析的ccpdi-impm方法 | |
CN102915552B (zh) | 一种可控的火焰动画生成方法及其系统 | |
CN109344450B (zh) | 基于pbf的流体凝固模拟方法及系统 | |
CN113947003A (zh) | 一种面向热流耦合场景的粒子型无网格仿真系统 | |
CN108269299A (zh) | 一种基于sph方法近似求解的粘性流体建模方法 | |
CN105006015A (zh) | 基于位置的流体模拟控制方法及系统 | |
CN103426196B (zh) | 一种流体环境下的关节动画建模方法 | |
CN106934192A (zh) | 一种参数优化的浅水方程模型水体建模方法 | |
CN103389649B (zh) | 一种基于球面拼接网格的飞行器机动运动模拟方法 | |
CN103839287B (zh) | 一种用湍流方程丰富火焰模拟细节的方法 | |
Zhang et al. | Position-based fluid control | |
Liu et al. | Modified ghost fluid method with acceleration correction (MGFM/AC) | |
Jang et al. | Multilevel vorticity confinement for water turbulence simulation | |
Feng et al. | Detail‐preserving SPH fluid control with deformation constraints | |
Wei et al. | Melting and flowing of viscous volumes | |
Im et al. | Visual simulation of rapidly freezing water based on crystallization | |
Zhang et al. | Hybrid modeling of multiphysical processes for particle‐based volcano animation | |
CN109063375B (zh) | 基于保密度和无散度的不可压缩流体的模拟方法及系统 | |
Liu et al. | Viscosity-based vorticity correction for turbulent sph fluids | |
CN105260619B (zh) | 一种改进的kgf-sph方法 | |
CN110929450A (zh) | 一种基于sph的真实感湍流模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170329 Termination date: 20200801 |