0
$\begingroup$

I have a time iteration function looked on a 2D surface like this.

enter image description here Since the numbers wee very small i.e. hbar=6.6260700404e-34./(2*pi), my professor told me to use our own "scaled unites" during the calculation. ie. instead of using $1e-9$ (meter), scale the unit and use $1e-6$ ($\mu m$) during the calculation. e.t.c.

It worked very well when I was only changing the space parameter, (m) However, when I adjust the space scale and time scale together. I encounter a problem, where the function stopped being evolving smoothly.

enter image description here

I checked my substitution, and they were correct. But what happened here? Can we really save time of calculation by using scaled units?

|cite|improve this question
$\endgroup$
  • $\begingroup$ If you're solving Schrödinger equation, I suggest that you use a system of equations were constants like $h$, $m$, and $e$ are all unitary. $\endgroup$ – nicoguaro♦ Apr 8 '18 at 2:46
  • $\begingroup$ I think that your marching time method becomes unstable. $\endgroup$ – HBR Apr 8 '18 at 15:51
  • $\begingroup$ @HBR Why it became unstable? $\endgroup$ – J C Apr 8 '18 at 16:32
  • $\begingroup$ These peaks result when some iterative method becomes unstable, like when using euler forwards scheme the Courant number is more than 1 for a pure convection equation. Be sure that this number (ratio between $\Delta t$ and $$\delta x$) is under 1. $\endgroup$ – HBR Apr 8 '18 at 16:35
  • 2
    $\begingroup$ And for the last question. No. Scaling your nondimensional numbers allows you to control better the numerical aspects. It is always better to treat with numbers of the order of unity. $\endgroup$ – HBR Apr 8 '18 at 16:53
2
$\begingroup$

To close this question, as @AntonMenshov said in a comment I will put two of mine toguether:

These peaks result when some iterative method becomes unstable, like when using euler forwards scheme the Courant number defined as: $$C = \frac{v \Delta t}{\Delta x}$$ where $\Delta x/v$ must be the minimum residence time in a computation cell defined as the characteristic length of this cell $\Delta x$ divided by its local velocity $v$ and this number must be compared with the simulation timestep $\Delta t$. For an explicit time-scheme the Courant number is usually under 1 for a pure convection equation.

Also, you can improve the handling by using nondimensional numbers: scaling your dimensional quantities such as the potential or the time allows you to control better the numerical aspects.

|cite|improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.