CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Dimesion error after switching from simpleFoam to rhoSimpleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 17, 2010, 07:06
Default Dimension error after switching from simpleFoam to rhoSimpleFoam
  #1
Member
 
Sebastian Saegeler
Join Date: Nov 2009
Location: Munich
Posts: 70
Rep Power: 8
sebastian is on a distinguished road
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:rintStack(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:perator=(Foam::dimensionSet const&) const in "/home/sebastian/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>:perator=(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:perator=(const dimensionSet& ds) const
[1] in file dimensionSet/dimensionSet.C at line 143.
[1]
FOAM parallel run aborting


Anybody an idea??


Best regards,
Sebastian

Last edited by sebastian; June 17, 2010 at 07:26.
sebastian is offline   Reply With Quote

Old   June 17, 2010, 07:24
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,693
Rep Power: 27
ngj will become famous soon enoughngj will become famous soon enough
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
ngj is online now   Reply With Quote

Old   June 17, 2010, 08:03
Default
  #3
Member
 
Sebastian Saegeler
Join Date: Nov 2009
Location: Munich
Posts: 70
Rep Power: 8
sebastian is on a distinguished road
Hi Niels!

Thanks a lot! Now it works!


Sebastian
sebastian is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Laminar simpleFoam and inviscid simpleFoam herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27


All times are GMT -4. The time now is 06:54.