# Multiphase flow - simulation crashing

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

 November 14, 2010, 14:06 Multiphase flow - simulation crashing #1 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 Hello I am doing a multiphase flow simulation, with water and air (main phase) as the phases. I have a duct. Air enters the duct (I have specified total temperature and total pressure as the inlet condition, setting air volume fraction = 1 and water volume fraction = 0). Inside the duct, I inject water through small holes (typically 2 mm in diameter). At the holes I give an inlet boundary condition of a lower total temperature but higher total pressure (water volume fraction = 1, air volume fraction = 0). Basically I want the water to cool the air. At the outlet, I have specified bulk mass flow rate. At the remaining surfaces I have the adiabatic wall (no slip for air, free slip for water droplets) conditions and the rotating periodic conditions. I am using k-epsilon turbulence model, buoyant flow, etc. The problem is that the simulation is crashing. Either it gives a notice of very high Mach number or sometimes it just crashes without displaying any error. Microsoft gives an error that solver.hpmpi has crashed. I have tried changing the mesh, boundary conditions, turbulence models, but to no avail. Can someone please advise what to do and what can be the possible reasons of failure? Thanks.

 November 14, 2010, 17:57 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,638 Rep Power: 98 Can you give a drawing and flow rates to describe what you are doing?

November 14, 2010, 22:19
#3
New Member

Yash
Join Date: Nov 2010
Posts: 8
Rep Power: 7
Hello
I have attached a pic of the geometry and BC.
The boundary conditions are as follows:

Inlet:
Type: inlet
Total pressure: 237.38 kPa
Total temp.: 456.8 K
Air vol. frac. = 1
Water vol. frac. = 0

Holes:
Type: inlet
Total pressure: 262 kPa
Total temp.: 370 K
Air vol. frac. = 0
Water vol. frac. = 1

Outlet:
Type: outlet
Bulk mass flow rate = 0.338 kg/s

Other BC: rotating periodic on the sides (as seen), remaining surfaces: adiabatic wall (no slip for air, free slip for water droplets).

The timescale is set to auto.
I have done this simulation earlier when instead of water, the holes were injecting cold air. It converged then and gave results.
I hope this helps.

Thanks.
Attached Images
 CFX_BC.jpg (36.8 KB, 97 views)

 November 15, 2010, 03:31 #4 Member   Ali Torbaty Join Date: Jul 2009 Location: Sydney, Australia Posts: 72 Rep Power: 9 did you try it without water spray? if it works that way then start with low pressure and use an expresion to increase the water pressure by time. another point, you have a fixed mass flow at outlet which should be equal to the total mass of water and air, are you sure about that? maybe using fixed mass flow at inlet and an open boundary at outlet also helps to get better convergence.

 November 15, 2010, 04:08 #5 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 Hey It works without the water spray. The problem is that I cannot reduce the water pressure too low. It cannot go below the pressure of air surrounding the holes or else the water will not be injected (liquid can go from higher to lower pressure only). As regards the outlet, I have also tried putting static pressure conditions at the outlet but to no avail. That too crashes. I'll try this opening condition though. However isn't specifying a mass flow rate at inlet and opening at outlet a less robust condition? I thought a pressure condition makes the simulation robust. Thanks for the help and please keep helping and advising. It's kinda urgent.

 November 15, 2010, 04:16 #6 Member   Ali Torbaty Join Date: Jul 2009 Location: Sydney, Australia Posts: 72 Rep Power: 9 if it works without water spray, then it is likely that the initial condition is responsible for crashes. did you try to run it without water spray and then use the results to initilise the main model? if this doesn't work, try to use low pressure and increase it gradualy.

 November 15, 2010, 04:21 #7 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 I am using the same initializing conditions (on velocity) as I have used for the simulation without the water spray.

 November 15, 2010, 04:24 #8 Member   Ali Torbaty Join Date: Jul 2009 Location: Sydney, Australia Posts: 72 Rep Power: 9 then use the result of converged simple model (no spray) for intilising the complex model.

 November 15, 2010, 06:29 #9 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 Is the error coming because of only wrong initialization? Or can there be some other reason? Why I ask this is because shouldn't a good steady state result be independent of the initialization conditions? Also currently I'm giving only the air Cartesian velocity as the initialization condition. Rest all conditions (pressure, temperature, water velocities, volume fractions) are set to automatic. Please help. I'm stuck

 November 15, 2010, 17:37 #10 Member   h-h Join Date: Sep 2009 Posts: 38 Rep Power: 9 when i use smaller physical timestep in my pre the mach problem solved.(0.0001) when mach number became large(6.47) solver display notice.how much is it?

November 15, 2010, 18:26
#11
Member

Ali Torbaty
Join Date: Jul 2009
Location: Sydney, Australia
Posts: 72
Rep Power: 9
Quote:
 Originally Posted by yschati Is the error coming because of only wrong initialization?
In a multiphase model you need proper initial condition. If you don't have an existing model result to use as initial condition, still you need to avoid automatic setup and provide your best approximate of pressure and velocity.

 November 16, 2010, 06:38 #12 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 Hi I tried giving the results of the basic model as initial conditions. What I did was that, in the solver manager, I gave the results file of the converged simulation without the water spray and then ran the simulation. However, it is again crashing, due to very high Mach numbers. My flow is not at all supersonic.

 November 16, 2010, 09:49 #13 New Member   Yash Join Date: Nov 2010 Posts: 8 Rep Power: 7 Hey I tried the following things. All failed. a. Initialization using the results without water spray. b. Reducing the physical timestep by 2 orders of magnitude (1e-8). c. Specifying opening (and giving static temperature and pressure conditions) in place of inlet conditions. d. Reducing the total pressure at which water is being injected (262 to 250 kPa). What to do now?????

 November 16, 2010, 17:35 #14 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,638 Rep Power: 98 Work through the hints listed here: http://www.cfd-online.com/Wiki/Ansys...gence_criteria and let us know how you go.

 November 16, 2010, 18:29 #15 Member   h-h Join Date: Sep 2009 Posts: 38 Rep Power: 9 change your air if it is ideal gas. before solver crash stop it and see where mach is high.

 Tags crashing, multiphase, simulation

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Michiel CFX 17 April 21, 2010 10:14 rockewan FLUENT 0 April 6, 2010 12:34 cfd_multiphyiscs CFX 2 March 10, 2010 14:43 SimonH. OpenFOAM 0 October 27, 2009 05:39 Tom FLUENT 8 January 18, 2006 11:54

All times are GMT -4. The time now is 20:55.