View Issue Details
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0003515||OpenFOAM||Contribution||public||2020-07-01 15:23||2020-07-01 15:56|
|Status||closed||Resolution||no change required|
|Fixed in Version|
|Summary||0003515: Support request about a proposal of fvOptions solidificationMeltingSource for irreversible phase change|
|Description||current fvOptions solidificationMeltingSource allows only reversible phase change namely at equilibrium temperatures without irreversible solidification undercooling.|
Thermal analysis highlights an undercooling as activation energy to get the solidification started, and a trend to restore equilibrium temperature.
Focusing on Sn-Pb alloy characteristic hypo, hyper and eutectic behavior it seems to me to be fine the results obtained considering appropriate calibration coefficients irrM (to control undercooling magnitude) and irrR (to control reheating rate), such that irrM is proportional to time average cooling rate and irrR is function of volume average melt fraction.
During development I was not able to consider the heating/cooling rate (dotT) in a smart way (as fvc::ddt(T) should allow), so I am primarily asking for support about it: please let see comment at row 151 and my workaround from row 200 of meltingSolidificationSource.C.
|Steps To Reproduce||1- compile fvOptions substituting attached files in fvOptions/sources/derived/solidificationMeltingSource/|
2- case definition requires a new additional syntax for fvOptions dictionary in order to consider solidification irreversibility:
i- optionally use irrMsol and/or irrMliq to calibrate undercooling magnitude respectively for eutectic and primary solidification;
ii- optionally use irrRsol and/or irrRliq to calibrate reheating(recalescence) rate respectively for eutectic and primary solidification;
iii- the attached implementation become unstable for heavy undercooling due to initial dotT calculation, so use Tstart in fvOption dictionary in order to delay the irreversible solidification calculation start.
|Additional Information|| +- I aspect Lirr < Lrev = L (latent heat)|
++- I’m testing a model to control nucleation delay
+++- I aspect irrM is a characteristic time proportional to: (dM^2)/(lambda/Cp*rho) | M=Volume/Area
|Tags||No tags attached.|
solidificationMeltingSource-v7-recalescence.tar.gz (6,811 bytes)
User support request, you will need to arrange a support contract to cover this.
|2020-07-01 15:23||lorenzotrevisan||New Issue|
|2020-07-01 15:23||lorenzotrevisan||File Added: solidificationMeltingSource-v7-recalescence.tar.gz|
|2020-07-01 15:56||henry||Assigned To||=> henry|
|2020-07-01 15:56||henry||Status||new => closed|
|2020-07-01 15:56||henry||Resolution||open => no change required|
|2020-07-01 15:56||henry||Note Added: 0011412|