We show that while Runge-Kutta methods cannot preserve polynomial invariants in general, they can preserve polynomials that are the energy invariant of canonical Hamiltonian systems.

We prove the existence of solutions to nonlinear parabolic problems of the following type: $$\left\{\begin{array}{cc}{\displaystyle \frac{\partial b\left(u\right)}{\partial t}}+A\left(u\right)=f+\mathrm{div}\left(\Theta (x;t;u)\right)\hfill & \text{in}Q,\hfill \\ u(x;t)=0\hfill & \text{on}\partial \Omega \times [0;T],\hfill \\ b\left(u\right)(t=0)=b\left(u_0\right)\hfill & \text{on}\Omega ,\hfill \end{array}\right.$$ where $b:ℝ\to ℝ$ is a strictly increasing function of class ${𝒞}^1$, the term $$A\left(u\right)=-\mathrm{div}\left(a\right(x,t,u,\nabla u\left)\right)$$ is an operator of Leray-Lions type which satisfies the classical Leray-Lions assumptions of Musielak type, $\Theta :\Omega \times [0;T]\times ℝ\to ℝ$ is a Carathéodory, noncoercive function which satisfies the following condition: ${sup}_{\left|s\right|\le k}\left|\Theta (·,·,s)\right|\in E_{\psi }\left(Q\right)$ for all $k>0$, where $\psi$ is the Musielak complementary function of $\Theta$, and the second term $f$ belongs to $L^1\left(Q\right)$.

For contractive interval functions $\left[g\right]$ we show that $\left[g\right]\left({\left[x\right]}_{ϵ}^{k_0}\right)\subseteq \int \left({\left[x\right]}_{ϵ}^{k_0}\right)$ results from the iterative process ${\left[x\right]}^{k+1}:=\left[g\right]\left({\left[x\right]}_{ϵ}^k\right)$ after finitely many iterations if one uses the epsilon-inflated vector ${\left[x\right]}_{ϵ}^k$ as input for $\left[g\right]$ instead of the original output vector ${\left[x\right]}^k$. Applying Brouwer’s fixed point theorem, zeros of various mathematical problems can be verified in this way.

In this work, the error behaviour of high-order exponential operator splitting methods for the time integration of nonlinear evolutionary Schrödinger equations is investigated. The theoretical analysis utilises the framework of abstract evolution equations on Banach spaces and the formal calculus of Lie derivatives. The general approach is substantiated on the basis of a convergence result for exponential operator splitting methods of (nonstiff) order p applied to the multi-configuration time-dependent...

A pointwise error estimate and an estimate in norm are obtained for a class of external methods approximating boundary value problems. Dependence of a superconvergence phenomenon on the external approximation method is studied. In this general framework, superconvergence at the knot points for piecewise polynomial external methods is established.

Galerkin reduced-order models for the semi-discrete wave equation, that preserve the second-order structure, are studied. Error bounds for the full state variables are derived in the continuous setting (when the whole trajectory is known) and in the discrete setting when the Newmark average-acceleration scheme is used on the second-order semi-discrete equation. When the approximating subspace is constructed using the proper orthogonal decomposition, the error estimates are proportional to the sums...