# Disturbed flow field at outlet boundary (Multiphase flow through pipe)

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 31, 2010, 05:01 Disturbed flow field at outlet boundary (Multiphase flow through pipe) #1 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 I'm want to analyse the transport of high density sand/water mixture throug a horizontal pipe with a diameter of 105 mm and a length of 3000 mm (3 m). The sand is defined as dispersed solid with a packing limit of 0.63 and a particle size of 0.2 mm. For the drag force i use the pre defined Gidaspow model. Bouyancy is activated. At the pipe inlet I have a uniform velocity which varies from 0.5 m/s up to 5 m/s. The volume fraction of sand is set to 0.4 and for water to 0.6. Before simulating the actual mixture I run a water only case which provides the initial flow field for the final analysis. The outlet is defined as average static pressure set to 1 bar. Reference pressure in the domain is also set to 1 bar. I inserted 2 pictures. On the picture of the outlet you can see the volume fraction of sand expressed in colours. Inconsistency at outlet boundary.jpg Because of the bouyancy the sand wants to concentrate on the bottom of the pipe. I would expact after 3 meters of pipe length, there should be a stable flow flied. But as you can see, in the region near the outlet there is an inconsistancy. Also the velocity vectors are showing a slightly donwards flow at the outlet boundary. How can this happen? The second picture is showing the inlet of the pipe. Flow field at inlet boundary.jpg It is clearly seen that the inlet has an uniform mixture. The sand tends to sink slowly as it gets further away from the inlet. After about 1 meter the flow flied reaches a stable situation, which remains there until it comes verry near the outlet. Why is the flow flied not stable at the outlet region? Did I define the outlet boundary (average static pressure, 1 bar) wrong? I tried several options, but i cannot get a proper result out of this...

 March 31, 2010, 07:05 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,721 Rep Power: 99 Have you tried using a velocity outlet and a pressure inlet? Just a guess.

 March 31, 2010, 07:15 #3 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Well, I thought I had. But I was just searching for the files, and i can't find them. It would be quite stupid to forget this option.... I will try this right now, and comment later.

 March 31, 2010, 09:28 #4 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Just finished the calculation... Results don't look good. Almost during the complete solution a wall was placed at the inlet, at the last iterations about 45% of the area was blocked. Convergence history looks like shit. Also the volume fraction of sand is I don't think this is the way to go... Most of the basic settings as I described in the first post came from the following model: http://works.bepress.com/sandip_lahiri/23/ The model I use has been successfully used by the author of this thesis. Do you have another suggestion?

 March 31, 2010, 10:30 #5 Senior Member     moein Join Date: Dec 2009 Posts: 132 Rep Power: 8 have you tried a Known pressure feild as outlet BC? i suggest at frist simulate a pipe without sand and investigate the velocity feild then try multiphase.

 March 31, 2010, 17:57 #6 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,721 Rep Power: 99 Also don't forget to try the pressure averaging options on the outlet BC. They might help.

 April 1, 2010, 02:52 #7 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Yesterday I couldn't fix it. Now i'm concentrating on the water only calculation. It seems the problem is not only present for the two phase flow. I have plotted the pressure contours near the outlet boundary for the water only flow: Pressure distribution near outlet.jpg It is clearly seen that the outlet boundary (right) has an influence on the pressure distribution. The velocity profile: Velocity distribution near outlet.jpg I would think the flow behaviour at the outlet should be the same as any cross section in the pipe where a fully developed flow field exists, or am I missing something?

 April 1, 2010, 06:30 #8 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,721 Rep Power: 99 As I said in my previous post, have you tried the pressure averaging options? I suspect this may be caused by the average static pressure option you are using.

 April 1, 2010, 10:44 #9 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Yes, all calculations performed with average static pressure outlet, the pressure averaging option was set to "average over whole outlet". I also tried the other option "circumferential" with the pipe center line difened as axis. This resulted in an error when the calculation reached the convergence target. Also tried a mass flow rate outlet, but still no improvement. I think I start all over again...

April 2, 2010, 06:09
#10
Super Moderator

Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,721
Rep Power: 99
Quote:
 all calculations performed with average static pressure outlet
That's my point. Have you tried turning pressure averaging off and specifying a pressure over the whole boundary?

 April 2, 2010, 08:04 #11 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Ow sorry, I misunderstood you. I also tried a static pressure outlet without averaging. The velocity profile looks better, but still there is a change near the boundary, see: Static pressure outlet velocity distribution.jpg

 April 3, 2010, 06:29 #12 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,721 Rep Power: 99 Are you sure your simulations are adequately converged and the mesh is of adequate quality? The bumpiness in your contour lines could be a sign things are not right there.

 April 7, 2010, 02:44 #13 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Solution is converged to target 1e-4. I also think there is something wrong with the mesh. Mesh at outlet: mesh at outlet.jpg That is what i meant by starting al over again; building a new mesh.

 April 7, 2010, 08:20 #14 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 I made a simple extrude mesh with hex elements. The velocity profile looks much better now: Static pressure outlet velocity distribution NEW.jpg There are no inflation layers, this is net yet the final mesh. I have to focus more on my meshing skills to get a proper mesh. So i will have some work for next week. Thank you for your comments!

 April 7, 2010, 18:38 #15 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,721 Rep Power: 99 Convergence to 1e-4 is generally loose convergence. I suspect your simulation is not fully converged yet. Do a sensitivity check on the convergence tolerance before changing anything else - there is not point working on a partly converged solution.

 April 8, 2010, 16:36 #16 Member   Freeman Adane Join Date: Apr 2010 Posts: 42 Rep Power: 8 hi Michiel, 1. I think using specifying 1bar pressure at the outlet means you're intrucing more errors into your solution filed. remember, is only pressure gradient is important in your governing equation. Therefore, setting reference pressure to 1bar, allows you to specify Zero at the outlet. In such a case, you reduce your error and therefore getting converged solutions. For a single phase flow (water) this shouldn't be a problem in CFX. If still, then I think you've more 'fundamental' problem than BC. 2. As for solid-liquid flow, is your flow laminar or Turbulent? My experience is that for turbulent you shouldn't be a problem if you run it in 'transient' mode. For laminar, I'm having a big challenge here. I had converged solution, however, my VF is wrong. CFX indicates that there's suspension/transportation of the sand, meanwhile, experiment indicate that there's no suspension of the sand and only a fraction of the sand is delivered. All the best

 April 21, 2010, 04:28 #17 Member   Michiel Join Date: Jul 2009 Location: The Netherlands Posts: 42 Rep Power: 9 Thank you both for your advice. I changed the unstructured mesh to a structured O-grid. The first calculations have no difficulty to converge to 1e-5. I will try to set the outlet pressure to 0 en see if this works out positive. Maybe I try a transient analysis. I still have a little difference in concentration profiles at the end of the pipe and just before the end of the pipe. Eventually i want to use the two phase model for turbulent pump calculations. The current analyses i use to get familiar with two phase cfd and to build a basic model which should be usable in slurry pump simulations. @freemankofi, I saw you postings about your problem. What you want to model is much more difficult I think. As I mentioned in a reaction on one of your postings, i think the settling of the solids is verry hard to model.

 April 21, 2010, 10:14 #18 Member   Freeman Adane Join Date: Apr 2010 Posts: 42 Rep Power: 8 Michiel, I wouldn't use the exact pipe length due to influence from outlet BC. If original pipe length is say, 10m, you want to add couple of meters to it. In such a case you don't extract your results from outlet but around 10m instead. This will avoid any outlet BC effect. If it is a fully-developed flow problem then, try to used normal single phase flow expression to estimate your entry length and here, add reasonable extras length to accommodate the outlet BC. If you have an idea about the outlet quantities such as volume fraction or concentration, then you might consider using "opening with p_spec=0" at the outlet. This BC is more robust and tend to give better convergence. Caution should also exercise when you're setting-up your final pump problems, especially the BCs. Here, you want to extend the domain such that any BCs will have no effect on the results.

 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 Saturn CFX 45 February 8, 2016 05:42 florian_krause OpenFOAM 22 June 13, 2013 21:25 saii CFX 2 September 18, 2009 08:07 venky FLUENT 1 April 5, 2006 04:09 John FLUENT 0 July 15, 2005 07:07

All times are GMT -4. The time now is 12:17.

 Contact Us - CFD Online - Top