June 20, 2012, 01:33
|
rhoSimpleFoam pipe flow
|
#1
|
New Member
Join Date: Mar 2012
Posts: 27
Rep Power: 14
|
Hi everyone,
I'm new to OpenFOAM and I am modeling a circular pipe using an axi-symmetric mesh for compressible turbulent flow. When I run the case using rhoSimpleFoam I get the following error message:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: no convergence criteria found. Calculations will run for 1000 steps.
Reading thermophysical properties
Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading field U
Reading/calculating face flux field phi
Creating turbulence model
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 1 average: 1
bounding epsilon, min: 0 max: 200 average: 200
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}
Starting time loop
Time = 1
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::tmp<Foam::Field<Foam::Vector<double> > > Foam:perator/<Foam::Vector<double> >(Foam::tmp<Foam::Field<Foam::Vector<double> > > const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#4 Foam::flowRateInletVelocityFvPatchVectorField::upd ateCoeffs() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#8 Foam::fv::laplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#9 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleTurbulenceModel.so"
#10 Foam::compressible::RASModels::kEpsilon::divDevRho Reff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libcompressibleRASModels.so"
#11
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoSimpleFoam"
#12 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#13
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/rhoSimpleFoam"
Floating point exception (core dumped)
I believe the cause of the error is in my boundary conditions, which I have attached, however I am not sure how to proceed from here. Any help that could be offered would be greatly appreciated.
|
|
|