
[Sponsors] 
July 5, 2012, 12:07 
A difficult challenge...

#1 
New Member
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 8 
Hi,
I'm trying to implement a new solver for simulate the flow inside a shot sleeve used in high pressure die casting (hpdc) process. The main features of this type of flow are:  dynamic mesh to simulate the boundary displacement  multiphase flow (liquid metal and air)  temperature dependent physical properties (such as density, etc...)  heat exchange with the solid wall  phase change in metal due to possible solidification caused by loss of temperature I think interDyMFoam is a good starting point but, for example, it doesn't resolve the temperature field... Also chtMultiRegionFoam could be a useful solver for determine the heat transfer between the solid and fluid region and to account for temperature dependent physical properties... There is a way to couple this two solvers at runtime or it is necessary to compile a new solver? The final target is to estimate the quantity of trapped air in the metal phase due to the waves system. Any suggestion is welcome. By the way, any progress in the development of this project will be post... 

July 6, 2012, 02:47 

#2 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,215
Rep Power: 23 
If you're happy with using interDyMFoam, add the temperature equation to it. You can find all the code you need on the forum for interFoam, it should be directly usable in interDyMFoam too.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

May 21, 2013, 17:17 

#3  
Member
Jace
Join Date: Oct 2012
Posts: 77
Rep Power: 7 
Quote:
The approach I would like to take is to replace continuity, momentum and energy equations in the fluid region of chtMultiRegionFoam to account for contribution from different phases (liquid, gas), but I'm stuck on the actual implementation. I guess my problem is I'm not sure which terms within the conservation equations I should remove and/or add to get combine heat transfer with phase change. Any suggestions will be appreciated. Thanks, 

May 23, 2013, 10:13 
working progress

#4 
New Member
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 8 
Hi, my work proceed and now I've implemented the temperature equation into the interDyMFoam solver.
To do this i modified the twoPhaseMixture source code to implement the possibility for varying properties of two phases (rho, cp, lambda and nu) in function of temperature through reading files that i write.
I create a simple case to verify the solver:  Geometry: prism with square section, one single outlet patch;  Initial conditions: Aluminum (T=1000 K) on the bottom of the prism that fill 50% of entire height of it, the remainder is occupy by air (T=300K);  T boundary condition: all the boundary sets the temperature to 300 K. If i try to simulate a simple case where the wall is moving slowly the solution is optimal and nothing suggest some errors in fact i can view the results of diffusion witch change the temperature distribution near the boundaries and, so important, the convective motion of the air. But, if i try to simulate a case where the wall is moved more rapidly, the solution give problems. Up to 0.05 s the simulation run correctly (i suppose it viewing the alpha1 field, viewing the start of wavefront) but then alpha1 go crazy. I suppose this is connected to the mesh motion but i don't understand how it is. To clarify all, i post two images that first at t=0.05 s and second at t=0.06s; both visualize only aluminum (note alpha1 threshold). The patch named "movingWall" represent the piston. I also attach part of the log file of simulation. Can anyone help me? All the suggestions are appreciated. Thanks a lot! Best regards, Enrico 

June 6, 2013, 05:41 

#5 
New Member
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 8 
I would notice that i resolved my question...I made a careless mistake because I imposed a boundary condition on all walls for temperature equal to 300 K, but this value is smaller than solidification temperature of aluminum and so the viscosity values became very hight. This causes the solver to crash.


May 19, 2017, 16:18 
interDyMFoam with temprature equation

#6 
New Member
S P Chauhan
Join Date: Feb 2017
Posts: 1
Rep Power: 0 
hii, i am also start to work with interDyMFoam and try to add temperature equation and my solver and model compiled but my results are not coming properly. so i request to anyone plz send you case file...
thanks 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
A challenge for you in cfd(pipe problem)  mahesh9  CFX  9  March 24, 2012 11:53 
GAMBITFLUENTBLOOD CHALLENGE  christoforos  FLUENT  0  October 16, 2008 10:24 
why it is difficult to get converge in dissipation  frank  FLUENT  0  January 9, 2006 06:24 
cfx is Difficult to learn?  SAM  CFX  3  October 31, 2004 18:05 
Open Challenge for Concept of Turbulence  K.K.J.Ranga Dinesh  Main CFD Forum  2  June 3, 2003 13:41 