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

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
 

Detailed Description

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)

Function Documentation

◆ xref()

robin.xref (   Bi)

Definition at line 13 of file robin.py.

◆ run()

robin.run (   ND,
  off,
  Bi 
)

Definition at line 17 of file robin.py.

Variable Documentation

◆ rng

robin.rng = rng.random((3, 2))

Definition at line 43 of file robin.py.

◆ NDs

list robin.NDs = [16, 32, 64, 128]

Definition at line 44 of file robin.py.

◆ ref

int robin.ref = xref(Bi) ** 2

Definition at line 46 of file robin.py.

◆ res

dict robin.res = {}

Definition at line 47 of file robin.py.

◆ line

str robin.line = None

Definition at line 53 of file robin.py.

◆ e

robin.e = np.sqrt(np.mean(np.square(res[(t, ND)])))

Definition at line 55 of file robin.py.

◆ prev

robin.prev = e

Definition at line 57 of file robin.py.