where N is the number of elements, and ϕ ϕ i ( x , y ) are the shape functions associated with each element. Substituting this expression into the Navier–Stokes equation and applying the finite element method, we can derive a system of ordinary differential equations: