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

Functions

 sdf (x, y, z)
 
 lap (x, y, z)
 
 face_open (N, axis, h)
 
 run (N)
 

Variables

float R = 0.1
 
 prev = None
 
 rms
 
 mx
 
 nc
 
 o = np.log2(prev[1]/mx) if prev else float('nan')
 

Detailed Description

A-priori truncation test of the FV viscous operator L_visc(u) = mu[ Σ_f o_f(u_i-u_nbr) - Σ_a W_a g_a ]
vs the exact ∫_fluidCV(-mu ∇²u) dV, for the manufactured field u = sdf (u=0 on the sphere, ∇u=n̂,
∇²u = 2/r). The wall gradient g_a = n_a is EXACT for u=sdf, so this isolates the FACE-flux term.
Compares (A) face-CENTRE two-point flux vs (B) face flux evaluated at the open-area centroid.

Function Documentation

◆ sdf()

fv_operator_truncation.sdf (   x,
  y,
  z 
)

Definition at line 7 of file fv_operator_truncation.py.

Referenced by run().

◆ lap()

fv_operator_truncation.lap (   x,
  y,
  z 
)

Definition at line 8 of file fv_operator_truncation.py.

Referenced by run().

◆ face_open()

fv_operator_truncation.face_open (   N,
  axis,
  h 
)

Definition at line 9 of file fv_operator_truncation.py.

Referenced by run().

◆ run()

fv_operator_truncation.run (   N)

Definition at line 20 of file fv_operator_truncation.py.

References face_open(), lap(), and sdf().

Variable Documentation

◆ R

float fv_operator_truncation.R = 0.1

Definition at line 6 of file fv_operator_truncation.py.

◆ prev

tuple fv_operator_truncation.prev = None

Definition at line 49 of file fv_operator_truncation.py.

◆ rms

fv_operator_truncation.rms

Definition at line 52 of file fv_operator_truncation.py.

◆ mx

fv_operator_truncation.mx

Definition at line 52 of file fv_operator_truncation.py.

◆ nc

fv_operator_truncation.nc

Definition at line 52 of file fv_operator_truncation.py.

◆ o

fv_operator_truncation.o = np.log2(prev[1]/mx) if prev else float('nan')

Definition at line 53 of file fv_operator_truncation.py.