dimensions error

Submitted by shokeyi on Tue, 03/04/2014 - 20:14

Helloeveryone,
There is an error for my test during coupling stage as follows,
[0] --> FOAM FATAL ERROR:
[0] [1]
[1]
[1] --> FOAM FATAL ERROR:
[1] incompatible dimensions for operation
[U[1 -2 -2 0 0 0 0] ] + [U[0 1 -2 0 0 0 0] ]
[1]
[1] From function checkMethod(const fvMatrix&, const fvMatrix&)
[1] in file /opt/openfoam211/src/finiteVolume/lnInclude/fvMatrix.C at line incompatible dimensions for operation
[U[1 -2 -2 0 0 0 0] ] + [U[0 1 -2 0 0 0 0] ]
[0]
[0] From function checkMethod(const fvMatrix&, const fvMatrix&)
[0] in file /opt/openfoam211/src/finiteVolume/lnInclude/fvMatrix.C at line 1316.
[1]
FOAM parallel run aborting
There is no dimension of [1 -2 -2 0 0 0 0] for my test, does anyone meet this error before or know how to modify it?

Thanks,
shokeyi

Forums:

alice's picture

alice | Wed, 03/05/2014 - 08:16

Good morning,
well I only came across this error when something with the dimensions was messed up. Is it possible that some wrong operation within your algorithm (Did you write a new solver?) or a boundary condition of the wrong dimensions causes this message? [1 -2 -2 ...] might as well result from an operation...
Cheers,
Alice

shokeyi | Wed, 03/05/2014 - 10:42

Hi Alice,
Yeah, I write a new solver for the test. Do you know where the error might occur for the algorithm?
Thanks.

alice's picture

alice | Wed, 03/05/2014 - 13:31

Hello,
well given that I have no idea of how your solver looks like that's a little hard to tell. Have a close look at all operations, operators, functions etc. you use. Seems like somewhere you tried to assign a velocity field with a quantity of wrong dimensions. Other sources are the setFields.H and the initial/boundary conditions files...
Cheers,
Alice

shokeyi | Thu, 03/06/2014 - 11:04

Hello Alice,
The dimensions error have been modified and could be run now. However, there is wrong calculation of alpha as follow
MULES: Solving for alpha1
Phase-1 volume fraction = 0.334474 Min(alpha1) = -0.00656116 Max(alpha1) = 1.20804
MULES: Solving for alpha1
Phase-1 volume fraction = 0.334323 Min(alpha1) = -2.11738e-25 Max(alpha1) = 1.45714
DICPCG: Solving for p_rgh, Initial residual = 2.28859e-08, Final residual = 2.28859e-08, No Iterations 0
time step continuity errors : sum local = 13176, global = -14.9368, cumulative = -14.9241
DICPCG: Solving for p_rgh, Initial residual = 2.30077e-08, Final residual = 2.30077e-08, No Iterations 0
time step continuity errors : sum local = 13322.1, global = -15.2523, cumulative = -30.1764
DICPCG: Solving for p_rgh, Initial residual = 2.3096e-08, Final residual = 2.3096e-08, No Iterations 0
time step continuity errors : sum local = 13421.4, global = 308.348, cumulative = 278.172
ExecutionTime = 2.16 s ClockTime = 2 s

Courant Number mean: 71615.9 max: 415336
Interface Courant Number mean: 26.7903 max: 13902.3
The alpha value should between 0 and 1. I gave these value in the source code but still not work. Could you let me know the setFields.H location? I could not find it in the source code.
Thanks.