interTrackFoam for rising bubble
I am working on adapting the hydrofoil tutorial case that comes with the interTrackFoam solver to a gas bubble rising in a stagnant liquid. What I've done until this point is the following:
- created a (2D) mesh with a 'hole' (the bubble) in the middle
- adjusted constant/polyMesh/boundary file to have a left, right, bottom and top and a freeSurface (the interface of the bubble) and of course the frontAndBackPlanes
- adjusted constant/faMesh/faMeshDefinition to have only frontAndBackPlanes attached to the freeSurface (unlike the hydrofoil case my freeSurface is not attached to other boundaries)
- adjusted 0/motionU to have slip on all sides and a fixedValue (000) on freeSurface
- adjusted 0/U to have zeroGradient on all sides and a fixedGradient (000) on freeSurface
- left 0/p as in the hydrofoil case
I can create the faMesh and interTrackFoam starts to run but in the first timestep the curvature gets out of hand and the solver crashes with 'floating point exception'.
Does anyone know what could be the reason for this?
And another question: does interTrackFoam remesh the domain while solving? (which I would say is needed for a rising bubble)
I've seen the new solver in OF 1.6-ext, but thanks for pointing it out anyway.
The rising bubble was actually just a testcase to try and get to grips with interTrackFoam. Eventually I would be interested in using it for a droplet or bubble in a partially wetting channel, but I think the boundary conditions for that are not yet implemented.
Ok, well, in that case you might be interested in the work of Sandeep Menon. He used interTrackFoam to model droplet breakup, if I remember right. You can find some info and even some downloadable code on Sandeep's website.
|All times are GMT -4. The time now is 15:55.|