View Issue Details
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0003518||OpenFOAM||Bug||public||2020-07-07 19:24||2020-07-07 19:24|
|Platform||Linux||OS||Scientific Linux||OS Version||6.3|
|Fixed in Version|
|Summary||0003518: Heat balance not closed in quarter of hollow cylinder with anisotropic thermal conductivity in cylindrical coordinate system.|
|Description||When calculating the steady state temperature distribution in a quarter of a hollow cylinder (end part) with an anisotropic thermal conductivity defined in a cylindrical coordinate system, connected to a cuboid (straight part) with anisotropic thermal conductivity defined in a Cartesian coordinate system, both with uniform volumetric heat generation, more heat is generated within the end part than physically possible (or given at the input).|
The cause is within the end part because, when changing the conductivity in the end part to isotropic or anisotropic with the same conductivity in the three directions (r, θ, z), the heat balance is closed. When using an anisotropic conductivity eg (1, 100, 1), probably the wrong kappa value is used when calculating the heat flux within the turbulentTemperatureCoupledBaffleMixed boundary condition.
The same issue was encountered before for the straight part, but changing the kappaMethod from solidThermo to directionalSolidThermo solved that issue for the straight part. For the end part however, this does not solve the issue.
|Steps To Reproduce||A test geometry was set up to reproduce the issue with three cases:|
- Test-case with anisotropic thermal conductivity → with the issue (Issue_Aniso_ani)
- Test-case with isotropic thermal conductivity → no issue (Issue_Aniso_iso1)
- Test-case with isotropic thermal conductivity, but defined as anisotropic with the same conductivity in the three directions → no issue (Issue_Aniso_iso3)
Open test case (Issue_Aniso_ani) and run the ‘Allrun’ file. Afterwards, run the ‘Allpostprocess’ file to check the heat balance:
- Within ‘Issue_Aniso_ani\postProcessing\end\wallHeatFlux\0\wallHeatFlux.dat’ or ‘Issue_Aniso_ani\postProcessing\straight\wallHeatFlux\0\ wallHeatFlux.dat’, it can be seen that the heat flux from ‘end’ to ‘straight’ is respectively 15.56157W or 15.53741W, which is 5.9284W higher than the input value 9.609W
- Within ‘Issue_Aniso_ani\postProcessing\straight\patchAverage(name=s1Conv,T)\0\ surfaceFieldValue.dat’, the area and average surface temperature of the cooled surfaces can be found. Calculating the heat flux with q=hA(T_avg-T_inf)=2000*0.001*( 317.9635-300)=35.927, which is 5.927W higher than the input value of 30W. This value is very close to the excess losses coming from the end part, which shows that the issue is in the end part.
A repetition of the previous steps for both other test-cases with isotropic thermal conductivity in the end part (Issue_Aniso_iso1 and Issue_Aniso_iso3), shows correct results (heat balance closed).
|Additional Information||This issue is related to issues:|
- 0002127: Volumetric heat source in solid with anisotropic conductivity
- 0002043: turbulentTemperatureCoupledBaffleMixed doesn't get heat fluxes balance in case of anisotropic conductivity
- 0002302: wallHeatFlux utility: wrong calculation of heatfluxes using anisotropic kappa
Some details of the test geometry: the geometry consists of a quarter of a hollow cylinder (end part) with:
- an anisotropic thermal conductivity (1,100,1) defined in a cylindrical coordinate system (r, θ, z)
- inner radius 10mm, outer radius 20mm and height 10mm
- heat generation of 9.6090W
connected to a cuboid (straight part) with:
- anisotropic thermal conductivity (1,1,100) defined in a Cartesian coordinate system (x,y,z)
- height and width of 10mm, length of 50mm
- heat generation of 20.3910W
- top and inner side cooled with convective boundary condition (h=2000W/m²K, T=300K)
All other surfaces are adiabatic.
- solver chtMultiRegionFoam
|Tags||No tags attached.|