fsi solver problem

版本为 foam-extend3.1

fluidStructureInterface 库下的 fluidStructureInterface.C 文件

void Foam::fluidStructureInterface::calcCurrentSolidZonePoints() const
{
    // Find global face zones
    if (currentSolidZonePointsPtr_)
    {
        FatalErrorIn
        (
            "void fluidStructureInterface::"
            "calcCurrentSolidZonePoints() const"
        )
            << "Current solid zone points alarady exist"
                << abort(FatalError);
    }

   currentSolidZonePointsPtr_ =
        new vectorField(stress().currentFaceZonePoints(solidZoneIndex()));
}

该函数首先判断固体网格节点是否存在,如果存在则报错,如果不存在则进行计算赋值。
最后的赋值语句不太清楚意思,globe face zones是全部固体网格的意思把。

2 posts were split to a new topic: boundaryMesh()[fluid2PatchIndex_].faceCentres() 什么意思

这段的大概意思是:

如果 currentSolidZonePointsPtr_ 存在,则报错;否则,对其进行赋值,其初始值是将 stress() 里面的值 copy 过来。