[AnishPal]

I want to implement a time varying vol. heat generation in chtmultiregionFoam. I have done so following the previous threads but my temperature is not increasing. Can anyone please help me. And sometimes while running the solver it is saying selection mode not define although i have defined it. I think the codedsource is not being properly read by the solver. I am attaching the case here as well.

[AnishPal]

Can anyone please help me?

[seregaxvm]

In all examples heat source value is added to the heSource field:

heSource += 1e3*(time.value())*V;

[seregaxvm]

After investigation using method described in this post I've learned that additional heat source value should be negative and must be divided among mesh cells. So I came up with this line:

heSource -= (heatW / (heSource.size()));

where heatW is a heat generation rate in [W].
Make sure that heatW is a floating point variable, otherwise you'll get rounded integer division.

[trailer]

Originally Posted by seregaxvm

After investigation using method described in this post I've learned that additional heat source value should be negative and must be divided among mesh cells. So I came up with this line:

heSource -= (heatW / (heSource.size()));

where heatW is a heat generation rate in [W].
Make sure that heatW is a floating point variable, otherwise you'll get rounded integer division.

Please correct me if I am wrong, but to implement a source term with fvOptions I think you should have the following:

const scalarField& V = mesh_.V();
scalarField& S = eqn.source();

forAll(S, cellI)
{
     // Calculation of your source term for each cell in the domain
     mySourceTermForCelllI = Something



     S[cellI] = -V[cellI] * mySourceTermForCelllI;
}

Additionally, do you know why the fvOptions dictionary is sometimes found on constant folder and other times on the system folder?
OpenFOAM documentation says that the file can "live" in both folders [https://www.openfoam.com/documentati...voptions.html]. This must be one of the only dictionaries with this type of flexibility ...

[seregaxvm]

Yes, I've skipped boilerplate code. Additionally, mesh_.V() is a cell volume list. You'll need it only if you define your heat source in [W/m^3] units.

[trailer]

Originally Posted by seregaxvm

Yes, I've skipped boilerplate code. Additionally, mesh_.V() is a cell volume list. You'll need it only if you define your heat source in [W/m^3] units.

Sorry if this is too basic, but do you know the reason why I am multiplying by the cell volume and the heat source is being given in [W/m^3].

[seregaxvm]

One other note. codedSourse doesn't seem to work when limitTemperature is on.

Originally Posted by trailer

Sorry if this is too basic, but do you know the reason why I am multiplying by the cell volume and the heat source is being given in [W/m^3].

Sorry, didn't understand the question.