precice / precice/tutorials

Perpendicular Flap crashes when using subcycling

Open
#575 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

student
Dominant language
C
Stars
142
Forks
167
Avg merge
20h 24m
Merged PRs (30d)
11

Description

Used library versions:

  • Adapter version/commit number/last download: v1.3.1
  • OpenFOAM version: 2406
  • preCICE version: 3.1.2
Error message of last time step: 
PIMPLE: iteration 1
smoothSolver:  Solving for cellDisplacementx, Initial residual = 0.820797, Final residual = 7.13498e-19, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0.799523, Final residual = 2.76883e-17, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.000672219, No Iterations 68
DICPCG:  Solving for pcorr, Initial residual = 0.0618846, Final residual = 8.14925e-09, No Iterations 88
time step continuity errors : sum local = 1.37376e-09, global = 3.56395e-11, cumulative = -8.73401e-07
smoothSolver:  Solving for Ux, Initial residual = 0.120197, Final residual = 8.62079e-07, No Iterations 39
smoothSolver:  Solving for Uy, Initial residual = 0.112921, Final residual = 9.31561e-07, No Iterations 34
DICPCG:  Solving for p, Initial residual = 0.302693, Final residual = 0.000293536, No Iterations 66
DICPCG:  Solving for p, Initial residual = 0.0386374, Final residual = 3.68061e-05, No Iterations 55
time step continuity errors : sum local = 7.24018e-06, global = -7.07944e-07, cumulative = -1.58135e-06
DICPCG:  Solving for p, Initial residual = 0.0220954, Final residual = 1.83289e-05, No Iterations 62
DICPCG:  Solving for p, Initial residual = 0.00480982, Final residual = 4.01741e-06, No Iterations 56
time step continuity errors : sum local = 8.28686e-07, global = -2.70158e-08, cumulative = -1.60836e-06
DICPCG:  Solving for p, Initial residual = 0.0142835, Final residual = 1.2279e-05, No Iterations 60
DICPCG:  Solving for p, Initial residual = 0.00274941, Final residual = 2.44604e-06, No Iterations 60
time step continuity errors : sum local = 5.1083e-07, global = 1.47916e-08, cumulative = -1.59357e-06
DICPCG:  Solving for p, Initial residual = 0.0118773, Final residual = 1.1042e-05, No Iterations 59
DICPCG:  Solving for p, Initial residual = 0.00217576, Final residual = 8.73635e-09, No Iterations 77
time step continuity errors : sum local = 1.83504e-09, global = 7.94249e-11, cumulative = -1.59349e-06
ExecutionTime = 3.92 s  ClockTime = 7 s

---[precice]  Mapping "Force" for t=0.04 from "Fluid-Mesh" to "Solid-Mesh"
---[precice]  Time window completed
---[precice]  Mapping "Displacement" for t=0.04 from "Solid-Mesh" to "Fluid-Mesh"
---[precice]  iteration: 1 of 50 (min 1), time-window: 5, time: 0.04 of 5, time-window-size: 0.01, max-time-step-size: 0.01, ongoing: yes, time-window-complete: yes, write-iteration-checkpoint 
---[preciceAdapter] The coupling timestep completed. Writing the updated results.
Courant Number mean: 1.10395 max: 7.67523
Time = 0.05

PIMPLE: iteration 1
smoothSolver:  Solving for cellDisplacementx, Initial residual = 0.000513781, Final residual = 5.88609e-18, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0.000501744, Final residual = 1.19195e-17, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.000924515, No Iterations 67
DICPCG:  Solving for pcorr, Initial residual = 0.061722, Final residual = 7.62549e-09, No Iterations 88
time step continuity errors : sum local = 1.4287e-09, global = -2.19138e-11, cumulative = -1.59351e-06
smoothSolver:  Solving for Ux, Initial residual = 0.342152, Final residual = 7.97886e-07, No Iterations 31
smoothSolver:  Solving for Uy, Initial residual = 0.484288, Final residual = 9.72157e-07, No Iterations 32
DICPCG:  Solving for p, Initial residual = 0.140689, Final residual = 0.000118379, No Iterations 66
DICPCG:  Solving for p, Initial residual = 0.056835, Final residual = 4.63498e-05, No Iterations 61
time step continuity errors : sum local = 1.29794e-05, global = -4.66981e-07, cumulative = -2.06049e-06
DICPCG:  Solving for p, Initial residual = 0.0313542, Final residual = 2.32046e-05, No Iterations 62
DICPCG:  Solving for p, Initial residual = 0.00546486, Final residual = 3.78607e-06, No Iterations 61
time step continuity errors : sum local = 1.04716e-06, global = -5.96731e-09, cumulative = -2.06646e-06
DICPCG:  Solving for p, Initial residual = 0.0161969, Final residual = 1.44598e-05, No Iterations 62
DICPCG:  Solving for p, Initial residual = 0.00295532, Final residual = 2.51032e-06, No Iterations 62
time step continuity errors : sum local = 6.90473e-07, global = 7.53007e-09, cumulative = -2.05893e-06
DICPCG:  Solving for p, Initial residual = 0.0135775, Final residual = 1.24775e-05, No Iterations 61
DICPCG:  Solving for p, Initial residual = 0.00229308, Final residual = 6.60812e-09, No Iterations 76
time step continuity errors : sum local = 1.80962e-09, global = -7.22364e-12, cumulative = -2.05894e-06
ExecutionTime = 4.13 s  ClockTime = 7 s

---[precice]  Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice]  Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice]  iteration: 2 of 50 (min 1), time-window: 5, time: 0.04 of 5, time-window-size: 0.01, max-time-step-size: 0.01, ongoing: yes, time-window-complete: no, read-iteration-checkpoint 
Courant Number mean: 1.10395 max: 7.67523
Time = 0.05

PIMPLE: iteration 1
smoothSolver:  Solving for cellDisplacementx, Initial residual = 0.826971, Final residual = 8.66142e-18, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0.81761, Final residual = 2.76045e-17, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.000860607, No Iterations 73
DICPCG:  Solving for pcorr, Initial residual = 0.40583, Final residual = 9.06378e-09, No Iterations 97
time step continuity errors : sum local = 2.75936e-08, global = -3.97749e-11, cumulative = -2.05898e-06
smoothSolver:  Solving for Ux, Initial residual = 0.903074, Final residual = 789.26, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.732706, Final residual = 2151.95, No Iterations 1000
DICPCG:  Solving for p, Initial residual = 0.999632, Final residual = 0.000697483, No Iterations 88
DICPCG:  Solving for p, Initial residual = 0.872759, Final residual = 0.000826547, No Iterations 82
time step continuity errors : sum local = 48.2994, global = -0.752954, cumulative = -0.752956
DICPCG:  Solving for p, Initial residual = 0.336066, Final residual = 0.000239103, No Iterations 83
DICPCG:  Solving for p, Initial residual = 0.964627, Final residual = 0.000917575, No Iterations 59
time step continuity errors : sum local = 5066.05, global = 3918.51, cumulative = 3917.76
DICPCG:  Solving for p, Initial residual = 0.478094, Final residual = 0.000424423, No Iterations 51
DICPCG:  Solving for p, Initial residual = 0.953852, Final residual = 0.00093829, No Iterations 60
time step continuity errors : sum local = 105551, global = 193420, cumulative = 197338
DICPCG:  Solving for p, Initial residual = 0.48052, Final residual = 0.000391051, No Iterations 51
DICPCG:  Solving for p, Initial residual = 0.952827, Final residual = 5.26409e-09, No Iterations 115
time step continuity errors : sum local = 46.4586, global = 2.09115, cumulative = 197340
ExecutionTime = 4.51 s  ClockTime = 7 s

---[precice]  Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice]  Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice]  iteration: 3 of 50 (min 1), time-window: 5, time: 0.04 of 5, time-window-size: 0.01, max-time-step-size: 0.01, ongoing: yes, time-window-complete: no, read-iteration-checkpoint 
Courant Number mean: 1.10395 max: 7.67523
Time = 0.05

PIMPLE: iteration 1
smoothSolver:  Solving for cellDisplacementx, Initial residual = 0.867691, Final residual = 6.23348e-18, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0.861202, Final residual = 4.80158e-18, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.000977818, No Iterations 71
DICPCG:  Solving for pcorr, Initial residual = 0.455588, Final residual = 6.71708e-09, No Iterations 95
time step continuity errors : sum local = 1.05801e-08, global = -1.17836e-10, cumulative = 197340
smoothSolver:  Solving for Ux, Initial residual = 0.920002, Final residual = 1.95588e+20, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.83389, Final residual = 1.85741e+20, No Iterations 1000
DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.000400409, No Iterations 82
DICPCG:  Solving for p, Initial residual = 0.745248, Final residual = 0.000709922, No Iterations 37
time step continuity errors : sum local = -4.03737e+18, global = -1.93588e+18, cumulative = -1.93588e+18
DICPCG:  Solving for p, Initial residual = 0.357979, Final residual = 0.000326437, No Iterations 69
DICPCG:  Solving for p, Initial residual = 0.961185, Final residual = 0.000784698, No Iterations 22
time step continuity errors : sum local = -1.19919e+20, global = -1.88691e+19, cumulative = -2.0805e+19
DICPCG:  Solving for p, Initial residual = 0.687401, Final residual = 0.000624804, No Iterations 77
DICPCG:  Solving for p, Initial residual = 0.986636, Final residual = 0.000771809, No Iterations 13
time step continuity errors : sum local = 7.17806e+22, global = -4.7452e+19, cumulative = -6.82569e+19
DICPCG:  Solving for p, Initial residual = 0.50931, Final residual = 0.000141059, No Iterations 80
DICPCG:  Solving for p, Initial residual = 0.983318, Final residual = 8.43859e-09, No Iterations 99
time step continuity errors : sum local = 2.81033e+20, global = -2.58029e+18, cumulative = -7.08372e+19
ExecutionTime = 4.87 s  ClockTime = 8 s

---[precice]  Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice]  Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice]  iteration: 4 of 50 (min 1), time-window: 5, time: 0.04 of 5, time-window-size: 0.01, max-time-step-size: 0.01, ongoing: yes, time-window-complete: no, read-iteration-checkpoint 
Courant Number mean: 1.10395 max: 7.67523
Time = 0.05

PIMPLE: iteration 1
smoothSolver:  Solving for cellDisplacementx, Initial residual = 1, Final residual = 5.62587e-18, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 1, Final residual = 3.74224e-17, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.00089826, No Iterations 64
DICPCG:  Solving for pcorr, Initial residual = 0.996469, Final residual = 6.35919e-09, No Iterations 89
time step continuity errors : sum local = 7.06987e-06, global = 1.28168e-07, cumulative = -7.08372e+19
[stack trace]
=============
#1  Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#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 /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#6  ? in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#7  Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#8  ? in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
#9  ? in /lib/x86_64-linux-gnu/libc.so.6
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in /usr/lib/openfoam/openfoam2406/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
=============
Floating point exception (core dumped)

To reproduce this issue apply the following changes in perpendicular-flap/solid-openfoam/system/controlDict:

- deltaT          0.01;
+ deltaT          0.001;

I also tested with the following additional modifications that lead to crashes:
In perpendicular-flap/fluid-openfoam/system/controlDict:

- deltaT          0.01;
+ deltaT          0.001;

Setting substeps=true in the precice-config.xml didn't change the behavior.

It seems like the OpenFOAM-OpenFOAM perpendicular flap example works if the participants' deltaT are at least the time-window-size.

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reproducing the Perpendicular Flap case with OpenFOAM 2406, preCICE 3.1.2, adapter v1.3.1, and subcycling enabled. Inspect the coupling iterations around time 0.05 and the OpenFOAM stack trace; done means the case completes without the floating-point exception.

Written by the indexing model from the issue text.

Assessment

Domain
backend
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
25/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.