[som]

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

[nimasam]

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

[som]

many thanks i'll try it asap

[fabian_roesler]

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

galliumErfMelting.zip
erfConvectiveMeltingPimpleFoam.zip

[som]

Originally Posted by fabian_roesler

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

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

[fabian_roesler]

Originally Posted by som

hello, Fabian

I read your artical. It's very interesting for me.
Now I try to solve the following case : moving spherical
heat source in a tube storrage. How do you think is it possble to apply your solver for this 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

[som]

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?

[aerospace1365]

Originally Posted by som

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 Som.
I want to cool a heat sink by butting it on the wax, just like your problem.
May I have your code?
would you please send it to my email?
aerospace1365@gmail.com

[zzj]

[QUOTE=fabian_roesler;329995][FONT=&quot]Hi

hello,fabian,your solver really works,but I have some questions about your solver.You define a variable named rhok,then add an annotation:Kinematic density for buoyancy force.
There is little information about the variable. I really can not figure out what kinematic density means,could you give me some suggestions?

[kera]

[QUOTE=zzj;637726]

Originally Posted by fabian_roesler

[FONT=&quot]Hi

hello,fabian,your solver really works,but I have some questions about your solver.You define a variable named rhok,then add an annotation:Kinematic density for buoyancy force.
There is little information about the variable. I really can not figure out what kinematic density means,could you give me some suggestions?

I hope this helps!

https://openfoamwiki.net/index.php/B...sinesqPisoFoam

Ricky

[zzj]

[QUOTE=kera;637754][QUOTE=zzj;637726]

hi, Ricky

thanks a lot,it really helps.

The essay gives a detailed explaination so that I can understand.

zzj

[som]

Originally Posted by nimasam

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

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

[nimasam]

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

be lucky

[AnjaMiehe]

Hello everybody,

Originally Posted by nimasam

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

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.

[som]

Originally Posted by AnjaMiehe

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?

[AnjaMiehe]

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

[som]

Originally Posted by AnjaMiehe

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?

[fabian_roesler]

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

[AnjaMiehe]

Originally Posted by som

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

Hi,

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

Regards, Anja

[salehda]

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