[nikoscham]

Hi community!

I would like to add a vector field (with constant value) in the momentum conservation equation.

I have the following equation form in the solver:

   fvVectorMatrix UEqn
    (
        fvm::ddt(U) + fvm::div(phi, U)
      + fvm::Sp(C*pow(1-lf, 2)/(rhoRef*(pow(lf, 3)+q)), U)
      + MRF.DDt(U)
      + turbulence->divDevReff(U)
     ==
        fvOptions(U)
     );

which is part of the meltingFoam solver (see here melting problem: looking for appropriate solvers). This solver is described in the paper https://www.researchgate.net/publica...th_experiments

I would like to subtract a constant vector field (Up) from the velocity vector U.

I changed the code as follows:

   fvVectorMatrix UEqn
    (
        fvm::ddt(U) + fvm::div(phi, U)
      + fvm::Sp(C*pow(1-lf, 2)/(rhoRef*(pow(lf, 3)+q)), U - Up) //modified
      + MRF.DDt(U)
      + turbulence->divDevReff(U)
     ==
        fvOptions(U)
     );

Moreover I added the following in the createFields.H

Info<< "Reading field Up\n" << endl;
volVectorField Up
(
    IOobject
    (
        "Up",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    ),
    mesh
 );

By compiling the solver, however I get errors. Part of the error log is the following:

/opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvmSup.C:142:1: note:   template argument deduction/substitution failed:
In file included from meltingPullFoam.C:97:0:
UEqn.H:12:63: note:   ‘Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >’ is not derived from ‘const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>’
       + fvm::Sp(C*pow(1-lf, 2)/(rhoRef*(pow(lf, 3)+q)), U - Up) //added
                                                               ^
In file included from /opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvmSup.H:165:0,
                 from /opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvm.H:49,
                 from /opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvCFD.H:10,
                 from meltingPullFoam.C:48:
/opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvmSup.C:156:1: note: candidate: template<class Type> Foam::tmp<Foam::fvMatrix<Type> > Foam::fvm::Sp(const dimensionedScalar&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&)
  Foam::fvm::Sp

It seems that the operation (U-Up) is not accepted. . .

Do you have any proposals or tips?

Thank you in advance!


Nikos