CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

A difficult challenge...

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 5, 2012, 12:07
Post A difficult challenge...
  #1
New Member
 
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 6
e_boesso is on a distinguished road
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...
e_boesso is offline   Reply With Quote

Old   July 6, 2012, 02:47
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
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.
*Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology
akidess is offline   Reply With Quote

Old   May 21, 2013, 17:17
Default
  #3
Member
 
Jace
Join Date: Oct 2012
Posts: 77
Rep Power: 5
zhengzh5 is on a distinguished road
Quote:
Originally Posted by akidess View Post
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.
Hi, I'm looking to add phase-change functionality to chtMultiRegionFoam such that I can determine if there's any phase change happening within individual fluid region during the heat transfer with neighboring solid regions.

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,
zhengzh5 is offline   Reply With Quote

Old   May 23, 2013, 10:13
Default working progress
  #4
New Member
 
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 6
e_boesso is on a distinguished road
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
Attached Images
File Type: jpg 0_05.jpg (43.3 KB, 31 views)
File Type: jpg 0_06.jpg (44.3 KB, 33 views)
Attached Files
File Type: zip log.zip (39.7 KB, 6 views)
e_boesso is offline   Reply With Quote

Old   June 6, 2013, 05:41
Default
  #5
New Member
 
Enrico Boesso
Join Date: Mar 2011
Location: Padova, Italy
Posts: 16
Rep Power: 6
e_boesso is on a distinguished road
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.
e_boesso is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


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
GAMBIT-FLUENT-BLOOD 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


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