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

interFoam more Axi-symmetric two phase flow ?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 30, 2010, 11:08
Default interFoam more Axi-symmetric two phase flow ?
  #1
Member
 
Javier Basurco
Join Date: Jan 2010
Location: Rio de Janeiro, Brazil
Posts: 33
Rep Power: 7
Ingenierias2003 is on a distinguished road
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:rintStack(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:perator*<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??? Also I post my case if anybody could be evaluate !!!


My deep regards

Respectfully
Ingenierias2003
Attached Images
File Type: jpg Captura_da_tela.jpg (44.4 KB, 45 views)
File Type: jpg Captura_da_tela-1.jpg (46.5 KB, 49 views)
Attached Files
File Type: gz AxiSym2D_Air-Water.tar.gz (3.4 KB, 15 views)
Ingenierias2003 is offline   Reply With Quote

Old   April 30, 2010, 15:34
Default
  #2
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
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
feijooos is offline   Reply With Quote

Old   May 10, 2010, 09:18
Default
  #3
Member
 
Javier Basurco
Join Date: Jan 2010
Location: Rio de Janeiro, Brazil
Posts: 33
Rep Power: 7
Ingenierias2003 is on a distinguished road
Dear Eelco Gehring,

Sorry but don't attend you, i retunr for a congress on multiphase flow. Do you revised my case and found any error in the implementation.

I look a forward to hearing from you!!
My deep regards

Respectfully
Ingenierias2003


Ingenierias2003 is offline   Reply With Quote

Old   May 10, 2010, 11:15
Default
  #4
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
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.
feijooos is offline   Reply With Quote

Old   May 10, 2010, 14:39
Default
  #5
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 272
Rep Power: 9
phsieh2005 is on a distinguished road
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
Attached Files
File Type: gz AxiSym2D_Air-WaterRev.tar.gz (4.0 KB, 35 views)
phsieh2005 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
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Two phase flow with condensation Sunnie FLUENT 0 May 31, 2009 10:02
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58
how 2 freeze 1 phase flow field & start lagrangian KK CFX 5 February 14, 2008 17:48
low speed compressible two phase flow?? cat CFX 0 November 15, 2005 08:59


All times are GMT -4. The time now is 18:24.