CFD中文网

@sungda 高手，默认为8，调节为4就好了 :quwan:

@txwwbz
OF5指南我从去年10 1就开始弄，一直很零散的弄，导致进度很慢。

2014年出OpenFOAM指南的时候，当时在读博士，无独立科研能力，时间很多，导师不给我活，我就闲着，索性翻译了。

现在有了独立科研能力，一直都在搞新的算法，指南一直没更新。

我看看下周五之前，发布OpenFOAM指南5.0

@李东岳 请问在哪里下载？东岳流体？

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

图在哪？
感觉描述的很有意思
再试试？？

看到自己提交的C-CFD了哈哈哈

明明是和zeroGradient等边界条件放在一起的，应该是个基础边界条件吧？为什么在使用的时候显示unknown呢

@izumi 感谢。

链接文本

参考页面左边栏。

小白试了你更新的代码，运行时出现下面的错误，请问是我编译位置错了吗0_1521630956325_error.png 图片地址)

@一朵苍凉 我也不太懂，可以看看控制领域的相频特性，也许有帮助:zoule:

可以吧，fluent给定的周期边界的案例就是类似的

请问燃料燃烧时这种‘理论计算产生的平均总温’应该怎么计算啊？翻了基本燃烧学的书都没这方面的定义。温升燃烧效率.png

自己定义了一个类，如图，通过该类求解出来的数据在进程间不连续，我应该怎么写进程间的通信，有没有大神碰见过和我一样的问题，或者openfoam内有没有可以参考的MPI案例，求解答

0_1521016902124_TIM截图20180314163901.png

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

0_1514836555114_xmas-ornaments-erick-2000-small.png

挺有意思

每个人都能出名15分钟

--安迪 奥霍尔

摩拜、简书、一直播、债无忧，乍一看这些近3年成熟的新产品毫无关系，但是他们的共同点是背后的资本推手。以共享单车为例，甚至在10年前，北京就存在“共享单车”，当然叫做公共自行车更为恰当。那么为什么共享单车做到了亿万市场，公共自行车缺黯然没落。这和背后的运营主体有直接关系。

以摩拜为例，以企业作为运营主体的摩拜2015年就获得千万元人民币的投资，并在1年后完成了近10亿人民币的C轮融资。资本催熟了共享单车的运营形式，只有在有了资本之后，才能投入“车满为患”的量级。只有路边随手就是共享单车，才能实现公共自行车的初衷。

民间资本1年投入的这10个亿，等于我国在CAE行业10年的投入。

北快手，南抖音，通过农村包围城市，成就乡村网红，随即转型名流。相对于大V云集的微博，“去中心化”，“激发普通人的创造力”是抖音短视频以及快手短视频的核心概念。后者已经获得3.5亿美元的战略融资。在抖音上，只要你发的作品不拘一格，有创造力，即可以获得数十万的点击，并获得上万的粉丝。以前那个名人效应不复存在，在“去中心化”的背景下，人人都可以成为网红，进而变现造富。

从人人都可以装B的知乎，到资本井喷的人脸识别，以及你女朋友用的小咖秀，到硕士起薪30万+的AI ，这些被资本催熟的产品，并不适用于CFD，或许不适合今天中国的CFD行业。

行业的发展永远摆脱不了资本的渗入。CFD的发展，更离不开资本的投入。只有资本投入，才能成就行业。但我们要十指交叉，在未来，总会有若干有远见的民间资本家，为我国CAE行业注入资本的血液，成就中国的工业制造自主涅槃。

各位大佬好~向大家请教一个问题。

式（1）是大涡模拟法夫尔滤波操作以后的连续性方程
\begin{equation}
\frac{\partial{\bar{\rho}}}{t}+\frac{\partial{\bar{\rho}\tilde{u}_i}}{x_i}=0
\end{equation}

如果对$\tilde{u}_i$取时均值，即$<\tilde{u}_i>$, 这样得到的是雷诺时均值的速度吗？如果想要取法夫尔时均值，是不是需要做$<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>$的处理？P.S(这里用<φ>表示对φ取时间平均)

第二个问题与第一个问题相似：基于法夫尔时均的可解尺度的湍动能的求解是否可以按照式（2）不可压缩流体的求解方式求解呢？我感觉式（2）的解更像是雷诺时均的解，因为可解尺度的湍流里面也包含了密度脉动。
\begin{equation}
k_{resolved}=\frac{1}{2}(<\tilde{u}_i\tilde{u}_i>-<\tilde{u}_i>^2)
\end{equation}

还是说需要用式（3）继续用法夫尔时均进行计算呢？
\begin{equation}
k_{resolved}=\frac{1}{2}\big(<\bar{\rho}\tilde{u}_i\tilde{u}_i>-<\bar{\rho}>(<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>)^2\big)/<\bar{\rho}>
\end{equation}

@zhanghan

科研环境感觉还是国外好。

但是国外好山好水好无聊。没有狐朋狗友也没有烤串火锅。语言也不通。同时毕竟是亚洲人，发展中国家，很多对中国人有歧视。日本人地位比中国人要高。

#include "udf.h"
DEFINE_EC_RATE(user_ec_rate, f, fthread, r, V, current, didV, Eeq)
{
double alpha_a = r->alpha_a, alpha_c = r->alpha_c;
double io = r->io
double T = F_T(f,fthread);
double arg1, arg2;
cxboolean tafelmethod = r->tafelmethod;
int i;
double eta;
if (tafelmethod)
{
alpha_a = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_a * FARADAY_CONSTANT);
alpha_c = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_c * FARADAY_CONSTANT);
}
Eeq = r->Eeq;
eta = V - Eeq;
for(i = 0; i<r->n_reactants; i++)
if( ABS( r->exp_reactant[i] ) > SMALL_S )
{
int ni = r->reactant[i];
io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_reactant[i]);
}
for(i = 0; i<r->n_products; i++)
if( ABS( r->exp_product[i] ) > SMALL_S )
{
int ni = r->product[i];
io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_product[i]);
}
arg1 = FARADAY_CONSTANT / (UNIVERSAL_GAS_CONSTANTT);
arg2 = arg1eta;
current = io( exp( arg2alpha_a ) - exp( -arg2alpha_c ) );
didV = io( arg1alpha_aexp( arg2alpha_a ) + arg1alpha_cexp( -arg2alpha_c ) );
/* If multiple electrochemical reactions are used, you can define rate for each reaction
using the following if-statement /
/
if (STREQ(r->name, "reaction-1"))
{
...
}
else if (STREQ(r->name, "reaction-2"))
{
...
}
*/
}