[quarkz]

Hi,

I've trying to modify the Turbulent kinetic energy equation.

Basically, for some reasons, I need to modify from e.g.:

fvm::laplacian(A,k)

to

fvm::laplacian(A/sqrt(k),sqrt(k))

where the full eqn is:

tmp<fvScalarMatrix> kEqn
    (
        fvm::ddt(alpha, rho, k)
     ...
      - fvm::laplacian(A/sqrt(k),sqrt(k))
...

The errors I got are:

error: no matching function for call to ‘laplacian(Foam::tmp<Foam::GeometricField<dou ble, Foam::fvPatchField, Foam::volMesh> >, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >)’

/opt/OpenFOAM/OpenFOAM-v1912/src/finiteVolume/lnInclude/fvmLaplacian.C:47:1: note: template argument deduction/substitution failed:
In file included from ../turbulenceModels/lnInclude/kOmegaSSTF.H:194:0,

If I change to explicit:

fvc::laplacian(A/sqrt(k),sqrt(k))

Then there is no problem.

But if I want it to be implicit, is there any way?

Thanks in advance.

[HPE]

to use fvm for sqrt(k), the governing equation should br sqrt(k) rather than k.

Is this a source term which can be put on the right hand side?

[quarkz]

Sorry, I think I made a mistake. I am solving for k, so I'm actually computing fvc::laplacian(A/sqrt(k),sqrt(k)) on the RHS as explicit.

Thanks anyway!