
[Sponsors] 
October 14, 2009, 09:17 
Incompatible dimensions for operation

#1 
Member
Pramod
Join Date: Jul 2009
Posts: 30
Rep Power: 9 
Hi friends,
I tried to combine solvers settlingFoam and bubbleFoam, I could compile successfully, bu when I ran my solver, which gives error, which states, incompatible dimensions for operation [alpha[0 0 1 0 0 0 0] ] + [alpha[1 3 1 0 0 0 0] ]#0 Foam::error:rintStack(Foam::Ostream&) in "/home/baburao/OpenFOAM/OpenFOAM1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/home/baburao/OpenFOAM/OpenFOAM1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::fvMatrix<double> const&, char const*) in "/home/baburao/OpenFOAM/OpenFOAM1.5/applications/bin/linux64GccDPOpt/comsol" #3 Foam::tmp<Foam::fvMatrix<double> > Foam:perator+<double>(Foam::tmp<Foam::fvMatrix<d ouble> > const&, Foam::tmp<Foam::fvMatrix<double> > const&) in "/home/baburao/OpenFOAM/OpenFOAM1.5/applications/bin/linux64GccDPOpt/comsol" #4 main in "/home/baburao/OpenFOAM/OpenFOAM1.5/applications/bin/linux64GccDPOpt/comsol" #5 __libc_start_main in "/lib/libc.so.6" #6 _start at /build/buildd/glibc2.9/csu/../sysdeps/x86_64/elf/start.S:116 From function checkMethod(const fvMatrix<Type>&, const fvMatrix<Type>&) in file /home/baburao/OpenFOAM/OpenFOAM1.5/src/finiteVolume/lnInclude/fvMatrix.C at line 1184. FOAM aborting Can anyone help with this errorsince I am very new to OpenFoam and C++, Thanks in advance, Anticipating for your kind reply, Pramod 

October 14, 2009, 09:25 

#2 
Member
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 9 
I guess there is problem with "rho" could u write the alphaEqn
Sachin 

October 14, 2009, 10:06 

#3 
Member
Pramod
Join Date: Jul 2009
Posts: 30
Rep Power: 9 
alphaEqn has been attached,
Thanks for your response Sachin Pramod 

October 14, 2009, 11:43 

#4 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 15 
You're obviously using the compressible formulation of 'phi', which includes the density, and that's why you have a dimensions mismatch.
Before 'sticking' codes together, make sure you know exactly what equations you are solving for. If you look at the difference between the alphaEqn in settlingFoam and bubbleFoam the reason it isn't working should be obvious.
__________________
Laurence R. McGlashan :: Website 

October 14, 2009, 11:54 

#5 
Member
Pramod
Join Date: Jul 2009
Posts: 30
Rep Power: 9 
yes , as I mentioned since I dont know C++ upto the mark , I tried to stick with the solvers. Sorry to tell this, can anyone help me regarding this,
Thanks in advance, 

October 14, 2009, 12:02 

#6 
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 15 
This doesn't have anything to do with C++. I am assuming that you are solving:
fvm::ddt(alpha) + fvm::div(phi,alpha) + whatever The problem you have is that your dimensions are different. alpha is dimensionless so the first term will have a dimension of . The second term has phi, which has dimensions of , and alpha which is dimensionless, so that term should also have dimensions of . You are calculating phi using something akin to rho*U, which means your phi has different dimensions to the phi in the equation above. Is your case compressible/incompressible?
__________________
Laurence R. McGlashan :: Website 

October 14, 2009, 13:44 

#7 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 27 
In bubbleFoam the flux is volumetric:
phia = fvc::interpolate(Ua) & mesh.Sf(); in settlingFoam phi is created with compressibeCreatePhi.H in createFields.H What is the purpose of combining these two codes?
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

October 15, 2009, 04:24 

#8 
Member
Pramod
Join Date: Jul 2009
Posts: 30
Rep Power: 9 
Quote:
I am solving with incompressible fluids, I am open for any suggesstion Thank you 

October 15, 2009, 05:00 

#9  
Member
Pramod
Join Date: Jul 2009
Posts: 30
Rep Power: 9 
Quote:
Any suggestions are welcome, Thank you 

Tags 
incompatible dimensions 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Different dimensions for FATAL ERROR  retech  OpenFOAM Running, Solving & CFD  2  August 14, 2007 10:17 
Dimensions of laplacian in PISO loop  kumar2  OpenFOAM Running, Solving & CFD  2  July 3, 2006 14:34 
How to change the dimensions in STARCD?  Particle  CDadapco  5  April 23, 2006 10:16 
Fluent: changing dimensions of a geometry genetaed  Madhukar Rapaka  FLUENT  3  October 12, 2005 11:40 
SOURCE DIMENSIONS  Andreas Abdon  CFX  1  March 13, 2000 11:51 