[linch]

Hi Foamers,

in compressible solvers continuity is solved for density:

solve(fvm::ddt(rho) + fvc::div(phi));

with mass flux phi.

But how can I solve continuity for velocity?

I already tried to solve:

    fvScalarMatrix UEqn
    (
        fvc::ddt(rho)
      + fvm::div(rho,U)
    );

And the compiler returned an error:

error: no matching function for call to ‘div(Foam::volScalarField&, Foam::volVectorField&)’

Then I tried: this one:

fvScalarMatrix UEqn
{
   fvc::ddt(rho)
+ fvm::div(fvc::interpolate(rho),U)
}

but then I got a type mismatch error:

error: no match for ‘operator+’ in ‘Foam::fvc::ddt(const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = double]() + Foam::fvm::div(const Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = Foam::Vector<double>](((Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&)(& U)))’

    fvScalarMatrix UEqn
    (
        fvc::ddt(rho)
      + fvm::div(rho,phi)
    );

->

error: no matching function for call to ‘div(Foam::volScalarField&, Foam::surfaceScalarField&)’

fvScalarMatrix UEqn
    (
        fvc::ddt(rho)
      + fvm::div(fvc::interpolate(rho),phi)
    );

->

error: no matching function for call to ‘div(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >, Foam::surfaceScalarField&)’

fvVectorMatrix UEqn
    (
        fvc::ddt(rho)
      + fvm::div(fvc::interpolate(rho),U)
    );

->

error: no match for ‘operator+’ in ‘Foam::fvc::ddt(const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = double]() + Foam::fvm::div(const Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = Foam::Vector<double>](((Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&)(& U)))’

I'm looking forward for an expert-answer

Best,