
[Sponsors] 
April 28, 2009, 04:49 
Solidification in OpenFoam

#1 
New Member
Luke Christ
Join Date: Apr 2009
Posts: 7
Rep Power: 8 
Dear all,
I'm going to use OpenFoam for simulating a continuous casting problem. I'm completely new to OpenFoam, and it seems to me that there is not any solidification model in OpenFoam. Should I create solidification model myself? As my case is continuous casting, I should consider a pull velocity for solid phase, too. Is it possible without great difficulties in OpenFoam? I will be thankful for you reply. 

January 30, 2013, 02:38 
Updates

#2 
Member
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5 
Are you still interested in solidification in OpenFOAM?


May 3, 2013, 10:25 

#3 
New Member
Mac
Join Date: May 2013
Posts: 2
Rep Power: 0 
Hello
I'm newbie with OpenFOAM, but I'm interested to model solidification of steel. If i have well understood, there is no builtin solver? Do you know where I can find one? Regards 

May 3, 2013, 11:00 
solidification problem

#4 
Senior Member
Fabian Roesler
Join Date: Mar 2009
Location: Bad Friedrichshall, Germany
Posts: 167
Rep Power: 8 
Hi
In the thread http://www.cfdonline.com/Forums/ope...gproblem.html you can find several solvers for melting and solidification using the enthalpy method and the enthalpy porosity method. Regards Fabian 

May 6, 2013, 01:45 

#5 
New Member
Mac
Join Date: May 2013
Posts: 2
Rep Power: 0 
Thanks a lot


June 9, 2013, 16:02 
Explanation for terms in TEqn.H?

#6 
Member
einat
Join Date: Jul 2012
Posts: 31
Rep Power: 5 
Hello all!
I am trying to use meltFOAM for my project on solidification of lavas. I am trying to understand the terms appearing in TEqn.H. I identified the dT/dt, advection and diffusion terms, and then see terms that involve latent heat, and what I figure is some sort of a melt fraction term, but I don't understand the details  what is the exponent, where did all the "4" factors come from, etc.? Here's the code: fvScalarMatrix TEqn ( fvm::ddt(cp, T) + fvm::div(phi*fvc::interpolate(cp), T) + hs*4.0*exp(pow(4.0*(TTmelt)/(TlTs),2))/Foam::sqrt(pi)/(TlTs)*fvm::ddt(T) + hs*4.0*exp(pow(4.0*(TTmelt)/(TlTs),2))/Foam::sqrt(pi)/(TlTs)*(U & fvc::grad(T))  fvm::laplacian(lambda/rho, T) ); Thanks!!! 

June 11, 2013, 08:58 

#7 
Senior Member
Fabian Roesler
Join Date: Mar 2009
Location: Bad Friedrichshall, Germany
Posts: 167
Rep Power: 8 
Have a look into my paper on the solver:
F. Rösler, D. Brüggemann (2011): Shellandtube type latent heat thermal energy storage: numerical analysis and comparison with experiments. Heat and Mass Transfer, Vol. 47 Issue 8 , 10271033, DOI: 10.1007/s0023101108669 http://www.springerlink.com/content/b1tp01k2u7q8j432/ Keep in mind that the solver does not account for nonlinearity like the linear source method by Voller. The erfMeltSolver just reduces the nonlinearity effects in the energy conservation equation. For simple problems, the results converge with a very small error for the energy conservation. With increasing convective transport, the error increases and some iteration to account for the nonlinearity have to be performed. In my previous work, I use a linear liquid fraction function and the linear source method proposed by Voller. Regards Fabian 

June 11, 2013, 11:15 

#8 
Member
einat
Join Date: Jul 2012
Posts: 31
Rep Power: 5 
Excellent paper Fabian! Thanks for pointing it. Just what I needed.


July 6, 2013, 02:24 
Stefan problem

#9  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Sir
i am also trying to analysize the phase change problem.I want to know is it possible to do it on Fluent and if yes which solver /model i should go for. thanks Quote:


July 10, 2013, 04:38 

#10 
Senior Member
Fabian Roesler
Join Date: Mar 2009
Location: Bad Friedrichshall, Germany
Posts: 167
Rep Power: 8 
Hi dinesh
Yes, Fluent offers an EnthalpyPorosityMethod for simulation of solid/liquid phase change. Unfortunately I never used Fluent for such simulations so go on and find out yourself. Good luck. Regards Fabian 

July 20, 2013, 11:59 

#11  
New Member
Mojtaba
Join Date: Nov 2012
Location: Tehran
Posts: 20
Rep Power: 4 
Quote:
As Fabian said, FLUENT can solve Solidification/Melting problems. I myself tested and validated the solver for my thesis. There are excellent tutorial and theory materials from Ansys Fluent. Using Google, you may find them with no problem. Have fun! 

August 30, 2013, 05:18 

#12  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Quote:
Ref "H. Shmueli et al. / International Journal of Heat and Mass Transfer 53 (2010) 4082–4091" this paper says "page number 4086 says: As for the pressure discretization, only PRESTO!and BodyForceWeighted schemes are available for the VOF and mixture multiphase models." Does this means that multiphase has to be used for my case. Enabling this i find that courant number can be inserted which the paper says to be kept around 0.5. what about the phase 1 and phase 2. see this post by me which flotus replied on 11july Ref "http://www.cfdonline.com/Forums/systemanalysis/120322stefanproblem.html#post439934" eagerly waiting for the reply 

August 30, 2013, 21:59 

#13  
New Member
Mojtaba
Join Date: Nov 2012
Location: Tehran
Posts: 20
Rep Power: 4 
Quote:
I think this open access paper can help you: A Numerical Study on TimeDependent Melting and Deformation Processes of Phase Change Material (PCM) Induced by Localized Thermal Input 

September 1, 2013, 08:34 

#14  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Quote:


September 3, 2013, 01:30 

#15 
New Member
Mojtaba
Join Date: Nov 2012
Location: Tehran
Posts: 20
Rep Power: 4 
You didn't answer my question about whether the solidification zone is a fully filled container or not? If there is no gas in contact, using VOF is not necessary at all.
When you plot the contour of liquid fraction (F) the meltsolid (or meltgas) interface can obviously be seen and extract. Refer to Kim's paper, Fig. 9 to 11. 

September 3, 2013, 02:53 

#16  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Quote:
Plz see the plot i got https://www.dropbox.com/s/gg3mqyum1u3mx5z/frac4000s.png How can i extract the data regarding the anount or percent of solid melted/or mushy zone volume. thanks 

September 3, 2013, 11:15 

#17  
New Member
Mojtaba
Join Date: Nov 2012
Location: Tehran
Posts: 20
Rep Power: 4 
Quote:
From Menu Bar, go to Report and choose Volume Integrals... to open Volume Integral Window. In Volume Integrals Window: Choose Fluid (or whatever the solidification zone's name is) from right menu (Cell Zones), Choose Volume Integral from left menu (Report Type), Under Field Variable, select Solidification/Melting... and then Liquid Fraction, Click Compute button... and it's done! Good Luck 

September 4, 2013, 08:56 

#18  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Quote:
Can you plz tell me about the discretization method. I used SIMPLE method with PRESTO scheme (without gravity being aplied) i got some result. Now i am using PISO with Body Force weighted with gravity added slowly from 1m/s2 to 9.81 m/s2 (as recommended by some user). I run the simulation for some time and then i get divergence either in epsilon or( x y z componenet of velocity) can you suggest some way to overcome this. 

September 4, 2013, 10:01 

#19 
New Member
Mojtaba
Join Date: Nov 2012
Location: Tehran
Posts: 20
Rep Power: 4 
What are Ra and Pr numbers related to your simulation?
I've done an unsteady simulation of melting Gallium in a 2D cavity with Ra=2.2x10^5. I used SIMPLE for pressurevelocity coupling, PRESTO! for pressure discretization and SST kw for turbulence modeling. You can obtain a good convergence with a small enough time step and fine enough grid. If the convergence problem is still bothering, bring up more info about your simulation. 

September 4, 2013, 14:03 

#20  
Member
Vjoess
Join Date: Oct 2012
Posts: 52
Rep Power: 4 
Quote:
Regarding Pr and Ra i have not calculated. My hot water flow velocity is 0.1m/s at 350 K which is used to melt parafin wax which initially is at 293K and melting temperature is 313316K. Hot fluid flows inside while the outer cylinder(0.6m dia 1m length) is enclosing the wax. 

Tags 
continuous casting 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Superlinear speedup in OpenFOAM 13  msrinath80  OpenFOAM Running, Solving & CFD  18  March 3, 2015 06:36 
Solidification simulation with OpenFOAM  James  OpenFOAM  5  June 26, 2012 02:33 
64bitrhel5 OF installation instructions  mirko  OpenFOAM Installation  2  August 12, 2008 18:07 
Adventure of fisrst openfoam installation on Ubuntu 710  jussi  OpenFOAM Installation  0  April 24, 2008 14:25 
OpenFOAM Debian packaging current status problems and TODOs  oseen  OpenFOAM Installation  9  August 26, 2007 13:50 