# Dynamic mesh Vector mesh error

 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 21, 2018, 04:55 Dynamic mesh Vector mesh error #1 New Member   Join Date: Jun 2017 Posts: 11 Rep Power: 8 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&, Foam::UList const&) at ??:? [3] #4 Foam::tmp > Foam::pow3(Foam::GeometricField const&) at ??:? [3] #5 Foam::RASModels::SpalartAllmaras >::fv1(Foam::GeometricField const&) const at ??:? [3] #6 Foam::RASModels::SpalartAllmaras >::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.

 June 21, 2018, 09:39 #2 New Member   Join Date: Jun 2017 Posts: 11 Rep Power: 8 **update* My flow field is crazy. The magnitude is reasonable but the direction is wack. So I thought it might have something to do with the boundary conditions. I've used freestream since it is a 2D model. Once again, Any help would be super appreciated.

 June 21, 2018, 10:21 #3 New Member   Join Date: Jun 2017 Posts: 11 Rep Power: 8 As an additional question is it possible to add a phase shift between 2 different types of motion in the multimotion parameters for the dynamic mesh?