“…For example, the nonlinear term u 2 in f u,x is linearised as U u, where U is the approximation of the interface velocity. The details regarding the iterative computation of the interface velocities using local BVPs are given in [4,5]. Thus, for the approximation of the flux f u,x i+1,j we solve the linearised local BVP…”