摘要: 1 控制方程 控制方程如下: \[ \frac{\part u}{\part t} + \frac{\part f(u)}{\part x} = s(x), \quad x \in [a, b] \] 我们通过各种数值方法来完成上述方程的求解。本软件采用有限差分法来离散方程,为了编程的方便,暂时采用 阅读全文
posted @ 2021-02-26 23:00 catfluid 阅读(342) 评论(0) 推荐(0)
摘要: 1. 第一层网格节点到壁面距离 1.1 雷诺数计算: \[ Re=\frac{\rho u d}{\mu} \] 1.2 壁面摩擦系数计算: 更多公式参考http://www.cfd-online.com/Wiki/Skin_friction_coeffcient The skin friction 阅读全文
posted @ 2022-03-21 20:28 catfluid 阅读(2561) 评论(0) 推荐(0)
摘要: 我们知道,现阶段数值模拟湍流的方法主要有三个层次: 雷诺平均(RANS)应用最广泛,但雷诺平均改变NS方程的瞬时性,湍流模型的封闭一直是其难点; 大涡模拟(LES)方法建立在人们对湍流的基本认识:大尺度涡由流场边界、流体粘性等决定,具有明显的各向异性和非均匀性,小尺度涡基本是均匀各向同性的。通过滤波 阅读全文
posted @ 2022-03-20 10:51 catfluid 阅读(1271) 评论(0) 推荐(0)
摘要: 一些常用的流量限制器/斜率限制器 CHARM [not 2nd order TVD] (Zhou, 1995) \[ \phi_{cm}(r)=\left\{ \begin{array}{ll} \frac{r\left(3r+1\right)}{\left(r+1\right)^{2}}, \qu 阅读全文
posted @ 2021-02-26 23:17 catfluid 阅读(999) 评论(0) 推荐(0)
摘要: 1. 重要概念 1.1 单调格式: 如果一个格式可以写成:$u_i^{n+1}=f(u_n,u_{i-p+1}n, \dots, un_{i+q})$的形式,并且对于每一个变元$u_kn,k=i-p,i-p+1,\dots,i+q$都有: \[ \frac{\part f(u_{i-p}^n,u_{ 阅读全文
posted @ 2021-02-26 22:54 catfluid 阅读(802) 评论(1) 推荐(1)
摘要: archlinux基本配置:更换清华源,配置fcitx5输入法。 阅读全文
posted @ 2020-08-07 14:01 catfluid 阅读(870) 评论(0) 推荐(0)
摘要: 布拉休斯方程如下: $$ \begin{equation} f f^{''}+2f^{'''}=0 \\ f(0)=f^{'}(0)=0;f^{''}(0)=1 \end{equation} $$ 这是一个非线性常微分方程,下面我们利用四阶龙格库塔方法来求解该方程。 我们引入新的变量: $$ \be 阅读全文
posted @ 2020-04-03 09:04 catfluid 阅读(1650) 评论(0) 推荐(0)
摘要: 0. 序言 在开始正文之前,请允许我先说明一下本文的目的和写作的动机,好让读者不惑。 我们知道,在Linux环境中,很多软件的组织都遵循GNU软件标准。不论是自己开发GNU软件还是阅读别人写好的源程序,能了解、熟悉GNU项目的构建方式,对我们的工作会起到事半功倍的效果。本文的目的,就是从零开始,告诉 阅读全文
posted @ 2019-09-09 15:12 catfluid 阅读(820) 评论(0) 推荐(2)
摘要: [TOC] 1. 什么是lxml lxml是干什么的?简单的说来,lxml是帮助我们解析HTML、XML文件,快速定位,搜索、获取特定内容的Python库。我们知道,对于纯文本的HTML文件的查找可以使用 、`BeautifulSoup`等完成。lxml也是对网页内容解析的一个库。 那么为什么要用l 阅读全文
posted @ 2019-08-16 16:51 catfluid 阅读(9523) 评论(0) 推荐(1)
摘要: 说明 plot3d格式是NASA制定并大量使用的CFD网格文件格式,在CFD编程过程中经常涉及到。本文利用Python语言编写一个读取plot3d文件,写入plot3d文件的通用类,并利用vtk完成网格的显示。 对于程序的说明 使用说明 在主函数中,写入如下代码: # 新建对象 mesh = Mes 阅读全文
posted @ 2019-07-15 21:55 catfluid 阅读(3749) 评论(0) 推荐(0)