[pramodopen4foam]

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/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/baburao/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::fvMatrix<double> const&, char const*) in "/home/baburao/OpenFOAM/OpenFOAM-1.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/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/comsol"
#4 main in "/home/baburao/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/comsol"
#5 __libc_start_main in "/lib/libc.so.6"
#6 _start at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116


From function checkMethod(const fvMatrix<Type>&, const fvMatrix<Type>&)
in file /home/baburao/OpenFOAM/OpenFOAM-1.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

[sachin]

I guess there is problem with "rho" could u write the alphaEqn

Sachin

[pramodopen4foam]

alphaEqn has been attached,

Thanks for your response Sachin

Pramod

[l_r_mcglashan]

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.

[pramodopen4foam]

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,

[l_r_mcglashan]

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?

[alberto]

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?

[pramodopen4foam]

Originally Posted by l_r_mcglashan

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?

Firstly thank you very much for your patient response Laurence, its helping me

I am solving with incompressible fluids,
I am open for any suggesstion

Thank you

[pramodopen4foam]

Originally Posted by alberto

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?

Thanks alberto for your response, I am trying to combine the codes since my case has both setlingFoam n bubbleFoam involved,

Any suggestions are welcome,
Thank you

[Miguel António]

Hi Pessola, I need help to resolve this error.

--> FOAM FATAL ERROR:
incompatible dimensions for operation
[CH4[1 -3 -1 0 0 0 0] ] + [CH4[-1 0 -1 0 1 0 0] ]

From function void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::fvMatrix<Type>&, const char*) [with Type = double]
in file /home/miguel/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude/fvMatrix.C at line 1256.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::fvMatrix<double> const&, char const*) in "/home/miguel/OpenFOAM/miguel-8/platforms/linux64GccDPInt32Opt/bin/myporousGasificationFoam"
#3 Foam::tmp<Foam::fvMatrix<double> > Foam::operator+<double>(Foam::tmp<Foam::fvMatrix<d ouble> > const&, Foam::tmp<Foam::fvMatrix<double> > const&) in "/home/miguel/OpenFOAM/miguel-8/platforms/linux64GccDPInt32Opt/bin/myporousGasificationFoam"
#4 ? in "/home/miguel/OpenFOAM/miguel-8/platforms/linux64GccDPInt32Opt/bin/myporousGasificationFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/home/miguel/OpenFOAM/miguel-8/platforms/linux64GccDPInt32Opt/bin/myporousGasificationFoam"
Abortado (imagem do núcleo gravada)

[Kar93]

Hi Miguel,
I think the error might have crept in when you did not change the units of pressure while trying to create your own case.

Recently, I have fixed this error !!! Let me know if I am correct on this and hope this should soothe you.

Regards,
Karthik.