[sErik]

Hi Foamers,

I'm currently working on a case with a turbulent flow in a tube. I want to compare the OpenFOAM results, with our previously ones, which we made with the Fluent "realizable k-epsilon" solver.

I thought about using the simpleFoam solver. Is the simpleFoam solver actually suitable for this case?
For initial condition, I have a velocity, peressure, temperature and the density of the fluid. But now I'm not sure about where to enter the temperature and the density for the calculation.

[sErik]

Now I'm using the buoyantSimpleFoam solver.

But when executing the solver, I recieve this error:

Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::surfaceInterpolation::makeWeights() const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::surfaceInterpolation::weights() const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::surfaceInterpolation::makeDeltaCoeffs() const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::surfaceInterpolation::deltaCoeffs() const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#7 Foam::fvPatch::deltaCoeffs() const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::fvPatchField<double>::snGrad() const in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/buoyantSimpleFoam"
#9 Foam::basicThermo::hBoundaryCorrection(Foam::Geome tricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#10 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#11 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > > >::New(Foam::fvMesh const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#12 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#13 main in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/buoyantSimpleFoam"
#14 __libc_start_main in "/lib/libc.so.6"
#15 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Gleitkomma-Ausnahme

I guess, something is wrong with my thermophysicalProperties file but I don't know what. This is the file:

thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;

mixture mixture 1 40.045 1000 0 1.8e-05 0.7;

pRef 100000;

My fluid is a mixture of water and another fluid. I think the the fault must be there, in that declaration. What are the valid entries for the mixture entry?

Any help is appreciated!