Solid/liquid phase change using large viscosity

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

 January 14, 2010, 11:01 Solid/liquid phase change using large viscosity #1 New Member   Gregor Join Date: Jan 2010 Location: Ljubljana, Slovenia Posts: 5 Rep Power: 8 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

 February 17, 2010, 10:41 #2 New Member   Gregor Join Date: Jan 2010 Location: Ljubljana, Slovenia Posts: 5 Rep Power: 8 I've found out that the PISO loop needs underrelaxation, now it is working fine . Regards, G

 February 18, 2010, 03:45 #3 Senior Member   matej forman Join Date: Mar 2009 Location: Brno, Czech Republic Posts: 104 Rep Power: 9 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

June 10, 2010, 06:31
#4
Member

Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 9
Quote:
 Originally Posted by gregor.vidmar153 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

Kind regards,
Jitao

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post fabian_roesler OpenFOAM 10 December 24, 2012 07:37 phdsantos FLUENT 0 March 20, 2009 11:19 qiulan CD-adapco 0 June 5, 2006 15:30 Sang-jin Lee CD-adapco 5 December 5, 2001 06:17 Yanhu Guo Main CFD Forum 4 January 24, 2001 00:16

All times are GMT -4. The time now is 14:46.