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

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.

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