|
|
T_of6.png (22,052 bytes) chtPlate.png (2,899 bytes) T_of7.png (22,271 bytes) T_of7mod.png (22,276 bytes) |
|
|
> I would have thought the enthalpy departure function for rhoConst should be zero. Why? When I evaluate the departure function for rhoConst I get p/rho. There are various ways it is presented in terms of the various forms of thermodynamic state variables but a convenient form for this purpose is on page 59 of http://web.nchu.edu.tw/pweb/users/cmchang/lesson/10174.pdf or specifically for liquids on page 17/18 of https://nptel.ac.in/content/storage2/courses/103101004/downloads/chapter-5.pdf If you evaluate the departure function from that what do you get? Incidentally what solution do you get when using internal energy in OpenFOAM-6 and 7: thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Do you get the same solution for both the internal energy and enthalpy forms in both versions of OpenFOAM? |
|
|
The problem is that the pressure in the solid is non-uniform and so the pressure contribution to the enthalpy is non-uniform, specifically the bulk pressure in the solid is 1e5 but it is specified as 0 on the heater-liquid boundary: heater_to_liquid { type calculated; value uniform 0; } when I change this to heater_to_liquid { type calculated; value uniform 100000; } the case gives the results you expected. |
|
|
Thank you very much! I can confirm that changing the entry as per your suggestion returns the good agreement seen with version 6, see attached T_of7_corr_hTL_p.png The values for the patches heater_to_liquid and liquid_to_heater are set (I think) during splitMeshRegions. I had assumed that the patches would inherit values from (say) 0/p, in which all the patches, apart from frontAndBack, have calculated values: #include "./system/parameters" dimensions [1 -1 -2 0 0 0 0]; internalField uniform $pInf; boundaryField { ".*" { type calculated; value $internalField; } frontAndBack { type empty; } } After splitMeshRegions is run, the patch value settings are modified using changeDictionary. For system/liquid/changeDictionaryDict, the values for p and p_rgh are set explicitly, whereas they are not modified in system/heater/changeDictionaryDict. Then, the entry for 0/liquid/p.liquid_to_heater is correctly set as 100000, whereas (as you noted) the value for 0/heater/p.heater_to_liquid is 0 Some other points: - when the corrected case - ie with the value for heater_to_liquid correctly set - is run using sensibleInternalEnergy, the agreement is good: the values are somewhat higher than previously, but the results are withiin ~ 0.4% of theory, see attached T_of7_sIE.png - when above run is repeated in version 6, there is an error: [3] --> FOAM FATAL ERROR: [3] Unknown rhoReactionThermo type thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Valid rhoReactionThermo types are: heRhoThermo homogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo homogeneousMixture const hConst perfectGas specie sensibleEnthalpy heRhoThermo homogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo homogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy heRhoThermo inhomogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo inhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy heRhoThermo inhomogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo inhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst adiabaticPerfectFluid specie sensibleInternalEnergy heRhoThermo multiComponentMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst incompressiblePerfectGas specie sensibleInternalEnergy heRhoThermo multiComponentMixture const hConst perfectFluid specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst perfectFluid specie sensibleInternalEnergy heRhoThermo multiComponentMixture const hConst perfectGas specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst perfectGas specie sensibleInternalEnergy heRhoThermo multiComponentMixture const hConst rhoConst specie sensibleEnthalpy heRhoThermo multiComponentMixture const hConst rhoConst specie sensibleInternalEnergy heRhoThermo multiComponentMixture polynomial hPolynomial icoPolynomial specie sensibleEnthalpy heRhoThermo multiComponentMixture polynomial hPolynomial icoPolynomial specie sensibleInternalEnergy heRhoThermo multiComponentMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo multiComponentMixture sutherland janaf incompressiblePerfectGas specie sensibleInternalEnergy heRhoThermo multiComponentMixture sutherland janaf perfectGas specie sensibleEnthalpy heRhoThermo multiComponentMixture sutherland janaf perfectGas specie sensibleInternalEnergy heRhoThermo pureMixture const hConst Boussinesq specie sensibleEnthalpy heRhoThermo pureMixture const hConst Boussinesq specie sensibleInternalEnergy heRhoThermo pureMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy heRhoThermo pureMixture const hConst adiabaticPerfectFluid specie sensibleInternalEnergy heRhoThermo pureMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy heRhoThermo pureMixture const hConst incompressiblePerfectGas specie sensibleInternalEnergy heRhoThermo pureMixture [1] T_of7_corr_hTL_p.png (21,982 bytes) T_of7_sIE.png (22,059 bytes) |
|
|
User error in the specification of the solid pressure. |
- Anonymous
