- #1

- 48

- 0

## Main Question or Discussion Point

Hi all!

I am coming back here after a long time. Last time I got the answer I was looking for, here. I hope that I will find it again.

I need the possible numerical methods for solving the following PDE: Its a nonlinear parabolic boundary value problem. I want a very stable algorithm as this equation is a part of a bigger code. I don’t want accumulated numerical errors as I am investigating an instability in the bigger code involving a lot of iterations.

The PDE is :

d(N(x,t))/dt=D*d2(N(x,t))/dx2-A*N(x,t)-B*N2(x,t)-C*N3(x,t)-S(x,t)+ R(x)

N(x=a,t)=N(x=b,t)=0 for all t. N(x,t=0)=f(x) . (f(x) is a Gaussian function)

where N2 is N-squared, N3 is N-cube and so on. S(x,t) can be a rapidly oscillating function in x and t. C is very small as compared to the other constants A, B and D (all positive).

Is it possible to use a relaxation method for solving the above PDE? I ask this because the corresponding steady state problem can be solved very effectively using a Relaxation method and I have a efficient code for that already.

Please also give me a reference with your suggestions so that I can refer it. Please ask me if you need any more information about the PDE.

Thanks for the help in advance and for having such a great forum.

Jam

I am coming back here after a long time. Last time I got the answer I was looking for, here. I hope that I will find it again.

I need the possible numerical methods for solving the following PDE: Its a nonlinear parabolic boundary value problem. I want a very stable algorithm as this equation is a part of a bigger code. I don’t want accumulated numerical errors as I am investigating an instability in the bigger code involving a lot of iterations.

The PDE is :

d(N(x,t))/dt=D*d2(N(x,t))/dx2-A*N(x,t)-B*N2(x,t)-C*N3(x,t)-S(x,t)+ R(x)

N(x=a,t)=N(x=b,t)=0 for all t. N(x,t=0)=f(x) . (f(x) is a Gaussian function)

where N2 is N-squared, N3 is N-cube and so on. S(x,t) can be a rapidly oscillating function in x and t. C is very small as compared to the other constants A, B and D (all positive).

Is it possible to use a relaxation method for solving the above PDE? I ask this because the corresponding steady state problem can be solved very effectively using a Relaxation method and I have a efficient code for that already.

Please also give me a reference with your suggestions so that I can refer it. Please ask me if you need any more information about the PDE.

Thanks for the help in advance and for having such a great forum.

Jam