Navier Stokes preconditioners

[luca-heltai]

• [x] Check that the current preconditioner for Stokes works with Linear Operators
• [x] Write a section in the Navier-Stokes interface, to select different preconditioners, and assemble the right one according to the parameter file (TODO: choose a list of preconditioners for the tests)
• [x] Make the grid and manifold descriptor for two benchmarks of Navier-Stokes:
  • [x] Lid Cavity problem
  • [x] Flow past a cylinder
• [x] Test all implemented preconditioners on the two problems above, checking scalability, for different Reynolds Numbers: 1, 10, 100, 1.000, 10.000 (not all preconditioners will work everywhere)
• [x] Repeat the tests above with a different selection of Finite Element system (Q(3)-Q(2)/Q(4)-Q(3))

[luca-heltai]

@ESeNonFossiIo , punto uno e' fatto? (se si: spunta..)

[ESeNonFossiIo]

Se basta che l'implementazione usata in step-32 funziona ok: assemblo i blocchi come linear operator e assemblo il precondizionatore usando questi blocchi. Senza usare back o forward.

[ESeNonFossiIo]

[ESeNonFossiIo]

$\nu=200$

[luca-heltai]

Figo!

[ESeNonFossiIo]

per $Flow past a cylinder$ avete una griglia già fatta da passarmi? ... o implementiamo il grid merge in grid generator?

[luca-heltai]

deve essere negli esempi di deal2lkit

[ESeNonFossiIo]

Da un test come si carica un file? non mi trova il path.

[ESeNonFossiIo]

con \nu = 10

[ESeNonFossiIo]

se metto \nu~0.1 non converge

[ESeNonFossiIo]

[ESeNonFossiIo]

Re ~ 400

[ESeNonFossiIo]

v=2*t per t<.5 v=1 per t>=.5 immagine presa a tempo t~1.2

[ESeNonFossiIo]

[ESeNonFossiIo]

t~2.2

[ESeNonFossiIo]

[ESeNonFossiIo]

t~3.2