[pruthvi1991]

Hey guys. I'm using openfoam230 to simulate a plunging airfoil. The simulation runs perfectly fine without the forceCoeffs function. The moment I add the function, it crashes. I get different error messages for different number of cores. Please help me out. For a certain case I managed to run the simulation for 3 cores. It crashed for all other choices including 1 core. I'm unable to figure out what is causing this error.

I'm giving the link to my github account where the case file exists exactly as I tried to run it. https://github.com/pruthvi1991/pimpl.../onemeterplate

CheckMesh log file :

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : checkMesh
Date : Oct 08 2014
Time : 03:24:59
Host : "ubuntu"
PID : 30815
Case : /home/jujja/GIT1/pimpleDyMFoam-Tutorials/onemeterplate/snappy
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: 702720
faces: 1924500
internal faces: 1760223
cells: 611384
faces per cell: 6.02686
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 605816
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 5568
Breakdown of polyhedra by number of faces:
faces number of cells
6 184
9 5296
12 87
15 1

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
topAndBottom 160 324 ok (non-closed singly connected)
inlet 962 1153 ok (non-closed singly connected)
outlet 962 1153 ok (non-closed singly connected)
symFront 78407 79395 ok (non-closed singly connected)
symBack 78410 79398 ok (non-closed singly connected)
wing 5376 5544 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-4 -6 -0.1) (14 6 0.1)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (3.98309e-17 -3.13155e-18 -9.03582e-16) OK.
Max cell openness = 3.69726e-16 OK.
Max aspect ratio = 5.59099 OK.
Minimum face area = 4.43241e-06. Maximum face area = 0.08302. Face area magnitudes OK.
Min volume = 3.45928e-08. Max volume = 0.0124572. Total volume = 43.1991. Cell volumes OK.
Mesh non-orthogonality Max: 45.6781 average: 3.68824
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.80463 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

Here are some pics of the Mesh





controlDict file:

libs
(
"libOpenFOAM.so"
);
application pimpleDyMFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 3000;

deltaT 1e-1;

writeControl adjustableRunTime;

writeInterval 4.1665;

purgeWrite 0;

writeFormat binary;

writePrecision 10;

writeCompression on;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 0.9;
functions
{
#include "forceCoeffs"
}

forceCoeffs file:

forceCoeffs
{
// type forceCoeffs;
type forces;

functionObjectLibs ( "libforces.so" );

outputControl timeStep;
timeInterval 1;

log yes;

patches ( "wing.*" );
pName p;
UName U;
rhoName rhoInf; // Indicates incompressible
log true;
rhoInf 1; // Redundant for incompressible
liftDir (0 1 0);
dragDir (1 0 0);
CofR (0.2285 0 0); // Axle midpoint on ground
pitchAxis (0 0 1);
magUInf 0.03854;
lRef 1; // Wheelbase length
Aref 0.0457; // Estimated

binData
{
nBin 20; // output data into 20 bins
direction (1 0 0); // bin direction
cumulative yes;
}
}

To prevent the post from getting too long I have attached the error logs for 1 and 8 cores. Please check them and help me fix this problem. Thanks. http://www.cfd-online.com/Forums/att...1&d=1412753744

[kkpal]

I also have this problem when I put binData in the controlDict when running pimpleDyMFoam.
binData works fine in the static cases where pimpleFoam is used.
Have you solved this problem?

[pruthvi1991]

Hello kai,

Which version of OpenFOAM are you using? This error is a bug fixed by a commit. If you have 2.3.X you should be fine. Just update your OpenFOAM if its older. Let me know if the bug persists.

Thanks,
Pru.