OK, from the comments I understand that the problem boils down to

$$\frac{\partial}{\partial t}\Phi_2(t,x_1,x_2)=\theta\bigl(F[\Phi_2(t,x_1,x_2)]-x_1\bigr),\;\;\Phi_2(0,x_1,x_2)=x_2,$$

where $\theta$ is the unit step function. Let me assume that $F$ is a nondecreasing function. Then the solution is

$$\Phi_2(t,x_1,x_2)=\begin{cases}
t+x_2&\text{if}\;\;F(x_2)>x_1\\
x_2&\text{if}\;\;F(x_2)\leq x_1.
\end{cases}
$$

A more general choice of $F$ can be readily accommodated, by piecing together increasing and decreasing segments. For any $F$, the function $\Phi_2(t,x_1,x_2)=x_2$ whenever $F(x_2)\leq x_1$, so we only need to consider regions in which $F(x_2)>x_1$ and $\Phi_2$ increases linearly with $t$ until $F(\Phi_2)$ becomes smaller than $x_1$.

as requested, the Mathematica code

```
Manipulate[Module[{sol =
NDSolve[{phi2'[t] == UnitStep[Sin[phi2[t]] - x1], phi2[0] == x2},
phi2, {t, 0, tfinal}]},
Plot[Evaluate[phi2[t] /. sol], {t, 0, tfinal}]],
{x1, 0, 2}, {x2, 0, 2}, {tfinal, 1, 10}]
```

will produce a plot of $\Phi_2(t)$ (for $F=\sin$) where you can vary $t_{\rm final}$, $x_1$, and $x_2$.

alternatively, the Mathematica code

```
sol = NDSolve[{D[phi2[t, x1, x2], t] ==
UnitStep[Sin[phi2[t, x1, x2]] - x1], phi2[0, x1, x2] == x2}, phi2,
{t, 0, 5}, {x1, 0, 2}, {x2, 0, 2}];
Manipulate[VectorPlot[{x1, Evaluate[phi2[t, x1, x2] /. sol]},
{x1, 0, 2}, {x2, 0, 2}], {t, 0, 5}]
Manipulate[StreamPlot[{x1, Evaluate[phi2[t, x1, x2] /. sol]},
{x1, 0, 2}, {x2, 0, 2}], {t, 0, 5}]
```

will produce a vector plot or a stream plot of $(\Phi_1,\Phi_2)$ in the $x$-$y$ plane, with $t$ as a parameter that you can vary.

2more comments