flow 1.3.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
Loading...
Searching...
No Matches
anderson_rev1_probes Namespace Reference

Functions

 build_case (flow, a, settings="production")
 
 inner (v, G=2)
 
 cmd_slowmode (flow, a)
 
 cmd_phase_a (flow, a)
 
 cmd_certify (flow, a)
 
 main ()
 

Detailed Description

Probes behind revision 1 of doc/steady_acceleration.md (an INSTRUMENT, not a gate).

Three subcommands, each on the oracle's cases (tests/study/anderson_oracle.py builds the solvers):

  slowmode   plain march on a case, per-step residual split (velocity / c_P-weighted pressure) and
             the full residual field at chosen steps; then the structure of the late residual:
             decay rates, alignment, energy shares, the fluid components of the face-aperture graph
             (sealed pockets) and the share of the pressure residual that sits off the main
             component at aperture thresholds 0 ... 0.2.                       (note §1.2)
  phase-a    phase A run unconditionally (acc.step(True) every call, no handover) with the monitor
             read after every map output: steps until the monitor stays within 1e-5 / 3e-5 of the
             reference and until the velocity / W residual first reaches 3e-7.  ("Revision 1" R1)
  certify    the full revision-1 driver with the monitor recorded after every step, then the
             certification blocks replayed: d, R and the remainder bound per block.   (§7, R4)

Run from the flow tree with the OpenMP pool bounded, e.g.

    OMP_NUM_THREADS=8 OMP_PROC_BIND=false PYTHONPATH=$PWD/build_omp \\
        python tests/study/anderson_rev1_probes.py phase-a --case sphere --scheme collocated --N 16 \\
        --window 5 --metric V --steps 300 --out e4.npz

The dense bed needs `--case bed --arrangement <npz>` (centres, radii, L; see the log for how it is
generated from the A1 scripts).

Function Documentation

◆ build_case()

anderson_rev1_probes.build_case (   flow,
  a,
  settings = "production" 
)
The oracle's sphere / bed cases (anderson_oracle.main's geometry), returning (solver,
oracle kwargs, settings, K-of-monitor).

Definition at line 37 of file anderson_rev1_probes.py.

Referenced by cmd_certify(), cmd_phase_a(), and cmd_slowmode().

◆ inner()

anderson_rev1_probes.inner (   v,
  G = 2 
)

Definition at line 71 of file anderson_rev1_probes.py.

Referenced by cmd_slowmode().

◆ cmd_slowmode()

anderson_rev1_probes.cmd_slowmode (   flow,
  a 
)

Definition at line 75 of file anderson_rev1_probes.py.

References build_case(), inner(), and mean().

◆ cmd_phase_a()

anderson_rev1_probes.cmd_phase_a (   flow,
  a 
)

Definition at line 165 of file anderson_rev1_probes.py.

References build_case(), and mean().

◆ cmd_certify()

anderson_rev1_probes.cmd_certify (   flow,
  a 
)

Definition at line 200 of file anderson_rev1_probes.py.

References build_case(), and mean().

◆ main()

anderson_rev1_probes.main ( )

Definition at line 238 of file anderson_rev1_probes.py.

References main().

Referenced by main().