Waveform relaxation methods for implicit differential equations

We apply a Runge-Kutta-based waveform relaxation method to initial-value problems for implicit differential equations. In the implementation of such methods, a sequence of nonlinear systems has to be solved iteratively in each step of the integration process. The size of these systems increases linearly with the number of stages of the underlying Runge-Kutta method, resulting in high linear algebra costs in the iterative process for high-order Runge-Kutta methods. In our earlier investigations of iterative solvers for implicit initial-value problems, we designed an iteration method in which the linear algebra costs are almost independent of the number of stages when implemented on a parallel computer system. In this paper, we use this parallel iteration process in the Runge-Kutta waveform relaxation method. In particular, we analyse the convergence of the method. The theoretical results are illustrated by a few numerical examples.     

Local solvability for a class of evolution equations

Motivated by the celebrated example of Y. Kannai of a linear partial differential operator which is hypoelliptic but not locally solvable, we consider a class of evolution operators with real-analytic coefficients and study their local solvability both in L² and in the weak sense. In order to do so we are led to propose a generalization of the Nirenberg-Treves condition (ψ) which is suitable to our study.     

On the solvability of implicit differential inclusions

In the paper, the set-valued covering mappings are studied. The statements on solvability, solution estimates, and well-posedness of inclusions with conditionally covering mappings are proved. The results obtained are applied to the investigation of differential inclusions unsolved for the unknown function. The statements on solvability, solution estimates, and well-posedness of these inclusions are derived.     

On the iterative solution of the algebraic equations in fully implicit Runge-Kutta methods

This paper deals with the iterative solution of stage equations which arise when some fully implicit Runge-Kutta methods, in particular those based on Gauss, Radau and Lobatto points, are applied to stiff ordinary differential equations. The error behaviour in the iterates generated by Newton-type and, particularly, by single-Newton schemes which are proposed for the solution of stage equations is studied. We consider stiff systems y'(t) = f(t,y(t)) which are dissipative with respect to a scalar product and satisfy a condition on the relative variation of the Jacobian of f(t,y) with respect to y, similar to the condition considered by van Dorsselaer and Spijker in [7] and [17]. We prove new convergence results for the single-Newton iteration and derive estimates of the iteration error that are independent of the stiffness. Finally, some numerical experiments which confirm the theoretical results are presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

Non-instantaneous impulsive fractional-order implicit differential equations with random effects

In this article, we study existence and stability of a class of non-instantaneous impulsive fractional-order implicit differential equations with random effects. First, we establish a framework to study impulsive fractional sample path associated with impulsive fractional L^p-problem, and present the relationship between them. We also derive the formula of the solution for inhomogeneous impulsive fractional L^p-problem and sample path. Second, we construct a sequence of Picard functions, which admits us to apply successive approximations method to seek the solution of impulsive fractional sample path. Further, we derive the existence of solutions to impulsive fractional L^p-problem. Third, the concepts of Ulam's type stability are introduced and sufficient conditions to guarantee Ulam–Hyers–Rassias stability are derived. Finally, an example is given to illustrate the theoretical results.     

Nonlocal boundary value problem for generalized Hilfer implicit fractional differential equations

In this paper, we derive the equivalent fractional integral equation to the nonlinear implicit fractional differential equations involving Ψ-Hilfer fractional derivative subject to nonlocal fractional integral boundary conditions. The existence of a solution, Ulam–Hyers, and Ulam–Hyers–Rassias stability have been acquired by means of an equivalent fractional integral equation. Our investigations depend on the fixed-point theorem due to Krasnoselskii and the Gronwall inequality involving Ψ-Riemann–Liouville fractional integral. Finally, examples are provided to show the utilization of primary outcomes.     

Generalized implicit Euler method for hyperbolic functional differential equations

Nonlinear hyperbolic functional differential equations with initial boundary conditions are considered. Theorems on the convergence of difference schemes and error estimates of approximate solutions are presented. The proof of the stability of the difference functional problem is based on a comparison technique. Nonlinear estimates of the Perron type with respect to the functional variable for given functions are used. Numerical examples are given (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Some remarks on planar Boussinesq equations

The main purpose of this paper is to prove the well-posedness of the two-dimensional Boussinesq equations when the initial vorticity ω 0 ∈L1 (R 2 ) (or the finite Radon measure space). Using the stream function form of the equations and the Schauder fixed-point theorem to get the new proof of these results, we get that when the initial vorticity is smooth, there exists a unique classical solutions for the Cauchy problem of the two dimensional Boussinesq equations.     

Integral equations, implicit functions, and fixed points

The problem is to show that (1) has a solution, where defines a contraction, , and defines a compact map, . A fixed point of would solve the problem. Such equations arise naturally in the search for a solution of where , but so that the standard conditions of the implicit function theorem fail. Now would be in the form for a classical fixed point theorem of Krasnoselskii if were a contraction. But fails to be a contraction for precisely the same reasons that the implicit function theorem fails. We verify that has enough properties that an extension of Krasnoselskii's theorem still holds and, hence, (1) has a solution. This substantially improves the classical implicit function theorem and proves that a general class of integral equations has a solution.

     

Local solvability for semilinear Fuchsian equations

We consider a semilinear elliptic operator P on a manifold B with a conical singular point. We assume P is Fuchs type in the linear part and has a non–linear lower order therms. Using the Schauder fixed point theorem, we prove the local solvability of P near the conical point in the weighted Sobolev spaces.     

Fourth-order diagonally implicit schemes for evolution equations

New fourth-order methods are proposed for solving both ordinary and partial differential equations. The derivation of the methods is based on the form of diagonally implicit schemes applied to stiff ordinary differential equations. The methods are absolutely and unconditionally stable. Test computations are presented.     

Solvability of partial differential equations by meshless kernel methods

This paper first provides a common framework for partial differential equation problems in both strong and weak form by rewriting them as generalized interpolation problems. Then it is proven that any well-posed linear problem in strong or weak form can be solved by certain meshless kernel methods to any prescribed accuracy. The work described in this paper was partially supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 101205). Robert Schaback’s research in Hong Kong was sponsored by DFG and City University of Hong Kong.     

Linearly implicit methods for nonlinear parabolic equations

We construct and analyze combinations of rational implicit and explicit multistep methods for nonlinear parabolic equations. The resulting schemes are linearly implicit and include as particular cases implicit-explicit multistep schemes as well as the combination of implicit Runge-Kutta schemes and extrapolation. An optimal condition for the stability constant is derived under which the schemes are locally stable. We establish optimal order error estimates.

     

Some remarks on extremal solutions of multivalued differential equations

In this paper, we consider extremal solutions of multivalued differential equations, i.e., solutions that steer to the boundary of the attainable set. Multivalued differential equations arise in a natural way from control systems governed by ordinary differential equations that have a variable control-constraint set. Extremal solutions of multi-valued differential equations are important in the study of the optimal control of such systems. We give conditions under which extremality of a solution at a certain time implies extremality of the solution at all previous times where it is defined. Necessary conditions for extremality are also obtained. We treat both the time-dependent case and the time-independent case.     

Stein implicit Runge-Kutta methods with high stage order for large-scale ordinary differential equations

The Runge-Kutta method is one of the most popular implicit methods for the solution of stiff ordinary differential equations. For large problems, the main drawback of such methods is the cost required at each integration step for computing the solution of a nonlinear system of equations. In this paper, we propose to reduce the cost of the computation by transforming the linear systems arising in the application of Newton's method to Stein matrix equations. We propose an iterative projection method onto block Krylov subspaces for solving numerically such Stein matrix equations. Numerical examples are given to illustrate the performance of our proposed method.     

Completeness of security markets and solvability of linear backward stochastic differential equations

For a standard Black-Scholes type security market, completeness is equivalent to the solvability of a linear backward stochastic differential equation (BSDE, for short). An ideal case is that the interest rate is bounded, there exists a bounded risk premium process, and the volatility matrix has certain surjectivity. In this case the corresponding BSDE has bounded coefficients and it is solvable leading to the completeness of the market. However, in general, the risk premium process and/or the interest rate could be unbounded. Then the corresponding BSDE will have unbounded coefficients. For this case, do we still have completeness of the market? The purpose of this paper is to discuss the solvability of BSDEs with possibly unbounded coefficients, which will result in the completeness of the corresponding market.