[alireza.glz]

Hi Dear foamers
I'm a new user of openfoam and I'm trying to simulate a pipe through which a fluid passes in laminar flow and is being heated.
for this problem I'm using buoyantboussinesqsimplefoam solver and my RASProperties file is as follows:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RASModel laminar;

turbulence off;

printCoeffs on;
/////////////////////////////////////////////////

but openfoam does not run my case and following errors appear:

Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Reading field T

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model laminar
Calculating field beta*(g.h)


Starting time loop

Time = 1



LHS and RHS of - have different dimensions
dimensions : [0 3 -2 0 0 0 0] - [1 0 -2 0 0 0 0]
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam:perator-(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#3 Foam::tmp<Foam::GeometricField<Foam::typeOfSum<dou ble, double>::type, Foam::fvsPatchField, Foam::surfaceMesh> > Foam:perator-<double, double, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField< double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantBoussinesqSimpleFoam"
#4 main in "/home/alireza/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantBoussinesqSimpleFoam"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122


From function operator-(const dimensionSet& ds1, const dimensionSet& ds2)
in file dimensionSet/dimensionSet.C at line 423.

FOAM aborting

Aborted
/////////////////////////////////////////////////////

unfortunately I don't understand the source of these errors.
could anybody here help me to correct my case?
Please give me your valuable advises, please!

[sixwp]

Hi Alireza,

the error tells you that you have wrongly defined the boundary unit.

if you compare:
[0 3 -2 0 0 0 0] with [1 0 -2 0 0 0 0]
m3.s-2 with kg.s-2 => therefore a difference in kg/m3

My guess is:
I bet your p/p_rgh initial value are define with [1 -1 -2 0 0 0 0]
Change it to: [0 2 -2 0 0 0 0] and you'll be fine.

The reason: depending the solver, density can be included or not

Best regards,

[alireza.glz]

Dear samuel
Thanks a lot
You are right
I made the changes and my case run
Thanks a milion for your help

[sixwp]

you're most welcome

[libia87]

Hi!

I have a similar simulation, and when I run the problem, an error ocurrs:

--> FOAM FATAL ERROR:
Invalid wall function specification
Patch type for patch top must be wall
Current patch type is patch



From function kappatJayatillekeWallFunctionFvPatchScalarField::c heckType()
in file derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C at line 56.

FOAM aborting


I change the geometry and everythings it's ok, but the error says that the surface "top" it's not a wall, and I don't know what I have to do,

some ideas??

[nimasam]

Invalid wall function specification
Patch type for patch top must be wall
Current patch type is patch

patch type is patch so change it into wall, for changing the patch type you need edit your blockMeshDict or if you import a file from gambit change your boundary file
look at constant/polyMesh

P.S. if you read the error carefully, the solution is in it !!!

[libia87]

Originally Posted by nimasam

patch type is patch so change it into wall, for changing the patch type you need edit your blockMeshDict or if you import a file from gambit change your boundary file
look at constant/polyMesh

P.S. if you read the error carefully, the solution is in it !!!

Ok, thanks!! I'm new with openfoam and it's difficult to me understand all the warnings that the program gives me. Thank you again! my simulation runs!

Another question, If I run my simulation and then I change some parameters and I run it again, paraFoam gives me the same results as the first run. How can I restart the simulation to obtain the new results? thanks!! (I change the writeInterval, but in my directory appears the same folders as before).