CFD中文网

谢谢各位的回复！

你最好还是找版本在2.x以上的例子吧

好吧，错误信息第一行显示在make/option文件里有问题，你看一下

你可以试试看。

@桥下卧龙 我说的是在整个文件夹里找，openfoam/...

文件搜索啊，linux和windows系统都有这个功能吧

你看错误信息就知道了啊，从 In file include from /opt/.... 那里就有提示了，在那个路径之下的h文件的第几行出现的代码，到最后提示找不到“Scalar.H”，所以你要找到这个文件在哪。

那个是1.3版本的吧，看错误信息是找不到h文件。也许是include方面的问题?

近来学习openFOAM，发现有fvoptions，但是在网上查不到详细的说明和较为粗浅的例子进行学习。请问

1. fvOptions是实现运行时添加源项的，但是如何实现？我看到的所有例子中 type 应该都是内建的，定义的文件在哪里？
2. 有没有可能通过fvOptions实现动态的源项？例如弹性的多孔介质什么的，孔的大小受到邻近流体流速场、压力场等的调节什么的，这样就构成互相影响的关系。如果有现成的类似例子就更好了，因为这和我的课题方向更类似一些。

好的，谢谢两位的指点！

看了给的帖子，以及 https://openfoamwiki.net/index.php/IcoFoam 的说明，虽然细节还不懂，但是有了粗糙的理解。
由于phiHbyA 对应 在face上的插值，对内部的面而言无所谓，但是在边界上就需要调整，保持一致性。那么phi 也是在面上的插值，由于它对应是 U 的插值，U 已经定义了边界条件，所以就不用调整了。 或者 adjustPhi(phi, U, p) 也是可以写的，但是只是在边界上重复一遍没什么用吧。
不知道这样的理解对不对。

谢谢了~~~，有没有发帖格式的说明？

东岳大神写了关于icoFoam solver的一些解释，不过还有一些细节问题不太清楚，特来请教 （OpenFOAM 3.0.1）：

1. 这段代码：

while (piso.correct())
        {
            volScalarField rAU(1.0/UEqn.A());
            .....
        },  

这里的UEqn.A和后面的UEqn.H是自动根据U的新值进行刷新的吗？
如果是固定不变的，那就没有必要写在while循环里；如果是自动更新的，那么不考虑什么非正交循环或特殊的算法，只是最简单的不断根据旧值求新值，solver是不是可以写成 下面这种？

    fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu, U)
        );
        while (......) {
        solve(UEqn == .....)
        }

但是在代码里写的是

            surfaceScalarField phiHbyA
            (
                "phiHbyA",
                (fvc::interpolate(HbyA) & mesh.Sf())
              + fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)
            );

            adjustPhi(phiHbyA, U, p);

为什么要加上 fvc::interpolate(rAU)*fvc::ddtCorr(U, phi) 这一项？以及在后面adjustPhi(phiHbyA, U, p) ？ 因为只按照图片中的方程显示，只要 fvc::interpolate(HbyA) & mesh.Sf()) 就够了。