[vitors]

Dear all,

I've been browsing through some threads asking/explaining how to use fvOptions to add, for example, a heat source to a solid.

Unfortunately I found nothing (and my apologies for opening a new thread if it is there) and since this specific solver has changed a lot on the last versions of OpenFOAM, I need to ask few questions to the experienced foamers.

{
1    for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
2    {
3        fvScalarMatrix hEqn
4        (
5          - fvm::laplacian(betav*alpha, h, "laplacian(alpha,h)")
6         ==
7            fvOptions(rho, h)
8        );
9

10        hEqn.relax();
11

12        fvOptions.constrain(hEqn);
13

14        hEqn.solve();
15

16        fvOptions.correct(h);
17    }
}

This is the solveSolid.H file of chtMultiRegionSimpleFoam.
Lines 7, 12 and 16 are puzzling me.

I know (by reading other threads and OpenFOAM documentation) that in line 7 the fvOptions(rho, h) works as a source-term, assuming the user has set the fvOptions file correctly.

But I'm stucked here. So, is this a volScalarField based in density(???) and enthalpy?

How this impact in the thermophysical model?

Lines 12 and 16 are also a mystery. How fvOptions constrains que equation if I do not set a fvOptions dictionary previously?

Sorry for the bunch of doubts, but former versions of this solver were simpler and I'm struggling to understand the meaning of fvOptions from solver's point of view.

Thank you very much in advance for any clarification.

Vitor

[zfaraday]

I've been having the same doubts as you for some time. Searching for answers in the forum about how to set up a heat source using fvOptions with no luck.

Let's see if we finally see the light!

Regards,

Alex