In some codes upwinding is available to decrease oscillations,

this accomplished through use of non-symmetric weighting functions on the concentration at each node

with:

where:

Wi	=	Ni	+		F(x)
Weight		linear basis function		modeler selects degree of upwinding ranging from = 0 to 1	a parabola or other function that fits the bill