COMSOL Forums: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/
Most recent forum messagesMon, 22 Oct 2012 11:43:59 +0000COMSOL Forums: Help modeling supercooling effect in water to ice transitionhttp://www.comsol.com/shared/images/logos/comsol_logo.gif
http://www.comsol.com/community/forums/general/thread/32034/
Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p88667
In the discussion about supercooling and freezing it is my feeling that a more apropiate model can be obtained by modelling the change of enthalpy and not the temperature. During cooling you continuously remove enthalpy and the temperature follows. The temperature is a function of the enthalpy content and the percentage of ice. The latter lags behind during super cooling.<br />
<br />
If you make a plot of H versus T for an equilibrium process, and take the derivative dH/dT, then you can smear the latent heat of fusion peak a bit and put it in a Cp(T) function. That is a 'normal' trick.<br />
But with supercooling this will not work. See the picture <a href="http://www.esveld.org/Comsol/HvsT%20and%20TvsH.png" title="www.esveld.org/Comsol/HvsT%20and%20TvsH.png">www.esveld.org/Comsol/HvsT%20and%20TvsH.png</a><br />
dH/dT is not a unique function of T. You can define dT/dH as function of H. Or even directly T als function of H(H)<br />
<br />
Thus you need to solve a system dH/dt = PDE (T) with T=F(H) or dT/dt = f(H). Will that work in Comsol?<br />
<br />
More physically correct would be to model the two phases ice and water and their transition reactions with entalpy change, like you can do for drying. This is a very stiff system as the local melt / freeze reactionrates are very fast. The system does not deviate from equilibrium, except when the crystallisation rate is low because of lack of nuclii.<br />
<br />
<br />
ErikMon, 22 Oct 2012 11:43:59 +000032034.1350906239.88667Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p88477
Hi<br />
<br />
a few suggestions:<br />
<br />
Use rather the analytical functions to define your full variable and material expressions, as then, with a click, you can check them, and make a plot to see how they behave, then call your variables or material parameters to use the analytical functions you have defined.<br />
<br />
A variable useful to help you set up your thermal models is the alpa heat diffusivity, as when you have a very different heat diffusivity its worth to have a fine mesh at the interfaces to resolve correctly the high gradients you develop here.<br />
<br />
define a model variable alpha = material.k11/material.rho/material.Cp<br />
<br />
in fact there is a relation to respect between average mesh size "h" alpha and the smallest time stepping dt<br />
<br />
then when you are ready to run your solver, start running only the initial conditions and set up a surface plot for alpha, at the common boundaries with large alpha difference, check that you have a dense(r) mesh<br />
<br />
Then only run your model<br />
<br />
In your case you consider your fluid is a "solid" no convection, which is often a perfectly valuable simplification hypothesis, <br />
<br />
for such diffusion problems often a power of 2 time stepping is good, but as you have changes in temperature and slope variations you must ensure that these are also well followed. Then you must set the solver time stepping node to "intermediate, or strict to ensure that COMSOL does not jump over one of your time transitions<br />
<br />
And while running long time solving, turn on the plot while solving FOR all Steps to see how it behaves, time is shorter like that and you learn how physics behave<br />
<br />
And have you checked that your total integrated Cp variation corresponds to the latent heat of melting/solidification ?<br />
<br />
--<br />
Good luck<br />
IvarFri, 19 Oct 2012 10:59:21 +000032034.1350644361.88477Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p88474
Hello again everybody! Here is a very simple model of what <br />
I am trying to model, a mass of water and a small plate of aluminum <br />
exposed to a temperature. Functions for the temperature profile and<br />
for the characteristics<br />
<br />
I am open to suggestions!!<br />
<br />
Thank you for your time!Fri, 19 Oct 2012 09:18:23 +000032034.1350638303.88474Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87658
Taking inspiration from the continuous casting model, the Cp function simulates, or is supposed to simulate,<br />
the latent heat of the water, 333kJ/kg, smoothed over a couple of degrees (I played with the dT value a bit from <br />
0.5 to 3 K)<br />
<br />
<br />
I will provide you with a simple model as soon as possible so that you can have a look.<br />
<br />
Tue, 09 Oct 2012 04:15:35 +000032034.1349756135.87658Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87621
Hello, Nikolaos,<br />
<br />
Sorry, but I still don't make an idea of your model and why you are unconfortable with the graphics T vs. t. If you could provide a simple (better 1D or 2D) version of the model, perhaps we would be able to help you more.<br />
<br />
Nevertheless, I have inspected your Cp function, and it has a sudden and strong peak at T about 273[K] (Cp reaches 1.92075e5 at T = 0[degC]!!!). Is that what you intend? How do you distinguish water from ice? Why don't you take into account the difference in k (thermal conductivity) and rho (density) between both phases? Could you go without any geometry, and thus use only an ODE?<br />
<br />
Jesus.Mon, 08 Oct 2012 13:56:46 +000032034.1349704606.87621Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87614
Hello Jesus and thank you for your reply. The water in this simulation acts only as a heat storage medium.<br />
The "problem" is that in during the climate chamber simulations the water is standing still and thus allowed<br />
to be supercooled before turning into ice. I might have to repeat the simulations trying to force the phase transition<br />
(by jiggling the container maybe?) around the 0°C point in order to get a smoother behavior. <br />
<br />
The water is cooled and warmed by a "Temperature" boundary condition applied on a plate <br />
underneath the container.<br />
<br />
By phase coexistence do you mean a "mushy" state? I don't think that is necessary at the moment.<br />
<br />
I have considered making the model a bit more complex to account for movement in the liquid but again<br />
I think that is too complex for what I am trying to simulate.<br />
<br />
Thank you for your time and advice!!Mon, 08 Oct 2012 09:53:23 +000032034.1349690003.87614Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87601
Hi, Nikolaos,<br />
<br />
Ivar's help has been very valuable, as always. And perhaps you can receive more help if you describe your model basically, or attach a simplified (if possible 2D) version of it: which physics do you use? (you mention "Heat transfer in solids", but you you try to simulate heat in a liquid?), which boundary conditions? (how do you heat the water), etc.<br />
Have you considered the coexistence of two phases?<br />
<br />
And about the general question: "how to simulate supercooling in water", have you taken into account that the phase trasition would be a function not only of the temperature, but also of the velocity magnitude? (and hence, you should simulate the fluid with a moving mesh, I think).<br />
<br />
Bye,<br />
<br />
Jesus.Mon, 08 Oct 2012 07:02:12 +000032034.1349679732.87601Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87596
I have concluded now that modeling the supercooling effect itself, isn't that useful for<br />
the scope of the simulation and an "averaged" behavior is more fitting.<br />
<br />
The question remains but this time for philosophical reasons.Mon, 08 Oct 2012 05:36:09 +000032034.1349674569.87596Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87543
Hi<br />
<br />
In COMSOL most entry fields accept full equations (including function calls) so you can define variable densties, temperatures etc. But depending on how you set these up, you might make your model very non-linear, hence the solver might start to have issues to converge.<br />
<br />
But be sure you are defining "physical" cases, else the results might be bizarre, or non physical too ;)<br />
<br />
--<br />
Good luck<br />
IvarSat, 06 Oct 2012 18:11:59 +000032034.1349547119.87543Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87486
Can I make the internal temperature of the mass of water rise that way?Fri, 05 Oct 2012 09:17:42 +000032034.1349428662.87486Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87481
Hi<br />
<br />
in V4 you can use the predefined "STEP()" which is a variant but perhaps with less parameters to tweak<br />
<br />
--<br />
Good luck<br />
IvarFri, 05 Oct 2012 08:55:13 +000032034.1349427313.87481Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87468
Thank you for your reply, I have used flc2sh to produce a smooth "step" but it<br />
gets smoothed or averaged over instead of having that relatively steep step. <br />
I don't mind so much at the moment because the final result is "close enough"<br />
but, you know how it is, wanting to get it precisely right ;)I tweaked the solver to <br />
take 0.5 sec steps around the step area but it did not do any difference.<br />
<br />
I am starting to believe that if the water is treated as a solid, the moment of<br />
crystallization isn't simulated correctly.<br />
<br />
Again, thank you for your time Mr. Kjelberg<br />
Fri, 05 Oct 2012 04:30:34 +000032034.1349411434.87468Re: Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87457
Hi<br />
<br />
I would suspect some difficulties to make it a true "step" but if you use a combination with the v4 "step" function and use smooth corners, I believe you could manage, but you would need to find a way to tweak your solver to detect the step and take "small" solver steps to resolve the steep gradients<br />
<br />
--<br />
Good luck<br />
IvarFri, 05 Oct 2012 03:18:47 +000032034.1349407127.87457Help modeling supercooling effect in water to ice transition
http://www.comsol.com/community/forums/general/thread/32034/#p87416
Hello everyone, relatively new to COMSOL here, I am trying to create a model<br />
describing a simulation that took place in a climate chamber. I have attached a <br />
graph describing the temperature on the water container surface. I want help modeling<br />
that step-up near the 1600s mark, how should I "make the water warmer" before it starts<br />
cooling again? I have looked into the continuous casting model for guidance and I have<br />
used some methods described in it but I am stuck...<br />
The physics I use is heat transfer in solids and the custom water I have<br />
has a Cp = 4200 - (2120*(T<T_freeze)) +exp(-(T-T_freeze)^2/((dT)^2))/sqrt(pi*(dT)^2)*333000<br />
<br />
where dT=1[K] and T_freeze = 0[degC]<br />
<br />
Any help would be welcome!<br />
<br />
Thank you!Thu, 04 Oct 2012 11:33:57 +000032034.1349350437.87416