Limit Theorem for Controlled Backward SDEs and Homogenization of Hamilton–Jacobi–Bellman Equations

We prove a convergence theorem for a family of value functions associated with stochastic control problems whose cost functions are defined by backward stochastic differential equations. The limit function is characterized as a viscosity solution to a fully nonlinear partial differential equation of second order. The key assumption we use in our approach is shown to be a necessary and sufficient assumption for the homogenizability of the control problem. The results generalize partially homogenization problems for Hamilton–Jacobi–Bellman equations treated recently by Alvarez and Bardi by viscosity solution methods. In contrast to their approach, we use mainly probabilistic arguments, and discuss a stochastic control interpretation for the limit equation.     

Low-storage Runge–Kutta methods for stochastic differential equations

Runge–Kutta methods that require only two memory locations per variable and have strong local order γ=1.5 for non-commutative systems of stochastic differential equations driven by one Wiener process are devised in this paper. A first step in the derivation is to extend existing deterministic methods to the commutative stochastic case, for which higher accuracy is also obtained. Numerical results are presented to validate the approach.     

On Convergence of One-Step Schemes for Weak Solutions of Quantum Stochastic Differential Equations

Several one-step schemes for computing weak solutions of Lipschitzian quantum stochastic differential equations (QSDE) driven by certain operator-valued stochastic processes associated with creation, annihilation and gauge operators of quantum field theory are introduced and studied. This is accomplished within the framework of the Hudson–Parthasarathy formulation of quantum stochastic calculus and subject to the matrix elements of solution being sufficiently differentiable. Results concerning convergence of these schemes in the topology of the locally convex space of solution are presented. It is shown that the Euler–Maruyama scheme,with respect to weak convergence criteria for Itô stochastic differential equation is a special case of Euler schemes in this framework. Numerical examples are given.     

Stochastic PDEs Driven by Nonlinear Noise and Backward Doubly SDEs

We study a new kind of backward doubly stochastic differential equations, where the nonlinear noise term is given by Itô–Kunita's stochastic integral. This allows us to give a probabilistic interpretation of classical and Sobolev's solutions of semilinear parabolic stochastic partial differential equations driven by a nonlinear space-time noise.     

Runge–Kutta–collocation methods for systems of functional–differential and functional equations

Systems of functional–differential and functional equations occur in many biological, control and physics problems. They also include functional–differential equations of neutral type as special cases. Based on the continuous extension of the Runge–Kutta method for delay differential equations and the collocation method for functional equations, numerical methods for solving the initial value problems of systems of functional–differential and functional equations are formulated. Comprehensive analysis of the order of approximation and the numerical stability are presented.     

Almost Sure Convergence of the Numerical Discretization of Stochastic Jump Diffusions

Based on the shuffle product expansion of exponential Lie series in terms of a Philip Hall basis for the stochastic differential equations of jump-diffusion type, we can establish Stratonovich–Taylor–Hall (STH) schemes. However, the STHr scheme converges only at order r in the mean-square sense. In order to have the almost sure Stratonovich–Taylor–Hall (ASTH) schemes, we have to include all the terms related to multiple Poissonian integrals as the moments of multiple Poissonian integrals always have lower orders of magnitudes as compared with those of multiple Brownian integrals.     

Robust Parameter Estimation for Stochastic Differential Equations

We consider the estimation of parameters in stochastic differential equations (SDEs). The problem is treated in the setting of nonlinear filtering theory with a degenerate diffusion matrix. A robust stochastic Feynman–Kac representation for solutions of SDEs of Zakai-type is derived. It is verified that these solutions are conditional densities for the conditional measures defined by degenerate filtering problems. We show that the corresponding estimator for the parameters is robust in the following sense: It depends continuously on both the measurement path and on the intensity of the measurement noise. An algorithm based on a Monte-Carlo approach is given for the practical application of the estimator, and numerical results are reported. Mathematics Subject Classifications (2000) Primary: 62M05, 62M20; secondary: 62F15.     

Max-Plus Stochastic Processes

This paper is concerned with processes which are max-plus counterparts of Markov diffusion processes governed by Ito sense stochastic differential equations. Concepts of max-plus martingale and max-plus stochastic differential equation are introduced. The max-plus counterparts of backward and forward PDEs for Markov diffusions turn out to be first-order PDEs of Hamilton–Jacobi–Bellman type. Max-plus additive integrals and a max-plus additive dynamic programming principle are considered. This leads to variational inequalities of Hamilton–Jacobi–Bellman type.     

Conservation laws for Camassa–Holm equation, Dullin–Gottwald–Holm equation and generalized Dullin–Gottwald–Holm equation

In this paper we construct the conservation laws for the Camassa–Holm equation, the Dullin–Gottwald–Holm equation (DGH) and the generalized Dullin–Gottwald–Holm equation (generalized DGH). The variational derivative approach is used to derive the conservation laws. Only first order multipliers are considered. Two multipliers are obtained for the Camassa–Holm equation. For the DGH and generalized DGH equations the variational derivative approach yields two multipliers; thus two conserved vectors are obtained.     

Exact travelling wave solutions for the generalized nonlinear Schrödinger (GNLS) equation with a source by Extended tanh–coth, sine–cosine and Exp-Function methods

The capability of Extended tanh–coth, sine–cosine and Exp-Function methods as alternative approaches to obtain the analytic solution of different types of applied differential equations in engineering mathematics has been revealed. In this study, the generalized nonlinear Schrödinger (GNLS) equation is solved by three different methods. To obtain the single-soliton solutions for the equation, the Extended tanh–coth and sine–cosine methods are used. Furthermore, for this nonlinear evolution equation the Exp-Function method is applied to derive various travelling wave solution. Results show that while the first two procedures easily provide a concise solution, the Exp-Function method provides a powerful mathematical means for solving nonlinear evolution equations in mathematical physics.     

Runge–Kutta methods for first-order periodic boundary value differential equations with piecewise constant arguments

In this paper we deal with the numerical solutions of Runge–Kutta methods for first-order periodic boundary value differential equations with piecewise constant arguments. The numerical solution is given by the numerical Green’s function. It is shown that Runge–Kutta methods preserve their original order for first-order periodic boundary value differential equations with piecewise constant arguments. We give the conditions under which the numerical solutions preserve some properties of the analytic solutions, e.g., uniqueness and comparison theorems. Finally, some experiments are given to illustrate our results.     

Fourier–Bessel-type series: The fourth-order case

The structured higher-order Bessel-type linear ordinary differential equations were first discovered in 1994. There is a denumerable infinity of these higher-order equations, all of then of even-order.These differential equations possess many of the properties of the classical second-order Bessel differential equation, but these higher-order cases bring remarkable new analytic structures. In many ways it is sufficient to study the properties of the fourth-order Bessel-type differential equation to be able to assess the corresponding properties of the sixth-and higher-order cases.This paper follows a number of earlier papers devoted to the study of the fourth-order case. These publications show the connections between the special function properties of solutions of the differential equation, and the properties of linear differential operators generated by the associated linear differential expression in certain weighted Lebesgue, and Lebesgue–Stieltjes function spaces.To follow the earlier papers on the study of the fourth-order Bessel-type differential equation, this present paper determines the form of the Fourier–Bessel-type series which best extends the classical theory of the second-order Fourier–Bessel series.In fact the Fourier–Bessel-type series are based on a new orthogonal system in terms of the regular eigensolutions of the fourth-order Bessel-type equation. The corresponding eigenvalues are obtained by restricting the spectral parameter to the zeros of an analytic function arising already in the Dini boundary conditions.     

Mixed Laguerre–Legendre Spectral Method for Incompressible Flow in an Infinite Strip

This paper concerns the mixed Laguerre–Legendre spectral approximation and its application to numerical simulation of incompressible flow in an infinite strip. Some approximation results in weighted Sobolev spaces are given. A Laguerre–Legendre spectral scheme for the stream function form of Navier–Stokes equations is constructed. The stability and the convergence of the proposed scheme are proved. The numerical experiments show the high accuracy of this method. The main techniques used in this paper are also applicable to other nonlinear partial differential equations in an infinite strip.     

Some properties of Lagrangian sequences

Here we discuss a sequence of Lagrangians and corresponding Euler–Lagrange equations and point out some interesting properties that this particular sequence holds. This is an extension on recent results obtained on sequences of differential equations to Lagrangians as it leads to differential equations on the application of the Variational Principle.     

Jacobi interpolation approximations and their applications to singular differential equations

Jacobi–Gauss-type interpolations are considered. Some approximation results in certain Hilbert spaces are established. They are used for numerical solutions of singular differential equations and other related problems. The numerical results are illustrated.     

On the approximation of solutions of stochastic differential inclusions

A sequence of approximating equations is constructed for stochastic differential inclusions, and the properties of the measures corresponding to solutions of the approximating equations are studied for the class of stochastic differential inclusions.Translated fromTeoriya Sluchaínykh Protsessov, Vol. 14, pp. 43–48, 1986.     

Anticipating stochastic differential systems with memory

This article establishes existence and uniqueness of solutions to two classes of stochastic systems with finite memory subject to anticipating initial conditions which are sufficiently smooth in the Malliavin sense. The two classes are semilinear stochastic functional differential equations (sfdes) and fully nonlinear sfdes with a sublinear drift term. For the semilinear case, we use Malliavin calculus techniques, existence of the stochastic semiflow and an infinite-dimensional substitution theorem. For the fully nonlinear case, we employ an anticipating version of the Itô–Ventzell formula due to Ocone and Pardoux [D. Ocone, E. Pardoux, A generalized Itô–Ventzell formula. Application to a class of anticipating stochastic differential equations, Annales de l’Institut Henri Poincaré. Probabilité s et Statistiques 25 (1) (1989) 39–71]. In both cases, the use of Malliavin calculus techniques is necessitated by the infinite dimensionality of the initial condition.     

Explicit Formulas for Generalized Action–Angle Variables in a Neighborhood of an Isotropic Torus and Their Application

Different versions of the Darboux–Weinstein theorem guarantee the existence of action–angle-type variables and the harmonic-oscillator variables in a neighborhood of isotropic tori in the phase space. The procedure for constructing these variables is reduced to solving a rather complicated system of partial differential equations. We show that this system can be integrated in quadratures, which permits reducing the problem of constructing these variables to solving a system of quadratic equations. We discuss several applications of this purely geometric fact in problems of classical and quantum mechanics.     

Partitioned Runge–Kutta Methods in Lie-Group Setting

We introduce partitioned Runge–Kutta (PRK) methods as geometric integrators in the Runge–Kutta–Munthe-Kaas (RKMK) method hierarchy. This is done by first noticing that tangent and cotangent bundles are the natural domains for the differential equations to be solved. Next, we equip the (co)tangent bundle of a Lie group with a group structure and treat it as a Lie group. The structure of the differential equations on the (co)tangent-bundle Lie group is such that partitioned versions of the RKMK methods are naturally introduced. Numerical examples are included to illustrate the new methods.     

Existence and global bounds for a fluid model of plasma display technology

This article considers the fluid model for the discharge of plasma particle species in display technology. The fluid equations are coupled with Poisson's equation, which describes the effect of the charged particles on the electric field. The diffusion and mobility coefficients for the positive ion particles depend on the electric field, while those for the electrons depend on the electron mean energy. The reaction rates are proportional to the products of the densities of the reacting particles involved in the particular ionization, conversion or recombination reactions. Moreover, the ionization coefficients are dependent on the electric field, which varies spatially and temporally. The main ionization and discharge reactions are described by an initial-boundary value problem for a system of coupled parabolic–elliptic partial differential equations. The system is first analyzed by upper–lower solution method. By means of the a priori bounds obtained for an arbitrary time, the existence of solution for the initial-boundary value problem is proved in an appropriate Hölder space.