|
[Sponsors] |
Disturbed flow field at outlet boundary (Multiphase flow through pipe) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 31, 2010, 06: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: 17 |
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, 08:05 |
|
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
Have you tried using a velocity outlet and a pressure inlet? Just a guess.
|
|
March 31, 2010, 08:15 |
|
#3 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 10:28 |
|
#4 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 11:30 |
|
#5 |
Senior Member
|
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, 18:57 |
|
#6 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
Also don't forget to try the pressure averaging options on the outlet BC. They might help.
|
|
April 1, 2010, 03:52 |
|
#7 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 07:30 |
|
#8 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
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, 11:44 |
|
#9 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 07:09 |
|
#10 | |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
Quote:
|
||
April 2, 2010, 09:04 |
|
#11 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 07:29 |
|
#12 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
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, 03:44 |
|
#13 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 09:20 |
|
#14 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 19:38 |
|
#15 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
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, 17:36 |
|
#16 |
Member
Freeman Adane
Join Date: Apr 2010
Posts: 42
Rep Power: 16 |
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, 05:28 |
|
#17 |
Member
Michiel
Join Date: Jul 2009
Location: The Netherlands
Posts: 42
Rep Power: 17 |
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, 11:14 |
|
#18 |
Member
Freeman Adane
Join Date: Apr 2010
Posts: 42
Rep Power: 16 |
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 | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind turbine simulation | Saturn | CFX | 60 | July 17, 2024 06:45 |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 06:21 |
setup problems - LES pipe flow with cyclic BC (1) and direct mapped inlet (2) | florian_krause | OpenFOAM | 22 | June 13, 2013 22:25 |
HOW TO GIVE FLOW RATE IN OUTLET BOUNDARY CONDITION | venky | FLUENT | 1 | April 5, 2006 05:09 |
Boundary condition for not fully-developed flow | John | FLUENT | 0 | July 15, 2005 08:07 |