CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

melting problem: looking for appropriate solvers

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

Like Tree146Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 8, 2015, 05:12
Default casting processes
  #201
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
Good morning SSSS,

I'll try to realize all your instructions over the day. The greates challenge seems to be the new convMeltFoam approach.

Attached you'll find the hole code. I also extended the immiscibleIncompressibleTwoPhaseMixture and incompressibleTwoPhaseMixture in order to store the thermophysical properties. You'll find the new transportModels attached as well.

Thanks,

ole
Attached Files
File Type: zip interTempFoam.zip (10.8 KB, 220 views)
File Type: zip transportModels.zip (54.3 KB, 135 views)
chaitanyaarige likes this.
kuechenrole is offline   Reply With Quote

Old   May 8, 2015, 08:27
Default simple test case for casting processes
  #202
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
Here is a simple test case for my solver attempt. It simulates a simple cube which is open to the atmosphere at the top and has an inlet at the bottom.
Attached Files
File Type: zip my_interFoam_fillingBox.zip (42.4 KB, 122 views)
chaitanyaarige likes this.
kuechenrole is offline   Reply With Quote

Old   May 8, 2015, 08:44
Default
  #203
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
A good introduction to the ideas, discretisation and solving strategies of interFoam seems to be the following pdf:
http://infofich.unl.edu.ar/upload/3b...7523c8ea52.pdf
kuechenrole is offline   Reply With Quote

Old   May 9, 2015, 04:10
Default casting processes
  #204
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
Here is a revised version of my 3 phase attempt. Following the hints of SSSS the main changes are:
  • reduced size, 2D domain inkl. symmetryPlane
  • rho and cp is equal for liquid and solid
  • rho becomes a volScalarField featuring the Boussinesq approach
  • added TEqn.relax()

The TEqn still doesn't work. Just cooling the cube results in higher temperatures than initiated and filling the isothermal cube results in floating Point error after raising and dropping T and liquid fraction.

Two ideas are still open:
  • Do anybody know an alternative approach to convMeltFoam, which doesn't use the convective term fvc::div( *** , fl )?
  • How can terms in the TEqn, which depend on the melt fraction (alpha1 = liquid + solid) be limited, other than just multiplying it with alpha1?

Thanks,

ole
Attached Files
File Type: zip solidiCubeCase.zip (14.8 KB, 60 views)
File Type: zip transportModels.zip (54.5 KB, 99 views)
File Type: zip interTempFoam.zip (11.3 KB, 83 views)
kuechenrole is offline   Reply With Quote

Old   May 11, 2015, 11:16
Default
  #205
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
New try with the approach from file:///home/konsole/Downloads/PhD_Thesis_ZSSALDI.pdf
Attached Files
File Type: zip interTempFoam.zip (11.5 KB, 206 views)
File Type: zip solidiCubeCase.zip (14.7 KB, 131 views)
geralt1997 likes this.
kuechenrole is offline   Reply With Quote

Old   May 22, 2015, 03:41
Default
  #206
New Member
 
Ole Richter
Join Date: Jun 2013
Posts: 11
Rep Power: 9
kuechenrole is on a distinguished road
Hello there,

does anybody know a common test case, which consider a phase change due to melting or solidification and the propagation of the PCM (phaseChangeMaterial). I want to benchmark my solver.

The frame is a casting process.
kuechenrole is offline   Reply With Quote

Old   May 28, 2015, 11:41
Default
  #207
Senior Member
 
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AUS
Posts: 136
Rep Power: 11
ahmmedshakil is on a distinguished road
Quote:
Originally Posted by kuechenrole View Post
Hello there,

does anybody know a common test case, which consider a phase change due to melting or solidification and the propagation of the PCM (phaseChangeMaterial). I want to benchmark my solver.

The frame is a casting process.
Hi Ole,
I did some validations for some extreme conditions of melting/solidification here: http://people.eng.unimelb.edu.au/ima...ngs/19/142.pdf

Cheers
shakil
chaitanyaarige likes this.
ahmmedshakil is offline   Reply With Quote

Old   July 16, 2015, 20:53
Default Species conservation equation in PimpleFoam
  #208
New Member
 
QP Dong
Join Date: Feb 2015
Posts: 7
Rep Power: 8
high_er is on a distinguished road
Hi there, I am a new former, which is described in attachment is species conservation equation . Now here is my problem, I don't know how to implement the species conservation equation in PimpleFoam accurately, if there were someone give me a guidance or show me how to do that, I will be very appreciate!species equation.png
high_er is offline   Reply With Quote

Old   September 2, 2015, 11:39
Default Regarding different versions of "convMeltFoam" solver
  #209
New Member
 
...
Join Date: Jun 2013
Posts: 19
Rep Power: 9
manalis is on a distinguished road
Hello phase-change foamers,

I downloaded the Anja's modified version of Fabian's "convMeltFoam" solver for OF2.3 (post #119) and also the latest (if I am not wrong) parallel version from Fabian's solver posted by himself (post #139). I tried the first one (let's call it "Anja's modified version" for clarity) and it runs normally but when I tried the other in serial case (let's call it "Fabian's latest parallel code") I get an error in Transport properties. Specifically, it asks for "mu0" and "muk" and from what I've seen they are declared this way in the respective "ReadTransportProperties.H"; when I put some values in the dictionary I receive an error about dimensions. I suppose that "mu" stands for dynamic viscosity with respective units [Pa*s]-->[kg/(m*s)]. On the contrary, "Anja's modified version" just asks for "nu" (kinematic viscosity) and the case runs without any problem (in OF 2.3.0). What could be the problem?

Thank you in advance!
manalis is offline   Reply With Quote

Old   September 3, 2015, 03:18
Default
  #210
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 15
fabian_roesler is on a distinguished road
Hi manalis

It's a long time ago since I posted this solver here. However, I remember that the viscosity was implemented temperature dependant.

Code:
mu = mu0*(scalar(1)-max(muk*(T-Tl),scalar(0)));
So mu0 has the same units as dynamic viscosity [kg/(m*s)] and muk has reciprocal temperature units [1/K].

Cheers

Fabian
chaitanyaarige likes this.
fabian_roesler is offline   Reply With Quote

Old   September 3, 2015, 06:26
Default
  #211
New Member
 
...
Join Date: Jun 2013
Posts: 19
Rep Power: 9
manalis is on a distinguished road
Hi Fabian and thanks for the immediate response!

After your suggestion and some minor changes in the "system" folder dictionaries I made it work! Up to now I have been using Fluent for my melting/solidification cases and very recently (actually last week) I decided to try the approach in OpenFOAM. So I am really new to this. Could you explain a little bit the expression for "mu" that you posted? What is the meaning of "scalar(0)", "scalar(1)" and "muk" in your code regarding the way that viscosity is changing with temperature? Have you tried to implement different properties for the 2 different phases with a high viscosity value for the solid phase?

Best regards
manalis is offline   Reply With Quote

Old   September 10, 2015, 09:34
Default convMeltFoam
  #212
New Member
 
...
Join Date: Jun 2013
Posts: 19
Rep Power: 9
manalis is on a distinguished road
Any thoughts/suggestions regarding viscosity equation? Apart from that, I would be very interested to know if someone has dealt with contact melting (e.g. melting inside a closed capsule) with this specific solver.

Best regards
manalis is offline   Reply With Quote

Old   October 7, 2015, 03:53
Exclamation solidificationMeltingSource
  #213
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 15
fabian_roesler is on a distinguished road
Hi folks. During a small solver research within OpenFOAM 2.4.x I found a new fvOption for solidification and melting:
Quote:
solidificationMeltingSource
Basically it uses the same approach as most solvers in this thread - the enthalpy porosity method by Voller and Prakash.
So just in case somebody is new to this thread or it's topic. This is one way to go besides the proposed solvers in this enlightening thread.

Cheers

Fabian
fabian_roesler is offline   Reply With Quote

Old   January 12, 2016, 16:13
Default
  #214
New Member
 
Paul
Join Date: May 2012
Posts: 23
Rep Power: 11
pmdelgado2 is on a distinguished road
Ole,
interTempFoam is created for what version of OF? The case files seem to imply version 2.3.0, but when I compile it, I get errors stating that the class 'immiscibleIncompressibleTwoPhaseMixture' (in alphaCourantNo.H) lacks a member named 'nearInterface', while UEqn.H and pEqn.H lack members named surfaceTensionForce. Do these members exist in some other version of OF's immiscibleIncompressibleTwoPhaseMixture?

Last edited by pmdelgado2; January 12, 2016 at 17:54. Reason: specified version number
pmdelgado2 is offline   Reply With Quote

Old   January 14, 2016, 11:00
Default Setting up melting problem with interTempFoam for 3 phases (ice, water, and air)
  #215
New Member
 
Paul
Join Date: May 2012
Posts: 23
Rep Power: 11
pmdelgado2 is on a distinguished road
Quote:
Originally Posted by kuechenrole View Post
Hi Ole,
I successfully compiled your solver based on Z. Saldi's thesis (see post #205). I'm trying to setup a problem with a block of ice immersed in heated air. I'm trying to use the 'setFieldsDict' to set the alpha1=0 in the air region and alpha1=1 in the water region.

Within the water region, I want to make part of it solid (ice) and the other part fluid (water). Should I set alphas=0 for the ice and alphas=1 for the liquid water region (only in the water region)? Or should I set alphas=0 for the ice region and alphas=1 everywhere else?

Please let me know how best to setup this problem with your solver.
Thanks,
Paul

Last edited by pmdelgado2; January 14, 2016 at 11:02. Reason: added reference post number
pmdelgado2 is offline   Reply With Quote

Old   January 14, 2016, 13:45
Default
  #216
New Member
 
Paul
Join Date: May 2012
Posts: 23
Rep Power: 11
pmdelgado2 is on a distinguished road
Hi Fabian,

Could you post your modified convMeltFoam solver that incorporates the third phase (air), as suggested by your thesis?

Quote:
Originally Posted by fabian_roesler View Post
Hi Folks

My PhD thesis on modeling and simulation of phase change processes in macro encapsulated latent heat thermal energy storage is now published. It is written in German and was published by the Logos Verlag Berlin.

Modellierung und Simulation der Phasenwechselvorgänge in makroverkapselten latenten thermischen Speichern
Thermodynamik: Energie - Umwelt - Technik, Bd. 24
Fabian Rösler
ISBN 978-3-8325-3787-6

http://www.logos-verlag.de/cgi-bin/e...87&lng=deu&id=

The convMeltFoam solver as well as the solver which takes unconstrained close contact melting and an additional gas phase into account are explained in detail. Experimental validation in an rectangular cavity is performed.
Hope you like it.

Cheers

Fabian
pmdelgado2 is offline   Reply With Quote

Old   February 19, 2016, 09:39
Default convMeltFoam OpenFoam-3.0.0
  #217
New Member
 
Alex Jarosch
Join Date: Dec 2015
Location: Austria
Posts: 24
Rep Power: 7
alexj is on a distinguished road
Dear phase-change Foamers,

did anyone convert the convMeltFoam solver from the post http://www.cfd-online.com/Forums/ope...tml#post484860 above to OpenFoam-3.0.0??

Cheers,
Alex
alexj is offline   Reply With Quote

Old   February 25, 2016, 08:59
Default
  #218
New Member
 
Alex Jarosch
Join Date: Dec 2015
Location: Austria
Posts: 24
Rep Power: 7
alexj is on a distinguished road
Quote:
Originally Posted by fabian_roesler View Post
Hi folks. During a small solver research within OpenFOAM 2.4.x I found a new fvOption for solidification and melting:

Basically it uses the same approach as most solvers in this thread - the enthalpy porosity method by Voller and Prakash.
So just in case somebody is new to this thread or it's topic. This is one way to go besides the proposed solvers in this enlightening thread.

Cheers

Fabian
I have uploaded a functional demo case of gallium melting here on the forum which might be helpful for people starting on phase change modeling using the fvOption "solidificationMeltingSource"

http://www.cfd-online.com/Forums/ope...tml#post586823

Best regards,
Alex
snak and Kummi like this.
alexj is offline   Reply With Quote

Old   March 16, 2016, 06:12
Default Melting in a body-fixed reference frame
  #219
New Member
 
Kai Schüller
Join Date: Jun 2015
Location: Germany
Posts: 8
Rep Power: 7
tetra-eder is on a distinguished road
Hi,
first I would like to thank for this very helpful thread.

I want to simulate a melting problem and I just thought that someone of you might help me.
A heat source is embedded into a phase change material and the heat source moves with a constant and given velocity. To solve for the phase change with convection, I use a solver based on Fabian Rösler's solver which was posted in this thread. I would like to use a fixed grid with a heat source fixed reference frame, so that the phase change material moves relative to the heat source. When using just an inlet and outlet condition for the velocity, the velocity changes within the solid phase, which is unphysical. I think the reason for this is mass conservation. Another possibility would be to add the velocity in the darcy term but this seems also not to work. Do you have any ideas how i can treat this relative motion phase change problem?

Thanks in advance!
Kummi likes this.
tetra-eder is offline   Reply With Quote

Old   March 16, 2016, 08:41
Default
  #220
New Member
 
Alex Jarosch
Join Date: Dec 2015
Location: Austria
Posts: 24
Rep Power: 7
alexj is on a distinguished road
@tetra-eder. Unfortunately I have not worked much with moving reference frame simulations yet in OpenFOAM so I don't have a answer to your question readily available for you.

Maybe if you would detail a bit more the simulation you want to achieve it will be easier for the forum to answer you.

If I understand you correctly, your heat source is moving in a static fluid? So you should be able to simulate that with an inlet, outlet boundary condition and let the fluid move instead of the heat source. However where is the solid phase in in your simulation? The term "heat source" implies that you want to melt something right? Or do you want to freeze the liquid onto the moving heat sink?

As you can see, a bit more detail might help.

One option which comes to mind is that you can also use a cyclic boundary condition for the inlet and outlet and use the fvOption "momentumSource" with a "meanVelocityForce" type to prescribe a fixed flow rate in your simulation for the fluid. However if you do that and keep the temperature field also cyclic, then you will effectively heat the fluid with your heat source as no thermal energy will escape the system.

Hope this helps somewhat.

Cheers,
Alex
Kummi likes this.
alexj is offline   Reply With Quote

Reply

Tags
melting openfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Melting and solidification with free surface problem? cqlwj123 CFX 6 July 25, 2013 02:46
Can I solve this problem by Fluent? Kai_kc FLUENT 1 October 27, 2010 05:29
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Melting Problem M FLUENT 0 April 29, 2007 16:07


All times are GMT -4. The time now is 15:22.