[Gennaro]

Dear all,

I'm solving, in my customized solver based on CHTMutliregionFoam, the following equation:

d/dt = alphaEff * laplacian (T)

which is a classic heat transfer equation, where alphaEff (Effective thermal diffusivity) is a scalar field. I'm solving it in the following way:

    fvScalarMatrix EEqn
    (
   fvm::ddt (T)
-  fvm::laplacian(alphaEff, T)
    );

Now I would like to solve

rho*Cp*d/dt = divergence (Keff * gradient(T))

where rho is the density, Cp is the specific heat and Keff is the effective thermal conductivity, i.e they are all scalar fields.

How can I write it in OpenFOAM?

Thanks in advance

Kind regards

Gennaro

[mohsen cheraghi]

Hi,

divergence of x is written like: dev(x)

cheers