|
flow 1.3.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
|
Functions | |
| xref (Bi) | |
| run (ND, off, Bi) | |
Variables | |
| rng = rng.random((3, 2)) | |
| list | NDs = [16, 32, 64, 128] |
| int | ref = xref(Bi) ** 2 |
| dict | res = {} |
| str | line = None |
| e = np.sqrt(np.mean(np.square(res[(t, ND)]))) | |
| prev = e | |
Robin eigenproblem on the disc: -lap c = mu c, D dc/dn + k c = 0 (n out of the fluid), D=R=1.
Exact: sqrt(mu) = x with x J1(x) = Bi J0(x). Variants:
RF : flux form - aperture FV, kappa mass, all kappa>0 cells, wall flux A_w c/(d_centroid + 1/k)
RFc : same, d measured from the cell CENTRE (sdf), clamped >= 0.05 h
RG : linear ghost on centre-fluid cells, c_Gamma eliminated through the normal Robin relation
(diag += k|n_a| / (h (1+k d_n)) per crossing link; k->0 = staircase Neumann)