Suppose I had a boundary value problem:
$$\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega$$ $$\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega$$ $$u=h \text{ in } \partial\Omega$$
My goal is to decompose the solution of this coupled problem into a sequence of uncoupled PDE's. To decouple the system, I'm applying a fixed point iteration over a sequence of approximations $(u^k,v^k)$ such that
$$\frac{d^2u^k}{dx^2} + \frac{dv^{k-1}}{dx}=f $$ $$\frac{du^{k-1}}{dx} +\frac{d^2v^k}{dx^2} =g $$
Theoretically, this would allow me to solve both equations as a purely elliptic PDE's. However, I have never seen fixed point iterations applied to PDE's in this way. I've seen fixed point iterations applied to the numerically discretized equations (finite difference method, finite element method, etc.), but never to the continuous equations directly.
Am I violating any blatant mathematical principle by doing this? Is this mathematically valid? Could I solve the coupled PDE as a sequence of uncoupled PDE's by using fixed point iteration applied to the CONTINUOUS variable problem, rather than the DISCRETE variable problem?
At this point, I'm not really concerned with whether it is practical to use this method, but rather whether it is theoretically plausible. Any feedback would be greatly appreciated!