[6863523]

Dear All,
I am trying to use a set of dimensionless equations in OpenFOAM, as in the post :
How does one turn off the dimensional checking
and in
http://www.cfd-online.com/Forums/ope...quation-2.html
I have run a case with dimension forms without any problem.
Attached please find the new solver I want to use. Everytime I run the case with the new solver, there exists an error as below. I have copied the controlDict out to my working directory and changed the dimensionSet from 1 to 0, but the same error shows up.
\*--------------------------------------------------------------------------------------------------*/
Build : 3.0.1-d8a290b55d28
Exec : myIcoFoam
Date: Dec 10 2016
Time: 17:53:57
Host:
PID:
Case:
nPorcs: 1
sigFpe: Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChcking: Monitoring run-time modified files using timeStampMaster
allowSystemOperation : Allowing user-supplied system call operations

//************************************************** ************************************************//
Create Time

Create mesh for time=0

PISO: Operating solver in PISO mode

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Starting time loop

cellSource functionObjects_allvolAverage:All():
total cells =103462
total volume = 0.0400007

Time=0.0125

Courant Number mean:0 max:0

-->FOAM FATAL IO ERROR:
incompatible dimensions for operation
[U[0 0 -1 0 0 0 0]]+ [U[0 -1 0 0 0 0 0]]

From function checkMethod (const fcMatrix<Type>&, const fvMatrix<Type>&)
in file /usr/local/OpenFOAM/OpenFOAM-3.0.1/src/finiteVolume/lnInclude/fcMatrix.C at line 1295.


FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
# Foam::IOerror::abort() at ??:?
#2 void Foam::checkMethod<Foam::Vector<double>>(Foam::fvMa trix<Foam::Vector<double>> const&, Foam::fvMatrix<Foam::Vector<double>> const&, char const*) at ??:?
#3 ? at ??:?
#4 ? at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? at ??:?

Aborted

I wonder how to make a change to the dimensionedType.C.
I wonder if any one could give me some help.
Thank you in advance.
Best Regards,
Bill

[6863523]

Hi,
I copied the dimensionedType.C out, as what I have done for controlDict; Even, I deleted the dimension checking syntax, there is still the same error, which means the dimensionedType.C was not changed at all. Or the problem is caused by something else?
Please help, I have stuck for several days.
Thank you in advance.
Best Regards,
Bill

[6863523]

Any suggestion would appreciated.
Thank you in advance.

[6863523]

The files have been updated after some modifications.
Thanks.

[6863523]

Just updated.
If I deleted the 'nu' in my new solver and the 'nu' definition in the transportProperties.
It will show that the nu is undefined in the transportProperties.
The problem seems to be clear, the newly defined solver may not be used at all by the program.
I will define the new solver again to test if it can be used or not.

[6863523]

Newly updated.
I have tried the new solver for other cases, the new solver can be used, showing that the undefined Re, which I have used in the new solver.
The problem seems that the dimension checking was not turned off at all by copying out the controlDict and switching the dimensionSet from 1 to 0.
Since I do not have the privilege to modify this file, I wonder if there exists any other solution, except copying the controlDict out.
Still, any suggestion would be appreicated.

[6863523]

Any information missed, please let me know.
Could any one try to turn off the dimension checking and let me know if it is possible?
Thank you in advance.

[6863523]

The dimensionless Navier-Stokes equations that I want to solve in the new created solver is:
fvVectorMatrix UEqn
(
Re*(fvm::ddt(U))
+ Re*(fvm::div(phi, U))
- fvm::laplacian(nu,U)
+ F
);

solve(UEqn == -fvc::grad(P));

The '' F '' is a unit vector, standing the gradient of the nondimensionalized mean pressure, as can be seen in the attached equations. I have it like an input (I will study my flow in function of given values of Re and F).

I have decomposed the traditional p into F+ grad(p), as in the attachment.
Thank you. Any information needed, please let me know.

[M. Nasr]

Originally Posted by 6863523

The dimensionless Navier-Stokes equations that I want to solve in the new created solver is:
fvVectorMatrix UEqn
(
Re*(fvm::ddt(U))
+ Re*(fvm::div(phi, U))
- fvm::laplacian(nu,U)
+ F
);

solve(UEqn == -fvc::grad(P));

The '' F '' is a unit vector, standing the gradient of the nondimensionalized mean pressure, as can be seen in the attached equations. I have it like an input (I will study my flow in function of given values of Re and F).

I have decomposed the traditional p into F+ grad(p), as in the attachment.
Thank you. Any information needed, please let me know.

Hello, Have you managed to develop the code to solve the Non-dimensional Navier-Stokes equations?