cht solid temperature distribution
Summary
Solid temperature distribution in steady cht solver with mixture models pureZoneMixture and pureMixture leads nonphysical results. The same case by simulated with separated regions give a different results than single region case.
In pureZoneMixture case temperature distribution look that heat goes from lower to higher temperature. In pureMixture case heat flux does not match between regions. Temperature distribution looks logical but temperature levels are much higher that single region case.
Steps to reproduce
Please find attached cases. By looking heat flux data from log files and temperature distribution is possible notice difference between mixture models.
Example case
What is the current bug behaviour?
Results are unreliable
What is the expected correct behavior?
In pureZoneMixture case temperature distribution can be "fix" by setting the same Cp values to all materials. I think that in solid solver is used alpha instead of kappa in temperature equation.
Relevant logs and/or images
T fixed wall0 to 300 and wall1 to 400 K. Heat source 10 W in middle of the model. pureZoneMixture heat flux:
min/max/integ(wall0) = -1322.838, -1322.838, -132.2838
min/max/integ(wall1) = 1222.838, 1222.838, 122.2838difference between in and out 10 W and it is OK. Respectively pureMixture case
min/max/integ(wall0) = -114.41282, -114.41282, -11.441282
min/max/integ(z1_to_z2) = 594.6436, 594.6436, 59.46436
min/max/integ(z2_to_z1) = -594.6436, -594.6436, -59.46436
min/max/integ(z2_to_z3) = 2408.3388, 2408.3388, 240.83388
min/max/integ(wall1) = -15889.687, -15889.687, -1588.9687
min/max/integ(z3_to_z2) = -2408.3388, -2408.3388, -240.83388wall0 and wall1 difference is much higher (1588 - 11) than 10 W.
Environment information
- OpenFOAM version : v2212
- Operating system : AlmaLinux 8
- Hardware info : Intel Core 8
- Compiler : gcc version 8.5.0 20210514 (Red Hat 8.5.0-18) (GCC)
Possible fixes
[chtCases.zip](/uploads/0b33599299701063e04e68c42a50fb47/chtCases.zip