Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    关于laplacian((1/A(U)),p)的疑惑与icoFoam.dep文件的疑惑

    OpenFOAM
    3
    7
    2197
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • D
      D.Benjamin last edited by

      今天我看了一个博客《OpenFOAM程序开发初步》,该博客作者向大家演示了将icoFoam修改成myicoFoam并给出了应用,在此,我遇到了两个问题,想向大家请教:

      1. 修改后的myicoFoam添加了温度方程,并用其求解cavity算例。如下图所示,我对laplacian((1/A(U)),p)离散格式中的A(U)有疑惑。我看到的东岳老师关于icoFoam中的解析以及源代码中的icoFoam.C中为下图所示,所以离散格式不应该为laplacian(rAU,p)嘛?
        0_1537606256306_4.PNG 0_1537606265373_6.PNG

      2. 博客中说道,修改求解器,需要删除依赖文件icoFoam.dep,如下图所说,进入icoFoam文件夹,删除icoFoam.dep。但是,我发现icoFoam文件夹中只包含icoFoam.C、createField.H和Make,并没有icoFoam.dep。我猜这有可能是OF的版本问题,也有可能根本没有icoFoam.dep文件,因此,我的问题是:我们在修改求解器的时候是否还需要删除该依赖文件,如果需要,它在哪个文件夹中。期待各位的回复,谢谢!
        0_1537606605985_5.PNG

      OpenFOAM初学者,希望和大家共同交流

      I W 2 Replies Last reply Reply Quote
      • I
        Izumi @D.Benjamin last edited by

        @d-benjamin

        我们在修改求解器的时候是否还需要删除该依赖文件,如果需要,它在哪个文件夹中。

        我使用of4、5,修改求解器时,直接复制原求解器文件夹、重命名、修改求解器及make文件,然后直接编译就行了

        D 1 Reply Last reply Reply Quote
        • W
          wwzhao 教授 @D.Benjamin last edited by

          @d-benjamin

          1. 离散格式指定的关键词与其名称有关。

          如果在声明变量的时候指定名称:

          volScalarField rAU(1.0/UEqn.A(), "rAU");
          

          那么可以在字典文件中用以下语句指定离散格式

          laplacian(rAU,p) Gauss linear corrected
          

          在声明rAU时不指定名称:

          volScalarField rAU(1.0/UEqn.A());
          

          那么rAU的名称通过量纲以及符号运算得到,1.0/UEqn.A()的名称为(1|A(U)),因此相应的离散格式应该写作:

          laplacian((1|A(U)),p) Gauss linear corrected
          
          1. 博客的例子OpenFOAM版本太老。OpenFOAM 2.x以前生成的dep依赖文件和源文件同目录,后续版本已经将这些文件移动到其他目录。
          D 2 Replies Last reply Reply Quote
          • D
            D.Benjamin @wwzhao last edited by

            @wwzhao 谢谢你的回复和帮助,让我受益良多,中秋节快乐!

            OpenFOAM初学者,希望和大家共同交流

            1 Reply Last reply Reply Quote
            • D
              D.Benjamin @Izumi last edited by

              @izumi 谢谢大佬的帮助,中秋节快乐

              OpenFOAM初学者,希望和大家共同交流

              1 Reply Last reply Reply Quote
              • D
                D.Benjamin @wwzhao last edited by

                @wwzhao 教授您好,我还想追问一个问题,我想请问一下UEqn.A()与UEqn.H()各代表着什么意思?A()与H()函数各能取得方程的哪个部分?期待您的回复,祝好!

                0_1537668280145_8.PNG
                0_1537668287241_9.PNG

                OpenFOAM初学者,希望和大家共同交流

                1 Reply Last reply Reply Quote
                • W
                  wwzhao 教授 last edited by

                  A是系数矩阵的对角部分,H是H操作符,是用于速度压力解耦的。

                  1 Reply Last reply Reply Quote
                  • First post
                    Last post

                  CFD中文网 | 东岳流体 | 京ICP备15017992号-2