《无痛苦N-S方程笔记》勘误
-
图中的$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 厉害厉害，最近我把壁面函数更新一下
-
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.
-
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) -
只能多看看书。
确实这次基本把湍流模型都撸了一遍，模型加代码。每次看这本笔记都收获很大。
学无止境，接着嗨
