2
$\begingroup$

We can solve some differential equation numerically using Runge-Kutta methods. i.e. for some function $f(x,y) = \frac{dy}{dx}$, then

$$ y_{n+1} = y_n + \sum_{i=1}^{s} c_i k_i $$

where e.g.

$$ k_2 = f(x_n+a_2 h, y_n + b_{21}k_1)$$

and $a_i$, $b_{ij}$ are some coefficeints e.g. Cash-Karp.


Now, I understand that there also exist symplectic Runge-Kutta methods. How would this method work? What would be the equivalent $y_{n+1}$ style equations for a symplectic Runge-Kutta? How are the coefficients chosen?

$\endgroup$
2

1 Answer 1

1
$\begingroup$

The following is sufficient and necessary for a Runge-Kutta method to be symplectic: $$ b_i a_{i,j} + b_j a_{j,i} - b_ib_j = 0 \, \quad 1 \leq i,j \leq s $$ where $s$ is the number of stages in your method. No explicit RK methods achieve this.

This condition is in addition to other consistency conditions too.

Also, don't forget the $dt$ in your method, and the convention is to use $b_i$ in the linear combination of stage derivatives $k_i$, $c_i$ multiplied by $h=dt$, and $a_{i,j}$ for the states.

source: http://www.sanzserna.org/pdf/56_physicaD.pdf

EDIT: As pointed out in the comments, this applies to symplectic Runge-Kutta methods, but generally methods may be symplectic but not Runge-Kutta.

$\endgroup$
2
  • $\begingroup$ This applies to the use of standard RK methods, like the midpoint method. A larger number of symplectic methods are constructed as partitioned methods, with different, but coupled stages for the position and impulse components of the system. $\endgroup$ Commented Mar 14, 2022 at 17:02
  • 1
    $\begingroup$ Indeed. But those methods are not Runge-Kutta per se, so I am answering the OP. $\endgroup$ Commented Mar 14, 2022 at 19:30

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .