[Roxxor91]

In the context of dynamicInkJetFvMesh, the motion of every points is given explicitly, as seen in its implementation of mesh.update:

bool Foam::dynamicInkJetFvMesh::update()
{
    scalar scalingFunction =
        0.5*
        (
            ::cos(constant::mathematical::twoPi*frequency_*time().value())
          - 1.0
        );

    Info<< "Mesh scaling. Time = " << time().value() << " scaling: "
        << scalingFunction << endl;

    pointField newPoints = stationaryPoints_;

    newPoints.replace
    (
        vector::X,
        stationaryPoints_.component(vector::X)*
        (
            1.0
          + pos0
            (
              - (stationaryPoints_.component(vector::X))
              - refPlaneX_
            )*amplitude_*scalingFunction
        )
    );

    fvMesh::movePoints(newPoints);

    lookupObjectRef<volVectorField>("U").correctBoundaryConditions();

    return true;
}

However, how faces and cells are updated ? More specifically, points are still belonging to the same faces, I can't figure out how face normals and cell volume are updated after point motion.

Moreover, I dont't understand how the method to get "mesh.phi" works. It is defined as

 const Foam::surfaceScalarField& Foam::fvMesh::phi() const
 {
     if (!phiPtr_)
     {
         FatalErrorInFunction
             << "mesh flux field does not exist, is the mesh actually moving?"
             << abort(FatalError);
     }
  
     // Set zero current time
     // mesh motion fluxes if the time has been incremented
     if (!time().subCycling() && phiPtr_->timeIndex() != time().timeIndex())
     {
         (*phiPtr_) = dimensionedScalar(dimVolume/dimTime, Zero);
     }
  
     phiPtr_->setOriented();
  
     return *phiPtr_;
 }

and it does not look to me that phiPtr_ points tothe desired surfaceScalarField..

Can you help me with that ?