CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

solidification and melting with VOF: how to sink solid PCM as it melts?

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

Like Tree5Likes
  • 1 Post By sakil2k3
  • 1 Post By msaeedsadeghi
  • 2 Post By Mahboobe365
  • 1 Post By thermal energy

Reply
 
LinkBack Thread Tools Display Modes
Old   January 24, 2013, 23:27
Talking solidification and melting with VOF: how to sink solid PCM as it melts?
  #1
New Member
 
mahmood
Join Date: Sep 2012
Posts: 14
Rep Power: 4
sakil2k3 is on a distinguished road
Hi all, I have an enclosure spherical glass that is contain with phase change material (PCM) that is in its solid state initially. Heat will be applied above the melting temperature of the PCM around the spherical glass. During the melting process of the PCM there will be 2 phases present in the enclosure, the liquid and solid state. The solid PCM will drop to the bottom due to the gravitation force.

A 2d mesh was created and exported out from Gambit.

In Fluent, I used Solidification and Melting function in my simulation with VOF. All Material properties of the PCM and glass have set accordingly. The Gravity was also enable with Y=-9.81 in the Operating Condition and the Variable-Density Parameters was also set to an operating density of PCM, that is 777kg/m^3. The Body force weighted under Discretization in Solution Control was also selected.

After iterating for 200 sec, the values for the Liquid fraction was calculated but the simulation did not show the solid PCM sink to the bottom of the enclosure.

How can I make the solid PCM to sink to the bottom as it melts?
Thanks.
thermal energy likes this.
sakil2k3 is offline   Reply With Quote

Old   January 26, 2013, 00:57
Default
  #2
Senior Member
 
SSL
Join Date: Oct 2012
Posts: 227
Rep Power: 5
msaeedsadeghi is on a distinguished road
I have done that. It's easy, you may have done a great mistake.
sakil2k3 likes this.
msaeedsadeghi is offline   Reply With Quote

Old   January 26, 2013, 03:52
Default
  #3
New Member
 
mahmood
Join Date: Sep 2012
Posts: 14
Rep Power: 4
sakil2k3 is on a distinguished road
Quote:
Originally Posted by msaeedsadeghi View Post
I have done that. It's easy, you may have done a great mistake.
Pls can u tell me where is the mistake. I m a M.Sc student and it is very urgent to complete my Thesis. Pls tell me how to do it. I will be ever grateful to you.
Thanx
sakil2k3 is offline   Reply With Quote

Old   January 26, 2013, 04:00
Default
  #4
Senior Member
 
SSL
Join Date: Oct 2012
Posts: 227
Rep Power: 5
msaeedsadeghi is on a distinguished road
It should be checked
msaeedsadeghi is offline   Reply With Quote

Old   January 26, 2013, 04:25
Talking
  #5
New Member
 
mahmood
Join Date: Sep 2012
Posts: 14
Rep Power: 4
sakil2k3 is on a distinguished road
Thank you for ur replay.
I have attached the full procedure and strategies.
I am eagerly waiting for ur replay.
Thanx
Attached Files
File Type: docx Simulation problem.docx (78.3 KB, 210 views)
sakil2k3 is offline   Reply With Quote

Old   January 28, 2013, 04:30
Default
  #6
New Member
 
mahmood
Join Date: Sep 2012
Posts: 14
Rep Power: 4
sakil2k3 is on a distinguished road
Please can anyone help me???
sakil2k3 is offline   Reply With Quote

Old   January 28, 2013, 10:23
Default
  #7
Member
 
Join Date: Dec 2012
Posts: 77
Rep Power: 4
beer is on a distinguished road
Hi
Of course the solid region doesn't move, the moment source holds it back. I would try to add a source term for gravity and use (U(n)-(n-1)) instead of U, where n is the timestep. I think with this modifications the source your solid reagion will sink to the bottom.

Greetings
beer is offline   Reply With Quote

Old   February 13, 2013, 16:56
Default solidification and vof
  #8
New Member
 
bouteldjamounir
Join Date: Feb 2011
Posts: 1
Rep Power: 0
coolmen_21 is on a distinguished road
I have a similar problem, can you send me your file . cas, to check on my side
coolmen_21 is offline   Reply With Quote

Old   January 14, 2014, 13:05
Default
  #9
New Member
 
Mahboobe Mahdavi
Join Date: Mar 2013
Posts: 22
Rep Power: 4
Mahboobe365 is on a distinguished road
Dear Saeed Sadeghi,

It seems that you are so happy knowing every thing and not helping others!
HyperNova and HamidIran like this.
Mahboobe365 is offline   Reply With Quote

Old   January 15, 2014, 04:17
Default
  #10
Senior Member
 
SSL
Join Date: Oct 2012
Posts: 227
Rep Power: 5
msaeedsadeghi is on a distinguished road
Quote:
Originally Posted by Mahboobe365 View Post
Dear Saeed Sadeghi,

It seems that you are so happy knowing every thing and not helping others!
Sorry, I am so busy.
msaeedsadeghi is offline   Reply With Quote

Old   March 14, 2014, 12:37
Default
  #11
Member
 
Join Date: Jan 2014
Location: Turkey
Posts: 36
Rep Power: 3
thermal energy is on a distinguished road
Hello dear friends,

I try to simulate melting process using VOF model. Pcm is solid phase initialy. Upper part of the solid there is air. They do not interpenetrate eachother. Natural convection is active. My primary phase air, secondary phase is pcm (phase change material). 85 % of volume is solid pcm initialy, 15% of total volume is air. I have some questions to be sure about my model.

First, I do not define interface between air and pcm. Mean, There are 3 cell zone: air, fin ,pcm but I do not create name selection for this interface in meshing. Is it ok?

Initial condition of Wholesystem temperature is 293K, but air temp is 300K. how to set air temp ?

For cell zone condition, all cell zone (air, pcm, fin ) show mixture. Also for example, if I change pcm to phase2, all of cell zone change to phase 2

Only for secondary phase (pcm), I patch within fluent by adapt-region-mark then patch. Ok? Value of volume fraction should be 1 or 0 while patching?? I am confused. Initialy pcm is solid.

For initial values phase-2 volume fraction is 0. Ok?

I did some iteration for a short time to check model. For results, in graphics-countors I can not see total temperature for phase 1 and 2 only entalphy values I see for phase1 and 2.

Also countors of solidification/melting it shows only mixture. Mixture show what?

Any interaction should I use like surface tension etc.

Thanks for your help. ı am in stuck.
zhangqinlian likes this.

Last edited by thermal energy; March 16, 2014 at 17:18.
thermal energy is offline   Reply With Quote

Old   October 9, 2014, 23:59
Default
  #12
New Member
 
Join Date: Oct 2014
Posts: 2
Rep Power: 0
HamidIran is on a distinguished road
A large mushy zone constant enforces the velocity in solid phase to go to zero (or to pull velocity). If you set it to zero the solid phase will move with a velocity obtained from the solution of the N-S equations. However, it will deform like a liquid.
HamidIran is offline   Reply With Quote

Old   December 13, 2014, 09:59
Default
  #13
Member
 
BijanKianoush
Join Date: May 2014
Location: Iran Tehran
Posts: 84
Rep Power: 3
HyperNova is on a distinguished road
Quote:
Originally Posted by Mahboobe365 View Post
Dear Saeed Sadeghi,

It seems that you are so happy knowing every thing and not helping others!
This who iranian are Mahboobe !!!
HyperNova is offline   Reply With Quote

Old   December 21, 2014, 10:25
Default
  #14
New Member
 
ZBC
Join Date: Dec 2014
Posts: 3
Rep Power: 2
zbc19910519 is on a distinguished road
Hi, would you please tell me more detail about how to acheive your idea to move the solid part, thx

Quote:
Originally Posted by beer View Post
Hi
Of course the solid region doesn't move, the moment source holds it back. I would try to add a source term for gravity and use (U(n)-(n-1)) instead of U, where n is the timestep. I think with this modifications the source your solid reagion will sink to the bottom.

Greetings
Hi can you tell some detail about how to accomplish your idea to move the solid part? thx

greetings

Quote:
Originally Posted by sakil2k3 View Post
Hi all, I have an enclosure spherical glass that is contain with phase change material (PCM) that is in its solid state initially. Heat will be applied above the melting temperature of the PCM around the spherical glass. During the melting process of the PCM there will be 2 phases present in the enclosure, the liquid and solid state. The solid PCM will drop to the bottom due to the gravitation force.

A 2d mesh was created and exported out from Gambit.

In Fluent, I used Solidification and Melting function in my simulation with VOF. All Material properties of the PCM and glass have set accordingly. The Gravity was also enable with Y=-9.81 in the Operating Condition and the Variable-Density Parameters was also set to an operating density of PCM, that is 777kg/m^3. The Body force weighted under Discretization in Solution Control was also selected.

After iterating for 200 sec, the values for the Liquid fraction was calculated but the simulation did not show the solid PCM sink to the bottom of the enclosure.

How can I make the solid PCM to sink to the bottom as it melts?
Thanks.
My friend , did you solve your problem?

Quote:
Originally Posted by coolmen_21 View Post
I have a similar problem, can you send me your file . cas, to check on my side
Hi, have you solved this problem?

Last edited by wyldckat; December 27, 2014 at 06:06. Reason: merge consecutive posts, posted a few minutes apart
zbc19910519 is offline   Reply With Quote

Old   January 2, 2015, 19:11
Default
  #15
Member
 
Join Date: Dec 2012
Posts: 77
Rep Power: 4
beer is on a distinguished road
Hi zbc19910519

Sure. Just use a very strong force, which is velocity dependent. For example 1e10*(vx-1) if you want your velocity component in x-direction to be 1. Just remember to tell the solver the derivative of this equation (I think it is called Source Coefficient in Fluent), otherwise the source will blow up your calculation. The solver uses the derivative to include the force in the implicit part of the solver (see Patankar).

Cheers
beer is offline   Reply With Quote

Old   January 3, 2015, 01:36
Smile
  #16
New Member
 
ZBC
Join Date: Dec 2014
Posts: 3
Rep Power: 2
zbc19910519 is on a distinguished road
Quote:
Originally Posted by beer View Post
Hi zbc19910519

Sure. Just use a very strong force, which is velocity dependent. For example 1e10*(vx-1) if you want your velocity component in x-direction to be 1. Just remember to tell the solver the derivative of this equation (I think it is called Source Coefficient in Fluent), otherwise the source will blow up your calculation. The solver uses the derivative to include the force in the implicit part of the solver (see Patankar).

Cheers
Thx, my friend, but if add a momentum source/sink in melting/solidification process, actually in which a huge momentum sink should be added to solid phase, thus, the solid phase is fixed, for the gravity is ignored compared to momentum sink, for example, if I initialize the case as the solid phase hangs in the air, and in next time step, it doesn't move any more and it is not mean to be, so how should I add this momentum sink properly to simulate the phase change process.

Thank you again !
zbc19910519 is offline   Reply With Quote

Old   January 3, 2015, 07:51
Default
  #17
Member
 
Join Date: Dec 2012
Posts: 77
Rep Power: 4
beer is on a distinguished road
I didn't really get what you're trying to say. Can you show a picture or something?
beer is offline   Reply With Quote

Old   March 17, 2015, 16:40
Default Question
  #18
Member
 
Mohamed
Join Date: Jan 2011
Location: Algeria
Posts: 43
Rep Power: 6
B.Hamada is on a distinguished road
how I can define the parameters of the liquid phase of PCM in Fluent initially solid.
Please help!!!!!!!!!!!!!!!!!!!!!!!!
B.Hamada is offline   Reply With Quote

Old   March 17, 2015, 17:24
Default
  #19
New Member
 
Join Date: Oct 2014
Posts: 2
Rep Power: 0
HamidIran is on a distinguished road
Hi Hamada,

Depending on the melting temperature that you define for your PCM in Materials tabs and the initial temperature you set in the "Solution Initialization" tab, the initial state of PCM is automatically determined. For example, if your PCM melting point is 309<Tm<311, if you set the initial temperature to 308, the PCM will start from solid phase. If your initial temperature is 312 the initial PCM phase will be liquid, and if the initial PCM temperature is 310, the initial volume fraction of liquid phase will be 0.5.

I hope this helps.

Hamid
HamidIran is offline   Reply With Quote

Old   March 17, 2015, 19:54
Default
  #20
Member
 
Mohamed
Join Date: Jan 2011
Location: Algeria
Posts: 43
Rep Power: 6
B.Hamada is on a distinguished road
Quote:
Originally Posted by HamidIran View Post
Hi Hamada,

Depending on the melting temperature that you define for your PCM in Materials tabs and the initial temperature you set in the "Solution Initialization" tab, the initial state of PCM is automatically determined. For example, if your PCM melting point is 309<Tm<311, if you set the initial temperature to 308, the PCM will start from solid phase. If your initial temperature is 312 the initial PCM phase will be liquid, and if the initial PCM temperature is 310, the initial volume fraction of liquid phase will be 0.5.

I hope this helps.

Hamid
Thank you
But i have a question, i have two data one for solid phase and the other for the liquid phase, for example there is a liquid density and a solid density, the same for the Specific heat capacity, and the Thermal conductivity.
but in fluent, I can enter the parameters of a single-phase.
please, what can i do?
B.Hamada is offline   Reply With Quote

Reply

Tags
vof phase change

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
VOF with melting of Phase change materials sakil2k3 FLUENT 3 March 17, 2015 16:26
How can I make the solid PCM to sink to the bottom as it melts? Lim FLUENT 6 December 22, 2014 03:29
PCM melting inside the Al cavity shoxan FLUENT 0 August 6, 2010 12:56
Laser vof melting solidification injection Michael FLUENT 0 February 20, 2007 22:28
vof + solidification melting James FLUENT 0 November 29, 2002 11:56


All times are GMT -4. The time now is 21:33.