Dimension error after switching from simpleFoam to rhoSimpleFoam
 

Hi everybody,

I got an dimension error after switching from simpleFoam to rhoSimpleFoam.

I created a starting solution with simpleFoam and everything worked fine. Now I have reached steady state and I want to switch to rhoSimpleFoam to solve for the compressible case.

Unfortunateley I received an error...

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

Create mesh for time = 25

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
--> Upgrading k to employ run-time selectable wall functions
Backup original k to k.old
Writing updated k
--> Upgrading epsilon to employ run-time selectable wall functions
Backup original epsilon to epsilon.old
Writing updated epsilon
--> Creating mut to employ run-time selectable wall functions
Writing new mut
--> Creating alphat to employ run-time selectable wall functions
Writing new alphat
[1]
[1]
[1] Different dimensions for =
dimensions : [1 -1 -1 0 0 0 0] = [0 2 -1 0 0 0 0]
#0 Foam::error::printStack(Foam::Ostream&) in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::dimensionSet::operator=(Foam::dimensionSet const&) const in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator=(Foam::tmp<Foam::Geometri cField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/rhoSimpleFoam"
#4 Foam::compressible::RASModels::kEpsilon::kEpsilon( Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libcompressibleRASModels.so"
#5 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kEpsilo n>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libcompressibleRASModels.so"
#6 Foam::compressible::RASModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libcompressibleRASModels.so"
#7 main in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/rhoSimpleFoam"
#8 __libc_start_main in "/lib/tls/libc.so.6"
#9 _start at ../sysdeps/i386/elf/start.S:122
[1]
[1]
[1]
[1]
[1] From function dimensionSet::operator=(const dimensionSet& ds) const
[1] in file dimensionSet/dimensionSet.C at line 143.
[1]
FOAM parallel run aborting

Anybody an idea??


Best regards,
Sebastian

Hi Sebastian

In simpleFoam the pressure is normalized with density, i.e. you need not to specify it, however in rhoSimpleFoam the density varies, hence you should modify

1. Dimensions in your pressure file
2. The magnitude of the pressure, i.e. multiply by rho.

Best regards,

Niels

Hi Niels!

Thanks a lot! Now it works!


Sebastian