[JohnMB]

I am having difficulties with the proper region to region boundaries using chtMultregionfoam. for heat transfer



To test I made a very simple model with 3 regions, in a stack, with the top boundary held at 800 and the bottom a symmetry plane. The middle region(buildBlock) has a high thermal conductivity compared to the upper (heater) and lower(filledVolume).


What I see is that the heat flux between the top and middle is very low, the high conductivity region is acting as a thermal barrier!

I'm assuming that I have set up the boundaries incorrectly somehow.


This is from the changeDictionaryDict of the top layer

T
{
internalField uniform 300;

boundaryField
{
".*"
{
type zeroGradient;
value uniform 300;
}
heater_to_buildBlock
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
Tnbr T;
kappaMethod solidThermo;
value uniform 300;
}
maxY{
type fixedValue;
value uniform 800;
}

}
}

The only difference in the thermophysical properties between the regions is in the value of kappa, which s 5 for the upper an lower regions and 50 for the middle region.

thermoType
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

mixture
{
specie
{
molWeight 12;
}

transport
{
kappa 5;
}

thermodynamics
{
Hf 0;
Cp 500;
}

equationOfState
{
rho 8000;
}
}

Can anyone suggest what i am getting wrong?

[JohnMB]

I think i may be misinterpreting the results.


and I cant figure out how to delete the thread

[pete20r2]

I'll just add that I have found turbulentTemperatureCoupledBaffleMixed to cause slow convergence.
I modified the laplacianFoam solver to take DT (thermal diffusivity) in as a field and I assign varying properties to different cellZones via setFields.
This way the whole problem is one mesh (for all solid) and no explicit coupling is required.

[JohnMB]

Originally Posted by pete20r2

I'll just add that I have found turbulentTemperatureCoupledBaffleMixed to cause slow convergence.
I modified the laplacianFoam solver to take DT (thermal diffusivity) in as a field and I assign varying properties to different cellZones via setFields.
This way the whole problem is one mesh (for all solid) and no explicit coupling is required.

Thanks that sounds like a useful method.

I may do that to validate some aspects of the model, however I ultimately need anisotropic heat conductivity and will probably also need some fluid regions for the final model, so I think I need to keep working with multiregion to get to the end goal.Unless you know a way of including anisotropic conductivity in LaplacianFoam.

[pete20r2]

Give me a couple of days, I'm pretty sure I can do it with tensors. I would like to make a cht style solver that uses the proper thermo classes, which make things like heat addition easy but I haven't got that far. I'm working on some problems like this at work now so I'll post here if I develop anything. Do you need arbitrary anisotropy or is it all orthogonal, i.e (condX condY condZ)?

[JohnMB]

Originally Posted by pete20r2

Give me a couple of days, I'm pretty sure I can do it with tensors. I would like to make a cht style solver that uses the proper thermo classes, which make things like heat addition easy but I haven't got that far. I'm working on some problems like this at work now so I'll post here if I develop anything. Do you need arbitrary anisotropy or is it all orthogonal, i.e (condX condY condZ)?

Orthogonal.
Thanks.