|
[Sponsors] |
April 19, 2021, 14:48 |
Error running AMI propeller simulation
|
#1 |
New Member
Join Date: Apr 2021
Posts: 3
Rep Power: 5 |
Hi Foamers!
I need some help with a propeller simulation (transient). As a quick note, I'm quite new to OpenFoam. A brief description of the project: I created the propeller geometry with Solidworks and then created the enclosure and the mesh with Salome. There are two regions: the main enclosure (that doesn't rotate) which represents the freestream flow and the propeller enclosure (which should rotate) that contains the propeller. The meshes of both domains are created separately and then merged using mergeMeshes. The rotor cellZone region is then set using setSet. I also use topoSetDict to add all AMI faces to a faceSet. Finally, due to the size of the simulation, I run it in parallel using 20 procs (using the singleProcessorFaceSets to keep all AMI faces on the same processor). The simulation blows after a few iterations. I must note that the simulation blows when it's run in series too. I've checked the mesh and it's okay (see attached log file). This is the log of a failed parallel simulation. Code:
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 4.1 Exec : pimpleDyMFoam -parallel Date : Apr 19 2021 Time : 15:56:34 Host : "93a46648863c" PID : 165 Case : /home/openfoam nProcs : 20 Slaves : 19 ( "93a46648863c.166" "93a46648863c.167" "93a46648863c.168" "93a46648863c.169" "93a46648863c.170" "93a46648863c.171" "93a46648863c.172" "93a46648863c.173" "93a46648863c.174" "93a46648863c.175" "93a46648863c.176" "93a46648863c.177" "93a46648863c.178" "93a46648863c.179" "93a46648863c.180" "93a46648863c.181" "93a46648863c.182" "93a46648863c.183" "93a46648863c.184" ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone rotor PIMPLE: no residual control data found. Calculations will employ 2 corrector loops Reading field p Reading field U Reading/calculating face flux field phi AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998918, 1.0014, 1 AMI: Patch target sum(weights) min/max/average = 0.998764, 1.00002, 0.999999 Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 -0.33; sigmak 1; sigmaEps 1.3; } No MRF models present Reading/calculating face velocity Uf No finite volume options present Courant Number mean: 3.12403e-06 max: 0.000570603 Starting time loop Courant Number mean: 3.0931e-06 max: 0.000564953 deltaT = 1.17647e-05 Time = 1.17647e-05 AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1 AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00635745, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.369898, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.36776, No Iterations 1000 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0071022, No Iterations 11 time step continuity errors : sum local = 5.33589e-08, global = -3.80218e-08, cumulative = -3.80218e-08 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.47015, Final residual = 0.0836701, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.265834, Final residual = 0.0436782, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.10333, Final residual = 0.00688327, No Iterations 1000 GAMG: Solving for p, Initial residual = 0.0941766, Final residual = 8.35714e-07, No Iterations 22 time step continuity errors : sum local = 5.98351e-09, global = 3.17688e-09, cumulative = -3.48449e-08 smoothSolver: Solving for epsilon, Initial residual = 0.994799, Final residual = 0.00754865, No Iterations 1000 smoothSolver: Solving for k, Initial residual = 1, Final residual = 2.46881e-06, No Iterations 1000 ExecutionTime = 857.17 s ClockTime = 866 s Courant Number mean: 0.00943062 max: 83392.3 deltaT = 2.82153e-10 Time = 1.1765e-05 AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1 AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.7144, Final residual = 0.0433931, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.775929, Final residual = 0.0502612, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.82386, Final residual = 0.0560173, No Iterations 2 GAMG: Solving for p, Initial residual = 0.910335, Final residual = 0.00880902, No Iterations 4 time step continuity errors : sum local = 3.13285e-10, global = 2.682e-12, cumulative = -3.48422e-08 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.530458, Final residual = 0.472351, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.437057, Final residual = 0.265858, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.371612, Final residual = 0.0227143, No Iterations 1000 GAMG: Solving for p, Initial residual = 1.58002e-08, Final residual = 1.58002e-08, No Iterations 0 time step continuity errors : sum local = 6.73282e-06, global = 3.84411e-12, cumulative = -3.48384e-08 smoothSolver: Solving for epsilon, Initial residual = 0.192127, Final residual = 7.48896e-07, No Iterations 4 smoothSolver: Solving for k, Initial residual = 0.827074, Final residual = 0.106168, No Iterations 1000 ExecutionTime = 1320.87 s ClockTime = 1331 s Courant Number mean: 0.000144388 max: 1942.94 deltaT = 2.9044e-13 Time = 1.1765e-05 AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1 AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.768038, Final residual = 0.0766915, No Iterations 118 smoothSolver: Solving for Uy, Initial residual = 0.81215, Final residual = 0.0808922, No Iterations 127 smoothSolver: Solving for Uz, Initial residual = 0.902063, Final residual = 0.0894033, No Iterations 22 GAMG: Solving for p, Initial residual = 0.880168, Final residual = 0.00271608, No Iterations 2 time step continuity errors : sum local = 0.000258367, global = 4.23745e-07, cumulative = 3.88907e-07 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.580994, Final residual = 0.0598996, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.720866, Final residual = 0.1011, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.615497, Final residual = 0.137957, No Iterations 1000 GAMG: Solving for p, Initial residual = 0.000621124, Final residual = 9.73629e-07, No Iterations 22 time step continuity errors : sum local = 4.05813e-06, global = -1.78319e-07, cumulative = 2.10588e-07 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.106181, No Iterations 1000 smoothSolver: Solving for k, Initial residual = 0.999862, Final residual = 0.154266, No Iterations 1000 ExecutionTime = 1944.37 s ClockTime = 1959 s Courant Number mean: 0.00645368 max: 6.77553e+07 deltaT = 8.5732e-21 Time = 1.1765e-05 AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1 AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.739661, Final residual = 0.0312613, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.873737, Final residual = 0.05126, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.800713, Final residual = 0.0290762, No Iterations 2 GAMG: Solving for p, Initial residual = 0.917227, Final residual = 0.00712336, No Iterations 5 time step continuity errors : sum local = 6.51386e-13, global = -6.0799e-15, cumulative = 2.10588e-07 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.58764, Final residual = 0.390951, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.510518, Final residual = 0.0398422, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.612049, Final residual = 0.0284529, No Iterations 1000 GAMG: Solving for p, Initial residual = 1.86722e-09, Final residual = 1.86722e-09, No Iterations 0 time step continuity errors : sum local = 8.88265e-09, global = -1.68866e-15, cumulative = 2.10588e-07 smoothSolver: Solving for epsilon, Initial residual = 0.203351, Final residual = 1.17043e-07, No Iterations 1 bounding epsilon, min: -2.63513e+29 max: 3.13976e+42 average: 2.11151e+37 smoothSolver: Solving for k, Initial residual = 0.996435, Final residual = 0.967348, No Iterations 1000 ExecutionTime = 2438.95 s ClockTime = 2456 s Courant Number mean: 3.76917e-08 max: 706.294 deltaT = 2.42766e-23 Time = 1.1765e-05 AMI: Creating addressing and weights between 165774 source faces and 170168 target faces AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1 AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.927617, Final residual = 0.0460816, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.484562, Final residual = 0.028569, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.983595, Final residual = 0.0748445, No Iterations 3 GAMG: Solving for p, Initial residual = 0.796851, Final residual = 0.00460583, No Iterations 1 time step continuity errors : sum local = 1.00549e-06, global = -1.28245e-10, cumulative = 2.1046e-07 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.473033, Final residual = 0.458308, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.609904, Final residual = 0.380841, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.565501, Final residual = 0.00530204, No Iterations 1000 GAMG: Solving for p, Initial residual = 1.06618e-10, Final residual = 1.06618e-10, No Iterations 0 time step continuity errors : sum local = 1.57834e-05, global = -1.45424e-10, cumulative = 2.10314e-07 [14] #0 Foam::error::printStack(Foam::Ostream&)[9] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [14] #1 Foam::sigFpe::sigHandler(int) at ??:? [9] #1 Foam::sigFpe::sigHandler(int) at ??:? [9] #2 ? at ??:? [14] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [9] #3 Foam::symGaussSeidelSmoother::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) in "/lib/x86_64-linux-gnu/libc.so.6" [14] #3 Foam::symGaussSeidelSmoother::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) at ??:? [9] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [9] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [14] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [9] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? [14] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [14] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? [9] #7 at ??:? [14] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? [9] #8 at ??:? [14] #8 Foam::fvMatrix<double>::solve()Foam::fvMatrix<double>::solve() at ??:? [9] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:? [14] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:? [9] #10 Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:? [14] #10 Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:? [9] #11 ? at ??:? [14] #11 ? at ??:? [9] #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [9] #13 at ??:? [14] #12 __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6" [14] #13 at ??:? [93a46648863c:00174] *** Process received signal *** [93a46648863c:00174] Signal: Floating point exception (8) [93a46648863c:00174] Signal code: (-6) [93a46648863c:00174] Failing at address: 0x3e9000000ae [93a46648863c:00174] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f3c98ee14b0] [93a46648863c:00174] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f3c98ee1428] [93a46648863c:00174] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f3c98ee14b0] [93a46648863c:00174] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x744)[0x7f3c9a134634] [93a46648863c:00174] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x7f3c9a134888] [93a46648863c:00174] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4f6)[0x7f3c9a12b8f6] [93a46648863c:00174] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f3c9ceef95b] [93a46648863c:00174] [ 7] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x470f41] [93a46648863c:00174] [ 8] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xd4)[0x471194] [93a46648863c:00174] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solveIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x2f)[0x7f3c9e601f3f] [93a46648863c:00174] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x9a9)[0x7f3c9e66e099] [93a46648863c:00174] [11] pimpleDyMFoam[0x425849] [93a46648863c:00174] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f3c98ecc830] [93a46648863c:00174] [13] pimpleDyMFoam[0x4286c9] [93a46648863c:00174] *** End of error message *** ? at ??:? [93a46648863c:00179] *** Process received signal *** [93a46648863c:00179] Signal: Floating point exception (8) [93a46648863c:00179] Signal code: (-6) [93a46648863c:00179] Failing at address: 0x3e9000000b3 [93a46648863c:00179] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7ff09b5174b0] [93a46648863c:00179] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7ff09b517428] [93a46648863c:00179] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7ff09b5174b0] [93a46648863c:00179] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x744)[0x7ff09c76a634] [93a46648863c:00179] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x7ff09c76a888] [93a46648863c:00179] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4f6)[0x7ff09c7618f6] [93a46648863c:00179] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7ff09f52595b] [93a46648863c:00179] [ 7] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x470f41] [93a46648863c:00179] [ 8] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xd4)[0x471194] [93a46648863c:00179] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solveIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x2f)[0x7ff0a0c37f3f] [93a46648863c:00179] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x9a9)[0x7ff0a0ca4099] [93a46648863c:00179] [11] pimpleDyMFoam[0x425849] [93a46648863c:00179] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7ff09b502830] [93a46648863c:00179] [13] pimpleDyMFoam[0x4286c9] [93a46648863c:00179] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 14 with PID 179 on node 93a46648863c exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- OpenFOAM-4.1(31) Group_1 is the outlet and Wall_noSlip is the wall of the outer domain (the name is misleading. It's actually a slip wall condition). I'm aware that this names are confusing and I'll fix it now. This is the link to the project files (I deleted /constant/faces to reduce it's size). https://drive.google.com/drive/folde...bE?usp=sharing Any help would be appreciated! |
|
Tags |
ami, boundaries condition, propeller |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to set mesh conditions for propeller simulation? | cfd2020 | FLUENT | 2 | May 11, 2020 04:48 |
Quadrotor helicopter propeller simulation | Attesz | CFX | 14 | December 4, 2017 16:45 |
3D CFD simulation of propeller blade - Turbulent viscosity limited to viscosity ratio | Quentin_C | FLUENT | 1 | April 4, 2017 06:32 |
Problems about propeller simulation | duguyoudi | CFX | 13 | July 3, 2014 01:50 |
Run A propeller case with AMI and not accelerated flow behind the blades... | reza1980 | OpenFOAM | 6 | September 11, 2013 04:09 |