CFD中文网

@wwzhao 谢谢。这个问题我现在已经搞清楚了，是为了增加矩阵的“对角占优”特性，提高收敛性。经过实验，假如修改成不加nu_t的公式，对于一些Co数较大的case就有可能计算发散。

@alex 看到你的网格了，感谢分享，感觉划分的非常不错，好厉害！学习了！

对于给定的NDF，划分为$i$个$N_{pp}$，对每个$i$上定义$k$阶矩$m_k^i$，给定$m_k^i$，可以计算第$i$区间的节点$d^i_0,d^i_1$以及权重$w^i_0,w^i_1$：
\begin{equation}
\begin{split}
w^i_0&=w^i_1=0.5
\\
d^i_0&=m_1^i-\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\\
d^i_1&=m_1^i+\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\end{split}
\end{equation}
对于仅考虑破碎的PBE：
\begin{equation}\label{pbe}
\frac{\p n(d)}{\p t}=\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'-g(d)n(d)
\end{equation}
对方程\eqref{pbe}在$i$上取$k$阶矩：
\begin{equation}\label{m}
\frac{\p m_k^i}{\p t}=\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')d^k\beta(d|d')n(d')\rd d'\rd d-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}
\begin{equation}
\begin{split}
\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'\rd d&=
\int_{d_{i-1/2}}^{d_{max}}g(d')n(d')\left(\int_{d_{i-1/2}}^{d'}\beta(d|d')\rd d\right)\rd d'
\\&=
\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)
\end{split}
\end{equation}
Therefore
\begin{equation}
\frac{\p m_k^i}{\p t}=\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}

那你先解释一下，为什么你认为这是个不可能事件

C1和C2差异很大，但Q3和w3几乎相同，请问如何解释

它违背了什么物理规律？和已有的什么结论相悖？

谢谢老师~~

应该还可以用这个办法来判断
if(THREAD_TYPE(t) == THRAD_F_WALL)

不是，这里很少有用OpenFOAM实现的

@李东岳 好的，谢谢李老师

邱老师现在玩的大项目，都是上亿网格的工业级别 :cheeky:

@CFD_xiaobai Luo破碎模型不是也是一个大公式么？

破碎这面我没太看中文的书，跟你的描述可能对不上。据我了解，破碎模型有俩种：

一种是breakge frequency(破碎频率)和daughtor size distribution (DSD)分开的， 一种是破碎频率+DSD写成一个大方程的；

比如Laakkonen是第一种，Luo是第二种，Liao也是第二种。是这样么？

@红豆沙 并没有。。。我们这边很神奇，几个学校的cluster是共享的，有各种CPU，5118，6142，EPYC 7501什么乱七八糟的，还有V100加速卡，然而我自己实验室给的一台破E5450，自己买了台Y540 打游戏。。。

@东岳 让孩子赢在起跑线上:high:

$\tau$这个雷诺应力，对角线的是normal stress，非对角线的是shear stress。不过你说的更详细，感谢分享！

请教一个比较初级的问题，查了帮助文档没有找到想得到的答案。
ICEM采用Octree方法划分非结构网格的时候，帮助文档和网上的一些教程中都写到：“建议全局网格参数设置中的Max element值要是2的幂”
这里2的幂是指：2的正整数幂？2的正数幂？2的实数幂？
实际测试过感觉比如填0.05啥的也可以划分，网格也确实要比2的密很多。所以想问下各位高手这里这个数值填写应该有什么经验或者规律呢？尤其是模型尺寸很小的话，有时候确实需要很小的尺寸啊~

@evensun 嗯，涡粘系数是二阶张量这个按我理解其实就是非线性湍流模型了

@禾斗石开 这个应该都是读取时候的提示，提示“正在从case中读取旧格式的信息”，能显示出来就说明读取成功了。

这个问题已经解决，麻烦大家了

@程迪 在 GPU VS CPU 中说：

隐式的得回头修改这个自由度那个自由度

你好，请问您所说的自由度是什么意思，能够稍微解释一二，不胜感激，谢谢！

@李东岳 好的谢谢您，我再试试

安装nodejs

curl -sL https://deb.nodesource.com/setup_12.x | sudo -E bash - sudo apt-get install -y nodejs

安装redis

sudo apt-get install software-properties-common sudo add-apt-repository ppa:chris-lea/redis-server sudo apt-get update sudo apt -y install redis-server

创建用户

adduser dyfluid

更改编辑器：

sudo update-alternatives --config editor

添加用户权限：

visudo

dyfluid ALL=(ALL:ALL) ALL

关闭停止：

/etc/init.d/redis-server stop /etc/init.d/redis-server start cd $HOME git clone https://github.com/NodeBB/NodeBB.git nodebb cd nodebb ./nodebb setup

清楚ssh

rm ~/.ssh/known_hosts

Nginx做ln -s注意使用绝对路径

ln -s /etc/nginx/sites-available/example.com.conf /etc/nginx/sites-enabled/example.com.conf systemctl stop nginx.service systemctl start nginx.service systemctl status nginx.service