Wong-Zakai method for stochastic differential equations in engineering

In this paper, Wong-Zakai approximation methods are presented for some stochastic differential equations in engineering sciences. Wong-Zakai approximate solutions of the equations are analyzed and the numerical results are compared with results from popular approximation schemes for stochastic differential equations such as Euler-Alartiyama and Milstein methods. Several differential equations from engineering problems containing stochastic noise are investigated as numerical examples. Results show that Wong-Zakai method is a reliable tool for studying stochastic differential equations and can he used as an alternative for the known approximation techniques for stochastic models

An intelligent processing environment for real-time simulation

The development of a highly efficient and thus truly intelligent processing environment for real-time general purpose simulation of continuous systems is described. Such an environment can be created by mapping the simulation process directly onto the University of Alamba's OPERA architecture. To facilitate this effort, the field of continuous simulation is explored, highlighting areas in which efficiency can be improved. Areas in which parallel processing can be applied are also identified, and several general OPERA type hardware configurations that support improved simulation are investigated. Three direct execution parallel processing environments are introduced, each of which greatly improves efficiency by exploiting distinct areas of the simulation process. These suggested environments are candidate architectures around which a highly intelligent real-time simulation configuration can be developed

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit

Implicit-explicit predictor-corrector methods combined with improved spectral methods for pricing European style vanilla and exotic options

In this paper we present a robust numerical method to solve several types of European style option pricing problems. The governing equations are described by variants of Black-Scholes partial differential equations (BS-PDEs) of the reaction-diffusion-advection type. To discretise these BS-PDEs numerically, we use the spectral methods in the asset (spatial) direction and couple them with a third-order implicit-explicit predictor-corrector (IMEX-PC) method for the discretisation in the time direction. The use of this high-order time integration scheme sustains the better accuracy of the spectral methods for which they are well-known. Our spectral method consists of a pseudospectral formulation of the BS-PDEs by means of an improved Lagrange formula. On the other hand, in the IMEX-PC methods, we integrate the diffusion terms implicitly whereas the reaction and advection terms are integrated explicitly. Using this combined approach, we first solve the equations for standard European options and then extend this approach to digital options, butterfly spread options, and European calls in the Heston model. Numerical experiments illustrate that our approach is highly accurate and very efficient for pricing financial options such as those described above