interFoam more Axi-symmetric two phase flow ?
 

Dear OpenFoam user's,

I am trying to setup a case in OpenFoam 1.6 for simulated a slug air - water two-phase flow in a micropipe but instead of a using all geometry; I trying to use an axi-symmetric approach. I build my mesh using a blockMesh in according to the UserGuidde.pdf from OpenFOAM using the wedge boundary for the front and back plane, still I am usage the setFields function to set a stratified flow to start the simulation.

Really, I don't know what's happend because my simulation stop suddenly after read p, U and alpha1 in the first step.

This is the message in the konsole:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-53b7f692aa41
Exec : interFoam
Date : Apr 30 2010
Time : 11:56:41
Host : javierbasurco-desktop
PID : 10345
Case : /home/javierbasurco/OpenFOAM/javierbasurco-1.6/run/tutorials/multiphase/interFoam/AxiSym2D_1
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading field p

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
time step continuity errors : sum local = 1.816956462e+279, global = 1.816956462e+279, cumulative = 1.816956462e+279
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.520219667e-11, No Iterations 333
time step continuity errors : sum local = 5.082309452e+270, global = 1.183822335e+268, cumulative = 1.816956462e+279
Courant Number mean: 9.800427968e-08 max: 1.018603267e-06

Starting time loop

Courant Number mean: 9.800427968e-08 max: 1.018603267e-06
Time = 1e-10

MULES: Solving for alpha1
Liquid phase volume fraction = -6.746125482e+276 Min(alpha1) = -5.096371243e+280 Max(alpha1) = 4.741194544e+271
#0 Foam::error::printStack(Foam::Ostream&) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::multiply<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/interFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator*<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/interFoam"
#6 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libfiniteVolume.so"
#7 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::fvc::flux<double>(Foam::GeometricField<doubl e, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/interFoam"
#8 main in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/interFoam"
#9 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#10 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/javierbasurco/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/interFoam"
Exception for flutuate point

Please, help me. I don't know where is the error or if is possible using interFoam with axi-symmetric mesh???:confused::confused::confused::( Also I post my case if anybody could be evaluate !!!


My deep regards

Respectfully
Ingenierias2003

I would check your boundary conditions and mesh, see this line:

time step continuity errors : sum local = 1.816956462e+279, global = 1.816956462e+279, cumulative = 1.816956462e+279

Sorry Javier, but you're going to have to look at the case by yourself for now. It is around that time where I am finishing up my thesis, so please check it out yourself and see if you can figure it out.

Hi, Javier,

I revised your blockMesh. I did not check your BCs. I did not get the initial error you encountered.

after untar the case, run ./runCase. It will do blockMesh, setFields, and then, run interFoam.

Pei