CFD Online Discussion Forums - chtMultiRegionFoam solid to solid boundaries

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - chtMultiRegionFoam solid to solid boundaries (https://www.cfd-online.com/Forums/openfoam-solving/210860-chtmultiregionfoam-solid-solid-boundaries.html)

chtMultiRegionFoam solid to solid boundaries

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

Quote:

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.

Quote:

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?

I think i may be misinterpreting the results.
:(

and I cant figure out how to delete the thread

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.

Quote:

Originally Posted by pete20r2 (Post 714414)

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.

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)?

Quote:

Originally Posted by pete20r2 (Post 714481)

Orthogonal.
Thanks.