CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent Multiphase

Fluid Flow into Still Air from Faucet

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 3, 2018, 11:02
Default Fluid Flow into Still Air from Faucet
New Member
Jackie Yik
Join Date: Dec 2017
Posts: 2
Rep Power: 0
bladder1996 is on a distinguished road
I am having difficulties simulating water jetting from a capillary or faucet into still air and breaking up. To me, this should not be so difficult to simulate, but it has proven very time consuming.

I want to simulate the first two regimes of a jet breakup stability curve as shown in the first two boxes of the first attachment (Dumouchel 2008) in 3D. I have simplified the geometry for this with just a cylinder acting as the free space. I have shown the top of the cylinder and meshing I am currently using in the attachments. At the top of the cylinder I have the orifice of the faucet/capillary which is the finer meshed area. I used a parabaloid UDF for the inlet velocity profile profile and set the peak velocity as 2.38m/s. Water is the material that is coming from the faucet and air is what fills the domain. With my current properties I have a Reynolds number = 711 and Weber number of the gas = 0.0286. This falls into the first box of the attachment, Reyleigh regime, but the current simulation (also in attachments) looks more like the second box, first wind-induced regime. It took 10 days to run which is ridiculous and I have yet to see a drop break off, not to mention that I suspect the solution is not correct.

Here is a little more information and settings I used:

Total Nodes: 941200
Total Elements: Hexahedrals: 927276; and Quadrilaterals: 26992
**According to ICEM CFD Pre-meshing information before I converted it into mesh**
Minimum Orthogonal Quality = 0.738859
Maximum Ortho Skew = 0.261141
Maximum Aspect Ratio = 16.5582

Pressure-based, Transient, gravity pulling fluid down
VOF, 2 phases
Level Set selected
K-epsilon turbulence, realizable, standard wall functions
Surface tension on (obviously) but no wall adhesion
Boundary Conditions
top of cylinder where orifice inlet is and sides set as free slip wall
bottom of cylinder set as pressure outlet
**Tried setting the free slip walls as pressure outlets as well,
but not much difference to solution**

PISO scheme
Gradient: Least Squares Cell Based
Pressure: PRESTO!
Momentum: Second Order Upwind
Volume Fraction: Geo-Reconstruct
Turbulent Kinetic Energy, Turbulent Dissipation Rate, and Level-set Function: Second Order Upwind
kept default; did change them before but simulation is too slow to see much difference
set to default 1e-3
entire cylinder filled with air
Calculation Timestep
Variable method with Global Courant Number = 1 (was set to 0.5 before)
**Courant Number = 2 seemed to diverge sometime after 5000 or so time steps and I ran it over the holidays so I did not want it to diverge.**
Time Step Size = ranges from 2e-7 to 7e-8
max Iteration/Time Step = 50

I have tried looking around to speed up this simulation, but nothings has worked. I have looked at other papers on simulations similar to this (in 2D however) and they usually use similar settings. My computer's processing power is not the issue here I think because I am doing this on a 64GB RAM and 2 Intel Xeon processors workstation.

My main problem right now is that it takes too long to run. Each run of the simulation takes so long I do not have time to diagnose it. I hope I do not have to run a simulation for two weeks and then check to see if it is right, change some minor detail, and run it again for two more weeks.

There are a few things I have considered while looking around but I am not too sure:
  • My meshing could be too fine leading to errors, so possibly increasing mesh size will lead to a better and faster solution. The problem with this is that I need a fine enough mesh to be able to visualize the thin thread of fluid that attaches the droplet to the stream before breakup.
  • Maybe Discrete Phase Model might be better? I have read it can model droplet breakups. Might be slower though because it tracks the particles and is useful when there is not a lot of the discrete phase compared to continuous phase. In this case, the flow is continuous and eventually the water will take up a lot more of the domain which is why I do not think it should be used.
  • Maybe Pseudo-Transient Steady instead of Transient might make things faster? Once the jet and breakup stabilizes, it should have somewhat an oscillating solution since the droplets will periodically fall off and the jet length will increase and decrease because of this. I essentially do not care about the beginning of the flow because I know how it will behave for the most part (gravity will pull the stream down and it will be a growing jet). The problem with this is that I do not know the documentation for it and every time I try to implement it, the jet seems to literally explode upon leaving the orifice and I have no idea why.
  • I could use two planes of symmetry for the cylinder cutting down the mesh drastically. This would work for the first regime where everything is indeed symmetrical, but for the second regime where the flow gets quite turbulent and not symmetrical, I do not think the physics of the problem elicits the use of symmetry planes.

I am in quite a pickle and I do not know what to do. The only license I have is an academic one which does not allow customer support for this, I think, and I cannot run more than one simulation at a time so I cannot try different things out at the same time. Unfortunately 2D is not an option as I need the results for some other application. If there are any Fluent experts out there that might be able to help or some other post/documentation/resource that would be useful (I swear I tried using the search bar first!) or someone I could contact that would be able to help that would be great! Please let me know or private message me; I have been trying to do this for 3 months with little progress.
Attached Images
File Type: jpg Picture3.jpg (14.9 KB, 17 views)
File Type: jpg Picture5.jpg (89.3 KB, 12 views)
File Type: jpg Picture4.jpg (52.1 KB, 17 views)
bladder1996 is offline   Reply With Quote

Old   January 5, 2018, 13:26
New Member
R. Kim Eadies's Avatar
Join Date: Jan 2018
Posts: 11
Rep Power: 8
R. Kim Eadies is on a distinguished road
Hello Jackie,

Have you experimented with Implicit vs Explicit Volume Fraction Formulation (in Multiphase Model dialog box)? Explicit will be best suited for accuracy, but implicit may offer time savings by allowing a larger time step. May be worth playing with if you haven't.

Also, I didn't see anything about serial/parallel settings. I believe the academic license still provides for parallel processing, which will speed up your runs substantially if you aren't familiar.
"The more you do what you're already doing, the more you'll get what you've already got."
R. Kim Eadies is offline   Reply With Quote

Old   January 8, 2018, 10:50
New Member
Jackie Yik
Join Date: Dec 2017
Posts: 2
Rep Power: 0
bladder1996 is on a distinguished road
I haven't experimented with Explicit/Implicit yet, maybe I will look into that, thanks.

I know about parallel computing vs serial and I am using parallel. What I meant was I don't know if I can run two or more independent simulations at the same time.
bladder1996 is offline   Reply With Quote


breakup length, droplet formation, jetting, reyleigh regime, still air

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Difficulty in calculating angular velocity of Savonius turbine simulation alfaruk CFX 14 March 17, 2017 06:08
Average Rising Velocity as boundary condition Floing CFX 11 August 21, 2016 05:47
Want Impeller Driven Fluid Flow: What Inlet and Outlet BC to use for Centrifugal Pump Zev Xavier FLUENT 3 May 9, 2016 06:42
ansys cfx solver exit with return code 1!!!!! mhabibnia CFX 7 August 19, 2013 03:53
Simulation of air flow in a chamber Issa CFX 3 November 23, 2009 17:16

All times are GMT -4. The time now is 11:42.