[论文解读] A Modified Center-of-Mass Conservation Law in Finite-Domain Simulations of the Zakharov--Kuznetsov Equation

Nobuyuki Sawado, Yuichiro Shimazaki|arXiv (Cornell University)|Jan 22, 2026

Nonlinear Waves and Solitons被引用 0

一句话总结

该论文表明有限域模拟中的边界通量会导致ZK方程的质心分量 I4x 的漂移，并引入修正的守恒量 I4x^mod 以恢复守恒并纠正质心动力学。

ABSTRACT

We investigate conservation laws of the two-dimensional Zakharov--Kuznetsov (ZK) equation, a natural higher-dimensional and non-integrable extension of the Korteweg--de Vries equation. The ZK equation admits three scalar conserved quantities -- mass, momentum, and energy -- represented as $I_1$, $I_2$, and $I_3$, as well as a vector-valued quantity $\bm{I}_4$. In high-accuracy numerical simulations on a finite double-periodic domain, most of these quantities are well preserved, while a systematic temporal drift is observed only in the $x$-component $I_{4x}$. We show that the nontrivial evolution of $I_{4x}$ originates from an explicit boundary-flux contribution, which is induced by fluctuations of the solution and its spatial derivatives at the domain boundaries. We successfully identify the source of the inaccuracy in the numerical solutions. Motivated by this analysis, we define a modified center-of-mass quantity $I_{4x}^{\mathrm{mod}}$ and demonstrate its conservation numerically for single-pulse configurations. The modified quantity thus provides a consistent conservation law for the ZK equation and yields an appropriate description of center-of-mass motion in finite-domain numerical simulations.

研究动机与目标

• Motivate and analyze conservation laws for the two-dimensional Zakharov–Kuznetsov (ZK) equation, a higher-dimensional generalization of KdV.
• Identify which conserved quantities are robust under finite-domain, double-periodic boundary conditions.
• Explain the boundary-flux mechanism responsible for the drift in the center-of-mass component I4x and propose a corrective invariant.
• Demonstrate that the modified center-of-mass quantity I4x^mod remains conserved in numerical experiments for localized pulses.

提出的方法

• Use high-accuracy Fourier pseudo-spectral discretization on a finite doubly periodic domain.
• Represent the ZK equation as u_t + ∂x Q = 0 with Q = 6u^2 + Δu to identify boundary flux through domain edges.
• Compute time derivatives of I4x and decompose into boundary terms to isolate flux contributions (A(t), B(t), C(t), D(t)).
• Show analytically that on finite domains, B, C, D vanish due to periodic sampling, leaving boundary flux A(t) as the source of drift.
• Define a modified center-of-mass invariant I4x^mod(t) = I4x(t) − ∫0^t A(τ) dτ and verify its conservation numerically for single-pulse solutions.
• Examine centroid coordinates and velocities derived from invariants to connect conserved quantities with physical centroid motion.

实验结果

研究问题

• RQ1 Does finite-domain periodic truncation preserve all ZK conservation laws, and if not, what causes deviations?
• RQ2 What is the explicit boundary-flux mechanism responsible for the drift of I4x in finite-domain simulations?
• RQ3 Can a modified conserved quantity be defined to restore exact center-of-mass conservation on finite domains, and does it yield physically meaningful centroid dynamics?
• RQ4 How does the modified centroid formulation affect inferred centroid position and especially centroid velocity for localized pulse solutions?

主要发现

• I1 (mass), I2 (momentum), I3 (energy), and I4y remain conserved to numerical accuracy in finite-domain simulations.
• The x-component of the center-of-mass invariant I4x exhibits a systematic temporal drift due to boundary-flux contributions from weak radiative tails.
• The analytically derived boundary flux A(t) exactly matches the observed drift in I4x over time.
• A modified center-of-mass quantity I4x^mod, defined by subtracting the accumulated boundary flux, remains nearly constant for single-pulse configurations.
• Centroid coordinates and velocities analyzed via invariants show that using I4x^mod yields a uniform centroid motion, aligning with intrinsic translational dynamics.
• Appendix results indicate similar boundary-flux effects can occur in the KdV equation, suggesting broader applicability of the modification.