Perpendicular Flap crashes when using subcycling
Aperta
I maintainer di solito rispondono entro 1 giorno
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 4/5
- Tempo stimato
- 3-5 giorni
- Idoneità per principianti
- 25/100
- Tipo di issue
- Bug
- Chiarezza
- Da chiarire
- Stato di attività
- Ferma
- Ambito
- backend
Direzione di ricerca
Inizia riproducendo il caso Perpendicular Flap con OpenFOAM 2406, preCICE 3.1.2, l’adapter v1.3.1 e il subcycling abilitato. Ispeziona le iterazioni di accoppiamento intorno al tempo 0.05 e lo stack trace di OpenFOAM; il lavoro è completato quando il caso termina senza l’eccezione in virgola mobile.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
student
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] [0m Mapping "Force" for t=0.04 from "Fluid-Mesh" to "Solid-Mesh"
---[precice] [0m Time window completed
---[precice] [0m Mapping "Displacement" for t=0.04 from "Solid-Mesh" to "Fluid-Mesh"
---[precice] [0m 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] [0m Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice] [0m Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice] [0m 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] [0m Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice] [0m Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice] [0m 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] [0m Mapping "Force" for t=0.05 from "Fluid-Mesh" to "Solid-Mesh"
---[precice] [0m Mapping "Displacement" for t=0.05 from "Solid-Mesh" to "Fluid-Mesh"
---[precice] [0m 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.
- Lingua principale
- C
- Stelle
- 143
- Fork
- 169
- Merge medio
- 1g 11h
- PR unite (30g)
- 20
Preparare l'ambiente
- Nessun Dockerfile né file Docker Compose
- Ha un modello di pull request
- Leggi la guida per i contributori
Come iniziare
- Leggi tutta la issue e poi la guida ai contributi del progetto.
- Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
- Fai un fork del repository e lavora su un branch.
- Apri una pull request che faccia riferimento al numero della issue.
Altre issue di precice/tutorials
-
good first issue
Difficoltà 2/5 1-2 giorni Idoneità per principianti 62/100
precice/tutorials#617 · 2 commenti · 1 reazione ·
I maintainer di solito rispondono entro 1 giorno
-
system-tests
Difficoltà 4/5 3-5 giorni Idoneità per principianti 55/100
I maintainer di solito rispondono entro 1 giorno
-
system-tests
Difficoltà 4/5 3-5 giorni Idoneità per principianti 45/100
precice/tutorials#933 · 1 commento ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 3/5 1-2 giorni Idoneità per principianti 66/100
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 4/5 3-5 giorni Idoneità per principianti 42/100
precice/tutorials#877 · 1 commento ·
I maintainer di solito rispondono entro 1 giorno
Tutte le issue di precice/tutorials
Issue simili
-
initramfs: -type f (#18686) skips the libcurl.so.4 symlink, libcurl no longer copied into initramfsAperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100
I maintainer di solito rispondono entro 2 giorni
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100
UNIVERSE-HPC/course-material#283 ·
-
bug C/C++ code
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 78/100
webarkit/WebARKitLib#85 ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 75/100
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 62/100
I maintainer di solito rispondono entro 1 giorno