CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Ice melting Simulation

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By ghorrocks

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 28, 2010, 17:22
Default Ice melting Simulation
  #1
Member
 
Pratik Mehta
Join Date: Mar 2009
Posts: 32
Rep Power: 17
pratikmehta is on a distinguished road
Hello Friends,

I am trying to perform a simulation for ice melting

Here is the Problem definition:

I have a aluminium hollow pipe ( 1mm thickness ) . Through this pipe water at 80 C is flowing at constant flow rate . Re number is about 20K
The ice deposition over this pipe is about 2mm thick .

My quest is to know how much time does it take for the ice at -5 C gets converted to water at 5 C . I just need to know when does water reach 5C.


What I have done so far:

I have got geometry and meshed it .I have the fluid domain ( water )and solid domain ( Ice ) .I am using CFX 12.1v. Material in fluid domain is water & material in Ice domain is user created material as CHT solid & thermodynamic state as solid.

My doubts,,,
1) How should I treat the interface ??. Should I keep it as conservative interface flux as the boundary condition and should I create a thin wall of 1 mm aluminium ,so as to save my meshes to aviod the extra physical 1mm thickness of the pipe.

2) Variables like Density, Thermal conductivty & Sp heat capacity as all function of temperature. Regarding the thermal conductivity & density I have used step function CEL expression that will suitably feed in the right number at the existence of ice or water. My biggest worry is to specify Cp variation as the ice start to heat up from -5C .

I need to know how should I specify Cp that will consider the Cp of Ice from -5C to -1C approx 2kJ/KgK , then it must consider Cp as 166kJ/KgK frmo -1C to +1 C and later from 1C onwards Cp as 4kJ/KgK .

I tried to use the step function in order to get this variation but could not get it correct.

So how would you suggest me to go forward .


Thanks for your suggestions

Best regards
Pratik
pratikmehta is offline   Reply With Quote

Old   September 28, 2010, 19:55
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Can't your problem be done as an analytical solution, heat input, assumed HTC, latent heat etc to get how much ice melts? What will CFD add to it?
ghorrocks is offline   Reply With Quote

Old   September 28, 2010, 19:57
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Actually, it probably won't be an analytical solution but sounds like a simple ODE so can be easily done in matlab or a hand written ODE solver.
abhay.uiet likes this.
ghorrocks is offline   Reply With Quote

Old   September 29, 2010, 04:13
Default
  #4
Member
 
Pratik Mehta
Join Date: Mar 2009
Posts: 32
Rep Power: 17
pratikmehta is on a distinguished road
Hi Glenn,

Ya actually this is what could be the solution but currently i am just tryin to get the melting process to work. Actually problem definition is that the water flows through the pipe at varied temperature over period of 20 mins ok , ,, so over the time of 20 mins of constant flow rate with varied temperature of water flowing through the pipe I need to melt certain substance called Ad blue for SCR ( having its material properties ).

So for a start I thought of setting up a simple case to understand the complications before i taken the above mentioned problem addition.

Well all I would think of that this simulation is pretty the same as the Deicing / Defrosting on a windshield which is a regular CFD problem.

Also my question is ,, to use either a step function for Cp or use a user rountine fortan but I have heard before that fortan sub routine cannot be run properly on a windows platform ???, well i am just speculating i am not sure.

thanks
pratikmehta is offline   Reply With Quote

Old   September 29, 2010, 09:26
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
There are many complexities in modelling freezing or melting. You still have not explained why CFD is needed here, why won't a simple ODE solution be good enough? You would have the answer after a few hours work rather than a few months by CFD.

If you insist on doing it by CFD then I would talk to CFX support and ask for some demo files on freezing or melting. I am pretty sure some example files are around.

Fortran runs fine on Windows, you just need to have the correct fortran compiler. Look in the documentation for details.
ghorrocks is offline   Reply With Quote

Old   September 29, 2010, 12:39
Default
  #6
Member
 
Pratik Mehta
Join Date: Mar 2009
Posts: 32
Rep Power: 17
pratikmehta is on a distinguished road
Hi Glenn,


Melting the substance under consideration while supplying the water which is responsible for heat transfer whose temperature at the inlet is sent in at different temperatures ( starting with 20 C upto 70 C within the time duration of 20 mins ) .Hence I thought doing this in CFD would be helpful as also we have no details abt the local heat transfer cofficient and assuming a constant heat transfer coefficient or by empherical value would cause errors. Over all geometry is also complex in reality. Also we are interested to know whether how much melting in liquid fraction has been obtained.

PS: Dont you think that deicing is very similar issue and people in industries do go for such CHT CFD simulation in order to know the fraction of liquid to solid and get more info.

let me know if you still think if this can be done using ODEs.

Well the CFX support have not given me any documentation as of yet.


Thanks for your reply

Best regards
Pratik
pratikmehta is offline   Reply With Quote

Old   September 29, 2010, 19:57
Default
  #7
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
we have no details abt the local heat transfer cofficient
Then do a simple single phase model to just get the HTC and use that in an ODE approach.

Quote:
Over all geometry is also complex in reality.
That may mean the ODE approach is intractable depending on the geometry and how the stuff melts. Can't say one way or the other without further details.

Quote:
PS: Dont you think that deicing is very similar issue and people in industries do go for such CHT CFD simulation in order to know the fraction of liquid to solid and get more info.
It is similar to deicing but many deicing simulations assume the ice is a thin layer on the surface and when it turns to liquid it disappears. This may not be appropriate for what you are doing. Can't be sure based on the sketchy details you have provided so far.
ghorrocks is offline   Reply With Quote

Old   September 30, 2010, 03:35
Default
  #8
Member
 
Pratik Mehta
Join Date: Mar 2009
Posts: 32
Rep Power: 17
pratikmehta is on a distinguished road
Hi ,

Ya may be in this case , I can assume all the liquid disappears as well as long as that is indicating the formation of liquid state. But never the less my query is to implement the correct Cp so as to observe this process in CFX. Assuming an average Cp would insert errors in overall melting process so I was looking for help for the same.

PS: I would be glad if anyone has done even the deicing simulation to direct me in how have they managed to put in Cp or inserted the latent heat in the data as I see that CFX does not allow to put this information directly. With conductivity and densitx changed can be managed with Step function.


Thank again for your kind help

Best Regards
Pratik
pratikmehta is offline   Reply With Quote

Old   September 30, 2010, 06:09
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
I do not understand what you are trying to do. Can you post some drawings and a description of how it works?
ghorrocks is offline   Reply With Quote

Old   February 29, 2012, 05:39
Default Freezing of an ice cream
  #10
New Member
 
KEYUR PATEL
Join Date: Feb 2012
Posts: 5
Rep Power: 14
10mmet21 is on a distinguished road
Dear all,

I am Working On an interesting M. Tech Project. Description is discussed bellow.

It is a Two phase problem in which convection takes place through air flowing outside and conduction in the ice cream.

1 litter rectangular Ice cream packet is placed in air blast freezer.
Air is flowing at the velocity of 2.5 m/s. from one side of rectangular duct and leave at the other end of the freezer. It flows over the Rectangular ice cream box. The Temperature of an air is (-30 C)

Initial Temperature of an ice cream is (-4 C)
initial freezing temperature is (-5.6 C)
specific heat above freezing point is 3.22 kJ/Kg K and bellow freezing Point is 2.74Kj/kg K
latent heat of fusion is 204 kj/kg
density 550 kg/cubic meter
conductivity of unfrozen ice cream is 0.518 w/m K and frozen ice cream is 1.8203 w/m k

I have to find that what time will it take to reach the center temperature of an ice cream to (-18 C)

Air properties are available. heat transfer coefficient h=15 w/m2*K

conductivity of an ice cream box which is made up of paper can be taken as 0.05 w/m k. Thickness of wraper is 0.5 mm.

If any one has related tutorial of freezing then please send it to me or upload on the forum.

I wanted to plot graph of Temperature vs time as an Output.

Please Help me.
10mmet21 is offline   Reply With Quote

Old   February 29, 2012, 07:02
Default
  #11
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Does the liquid ice cream move at all?

Are there any moulds, pipes, formers or anything else to hold things in place? I assume they are made of metal or plastic or something. Or is the liquid ice cream in a paper box?
ghorrocks is offline   Reply With Quote

Old   February 29, 2012, 08:25
Default
  #12
New Member
 
KEYUR PATEL
Join Date: Feb 2012
Posts: 5
Rep Power: 14
10mmet21 is on a distinguished road
Thank you very much for taking an interest in my problem.

The ice cream is stationary in 1 liter rectangular family pack. (It is Kept in Rectangular Box made of paper card board... You have Already seen it)

This Box is placed in air blast freezer.

One more thing is Ice cream is in paste form (Not a liquid Not a solid) so i am confused that what should i consider it initially. finally it will became solid after freezing.

Another thing is Ice cream properties like specific heat, density and conductivity are Changing from paste form to solid form.

if i use solidification and melting model software ask me about solidus temperature and liquidious temperature. what value should i take. the freezing temperature is (-5.6 degree centigrade)

so cold air of (-30 C) at velocity 2.5 m/s (turbulent flow) is blown on (-4 C) temperature ice cream and air continuous to flow till the ice cream center temperature reaches to (-18 C)

please help me if you can do some thing.

waiting for your positive reply.
10mmet21 is offline   Reply With Quote

Old   February 29, 2012, 18:14
Default
  #13
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Sounds like the simplest approach is to model the ice cream as a solid with a variable specific heat as a function of temperature, and include in it a jump which represents the latent heat of freezing. Likewise the thermal conductivity is simply a function of temperature.

The paper is a thermal resistance on the outside of the ice cream solid.

Then you can do a CHT simulation to couple it to the fluid flow. But I suspect the coupling will be weak and you can probably do a fluid simulation to get heat transfer coefficients over the ice cream and then just do an isolated fluid flow model and isolated solid temperature model.
ghorrocks is offline   Reply With Quote

Old   February 29, 2012, 22:14
Default
  #14
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,162
Rep Power: 23
evcelica is on a distinguished road
Looks like you were given the convection coefficient of 15W/m^2*K?
If so forget CFD and just do a transient solid model and apply the convective boundary condition and temperature to the outside. Just increase the Cp over a small range like you were already doing to represent the latent heat.


A solid model like this shouldn't take more than 30 minutes to setup, mesh, and solve.
evcelica is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
simulation of phase change material melting rosol FLUENT 4 December 8, 2015 08:56
ice melting simulation pratikmehta ANSYS 3 April 17, 2013 10:52
Continuous vs interrupted simulation sega OpenFOAM Running, Solving & CFD 4 November 3, 2008 15:29
can CFX use for simulation in ICE? hurin CFX 1 December 30, 2005 09:30
3-D Contaminant Dispersal Simulation Apple L S Chan Main CFD Forum 1 December 23, 1998 11:06


All times are GMT -4. The time now is 00:16.