My question goes along a similar line as this other one that was left unanswered because the OP did not have enough equations to produce a solution.
I want to solve for the vector field $\mathbf f(\mathbf r)$, defined over $V\subset\mathbb R^3$, that satisfies
$$\boldsymbol\nabla\cdot \mathbf f-\mathbf g(\mathbf r)\cdot\mathbf f=h(\mathbf r), \text{ and}$$ $$\boldsymbol\nabla\times \mathbf f=\mathbf 0.$$
Also, we assume that $\mathbf f|_{\partial V}$ is known.
How would one go about solving it? I saw this webpage that showed how to solve a similar kind of equations but with scalar field, and I am at a loss on how one would generalize that method to a vector field.
I've thought that one could possibly exploit the fact $\mathbf f$ is conservative to write the equations instead as
$$\nabla^2\Phi-\mathbf g\cdot\boldsymbol\nabla\Phi=h$$
but I don't see how that would help.
and $$\frac{\partial u}{\partial y}=\frac{\partial v}{\partial x}, \frac{\partial v}{\partial z}=\frac{\partial w}{\partial y}, \frac{\partial u}{\partial z}=\frac{\partial w}{\partial x}$$
The page only shows how to solve when there is 1 unknown $u$, but here we have 3 coupled ones.
– Chaotic May 02 '22 at 18:39