Skip to content

upwind_1d(x, F)

First-order accurate, unconditionally stable, non-conservative upwind advection scheme.

Source code in openterrace/convection_schemes/upwind_1d.py 
 4
 5
 6
 7
 8
 9
10
11
12
13
14
@nb.njit
def upwind_1d(x, F):
    """First-order accurate, unconditionally stable, non-conservative upwind advection scheme.
    """
    _out = np.zeros_like(x)
    for j in range(0, x.shape[0]):
        for i in range(1, x.shape[1]-1):
            _out[j,i] = x[j,i+1]*(-np.minimum(F[0,j,i],0))\
                + x[j,i-1]*(np.maximum(F[1,j,i],0))\
                + x[j,i]*(np.minimum(F[0,j,i],0)-np.maximum(F[1,j,i],0))
    return _out