|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
Hi, everybody. I am new in Openfoam. Recently I am using interFoam to simulate a falling droplet under gravity. However, the droplet falls much more slower than expected.
I have been stuck here for one month; really hope some one can help me out. ![]() Attached is the p_rgh and U file. Thanks! Screenshot from 2014-08-28 10:45:39.jpg Screenshot from 2014-08-28 10:46:12.jpg |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Vignesh
Join Date: Oct 2012
Location: Darmstadt, Germany
Posts: 66
Rep Power: 13 ![]() |
Hi Cara !!
Did you switch on gravity ? There's a file named g in constant directory, Check inside the file whether g is set as value ( 0 -9.81 0 );
__________________
Thanks and Regards Vignesh |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
Quote:
Yeah, I believe I did edit the g file in constant. I am guessing that maybe my boundary conditions for the case is not right. I set all the surrounding boundaries to be patches, and set pressure and velocity boundary conditions accordingly. Thanks, Cara |
||
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
Hi Cara.
Can you specify more precisely what you mean by "slower than expected?". |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
I think further information about the case is necessary to help you. What does your mesh look like? How did you set the alpha field? Maybe you can upload your dictionaries? As is saw you based your simulation on the damBreak tutorial, is that right? If so, what did you change and what is standard?
Regards, Sebastian Last edited by Googlestuff; September 7, 2014 at 06:10. |
|
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
Yes, the case is based on the dambreak case. In the damBreak case, there is only one patch for the side to atmosphere, the others are walls; but in my case (3D), the droplet falls in the atmosphere, so all the boundaries are patches except for the groud it falls on. And the boundary conditions are the same with that in damBreak.
p_rgh: total pressure U: pressureInletOutletVelocty The mesh is not the directory (too large to upload successfully). But all the other stuff are in it. Thanks a lot for your help. By the way, I am also working on a case about channel flow with fixed heat flux into it, but I didn't find any similar case in the tutorial. I want to set the heat flux directly instead of set a temperature gradient. Do you have any suggestion?? Thanks, Best, Cara |
|
![]() |
![]() |
![]() |
![]() |
#7 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
Okay. I see you have five areas you assign boundary conditions to. Inlet, Outlet and Walls 1-3. But if you have a 3D Case there must be a sixth face. Maybe you forgot to assign a BC to that? Your dictionaries look just fine to me.
I am also simulating a falling drop under similar conditions, but i assign a start velocity to the drop so i don't have to wait for too long for the splash. Maybe that's an option for you too? Regards! |
|
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
Dear Sebastian,
Thanks for your quick reply. I group the walls at the same level, and name it wall 1,2,3. And then set boundary condition to these walls. I think it will be easier if the droplet is assigned a initial velocity, but I cannot see the deformation(might be very small) at the beginning. Thanks, Cara |
|
![]() |
![]() |
![]() |
![]() |
#9 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
Okay maybe i missunderstand. If you did group the walls, why did you assign the BC's to them individually? Maybe a picture of your mesh might help me to understand your geometry.
Code:
boundaryField { INLET { type pressureInletOutletVelocity; value uniform (0 0 0); } WALL_1 { type pressureInletOutletVelocity; value uniform (0 0 0); } WALL_2 { type pressureInletOutletVelocity; value uniform (0 0 0); } WALL_3 { type pressureInletOutletVelocity; value uniform (0 0 0); } OUTLET { type fixedValue; value uniform (0 0 0); } } Best regards, Sebastian |
|
![]() |
![]() |
![]() |
![]() |
#10 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
Hi, Sebastian,
I have attatched the pictures of the geometry, each picture for each boundary, except the inlet surface(because the maximam number of files is 5). I have ever checked the presssure and velocity fields in paraFoam, and then became more confused ![]() Best, Cara |
|
![]() |
![]() |
![]() |
![]() |
#11 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
Ah ok. Now I get it
![]() So with that configuration your BC's look reasonable. Why were you confused looking at p and U? I see your outlet is also a patch. I would recommend to define the type as "wall" in the blockMeshDict. I think in the tutorial they also use wall, because with this type the wall-distance is treated in a special way. |
|
![]() |
![]() |
![]() |
![]() |
#12 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
The outlet is wall; I think in paraView it set all boundaries to be patch.
Thanks, Best, Cara |
|
![]() |
![]() |
![]() |
![]() |
#13 |
New Member
Sebastian
Join Date: Nov 2013
Location: Germany
Posts: 14
Rep Power: 12 ![]() |
So what confuses you looking at the pressure and velocity fields?
|
|
![]() |
![]() |
![]() |
![]() |
#14 |
New Member
Fei Gao
Join Date: Aug 2014
Posts: 12
Rep Power: 12 ![]() |
It's like the droplet is in the areospace, the gravity is 10 times smaller. And the the accelation happens really slow.
Best, Cara |
|
![]() |
![]() |
![]() |
![]() |
#15 |
Senior Member
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 16 ![]() |
Hi,
I am also simulating a falling droplet in OpenFOAM. Although the droplet is just under the action of gravity and it is in the air, but it tends to right or left depending on the max Courant number when it is falling . Can anybody help me to solve it? |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Modelling falling solid sphere using interFoam VOF model | eelcovv | OpenFOAM Running, Solving & CFD | 6 | August 7, 2021 21:52 |
2D water droplet fall in InterFoam | muhanad | OpenFOAM Pre-Processing | 3 | May 3, 2019 03:48 |
Falling Droplet using InterFoam | cosbergel | OpenFOAM Running, Solving & CFD | 5 | May 16, 2018 05:21 |
droplet falling - VOF | bohis | FLUENT | 1 | July 10, 2013 04:28 |
sphere falling down droplet | naderafshar | FLUENT | 1 | December 24, 2011 10:30 |