finding the potential v(x,y)

1.3k Views Asked by At

Consider the system

$\dot{x}=3x^2-1-e^{2y}, \dot{y}=-2xe^{2y}$

1)Show that $\frac {\partial{f}}{\partial{y}}=\frac {\partial{g}}{\partial{x}}$

2)Find the potential $V(x,y)$

3)Show trajectories always cross the equipotentials at right angles.

$\dot{x}=f(x,y)$ and $\dot{y}=g(x,y)$

I showed the first partial derivative part but i couldn't second and third part

2

There are 2 best solutions below

0
On

Given:

$$x'= f(x,y) = 3x^2-1-e^{2y}, y'= g(x,y) = -2xe^{2y}$$

1) Show that $\dfrac {\partial{f}}{\partial{y}}=\dfrac {\partial{g}}{\partial{x}}$

  • $\dfrac{\partial f}{\partial y} = -2 e^{2y}$
  • $\dfrac{\partial g}{\partial x} = -2 e^{2y}$

2) Find the potential function $V(x,y)$

  • $x' = 3x^2 - 1 - e^{2y} = -V_x$, so $$ V = - \int (3x^2 - 1 - e^{2y})~ dx + F(y) = -x^3 + x + xe^{2y} + F(y)$$
  • $y' = -2xe^{2y} = -V_y$, so $$V = - \int (-2xe^{2y})~dy + G(x) = xe^{2y} + G(x)$$
  • From the needed conditions for a gradient system, we can take $F(y) = 0$ and $G(x) = -x^3 + x$, hence

$$V(x,y) = -x^3 + x + xe^{2y}$$

3) Show trajectories always cross the equipotentials at right angles.

Please add what you have done here. If they only want you to show it, the following phase portrait may suffice. However, if they want a proof, that is harder.

enter image description here

5
On

Nota Bene: In preparing this answer, I in my haste used $\phi$ for the potential instead of the OP's $V(x, y)$; I have let it so stand to avoid the editing a change would entail. Ah well, as Treebeard the Ent might say, "hrooooom, hrooooom, let us not be hasty!" Sorry for any confusion, and remember, here

$\phi = V. \tag{0}$

This being said:

According to our OP hello, we can take as given the evident fact that

$\dfrac{\partial f}{\partial y} = \dfrac{\partial g}{\partial x}. \tag{1}$

Denoting the vector field $(f, g)^T$ by $X$:

$X = \begin{pmatrix} f(x, y) \\ g(x, y) \end{pmatrix}, \tag{2}$

we see that (1) implies that

$\nabla \times X = 0, \tag{3}$

which in turn implies the existence of a function $\phi(x, y)$ such that

$X = \nabla \phi. \tag{4}$

(3) also implies, via Stoke's theorem, that around any (sufficiently smooth) closed path $\Gamma$, we have

$\int_\Gamma X \cdot ds = \int_\Omega (\nabla \times X) \cdot dA = 0, \tag{5}$

where $\Omega$ is the region bounded by $\Gamma$ and $dA$ is the corresponding area element. (5) in turn shows that the function $\phi(x, y)$ may be defined, up to an arbitrary constant, by

$\phi(x, y) = \int_\gamma X \cdot ds \tag{6}$

for any path $\gamma(s)$ joining some fixed point $(x_0, y_0)$ and $(x, y)$, since (5) in fact implies such integrals are independent of the particular path chosen. We are thus free to evaluate $\phi$ by using a path of exceptional convenience. Since $X$ is defined on the entire $x$-$y$ plane $\Bbb R^2$, we may choose $(x_0, y_0) = (0, 0)$ and take $\gamma$ to run from $(0, 0)$ to $(x, 0)$ along the $x$-axis, then from $(x, 0)$ to $(x, y)$ along the segment parallel to the $y$-axis. We thus have

$\phi(x, y) = \int_{(0, 0)}^{(x,0)} (3s^2 - 2) ds + \int_{(x,0)}^{(x, y)} (-2xe^{2s}) ds; \tag{7}$

in the first of these integrals the $y$-component of $X$, $g(x, y)$, does not appear since the path is solely in the $x$-direction; $X \cdot \gamma'(s) = f(x, y)$ along this path component. And $f(x, y)$ does not appear in the second integral for the same reason, with the roles of $x$ and $y$ reversed. Evaluating these integrals, we find that

$\phi(x, y) = (x^3 - 2x) - x(e^{2y} - 1) = x^3 - x - xe^{2y}; \tag{8}$

it is now easy to see that

$\nabla \phi = (f, g)^T = X, \tag{9}$

that is, $\phi$ is a potential function for $X$. To see that the trajectories of $X$ are normal to the equipotentials of $\phi$, it is perhaps simplest at this point to observe that since the equipotentials are the curves

$\phi(x, y) = c \tag{10}$

for constant $c$, implicit differentiation of (10) with respect to $x$ yields

$f + gy'(x) = 0, \tag{11}$

which may in fact be written as

$X \cdot (1, y'(x)) = 0, \tag{12}$

and then use the fact that $(1, y'(x))$ is the tangent to the curve $(x, y(x))$, a parametric representation of the equipotential with parameter $x$ itself. This is valid everywhere except where $g(x, y) = 0$, that is, where $x \ne 0$; but it is easy to see that when $g = 0$ we are on the $y$-axis, and there $X$ points in the $x$-direction, so orthogonality of $X$ to the equipotentials holds there as well.

Hope this helps. Happy New Year,

and as always,

Fiat Lux!!!