melting problem: looking for appropriate solvers
 

Did anybody solve melting problem with OpenFoam. Are there any appropriate solvers exist?

meltFoam
 

1) here is a solver i found some months ago in forum so if you search in forum you maybe find more

2) to use this solver, compile it with wmake command

3) it has a test case ;) too

many thanks i'll try it asap

New non-isothermal melting solver
 

Hi

I was writing a new melting solver for non-isothermal phase change. It can be used for isothermal melting as well. In the near future I want to combine meltFoam with my new solver which is significantly faster than the old solver. I although wrote an article about my work. Have a look here:

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/

What do you plan to do?

Regards

Fabian

Attachment 9805
Attachment 9806

It's been written for openFoam 1.3. I've managed to compile it after some minor changes for openfoam v2.01 .It works fine. Thanks again, Serge

yourwelcome, you can pay it forward :D
by putting your modified solver here for other users use ;)

be lucky

Hello. Your solver also works after some minore changes. Very fast for testcase geometry

Re: Melting problem
 

Hi som

I don't quite know what you mean with moving spherical heat source. But it sounds to me like a boundary condition, right? So it must be possible to attach my solver to this problem.
Glade you like my solver. I'm still working on it to allow moving rigid body (solid phase) and volume expansion of the melt into a gaseous phase with volume-of-fluid method.

Regards

Fabian

I have a ball with fixed T gradient on the boundary. this ball is placed into tube full of wax. Your solver acts very well with such problem.

Now i have new problem. I want make the ball moving in the tube.
I've tried to setup the boundary of the ball to type movingWallVelocity and run the case. Everything is good but a cant see any boundary displacement in ParaFoam.
Any ideas?

Hello everybody,

I quite agree and as I succeeded in compiling it in 2.1.0, see for the adapted solver below. The testcase given earlier by Fabian works and a Stefan-Problem testcase shows very good agreement between the numerical and the analytical solution.

What Stefan problem testcase do you mean? Where can i find it?

Hello Serge,

and welcome back to the thread. The Stefan Problem or Stefan Task is a classical testcase of solidification or melting, in fact where a moving boundary is met. The problem has an analytical solution depending on the various possible boundary conditions and as it is but 1 dimensional, it is very fast to validate parts of the code.

Try the key words "Stefan task solidification melting" in google or google scholar. I have just tried it and received lots of good hits.

Have a nice evening, Anja

yes, i know about it. I asked you about testcase for openfoam. Did you mean galliumErfMelting testcase or something else?

Stefan Testcase
 

Hi

I used the gallium melting test case by Gau and Viskanta [Gau C. and Viskanta R.: Melting and Solidification of a Pure Metal on a Vertical Wall, Journal of Heat Transfer, 1986, 108(1), pp. 174–181] for validation. However there are two Stefan problem test cases that are also often used when it comes to melting simulation. The one phase or two phase Stefan problem. Have a look into Cranks book [Crank J.: Free and Moving Boundary Problems, 1984, Oxford University Press, Oxford].

Regards

Fabian

Hi,

I used the testcase of Fabian, the gallium melting and the Stefan Problem.

Regards, Anja

Hi
I am trying to Compile the meltFoam.C , but I wasn't succesful because it was made for an older version - I had OpenFOAM-2.1.0 .
and while compiling it gave me an error message related to
readSIMPLEControls.H: No such file or directory
Any suggestions how to solve that please ?

Thanks
David

Hello David,

this I not trivial for me either as the solver was coded in 1.4.1, or at least that is the tag in the test case. Also, the test case proposes icoFoam, a PISO based solver, to solve the test case. I tried to rewrite meltFoam in 2.1.0 using the SIMPLE Algorithm but I found far more similarities in coding with the icoFoam solver in 2.1.0.

Therefore, and as my try for the meltSimpleFoam does compile but not solve the test case, here is my idea for the meltIcoFoam in 2.1.0 . For me, it compiles only with warnings and the rough test case given earlier in this thread works alright. I can't guarantee for anything.


I hope this is a help for you,
Regards, Anja

compilation meltFoam
 

Hi

The meltFoam solver is based on the old buoyantBoussinesqPisoFoam solver as buoyancy is taken into account. You should be able to translate the solver to OpenFOAM 2.1.x. I did that few month ago but without comments and bad syntax. Maybe I can post the solver for OF_2.1.x after a small revision in the next days if you still need it :).

Regards

Fabian

new MeltFoam solver for OF 2.1
 

Hi all

I posted that it might take some time to review the meltFoam solver and port it from OF_1.7 to OF_2.1 but I felt like doing it instantly. So here you go :). I hope you enjoy it.

Regards

Fabian

Hello,
thanks for the code Fabian, nice to see that I did a good job with the changings of erfConvectiveMeltingPimpleFoam posted here in #10 earlier ;-)
As far as my post on meltIcoFoam goes, I thought David was asking for an adapted coding of the solver in post #2 and that one is based on 1.4.1.
Thanks again for your code in general, I am having fun with it.

Anja