CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solid/liquid phase change using large viscosity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2010, 10:01
Default Solid/liquid phase change using large viscosity
  #1
New Member
 
Gregor
Join Date: Jan 2010
Location: Ljubljana, Slovenia
Posts: 5
Rep Power: 6
gregor.vidmar153 is on a distinguished road
Hello!

Does anybody has some experiences on the modeling of the phase change from solid to liquid taking large viscosity for solid. I know that the problem is stiff but ...

I'm modeling it with a corrected version of the buoyantBoussinesqPisoFoam. I have a normal wax taking 100 000x bigger viscosity for solid than for liquid (smaller it does not bring the correct heat transfer in the solid). I model it for the 2D buoyancy driven cavity starting from the full solid and heating it from one side. The viscosity around the phase change is changing linearly or even quadraticaly in an interval of 1K, the difference in left and right side wall temperatures is 15K.

Running the model with pure "solid" or pure liquid works nice. The problem is, that when the phase change begins, the time step decreases for the same factor as is the ratio between solid and liquid phases, thus model becomes unusable. I'm sure that the problem is somewhere in change of the viscosity, while the time step decrease happens even if a leave all the constats except viscosity unchaged.

After trying everithing falling to my mind (shemes, mesh, better v-T copling, ...) I'm close to the conclusion that the problem is in the model itself. Am I right???

Gregor
gregor.vidmar153 is offline   Reply With Quote

Old   February 17, 2010, 09:41
Default
  #2
New Member
 
Gregor
Join Date: Jan 2010
Location: Ljubljana, Slovenia
Posts: 5
Rep Power: 6
gregor.vidmar153 is on a distinguished road
I've found out that the PISO loop needs underrelaxation, now it is working fine .

Regards, G
gregor.vidmar153 is offline   Reply With Quote

Old   February 18, 2010, 02:45
Default
  #3
Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 91
Rep Power: 7
matejfor is on a distinguished road
Hi Gregor.

I would suggest you think a bit what the under-relaxation is doing with PISO.
Pure PISO is getting the result in one step with several correction (you set the number of corrections). When you under-relax variables in PISO, you're relaxing the physics. Like for temperature it would mean adding heat capacity - thus slowing down the temperature propagation. If you have problems with convergence with large time steps, you should keep the time step small, or use PIMPLE, which is basicly combination of PISO with inner SIMPLE loops, which you can under-relax as you wish.
This book is good reference: J. H. Ferziger, M. Peric, Computational Methods for Fluid Dynamics, Springer, 3rd Ed., 2001.
I hope people will correct me, if I'm wrong.

Good luck
Matej
matejfor is offline   Reply With Quote

Old   June 10, 2010, 06:31
Default
  #4
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 7
awacs is on a distinguished road
Quote:
Originally Posted by gregor.vidmar153 View Post
Hello!

Does anybody has some experiences on the modeling of the phase change from solid to liquid taking large viscosity for solid. I know that the problem is stiff but ...

I'm modeling it with a corrected version of the buoyantBoussinesqPisoFoam. I have a normal wax taking 100 000x bigger viscosity for solid than for liquid (smaller it does not bring the correct heat transfer in the solid). I model it for the 2D buoyancy driven cavity starting from the full solid and heating it from one side. The viscosity around the phase change is changing linearly or even quadraticaly in an interval of 1K, the difference in left and right side wall temperatures is 15K.

Running the model with pure "solid" or pure liquid works nice. The problem is, that when the phase change begins, the time step decreases for the same factor as is the ratio between solid and liquid phases, thus model becomes unusable. I'm sure that the problem is somewhere in change of the viscosity, while the time step decrease happens even if a leave all the constats except viscosity unchaged.

After trying everithing falling to my mind (shemes, mesh, better v-T copling, ...) I'm close to the conclusion that the problem is in the model itself. Am I right???

Gregor
Dear Gregor,

I encountered a similar problem. I want to simulate a two phase (high viscous liquid and air) flow problem consiering soliderification. The phase change from liquid to solid can be achieved by taking large viscosity for solid.

I want to create a new piecewise viscosity moedl:

(1) when T<Tg, viscosity is equal to a constant:

nu = nu0;

(2) when T>=Tg, viscosity varies with temperature, pressure and strain rate:

nu= (B_*exp(B1_*pd)*exp(Tb_/T))/(RHO_*(scalar(1)+pow((B_*exp(B1_*pd)*exp(Tb_/T))*strainRate()/t1_, scalar(1)-n_)));

where B, B1, Tb, t1, n are material constants.

When compiling this model using wmake, I encountered the following problems: http://www.cfd-online.com/Forums/openfoam-solving/76995-piecewise-viscosity-model.html
Please give me some suggestions. Thanks in advance.

Kind regards,
Jitao
awacs is offline   Reply With Quote

Reply

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
Solid/liquid phase change fabian_roesler OpenFOAM 10 December 24, 2012 06:37
Phase change phdsantos FLUENT 0 March 20, 2009 11:19
Solid Phase change material qiulan CD-adapco 0 June 5, 2006 15:30
Phase change process Sang-jin Lee CD-adapco 5 December 5, 2001 05:17
help needed about phase change Yanhu Guo Main CFD Forum 4 January 23, 2001 23:16


All times are GMT -4. The time now is 18:12.