CFD中文网

写jou.
设置并行

比如CSDN 上提供的解决方案

/define/user-defined/compiled-functions/compile
        libraryName
        yes
        sourcefilename.c

/define/user-defined/compiled-functions/load
       libraryName

;#编译并载入UDF

/file read-case
       yourcasename.cas
       /file read-data
       yourdatname.dat

;#读入调试好的case、data文件

/define/user-defined/function-hooks/adjust "youradjustname::libraryName"""

;#hook像诸如init adjust这些需要hook的udf,不知道为何除了udf里定义的DEFINE_ADJUST(youradjustname,d)外，还会有一个名字""为空的，并且三个"要连写，两者之间不能有空格，不知道代表什么

/solve/patch

phase-2

12

mp

1

;#我计算的是一个液面晃荡的问题，所以是两相流，要Patch一下，选择phase-2作为我们的domain,然后我要patch的zone的编号是12，然后要空格一下进入下一步,然后选择mp即体积分数，令其为1
       /file/auto-save/root-name yourrootname
       /file/auto-save/data-frequency 10000

/file/auto-save/case-frequency 10000

;#自动保存data，case,每10000个时间步或是10000迭代步
       /solve/iterate
       200000

;#如果是稳态问题，迭代200000次

/solve/set/time-step 0.0001
       /solve/dual-time-iterate 
       50000
       200

;#如果是非稳态问题，时间步为0.0001，时间步数50000，每步迭代200次
       /file/write-case-data
       yourfinalname.out
       exit
       yes

;#最后导出最后的case与data文件

然后最后提交任务可以用lsf、pbs啥的，但是我不会，不过直接用命令提交貌似也麻烦不了多少。

新建一个文件new

然后文件里面写

node9

node9

node9

....

node9

这样共n列代表你使用节点9，里面的n个核。

nohup fluent 3d -t12 -cnf=new -g -i yourjouname.jou >&youroutname.out&

其实-t12表示一共12个核，-cnf=new使用node9的12核，然后运行过程输出到youroutname.out中，这是使用一个节点的命令，使用多个节点只需对命令作些许修改如下 只需加入-pethernet 和 -ssh

nohup fluent 3ddp -pethernet -cnf=myhosts -t36 -ssh -g -i f30.jou>&out.out &

作者：aleigus
来源：CSDN
原文：https://blog.csdn.net/aleigus/article/details/77161181
版权声明：本文为博主原创文章，转载请附上博文链接！

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/HrvojeJasak/DiscretisationBestPractice.pdf

请自己看手册。。。
祝早日脱离苦海

This questions is not clear. What is 液力变矩器?

@李东岳 车灯热分析。

@李东岳 谢谢东岳大神。祝论坛越来越好。

就是下面这个
$$
\frac{\p T_{ij}}{\p x_i\p x_j}
$$
其实$T_{ij}=u_iu_j$，所以演变为
$$
\frac{\p u_iu_j}{\p x_i\p x_j}
$$
这个公式是可以写成$\nabla\cdot(\nabla\cdot(\bfU\bfU))$这种形式，$\frac{\p T_{ij}}{\p x_i\p x_j}$好像写不出来 :134:

大佬们看看？

@李东岳 在 有关侵入边界法。 中说：

Hi，

OpenFOAM-extend里面有IBM方法，Extend和官方的OpenFOAM侧重不同，几个我了解到的不同主要是IBM，粘弹流体，动网格表面跟踪。

还有可压缩流动高分辨率格式。https://github.com/Unofficial-Extend-Project-Mirror/Unofficial-Extend-Project-Mirror.github.io/wiki

@simthere 好久没有画结构网格了，回归下，练练手，暂时还能画出来，哈哈0_1503487106341_l1.png 0_1503487116733_l2.png

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。

硕士毕业找到了一些研究所的工作，还是拒了选择读博。因为自己总觉得还有很多东西值得深入地学习。读博是个难得的提升自我的机会，比较enjoy独处的时候可以静下心来思考。生活也不算惨淡，有funding可以不愁吃喝。

@xpqiu 十分感谢！我会仔细根据您提的意见学习的，谢谢！

@李东岳 谢谢版主:laughing:

0_1472638385877_捕获.JPG
2_1472638455925_捕获4.JPG 1_1472638455924_捕获3.JPG 0_1472638455924_捕获2.JPG

Turbulence
CEFRC Combustion Summer School
Prof. Dr.-Ing. Heinz Pitsch
2014

List 13 desolate
荒凉的 despair n. vi. desperate adj.
a feeling that you have no hope at all
drive sb. to despair
对某人感到绝望 despair of sb.
My teacher began to despair of me dispatch
formal 为了特定的目的send某人或某物到某地
Dispatch several officers to Beijing despise
蔑视
I despise all cowards detain
拘留 deteriorate
恶化，变坏
cause the leather to deteriorate deviate
偏离
deviate from the truth
deviate from the original plan dictate dictation
口述，使听写
The boss dictate the plan to us
She is dictating a plan to her secretary diminish
reduce discern
to notice or understand something by thinking about it carefully
Officials were keen to discern how much public support there was
I would discern what his heart and soul were like discourse
演讲
I felt at home in this kind of discourse
a discourse on art discreet
小心谨慎的
I stood back at a discreet distance
Reporters asked some discreet questions discrepancy
a difference between two amounts, details, reports etc that should be the same
Police found discrepancies in the two men's reports

很有帮助。

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

有个方向，叫混沌之类的。详细可以了解下蝴蝶效应和非线性的一些内容。

我导收了一个湘潭大学的搅拌器的英文sci。实验+cfd，用openfoam搞的。我针对算例设置这一段提出了1个不痛不痒的问题，不能一个不提啊。目前已返回我导。估计明后天就返回编辑吧？哈哈。我导没有抵制，不过提了2页问题...anyway，问题不大。:sunglasses:

z 和 x 用得最顺。

其他的还在熟悉中。