... (as in Fig. 2b ). The system is then solved using the ILU(0) preconditioned GMRES method (see Ref. 23 ) with adaptive subspace selection. For finite element discretization, we use the streamline upwind Petrov–Galerkin formulation proposed by Brooks and Hughes ( 20 ). Using linear triangular elements...