CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   DecomposeParDict (http://www.cfd-online.com/Forums/openfoam-solving/60280-decomposepardict.html)

hsieh March 20, 2006 16:18

Hi, I have a workstation wi
 
Hi,

I have a workstation with 2 CPUs. I ran the decomposePar utility to decompose the case into two. In the original time=0 directory, gamma field has been initialized to have non-uniform distribution.

After running the decomposePar utility, I got processor0 and processor1. However, the gamma field on processor0/0 directory does not have non-uniform gamma field. Can someone suggest a solution? Thanks!

pei

gschaider March 20, 2006 17:05

Have you checked how the domai
 
Have you checked how the domain has been decomposed? Maybe in processor0's part the gamma-field IS uniform. For instance in the famous damBreak case this would be the case if processor0 gets the right part.

hsieh March 21, 2006 10:15

Hi, Bernhard, Thanks for th
 
Hi, Bernhard,

Thanks for the reply. You are 100% correct. I found a problem (sort of). I compiled a modified version of rasInterFoam, called rasInterFoamInlet. I have no problem running the serial version of rasInterFoamInlet, however, when I ran rasInterFoamInlet in parallel, I got the following error:

phsieh@sun:~/OpenFOAM/phsieh-1.2/run/tutorials/rasInterFoam> mpirun -np 2 rasInterFoamInlet . CP-probeB -parallel
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.2 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : rasInterFoamInlet . CP-probeB -parallel
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.2 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : rasInterFoamInlet . CP-probeB -parallel
[1] Date : Mar 21 2006
[0] Date : Mar 21 2006
[0] Time : 09:30:22
[0] Host : sun
[0] PID : 8997
[1] Time : 09:30:22
[1] Host : sun
[1] PID : 8998
[1] Root : /home/phsieh/OpenFOAM/phsieh-1.2/run/tutorials/rasInterFoam
[0] Root : /home/phsieh/OpenFOAM/phsieh-1.2/run/tutorials/rasInterFoam
[0] Case : CP-probeB
[0] Nprocs : 2
[0] Slaves :
[0] 1
[0] (
[0] sun.8998
[0] )
[0]
Create time

[1] Case : CP-probeB
[1] Nprocs : 2
Create mesh, no clear-out for time = 0


Reading environmentalProperties
Reading field pd

Reading field gamma

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Calculating field g.h

Selecting turbulence model kEpsilon
Mean and max Courant Numbers = 0 0.02968
ICCG: Solving for pcorr, Initial residual = 1, Final residual = 9.92634e-11, No Iterations 205
time step continuity errors : sum local = 6.33051e-16, global = -6.8148e-17, cumulative = -6.8148e-17

Starting time loop

Mean and max Courant Numbers = 0.000532038 0.845396
deltaT = 1.5873e-06
Time = 1.5873e-06

-----------------------------------------------------------------------------
One of the processes started by mpirun has exited with a nonzero exit
code. This typically indicates that the process finished in error.
If your process did not finish in error, be sure to include a "return
0" or "exit(0)" in your C code before exiting the application.

PID 8997 failed on node n0 (192.168.0.10) due to signal 11.
-----------------------------------------------------------------------------
phsieh@sun:~/OpenFOAM/phsieh-1.2/run/tutorials/rasInterFoam>

******************
The following is the section I added to the original rasInterFoam:
------

// Get index of patch
label inletPatchID = mesh.boundaryMesh().findPatchID("inlet");

// Get reference to boundary value
vectorField& inletU = U.boundaryField()[inletPatchID];

// Manipulate value
velInlet = Foam::sin(3.14159/0.3158 * runTime.value());
velInlet = -150.0/0.3158/160.0 * velInlet;
if(runTime.value() >= 0.3158) velInlet=0.0;

inletU = vector(velInlet,0,0);
Info<< "inletU:" << min(inletU) << endl;

**********************

I have no problem running the original rasInterFoam in parallel.

Any suggestion will be highly appreciated.

Pei


All times are GMT -4. The time now is 03:26.