CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Dimesion error after switching from simpleFoam to rhoSimpleFoam (http://www.cfd-online.com/Forums/openfoam-bugs/77232-dimesion-error-after-switching-simplefoam-rhosimplefoam.html)

sebastian June 17, 2010 07:06

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

ngj June 17, 2010 07:24

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

sebastian June 17, 2010 08:03

Hi Niels!

Thanks a lot! Now it works!


Sebastian


All times are GMT -4. The time now is 21:39.