关于buoyangBoussinesqSimpleFoam添加新的边界条件



  • 0_1543624781197_6cb4ae83-65e3-4b61-a111-b000bccf7e68-image.png
    旧的边界条件位置如图
    打开alphatJayatillekeWallFunction.C文件,可以看到其分别对nut/Prt,nut/E(nut/E可能不对)做出设置求解方法。现在我想在内部添加一个nut/Sct的计算(Sct也是和Prt一样的常数),那么我该对这个文件的.C和.H文件怎样修改,如何编译?
    小弟初学乍到,描述不周,欢迎质疑讨论。在此拜谢,不胜感激!!!
    PS:我用的OpenFOAM版本为2.4.0



  • 因为nut/Sct这个数的设置体现了浓度扩散系数ν的系数,也就是说我在buoyantBoussinesqSimpleFoam中又基于温度方程增加了一个浓度方程。如果不去改动tutorials中0文件夹下的alphat文件(alphat里的边界条件就是alphatJayatillekeWallFunction,作用是设置nut/Sct+nu/Sc),只是复制alphat一个作为nut/Sct的边界条件文件夹,算出来的温度场呈浓度场样貌。故提出1L所问,请问我的问法或者说思路是否有毛病?



  • 也就是说我在buoyantBoussinesqSimpleFoam中又基于温度方程增加了一个浓度方程

    等下,为何不在controlDict文件下直接添加一个浓度方程?



  • @东岳 那我去试试,但是添加浓度方程的话一样需要将alphat里的文件中nut/Prt+nu/Pr更改成nut/Sct+nu/Sc,这个新的Sc常数我在transportproprities已经设置好了。但是在0文件夹下的alphatJayatillekeWallFunction对应做更改怎么实现?恳请东岳老师赐教



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



  • 在你的controlDict下面添加

    functions
    {
        #includeFunc scalarTransport
    }
    

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

    
    type            scalarTransport;
    libs            ("libsolverFunctionObjects.so");
    
    field           T;
    schemesField    T;
    D               1e-09;//扩散率
    

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



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



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



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



  • @LiuZhuang

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

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



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



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

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



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

    替代文字



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



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

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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2