# Solve Linear PDE via a change of variables

Show 40 post(s) from this thread on one page
Page 1 of 2 12 Last
• Feb 11th 2011, 08:06 AM
bugatti79
Solve Linear PDE via a change of variables
$\displaystyle \omega(r,\theta)=u(r cos (\theta),r sin (\theta))$ for some function u(x,y)

$\displaystyle \displaystyle y \frac{\partial u}{\partial x}-x\frac{\partial u}{\partial y}=1$

Rewrite the PDE and IC's in terms of $\displaystyle \omega,r, \theta$ by expressing

$\displaystyle \displaystyle \frac{\partial \omega}{\partial r}$ and $\displaystyle \displaystyle \frac{\partial \omega}{\partial \theta}$in terms of $\displaystyle \displaystyle \frac{\partial u}{\partial x}$ and $\displaystyle \displaystyle \frac{\partial u}{\partial y}$ and solve.

My attempt: I believe the first equation is of the form

$\displaystyle \displaystyle \omega(r,\theta) = u(x(r,\theta),y(r,\theta))$ therefore$\displaystyle \begin{cases} \frac{\partial \omega}{\partial r} = \frac{\partial \omega}{\partial x} \frac{\partial x}{\partial r} + \frac{\partial \omega}{\partial y} \frac{\partial y}{\partial r}\\\frac{\partial \omega}{\partial \theta} = \frac{\partial \omega}{\partial x} \frac{\partial x}{\partial \theta} + \frac{\partial \omega}{\partial y} \frac{\partial y}{\partial \heta}\end{cases}$ but they are not in terms of $\displaystyle \displaystyle \frac{\partial u}{\partial x}$ and $\displaystyle \displaystyle \frac{\partial u}{\partial y}$

Have I gone one step too far?
• Feb 11th 2011, 08:11 AM
Ackbeet
Could be way off here, but couldn't you just use this:

$\displaystyle \displaystyle\begin{cases} \frac{\partial \omega}{\partial r} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial r} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial r}\\\frac{\partial \omega}{\partial \theta} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial \theta} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial \theta}\end{cases}?$
• Feb 11th 2011, 08:21 AM
bugatti79
Quote:

Originally Posted by Ackbeet
Could be way off here, but couldn't you just use this:

$\displaystyle \displaystyle\begin{cases} \frac{\partial \omega}{\partial r} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial r} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial r}\\\frac{\partial \omega}{\partial \theta} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial \theta} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial \theta}\end{cases}?$

Hi Ackbeet,

You could be right, I noticed my error...the dependant variable on the RHS should be u not omega!...I think
• Feb 11th 2011, 08:43 AM
Ackbeet
Right. And, of course, the

$\displaystyle \displaystyle\frac{\partial x}{\partial r},\;\frac{\partial x}{\partial \theta},\;\frac{\partial y}{\partial r},\;\frac{\partial y}{\partial \theta}$

are all easy to compute.
• Feb 11th 2011, 08:57 AM
bugatti79
Quote:

Originally Posted by bugatti79
Hi Ackbeet,

You could be right, I noticed my error...the dependant variable on the RHS should be u not omega!...I think

$\displaystyle \displaystyle\begin{cases} \frac{\partial \omega}{\partial r} = \frac{\partial u}{\partial x} Cos(\theta) + \frac{\partial u}{\partial y} Sin(\theta)\\\frac{\partial \omega}{\partial \theta} = -\frac{\partial u}{\partial x} r Sin(\theta) + \frac{\partial u}{\partial y} r Cos (\theta)\end{cases}$

The task is to rewrite the PDE and IC's in terms of $\displaystyle \omega,r, \theta$ but $\displaystyle \dispalysytle \frac{\partial u}{\partial x}$ and $\displaystyle \dispalysytle \frac{\partial u}{\partial y}$ both appear twice in the above equations. My idea was to make $\displaystyle \dispalysytle \frac{\partial u}{\partial x}$ the subject of both equations and then equate...I dont think that would work...I have ran out of ideas on how to solve this!
• Feb 11th 2011, 09:03 AM
Ackbeet
Quote:

Originally Posted by bugatti79
$\displaystyle \displaystyle\begin{cases} \frac{\partial \omega}{\partial r} = \frac{\partial u}{\partial x} Cos(\theta) + \frac{\partial u}{\partial y} Sin(\theta)\\\frac{\partial \omega}{\partial \theta} = -\frac{\partial u}{\partial x} r Sin(\theta) + \frac{\partial u}{\partial y} r Cos (\theta)\end{cases}$

The task is to rewrite the PDE and IC's in terms of $\displaystyle \omega,r, \theta$ but $\displaystyle \dispalysytle \frac{\partial u}{\partial x}$ and $\displaystyle \dispalysytle \frac{\partial u}{\partial y}$ both appear twice in the above equations. My idea was to make $\displaystyle \dispalysytle \frac{\partial u}{\partial x}$ the subject of both equations and then equate...I dont think that would work...I have ran out of ideas on how to solve this!

Yeah, I think it's going to take some tricks. How about this (since I see a rotation matrix present):

$\displaystyle \begin{bmatrix}\omega_{r}\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}\cos(\theta) &\sin(\theta)\\ -\sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}u_{x}\ \ u_{y}\end{bmatrix}.$

Inverting a rotation matrix is straight-forward: just transpose it! We get

$\displaystyle \begin{bmatrix}\cos(\theta) &-\sin(\theta)\\ \sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}\omega_{r }\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}u_{x}\\ u_{y}\end{bmatrix}.$

How's that?
• Feb 11th 2011, 10:34 AM
bugatti79
Quote:

Originally Posted by Ackbeet
Yeah, I think it's going to take some tricks. How about this (since I see a rotation matrix present):

$\displaystyle \begin{bmatrix}\omega_{r}\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}\cos(\theta) &\sin(\theta)\\ -\sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}u_{x}\ \ u_{y}\end{bmatrix}.$

Inverting a rotation matrix is straight-forward: just transpose it! We get

$\displaystyle \begin{bmatrix}\cos(\theta) &-\sin(\theta)\\ \sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}\omega_{r }\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}u_{x}\\ u_{y}\end{bmatrix}.$

How's that?

Wow, that looks good I must say...never saw this implemented in pde's before. Transposing of the matrix involves making the first row the first column and the second row the scond column. so if im right then

$\displaystyle u_x=-Sin(\theta)\omega_r+Cos(\theta)\omega_\theta /r$

$\displaystyle u_y=Cos(\theta)\omega_r+Sin(\theta)\omega_\theta /r$. Rewriting the pde becomes

$\displaystyle \displaystyle r Sin (\theta)[-Sin(\theta)\omega_r+Cos(\theta)\omega_\theta /r]-r Cos(\theta)[Cos(\theta)\omega_r+Sin(\theta)\omega_\theta /r]$=1

The IC being $\displaystyle \omega(r Cos(\theta),0)=0$ for $\displaystyle 0< r Cos(\theta)< \infty$

Is this correct...Can this be solved using the method of characteristics? ill try..
• Feb 11th 2011, 10:45 AM
Ackbeet
From

$\displaystyle \begin{bmatrix}\cos(\theta) &-\sin(\theta)\\ \sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}\omega_{r }\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}u_{x}\\ u_{y}\end{bmatrix},$

you should do regular matrix multiplication to arrive at

$\displaystyle u_{x}=\omega_{r}\cos(\theta)-\omega_{\theta}\sin(\theta)/r,$

$\displaystyle u_{y}=\sin(\theta)\omega_{r}+\omega_{\theta}\cos(\ theta)/r.$

Where does that take you?
• Feb 11th 2011, 10:53 AM
bugatti79
Quote:

Originally Posted by Ackbeet
From

$\displaystyle \begin{bmatrix}\cos(\theta) &-\sin(\theta)\\ \sin(\theta) &\cos(\theta)\end{bmatrix}\begin{bmatrix}\omega_{r }\\ \omega_{\theta}/r\end{bmatrix}=\begin{bmatrix}u_{x}\\ u_{y}\end{bmatrix},$

you should do regular matrix multiplication to arrive at

$\displaystyle u_{x}=\omega_{r}\cos(\theta)-\omega_{\theta}\sin(\theta)/r,$

$\displaystyle u_{y}=\sin(\theta)\omega_{r}+\omega_{\theta}\cos(\ theta)/r.$

Where does that take you?

Sorry I multiplied from the wrong side

$\displaystyle \displaystyle r Sin (\theta)[Cos (\theta)\omega_r-Sin(\theta)\omega_\theta /r]-r Cos(\theta)[Sin(\theta)\omega_r+Cos(\theta)\omega_\theta /r]$=1

I am not sure if this can be tackled using method of characteristics..I will try later
• Feb 11th 2011, 10:56 AM
Ackbeet
Assuming you've plugged everything in correctly, doesn't that simplify down quite a bit? You've got a cancellation, and if I'm not mistaken, the identity $\displaystyle \sin^{2}(\theta)+\cos^{2}(\theta)=1$ might be useful.
• Feb 11th 2011, 12:48 PM
bugatti79
Quote:

Originally Posted by Ackbeet
Assuming you've plugged everything in correctly, doesn't that simplify down quite a bit? You've got a cancellation, and if I'm not mistaken, the identity $\displaystyle \sin^{2}(\theta)+\cos^{2}(\theta)=1$ might be useful.

I think reduces quite considerably, I get

$\displaystyle \displaystyle \frac{\partial \omega}{\partial \theta}[Sin^2 \theta +Cos^2 \theta]+\frac{\partial \omega}{\partial r}[rSin \theta Cos \theta-rSin \theta Cos \theta]=-1$ therefore

$\displaystyle \omega (\theta)=-\theta+C$ would this translate to

$\displaystyle \omega (\theta)=u(x(\theta),y(\theta))$...?
• Feb 11th 2011, 12:50 PM
Ackbeet
No, I think when you integrate $\displaystyle \omega_{\theta}=-1,$ you're going to have to have

$\displaystyle \omega=-\theta+f(r),$ right?
• Feb 11th 2011, 01:01 PM
bugatti79
Quote:

Originally Posted by Ackbeet
No, I think when you integrate $\displaystyle \omega_{\theta}=-1,$ you're going to have to have

$\displaystyle \omega=-\theta+f(r),$ right?

Ok is that because the initial statement says that omega is a function of 2 variables even though $\displaystyle \frac{\partial\omega}{\partial \theta}=-1$ is explicitly a function of theta only?

From your equation using the IC $\displaystyle \omega(r Cos(\theta),0)=0$ i get

$\displaystyle 0=-rCos \theta+f(r)$ therefore $\displaystyle \omega(r, \theta)=- \theta+rCos \theta$

That theta doesnt look right on its own!...
• Feb 11th 2011, 01:05 PM
Ackbeet
It's because the integration you did was sort of "partial integration w.r.t. $\displaystyle y$", not a total integration like you'd have in an ODE setting. You remember doing exact ODE's? When you integrated w.r.t. $\displaystyle y$, you'd get an arbitrary function of $\displaystyle x$. In your case, $\displaystyle \omega=\omega(r,\theta),$ period. It is a function of two variables. If you doubt the extra function $\displaystyle f(r)$, then just do a partial differentiation and see if you get back to the original DE. You do, right?

What is the IC, exactly?
• Feb 11th 2011, 01:15 PM
bugatti79
Quote:

Originally Posted by Ackbeet
It's because the integration you did was sort of "partial integration w.r.t. $\displaystyle y$", not a total integration like you'd have in an ODE setting. You remember doing exact ODE's? When you integrated w.r.t. $\displaystyle y$, you'd get an arbitrary function of $\displaystyle x$. In your case, $\displaystyle \omega=\omega(r,\theta),$ period. It is a function of two variables. If you doubt the extra function $\displaystyle f(r)$, then just do a partial differentiation and see if you get back to the original DE. You do, right?

What is the IC, exactly?

No I agree with you..I actually had your expression intially but I somehow changed my mind

Its $\displaystyle u(x,0)=0 for 0<x<\infty$ which I believe translates to $\displaystyle \omega(r Cos(\theta),0)=0$ for $\displaystyle 0< r Cos(\theta)< \infty$
Show 40 post(s) from this thread on one page
Page 1 of 2 12 Last