|
[Sponsors] |
reading transportProperties of multiphaseInterFoam Error |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 3, 2016, 17:08 |
reading transportProperties of multiphaseInterFoam Error
|
#1 |
New Member
Abbas Rahimi
Join Date: Jan 2013
Posts: 28
Rep Power: 13 |
Hello everyone,
I'm trying to run multiphaseInterFoam for a case and I get this error. I wonder if anyone here can help me to overcome this error? Thank you. Create mesh for time = 0 Reading field p_rgh Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib64/libc.so.6" #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:? #5 Foam::multiphaseMixture::multiphaseMixture(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:? #6 at ??:? #7 __libc_start_main in "/lib64/libc.so.6" #8 at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116 Floating exception |
|
February 4, 2016, 05:27 |
|
#2 |
Senior Member
|
Hi,
Answering your question: use correct initial conditions (or properties, maybe you have made a mistake there). More detailed answer would be: If you look at the error, you can find key word "divide": Code:
#2 in "/lib64/libc.so.6" #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:? Code:
Foam::multiphaseMixture::multiphaseMixture ( const volVectorField& U, const surfaceScalarField& phi ) : ... nu_ ( IOobject ( "nu", mesh_.time().timeName(), mesh_ ), mu()/rho() ), ... Code:
Foam::tmp<Foam::volScalarField> Foam::multiphaseMixture::rho() const { PtrDictionary<phase>::const_iterator iter = phases_.begin(); tmp<volScalarField> trho = iter()*iter().rho(); for (++iter; iter != phases_.end(); ++iter) { trho() += iter()*iter().rho(); } return trho; } |
|
February 4, 2016, 05:41 |
|
#3 |
New Member
Abbas Rahimi
Join Date: Jan 2013
Posts: 28
Rep Power: 13 |
Alexey,
Thank you for response. Your answer totally make sense but I double checked my transportProperties and couldn't find any rho=0 anywhere. /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // phases ( air { transportModel Newtonian; nu nu [ 0 2 -1 0 0 0 0 ] 29.522; rho rho [ 1 -3 0 0 0 0 0 ] 1129; BirdCarreauCoeffs { nu0 nu0 [ 0 2 -1 0 0 0 0 ] 378.7336; nuInf nuInf [ 0 2 -1 0 0 0 0 ] 8.9686e-07; k k [ 0 0 1 0 0 0 0 ] 10; n n [ 0 0 0 0 0 0 0 ] 0.2407; } } base { transportModel Newtonian; nu nu [ 0 2 -1 0 0 0 0 ] 295.22; rho rho [ 1 -3 0 0 0 0 0 ] 1129; BirdCarreauCoeffs { nu0 nu0 [ 0 2 -1 0 0 0 0 ] 262.65; nuInf nuInf [ 0 2 -1 0 0 0 0 ] 0.000886; k k [ 0 0 1 0 0 0 0 ] 10; n n [ 0 0 0 0 0 0 0 ] 0.3005; } } tcap { transportModel Newtonian; nu nu [ 0 2 -1 0 0 0 0 ] 295.22; rho rho [ 1 -3 0 0 0 0 0 ] 1129; BirdCarreauCoeffs { nu0 nu0 [ 0 2 -1 0 0 0 0 ] 262.65; nuInf nuInf [ 0 2 -1 0 0 0 0 ] 0.000886; k k [ 0 0 1 0 0 0 0 ] 10; n n [ 0 0 0 0 0 0 0 ] 0.3005; } } bcap { transportModel Newtonian; nu nu [ 0 2 -1 0 0 0 0 ] 295.22; rho rho [ 1 -3 0 0 0 0 0 ] 1129; BirdCarreauCoeffs { nu0 nu0 [ 0 2 -1 0 0 0 0 ] 262.65; nuInf nuInf [ 0 2 -1 0 0 0 0 ] 0.000886; k k [ 0 0 1 0 0 0 0 ] 10; n n [ 0 0 0 0 0 0 0 ] 0.3005; } } skirt { transportModel Newtonian; nu nu [ 0 2 -1 0 0 0 0 ] 295.22; rho rho [ 1 -3 0 0 0 0 0 ] 1129; BirdCarreauCoeffs { nu0 nu0 [ 0 2 -1 0 0 0 0 ] 262.65; nuInf nuInf [ 0 2 -1 0 0 0 0 ] 0.000886; k k [ 0 0 1 0 0 0 0 ] 10; n n [ 0 0 0 0 0 0 0 ] 0.3005; } } ); refPhase tcap; sigmas ( (air tcap) 0.0 (air bcap) 0.0 (air skirt) 0.0 (air base) 0.0 (tcap bcap) 0.0 (tcap skirt) 0.0 (tcap base) 0.0 (bcap base) 0.0 (bcap skirt) 0.0 (base skirt) 0.0 ); // ************************************************** *********************** // |
|
February 4, 2016, 05:49 |
|
#4 |
Senior Member
|
OK. You did not find it in transportProperties, search in 0/alpha.* files. Since you did not post your case, other people can only GUESS what is wrong. And, in general, guessing game is rather annoying.
|
|
Tags |
multi phase flow, multiphaseinterfoam, transportproperties |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM] an error in Calculator's equation | immortality | ParaView | 12 | June 29, 2021 01:10 |
decomposePar problem: Cell 0contains face labels out of range | vaina74 | OpenFOAM Pre-Processing | 37 | July 20, 2020 06:38 |
How to install CGNS under windows xp? | lzgwhy | Main CFD Forum | 1 | January 11, 2011 19:44 |
Installation OF1.5-dev | ttdtud | OpenFOAM Installation | 46 | May 5, 2009 03:32 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |