CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (http://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   dsmcFoam (http://www.cfd-online.com/Forums/openfoam-pre-processing/109811-dsmcfoam.html)

dark lancer November 27, 2012 09:21

dsmcFoam
 
hi
I try to run a problem with dsmc and mesh my geometry in gambit then convert to OF.
I excute the checkMesh All of them its Ok but at the end line :

Quote:

moosaie@:~/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical> checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : checkMesh
Date : Nov 17 2012
Time : 12:53:01
Host : ".site"
PID : 6467
Case : /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 1489344
faces: 4387940
internal faces: 4308820
cells: 1449460
boundary patches: 5
point zones: 0
face zones: 1
cell zones: 1

Overall number of cells of each type:
hexahedra: 1449460
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
body 17480 17694 ok (non-closed singly connected)
symmetry 34040 34576 ok (non-closed singly connected)
outlet_from_side 13524 13764 ok (non-closed singly connected)
outlet 10902 11079 ok (non-closed singly connected)
inlet 3174 3267 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-262.5 -262.5 -0.4142099593) (1375 262.5 262.5)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (1.756589626e-16 1.121895293e-15 5.86943899e-15) OK.
Max cell openness = 1.916326799e-15 OK.
Max aspect ratio = 61.10879384 OK.
Minumum face area = 0.01039895189. Maximum face area = 280.93608. Face area magnitudes OK.
Min volume = 0.001605026005. Max volume = 2372.742765. Total volume = 167688127. Cell volumes OK.
Mesh non-orthogonality Max: 86.44406365 average: 27.94667529
*Number of severely non-orthogonal faces: 20108.
Non-orthogonality check OK.
<<Writing 20108 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 4.752768679, 879 highly skew faces detected which may impair the quality of the results
<<Writing 879 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End

and I use dsmcFoam I see this Error:
Quote:

moosaie@:~/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical> dsmcFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : dsmcFoam
Date : Nov 17 2012
Time : 12:37:32
Host : ".site"
PID : 6370
Case : /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Constructing dsmcCloud


--> FOAM FATAL IO ERROR:
cannot find file

file: /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical/0/dsmcSigmaTcRMax at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting

moosaie@:~/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical>
I use this order:

Quote:

./Allrun
I see this error:
Quote:

moosaie@:~/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical> ./Allrun
Running blockMesh on /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical
Running decomposePar on /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical
Running dsmcInitialise in parallel on /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical using 4 processes
Running dsmcFoam in parallel on /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical using 4 processes
Running reconstructPar on /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical

and I don't have the blockMesh
.

for example when run the cavity
I type this order and now that I don't have the blockMesh what should I do?
cd
blockMesh
icoFoam
paraFoam

is this correct:
cd
.....????
dsmcFoam
paraFoam

nimasam November 30, 2012 02:43

dear hadi
you mixedup several things
1- your import mesh has a skewness problem
2- dsmcSigmaTcRMax can not be find in your 0 folder
3-Allrun executes commands which are written in it successively,
so i suggest to type command manually, or edit your Allrun file

dark lancer December 15, 2012 03:25

dear nimasam
I solved my problem mesh.
I have a bit confused.
now for run my problem
I have one pc for run my problem then I don;t excute (./Allrun) in my terminal because this order run in parallel,is it correct?
and after check mesh I type this in my terminal
Quote:

moosaie@:~/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical> dsmcFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : dsmcFoam
Date : Nov 17 2012
Time : 12:37:32
Host : ".site"
PID : 6370
Case : /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Constructing dsmcCloud


--> FOAM FATAL IO ERROR:
cannot find file

file: /home/moosaie/OpenFOAM/OpenFOAM-2.1.1/tutorials/discreteMethods/dsmcFoam/hemispherical/0/dsmcSigmaTcRMax at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting
I use openFOAM 2.1.1 and in the 0 folder of dsmc code don't exist the dsmcSigmaTcRMax

in my 0 folder
bounder T
boundery U
fD
iDof
q
internal KE
internalE
mommentum
RhoN
RhoM
dsmcRhoN

nimasam December 15, 2012 03:30

Dear hadi
as you correctly point there is no dsmcSigmaTcRMax file in folder 0, so you can add it to folder 0 :)

dark lancer December 15, 2012 03:35

thanks
where can I give it?

nimasam December 15, 2012 06:32

copy paste another file and rename it
"dsmcSigmaTcRMax"
then you need to assign appropriate dimension to it ;)

dark lancer December 16, 2012 14:26

thanks nimasam

I have a problem for run my case.we run a code that create with fortran and we decide run my case together because we use one part of cpu for run fortran but this error can't let us:
Quote:

error while loading shared libraries: libopen-rte.so.0: cannot open shared object file: No such file or directory
I do this
Quote:

export LD_LIBRARY_PATH=/usr/lib64/mpi/gcc/openmpi/lib64:$LD_LIBRARY_PATH
but not answer

nimasam December 16, 2012 15:35

may be you need to add this lib in your controldict to make it work,
but i suggest you to open new post and declare your problem there


All times are GMT -4. The time now is 12:52.