《无痛苦N-S方程笔记》勘误

[bestucan]

没搜到这个笔记的勘误帖子
那我就开了～

[bestucan]

第二页 第二段 倒数第三行 “大的颗粒移动速度较慢，小的颗粒移动速度较慢”->“小的颗粒移动速度较快”
第二页 第三段 倒数第四行 “多相流体力学主要通天实验室研究”->“通过”
第四页 (1.31)式第一个偏倒数应该是偏t
第十一页 倒数第二行 “如图2”->“如图2.2”
第十二页 第四行 “结合方程(2.8)和(2.30)”->“结合方程(2.30)和(2.32)”
第十二页 2.4.3节 上面一段 “从方程(2.9)推导到(2.34)”“因此积分形式的方程(2.9)”(2.9)好像引错了
第十六页 倒数第二行 “图8右图表示”->“图2.14”
第十七页 倒数第一段 第六行 “图10中的平面efgh”->“图2.16”
第十七页 倒数第一段 第九行 “图10中其他剪切应力”->“图2.16”

[李东岳]

大侠若留真名，我更新在书里，要不就用你的ID了

咦不对啊，怎么跟我的这个版本对不上? https://www.jianguoyun.com/p/DdvrhJMQ9s3ZBhjw7p4B

[bestucan]

@东岳 李老师客气了，我叫 凃灿。
我这您这里学到很多，做这些心甘情愿。
我是从东岳流体的主页上点进去下载的。也是坚果云。

[李东岳]

@bestucan 在 《无痛苦N-S方程笔记》勘误 中说：

第二页 第二段 倒数第三行 “大的颗粒移动速度较慢，小的颗粒移动速度较慢”->“小的颗粒移动速度较快”

比如这个， 没找到呢。其他很多也是。
第二页应该在介绍张量标识法，你的是么

[bestucan]

@东岳 李老师，不好意思，记笔记记窜了。
前面几个错误是《多相流及群体平衡模型的数值计算方法》里的。后面的是“笔记”里的。

[bestucan]

《多相流及群体平衡模型的数值计算方法》看完了，就又发现了一个错误，不值当再开新帖了。
第十一页 倒数第一段 第二行 两处“动理学方程”->“动力学方程”

[李东岳]

等下，一个一个来，无痛苦N-S方程笔记这个主要是哪几个？我没对上

第二页 第二段 倒数第三行 “大的颗粒移动速度较慢，小的颗粒移动速度较慢”->“小的颗粒移动速度较快”

第二页 第三段 倒数第四行 “多相流体力学主要通天实验室研究”->“通过”

第四页 (1.31)式第一个偏倒数应该是偏t

第十一页 倒数第二行 “如图2”->“如图2.2”

第十二页 第四行 “结合方程(2.8)和(2.30)”->“结合方程(2.30)和(2.32)”

第十二页 2.4.3节 上面一段 “从方程(2.9)推导到(2.34)”“因此积分形式的方程(2.9)”(2.9)好像引错了

第十六页 倒数第二行 “图8右图表示”->“图2.14”

第十七页 倒数第一段 第六行 “图10中的平面efgh”->“图2.16”

第十七页 倒数第一段 第九行 “图10中其他剪切应力”->“图2.16”

[bestucan]

@东岳 这个一大堆，除了前两个是多相流，后面的都是ns笔记的

[李东岳]

妥了，晚上吃完饭改一下，然后at你

[李东岳]

已勘误，多谢凃灿

《多相流及群体平衡模型的数值计算方法》这个还没写完，我捉摸要不要申请个国家的出版物资助 :w:

[bestucan]

不阔气:w:

[lin_lia]

老师您好，关于无痛苦N-S方程笔记中动量方程（2.68）的写法，是否能够写为$\rho \frac{DU}{Dt}= \triangledown \cdot \tau$这种形式。
是否这种形式这样更有助于理解本构关系

[李东岳]

@lin_lia 已更新，多谢建议

[lin_lia]

[bestucan]

@bestucan 在 《无痛苦N-S方程笔记》勘误 中说：

《多相流及群体平衡模型的数值计算方法》看完了，就又发现了一个错误，不值当再开新帖了。
第十一页 倒数第一段 第二行 两处“动理学方程”->“动力学方程”

这个不是错误。是我孤陋寡闻了

[.J.]

@bestucan 学习了，我们做颗粒动理学，只知道老师这么叫我们也跟着这么叫，原来还有这样一段故事。

[Frank0514]

李老师：今天学习您的《无痛苦N-S方程笔记》发现有点小笔误：
第10页公式（2.20）下面一句话：“现在看图4中的无穷小微团”，应该是图2.5吧；
第16页、17页、18页中有几个地方写的是“收到xxx力”，应该是“受到”才对。

[李东岳]

@Frank0514 已更新，谢谢！！

[FluidGao]

方程（1.4），是不是漏了一个转置符号？

[FluidGao]

@FluidGao 在 《无痛苦N-S方程笔记》勘误 中说：

方程（1.4），是不是漏了一个转置符号？

最新版的笔记已经更新

[李东岳]

@FluidGao 啊？我才看见

[wangfei9088]

最近几天再看这半部书，感觉收获很多，不只是加深印象，对很多问题也有了新的认识。基于2020年12月29日修订版，也就是东岳流体今天下载的那个版本有一些typo。说得不对的地方，还请东岳老师包涵。

第一页倒数第二行：
初学者也不看的不明不白。
是不是多了个“不”？

第二页第一行和第二行：
GCFD 课程因此分为两步，第一步是基本的N-S 方程入门，第二部是N-S 方程求解。第一步需要同学们通过本笔记进行预习。第二部将在GCFD 课程上讲授。
“部”应该是“步”？

第三十四页4.1.4节第二行：
另但是大量的文献依旧通过二维大涡模拟、甚至二维直接模拟进行相关研究。例如，Bouris and Bergeles 的文章中表示[7]
多了个“另”？

第三十五页第二段引用：
我们的研究表示二维的大涡模拟，甚至直接模拟完全是扯淡的。
研究表明？

第三十五页4.2节第一段倒数第二行：
定于网格节点与面相连接的序号等。
这个是“定义”？

第三十六页4.2.1节第二段第二行和第四行：
左侧速度的方向与面试量为近似同向，因此Sf · Uf 为正值。同时，owner 单元的流体流出，neighbour 单元有流体流入。右侧度的方向与面试量为近似反向，因此Sf · Uf 为负值。
两个应该都是“面矢量”？

第三十九页倒数第十一行：
在这种情况下，方程(4.29)变为：
这里应该是“方程(4.24)”？

第四十六页：

1. 方程4.69应该是最后一个“dβ”应该是“dγ”？

2. 然后对其求解散度项中代入到动量方程中进行显性离散计算。
   这个“中”是不是要换成“再”？

2021年第一个回帖，祝各位老师新的一年计算收敛，多发paper！新年快乐！！！

[李东岳]

@wangfei9088 非常感谢！已更新。最新版本的方程4.29、以及四十六页的没太对上，你看新版本的是哪个不？

[wangfei9088]

@李东岳 啊，东岳老师，刚重新对了一下东岳流体上的文档。不好意思，是我写错了。应该是这样的：

第四十一页：
方程4.43多了“Am”？

第四十一页倒数第八行：
在这种情况下，方程4.47变为：
这里应该是“方程(4.42)”？

第四十六页：
方程4.79最后一个“dβ”应该是“dγ”？

第四十七页第二段最后一行：
然后对其求解散度项中代入到动量方程中进行显性离散计算。
这个“中”是不是要换成“再”？

我写的勘误不一定都对，请东岳老师斟酌。

[李东岳]

@wangfei9088 非常感谢！全部都更正了！

顺便又更新了一下

• 麦克斯韦分布与高斯分布

• 矩方程、五矩模型、十矩模型以及欧拉方程

[wangfei9088]

@李东岳 作为伸手党，应该是我们感谢您！
其实这两个人都是我

[CUGB_Zhou]

受教了，谢谢

[Kingkong]

最近重看笔记，刚看了前言发现两个小笔误。
（1）前言第二段：Computational译为“计算”。Fluid Dynamics则为流体动力学。 计算后的句号是否应为逗号？
（2）IV页“总而言之”这段第5行“玩弄于掌股”是否为“玩弄于股掌”？
两处无关紧要的笔误(*￣︶￣)

[李东岳]

@kingkong 谢谢，已更新，你看下

[wangfei9088]

@李东岳
第95页（页码是83页）公式（7.27）中的
$$
\frac{\alpha}{m_0} ^ {2/3}
$$
加个括号应该会更好吧。
$$
(\frac{\alpha}{m_0}) ^ {2/3}
$$

[wangfei9088]

更新一个笔误。图中的$y^{*}$应该是
$$y^*=\frac{\sqrt{k} y}{\nu}$$

[李东岳]

谢谢！已更新

[wangfei9088]

图中的$r$应该是
\begin{equation}
r=\min\left(\frac{\nu+\nu_t}{|\nabla \mathbf{U}| \kappa^2 \tilde{d}^2}, 10\right)
\end{equation}
后面$\tilde{d}$更严格的写法应该是
\begin{equation}
\tilde{d}=\max\left[y-\left(1- \tanh\left(\left(8r\right)^3\right)\right) \max\left(y-C_{DES}\Delta, 0\right), 0\right]
\end{equation}

请东岳老师斟酌。

[李东岳]

@wangfei9088 多谢大佬检查。稍后更新。大佬把湍流模型捋了个遍么。做什么研究方向的

[wangfei9088]

@李东岳 哈哈，东岳老师，我不是大佬。小青铜科研找不到方向和突破口，只能多看看书。 确实这次基本把湍流模型都撸了一遍，模型加代码。每次看这本笔记都收获很大。

[李东岳]

@wangfei9088 厉害厉害，最近我把壁面函数更新一下

[wangfei9088]

@李东岳 学无止境，接着嗨

[袁宝强]

LienLeschziner低雷诺数湍流模型勘误。

Openfoam *.c文件中源码和出处

tmp<volScalarField> LienLeschziner::fMu() const
{
    const volScalarField yStar(sqrt(k_)*y_/nu());

    return
        (scalar(1) - exp(-Anu_*yStar))
       /((scalar(1) + SMALL) - exp(-Aeps_*yStar));
}

以下是原论文出处
Lien, F., and Leschziner, M. A. (December 1, 1993). "A Pressure-Velocity Solution Strategy for Compressible Flow and Its Application to Shock/Boundary-Layer Interaction Using Second-Moment Turbulence Closure." ASME. J. Fluids Eng. December 1993; 115(4): 717–725.

[李东岳]

@袁宝强 已更新，多谢老铁

[ESAYUAN]

@李东岳 东岳老师 我在学习过程中跟着公式推导 发现3.4.4条中的最后一步应该是将方程（3.51）带入到方程（3.50），是笔误吧？

[Tong]

hello李老师~ 最近看了看这本笔记，发现有的地方有点小问题~~~。
（1）多孔介质部分，P115页，公式7.58的以及后文的注解中，您将D解释为渗透率(permeability or k)，这应该是有些不准确的，根据达西公式描述的压力速度关系，渗透率的量纲为m^2 (https://www.calculator.org/properties/permeability.html)，通常左乘压力梯度(可以参考[1]中达西公式描述(Equation-(1)))，OF中的多孔介质模型实现是在动量方程添加阻力源项来实现的，与达西公式不同的是：在公式7.59中，D左乘了速度(U)，大小为渗透率的倒数，量纲为m^-2，为了避免误解，在这里解释D为达西阻力系数，F为 Forchheimer 阻力系数更合适，多孔介质渗透率(k)越大对应的达西阻力系数（D）就越小，反之亦然。个人观点：将” 沉降项 “翻译成”阻力源项“更好理解，也和Fluent的手册里对这种实现方法的解释一致。
(2)同时在 P125页中 7.9.6 setReference()， 对公式7.101的解释中，‘p无解’。实际上该方程的秩（rank）为 2（矩阵的第1、2行与第3行 线性相关），而未知数的个数为3，这种情况下，p有无穷多个线性相关解，这种情况下设置参考值能够将这无穷多个线性相关解退化到该参考值下的定解。同时，建议将rank翻译成’秩‘，符合国内线代教材对矩阵所在的线性空间的维数这一定义，用‘阶’的话容易和‘m*n阶’矩阵这种叫法中的‘阶’混淆造成误解。

[1]T. Wang and J. Yao, An improved embedded discrete fracture model and domain connectivity algorithms on 3D
unstructured grids, Journal of Computational Physics, 111142, doi: https://doi.org/10.1016/j.jcp.2022.111142. (求引用QAQ)

[李东岳]

@tong 非常感谢老铁！已更新。老铁留个地址，给你有几个CFD记事本 http://dyfluid.com/profs.html

[李东岳]

@tong 哎我看下你的文章 我刚发现你发了一个自己的JCP 我下载看一下

[李东岳]

@tong 明天给你处理，前几天疫情，这面快递都停发了，不好意思老铁

[Tong]

@李东岳
不要紧~ 啥时候都行

[Harvest]

李老师，我个人觉得公式3.7，3.8的表述有点问题，如果在某一个位置（固定坐标）求导数，那么就是偏导数，而不是全导数，全导数下面是不应该加坐标限制的，也许是理解不太一样