CFD中文网

给你一个搜各种书的网站，不需外网。一直用，偶尔令人失望。
https://vdoc.pub/

<a href="#" class="icon icon-home"> 东岳流体|OpenFOAM|CFD <img src="_static/logo.svg" class="logo" alt="Logo"/> </a>

img一行代码缺失，需要加进去

@李东岳 后面我测试了粒径更小的颗粒，发现同等条件下更小粒径的颗粒不会出现上述情况，应当是MPPICInterFoam要求网格体积和颗粒体积的比值比DPMFoam的更高

http://blog.sina.com.cn/s/blog_599d8faa010121yc.html
参考一下流沙大哥的

getdata取点；
给原作者发个邮件问一下；
这个图是excel 做的吧

你可以用excel 也可以用origin

@bestucan 好嘞感谢大佬的回复:146:

可以加载多个

@cfd-china 写文章，做实验。。

我们都知道，速度可以做Helmholtz分解，用速度势的梯度和流函数旋度叠加的形式，求解二维粘性可压缩翼型绕流中的速度势和流函数相当于求两个泊松方程，但是需要给出边界条件，一般情况下，在速度场已知的情况下，即涡量场和胀量场已知，该如何分别给速度势和流函数的物面和远场边界条件？

@李东岳 找到了 https://vimeo.com/148670249

你查一下这方面的实验论文吧？
壁面速度低，容易造成管壁聚集

@xuanze 3374769460@qq.com

就是两种液体接触的界面部分就没有办法加密啊？

不就是界面（表面）么？

@chpjz0391
非常感谢，已更新。:sunglasses:

Cd 变化，物理上就是这样吧

@liqiwen
Hello, 欢迎加入！

流体理论中层流中的边界层的确定使用的是普朗特的理论，而ansys flunt中加入边界层只需要把流体域边界处的网格加密，无需其他设定就可以计算了，请问普朗特的理论在cfd中是怎么引入的；

普朗特什么理论？虽然边界层分为层流边界层和湍流边界层。但是在计算中主要调查湍流边界层。CFD软件把网格加密或者加密到合适的高度，是因为边界层内速度和网格距离存在一定的关系。这个关系成为壁面法则。要实施壁面法则，壁面处网格需要满足距离上的一些要求，因此需要加密网格。湍流神秘y+

湍流和层流的边界层理论在flunt中的设定是否一致？

湍流可以使用壁面函数。层流边界层只是速度变化比较大。不能使用壁面函数。因此需要细化网格。

百度 google也搜不到啥。我就写在这里给大家提个醒吧。
之前尝试从华夏银行APP汇入法国账户1万欧。
应该是账户写错了。
折腾几天系统自动退回来了。
扣了90欧。
华夏银行这面没有扣钱。
就是在欧洲不知道是哪个环节扣了90欧。
大家知道就好。以后要跨境汇款心里有个数。

微信图片_20240214134432.jpg

@User0221 你具体做的什么东西我不太清楚，不好评价，但是个人觉得和网格细化没有多大的关系。

#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"))
{
...
}
*/
}