|
flow 1.3.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
|
Functions | |
| arg (name, default) | |
| col (path, ncol) | |
| main () | |
Variables | |
| float | MU = 3.33333333333333e-4 |
| float | RET = 127.3 |
| float | UTAU = RET * MU |
| float | LY = 2.0 |
Convert one TBFsolver output folder's statistics into the same NPZ `run_channel_18.py` writes.
python tbf_profiles.py <TBFsolver output folder> [--ny 160] [--out tbf_profiles.npz]
TBFsolver's statistics (`src/statistics/statistics.f90`) are TIME averages of PLANE integrals,
divided by the slice volume — so for the region `R` (`_W` whole, `_G` gas, `_L` liquid, with the
weight `cv` = 1, `c`, `1 − c` and `c` the GAS fraction, the same convention as our transcription):
stats_c <c>(y) 1 column
stats_u<R> <cv ux>, <cv uy>, <cv uz>, 9 columns
<cv uxux>, <cv uxuy>, <cv uxuz>,
<cv uyuy>, <cv uyuz>, <cv uzuz>
tauw the time-averaged wall shear stress
They are NOT conditioned averages: the liquid mean velocity is `<(1−c)u>/<1−c>`, which is what this
script forms, exactly as `run_channel_18.py` does with its own `L`-weighted sums. Both codes'
NPZ therefore carry the same quantities under the same names and the plotting script can overlay
them without further arithmetic.
| tbf_profiles.arg | ( | name, | |
| default | |||
| ) |
Definition at line 32 of file tbf_profiles.py.
Referenced by main().
| tbf_profiles.col | ( | path, | |
| ncol | |||
| ) |
Definition at line 36 of file tbf_profiles.py.
Referenced by main().
| tbf_profiles.main | ( | ) |
| float tbf_profiles.MU = 3.33333333333333e-4 |
Definition at line 26 of file tbf_profiles.py.
| float tbf_profiles.RET = 127.3 |
Definition at line 27 of file tbf_profiles.py.
Definition at line 28 of file tbf_profiles.py.
| float tbf_profiles.LY = 2.0 |
Definition at line 29 of file tbf_profiles.py.