|
[Sponsors] |
Run time post processing is crashing pimpleDyMFoam for a plunging airfoil |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 8, 2014, 04:51 |
Run time post processing is crashing pimpleDyMFoam for a plunging airfoil
|
#1 |
Member
Pruthvi
Join Date: Feb 2014
Posts: 41
Rep Power: 12 |
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 |
|
August 31, 2015, 14:44 |
|
#2 |
Senior Member
Join Date: Jan 2013
Posts: 135
Rep Power: 13 |
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? |
|
September 4, 2015, 01:31 |
|
#3 |
Member
Pruthvi
Join Date: Feb 2014
Posts: 41
Rep Power: 12 |
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. |
|
Tags |
forcecoeffs, openfoam2.3.x, pimpledymfoam, post procesing, run time |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SU2 AOA optimization | 454514566@qq.com | SU2 | 9 | March 7, 2022 17:17 |
How to export time series of variables for one point? | mary mor | OpenFOAM Post-Processing | 8 | July 19, 2017 11:54 |
Star cd es-ice solver error | ernarasimman | STAR-CD | 2 | September 12, 2014 01:01 |
Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 07:20 |
Unstabil Simulation with chtMultiRegionFoam | mbay101 | OpenFOAM Running, Solving & CFD | 13 | December 28, 2013 14:12 |