Solidification in OpenFoam
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. |
Updates
Are you still interested in solidification in OpenFOAM?
|
Hello
I'm newbie with OpenFOAM, but I'm interested to model solidification of steel. If i have well understood, there is no built-in solver? Do you know where I can find one? Regards |
solidification problem
Hi
In the thread http://www.cfd-online.com/Forums/ope...g-problem.html you can find several solvers for melting and solidification using the enthalpy method and the enthalpy porosity method. Regards Fabian |
Thanks a lot
|
Explanation for terms in TEqn.H?
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*(T-Tmelt)/(Tl-Ts),2))/Foam::sqrt(pi)/(Tl-Ts)*fvm::ddt(T) + hs*4.0*exp(-pow(4.0*(T-Tmelt)/(Tl-Ts),2))/Foam::sqrt(pi)/(Tl-Ts)*(U & fvc::grad(T)) - fvm::laplacian(lambda/rho, T) ); Thanks!!! |
Have a look into my paper on the solver:
F. Rösler, D. Brüggemann (2011): Shell-and-tube type latent heat thermal energy storage: numerical analysis and comparison with experiments. Heat and Mass Transfer, Vol. 47 Issue 8 , 1027-1033, DOI: 10.1007/s00231-011-0866-9 http://www.springerlink.com/content/b1tp01k2u7q8j432/ Keep in mind that the solver does not account for non-linearity like the linear source method by Voller. The erfMeltSolver just reduces the non-linearity 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 non-linearity 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 |
Excellent paper Fabian! Thanks for pointing it. Just what I needed.
|
Stefan problem
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:
|
Hi dinesh
Yes, Fluent offers an Enthalpy-Porosity-Method 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 |
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! |
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 Body-Force-Weighted 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.cfd-online.com/Forums/system-analysis/120322-stefan-problem.html#post439934" eagerly waiting for the reply |
Quote:
I think this open access paper can help you: A Numerical Study on Time-Dependent Melting and Deformation Processes of Phase Change Material (PCM) Induced by Localized Thermal Input |
Quote:
|
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 melt-solid (or melt-gas) interface can obviously be seen and extract. Refer to Kim's paper, Fig. 9 to 11. |
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 |
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 |
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. |
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 pressure-velocity coupling, PRESTO! for pressure discretization and SST k-w 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. |
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 313-316K. Hot fluid flows inside while the outer cylinder(0.6m dia 1m length) is enclosing the wax. |
All times are GMT -4. The time now is 16:28. |