NANov 8, 2016
Numerical validation of blow-up solutions of ordinary differential equationsAkitoshi Takayasu, Kaname Matsue, Takiko Sasaki et al.
This paper focuses on blow-up solutions of ordinary differential equations (ODEs). We present a method for validating blow-up solutions and their blow-up times, which is based on compactifications and the Lyapunov function validation method. The necessary criteria for this construction can be verified using interval arithmetic techniques. Some numerical examples are presented to demonstrate the applicability of our method.
NAFeb 5, 2019
Rigorous numerics of blow-up solutions for ODEs with exponential nonlinearityKaname Matsue, Akitoshi Takayasu
Our concerns here are blow-up solutions for ODEs with exponential nonlinearity from the viewpoint of dynamical systems and their numerical validations. As an example, the finite difference discretization of $u_t = u_{xx} + e^{u^m}$ with the homogeneous Dirichlet boundary condition is considered. Our idea is based on compactification of phase spaces and time-scale desingularization as in previous works. In the present case, treatment of exponential nonlinearity is the main issue. Fortunately, under a kind of exponential homogeneity of vector field, we can treat the problem in the same way as polynomial vector fields. In particular, we can characterize and validate blow-up solutions with their blow-up times for differential equations with such exponential nonlinearity in the similar way to previous works. A series of technical treatments of exponential nonlinearity in blow-up problems is also shown with concrete validation examples.
NAJul 23, 2018
Rigorous numerical computations for 1D advection equations with variable coefficientsAkitoshi Takayasu, Suro Yoon, Yasunori Endo
This paper provides a methodology of verified computing for solutions to 1-dimensional advection equations with variable coefficients. The advection equation is typical partial differential equations (PDEs) of hyperbolic type. There are few results of verified numerical computations to initial-boundary value problems of hyperbolic PDEs. Our methodology is based on the spectral method and semigroup theory. The provided method in this paper is regarded as an efficient application of semigroup theory in a sequence space associated with the Fourier series of unknown functions. This is a foundational approach of verified numerical computations for hyperbolic PDEs. Numerical examples show that the rigorous error estimate showing the well-posedness of the exact solution is given with high accuracy and high speed.
NAJul 30, 2017
Accurate method of verified computing for solutions of semilinear heat equationsAkitoshi Takayasu, Makoto Mizuguchi, Takayuki Kubo et al.
We provide an accurate verification method for solutions of heat equations with a superlinear nonlinearity. The verification method numerically proves the existence and local uniqueness of the exact solution in a neighborhood of a numerically computed approximate solution. Our method is based on a fixed-point formulation using the evolution operator, an iterative numerical verification scheme to extend a time interval in which the validity of the solution can be verified, and rearranged error estimates for avoiding the propagation of an overestimate. As a result, compared with the previous verification method using the analytic semigroup, our method can enclose the solution for a longer time. Some numerical examples are presented to illustrate the efficiency of our verification method.
NAJul 19, 2017
Numerical validation of blow-up solutions with quasi-homogeneous compactificationsKaname Matsue, Akitoshi Takayasu
We provide a numerical validation method of blow-up solutions for finite dimensional vector fields admitting asymptotic quasi-homogeneity at infinity. Our methodology is based on quasi-homogeneous compactifications containing a new compactification, which shall be called a quasi-parabolic compactification. Divergent solutions including blow-up solutions then correspond to global trajectories of associated vector fields with appropriate time-variable transformation tending to equilibria on invariant manifolds representing infinity. We combine standard methodology of rigorous numerical integration of differential equations with Lyapunov function validations around equilibria corresponding to divergent directions, which yields rigorous upper and lower bounds of blow-up times as well as rigorous profile enclosures of blow-up solutions.