kinbean |
June 21, 2018 04:55 |
Dynamic mesh Vector mesh error
Hey, I am have a bit of trouble running a simulation with a dynamic mesh where the vector field for movement is defined. It was working fine with oscillating linear and rotating motion. But not with tabulated6DoFMotion. It runs for a couple hours then throws an error message before completing the time the vectors are defined for.
This is the error message
Code:
solidBodyMotionFunctions::tabulated6DoFMotion::transformation(): Time = 0.00964224 transformation: ((0.00453319 0.0616797 0) (0.95434 (0 0 0.298724)))
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.72875e-06, No Iterations 146
DICPCG: Solving for pcorr, Initial residual = 0.00647977, Final residual = 9.61029e-06, No Iterations 92
time step continuity errors : sum local = 5.32607e-09, global = -3.10113e-16, cumulative = 7.40874e-16
PIMPLE: iteration 1
DICPCG: Solving for p, Initial residual = 0.898125, Final residual = 0.044323, No Iterations 17
DICPCG: Solving for p, Initial residual = 0.0240248, Final residual = 0.00119808, No Iterations 32
time step continuity errors : sum local = 5.02539e-10, global = -3.10115e-16, cumulative = 4.30759e-16
DICPCG: Solving for p, Initial residual = 0.502555, Final residual = 0.0247756, No Iterations 17
DICPCG: Solving for p, Initial residual = 0.0168104, Final residual = 0.000839935, No Iterations 32
time step continuity errors : sum local = 3.59747e-10, global = -3.10115e-16, cumulative = 1.20645e-16
PIMPLE: iteration 2
DICPCG: Solving for p, Initial residual = 0.729046, Final residual = 0.0351105, No Iterations 27
DICPCG: Solving for p, Initial residual = 0.0886994, Final residual = 0.00392495, No Iterations 38
time step continuity errors : sum local = 4.75911e-10, global = -3.10114e-16, cumulative = -1.89469e-16
DICPCG: Solving for p, Initial residual = 0.00881001, Final residual = 0.000431365, No Iterations 78
DICPCG: Solving for p, Initial residual = 0.000629799, Final residual = 9.7113e-06, No Iterations 91
time step continuity errors : sum local = 1.03688e-12, global = -3.10114e-16, cumulative = -4.99583e-16
DILUPBiCG: Solving for nuTilda, Initial residual = 2.43288e-05, Final residual = 4.74751e-10, No Iterations 3
bounding nuTilda, min: -2.48209e-06 max: 0.29104 average: 0.0287184
PIMPLE: converged in 2 iterations
ExecutionTime = 537.26 s ClockTime = 547 s
forceCoeffs forceCoeffs output:
Cm = 44282.2
Cd = -25037.2
Cl = 36030.6
Cl(f) = 62297.5
Cl(r) = -26266.9
fieldMinMax fieldMinMax output:
min(mag(U)) = 0.00194439 at location (0.029796 0.0744064 0.00381) on processor 0
max(mag(U)) = 35.1152 at location (0.004807 0.0610662 0.00381) on processor 0
min(p) = -1665.48 at location (0.04586 -0.803495 0.00381) on processor 3
max(p) = 1228.65 at location (0.837063 0.508323 0.00381) on processor 3
Courant Number mean: 0.000629294 max: 0.899674
deltaT = 3.31752e-06
Time = 0.0096455617
solidBodyMotionFunctions::tabulated6DoFMotion::transformation(): Time = 0.00964556 transformation: ((0.00453323 0.0616796 0) (0.954339 (0 0 0.298725)))
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.38253e-06, No Iterations 146
DICPCG: Solving for pcorr, Initial residual = 0.00637304, Final residual = 9.45022e-06, No Iterations 92
time step continuity errors : sum local = 4.62589e-09, global = 5.28549e-16, cumulative = 2.89663e-17
PIMPLE: iteration 1
DICPCG: Solving for p, Initial residual = 0.9042, Final residual = 0.042896, No Iterations 26
DICPCG: Solving for p, Initial residual = 0.0234949, Final residual = 0.00105508, No Iterations 38
time step continuity errors : sum local = 4.41989e-10, global = 5.28549e-16, cumulative = 5.57515e-16
DICPCG: Solving for p, Initial residual = 0.514141, Final residual = 0.0244435, No Iterations 26
DICPCG: Solving for p, Initial residual = 0.0164958, Final residual = 0.000741193, No Iterations 38
time step continuity errors : sum local = 3.10885e-10, global = 5.28549e-16, cumulative = 1.08606e-15
PIMPLE: iteration 2
DICPCG: Solving for p, Initial residual = 0.721387, Final residual = 0.0332969, No Iterations 23
DICPCG: Solving for p, Initial residual = 0.0848257, Final residual = 0.00383372, No Iterations 31
time step continuity errors : sum local = 4.85069e-10, global = 5.2855e-16, cumulative = 1.61461e-15
DICPCG: Solving for p, Initial residual = 0.0090906, Final residual = 0.00039283, No Iterations 84
DICPCG: Solving for p, Initial residual = 0.000784486, Final residual = 9.5523e-06, No Iterations 91
time step continuity errors : sum local = 1.02909e-12, global = 5.28549e-16, cumulative = 2.14316e-15
DILUPBiCG: Solving for nuTilda, Initial residual = 2.43002e-05, Final residual = 4.73755e-10, No Iterations 3
bounding nuTilda, min: -4.80249e-06 max: 0.290994 average: 0.0287156
PIMPLE: converged in 2 iterations
ExecutionTime = 537.34 s ClockTime = 547 s
forceCoeffs forceCoeffs output:
Cm = 49281.2
Cd = -27272.9
Cl = 41439.7
Cl(f) = 70001.1
Cl(r) = -28561.4
fieldMinMax fieldMinMax output:
min(mag(U)) = 0.00194445 at location (0.029796 0.0744065 0.00381) on processor 0
max(mag(U)) = 35.1136 at location (0.00480706 0.0610661 0.00381) on processor 0
min(p) = -561.208 at location (0.837061 0.508326 0.00381) on processor 3
max(p) = 2084.84 at location (0.0458629 -0.803495 0.00381) on processor 3
Courant Number mean: 0.000629519 max: 0.900107
deltaT = 3.3171e-06
Time = 0.0096488788
solidBodyMotionFunctions::tabulated6DoFMotion::transformation(): Time = 0.00964888 transformation: ((0.00453328 0.0616795 0) (0.954339 (0 0 0.298727)))
^C
parallels@ubuntu:~/Desktop/casefixvect$ [3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::pow3(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[3] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow3<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[3] #5 Foam::RASModels::SpalartAllmaras<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::fv1(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[3] #6 Foam::RASModels::SpalartAllmaras<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[3] #7 ? at ??:?
[3] #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #9 ? at ??:?
[ubuntu:11141] *** Process received signal ***
[ubuntu:11141] Signal: Floating point exception (8)
[ubuntu:11141] Signal code: (-6)
[ubuntu:11141] Failing at address: 0x3e800002b85
[ubuntu:11141] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7f5577282cb0]
[ubuntu:11141] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f5577282c37]
[ubuntu:11141] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7f5577282cb0]
[ubuntu:11141] [ 3] /opt/openfoam30/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4pow3ERNS_5FieldIdEERKNS_5UListIdEE+0x56) [0x7f55785bcb16]
[ubuntu:11141] [ 4] /opt/openfoam30/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam4pow3INS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_+0x153) [0x7f557bcc0bf3]
[ubuntu:11141] [ 5] /opt/openfoam30/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZNK4Foam9RASModels15SpalartAllmarasINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE3fv1ERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x2b) [0x7f557bcc0e2b]
[ubuntu:11141] [ 6] /opt/openfoam30/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels15SpalartAllmarasINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0xad) [0x7f557bd0fedd]
[ubuntu:11141] [ 7] pimpleDyMFoam() [0x420689]
[ubuntu:11141] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f557726df45]
[ubuntu:11141] [ 9] pimpleDyMFoam() [0x421b7e]
[ubuntu:11141] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 11141 on node ubuntu exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
This is my vector data
Code:
51 (
( 0 ( ( 0 0.0762 0 ) ( 0 0 34.2014591191 ) ) )
( 0.1428571429 ( ( 0 0.0755991402 0 ) ( 0 0 42.4843555122 ) ) )
( 0.2857142857 ( ( 0 0.0738060369 0 ) ( 0 0 50.09724824 ) ) )
( 0.4285714286 ( ( 0 0.0708489682 0 ) ( 0 0 56.9200774363 ) ) )
( 0.5714285714 ( ( 0 0.066774569 0 ) ( 0 0 62.8452430089 ) ) )
( 0.7142857143 ( ( 0 0.061647095 0 ) ( 0 0 67.7793015574 ) ) )
( 0.8571428572 ( ( 0 0.0555474094 0 ) ( 0 0 71.6444400309 ) ) )
( 1 ( ( 0 0.048571708 0 ) ( 0 0 74.3797028867 ) ) )
( 1.1428571429 ( ( 0 0.0408300018 0 ) ( 0 0 75.9419533957 ) ) )
( 1.2857142857 ( ( 0 0.032444382 0 ) ( 0 0 76.306553934 ) ) )
( 1.4285714286 ( ( 0 0.023547095 0 ) ( 0 0 75.4677545335 ) ) )
( 1.5714285715 ( ( 0 0.0142784562 0 ) ( 0 0 73.4387835616 ) ) )
( 1.7142857143 ( ( 0 0.0047846376 0 ) ( 0 0 70.2516391028 ) ) )
( 1.8571428572 ( ( 0 -0.0047846376 0 ) ( 0 0 65.9565843289 ) ) )
( 2 ( ( 0 -0.0142784562 0 ) ( 0 0 60.6213548196 ) ) )
( 2.1428571429 ( ( 0 -0.023547095 0 ) ( 0 0 54.3300903313 ) ) )
( 2.2857142858 ( ( 0 -0.032444382 0 ) ( 0 0 47.1820078632 ) ) )
( 2.4285714286 ( ( 0 -0.0408300018 0 ) ( 0 0 39.2898369459 ) ) )
( 2.5714285715 ( ( 0 -0.048571708 0 ) ( 0 0 30.7780418294 ) ) )
( 2.7142857143 ( ( 0 -0.0555474094 0 ) ( 0 0 21.7808586072 ) ) )
( 2.8571428572 ( ( 0 -0.061647095 0 ) ( 0 0 12.4401782336 ) ) )
( 3.0000000001 ( ( 0 -0.066774569 0 ) ( 0 0 2.9033088178 ) ) )
( 3.1428571429 ( ( 0 -0.0708489682 0 ) ( 0 0 -6.6793475124 ) ) )
( 3.2857142858 ( ( 0 -0.0738060369 0 ) ( 0 0 -16.1566665423 ) ) )
( 3.4285714286 ( ( 0 -0.0755991402 0 ) ( 0 0 -25.3791852893 ) ) )
( 3.5714285715 ( ( 0 -0.0762 0 ) ( 0 0 -34.2014591234 ) ) )
( 3.7142857144 ( ( 0 -0.0755991402 0 ) ( 0 0 -42.4843555162 ) ) )
( 3.8571428572 ( ( 0 -0.0738060369 0 ) ( 0 0 -50.0972482436 ) ) )
( 4.0000000001 ( ( 0 -0.0708489682 0 ) ( 0 0 -56.9200774395 ) ) )
( 4.1428571429 ( ( 0 -0.066774569 0 ) ( 0 0 -62.8452430117 ) ) )
( 4.2857142858 ( ( 0 -0.061647095 0 ) ( 0 0 -67.7793015596 ) ) )
( 4.4285714287 ( ( 0 -0.0555474094 0 ) ( 0 0 -71.6444400325 ) ) )
( 4.5714285715 ( ( 0 -0.048571708 0 ) ( 0 0 -74.3797028878 ) ) )
( 4.7142857144 ( ( 0 -0.0408300018 0 ) ( 0 0 -75.9419533962 ) ) )
( 4.8571428572 ( ( 0 -0.032444382 0 ) ( 0 0 -76.3065539339 ) ) )
( 5.0000000001 ( ( 0 -0.023547095 0 ) ( 0 0 -75.4677545328 ) ) )
( 5.142857143 ( ( 0 -0.0142784562 0 ) ( 0 0 -73.4387835603 ) ) )
( 5.2857142858 ( ( 0 -0.0047846376 0 ) ( 0 0 -70.2516391009 ) ) )
( 5.4285714287 ( ( 0 0.0047846376 0 ) ( 0 0 -65.9565843265 ) ) )
( 5.5714285715 ( ( 0 0.0142784562 0 ) ( 0 0 -60.6213548167 ) ) )
( 5.7142857144 ( ( 0 0.023547095 0 ) ( 0 0 -54.3300903279 ) ) )
( 5.8571428573 ( ( 0 0.032444382 0 ) ( 0 0 -47.1820078594 ) ) )
( 6.0000000001 ( ( 0 0.0408300018 0 ) ( 0 0 -39.2898369418 ) ) )
( 6.142857143 ( ( 0 0.048571708 0 ) ( 0 0 -30.778041825 ) ) )
( 6.2857142858 ( ( 0 0.0555474094 0 ) ( 0 0 -21.7808586026 ) ) )
( 6.4285714287 ( ( 0 0.061647095 0 ) ( 0 0 -12.4401782288 ) ) )
( 6.5714285716 ( ( 0 0.066774569 0 ) ( 0 0 -2.903308813 ) ) )
( 6.7142857144 ( ( 0 0.0708489682 0 ) ( 0 0 6.6793475172 ) ) )
( 6.8571428573 ( ( 0 0.0738060369 0 ) ( 0 0 16.156666547 ) ) )
( 7.0000000001 ( ( 0 0.0755991402 0 ) ( 0 0 25.3791852938 ) ) )
( 7.142857143 ( ( 0 0.0762 0 ) ( 0 0 34.2014591277 ) ) )
)
Dynamic Mesh file
Code:
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dynamicFvMesh solidBodyMotionFvMesh;
solidBodyMotionFvMeshCoeffs
{
solidBodyMotionFunction tabulated6DoFMotion;
tabulated6DoFMotionCoeffs
{
CofG (0.0254 0 0 );
timeDataFileName "/home/parallels/Desktop/casefixvect/constant/datafile.dat";
}
}
// ************************************************************************* //
Any help would be appreciated.
|