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

wangfei9088

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

李东岳

谢谢！已更新

wangfei9088

图中的 $r$ 应该是
$\begin{matrix} (1) & r = m i n (\frac{ν + ν_{t}}{| \nabla U | κ^{2} {\tilde{d}}^{2}}, 10) \end{matrix}$
后面 $\tilde{d}$ 更严格的写法应该是
$\begin{matrix} (2) & \tilde{d} = m a x [y - (1 - \tanh ({(8 r)}^{3})) m a x (y - C_{D E S} Δ, 0), 0] \end{matrix}$

请东岳老师斟酌。

李东岳

@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的表述有点问题，如果在某一个位置（固定坐标）求导数，那么就是偏导数，而不是全导数，全导数下面是不应该加坐标限制的，也许是理解不太一样

知行

@李东岳李老师，基于230105修改版有几处不明确的地方请您过目:
1.p96 第3个公式右端转置符号是否应该为正体的 $T$ ，类似情况在公式(2.4）中也有出现
2.p102 第二个公式中间第三项 $ϵ$ 是否应该为 $ε$
3.pIII(前言) 第二行,“来预测天气预报” $\to$ “来预测天气” or "来进行天气预报"

知行

@知行 4.p103 kOmega部分:第四行公式 $D_{ε}$ 是否应为 $D_{ω}$

知行

@知行 5.p63、p96~p107 每页的转置符号 $\to$ $T$

疏影横斜水清浅

81页：方程（4.179）应该改为方程（4.170）

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

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

CFD中文网

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