关于buoyangBoussinesqSimpleFoam添加新的边界条件

LiuZhuang

@东岳我是一个OF的小白丁，如何在controlDict中添加新的方程，以前没听说过，有什么资料可供参考？

李东岳

在你的controlDict下面添加

functions
{
    #includeFunc scalarTransport
}

然后在controlDict同级目录下创建字典文件scalarTransport，输入下面的信息


type            scalarTransport;
libs            ("libsolverFunctionObjects.so");

field           T;
schemesField    T;
D               1e-09;//扩散率

然后去你的0文件夹下给T的边界条件，并设置T的fvScheme参数

LiuZhuang

@东岳理解了很多，东岳老师，第二张图片中的D扩散率我是根据nut/Sct+nu/Sc进行添加，这就要引入一个新的边界条件。参考于alphat中的alphatJayatillekeWallFunction中的nut/Prt+nu/Pr,增加一个nut/Sct+nu/Sc。我该怎么做？

LiuZhuang

@东岳按照您的设置，我加入这个方程之后和原来的buoyangBoussinesqSimpleFoam中的温度方程互不影响么？其T表达的意义可以是浓度的求解感觉么？然后扩散率的设置变成了nut/Prt+nu/Pr TO nut/Sct+nu/Sc的形式么？这几个迷惑之处，希望东岳老师指教，拜谢！！！

LiuZhuang

@东岳，东岳老师，我做这个求解器的浓度方程添加进去了，程序也能算了。浓度的效果还不错。但是我现在遇到的问题是温度的变化也呈现出一种类似于浓度扩散的变化，请问是为什么？我需要把哪些代码发上去给您看看？？？

李东岳

温度的变化也呈现出一种类似于浓度扩散的变化

把你的扩散系数降低10倍100倍试试

LiuZhuang

@东岳原先的温度变化比较合理，自从我加入了自己编写的浓度方程和新的nut/Sct+nu/Sc的边界条件后，温度场就呈现出了浓度扩散那种一圈一圈的变化

李东岳

加入了自己编写的浓度方程

把代码复制到这里，并高亮看一下

LiuZhuang

这个是按照温度方程改变的浓度方程

替代文字

LiuZhuang

@东岳 alphatJayatillekeWallFunction;这是0文件alphat文件夹下的边界条件，我仿照这个设置了A边界条件，其意义是将nut/Prt+nu/Pr的热扩散系数更换为nut/Sct+nu/Sc的浓度扩散系数给出浓度扩散的边界条件。不知我的思路是否正确

李东岳

我仿照这个设置了A边界条件，

浓度扩散系数也应该用壁面函数么？

LiuZhuang

@东岳那个我添加成了，可能是我考虑多了，加和不加区别不大

CFD中文网