Jalued |
July 16, 2014 05:52 |
Arithmetic exception
2 Attachment(s)
For some reason, I am continuing to receive an "Arithmetric exception" error message when I attempt to run any 2D simulation after importing from Pointwise. In order to find the culprit I have gone to the simplest case I could think of: a 2D simulation that mimics the 2DAirfoil tutorial but without the airfoil.
I used Pointwise to construct the parameters which consist of a box that was extruded by 0.05m and given a simple structured mesh, passing checkMesh with flying colours:
PHP Code:
Checking geometry... Overall domain bounding box (1 1 0) (101 101 0.05) Mesh (non-empty, non-wedge) directions (1 1 0) Mesh (non-empty) directions (1 1 0) All edges aligned with or perpendicular to non-empty directions. Boundary openness (-3.46598e-21 3.46598e-21 9.19282e-15) OK. Max cell openness = 1.08813e-16 OK. Max aspect ratio = 1 OK. Minimum face area = 0.0505051. Maximum face area = 1.0203. Face area magnitudes OK. Min volume = 0.0510152. Max volume = 0.0510152. Total volume = 500. Cell volumes OK. Mesh non-orthogonality Max: 0 average: 0 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.01608e-13 OK. Coupled point location match (average 0) OK.
Mesh OK.
The 0 and system folders are identical to the 2DAirfoil except for the removal of the wall.
I ran the code through gdb simpleFoam and it exploded within the first iteration:
PHP Code:
smoothSolver: Solving for Ux, Initial residual = 0.343067, Final residual = 1.85346e+13, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.00320924, Final residual = 2.19375e+11, No Iterations 1000 GAMG: Solving for p, Initial residual = 1, Final residual = 5.19741e+118, No Iterations 1000 time step continuity errors : sum local = 2.01405e+117, global = 3.29242e+101, cumulative = 3.29242e+101
Program received signal SIGFPE, Arithmetic exception. 0x00007ffff4ab9575 in Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) () from /opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so
I have no idea what could be wrong with this, as the 2DAirfoil example works perfectly and a similar "laminar flow through a straight pipe" Pointwise to OpenFoam tutorial also works with no issues.
I was unable to attach the pw file or constant folder due to space constraints
Thanks,
Ed
|