[sega]

Hello.

Now I'm getting to my next problem:
The rise of sperical bubble to it's terminal velocity.

I used the test-setup described in http://test.interface.free.fr/Case01.pdf

Unfortunately there is nothing about the boundary conditions in this case.
As the boundary should not have an influence on the shape and even on the terminal veloctiy (if the computational domain is large enought), i chosed the following:

left: symmetryPlane for all values (U,pd,gamma) (I use an axisymmetric description
top: totalPressure (0) for pd, uniform (0 0 0) for U, zeroGradient for gamma
left & bottom: zeroGradient for pd and gamma, uniform (0 0 0) for U

To initialize the bubble I used funkySetFields to create a gamma-field at the center of the domain. Here is the image:



I used the same function to initialize the pressure inside the bubble to the value defined by sigma/radius:



Transport Properties were taken from the test-case description mentioned above.

Sounds bad, but it's real: the solver is doing nothing:


Exec : interFoam . .
Date : May 20 2008
Time : 15:16:11
Host : M1530
PID : 6392
Root : /home/sega/OpenFOAM/sega-1.4.1/run/rise
Case : .
Nprocs : 1
Create time

Create mesh for time = 0


Reading environmentalProperties
Reading field pd

Reading field gamma

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Calculating field g.h

time step continuity errors : sum local = 0, global = 0, cumulative = 0
#0 Foam::error::printStack(Foam:stream&)

So, there must be something wrong.
Any ideas?

My OpenFOAM-case is located here:
http://therealsega.th.funpic.de/openfoam/rise.tar.gz