Solving second order PDE by separation of variables (getting 2 ODE's)

[mathy_girl]

Hi all,

For my thesis I would like to solve the following second order nonlinear PDE for V(x,\sigma,t):

\frac{1}{2}\sigma^2\frac{\partial^2 V}{\partial x^2}+\frac{1}{2}B^2\frac{\partial^2 V}{\partial \sigma^2}+a\frac{\partial V}{\partial \sigma}=0,

subject to the following boundary conditions:
V(x,0,t)=\epsilon^{3/2}K \max(x,0)e^{-r(T-t)} and
V(x,\infty,t)=K (1+\epsilon^{3/2}x)

and with terminal condition
<br /> V(x,\sigma,T)=\epsilon^{3/2} K \max(x,0).<br />

I've tried separation of variables, by writing V=X(x)Y(\sigma)
which gives 2 ODEs
X&#039;&#039;(x)=kX(x) (which is nicely solvable) and
B^2Y&#039;&#039;(\sigma)+2aY&#039;(\sigma)+k\sigma^2Y(\sigma)=0, which is quite difficult, because it's still nonlinear.

Can anyone help me solve this? I don't know if this is the right way to do it, so other suggestions are also welcome!

Thanks!
Mathy_girl

PS: Help, I don't get the tex-code in my post in the right way... how does this work?

 

[HallsofIvy]

End with [/tex], not [\tex].

 

[HallsofIvy]

Hi all,

For my thesis I would like to solve the following second order nonlinear PDE for V(x,\sigma,t):

\begin{equation}\frac{1}{2}\sigma^2\frac{\partial^2 V}{\partial x^2}+\frac{1}{2}B^2\frac{\partial^2 V}{\partial \sigma^2}+a\frac{\partial V}{\partial \sigma}=0,\end{equation}

subject to the following boundary conditions:
\begin{equation} V(x,0,t)=\epsilon^{3/2}K \max(x,0)e^{-r(T-t)}\end{equation} and
\begin{equation}V(x,\infty,t)=K (1+\epsilon^{3/2}x)\end{equation}

and with terminal condition
<br /> \begin{equation}V(x,\sigma,T)=\epsilon^{3/2} K \max(x,0).\end{equation}<br />

I've tried separation of variables, by writing V=X(x)Y(\sigma)
which gives 2 ODEs
X&#039;&#039;(x)=kX(x) (which is nicely solvable) and
B^2Y&#039;&#039;(\sigma)+2aY&#039;(\sigma)+k\sigma^2Y(\sigma)=0, which is quite difficult, because it's still nonlinear.

No, that's a linear equation. It just has a variable coefficient. You might be able to do it with a series solution.

Can anyone help me solve this? I don't know if this is the right way to do it, so other suggestions are also welcome!

Thanks!
Mathy_girl

PS: Help, I don't get the tex-code in my post in the right way... how does this work?

 

[mathy_girl]

Thanks for your help with the tex-parts.

Unforunately, I already tried to construct a power series solution, and I didn't succeed so far.
This is what I tried:

We have to solve
\begin{equation}B^2Y&#039;&#039;(\sigma)+2aY&#039;(\sigma)+k\sigma^2Y(\sigma)=0.\label{ODE_y}\end{equation}

Expand Y(\sigma)=\sum_{n=0}^\infty c_n \sigma^n, and substitute this into (1) to get

B^2\sum_{n=2}^\infty n(n-1)c_n \sigma^{n-2}+2a \sum_{n=1}^\infty n c_n \sigma^{n-1}+k \sigma^2 \sum_{n=0}^\infty c_n \sigma^n = 0.

If we let all sums start at n=0, we get
\sum_{n=0}^\infty \left( B^2(n+2)(n+1)c_{n+2} + 2a (n+1) c_n + k \sigma^2 c_n \right) \sigma^n = 0.

Now how do I solve this? I need to find a recurrency relation for the coefficients c_n, but I don't know what to do with the \sigma^2 in the last term on the LHS.

 

[AiRAVATA]

The last term is not well stated, the equation for the powers should read

B^2\sum_{n=2}^\infty n(n-1)c_n \sigma^{n-2}+2a \sum_{n=1}^\infty n c_n \sigma^{n-1}+k \sum_{n=0}^\infty c_n \sigma^{n+2} = 0.

Now, changing coefficients so all powers are the same, we have

\sum_{n=0}^\infty \bigl( (n+2) B^2 c_{n+2} + 2a c_n\bigr)(n+1)\sigma^n + \sum_{n=2}^\infty k c_{n-2} \sigma^n= 0.

Can you take it from here?

 

[egust]

If a is zero, then your equation for Y is closely related to the quantum-mechanical simple harmonic oscillator in one dimension, and its polynomial solutions should be something like Gaussians multiplied by Hermite polynomials. Perhaps this can get you further along.

 

[mathy_girl]

The last term is not well stated, the equation for the powers should read

B^2\sum_{n=2}^\infty n(n-1)c_n \sigma^{n-2}+2a \sum_{n=1}^\infty n c_n \sigma^{n-1}+k \sum_{n=0}^\infty c_n \sigma^{n+2} = 0.

Now, changing coefficients so all powers are the same, we have

\sum_{n=0}^\infty \bigl( (n+2) B^2 c_{n+2} + 2a c_n\bigr)(n+1)\sigma^n + \sum_{n=2}^\infty k c_{n-2} \sigma^n= 0.

Can you take it from here?

When I take out the first two terms (for n=0 and n=1), we have
2B^2c_2 + 2a c_0 =0 and
(6B^62 c_3 + 4a c_1) \sigma =0,
such that
c_2 = -\frac{a}{B^2}c_0 and
c_3 = -\frac{2a}{3B^2}c_1.
Is that right?

So now both summation terms start at n=2, and can be combined to one summation as follows:

\sum_{n=2}^\infty ( ((n+2)B^2 c_{n+2} + 2a c_n)(n+1) + kc_{n-2} ) \sigma^n = 0

Next, I have to find this recurrency relation for all terms with n>1.. but how?
Is it given by

c_{n+2}=\frac{-2a(n+1) c_n - k c_{n-2}}{(n+1)(n+2)B^2}?

And how do I get my boundary conditions for X(x) and Y(\sigma)? Because we need them to determine the constants, right?

 

[AiRAVATA]

c_2 = -\frac{a}{B^2}a_0,

c_4 = \frac{2 \cdot 3 a^2-k}{3 \cdot 4 B^4}a_0,

and so on.

For the constants, I believe the best way is to express the boundary conditions in power series, but I'm not sure.

 

[mathy_girl]

c_2 = -\frac{a}{B^2}a_0,

c_4 = \frac{2 \cdot 3 a^2-k}{3 \cdot 4 B^4}a_0,

and so on.

For the constants, I believe the best way is to express the boundary conditions in power series, but I'm not sure.

Do you mean c_0 instead of a_0?

I understand how you got the first expression, for c_2, but how did you get the one for c_4??

 

[mathy_girl]

Anyone who can help me out with this? I'm still stuck at the same point...

 

[AiRAVATA]

Do you mean c_0 instead of a_0?

I understand how you got the first expression, for c_2, but how did you get the one for c_4??

Sorry for the late reply, didn't saw the post.

First of all, you are right, I meant c_0. Second, there is a little mistake in my algebra.

If

<br /> \sum_{n=0}^\infty \bigl( (n+2) B^2 c_{n+2} + 2a c_n\bigr)(n+1)\sigma^n + \sum_{n=2}^\infty k c_{n-2} \sigma^n= 0,<br />

then

<br /> c_2 = -\frac{a}{B^2} c_0,<br />

<br /> c_3 = -\frac{2a}{3B^2}c_1,<br />

and

<br /> c_{n+2}=- \frac{2 (n+1) a c_n + k c_{n-2}}{(n + 1) (n + 2) B^2}.<br />

So far so good. Now, For c_4,

<br /> c_4 = c_{2+2} = - \frac{2 (2 + 1) a c_2 + k c_0}{(2 + 1) (2 + 2) B^2} = - \frac{ 2 \cdot 3 a c_2 + k c_0}{ 3 \cdot 4 B^2}.<br />

Substituting c_2, then

<br /> c_4 = \frac{2 \cdot 3 a^2 - k B^2}{3 \cdot 4 B^4} c_0.<br />

In the same way you can calculate c_5,\,c_6,.... The reason why I let the constants expressed as multiplications is because it might be easier to spot factorials in the recurrence relation.

What is clear is that the general expression for c_{2n} is going to be complicated, as it will be a sum of n terms, while c_{2n+1} will be a sum of n-1 terms (i think).

Now that is clear how to calculate the coefficients, have you figured out how to deal with the boundary conditions?

P. S. Please redo all the math, just to make sure there's no errors this time.

 

[mathy_girl]

Thanks for your reply :) In the mean time, I indeed figured out how to calculate the coefficients, using the recursion, so on that I agree with you (or at least we do it the same way, tomorrow I'll check the outcomes with my own paperwork).

About the boundary conditions: still have no idea how to do this!

And about the recurrency: I don't expect to see a nice (well known) series or something in that. Do you think it's a logical step to just take the first x (some reasonable number) terms and use that as an approximation of the solution?

 

[AiRAVATA]

Well, how good the approximation is depends on how big the error of taking the truncated series. Remember that the solution is only valid whitin the radius of convergence of the series.

Let me ask you a question, is \epsilon a small parameter?

 

[mathy_girl]

Let me ask you a question, is \epsilon a small parameter?

Yes, it is: \epsilon is a small scaling parameter.

 

[AiRAVATA]

Then use asymptotic theory. The canonical reference would be https://www.amazon.com/dp/0387942025/?tag=pfamazon01-20, just to name a few.

 

[mathy_girl]

Then use asymptotic theory. The canonical reference would be https://www.amazon.com/dp/0387942025/?tag=pfamazon01-20, just to name a few.

I know, that's exactly what I've been doing. By applying this asymptotic theory (using the method of asymptotic expansions to solve a singular perturbation problem) I got to the problem I stated here... So applying it twice would be a bit strange, wouldn't it?

About the books: Holmes indeed is a good one, I've been using it. The other ones I don't really know, but thanks for the advice: maybe I'll take a look at them too...

 

[AiRAVATA]

The whole point of perturbation techniques is to simplify a problem by using the small parameter. How come you ended up with a simpler problem that is still dependent of the small parameter?

 

[mathy_girl]

Let me show you the complete problem then. Maybe that will make things a bit more clear. In the attachment you can find a part (scrap version of a chapter) of my report so far. Here I describe the problem completely. It's an application of perturbation methods on a financial model.

Section 1.3 is what this problem is about. Here we are doing some boundary layer analysis.

Hope you can help me out with this. If you know any other method to solve the inner problem: that is also fine!

 

[mathy_girl]

I'll try a Fourier transform on the original problem.. maybe that will help