Right answer is (c) Picard iteration
To explain I would say: The usual solution of non-linear coupled system is the sequential decoupled method. For this, the non-linear terms are linearized using the Picard iteration method. This is applied for convection and source terms of the equation.