|
[Sponsors] |
July 21, 2021, 15:06 |
Multiple AMI zones not working
|
#1 |
New Member
Pedro
Join Date: Jul 2021
Posts: 11
Rep Power: 5 |
Hello everyone!
I'm trying to study how the wake produced by a wind turbine affects to another wind turbine placed next to the first one (Check attached images for a visual). For that purpose, I set a 2D case with 2 vertical axys wind turbines with 1DoF each and two AMI zones, so I can get a flow induced movement. (I'm using this Holzmann tutorial as a base case: https://holzmann-cfd.com/community/t...l-wind-turbine). Also, I'm using OpenFOAM2012, as long as OpenFOAM8 doesn't support multiple 6DoF body motion. Everything works fine when I give 1DoF to only one turbine (doesn't matter which one) and the other stays as a wall, but when I try to give 1DoF to both of them, the solver crashes after a few iterations. After checking the logSolve file and the mesh in paraview, it seens that the second AMI zone destroys the mesh (see detailed images). Here is a extract of the logSolve file, where you can see how Patch target sum weights (Red line) it is not close to 1 anymore and in consecuence, p diverges and so time steps continuity error gets huge just before simulation crashes: Code:
PIMPLE: iteration 1 forces forces: rho: rhoInf Freestream density (rhoInf) set to 1.2 Not including porosity effects 6-DoF rigid body motion Centre of rotation: (0 0 0.025) Centre of mass: (0 0 0.025) Orientation: (1 -3.38655e-07 0 3.38655e-07 1 0 0 0 1) Linear velocity: (0 0 0) Angular velocity: (0 0 0.00257071) forces forces: rho: rhoInf Freestream density (rhoInf) set to 1.2 Not including porosity effects 6-DoF rigid body motion Centre of rotation: (2.5 0 0.025) Centre of mass: (2.5 0 0.025) Orientation: (1 -3.37314e-07 0 3.37314e-07 1 0 0 0 1) Linear velocity: (-0 -0 -0) Angular velocity: (0 0 0.00244335) AMI: Creating addressing and weights between 252 source faces and 256 target faces AMI: Patch source sum(weights) min:0.999188 max:1.00063 average:0.999958 AMI: Patch target sum(weights) min:0.971359 max:1.00632 average:0.999636 AMI: Creating addressing and weights between 252 source faces and 256 target faces AMI: Patch source sum(weights) min:0.981669 max:1.00147 average:0.999634 AMI: Patch target sum(weights) min:0.776676 max:1.31567 average:1.00391 FDICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.9477e-07, No Iterations 677 time step continuity errors : sum local = 4.12212e-12, global = 2.76544e-14, cumulative = 7.04299e-09 FDICPCG: Solving for p, Initial residual = 0.672125, Final residual = 2.1746, No Iterations 4000 time step continuity errors : sum local = 0.000157631, global = 3.22905e-07, cumulative = 3.29948e-07 FDICPCG: Solving for p, Initial residual = 0.783422, Final residual = 0.000685805, No Iterations 128 time step continuity errors : sum local = 4.11302e-07, global = 1.879e-07, cumulative = 5.17848e-07 FDICPCG: Solving for p, Initial residual = 0.985528, Final residual = 0.488492, No Iterations 4000 time step continuity errors : sum local = 0.0109226, global = -3.73584e-07, cumulative = 1.44265e-07 DILUPBiCGStab: Solving for omega, Initial residual = 0.10746, Final residual = 6.52578e-05, No Iterations 2 bounding omega, min: -11.2342 max: 349628 average: 176.709 DILUPBiCGStab: Solving for k, Initial residual = 0.0128361, Final residual = 1.64071e-07, No Iterations 3 bounding k, min: -3.31861e-05 max: 0.3456 average: 0.000188273 PIMPLE: iteration 2 FDICPCG: Solving for p, Initial residual = 0.997817, Final residual = 0.000982528, No Iterations 94 time step continuity errors : sum local = 0.00435951, global = -0.000702525, cumulative = -0.000702381 FDICPCG: Solving for p, Initial residual = 0.364384, Final residual = 0.000334422, No Iterations 72 time step continuity errors : sum local = 0.0318397, global = -0.000639627, cumulative = -0.00134201 FDICPCG: Solving for p, Initial residual = 0.157277, Final residual = 7.55113e-06, No Iterations 4000 time step continuity errors : sum local = 0.0143722, global = 0.00066907, cumulative = -0.000672938 DILUPBiCGStab: Solving for omega, Initial residual = 0.361791, Final residual = 0.000272681, No Iterations 3 DILUPBiCGStab: Solving for k, Initial residual = 0.352445, Final residual = 0.000127851, No Iterations 3 PIMPLE: iteration 3 FDICPCG: Solving for p, Initial residual = 0.48606, Final residual = 0.316725, No Iterations 4000 time step continuity errors : sum local = 2718.75, global = 7.53898, cumulative = 7.53831 FDICPCG: Solving for p, Initial residual = 0.00660017, Final residual = 0.000568018, No Iterations 4000 time step continuity errors : sum local = 2950.82, global = -22.8372, cumulative = -15.2989 FDICPCG: Solving for p, Initial residual = 0.000567995, Final residual = 0.000391576, No Iterations 4000 time step continuity errors : sum local = 167581, global = -3837.96, cumulative = -3853.26 DILUPBiCGStab: Solving for omega, Initial residual = 0.261805, Final residual = 0.000189136, No Iterations 2 DILUPBiCGStab: Solving for k, Initial residual = 0.289841, Final residual = 7.92111e-05, No Iterations 3 PIMPLE: iteration 4 Could someone please shed a light about this issue? I would be really grateful for any kind of suggestion, I've been weeks stuck on this. If you need any more info, please ask for it, I'm quite new to OpenFOAM and I don't know if I have explained myself correctly. Thank you for reading. Pedro. |
|
July 23, 2021, 03:31 |
|
#2 |
New Member
Pedro
Join Date: Jul 2021
Posts: 11
Rep Power: 5 |
I'm leaving here some of the important files, so you have more info.
dynamicMeshDict Code:
dynamicFvMesh dynamicMotionSolverListFvMesh; //dynamicMotionSolverFvMesh solvers { VAWT { motionSolverLibs ("libsixDoFRigidBodyMotion.so"); solver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { patches (VAWT); innerDistance 0.6; outerDistance 0.7; rho rhoInf; rhoInf 1.2; mass 0.691737; momentOfInertia (0.01035 0.00532 0.005317); centreOfMass (0 0 0.025); orientation ( 1 0 0 0 1 0 0 0 1 ); velocity (0 0 0); acceleration (0 0 0); angularMomentum (0 0 0); torque (0 0 0); report on; accelerationRelaxation 0.7; solver { type Newmark; } constraints { zAxis { sixDoFRigidBodyMotionConstraint axis; axis (0 0 1); } fixedPt { sixDoFRigidBodyMotionConstraint point; centreOfRotation (0 0 0.025); } } } } VAWT2 { motionSolverLibs ("libsixDoFRigidBodyMotion.so"); solver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { patches (VAWT2); innerDistance 0.6; outerDistance 0.7; rho rhoInf; rhoInf 1.2; mass 0.691737; momentOfInertia (0.01035 0.00532 0.005317); centreOfMass (2.5 0 0.025); orientation ( 1 0 0 0 1 0 0 0 1 ); velocity (0 0 0); acceleration (0 0 0); angularMomentum (0 0 0); torque (0 0 0); report on; accelerationRelaxation 0.7; solver { type Newmark; } constraints { zAxis { sixDoFRigidBodyMotionConstraint axis; axis (0 0 1); } fixedPt { sixDoFRigidBodyMotionConstraint point; centreOfRotation (2.5 0 0.025); } } } } } Code:
ddtSchemes { default Euler; } gradSchemes { default none; grad(U) Gauss linear; grad(pcorr) Gauss linear; grad(p) Gauss linear; grad(S) cellLimited Gauss linear 1; grad(omega) cellLimited Gauss linear 1; grad(k) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss linearUpwindV grad(U); div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,S) Gauss limitedLinear01 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; laplacian(diffusivity,cellMotionU) Gauss linear uncorrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } Code:
solvers { p { solver PCG; preconditioner FDIC; tolerance 1e-8; relTol 0.01; minIter 1; maxIter 400; } pFinal { $p; tolerance 1e-08; relTol 0; } pcorr { $p relTol 0.01; } pcorrFinal { $pcorr tolerance 1e-06; relTol 0; } "(k|omega|U|S)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-8; relTol 0.1; minIter 1; } "(k|omega|U|S)Final" { $U; tolerance 1e-08; relTol 0; minIter 1; } cellMotionUx { solver PCG; preconditioner FDIC; tolerance 1e-06; relTol 0; minIter 1; } } PIMPLE { correctPhi true; momentumPredictor false; turbOnFinalIterOnly false; consistent true; nNonOrthogonalCorrectors 0; nOuterCorrectors 15; nCorrectors 3; residualControl { p { $p; tolerance 1e-5; relTol 0.01; } U { $U; tolerance 4e-6; relTol 0.01; } } outerCorrectorResidualControl { p { relTol 0; tolerance 1e-3; } U { relTol 0; tolerance 1e-3; } } } relaxationFactors { fields { p 0.85; pFinal 1; } equations { "(k|S|omega|U)" 0.8; "(k|S|omega|U)Final" 1; } } Code:
application pimpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 30.5; deltaT 1e-05; writeControl adjustableRunTime; writeInterval 0.005; purgeWrite 0; writeFormat binary; writePrecision 6; writeCompression true; timeFormat general; timePrecision 6; runTimeModifiable true; adjustTimeStep yes; maxCo 0.8; maxDeltaT 0.001; libs ( "libsolverFunctionObjects.so" ); Error message Code:
[0] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [0] #1 Foam::sigFpe::sigHandler(int) at ??:? [0] #2 ? at ??:? in /lib/x86_64-linux-gnu/libpthread.so.0 [1] #1 Foam::sigFpe::sigHandler(int)[0] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:? [0] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #6 at ??:? [1] #2 ?? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #7 __libc_start_main in /lib/x86_64-linux-gnu/libpthread.so.0 [1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in /lib/x86_64-linux-gnu/libc.so.6 [0] #8 ? at ??:? [1] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #5 in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [pedro-HP-Laptop-15s-fq1xxx:09771] *** Process received signal *** [pedro-HP-Laptop-15s-fq1xxx:09771] Signal: Floating point exception (8) [pedro-HP-Laptop-15s-fq1xxx:09771] Signal code: (-6) [pedro-HP-Laptop-15s-fq1xxx:09771] Failing at address: 0x3e80000262b [pedro-HP-Laptop-15s-fq1xxx:09771] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f58dbccb3c0] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xcb)[0x7f58dbccb24b] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f58dbccb3c0] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 3] /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0xb7)[0x7f58dc897e17] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 4] pimpleFoam(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKNS_11dimensionedIdEERKS6_+0x40)[0x560fc9da4950] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 5] pimpleFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x294)[0x560fc9dd9bd4] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 6] pimpleFoam(+0x2f5ce)[0x560fc9d895ce] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f58dbaeb0b3] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 8] pimpleFoam(+0x3160e)[0x560fc9d8b60e] [pedro-HP-Laptop-15s-fq1xxx:09771] *** End of error message *** Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #6 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 [1] #8 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [pedro-HP-Laptop-15s-fq1xxx:09772] *** Process received signal *** [pedro-HP-Laptop-15s-fq1xxx:09772] Signal: Floating point exception (8) [pedro-HP-Laptop-15s-fq1xxx:09772] Signal code: (-6) [pedro-HP-Laptop-15s-fq1xxx:09772] Failing at address: 0x3e80000262c [pedro-HP-Laptop-15s-fq1xxx:09772] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f94b97b83c0] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xcb)[0x7f94b97b824b] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f94b97b83c0] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 3] /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0xb7)[0x7f94ba384e17] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 4] pimpleFoam(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKNS_11dimensionedIdEERKS6_+0x40)[0x556c7c64c950] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 5] pimpleFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x294)[0x556c7c681bd4] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 6] pimpleFoam(+0x2f5ce)[0x556c7c6315ce] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f94b95d80b3] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 8] pimpleFoam(+0x3160e)[0x556c7c63360e] [pedro-HP-Laptop-15s-fq1xxx:09772] *** End of error message *** -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun noticed that process rank 0 with PID 0 on node pedro-HP-Laptop-15s-fq1xxx exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Please, if you need more information or a futher explanation, ask for it. I'm totally desperate about this issue. Thank you. |
|
July 29, 2021, 03:21 |
|
#3 |
New Member
Norge
Join Date: Dec 2020
Posts: 7
Rep Power: 5 |
Hi,
Does the Courant number increase during the steps before the crash? If the courant number increases above 1, the simulation will crash. Would it be an idea to create a refinement area between the interacting wind mills? |
|
July 29, 2021, 05:26 |
|
#4 |
New Member
Pedro
Join Date: Jul 2021
Posts: 11
Rep Power: 5 |
Hello abeda, thanks for your reply.
Courant number is limited to 0.8, as you can see in my controlDict file. I have tried lowering Co to 0.3, but it was useless. Maybe I can try a finer mesh, but I won't be able to make a new mesh until next week. Best regards. |
|
July 29, 2021, 06:42 |
|
#5 |
New Member
Norge
Join Date: Dec 2020
Posts: 7
Rep Power: 5 |
Hi,
I didn't mean the courant number given in controlDict file. Sorry for being so unclear. I believe pimpleFoam also prints out the courant number between each time step. So, what is the courant numbers printed out by pimpleFoam right before the crash message given below? When you reduced the courant number in controlDict file to 0.3, did it crash sooner than when the courant number was 0.8? Error message: Code:
[0] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [0] #1 Foam::sigFpe::sigHandler(int) at ??:? [0] #2 ? at ??:? in /lib/x86_64-linux-gnu/libpthread.so.0 [1] #1 Foam::sigFpe::sigHandler(int)[0] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:? [0] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #6 at ??:? [1] #2 ?? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [0] #7 __libc_start_main in /lib/x86_64-linux-gnu/libpthread.so.0 [1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in /lib/x86_64-linux-gnu/libc.so.6 [0] #8 ? at ??:? [1] #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #5 in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [pedro-HP-Laptop-15s-fq1xxx:09771] *** Process received signal *** [pedro-HP-Laptop-15s-fq1xxx:09771] Signal: Floating point exception (8) [pedro-HP-Laptop-15s-fq1xxx:09771] Signal code: (-6) [pedro-HP-Laptop-15s-fq1xxx:09771] Failing at address: 0x3e80000262b [pedro-HP-Laptop-15s-fq1xxx:09771] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f58dbccb3c0] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xcb)[0x7f58dbccb24b] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f58dbccb3c0] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 3] /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0xb7)[0x7f58dc897e17] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 4] pimpleFoam(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKNS_11dimensionedIdEERKS6_+0x40)[0x560fc9da4950] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 5] pimpleFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x294)[0x560fc9dd9bd4] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 6] pimpleFoam(+0x2f5ce)[0x560fc9d895ce] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f58dbaeb0b3] [pedro-HP-Laptop-15s-fq1xxx:09771] [ 8] pimpleFoam(+0x3160e)[0x560fc9d8b60e] [pedro-HP-Laptop-15s-fq1xxx:09771] *** End of error message *** Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #6 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [1] #7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 [1] #8 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [pedro-HP-Laptop-15s-fq1xxx:09772] *** Process received signal *** [pedro-HP-Laptop-15s-fq1xxx:09772] Signal: Floating point exception (8) [pedro-HP-Laptop-15s-fq1xxx:09772] Signal code: (-6) [pedro-HP-Laptop-15s-fq1xxx:09772] Failing at address: 0x3e80000262c [pedro-HP-Laptop-15s-fq1xxx:09772] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f94b97b83c0] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xcb)[0x7f94b97b824b] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x153c0)[0x7f94b97b83c0] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 3] /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0xb7)[0x7f94ba384e17] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 4] pimpleFoam(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKNS_11dimensionedIdEERKS6_+0x40)[0x556c7c64c950] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 5] pimpleFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x294)[0x556c7c681bd4] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 6] pimpleFoam(+0x2f5ce)[0x556c7c6315ce] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f94b95d80b3] [pedro-HP-Laptop-15s-fq1xxx:09772] [ 8] pimpleFoam(+0x3160e)[0x556c7c63360e] [pedro-HP-Laptop-15s-fq1xxx:09772] *** End of error message *** -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun noticed that process rank 0 with PID 0 on node pedro-HP-Laptop-15s-fq1xxx exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- |
|
August 3, 2021, 04:44 |
|
#6 |
New Member
Pedro
Join Date: Jul 2021
Posts: 11
Rep Power: 5 |
Hello abeda, sorry for the late reply.
Here is what I get when I use maxCo = 0.8 (this is the last piece of info about Co, since in the next iteration the simulation crashes): Courant Number mean: 0.0197082 max: 0.831334 deltaT = 3.81029e-05 Time = 0.000922992 And here MaxCo = 0.3 : Courant Number mean: 0.00579354 max: 0.332099 deltaT = 1.05022e-05 Time = 0.000232002 Looks like in both cases the simulation crashed when maxCo is surpassed. Also, I have notice that both of them crashed at the same iteration (PIMPLE: iteration 4) even thought time step is smaller where maxCo = 0.3 Thank you for replying and sorry again for my late reply. Best regards. |
|
September 10, 2021, 07:03 |
|
#7 |
New Member
Pedro
Join Date: Jul 2021
Posts: 11
Rep Power: 5 |
Please, does anyone have any experience using multiple AMI zones? I am still stuck on this, any piece of advice would be really welcomed.
Thank you |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
pimpleDyMfoam , AMI -- The mesh has multiple regions not connected by any face | coolcrasher | OpenFOAM Running, Solving & CFD | 3 | May 16, 2019 05:58 |
Splitting Geometry into Multiple Zones | max.kozak | FLUENT | 1 | April 13, 2018 04:17 |
Replacing a mesh with multiple zones | Gert-Jan | FLUENT | 0 | March 23, 2018 16:26 |
Isosurface using multiple zones | FabianDe | Tecplot | 0 | July 18, 2016 08:43 |
Mixing Planes with Multiple Downstream Zones | Will Anderson | FLUENT | 0 | November 10, 2010 17:26 |