|View Issue Details|
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0002907||OpenFOAM||Contribution||public||2018-04-21 10:39||2018-05-01 16:08|
|Target Version||Fixed in Version||dev|
|Summary||0002907: a proposal of fvOptions solidificationMeltingSource for not isotherm phase change|
|Description||current fvOptions solidificationMeltingSource allows only isotherm phase change |
namely for pure substance or eutectic composition mixture.
Following Voller's indication, it should be possible to take into account
the not isotherm miscible systems phase change and
the not isotherm portion of not eutectic systems phase change
introducing an at least linear relation between melt fraction alpha1 and temperature.
Focusing on cast iron portion of Fe-C phase diagram it seems to me to be fine the results obteined considering Tmelt_-> Tsol_idus-eutectic, a new Tliq_uidus > Tsol and a new alpha1e_ max eutectic melt fraction as proposed in attached files. I noted some instability on temperature results for alpha1e >= 0.9.
|Steps To Reproduce||1- compile fvOptions substituting attached files in fvOptions/sources/derived/solidificationMeltingSource/|
2- case definition requires a new syntax for fvOptions dictionary:
i- use Tsol instead Tmelt as previously to define melting temperature in isotherm phase change (for pure substance or eutectic mixture -> Tsol=Tliq);
ii- optionally define new Tliq > Tsol to consider liquidus temperature in not isotherm phase change (for miscible mixture), where previous defined Tsol defines solidus temperature;
iii- optionally define also alpha1e to consider max eutectic melt fraction (that should be the percentage of solvent phase changed from initial to eutectic liquid concentration) in partially isotherm (at Tsol=Teutectic) and not
isotherm (from Tsol=Teutectic to Tliq) phase change (for solid not miscible mixture) [alpha1e=0 -> pure substance; alpha1e=1 -> eutectic mixture that is strictely not permitted[.
|Additional Information||the attached implementation has only eutectic linear capability (lever rule) in the hope that it will be considered; |
i belive a better fit to real thermal behavior should be possible introducing:
+- a negative exponent n for the factor (alpha1[celli]-alpha1e)/(1-alpha1e) for
eutectic power capability (n=1->linear, n=1/(1-k)->Sheil)
++- a control to consider a not isotherm phase change arrest at alpha1e proportional to the cooling rate (undercooling) due to solidification heat release (recalescence of irreversibility), and probably related to the exponent n too.
Thanks for the contribution.
Integrated by commit bad26bee71be1b56c0cc41e4bedbf491b5832567
|2018-04-21 10:39||lorenzotrevisan||New Issue|
|2018-04-21 10:39||lorenzotrevisan||File Added: solidificationMeltingSource.tar.gz|
|2018-04-21 10:39||lorenzotrevisan||Tag Attached: solidificationMelting|
|2018-05-01 16:08||henry||Assigned To||=> henry|
|2018-05-01 16:08||henry||Status||new => resolved|
|2018-05-01 16:08||henry||Resolution||open => fixed|
|2018-05-01 16:08||henry||Fixed in Version||=> dev|
|2018-05-01 16:08||henry||Note Added: 0009532|