CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Bubble in zero gravity conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 10, 2008, 11:52
Default Now that I did the "blocked" m
  #21
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Now that I did the "blocked" mesh I have a liquid phase fraction which is exactly as much as described in Case 10: http://test.interface.free.fr/Case10.pdf

The initial gamma-field is sharp if I take a look at the cell-gamma.

Unfortunately the pressure inside the bubble is still too low. It reaches a level of 0.943 N/m².
As described in the test-setup above it should be 1.04 N/m²

Now we should consider the next two options suggested:
2: No circular state. As I'm still not able to export a contour plot or similar to another post-processing tool, I can't get quantitative results of how the actual shape differs from a circle

3: Surface Tension and snGrad(). In Case10 there is a table listing the numerical error with different degrees of discretization.
In these terms (radius / gridspace) my grid has a value of 36.16 which is far more than the largest value (20) in Case10.
Concerning this the error of my model should not be that big!

So, there is still some thinking to do ...

What about these parasitic currents?
This whole setup is mainly about this.
Is there influenc that strong?

BTW: How can I compute the magnitude of these currents INSIDE the bubble?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 10, 2008, 16:21
Default Hi Sebastian, Could you pl
  #22
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sebastian,

Could you please upload your test case, then I will try to have a look at it in the beginning of next week.
If it is to big, do not hesitate sending it to my email.
Until then, enjoy yourself

- Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 11, 2008, 08:31
Default Hey. I send the files to yo
  #23
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hey.

I send the files to your account.
For all the other people who want to toy around with this problem: Here are the files attached:


__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 11, 2008, 08:34
Default Hey. I send the files to yo
  #24
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hey.

I send the files to your account.
For all the other people who want to toy around with this problem: Here are the files attached:

http://therealsega.th.funpic.de/open...ionSend.tar.gz
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 13, 2008, 12:16
Default Hi all Previously, I upload
  #25
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi all

Previously, I uploaded a small movie with a stick which turns into 3 spheres. Now I have added gravity just for the fun of it, and it looks so amazing that I wanted to share it with you.

Especially look at the symmetry which remains even though the flow is very complex. The movie shows the contour surface of gamma = 0.5 and all sides are walls.

www.student.dtu.dk/~s001581/OF_upload/cubeStickGravity.mpg

Have fun

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 13, 2008, 15:22
Default Hi Sebastian I have now loo
  #26
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sebastian

I have now looked at your setup, and I cannot say anything but: "I have run out of ideas.".

I have run your model and tried different things, and finally I tried to run it with double the discretization and the result is given below as "cell pd".



Suddenly large pressure gradients at the interface occur. Could it be that there is missing a term in the pressure or momentum equations? Otherwise these very large oscillations in the 'convergence' of the result is quite strange.

Please anyone who might have an idea, come forth

Best regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 14, 2008, 02:04
Default Hi The present scheme used
  #27
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi

The present scheme used for div(phirb,gamma) is interfaceCompression. I have made a small and fast test, and it looks like the pressure increases if another scheme is chosen.

I used vanLeer, but wouldn't know if that is a clever choice.

Best regads,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 14, 2008, 03:10
Default Well, ok. I think I will have
  #28
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Well, ok. I think I will have a talk about this divScheme with the guys at the university.

But first I have to know what div(phirb,gamma) stands for ...
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 14, 2008, 11:03
Default Is it possible that we are com
  #29
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Is it possible that we are comparing dynamic with static pressure?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 14, 2008, 11:43
Default Two things: 1. There seems
  #30
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Two things:

1. There seems to be a tool called interFoamPressure to calculate the static pressure instead of the dynamic.
Maybe this can settle the error?

2. How can I get a value for the magnitude of these currents inside the bubble?

Thanks so far. S.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 15, 2008, 14:47
Default Somehow I could make interFoam
  #31
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Somehow I could make interFoamPressure work and ran it over my results.

So now the pressure (p) is 0.85 N/m² ... 0.1 N/m² less than pd
Wrong direction as it should be 1.04 N/m² ...

Any ideas?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 16, 2008, 05:37
Default Hi Sebastian Sorry for not
  #32
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sebastian

Sorry for not replying, but I have been in the lab all week ... quite fun some times to see what reality actually looks like

div(phirb,gamma) is the convective term of the transport equation for gamma, where phirb is a surface normal flux. Changing from interface conpression to vanLeer gives a changes in pd from 0.944 to 0.957, though still not enough to reach 1.04.

I have not been able to find interFoamPressure, so could you please direct me to the source code?

There is a tool in paraFoam where you can see parts of the domain based on a threshold on one of the parameters ... thus using gamma > 0.5 would give you the interior of the bubble. You can find it in the toolbar next to the isosurface (the the left hand side).

Enjoy,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 16, 2008, 07:32
Default Hey Niels. I stumbled over
  #33
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hey Niels.

I stumbled over interFoamPressure in this post:
http://www.cfd-online.com/OpenFOAM_D...es/1/5619.html
The link (svn checkout) is in the third post ...

Greetings - S.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 16, 2008, 07:58
Default Hi Sebastian At first glanc
  #34
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sebastian

At first glance I see that the UEqn does not look like the one in interFoam, thus it is not surprising that the results deviate. E.g. in the one fvc:snGrad(gamma) is used and in the other fvc:grad(gamma) is used. The deviation seems odd and it is not the only one!

- Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 19, 2008, 09:45
Default Hey. I'm still playing with
  #35
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hey.

I'm still playing with this problem.

I did the following:
Defining a spherical gamma-field with this setup:



It's a quarter of a circel (due to symmetrical reasons).

I'm just starting the problem to have a look at the pressure inside the bubble.

Unfortunately the simulation is running at (adjustable) timesteps of 5e-6 by now.
This will take until eternity.

May it be that there is something wrong with the setup, because the solver is "stuck" like this?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 19, 2008, 09:53
Default Of course the file was too big
  #36
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Of course the file was too big.

Now it's here:
http://therealsega.th.funpic.de/open...10aSend.tar.gz
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 19, 2008, 10:44
Default Hi Sebastian Your velocity
  #37
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,595
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sebastian

Your velocity field looks something like this:



thus somethings definitely wrong at the symmetry condition.

1. What happens is you first try to make it a wall? Just to see if that removes the large velocities, it should

2. You have mapped gamma. Could you also map U and pd to get a initial estimate of the fields which are more or less in equilibrium in your symmetric case?

- Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   May 19, 2008, 14:44
Default Hey. Mapping (funkySetField
  #38
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hey.

Mapping (funkySetFields) the pd-field brought some even more disturbing results.

What about U? I think it should be zero all over the geometry, thus I leave it uniform (0 0 0).

I think I will have to do something about the right and top faces, concerning the pressure.
Setting it zero may not be the right approach...
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 20, 2008, 10:00
Default So, now I made some changes to
  #39
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
So, now I made some changes to the boundary conditions:

left & bottom: symmetryPlane for all variable (U,pd,gamma)

top: 0 for U, totalPressure for pd

right: 0 for U, zeroGradient for pd

I initialized the gamma and pd fields using funkySetFields so the values of a circle and the pressure inside the bubble due to surface tension.




The setting of pd with funkySetFields is destorying my boundary conditions totalPressure every time, making the solver demand a pdRefCell and Value. So I insered the totalPressure BC by hand to the nonuniform-pd-file.

The actual setup is here:
http://therealsega.th.funpic.de/open...ase10_a.tar.gz

The solver is stuck with timesteps of 10e-6 from the first 0.1 seconds. For this timestep there is "good" looking gamma-field, but a strange pressure-field:



So, I think there is some problem with the pressure boundary conditions.

Any more ideas?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   May 20, 2008, 10:10
Default Meanwhile I had a talk with my
  #40
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Meanwhile I had a talk with my academic advisor at the faculty I'm writing my student research project.

Obviously there is another student who is doing the same test-case like me.
He is investigating the square liquid volume under zero gravity conditions and looks after the final shape and pressure inside the bubble.
He is doing this with Fluent!

As a matter of fact he has the same problems!!!

On the same sized grid with equal amount of cells he is getting an error with the pressure jump of about 20%.

A guy from the faculty staff used the interFoam-solver with a (self-build) level-set-method and got even worse numbers.

So I think the problem I'm experiencing is not that trivial.
It even looks like my error of about 9% in the pressure jump is a pretty good value.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
GRAVITY Lorenzo FLUENT 0 September 19, 2007 10:10
gravity-table.scm scheme file for trans gravity Shankar FLUENT 0 May 31, 2006 04:30
Turbulent boundary conditions for bubble column anjai FLUENT 12 October 17, 2005 06:34
gravity driven flow- boundary conditions? Shankar FLUENT 4 November 26, 2003 16:45
Boundary Conditions for Bubble Column cfd-novice FLUENT 0 April 14, 2003 01:50


All times are GMT -4. The time now is 14:01.