Error in Numerical Solution of ODE by Euler Method - Patrick

patrick1990
Messages
2
Reaction score
0
Hi,
I recently need to do some numerical simulation by Euler method to solve a PDE.
However, I noticed that there are some errors which are obtained with bigger numerical steps, when applying Euler scheme.
Since my major is not mathematics, I do not know what this phenomenon is called. I have read it somewhere else (from Wikipedia ?), but unfortunately I cannot recall it at all.
Anyone knows the name of this ?
Thank you so much !
patrick
 

Attachments

  • euler.png
    euler.png
    2.5 KB · Views: 619
Physics news on Phys.org


It is generally called numerical instabilities or divergence and occurs because the state of discrete system (the solver) does not converge to the solution of the continuous system. The Euler method is a very simple method with a very small region of absolute stability. In general one usually achieves much faster or more precise (or both) solutions by integrating using a more capable method.

For more information you may want to search your references for "Euler region of absolute stability".
 


Thank you so much !
:)
 
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top