关于压力场边界更新的实现机制的疑问

李东岳

correctBoundaryConditions这个函数的作用非常简单。比如你通过下面的代码计算速度场：

U=Ua+Ub;

在这种情况下并没有计算U的边界条件，因此你需要

U.correctBoundaryConditions()

但是如果你通过UEqn.solve()的方式计算，solve()函数里面自动调用了correctBoundaryConditions。

多相流求解器里面也一样，比如

// Correct p_rgh for consistency with p and the updated densities
    p_rgh = p - rho*gh;
    p_rgh.correctBoundaryConditions();

是因为前面存在赋值 p_rgh = p - rho*gh;

因此，我的个人习惯是，如果是solve()出来的，不需要调用，如果是=赋值出来的，都需要correctBoundaryConditions。

千里之行

@李东岳哦，明白了，谢谢东岳老师，我再去好好看看相关的代码。那 "calculated" 类型的边界条件对应的类是 calculatedFvPatchField ，但是并没有重新定义虚函数 evaluate，肯定是基类 fvPatchField 的 evaluate 函数被调用了，可是那个 evaluate 函数并没有做什么计算，所以即使调用了场对象的 correctBoundaryConditions 函数，在 calculated 类型的边界也没有发生计算，这样理解对吗？
另外，经常会遇到这样的边界条件：

inlet
{
    type        fixedValue;
    value       $internalField;
}

这里的 $internalField 是什么意思呢，我认为是这块边界的值取邻近内场的值，但是一直也没有找到代码的依据，并且这样的边界条件和 zeroGradient 又有什么区别呢？

dzw05

@千里之行这是一个宏，意思是该处的value采用的是“internalField”这个变量的值。你在这个文件里开头的部分，会有类似internalfield uniform 100这样的语句。

千里之行

@dzw05 哦，懂了，原来如此，十分感谢！！

浪迹天大

@东岳
我也想了解关于 correctBoundaryConditions 使用的知识。
但是我发现等于号操作符也更新了边界场呀：

template<class Type, template<class> class PatchField, class GeoMesh>
void Foam::GeometricField<Type, PatchField, GeoMesh>::operator=
(
    const GeometricField<Type, PatchField, GeoMesh>& gf
)
{
    if (this == &gf)
    {
        FatalErrorInFunction
            << "attempted assignment to self"
            << abort(FatalError);
    }

    checkField(*this, gf, "=");

    // Only assign field contents not ID

    ref() = gf();
    boundaryFieldRef() = gf.boundaryField();
}

李东岳

执行下面的代码后，

U = Ua+Ub;

Ua+Ub的边界条件为calculated，然后U也为calculated，若需要更新U的边界条件，需要继续update

类似讨论： http://www.cfd-china.com/topic/1620

浪迹天大

就是说，场进行加减乘除操作之后，边界场虽然也会更新数据，但是此时的值可能违反它本身边界条件的约束，所以需要 correct 约束一下，是这样吗？

李东岳

是这个意思

evensun

@李东岳
老师，请问对于定义在体心上的物理量，这里correct的边界场，是指边界上的一层网格吗

李东岳

@evensun 是boundary的patch上的值，不是边界第一层网格的值

CFD中文网

关于压力场边界更新的实现机制的疑问