[willzyba]

All,

Within iterDyMFoam you'll find on line 128

gh = (g & mesh.C()) - ghRef;
ghf = (g & mesh.Cf()) - ghRef;

What does ghRef actually do?

We're created as an enhancement to the package waveFoam and waveDyMFoam. Simply merging the two. But in waveDyMFoam we have to write:

gh = (g & mesh.C());
ghf = (g & mesh.Cf());

Because in create fields we have to swap out

#include "readGravitationalAcceleration.H"
#include "readhRef.H"
#include "gh.H"

with

volScalarField gh("gh", g & (mesh.C() - referencePoint));
surfaceScalarField ghf("ghf", g & (mesh.Cf() - referencePoint));

to avoid various double assignment of variables in the waveFoam.

Obviously the solution is to sort though headers, but before I get that far I'd love to know the importance of ghRef.

[phmusi]

Hello Manoj and Agustino,

I am struggling at the same step.

I found this:
http://www.cfd-online.com/Forums/ope...c-meshing.html

And proceed to the opposite: in createFields.H
uncomment the two #include and comment the two solScalar and surfaceSaclar:

original:

//include "readhRef.H"
//include "gh.H"
#volScalarField gh("gh", g & (mesh.C() - referencePoint));
#surfaceScalarField ghf("ghf", g & (mesh.Cf() - referencePoint));

new:

#include "readhRef.H"
#include "gh.H"
//volScalarField gh("gh", g & (mesh.C() - referencePoint));
//surfaceScalarField ghf("ghf", g & (mesh.Cf() - referencePoint));

At least now it compiles and the floating run (crash after but might be another issue ^^)

The other solution I found is to copy-paste the definition of gh.H just before it is used in waveDyMFoam, but in removing the two last lines.

if (mesh.changing())
                {
                    Info<< "Execution time for mesh.update() = "
                        << runTime.elapsedCpuTime() - timeBeforeMeshUpdate
                        << " s" << endl;

        dimensionedScalar ghRef
            (
                mag(g.value()) > SMALL
              ? g & (cmptMag(g.value())/mag(g.value()))*hRef
              : dimensionedScalar("ghRef", g.dimensions()*dimLength, 0)
            );
            //volScalarField gh("gh", (g & mesh.C()) - ghRef);
            //surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);

                    gh = (g & mesh.C()) - ghRef;
                    ghf = (g & mesh.Cf()) - ghRef;
                }

AND by adding

#include "readhRef.H"

in the main (just after #include "readWaveProperties.H")

This way i am not altering the createFields.H

The simulation runs and then crash the same (hopefully meaning it does not come from here)

Cheers

Pierre-Henri