![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
cfd
文章平均质量分 81
lizhengjiang
这个作者很懒,什么都没留下…
展开
-
空间划分与粒子相互作用(interaction) 实例
SPHysics 是采用SPH(光滑粒子流)方法研究流体动力学的开源软件。源码中用到了均匀空间划分方法和邻居粒子碰撞搜索(searching & interaction),作为完善前两篇的实例。本文学习SPHysics这部分代码(FORTRAN)。仅供参考。首先,SPH方法程序实现大致有如下几个过程:空间划分、构造粒子链表;碰撞(interaction)对搜索、粒子变量更新(加速度、速原创 2013-04-19 09:10:25 · 1139 阅读 · 0 评论 -
算子分裂算法介绍(Introduction to Operator Splitting algorithms)
算子分裂(operator splitting algorithm). 考虑 不可压缩NS方程, 包含非线性对流项,粘性扩散项,压力梯度项。如果采用一个算子描述这三项,那么该算子会有非常多限制,不便计算。而考虑将这三项拆开,那么将得到三个基本方程(对流方程,波松方程,拉普拉斯方程),而且基本方程的数值算法可选项非常丰富。这就是算子分裂的意义。当然,拆分后的算子与原算子到底在多大程度上一致,也是该类原创 2014-03-25 01:46:17 · 9138 阅读 · 0 评论 -
Introduction to FEM based CFD algorithms ( 基于有限元的CFD算法)
recommand materials: 1 Numerical Methods for Incompressible Flows, by Martin Burger : 介绍非线性单步迭代算法和时间步积分/分裂算法; 介绍NS方程弱解的矩阵形式 2 Numerical Methods for Incompressible VIscous Flow, by Ragnar Winthe原创 2014-03-24 01:26:06 · 1335 阅读 · 1 评论 -
计算不可压缩流体- NS方程求解算法
在上一篇介绍差分格式可以看到求解不可压缩NS方程,主要是求解速度的动量方程和压力的波松方程。本文介绍一些历史中的算法,下一文介绍目前常用的算法。一 Marker and Celler (MAC) 算法在自然差分格式中给出PPE如下: laplace ( p) = laplace ( body_force) - rho * div ( u * grad_u) (1) 。原创 2014-03-20 02:11:59 · 15127 阅读 · 2 评论 -
计算不可压缩流体--差分格式的两个问题
一 Cell - Re 问题 举例, 稳态BURGER'S 方程 u du/ dx = miu * d^2 u/ d x^2 , 中心差分离散 Re(u(i + 1) - u(i) )/ 2 - ( u(i - 1) - 2 u(i) + u(i+1) ) = 0 定义单元雷诺数, Re = U * h/ miu, % h 网格尺寸原创 2014-03-20 00:04:49 · 1781 阅读 · 0 评论 -
计算不可压缩流体 -- 数学基础
推荐教材: a lecture in computational fluid dynamics of incompressible flow: mathematics, algorithms and implementations. Author: J. M. McDonough, Depart. of Mechanical Engineering and Mathematics,原创 2014-03-19 02:25:54 · 3185 阅读 · 0 评论 -
level set method reinitialization (水平集方法的重初始化)
1 回顾 水平集方法LSM使用了背景网格。不同于VOF,MPM等采用背景网格节点插值界面的构造方式,在于LSM采用 phi 函数隐式描述界面, phi函数是从微分几何过来的,图像学里面叫,曲线演化。LSM对于几何信息(诸如法向量,曲率)描述精确。界面演化由全局演化方程描述,(即上一篇中使用的演化方程,并不局限解界面演化,实际上是用来更新全域演化信息。所以它需要将已知在界面上的速度原创 2014-03-16 05:23:30 · 3918 阅读 · 0 评论 -
level set method 水平集方法基本问题
版里搜了下,呵呵,发现level set method应该是从图像识别领域借过来的。 基本思想: 令(n+1)维函数 t = phi(x(t)) 等于常数, 得到n维曲线,即水平集。 基本问题: 1 初始化 定义 level set 函数 phi(x) 为符号距离函数,界面线上距离为0。一般为隐函数。原创 2014-03-13 13:55:58 · 4888 阅读 · 0 评论 -
碰撞测试 (Real Time Collision Detection Chapter 7)
7.1.6.1 one test at a time考虑当一个对象完整"落入"某单元里。对其碰撞测试不仅要包含同一单元内的其他对象还有该单元邻居单元内的对象。哪些邻居单元需要测试取决于对象是如何“占据”邻居单元的,以及采用对象的何种特征定义对象“占据”(在..内)某个单元。可选的对象特征:对象包围球球心(object bounding sphere center)以及轴对齐包围盒的最小顶点值(翻译 2013-04-17 17:08:31 · 1081 阅读 · 0 评论 -
SPHysics浮体模型及对应边界处理
SPH方法的特点是采用粒子求和来近似场变量。有效的粒子求和近似依赖空间划分、邻近粒子搜索等算法,如前文介绍。本文主要学习SPHysics中的浮体处理方法。 首先需要了解SPH与结构的耦合(刚体、弹性体)都是通过界面粒子实现的。界面粒子当然也属于一类边界粒子,文献中称为“free moving boundary particles”(FM)。相对于固定在壁面上为防止流体粒子四处乱跑而设计的M原创 2013-04-24 14:55:59 · 1831 阅读 · 0 评论 -
空间划分(real time collision detection chapter 7)
本章介绍三种空间划分方法:均匀网格划分、分层网格划分、树划分7.1 均匀网格划分7.1.1 链表存储 定义一个数组,然后将网格单元(grid cell)与数组元素一一对应。当一个单元中有对象(object)时,数组元素就是一个链表,保存在该单元内的所有对象,否则指NULL。最简单维持对象链表的办法是,在对象类型(结构体)中包含链表成员变量。该方法既可以提高数据缓翻译 2013-04-17 16:20:02 · 1285 阅读 · 0 评论 -
su2 上手测试
su2 项目 及 与其他开源cfd比较原创 2015-01-26 06:17:43 · 3315 阅读 · 0 评论