美国人编制的溃坝的sph源程序，非常经典，适合初学者入门使用！

使用opengl对水流进行模拟，由于初学，水流效果很粗糙。本人电脑i5 4核，内存4G。水流粒子数5000，没有进行任何加速（没有使用领域搜索法），电脑跑起来只有10fps左右。初学opengl或者光滑粒子流体动力学(sph)的大神可以下载下来看一看。不过要自己配置opengl库对代码进行编译...

干货！这是本人基于光滑粒子流体动力学做的一个杯中水的模拟，比较成功！两年的sph学习不容易，今天特意拿出来与你们分享——sph方法的编程核心技巧！

一本目前为止最好的fluent学习书本 第一章 流体力学基础与FLUENT简介 第一节 概论 一、流体的密度、重度和比重 二、流体的黏性——牛顿流体与非牛顿流体 三、流体的压缩性——可压缩与不可压缩流体 四、液体的表面张力 第二节 流体力学中的力与压强 一、质量力与表面力 二、绝对压强、相对压强与真空度 三、液体的汽化压强 四、静压、动压和总压 第三节 能量损失与总流的能量方程 一、沿程损失与局部损失 二、总流的伯努里方程 三、人口段与充分发展段 第四节 流体运动的描述 一、定常流动与非定常流动 二、流线与迹线 三、流量与净通量 四、有旋流动与有势流动 五、层流与湍流 第五节 亚音速与超音速流动 一、音速与流速 二、马赫数与马赫锥 三、速度系数与临界参数 四、可压缩流动的伯努里方程 五、等熵滞止关系式 第六节 正激波与斜激波 一、正激波 二、斜激波 第七节 流体多维流动基本控制方程 一、物质导数 二、连续性方程 三、N—S方程 第八节 边界层与物体阻力 一、边界层及基本特征 二、层流边界层微分方程 三、边界层动量积分关系式 四、物体阻力 第九节 湍流模型 第十节 FLUENT简介 一、程序的结构 二、FLUENT程序可以求解的问题 三、用FLUENT程序求解问题的步骤 四、关于FLUENT求解器的说明 五、FLUENT求解方法的选择 六、边界条件的确定 第二章 二维流动与传热的数值计算 第一节 冷、热水混合器内部二维流动 一、前处理——利用GAMBIT建立计算模型 第1步 确定求解器 第2步 创建坐标网格图 第3步 由节点创建直线 第4步 创建圆弧边 第5步 创建小管嘴 第6步 由线组成面 第7步 确定边界线的内部节点分布并创建结构化网格 第8步 设置边界类型 第9步 输出网格并保存会话 二、利用FLUENT进行混合器内流动与热交换的仿真计算 第1步 与网格相关的操作 第2步 建立求解模型 第3步 设置流体的物理属性 第4步 设置边界条件 第5步 求解 第6步 显示计算结果 第7步 使用二阶离散化方法重新计算 第8步 自适应性网格修改功能 小结 课后练习 第二节 喷管内二维非定常流动 一、利用GAMBIT建立计算模型 第1步 确定求解器 第2步 创建坐标网格图和边界线的节点 第3步 由节点创建直线 第4步 利用圆角功能对I点处的角倒成圆弧 第5步 由边线创建面 第6步 定义边线上的节点分布 第7步 创建结构化网格 第8步 设置边界类型 第9步 输出网格并保存会话 二、利用FLUENT进行喷管内流动的仿真计算 第1步 与网格相关的操作 第2步 确定长度单位 第3步 建立求解模型 第4步 设置流体属性 第5步 设置工作压强为0 atm 第6步 设置边界条件 第7步 求解定常流动 第8步 非定常边界条件设置以及非定常流动的计算 第9步 求解非定常流 第10步 对非定常流动计算数据的保存与后处理 小结 课后练习 第三节 三角翼的可压缩外部绕流 一、利用GAMBIT建立计算模型 第1步 启动Gambit，并选择求解器为FLUENT5／6 第2步 创建节点 第3步 由节点连成线 第4步 由边线创建面 第5步 创建网格 第6步 设置边界类型 第7步 输出网格文件 二、利用FLUENT进行仿真计算 第1步 启动FLUENT 2D求解器并读入网格文件 第2步 网格检查与确定长度单位 第3步 建立计算模型 第4步 设置流体材料属性 第5步 设置工作压强 第6步 设置边界条件 第7步 利用求解器进行求解 第8步 计算结果的后处理 小结 课后练习 第四节 三角翼不可压缩的外部绕流(空化模型应用) 第1步 启动FLUENT 2D求解器并读入网格文件 第2步 网格检查与确定长度单位 第3步 设置求解器 第4步 设置流体材料及其物理性质 第5步 设置流体的流相 第6步 设置边界条件 第7步 求解 第8步 对计算结果的后处理 小结 课后练习 第五节 VOF模型的应用 一、利用GAMBIT建立计算模型 第1步 启动GAMBIT并选择FLUENT5／6求解器 第2步 建立坐标网格并创建节点 第3步 由节点连成直线段 第4步 创建圆弧 第5步 创建线段的交点G 第6步 将两条线在G点处分别断开 第7步 删除DG直线和FG弧线 第8步 由边创建面 第9步 定义边线上的节点分布 第10步 在面上创建结构化网格 第11步 设置边界类型 第12步 输出网格文件并保存会话 二、利用FLUENT 2D求解器进行求解 第1步 读入、显示网格并设置长度单位 第2步 设置求解器 第3步 设置流体材料及属

《近海环境流体动力学数值模型》比较系统地介绍了设计和使用近海环境流体动力学数值模型的相关知识，内容包括海洋流体动力学若干基础模型和完整方程组及定解条件；离散动力学方程的一些数值方法和差分格式以及格式与物理规律的关系；展示了从二维至三维，由均匀海洋至非均匀海洋，由固定岸界至运动岸界，由规则矩形网格至无结构网格，由正问题至逆问题等诸多模型或设计模型的特殊方法。书中强调了对近海环境至关重要的近海天文潮和近海环流模型的特征；对与海洋环境动力学直接相关的物质输运模型，也做了普遍性的介绍。　《近海环境流体动力学数值模型》可作为海洋科学和环境科学相关专业的研究生与本科生的参考书，也可供开始涉足海洋环境动力学的科技工作者参考使用

《光滑粒子流体动力学——一种无网格粒子法》 2005年 湖南大学出版社 作者：韩旭

SPH光滑粒子流体动力学中英文都有，中文版本以及英文版的都有，拿去参考吧。光滑粒子流体动力学-一种无网格粒子法 第1章 绪论 1.1 数值模拟 1.1.1 数值模拟的作用 1.1.2 一般数值模拟的求解过程 1.2 基于网格的方法 1.2.1 拉格朗日网格 1.2.2 欧拉网格 1.2.3 拉格朗日网格和欧拉网格的结合 1.2.4 基于网格的数值方法的局限性 1.3 无网格法 1.4 无网格粒子法(MPMS) 1.5 MPMs的求解策略 1.5.1 粒子描述法 1.5.2 粒子近似 1.5.3 MPMS的求解过程 1.6 光滑粒子流体动力学(SPH) 1.6.1 SPH方法 1.6.2 SPH方法简史 1.6.3 本书中的SPH方法 第2章 SPH的概念和基本方程 2.1 SPH的基本思想 2.2 SPH的基本方程 2.2.1 函数的积分表示法 2.2.2 函数的导数积分表示法 2.2.3 粒子近似法 2.2.4 推导SPH公式的一些技巧 2.3 其他基本概念 2.3.1 支持域和影响域 2.3.2 物理影响域 2.3.3 particle—in-cell(PIC)方法 2.4 结论 第3章 光滑函数的构造 3.1 引言 3.2 构造光滑函数的条件 3.2.1 场函数的近似 3.2.2 场函数导数的近似 3.2.3 核近似的连续性 3.2.4 粒子近似的连续性 3.3 构造光滑函数 3.3.1 构造多项式光滑函数 3.3.2 一些相关的问题 3.3.3 光滑函数构造举例 3.4 数值测试 3.5 结论 第4章 SPH方法在广义流体动力学问题中的应用 4.1 引言 4.2 拉格朗日型的Navier—Stokes方程 4.2.1 有限控制体与无穷小流体单元 4.2.2 连续性方程 4.2.3 动量方程 4.2.4 能量方程 4.2.5 Navier-Stokes方程 4.3 用SPH公式解Navier-Stokes方程组 4.3.1 密度的粒子近似法 4.3.2 动量方程的粒子近似法 4.3.3 能量方程的粒子近似法 4.4 流体动力学的SPH数值相关计算 4.4.1 人工粘度 4.4.2 人工热量 4.4.3 物理粘度 4.4.4 可变光滑长度 4.4.5 粒子间相互作用的对称化 4.4.6 零能模式 4.4.7 人工压缩率 4.4.8 边界处理 4.4.9 时间积分 4.5 粒子的相互作用 4.5.1 最近相邻粒子搜索法(NNPS) 4.5.2 粒子对的相互作用 4.6 数值算例 4.6.1 在不可压缩流的应用 4.6.2 在自由表面流的应用 4.6.3 SPH对可压缩流的应用 4.7结论 第5章 非连续的SPH(DSPH) 5.1 引言 5.2修正光滑粒子法 5.2.1一维情况 5.2.2 多维情况 5.3 模拟非连续现象的DSPH公式 5.3.1 DSPH公式 5.3.2 非连续的确定 5.4 数值性能研究 5.5 冲击波的模拟 5.6 结论 第6章 SPH在爆炸模拟中的应用 6.1 引言 6.2 HE爆炸和控制方程 6.2.1 爆炸过程 6.2.2 HE的稳态爆轰 6.2.3 控制方程 6.3 SPH公式 6.4 光滑长度 6.4.1 粒子的初始分布 6.4.2 光滑长度的更新 6.4.3 优化和松弛过程 6.5 数值算例 6.6 应用SPH方法模拟锥孔炸药 6.7 结论 第7章 SPH在水下爆炸冲击模拟中的应用 7.1 引言 7.2 水下爆炸和控制方程 7.2.1 水下爆炸冲击的物理特性 7.2.2 控制方程 7.3 SPH公式 7.4 交界面处理 7.5 数值算例 7.6 真实爆炸模型与人工爆炸模型的比较研究 7.7 水介质缓冲模拟 7.7.1 背景 7.7.2 模拟设置 7.7.3 模拟结果 7.7.4 小结 7.8 结论 第8章 SPH方法在具有材料强度的动力学中的应用 8.1 引言 8.2 具有材料强度的动力学 8.2.1 控制方程 8.2.2 本构模型 8.2.3 状态方程 8.2.4 温度 8.2.5 声速 8.3 具有材料强度的动力学SPH公式 8.4 张力不稳定问题 8.5 自适应光滑粒子流体动力学(ASPH) 8.5.1 为什么需要ASPH方法 8.5.2 ASPH的主要思想 8.6 对具有材料强度的动力学的应用 8.7 结论 第9章 与分子动力学耦合的多尺度模拟 9.1 引言 9.2 分子动力学 9.2.1 分子动力学的基本原理 9.2.2 经典分子动力学 9.2.3 经典MD模拟 9.2.4 Poiseuille流的MD模拟 9.3 MD与FEM和FDM的耦合 9.4 MD与SPH的耦合 9.4.1 模型I：双重功能(具有重叠区域的模型) 9.4.2 模型Ⅱ：力桥(没有重叠区域的模型) 9.4.3

Dr. Liu 在《光滑粒子流体动力学--一种无网格方法》一书中提供的光滑粒子流体动力学（SPH）源程序，采用Fortran语言编写，用于计算可压缩流体力学问题。