[sasanghomi]

Hi Dear Foamers,

I am trying to understand calculating of motion continuity error.
I think first line in below code is creating an object without any dimension and I don't understand the meaning of [d/dt ( Object ) ] . Also I couldn't understand what do meshPhi(U) do exactly .

volScalarField motionContErr =
        fvc::ddt(dimensionedScalar("1", dimless, 1.0), mesh)
      - fvc::div(fvc::meshPhi(U));

Can anybody help me for understanding this code?
I appreciate any help from you.

Thanks and best regards,
Sasan.

[sharonyue]

Hi,

This is for space conservative law equation, check it out in Peric:

[thomasArk47]

Hello,

Sharonyue gives you good answer but maybe not so easy to understand if you are "new" to CFD discretization techniques If it is the case, maybe it is difficult to understand what means the time derivative of 1. Isn't it equalt to zero???

In fact, if you want to really understand the meshPhi(U) goal, you should read some courses on the topic "ALE technique" (for Arbitrary Lagrangian Eulerian) which allows the computation of flows on moving meshes. One difficulty in doing that is to respect a "hidden" conservation law which is called in the litterature the "GCL" (for Geometric Conservation Law). Roughly speaking, the time evolution of the mesh must be done in such a way that one respect a kind of compatibility condition( the so-called GCL) between the temporal evolution of the volume of a cell and the face flux resulting from the velocity associated with the mesh motion. This is precisely the meaning of the equation writtent by Sharonyue.