Shun Sato

Daisuke Furihata, Shun Sato, Takayasu Matsuo

We consider structure-preserving methods for conservative systems, which rigorously replicate the conservation property yielding better numerical solutions. There, corresponding to the skew-symmetry of the differential operator, that of difference operators is essential to the discrete conservation law. Unfortunately, however, when we employ the standard central difference operator, the simplest one, the numerical solutions often suffer from undesirable spatial oscillations. In this letter, we propose a novel "average-difference method," which is tougher against such oscillations, and combine it with an existing conservative method. Theoretical and numerical analysis in the linear case show the superiority of the proposed method.

Shun Sato, Takayasu Matsuo

Recently, various evolutionary partial differential equations (PDEs) with a mixed derivative have been emerged and drawn much attention. Nonetheless, their PDE-theoretical and numerical studies are still in their early stage. In this paper, we aim at the unified framework of numerical methods for such PDEs. However, due to the presence of the mixed derivative, we cannot discuss numerical methods without some appropriate reformulation, which is mathematically challenging itself. Therefore, we first propose a novel procedure for the reformulation of target PDEs into a standard form of evolutionary equations. This contribution may become an important basis not only of numerical analysis, but also of PDE-theory. In order to illustrate this point, we establish the global well-posedness of the sine-Gordon equation. After that, we classify and discuss the spatial discretizations based on the proposed reformulation technique. As a result, we show the average-difference method is suitable for the discretization of the mixed derivative.

Kaito Sato, Shun Sato, Takayasu Matsuo

Multi-symplectic diamond schemes proposed by McLachlan and Wilkins (2015) provide a framework for the numerical integration of Hamiltonian partial differential equations, combining local implicitness with high-order accuracy and discrete multi-symplectic conservation laws. Despite these advantages, their behavior beyond a limited class of model equations remains poorly understood, and numerical difficulties may arise depending on the underlying multi-symplectic formulation. In this paper, we present a systematic stability analysis framework for diamond schemes applied to general multi-symplectic PDEs. The approach consists of three stages. First, we identify structural inconsistency of the local diamond update using Dulmage--Mendelsohn decomposition, revealing cases in which the scheme is intrinsically unsolvable. Second, we introduce a graph-based error-propagation analysis that yields a necessary stability condition by detecting negative cycles in a weighted directed graph. Third, for equations that pass the preliminary tests, we derive eigenvalue-based timestep restrictions providing sufficient conditions for stability. The analysis leads to a comprehensive classification of multi-symplectic PDEs according to whether diamond schemes are structurally inconsistent, unconditionally unstable, or conditionally stable. In particular, we show that benchmark equations such as the Korteweg--de Vries equation are intrinsically incompatible with the diamond update, while systems including the nonlinear Dirac and ``good'' Boussinesq equations admit stability regimes under mild timestep scaling. Extensive numerical experiments confirm the theoretical predictions and demonstrate the practical implications of the proposed framework. Our results clarify fundamental limitations of diamond schemes and provide practical guidelines for their reliable application to new PDE models.

S. Sato, K. Oguma, T. Matsuo et al.

The short pulse equation was introduced by Schaefer--Wayne (2004) for modeling the propagation of ultrashort optical pulses. While it can describe a wide range of solutions, its ultrashort pulse solutions with a few cycles, which the conventional nonlinear Schroedinger equation does not possess, have drawn much attention. In such a region, existing numerical methods turn out to require very fine numerical mesh, and accordingly are computationally expensive. In this paper, we establish a new efficient numerical method by combining the idea of the hodograph transformation and the structure-preserving numerical methods. The resulting scheme is a self-adaptive moving mesh scheme that can successfully capture not only the ultrashort pulses but also exotic solutions such as loop soliton solutions.