A class of refined implicit-explicit Runge-Kutta methods with robust time adaptability and unconditional convergence for the Cahn-Hilliard model

One of main obstacles in verifying the energy dissipation laws of implicit-explicit Runge-Kutta (IERK) methods for phase field equations is to establish the uniform boundedness of stage solutions without the global Lipschitz continuity assumption of nonlinear bulk. With the help of discrete orthogonal convolution kernels, an updated time-space splitting technique is developed to establish the uniform boundedness of stage solutions for a refined class of IERK methods in which the associated differentiation matrices and the average dissipation rates are always independent of the time-space discretization meshes. This makes the refined IERK methods highly advantageous in self-adaptive time-stepping procedures as some larger adaptive step sizes in actual simulations become possible. From the perspective of optimizing the average dissipation rate, we construct some parameterized refined IERK methods up to third-order accuracy, in which the involved diagonally implicit Runge-Kutta methods for the implicit part have an explicit first stage and allow a stage-order of two such that they are not necessarily algebraically stable. Then we are able to establish, for the first time, the original energy dissipation law and the unconditional <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper L squared"> <mml:semantics> <mml:msup> <mml:mi>L</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">L^2</mml:annotation> </mml:semantics> </mml:math> </inline-formula> norm convergence. Extensive numerical tests are presented to support our theory.