|
|
model_scheme.PNG (54,986 bytes) comparison.png (46,910 bytes) comparison2.PNG (139,656 bytes) |
|
|
The handling of anisotropic thermal conductivity is incomplete in OpenFOAM-10 and before, it assumes alignment with the boundary. I have recently completed a rewrite of the entire solid energy transport in OpenFOAM-dev which might/should resolve this issue: commit 01d2c6ec52bf0c5056d06d26245f0f46af1eefe3 Author: Henry Weller <http://cfd.direct> Date: Mon Oct 24 19:01:52 2022 +0100 solidThermophysicalTransportModel::anisotropic: Added heat-flux correction qCorr to handle the wall heat-flux generated by the tangential components of the temperature gradient when the thermal conductivity tensor Kappa does not align with the boundary. This is an uncommon situation but for cases where it is important it must be handled correctly; previously this term was ignored. For efficiency the temperature gradient field is automatically cached so that it is not evaluated for every patch for which Kappa is not aligned with the boundary, and the correction is evaluated and applied only for those patches. The heat-flux correction qCorr is used in the coupled and external heat-flux boundary conditions coupledTemperature and externalTemperature. |
|
|
You will need to check very carefully the setup of the coordinate system, it looks like from your results that the first component of kappa is not in the radial direction. In the OpenFOAM-dev implementation you can write out the Kappa tensor field and study it. |
|
|
I could not unpack the test-case, RAR is not a standard Linux compression format. Try your case with (28 1.3 28) and (28 28 1.3) |
|
|
Sorry, I have uploaded the case in .tar.gz format now. Regarding your suggestions: - I have already tried running the case with (28 1.3 28) and (28 28 1.3) using OpenFOAM 10. Tried the same approach again with OpenFOAM-dev and uploaded the comparison of the results. The behaviour in OpenFOAM-dev is different - all the results were the same. Also, I have tried to use (1.3 0.1 0.1) with OpenFOAM-dev and got max temp equal to 308.44 K, which is close to the value obtained using isotropic properties. - I have run the same case using freshly installed OpenFOAM-dev - the inner radius temperature decreased to 300.279 K. - I don't know how to get a kappa tensor using OpenFOAM-dev but I managed to obtain a KappaLocal tensor from OpenFOAM 10 using the debugger. The values were KappaLocal[0] = 9.73, [1] = 12.41, [2] = 0, [3] = 19.57, [4] = 0, [5] = 28. I presume that this is just an example value for 1 point, but I hope it might help. comparison_10_vs_dev.PNG (185,279 bytes) OpenFOAM-dev_1.3-0.1-0.1.PNG (75,429 bytes) |
|
|
I set the case up in dev, reorganising as a solid only case running with foamRun -solver solid setting the coordinate system in thermophysicalTransport model anisotropic; coordinateSystem { type cartesian; origin (0 0 0); coordinateRotation { type cylindrical; axis (0 0 1); } } and used a more accurate grad scheme because of the horrible tet-mesh: gradSchemes { default leastSquares; } and I get the same results for constIsoSolid with kappa 1.3 constAnIsoSolid with kappa (28 28 28) constAnIsoSolid with kappa (1.3 28 28) There is no bug, this is a user-support request, if you need further assistance you will need to arrange for a support contract. |
|
|
User error and support request. |
- Anonymous
