A difficult challenge...
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...
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.
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.
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,
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.
|All times are GMT -4. The time now is 11:26.|