High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (meansquare stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

Stability of numerical method for semi-linear stochastic pantograph differential equations

Directory of Open Access Journals (Sweden)

Yu Zhang

2016-01-01

Full Text Available Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results.

Stochastic fractional differential equations: Modeling, method and analysis

International Nuclear Information System (INIS)

Pedjeu, Jean-C.; Ladde, Gangaram S.

2012-01-01

By introducing a concept of dynamic process operating under multi-time scales in sciences and engineering, a mathematical model described by a system of multi-time scale stochastic differential equations is formulated. The classical Picard–Lindelöf successive approximations scheme is applied to the model validation problem, namely, existence and uniqueness of solution process. Naturally, this leads to the problem of finding closed form solutions of both linear and nonlinear multi-time scale stochastic differential equations of Itô–Doob type. Finally, to illustrate the scope of ideas and presented results, multi-time scale stochastic models for ecological and epidemiological processes in population dynamic are outlined.

Numerical methods for stochastic partial differential equations with white noise

CERN Document Server

Zhang, Zhongqiang

2017-01-01

This book covers numerical methods for stochastic partial differential equations with white noise using the framework of Wong-Zakai approximation. The book begins with some motivational and background material in the introductory chapters and is divided into three parts. Part I covers numerical stochastic ordinary differential equations. Here the authors start with numerical methods for SDEs with delay using the Wong-Zakai approximation and finite difference in time. Part II covers temporal white noise. Here the authors consider SPDEs as PDEs driven by white noise, where discretization of white noise (Brownian motion) leads to PDEs with smooth noise, which can then be treated by numerical methods for PDEs. In this part, recursive algorithms based on Wiener chaos expansion and stochastic collocation methods are presented for linear stochastic advection-diffusion-reaction equations. In addition, stochastic Euler equations are exploited as an application of stochastic collocation methods, where a numerical compa...

Exponential p-stability of impulsive stochastic differential equations with delays

International Nuclear Information System (INIS)

Yang Zhiguo; Xu Daoyi; Xiang Li

2006-01-01

In this Letter, we establish a method to study the exponential p-stability of the zero solution of impulsive stochastic differential equations with delays. By establishing an L-operator inequality and using the properties of M-cone and stochastic analysis technique, we obtain some new conditions ensuring the exponential p-stability of the zero solution of impulsive stochastic differential equations with delays. Two illustrative examples have been provided to show the effectiveness of our results

Numerical Analysis for Stochastic Partial Differential Delay Equations with Jumps

OpenAIRE

Li, Yan; Hu, Junhao

2013-01-01

We investigate the convergence rate of Euler-Maruyama method for a class of stochastic partial differential delay equations driven by both Brownian motion and Poisson point processes. We discretize in space by a Galerkin method and in time by using a stochastic exponential integrator. We generalize some results of Bao et al. (2011) and Jacob et al. (2009) in finite dimensions to a class of stochastic partial differential delay equations with jumps in infinite dimensions.

Stochastic differential equations and diffusion processes

CERN Document Server

Ikeda, N

1989-01-01

Being a systematic treatment of the modern theory of stochastic integrals and stochastic differential equations, the theory is developed within the martingale framework, which was developed by J.L. Doob and which plays an indispensable role in the modern theory of stochastic analysis.A considerable number of corrections and improvements have been made for the second edition of this classic work. In particular, major and substantial changes are in Chapter III and Chapter V where the sections treating excursions of Brownian Motion and the Malliavin Calculus have been expanded and refined. Sectio

Parameter estimation in stochastic differential equations

CERN Document Server

Bishwal, Jaya P N

2008-01-01

Parameter estimation in stochastic differential equations and stochastic partial differential equations is the science, art and technology of modelling complex phenomena and making beautiful decisions. The subject has attracted researchers from several areas of mathematics and other related fields like economics and finance. This volume presents the estimation of the unknown parameters in the corresponding continuous models based on continuous and discrete observations and examines extensively maximum likelihood, minimum contrast and Bayesian methods. Useful because of the current availability of high frequency data is the study of refined asymptotic properties of several estimators when the observation time length is large and the observation time interval is small. Also space time white noise driven models, useful for spatial data, and more sophisticated non-Markovian and non-semimartingale models like fractional diffusions that model the long memory phenomena are examined in this volume.

Asymptotic problems for stochastic partial differential equations

Science.gov (United States)

Salins, Michael

Stochastic partial differential equations (SPDEs) can be used to model systems in a wide variety of fields including physics, chemistry, and engineering. The main SPDEs of interest in this dissertation are the semilinear stochastic wave equations which model the movement of a material with constant mass density that is exposed to both determinstic and random forcing. Cerrai and Freidlin have shown that on fixed time intervals, as the mass density of the material approaches zero, the solutions of the stochastic wave equation converge uniformly to the solutions of a stochastic heat equation, in probability. This is called the Smoluchowski-Kramers approximation. In Chapter 2, we investigate some of the multi-scale behaviors that these wave equations exhibit. In particular, we show that the Freidlin-Wentzell exit place and exit time asymptotics for the stochastic wave equation in the small noise regime can be approximated by the exit place and exit time asymptotics for the stochastic heat equation. We prove that the exit time and exit place asymptotics are characterized by quantities called quasipotentials and we prove that the quasipotentials converge. We then investigate the special case where the equation has a gradient structure and show that we can explicitly solve for the quasipotentials, and that the quasipotentials for the heat equation and wave equation are equal. In Chapter 3, we study the Smoluchowski-Kramers approximation in the case where the material is electrically charged and exposed to a magnetic field. Interestingly, if the system is frictionless, then the Smoluchowski-Kramers approximation does not hold. We prove that the Smoluchowski-Kramers approximation is valid for systems exposed to both a magnetic field and friction. Notably, we prove that the solutions to the second-order equations converge to the solutions of the first-order equation in an Lp sense. This strengthens previous results where convergence was proved in probability.

Stochastic differential equations and a biological system

DEFF Research Database (Denmark)

Wang, Chunyan

1994-01-01

The purpose of this Ph.D. study is to explore the property of a growth process. The study includes solving and simulating of the growth process which is described in terms of stochastic differential equations. The identification of the growth and variability parameters of the process based...... on experimental data is considered. As an example, the growth of bacteria Pseudomonas fluorescens is taken. Due to the specific features of stochastic differential equations, namely that their solutions do not exist in the general sense, two new integrals - the Ito integral and the Stratonovich integral - have...... description. In order to identify the parameters, a Maximum likelihood estimation method is used together with a simplified truncated second order filter. Because of the continuity feature of the predictor equation, two numerical integration methods, called the Odeint and the Discretization method...

Modeling animal movements using stochastic differential equations

Science.gov (United States)

Haiganoush K. Preisler; Alan A. Ager; Bruce K. Johnson; John G. Kie

2004-01-01

We describe the use of bivariate stochastic differential equations (SDE) for modeling movements of 216 radiocollared female Rocky Mountain elk at the Starkey Experimental Forest and Range in northeastern Oregon. Spatially and temporally explicit vector fields were estimated using approximating difference equations and nonparametric regression techniques. Estimated...

Ambit processes and stochastic partial differential equations

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole; Benth, Fred Espen; Veraart, Almut

Ambit processes are general stochastic processes based on stochastic integrals with respect to Lévy bases. Due to their flexible structure, they have great potential for providing realistic models for various applications such as in turbulence and finance. This papers studies the connection betwe...... ambit processes and solutions to stochastic partial differential equations. We investigate this relationship from two angles: from the Walsh theory of martingale measures and from the viewpoint of the Lévy noise analysis....

Symmetries of th-Order Approximate Stochastic Ordinary Differential Equations

OpenAIRE

Fredericks, E.; Mahomed, F. M.

2012-01-01

Symmetries of $n$ th-order approximate stochastic ordinary differential equations (SODEs) are studied. The determining equations of these SODEs are derived in an Itô calculus context. These determining equations are not stochastic in nature. SODEs are normally used to model nature (e.g., earthquakes) or for testing the safety and reliability of models in construction engineering when looking at the impact of random perturbations.

Modelling Evolutionary Algorithms with Stochastic Differential Equations.

Science.gov (United States)

Heredia, Jorge Pérez

2017-11-20

There has been renewed interest in modelling the behaviour of evolutionary algorithms (EAs) by more traditional mathematical objects, such as ordinary differential equations or Markov chains. The advantage is that the analysis becomes greatly facilitated due to the existence of well established methods. However, this typically comes at the cost of disregarding information about the process. Here, we introduce the use of stochastic differential equations (SDEs) for the study of EAs. SDEs can produce simple analytical results for the dynamics of stochastic processes, unlike Markov chains which can produce rigorous but unwieldy expressions about the dynamics. On the other hand, unlike ordinary differential equations (ODEs), they do not discard information about the stochasticity of the process. We show that these are especially suitable for the analysis of fixed budget scenarios and present analogues of the additive and multiplicative drift theorems from runtime analysis. In addition, we derive a new more general multiplicative drift theorem that also covers non-elitist EAs. This theorem simultaneously allows for positive and negative results, providing information on the algorithm's progress even when the problem cannot be optimised efficiently. Finally, we provide results for some well-known heuristics namely Random Walk (RW), Random Local Search (RLS), the (1+1) EA, the Metropolis Algorithm (MA), and the Strong Selection Weak Mutation (SSWM) algorithm.

Probabilistic Forecasts of Solar Irradiance by Stochastic Differential Equations

DEFF Research Database (Denmark)

Iversen, Jan Emil Banning; Morales González, Juan Miguel; Møller, Jan Kloppenborg

2014-01-01

approach allows for characterizing both the interdependence structure of prediction errors of short-term solar irradiance and their predictive distribution. Three different stochastic differential equation models are first fitted to a training data set and subsequently evaluated on a one-year test set...... included in probabilistic forecasts may be paramount for decision makers to efficiently make use of this uncertain and variable generation. In this paper, a stochastic differential equation framework for modeling the uncertainty associated with the solar irradiance point forecast is proposed. This modeling...

STABILITY OF SOME KIND OF STOCHASTIC DIFFERENTIAL EQUATION

Institute of Scientific and Technical Information of China (English)

无

2011-01-01

In this paper,a kind of stochastic differential equation is investigated and the almost sure exponential stability of the equation is obtained using Gronwall's inequality.Further,we also give other noise intensity function to keep the stability of the system.

Numerical solution of second-order stochastic differential equations with Gaussian random parameters

Directory of Open Access Journals (Sweden)

Rahman Farnoosh

2014-07-01

Full Text Available In this paper, we present the numerical solution of ordinary differential equations (or SDEs, from each orderespecially second-order with time-varying and Gaussian random coefficients. We indicate a complete analysisfor second-order equations in specially case of scalar linear second-order equations (damped harmonicoscillators with additive or multiplicative noises. Making stochastic differential equations system from thisequation, it could be approximated or solved numerically by different numerical methods. In the case oflinear stochastic differential equations system by Computing fundamental matrix of this system, it could becalculated based on the exact solution of this system. Finally, this stochastic equation is solved by numericallymethod like E.M. and Milstein. Also its Asymptotic stability and statistical concepts like expectationand variance of solutions are discussed.

Population stochastic modelling (PSM)-An R package for mixed-effects models based on stochastic differential equations

DEFF Research Database (Denmark)

Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode

2009-01-01

are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE1 approximation......The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model...... development, J. Pharmacokinet. Pharmacodyn. 32 (February(l)) (2005) 109-141; C.W. Tornoe, R.V Overgaard, H. Agerso, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8...

PC analysis of stochastic differential equations driven by Wiener noise

KAUST Repository

Le Maitre, Olivier

2015-03-01

A polynomial chaos (PC) analysis with stochastic expansion coefficients is proposed for stochastic differential equations driven by additive or multiplicative Wiener noise. It is shown that for this setting, a Galerkin formalism naturally leads to the definition of a hierarchy of stochastic differential equations governing the evolution of the PC modes. Under the mild assumption that the Wiener and uncertain parameters can be treated as independent random variables, it is also shown that the Galerkin formalism naturally separates parametric uncertainty and stochastic forcing dependences. This enables us to perform an orthogonal decomposition of the process variance, and consequently identify contributions arising from the uncertainty in parameters, the stochastic forcing, and a coupled term. Insight gained from this decomposition is illustrated in light of implementation to simplified linear and non-linear problems; the case of a stochastic bifurcation is also considered.

Numerical analysis of systems of ordinary and stochastic differential equations

CERN Document Server

Artemiev, S S

1997-01-01

This text deals with numerical analysis of systems of both ordinary and stochastic differential equations. It covers numerical solution problems of the Cauchy problem for stiff ordinary differential equations (ODE) systems by Rosenbrock-type methods (RTMs).

Population stochastic modelling (PSM)--an R package for mixed-effects models based on stochastic differential equations.

Science.gov (United States)

Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode; Overgaard, Rune Viig; Madsen, Henrik

2009-06-01

The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model development, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 109-141; C.W. Tornøe, R.V. Overgaard, H. Agersø, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8)) (2005) 1247-1258; R.V. Overgaard, N. Jonsson, C.W. Tornøe, H. Madsen, Non-linear mixed-effects models with stochastic differential equations: implementation of an estimation algorithm, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 85-107; U. Picchini, S. Ditlevsen, A. De Gaetano, Maximum likelihood estimation of a time-inhomogeneous stochastic differential model of glucose dynamics, Math. Med. Biol. 25 (June(2)) (2008) 141-155]. PK/PD models are traditionally based ordinary differential equations (ODEs) with an observation link that incorporates noise. This state-space formulation only allows for observation noise and not for system noise. Extending to SDEs allows for a Wiener noise component in the system equations. This additional noise component enables handling of autocorrelated residuals originating from natural variation or systematic model error. Autocorrelated residuals are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE(1) approximation to the population likelihood which is generated from the individual likelihoods that are approximated using the Extended Kalman Filter's one-step predictions.

On stochastic differential equations with random delay

International Nuclear Information System (INIS)

Krapivsky, P L; Luck, J M; Mallick, K

2011-01-01

We consider stochastic dynamical systems defined by differential equations with a uniform random time delay. The latter equations are shown to be equivalent to deterministic higher-order differential equations: for an nth-order equation with random delay, the corresponding deterministic equation has order n + 1. We analyze various examples of dynamical systems of this kind, and find a number of unusual behaviors. For instance, for the harmonic oscillator with random delay, the energy grows as exp((3/2) t 2/3 ) in reduced units. We then investigate the effect of introducing a discrete time step ε. At variance with the continuous situation, the discrete random recursion relations thus obtained have intrinsic fluctuations. The crossover between the fluctuating discrete problem and the deterministic continuous one as ε goes to zero is studied in detail on the example of a first-order linear differential equation

Approximate Controllability for Linear Stochastic Differential Equations in Infinite Dimensions

International Nuclear Information System (INIS)

Goreac, D.

2009-01-01

The objective of the paper is to investigate the approximate controllability property of a linear stochastic control system with values in a separable real Hilbert space. In a first step we prove the existence and uniqueness for the solution of the dual linear backward stochastic differential equation. This equation has the particularity that in addition to an unbounded operator acting on the Y-component of the solution there is still another one acting on the Z-component. With the help of this dual equation we then deduce the duality between approximate controllability and observability. Finally, under the assumption that the unbounded operator acting on the state process of the forward equation is an infinitesimal generator of an exponentially stable semigroup, we show that the generalized Hautus test provides a necessary condition for the approximate controllability. The paper generalizes former results by Buckdahn, Quincampoix and Tessitore (Stochastic Partial Differential Equations and Applications, Series of Lecture Notes in Pure and Appl. Math., vol. 245, pp. 253-260, Chapman and Hall, London, 2006) and Goreac (Applied Analysis and Differential Equations, pp. 153-164, World Scientific, Singapore, 2007) from the finite dimensional to the infinite dimensional case

Short-term Probabilistic Forecasting of Wind Speed Using Stochastic Differential Equations

DEFF Research Database (Denmark)

Iversen, Jan Emil Banning; Morales González, Juan Miguel; Møller, Jan Kloppenborg

2016-01-01

and uncertain nature. In this paper, we propose a modeling framework for wind speed that is based on stochastic differential equations. We show that stochastic differential equations allow us to naturally capture the time dependence structure of wind speed prediction errors (from 1 up to 24 hours ahead) and......It is widely accepted today that probabilistic forecasts of wind power production constitute valuable information for both wind power producers and power system operators to economically exploit this form of renewable energy, while mitigating the potential adverse effects related to its variable......, most importantly, to derive point and quantile forecasts, predictive distributions, and time-path trajectories (also referred to as scenarios or ensemble forecasts), all by one single stochastic differential equation model characterized by a few parameters....

Introduction to stochastic analysis integrals and differential equations

CERN Document Server

Mackevicius, Vigirdas

2013-01-01

This is an introduction to stochastic integration and stochastic differential equations written in an understandable way for a wide audience, from students of mathematics to practitioners in biology, chemistry, physics, and finances. The presentation is based on the naïve stochastic integration, rather than on abstract theories of measure and stochastic processes. The proofs are rather simple for practitioners and, at the same time, rather rigorous for mathematicians. Detailed application examples in natural sciences and finance are presented. Much attention is paid to simulation diffusion pro

Efficient Estimating Functions for Stochastic Differential Equations

DEFF Research Database (Denmark)

Jakobsen, Nina Munkholt

The overall topic of this thesis is approximate martingale estimating function-based estimationfor solutions of stochastic differential equations, sampled at high frequency. Focuslies on the asymptotic properties of the estimators. The first part of the thesis deals with diffusions observed over...

COMPARISON THEOREM OF BACKWARD DOUBLY STOCHASTIC DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

This paper is devoted to deriving a comparison theorem of solutions to backward doubly stochastic differential equations driven by Brownian motion and backward It-Kunita integral. By the application of this theorem, we give an existence result of the solutions to these equations with continuous coefficients.

Accelerated Genetic Algorithm Solutions Of Some Parametric Families Of Stochastic Differential Equations

Directory of Open Access Journals (Sweden)

Eman Ali Hussain

2015-01-01

Full Text Available Absract In this project A new method for solving Stochastic Differential Equations SDEs deriving by Wiener process numerically will be construct and implement using Accelerated Genetic Algorithm AGA. An SDE is a differential equation in which one or more of the terms and hence the solutions itself is a stochastic process. Solving stochastic differential equations requires going away from the recognizable deterministic setting of ordinary and partial differential equations into a world where the evolution of a quantity has an inherent random component and where the expected behavior of this quantity can be described in terms of probability distributions. We applied our method on the Ito formula which is equivalent to the SDE to find approximation solution of the SDEs. Numerical experiments illustrate the behavior of the proposed method.

Modeling and Prediction Using Stochastic Differential Equations

DEFF Research Database (Denmark)

Juhl, Rune; Møller, Jan Kloppenborg; Jørgensen, John Bagterp

2016-01-01

Pharmacokinetic/pharmakodynamic (PK/PD) modeling for a single subject is most often performed using nonlinear models based on deterministic ordinary differential equations (ODEs), and the variation between subjects in a population of subjects is described using a population (mixed effects) setup...... deterministic and can predict the future perfectly. A more realistic approach would be to allow for randomness in the model due to e.g., the model be too simple or errors in input. We describe a modeling and prediction setup which better reflects reality and suggests stochastic differential equations (SDEs...

Stochastic differential equations used to model conjugation

DEFF Research Database (Denmark)

Philipsen, Kirsten Riber; Christiansen, Lasse Engbo

Stochastic differential equations (SDEs) are used to model horizontal transfer of antibiotic resis- tance by conjugation. The model describes the concentration of donor, recipient, transconjugants and substrate. The strength of the SDE model over the traditional ODE models is that the noise can...

General Large Deviations and Functional Iterated Logarithm Law for Multivalued Stochastic Differential Equations

OpenAIRE

Ren, Jiagang; Wu, Jing; Zhang, Hua

2015-01-01

In this paper, we prove a large deviation principle of Freidlin-Wentzell's type for the multivalued stochastic differential equations. As an application, we derive a functional iterated logarithm law for the solutions of multivalued stochastic differential equations.

PC analysis of stochastic differential equations driven by Wiener noise

KAUST Repository

Le Maitre, Olivier; Knio, Omar

2015-01-01

A polynomial chaos (PC) analysis with stochastic expansion coefficients is proposed for stochastic differential equations driven by additive or multiplicative Wiener noise. It is shown that for this setting, a Galerkin formalism naturally leads

Noise Analysis of Single-Ended Input Differential Amplifier using Stochastic Differential Equation

OpenAIRE

Tarun Kumar Rawat; Abhirup Lahiri; Ashish Gupta

2008-01-01

In this paper, we analyze the effect of noise in a single- ended input differential amplifier working at high frequencies. Both extrinsic and intrinsic noise are analyzed using time domain method employing techniques from stochastic calculus. Stochastic differential equations are used to obtain autocorrelation functions of the output noise voltage and other solution statistics like mean and variance. The analysis leads to important design implications and suggests changes in the device parame...

Causal interpretation of stochastic differential equations

DEFF Research Database (Denmark)

Sokol, Alexander; Hansen, Niels Richard

2014-01-01

We give a causal interpretation of stochastic differential equations (SDEs) by defining the postintervention SDE resulting from an intervention in an SDE. We show that under Lipschitz conditions, the solution to the postintervention SDE is equal to a uniform limit in probability of postintervention...... structural equation models based on the Euler scheme of the original SDE, thus relating our definition to mainstream causal concepts. We prove that when the driving noise in the SDE is a Lévy process, the postintervention distribution is identifiable from the generator of the SDE....

2–stage stochastic Runge–Kutta for stochastic delay differential equations

Energy Technology Data Exchange (ETDEWEB)

Rosli, Norhayati; Jusoh Awang, Rahimah [Faculty of Industrial Science and Technology, Universiti Malaysia Pahang, Lebuhraya Tun Razak, 26300, Gambang, Pahang (Malaysia); Bahar, Arifah; Yeak, S. H. [Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor (Malaysia)

2015-05-15

This paper proposes a newly developed one-step derivative-free method, that is 2-stage stochastic Runge-Kutta (SRK2) to approximate the solution of stochastic delay differential equations (SDDEs) with a constant time lag, r > 0. General formulation of stochastic Runge-Kutta for SDDEs is introduced and Stratonovich Taylor series expansion for numerical solution of SRK2 is presented. Local truncation error of SRK2 is measured by comparing the Stratonovich Taylor expansion of the exact solution with the computed solution. Numerical experiment is performed to assure the validity of the method in simulating the strong solution of SDDEs.

A concise course on stochastic partial differential equations

CERN Document Server

Prévôt, Claudia

2007-01-01

These lectures concentrate on (nonlinear) stochastic partial differential equations (SPDE) of evolutionary type. All kinds of dynamics with stochastic influence in nature or man-made complex systems can be modelled by such equations. To keep the technicalities minimal we confine ourselves to the case where the noise term is given by a stochastic integral w.r.t. a cylindrical Wiener process.But all results can be easily generalized to SPDE with more general noises such as, for instance, stochastic integral w.r.t. a continuous local martingale. There are basically three approaches to analyze SPDE: the "martingale measure approach", the "mild solution approach" and the "variational approach". The purpose of these notes is to give a concise and as self-contained as possible an introduction to the "variational approach". A large part of necessary background material, such as definitions and results from the theory of Hilbert spaces, are included in appendices.

Improved stochastic approximation methods for discretized parabolic partial differential equations

Science.gov (United States)

Guiaş, Flavius

2016-12-01

We present improvements of the stochastic direct simulation method, a known numerical scheme based on Markov jump processes which is used for approximating solutions of ordinary differential equations. This scheme is suited especially for spatial discretizations of evolution partial differential equations (PDEs). By exploiting the full path simulation of the stochastic method, we use this first approximation as a predictor and construct improved approximations by Picard iterations, Runge-Kutta steps, or a combination. This has as consequence an increased order of convergence. We illustrate the features of the improved method at a standard benchmark problem, a reaction-diffusion equation modeling a combustion process in one space dimension (1D) and two space dimensions (2D).

STABILITY OF STOCHASTIC DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper,we obtain suffcient conditions for the stability in p-th moment of the analytical solutions and the mean square stability of a stochastic differential equation with unbounded delay proposed in [6,10] using the explicit Euler method.

Stochastic Differential Equation-Based Flexible Software Reliability Growth Model

Directory of Open Access Journals (Sweden)

P. K. Kapur

2009-01-01

Full Text Available Several software reliability growth models (SRGMs have been developed by software developers in tracking and measuring the growth of reliability. As the size of software system is large and the number of faults detected during the testing phase becomes large, so the change of the number of faults that are detected and removed through each debugging becomes sufficiently small compared with the initial fault content at the beginning of the testing phase. In such a situation, we can model the software fault detection process as a stochastic process with continuous state space. In this paper, we propose a new software reliability growth model based on Itô type of stochastic differential equation. We consider an SDE-based generalized Erlang model with logistic error detection function. The model is estimated and validated on real-life data sets cited in literature to show its flexibility. The proposed model integrated with the concept of stochastic differential equation performs comparatively better than the existing NHPP-based models.

On the definition of an admitted Lie group for stochastic differential equations with multi-Brownian motion

International Nuclear Information System (INIS)

Srihirun, B; Meleshko, S V; Schulz, E

2006-01-01

The definition of an admitted Lie group of transformations for stochastic differential equations has been already presented for equations with one-dimensional Brownian motion. The transformation of the dependent variables involves time as well, and it has been proven that Brownian motion is transformed to Brownian motion. In this paper, we will discuss this concept for stochastic differential equations involving multi-dimensional Brownian motion and present applications to a variety of stochastic differential equations

Modelling conjugation with stochastic differential equations

DEFF Research Database (Denmark)

Philipsen, Kirsten Riber; Christiansen, Lasse Engbo; Hasman, Henrik

2010-01-01

Enterococcus faecium strains in a rich exhaustible media. The model contains a new expression for a substrate dependent conjugation rate. A maximum likelihood based method is used to estimate the model parameters. Different models including different noise structure for the system and observations are compared......Conjugation is an important mechanism involved in the transfer of resistance between bacteria. In this article a stochastic differential equation based model consisting of a continuous time state equation and a discrete time measurement equation is introduced to model growth and conjugation of two...... using a likelihood-ratio test and Akaike's information criterion. Experiments indicating conjugation on the agar plates selecting for transconjugants motivates the introduction of an extended model, for which conjugation on the agar plate is described in the measurement equation. This model is compared...

Reflected backward stochastic differential equations in an orthant

Indian Academy of Sciences (India)

R. Narasimhan (Krishtel eMaging) 1461 1996 Oct 15 13:05:22

Since backward stochastic differential equations were introduced about a decade back there has been a lot .... is given in ([10] pp. 76–80) in connection with financial networks; these authors call ..... Applying Theorem 3.2 to (̂Y. (k) i. , ̂Z. (k).

A Simple Stochastic Differential Equation with Discontinuous Drift

DEFF Research Database (Denmark)

Simonsen, Maria; Leth, John-Josef; Schiøler, Henrik

2013-01-01

In this paper we study solutions to stochastic differential equations (SDEs) with discontinuous drift. We apply two approaches: The Euler-Maruyama method and the Fokker-Planck equation and show that a candidate density function based on the Euler-Maruyama method approximates a candidate density...... function based on the stationary Fokker-Planck equation. Furthermore, we introduce a smooth function which approximates the discontinuous drift and apply the Euler-Maruyama method and the Fokker-Planck equation with this input. The point of departure for this work is a particular SDE with discontinuous...

Stochastic integration and differential equations

CERN Document Server

Protter, Philip E

2003-01-01

It has been 15 years since the first edition of Stochastic Integration and Differential Equations, A New Approach appeared, and in those years many other texts on the same subject have been published, often with connections to applications, especially mathematical finance. Yet in spite of the apparent simplicity of approach, none of these books has used the functional analytic method of presenting semimartingales and stochastic integration. Thus a 2nd edition seems worthwhile and timely, though it is no longer appropriate to call it "a new approach". The new edition has several significant changes, most prominently the addition of exercises for solution. These are intended to supplement the text, but lemmas needed in a proof are never relegated to the exercises. Many of the exercises have been tested by graduate students at Purdue and Cornell Universities. Chapter 3 has been completely redone, with a new, more intuitive and simultaneously elementary proof of the fundamental Doob-Meyer decomposition theorem, t...

Attempts at a numerical realisation of stochastic differential equations containing Preisach operator

International Nuclear Information System (INIS)

McCarthy, S; Rachinskii, D

2011-01-01

We describe two Euler type numerical schemes obtained by discretisation of a stochastic differential equation which contains the Preisach memory operator. Equations of this type are of interest in areas such as macroeconomics and terrestrial hydrology where deterministic models containing the Preisach operator have been developed but do not fully encapsulate stochastic aspects of the area. A simple price dynamics model is presented as one motivating example for our studies. Some numerical evidence is given that the two numerical schemes converge to the same limit as the time step decreases. We show that the Preisach term introduces a damping effect which increases on the parts of the trajectory demonstrating a stronger upwards or downwards trend. The results are preliminary to a broader programme of research of stochastic differential equations with the Preisach hysteresis operator.

Linear stochastic differential equations with anticipating initial conditions

DEFF Research Database (Denmark)

Khalifa, Narjess; Kuo, Hui-Hsiung; Ouerdiane, Habib

In this paper we use the new stochastic integral introduced by Ayed and Kuo (2008) and the results obtained by Kuo et al. (2012b) to find a solution to a drift-free linear stochastic differential equation with anticipating initial condition. Our solution is based on well-known results from...... classical Itô theory and anticipative Itô formula results from Kue et al. (2012b). We also show that the solution obtained by our method is consistent with the solution obtained by the methods of Malliavin calculus, e.g. Buckdahn and Nualart (1994)....

EXISTENCE AND UNIQUENESS OF SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATION WITH RANDOM COEFFICIENTS

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

This paper mainly deals with a stochastic differential equation (SDE) with random coefficients. Sufficient conditions which guarantee the existence and uniqueness of solutions to the equation are given.

Stability analysis for neutral stochastic differential equation of second order driven by Poisson jumps

Science.gov (United States)

Chadha, Alka; Bora, Swaroop Nandan

2017-11-01

This paper studies the existence, uniqueness, and exponential stability in mean square for the mild solution of neutral second order stochastic partial differential equations with infinite delay and Poisson jumps. By utilizing the Banach fixed point theorem, first the existence and uniqueness of the mild solution of neutral second order stochastic differential equations is established. Then, the mean square exponential stability for the mild solution of the stochastic system with Poisson jumps is obtained with the help of an established integral inequality.

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer from the step size reduction faced by standard explicit methods. The family is based on the standard second order orthogonal Runge-Kutta-Chebyshev (ROCK2) methods for deterministic problems. The convergence, meansquare, and asymptotic stability properties of the methods are analyzed. Numerical experiments, including applications to nonlinear SDEs and parabolic stochastic partial differential equations are presented and confirm the theoretical results. © 2013 Society for Industrial and Applied Mathematics.

Numerical Solution of Stochastic Nonlinear Fractional Differential Equations

KAUST Repository

El-Beltagy, Mohamed A.

2015-01-07

Using Wiener-Hermite expansion (WHE) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. WHE is the only known expansion that handles the white/colored noise exactly. This work introduces a numerical estimation of the stochastic response of the Duffing oscillator with fractional or variable order damping and driven by white noise. The WHE technique is integrated with the Grunwald-Letnikov approximation in case of fractional order and with Coimbra approximation in case of variable-order damping. The numerical solver was tested with the analytic solution and with Monte-Carlo simulations. The developed mixed technique was shown to be efficient in simulating SPDEs.

Numerical Solution of Stochastic Nonlinear Fractional Differential Equations

KAUST Repository

El-Beltagy, Mohamed A.; Al-Juhani, Amnah

2015-01-01

Using Wiener-Hermite expansion (WHE) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. WHE is the only known expansion that handles the white/colored noise exactly. This work introduces a numerical estimation of the stochastic response of the Duffing oscillator with fractional or variable order damping and driven by white noise. The WHE technique is integrated with the Grunwald-Letnikov approximation in case of fractional order and with Coimbra approximation in case of variable-order damping. The numerical solver was tested with the analytic solution and with Monte-Carlo simulations. The developed mixed technique was shown to be efficient in simulating SPDEs.

Stability of the trivial solution for linear stochastic differential equations with Poisson white noise

International Nuclear Information System (INIS)

Grigoriu, Mircea; Samorodnitsky, Gennady

2004-01-01

Two methods are considered for assessing the asymptotic stability of the trivial solution of linear stochastic differential equations driven by Poisson white noise, interpreted as the formal derivative of a compound Poisson process. The first method attempts to extend a result for diffusion processes satisfying linear stochastic differential equations to the case of linear equations with Poisson white noise. The developments for the method are based on Ito's formula for semimartingales and Lyapunov exponents. The second method is based on a geometric ergodic theorem for Markov chains providing a criterion for the asymptotic stability of the solution of linear stochastic differential equations with Poisson white noise. Two examples are presented to illustrate the use and evaluate the potential of the two methods. The examples demonstrate limitations of the first method and the generality of the second method

Stochastic TDHF and the Boltzman-Langevin equation

International Nuclear Information System (INIS)

Suraud, E.; Reinhard, P.G.

1991-01-01

Outgoing from a time-dependent theory of correlations, we present a stochastic differential equation for the propagation of ensembles of Slater determinants, called Stochastic Time-Dependent Hartree-Fock (Stochastic TDHF). These ensembles are allowed to develop large fluctuations in the Hartree-Fock mean fields. An alternative stochastic differential equation, the Boltzmann-Langevin equation, can be derived from Stochastic TDHF by averaging over subensembles with small fluctuations

Stochastic differential equation model for linear growth birth and death processes with immigration and emigration

International Nuclear Information System (INIS)

Granita; Bahar, A.

2015-01-01

This paper discusses on linear birth and death with immigration and emigration (BIDE) process to stochastic differential equation (SDE) model. Forward Kolmogorov equation in continuous time Markov chain (CTMC) with a central-difference approximation was used to find Fokker-Planckequation corresponding to a diffusion process having the stochastic differential equation of BIDE process. The exact solution, mean and variance function of BIDE process was found

Stochastic differential equation model for linear growth birth and death processes with immigration and emigration

Energy Technology Data Exchange (ETDEWEB)

Granita, E-mail: granitafc@gmail.com [Dept. Mathematical Education, State Islamic University of Sultan Syarif Kasim Riau, 28293 Indonesia and Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310,Johor (Malaysia); Bahar, A. [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310,Johor Malaysia and UTM Center for Industrial and Applied Mathematics (UTM-CIAM) (Malaysia)

2015-03-09

This paper discusses on linear birth and death with immigration and emigration (BIDE) process to stochastic differential equation (SDE) model. Forward Kolmogorov equation in continuous time Markov chain (CTMC) with a central-difference approximation was used to find Fokker-Planckequation corresponding to a diffusion process having the stochastic differential equation of BIDE process. The exact solution, mean and variance function of BIDE process was found.

Application of Stochastic Partial Differential Equations to Reservoir Property Modelling

KAUST Repository

Potsepaev, R.

2010-09-06

Existing algorithms of geostatistics for stochastic modelling of reservoir parameters require a mapping (the \\'uvt-transform\\') into the parametric space and reconstruction of a stratigraphic co-ordinate system. The parametric space can be considered to represent a pre-deformed and pre-faulted depositional environment. Existing approximations of this mapping in many cases cause significant distortions to the correlation distances. In this work we propose a coordinate free approach for modelling stochastic textures through the application of stochastic partial differential equations. By avoiding the construction of a uvt-transform and stratigraphic coordinates, one can generate realizations directly in the physical space in the presence of deformations and faults. In particular the solution of the modified Helmholtz equation driven by Gaussian white noise is a zero mean Gaussian stationary random field with exponential correlation function (in 3-D). This equation can be used to generate realizations in parametric space. In order to sample in physical space we introduce a stochastic elliptic PDE with tensor coefficients, where the tensor is related to correlation anisotropy and its variation is physical space.

Simulation of quantum dynamics based on the quantum stochastic differential equation.

Science.gov (United States)

Li, Ming

2013-01-01

The quantum stochastic differential equation derived from the Lindblad form quantum master equation is investigated. The general formulation in terms of environment operators representing the quantum state diffusion is given. The numerical simulation algorithm of stochastic process of direct photodetection of a driven two-level system for the predictions of the dynamical behavior is proposed. The effectiveness and superiority of the algorithm are verified by the performance analysis of the accuracy and the computational cost in comparison with the classical Runge-Kutta algorithm.

Multivalued stochastic delay differential equations and related ...

African Journals Online (AJOL)

We study the existence and uniqueness of a solution for the multivalued stochastic differential equation with delay (the multivalued term is of subdifferential type):. dX(t) + aφ (X(t))dt ∍ b(t,X(t), Y(t), Z(t)) dt. ⎨ +σ (t, X (t), Y (t), Z (t)) dW (t), t ∈ (s, T). X(t) = ξ (t - s), t ∈ [s - δ, s]. Specify that in this case the coefficients at time t ...

Fractional Stochastic Differential Equations Satisfying Fluctuation-Dissipation Theorem

Science.gov (United States)

Li, Lei; Liu, Jian-Guo; Lu, Jianfeng

2017-10-01

We propose in this work a fractional stochastic differential equation (FSDE) model consistent with the over-damped limit of the generalized Langevin equation model. As a result of the `fluctuation-dissipation theorem', the differential equations driven by fractional Brownian noise to model memory effects should be paired with Caputo derivatives, and this FSDE model should be understood in an integral form. We establish the existence of strong solutions for such equations and discuss the ergodicity and convergence to Gibbs measure. In the linear forcing regime, we show rigorously the algebraic convergence to Gibbs measure when the `fluctuation-dissipation theorem' is satisfied, and this verifies that satisfying `fluctuation-dissipation theorem' indeed leads to the correct physical behavior. We further discuss possible approaches to analyze the ergodicity and convergence to Gibbs measure in the nonlinear forcing regime, while leave the rigorous analysis for future works. The FSDE model proposed is suitable for systems in contact with heat bath with power-law kernel and subdiffusion behaviors.

On the Existence and the Applications of Modified Equations for Stochastic Differential Equations

KAUST Repository

Zygalakis, K. C.

2011-01-01

In this paper we describe a general framework for deriving modified equations for stochastic differential equations (SDEs) with respect to weak convergence. Modified equations are derived for a variety of numerical methods, such as the Euler or the Milstein method. Existence of higher order modified equations is also discussed. In the case of linear SDEs, using the Gaussianity of the underlying solutions, we derive an SDE which the numerical method solves exactly in the weak sense. Applications of modified equations in the numerical study of Langevin equations is also discussed. © 2011 Society for Industrial and Applied Mathematics.

A class of degenerate stochastic differential equations with non ...

Indian Academy of Sciences (India)

Introduction. In this article we consider (possibly degenerate) stochastic differential equations (SDEs) with non-Lipschitz coefficients. If the coefficients are Lipschitz, we can prove the existence of a unique strong solution (see [9]). But uniqueness fails in the case of non-Lipschitz coefficients. The literature on this topic is not ...

Optimal Control for Stochastic Delay Evolution Equations

Energy Technology Data Exchange (ETDEWEB)

Meng, Qingxin, E-mail: mqx@hutc.zj.cn [Huzhou University, Department of Mathematical Sciences (China); Shen, Yang, E-mail: skyshen87@gmail.com [York University, Department of Mathematics and Statistics (Canada)

2016-08-15

In this paper, we investigate a class of infinite-dimensional optimal control problems, where the state equation is given by a stochastic delay evolution equation with random coefficients, and the corresponding adjoint equation is given by an anticipated backward stochastic evolution equation. We first prove the continuous dependence theorems for stochastic delay evolution equations and anticipated backward stochastic evolution equations, and show the existence and uniqueness of solutions to anticipated backward stochastic evolution equations. Then we establish necessary and sufficient conditions for optimality of the control problem in the form of Pontryagin’s maximum principles. To illustrate the theoretical results, we apply stochastic maximum principles to study two examples, an infinite-dimensional linear-quadratic control problem with delay and an optimal control of a Dirichlet problem for a stochastic partial differential equation with delay. Further applications of the two examples to a Cauchy problem for a controlled linear stochastic partial differential equation and an optimal harvesting problem are also considered.

Analytical Solutions for Multi-Time Scale Fractional Stochastic Differential Equations Driven by Fractional Brownian Motion and Their Applications

Directory of Open Access Journals (Sweden)

Xiao-Li Ding

2018-01-01

Full Text Available In this paper, we investigate analytical solutions of multi-time scale fractional stochastic differential equations driven by fractional Brownian motions. We firstly decompose homogeneous multi-time scale fractional stochastic differential equations driven by fractional Brownian motions into independent differential subequations, and give their analytical solutions. Then, we use the variation of constant parameters to obtain the solutions of nonhomogeneous multi-time scale fractional stochastic differential equations driven by fractional Brownian motions. Finally, we give three examples to demonstrate the applicability of our obtained results.

Stationary solutions of linear stochastic delay differential equations: applications to biological systems.

Science.gov (United States)

Frank, T D; Beek, P J

2001-08-01

Recently, Küchler and Mensch [Stochastics Stochastics Rep. 40, 23 (1992)] derived exact stationary probability densities for linear stochastic delay differential equations. This paper presents an alternative derivation of these solutions by means of the Fokker-Planck approach introduced by Guillouzic [Phys. Rev. E 59, 3970 (1999); 61, 4906 (2000)]. Applications of this approach, which is argued to have greater generality, are discussed in the context of stochastic models for population growth and tracking movements.

Some Additional Remarks on the Cumulant Expansion for Linear Stochastic Differential Equations

NARCIS (Netherlands)

Roerdink, J.B.T.M.

1984-01-01

We summarize our previous results on cumulant expansions for linear stochastic differential equations with correlated multipliclative and additive noise. The application of the general formulas to equations with statistically independent multiplicative and additive noise is reconsidered in detail,

Some additional remarks on the cumulant expansion for linear stochastic differential equations

NARCIS (Netherlands)

Roerdink, J.B.T.M.

1984-01-01

We summarize our previous results on cumular expasions for linear stochastic differential equations with correlated multipliclative and additive noise. The application of the general formulas to equations with statistically independent multiplicative and additive noise is reconsidered in detail,

Reflected stochastic differential equation models for constrained animal movement

Science.gov (United States)

Hanks, Ephraim M.; Johnson, Devin S.; Hooten, Mevin B.

2017-01-01

Movement for many animal species is constrained in space by barriers such as rivers, shorelines, or impassable cliffs. We develop an approach for modeling animal movement constrained in space by considering a class of constrained stochastic processes, reflected stochastic differential equations. Our approach generalizes existing methods for modeling unconstrained animal movement. We present methods for simulation and inference based on augmenting the constrained movement path with a latent unconstrained path and illustrate this augmentation with a simulation example and an analysis of telemetry data from a Steller sea lion (Eumatopias jubatus) in southeast Alaska.

Malliavin Calculus With Applications to Stochastic Partial Differential Equations

CERN Document Server

Sanz-Solé, Marta

2005-01-01

Developed in the 1970s to study the existence and smoothness of density for the probability laws of random vectors, Malliavin calculus--a stochastic calculus of variation on the Wiener space--has proven fruitful in many problems in probability theory, particularly in probabilistic numerical methods in financial mathematics.This book presents applications of Malliavin calculus to the analysis of probability laws of solutions to stochastic partial differential equations driven by Gaussian noises that are white in time and coloured in space. The first five chapters introduce the calculus itself

Green function of the double-fractional Fokker-Planck equation: Path integral and stochastic differential equations

Science.gov (United States)

Kleinert, H.; Zatloukal, V.

2013-11-01

The statistics of rare events, the so-called black-swan events, is governed by non-Gaussian distributions with heavy power-like tails. We calculate the Green functions of the associated Fokker-Planck equations and solve the related stochastic differential equations. We also discuss the subject in the framework of path integration.

A stochastic differential equation framework for the turbulent velocity field

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole Eiler; Schmiegel, Jürgen

We discuss a stochastic differential equation, as a modelling framework for the turbulent velocity field, that is capable of capturing basic stylized facts of the statistics of velocity increments. In particular, we focus on the evolution of the probability density of velocity increments...

Improving multilevel Monte Carlo for stochastic differential equations with application to the Langevin equation.

Science.gov (United States)

Müller, Eike H; Scheichl, Rob; Shardlow, Tony

2015-04-08

This paper applies several well-known tricks from the numerical treatment of deterministic differential equations to improve the efficiency of the multilevel Monte Carlo (MLMC) method for stochastic differential equations (SDEs) and especially the Langevin equation. We use modified equations analysis as an alternative to strong-approximation theory for the integrator, and we apply this to introduce MLMC for Langevin-type equations with integrators based on operator splitting. We combine this with extrapolation and investigate the use of discrete random variables in place of the Gaussian increments, which is a well-known technique for the weak approximation of SDEs. We show that, for small-noise problems, discrete random variables can lead to an increase in efficiency of almost two orders of magnitude for practical levels of accuracy.

Analytical Solutions for Multi-Time Scale Fractional Stochastic Differential Equations Driven by Fractional Brownian Motion and Their Applications

OpenAIRE

Xiao-Li Ding; Juan J. Nieto

2018-01-01

In this paper, we investigate analytical solutions of multi-time scale fractional stochastic differential equations driven by fractional Brownian motions. We firstly decompose homogeneous multi-time scale fractional stochastic differential equations driven by fractional Brownian motions into independent differential subequations, and give their analytical solutions. Then, we use the variation of constant parameters to obtain the solutions of nonhomogeneous multi-time scale fractional stochast...

Controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with delay and Poisson jumps

Directory of Open Access Journals (Sweden)

Diem Dang Huan

2015-12-01

Full Text Available The current paper is concerned with the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps in Hilbert spaces. Using the theory of a strongly continuous cosine family of bounded linear operators, stochastic analysis theory and with the help of the Banach fixed point theorem, we derive a new set of sufficient conditions for the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps. Finally, an application to the stochastic nonlinear wave equation with infinite delay and Poisson jumps is given.

Mild Solutions of Neutral Stochastic Partial Functional Differential Equations

Directory of Open Access Journals (Sweden)

T. E. Govindan

2011-01-01

Full Text Available This paper studies the existence and uniqueness of a mild solution for a neutral stochastic partial functional differential equation using a local Lipschitz condition. When the neutral term is zero and even in the deterministic special case, the result obtained here appears to be new. An example is included to illustrate the theory.

The Application of backward stochastic differential equation with stopping time in hedging American contingent claims

International Nuclear Information System (INIS)

Wang Bo; Song Ruili

2009-01-01

We consider a more general wealth process with a drift coefficient which is Lipschitz continuous and the portfolio process with convex constraint. We convert the problem of hedging American contingent claims into the problem of minimal solution of backward stochastic differential equation with stopping time. We adopt the penalization method for constructing the minimal solution of stochastic differential equations and obtain the upper hedging price of American contingent claims.

A stochastic differential equation framework for the timewise dynamics of turbulent velocities

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole Eiler; Schmiegel, Jürgen

2008-01-01

We discuss a stochastic differential equation as a modeling framework for the timewise dynamics of turbulent velocities. The equation is capable of capturing basic stylized facts of the statistics of temporal velocity increments. In particular, we focus on the evolution of the probability density...

Application of Legendre spectral-collocation method to delay differential and stochastic delay differential equation

Science.gov (United States)

Khan, Sami Ullah; Ali, Ishtiaq

2018-03-01

Explicit solutions to delay differential equation (DDE) and stochastic delay differential equation (SDDE) can rarely be obtained, therefore numerical methods are adopted to solve these DDE and SDDE. While on the other hand due to unstable nature of both DDE and SDDE numerical solutions are also not straight forward and required more attention. In this study, we derive an efficient numerical scheme for DDE and SDDE based on Legendre spectral-collocation method, which proved to be numerical methods that can significantly speed up the computation. The method transforms the given differential equation into a matrix equation by means of Legendre collocation points which correspond to a system of algebraic equations with unknown Legendre coefficients. The efficiency of the proposed method is confirmed by some numerical examples. We found that our numerical technique has a very good agreement with other methods with less computational effort.

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr; Vilmart, Gilles; Zygalakis, Konstantinos C.

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer

Parameter estimation in a simple stochastic differential equation for phytoplankton modelling

DEFF Research Database (Denmark)

Møller, Jan Kloppenborg; Madsen, Henrik; Carstensen, Jacob

2011-01-01

The use of stochastic differential equations (SDEs) for simulation of aquatic ecosystems has attracted increasing attention in recent years. The SDE setting also provides the opportunity for statistical estimation of ecosystem parameters. We present an estimation procedure, based on Kalman...

From quantum stochastic differential equations to Gisin-Percival state diffusion

Science.gov (United States)

Parthasarathy, K. R.; Usha Devi, A. R.

2017-08-01

Starting from the quantum stochastic differential equations of Hudson and Parthasarathy [Commun. Math. Phys. 93, 301 (1984)] and exploiting the Wiener-Itô-Segal isomorphism between the boson Fock reservoir space Γ (L2(R+ ) ⊗(Cn⊕Cn ) ) and the Hilbert space L2(μ ) , where μ is the Wiener probability measure of a complex n-dimensional vector-valued standard Brownian motion {B (t ) ,t ≥0 } , we derive a non-linear stochastic Schrödinger equation describing a classical diffusion of states of a quantum system, driven by the Brownian motion B. Changing this Brownian motion by an appropriate Girsanov transformation, we arrive at the Gisin-Percival state diffusion equation [N. Gisin and J. Percival, J. Phys. A 167, 315 (1992)]. This approach also yields an explicit solution of the Gisin-Percival equation, in terms of the Hudson-Parthasarathy unitary process and a randomized Weyl displacement process. Irreversible dynamics of system density operators described by the well-known Gorini-Kossakowski-Sudarshan-Lindblad master equation is unraveled by coarse-graining over the Gisin-Percival quantum state trajectories.

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Angstmann, C.N.; Donnelly, I.C. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Henry, B.I., E-mail: B.Henry@unsw.edu.au [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Jacobs, B.A. [School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050 (South Africa); DST–NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS) (South Africa); Langlands, T.A.M. [Department of Mathematics and Computing, University of Southern Queensland, Toowoomba QLD 4350 (Australia); Nichols, J.A. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia)

2016-02-15

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also show that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.

Polynomial asymptotic stability of damped stochastic differential equations

Directory of Open Access Journals (Sweden)

John Appleby

2004-08-01

Full Text Available The paper studies the polynomial convergence of solutions of a scalar nonlinear It\\^{o} stochastic differential equation\\[dX(t = -f(X(t\\,dt + \\sigma(t\\,dB(t\\] where it is known, {\\it a priori}, that $\\lim_{t\\rightarrow\\infty} X(t=0$, a.s. The intensity of the stochastic perturbation $\\sigma$ is a deterministic, continuous and square integrable function, which tends to zero more quickly than a polynomially decaying function. The function $f$ obeys $\\lim_{x\\rightarrow 0}\\mbox{sgn}(xf(x/|x|^\\beta = a$, for some $\\beta>1$, and $a>0$.We study two asymptotic regimes: when $\\sigma$ tends to zero sufficiently quickly the polynomial decay rate of solutions is the same as for the deterministic equation (when $\\sigma\\equiv0$. When $\\sigma$ decays more slowly, a weaker almost sure polynomial upper bound on the decay rate of solutions is established. Results which establish the necessity for $\\sigma$ to decay polynomially in order to guarantee the almost sure polynomial decay of solutions are also proven.

Modelling biochemical reaction systems by stochastic differential equations with reflection.

Science.gov (United States)

Niu, Yuanling; Burrage, Kevin; Chen, Luonan

2016-05-07

In this paper, we gave a new framework for modelling and simulating biochemical reaction systems by stochastic differential equations with reflection not in a heuristic way but in a mathematical way. The model is computationally efficient compared with the discrete-state Markov chain approach, and it ensures that both analytic and numerical solutions remain in a biologically plausible region. Specifically, our model mathematically ensures that species numbers lie in the domain D, which is a physical constraint for biochemical reactions, in contrast to the previous models. The domain D is actually obtained according to the structure of the corresponding chemical Langevin equations, i.e., the boundary is inherent in the biochemical reaction system. A variant of projection method was employed to solve the reflected stochastic differential equation model, and it includes three simple steps, i.e., Euler-Maruyama method was applied to the equations first, and then check whether or not the point lies within the domain D, and if not perform an orthogonal projection. It is found that the projection onto the closure D¯ is the solution to a convex quadratic programming problem. Thus, existing methods for the convex quadratic programming problem can be employed for the orthogonal projection map. Numerical tests on several important problems in biological systems confirmed the efficiency and accuracy of this approach. Copyright © 2016 Elsevier Ltd. All rights reserved.

Modelling the heat dynamics of a building using stochastic differential equations

DEFF Research Database (Denmark)

Andersen, Klaus Kaae; Madsen, Henrik; Hansen, Lars Henrik

2000-01-01

estimation and model validation, while physical knowledge is used in forming the model structure. The suggested lumped parameter model is thus based on thermodynamics and formulated as a system of stochastic differential equations. Due to the continuous time formulation the parameters of the model...

Polynomial chaos functions and stochastic differential equations

International Nuclear Information System (INIS)

Williams, M.M.R.

2006-01-01

The Karhunen-Loeve procedure and the associated polynomial chaos expansion have been employed to solve a simple first order stochastic differential equation which is typical of transport problems. Because the equation has an analytical solution, it provides a useful test of the efficacy of polynomial chaos. We find that the convergence is very rapid in some cases but that the increased complexity associated with many random variables can lead to very long computational times. The work is illustrated by exact and approximate solutions for the mean, variance and the probability distribution itself. The usefulness of a white noise approximation is also assessed. Extensive numerical results are given which highlight the weaknesses and strengths of polynomial chaos. The general conclusion is that the method is promising but requires further detailed study by application to a practical problem in transport theory

Stochastic differential equations for quantum dynamics of spin-boson networks

International Nuclear Information System (INIS)

Mandt, Stephan; Sadri, Darius; Houck, Andrew A; Türeci, Hakan E

2015-01-01

A popular approach in quantum optics is to map a master equation to a stochastic differential equation, where quantum effects manifest themselves through noise terms. We generalize this approach based on the positive-P representation to systems involving spin, in particular networks or lattices of interacting spins and bosons. We test our approach on a driven dimer of spins and photons, compare it to the master equation, and predict a novel dynamic phase transition in this system. Our numerical approach has scaling advantages over existing methods, but typically requires regularization in terms of drive and dissipation. (paper)

Fractional Order Stochastic Differential Equation with Application in European Option Pricing

Directory of Open Access Journals (Sweden)

Qing Li

2014-01-01

Full Text Available Memory effect is an important phenomenon in financial systems, and a number of research works have been carried out to study the long memory in the financial markets. In recent years, fractional order ordinary differential equation is used as an effective instrument for describing the memory effect in complex systems. In this paper, we establish a fractional order stochastic differential equation (FSDE model to describe the effect of trend memory in financial pricing. We, then, derive a European option pricing formula based on the FSDE model and prove the existence of the trend memory (i.e., the mean value function in the option pricing formula when the Hurst index is between 0.5 and 1. In addition, we make a comparison analysis between our proposed model, the classic Black-Scholes model, and the stochastic model with fractional Brownian motion. Numerical results suggest that our model leads to more accurate and lower standard deviation in the empirical study.

Analytic Approximation of the Solutions of Stochastic Differential Delay Equations with Poisson Jump and Markovian Switching

Directory of Open Access Journals (Sweden)

Hua Yang

2012-01-01

Full Text Available We are concerned with the stochastic differential delay equations with Poisson jump and Markovian switching (SDDEsPJMSs. Most SDDEsPJMSs cannot be solved explicitly as stochastic differential equations. Therefore, numerical solutions have become an important issue in the study of SDDEsPJMSs. The key contribution of this paper is to investigate the strong convergence between the true solutions and the numerical solutions to SDDEsPJMSs when the drift and diffusion coefficients are Taylor approximations.

Stochastic partial differential equations a modeling, white noise functional approach

CERN Document Server

Holden, Helge; Ubøe, Jan; Zhang, Tusheng

1996-01-01

This book is based on research that, to a large extent, started around 1990, when a research project on fluid flow in stochastic reservoirs was initiated by a group including some of us with the support of VISTA, a research coopera­ tion between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap A.S. (Statoil). The purpose of the project was to use stochastic partial differential equations (SPDEs) to describe the flow of fluid in a medium where some of the parameters, e.g., the permeability, were stochastic or "noisy". We soon realized that the theory of SPDEs at the time was insufficient to handle such equations. Therefore it became our aim to develop a new mathematically rigorous theory that satisfied the following conditions. 1) The theory should be physically meaningful and realistic, and the corre­ sponding solutions should make sense physically and should be useful in applications. 2) The theory should be general enough to handle many of the interesting SPDEs that occur in r...

XMDS2: Fast, scalable simulation of coupled stochastic partial differential equations

Science.gov (United States)

Dennis, Graham R.; Hope, Joseph J.; Johnsson, Mattias T.

2013-01-01

XMDS2 is a cross-platform, GPL-licensed, open source package for numerically integrating initial value problems that range from a single ordinary differential equation up to systems of coupled stochastic partial differential equations. The equations are described in a high-level XML-based script, and the package generates low-level optionally parallelised C++ code for the efficient solution of those equations. It combines the advantages of high-level simulations, namely fast and low-error development, with the speed, portability and scalability of hand-written code. XMDS2 is a complete redesign of the XMDS package, and features support for a much wider problem space while also producing faster code. Program summaryProgram title: XMDS2 Catalogue identifier: AENK_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 2 No. of lines in distributed program, including test data, etc.: 872490 No. of bytes in distributed program, including test data, etc.: 45522370 Distribution format: tar.gz Programming language: Python and C++. Computer: Any computer with a Unix-like system, a C++ compiler and Python. Operating system: Any Unix-like system; developed under Mac OS X and GNU/Linux. RAM: Problem dependent (roughly 50 bytes per grid point) Classification: 4.3, 6.5. External routines: The external libraries required are problem-dependent. Uses FFTW3 Fourier transforms (used only for FFT-based spectral methods), dSFMT random number generation (used only for stochastic problems), MPI message-passing interface (used only for distributed problems), HDF5, GNU Scientific Library (used only for Bessel-based spectral methods) and a BLAS implementation (used only for non-FFT-based spectral methods). Nature of problem: General coupled initial-value stochastic partial differential equations. Solution method: Spectral method

Parametric inference for stochastic differential equations: a smooth and match approach

NARCIS (Netherlands)

Gugushvili, S.; Spreij, P.

2012-01-01

We study the problem of parameter estimation for a univariate discretely observed ergodic diffusion process given as a solution to a stochastic differential equation. The estimation procedure we propose consists of two steps. In the first step, which is referred to as a smoothing step, we smooth the

Exercise effects in a virtual type 1 diabetes patient: Using stochastic differential equations for model extension

DEFF Research Database (Denmark)

Duun-Henriksen, Anne Katrine; Schmidt, S.; Nørgaard, K.

2013-01-01

extension incorporating exercise effects on insulin and glucose dynamics. Our model is constructed as a stochastic state space model consisting of a set of stochastic differential equations (SDEs). In a stochastic state space model, the residual error is split into random measurement error...

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-09-19

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-01-01

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Semilinear Kolmogorov Equations and Applications to Stochastic Optimal Control

International Nuclear Information System (INIS)

Masiero, Federica

2005-01-01

Semilinear parabolic differential equations are solved in a mild sense in an infinite-dimensional Hilbert space. Applications to stochastic optimal control problems are studied by solving the associated Hamilton-Jacobi-Bellman equation. These results are applied to some controlled stochastic partial differential equations

EXISTENCE OF SOLUTION TO NONLINEAR SECOND ORDER NEUTRAL STOCHASTIC DIFFERENTIAL EQUATIONS WITH DELAY

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

This paper is concerned with nonlinear second order neutral stochastic differential equations with delay in a Hilbert space. Sufficient conditions for the existence of solution to the system are obtained by Picard iterations.

Statistical inference for discrete-time samples from affine stochastic delay differential equations

DEFF Research Database (Denmark)

Küchler, Uwe; Sørensen, Michael

2013-01-01

Statistical inference for discrete time observations of an affine stochastic delay differential equation is considered. The main focus is on maximum pseudo-likelihood estimators, which are easy to calculate in practice. A more general class of prediction-based estimating functions is investigated...

An estimator for the relative entropy rate of path measures for stochastic differential equations

Energy Technology Data Exchange (ETDEWEB)

Opper, Manfred, E-mail: manfred.opper@tu-berlin.de

2017-02-01

We address the problem of estimating the relative entropy rate (RER) for two stochastic processes described by stochastic differential equations. For the case where the drift of one process is known analytically, but one has only observations from the second process, we use a variational bound on the RER to construct an estimator.

Estimating the parameters of stochastic differential equations using a criterion function based on the Kolmogorov-Smirnov statistic

OpenAIRE

McDonald, A. David; Sandal, Leif Kristoffer

1998-01-01

Estimation of parameters in the drift and diffusion terms of stochastic differential equations involves simulation and generally requires substantial data sets. We examine a method that can be applied when available time series are limited to less than 20 observations per replication. We compare and contrast parameter estimation for linear and nonlinear first-order stochastic differential equations using two criterion functions: one based on a Chi-square statistic, put forward by Hurn and Lin...

Mixed Effects Modeling Using Stochastic Differential Equations: Illustrated by Pharmacokinetic Data of Nicotinic Acid in Obese Zucker Rats.

Science.gov (United States)

Leander, Jacob; Almquist, Joachim; Ahlström, Christine; Gabrielsson, Johan; Jirstrand, Mats

2015-05-01

Inclusion of stochastic differential equations in mixed effects models provides means to quantify and distinguish three sources of variability in data. In addition to the two commonly encountered sources, measurement error and interindividual variability, we also consider uncertainty in the dynamical model itself. To this end, we extend the ordinary differential equation setting used in nonlinear mixed effects models to include stochastic differential equations. The approximate population likelihood is derived using the first-order conditional estimation with interaction method and extended Kalman filtering. To illustrate the application of the stochastic differential mixed effects model, two pharmacokinetic models are considered. First, we use a stochastic one-compartmental model with first-order input and nonlinear elimination to generate synthetic data in a simulated study. We show that by using the proposed method, the three sources of variability can be successfully separated. If the stochastic part is neglected, the parameter estimates become biased, and the measurement error variance is significantly overestimated. Second, we consider an extension to a stochastic pharmacokinetic model in a preclinical study of nicotinic acid kinetics in obese Zucker rats. The parameter estimates are compared between a deterministic and a stochastic NiAc disposition model, respectively. Discrepancies between model predictions and observations, previously described as measurement noise only, are now separated into a comparatively lower level of measurement noise and a significant uncertainty in model dynamics. These examples demonstrate that stochastic differential mixed effects models are useful tools for identifying incomplete or inaccurate model dynamics and for reducing potential bias in parameter estimates due to such model deficiencies.

Approximate controllability of Sobolev type fractional stochastic nonlocal nonlinear differential equations in Hilbert spaces

Directory of Open Access Journals (Sweden)

Mourad Kerboua

2014-12-01

Full Text Available We introduce a new notion called fractional stochastic nonlocal condition, and then we study approximate controllability of class of fractional stochastic nonlinear differential equations of Sobolev type in Hilbert spaces. We use Hölder's inequality, fixed point technique, fractional calculus, stochastic analysis and methods adopted directly from deterministic control problems for the main results. A new set of sufficient conditions is formulated and proved for the fractional stochastic control system to be approximately controllable. An example is given to illustrate the abstract results.

Efficient Estimating Functions for Stochastic Differential Equations

DEFF Research Database (Denmark)

Jakobsen, Nina Munkholt

The overall topic of this thesis is approximate martingale estimating function-based estimationfor solutions of stochastic differential equations, sampled at high frequency. Focuslies on the asymptotic properties of the estimators. The first part of the thesis deals with diffusions observed over...... a fixed time interval. Rate optimal and effcient estimators areobtained for a one-dimensional diffusion parameter. Stable convergence in distribution isused to achieve a practically applicable Gaussian limit distribution for suitably normalisedestimators. In a simulation example, the limit distributions...... multidimensional parameter. Conditions for rate optimality and effciency of estimatorsof drift-jump and diffusion parameters are given in some special cases. Theseconditions are found to extend the pre-existing conditions applicable to continuous diffusions,and impose much stronger requirements on the estimating...

Numerical approximations of stochastic differential equations with non-globally Lipschitz continuous coefficients

CERN Document Server

Hutzenthaler, Martin

2015-01-01

Many stochastic differential equations (SDEs) in the literature have a superlinearly growing nonlinearity in their drift or diffusion coefficient. Unfortunately, moments of the computationally efficient Euler-Maruyama approximation method diverge for these SDEs in finite time. This article develops a general theory based on rare events for studying integrability properties such as moment bounds for discrete-time stochastic processes. Using this approach, the authors establish moment bounds for fully and partially drift-implicit Euler methods and for a class of new explicit approximation method

Stability of Nonlinear Neutral Stochastic Functional Differential Equations

Directory of Open Access Journals (Sweden)

Minggao Xue

2010-01-01

Full Text Available Neutral stochastic functional differential equations (NSFDEs have recently been studied intensively. The well-known conditions imposed for the existence and uniqueness and exponential stability of the global solution are the local Lipschitz condition and the linear growth condition. Therefore, the existing results cannot be applied to many important nonlinear NSFDEs. The main aim of this paper is to remove the linear growth condition and establish a Khasminskii-type test for nonlinear NSFDEs. New criteria not only cover a wide class of highly nonlinear NSFDEs but they can also be verified much more easily than the classical criteria. Finally, several examples are given to illustrate main results.

Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations

DEFF Research Database (Denmark)

Tornøe, Christoffer Wenzel; Overgaard, Rune Viig; Agerso, H.

2005-01-01

of noise: a measurement and a system noise term. The measurement noise represents uncorrelated error due to, for example, assay error while the system noise accounts for structural misspecifications, approximations of the dynamical model, and true random physiological fluctuations. Since the system noise...... degarelix. Conclusions. The EKF-based algorithm was successfully implemented in NONMEM for parameter estimation in population PK/PD models described by systems of SDEs. The example indicated that it was possible to pinpoint structural model deficiencies, and that valuable information may be obtained......Purpose. The objective of the present analysis was to explore the use of stochastic differential equations (SDEs) in population pharmacokinetic/pharmacodynamic (PK/PD) modeling. Methods. The intra-individual variability in nonlinear mixed-effects models based on SDEs is decomposed into two types...

Model identification using stochastic differential equation grey-box models in diabetes.

Science.gov (United States)

Duun-Henriksen, Anne Katrine; Schmidt, Signe; Røge, Rikke Meldgaard; Møller, Jonas Bech; Nørgaard, Kirsten; Jørgensen, John Bagterp; Madsen, Henrik

2013-03-01

The acceptance of virtual preclinical testing of control algorithms is growing and thus also the need for robust and reliable models. Models based on ordinary differential equations (ODEs) can rarely be validated with standard statistical tools. Stochastic differential equations (SDEs) offer the possibility of building models that can be validated statistically and that are capable of predicting not only a realistic trajectory, but also the uncertainty of the prediction. In an SDE, the prediction error is split into two noise terms. This separation ensures that the errors are uncorrelated and provides the possibility to pinpoint model deficiencies. An identifiable model of the glucoregulatory system in a type 1 diabetes mellitus (T1DM) patient is used as the basis for development of a stochastic-differential-equation-based grey-box model (SDE-GB). The parameters are estimated on clinical data from four T1DM patients. The optimal SDE-GB is determined from likelihood-ratio tests. Finally, parameter tracking is used to track the variation in the "time to peak of meal response" parameter. We found that the transformation of the ODE model into an SDE-GB resulted in a significant improvement in the prediction and uncorrelated errors. Tracking of the "peak time of meal absorption" parameter showed that the absorption rate varied according to meal type. This study shows the potential of using SDE-GBs in diabetes modeling. Improved model predictions were obtained due to the separation of the prediction error. SDE-GBs offer a solid framework for using statistical tools for model validation and model development. © 2013 Diabetes Technology Society.

A stochastic differential equation analysis of cerebrospinal fluid dynamics.

Science.gov (United States)

Raman, Kalyan

2011-01-18

Clinical measurements of intracranial pressure (ICP) over time show fluctuations around the deterministic time path predicted by a classic mathematical model in hydrocephalus research. Thus an important issue in mathematical research on hydrocephalus remains unaddressed--modeling the effect of noise on CSF dynamics. Our objective is to mathematically model the noise in the data. The classic model relating the temporal evolution of ICP in pressure-volume studies to infusions is a nonlinear differential equation based on natural physical analogies between CSF dynamics and an electrical circuit. Brownian motion was incorporated into the differential equation describing CSF dynamics to obtain a nonlinear stochastic differential equation (SDE) that accommodates the fluctuations in ICP. The SDE is explicitly solved and the dynamic probabilities of exceeding critical levels of ICP under different clinical conditions are computed. A key finding is that the probabilities display strong threshold effects with respect to noise. Above the noise threshold, the probabilities are significantly influenced by the resistance to CSF outflow and the intensity of the noise. Fluctuations in the CSF formation rate increase fluctuations in the ICP and they should be minimized to lower the patient's risk. The nonlinear SDE provides a scientific methodology for dynamic risk management of patients. The dynamic output of the SDE matches the noisy ICP data generated by the actual intracranial dynamics of patients better than the classic model used in prior research.

Probabilistic Forecast of Wind Power Generation by Stochastic Differential Equation Models

KAUST Repository

Elkantassi, Soumaya

2017-04-01

Reliable forecasting of wind power generation is crucial to optimal control of costs in generation of electricity with respect to the electricity demand. Here, we propose and analyze stochastic wind power forecast models described by parametrized stochastic differential equations, which introduce appropriate fluctuations in numerical forecast outputs. We use an approximate maximum likelihood method to infer the model parameters taking into account the time correlated sets of data. Furthermore, we study the validity and sensitivity of the parameters for each model. We applied our models to Uruguayan wind power production as determined by historical data and corresponding numerical forecasts for the period of March 1 to May 31, 2016.

Higher-order stochastic differential equations and the positive Wigner function

Science.gov (United States)

Drummond, P. D.

2017-12-01

General higher-order stochastic processes that correspond to any diffusion-type tensor of higher than second order are obtained. The relationship of multivariate higher-order stochastic differential equations with tensor decomposition theory and tensor rank is explained. Techniques for generating the requisite complex higher-order noise are proved to exist either using polar coordinates and γ distributions, or from products of Gaussian variates. This method is shown to allow the calculation of the dynamics of the Wigner function, after it is extended to a complex phase space. The results are illustrated physically through dynamical calculations of the positive Wigner distribution for three-mode parametric downconversion, widely used in quantum optics. The approach eliminates paradoxes arising from truncation of the higher derivative terms in Wigner function time evolution. Anomalous results of negative populations and vacuum scattering found in truncated Wigner quantum simulations in quantum optics and Bose-Einstein condensate dynamics are shown not to occur with this type of stochastic theory.

Backward-stochastic-differential-equation approach to modeling of gene expression.

Science.gov (United States)

Shamarova, Evelina; Chertovskih, Roman; Ramos, Alexandre F; Aguiar, Paulo

2017-03-01

In this article, we introduce a backward method to model stochastic gene expression and protein-level dynamics. The protein amount is regarded as a diffusion process and is described by a backward stochastic differential equation (BSDE). Unlike many other SDE techniques proposed in the literature, the BSDE method is backward in time; that is, instead of initial conditions it requires the specification of end-point ("final") conditions, in addition to the model parametrization. To validate our approach we employ Gillespie's stochastic simulation algorithm (SSA) to generate (forward) benchmark data, according to predefined gene network models. Numerical simulations show that the BSDE method is able to correctly infer the protein-level distributions that preceded a known final condition, obtained originally from the forward SSA. This makes the BSDE method a powerful systems biology tool for time-reversed simulations, allowing, for example, the assessment of the biological conditions (e.g., protein concentrations) that preceded an experimentally measured event of interest (e.g., mitosis, apoptosis, etc.).

ADAPTIVE METHODS FOR STOCHASTIC DIFFERENTIAL EQUATIONS VIA NATURAL EMBEDDINGS AND REJECTION SAMPLING WITH MEMORY.

Science.gov (United States)

Rackauckas, Christopher; Nie, Qing

2017-01-01

Adaptive time-stepping with high-order embedded Runge-Kutta pairs and rejection sampling provides efficient approaches for solving differential equations. While many such methods exist for solving deterministic systems, little progress has been made for stochastic variants. One challenge in developing adaptive methods for stochastic differential equations (SDEs) is the construction of embedded schemes with direct error estimates. We present a new class of embedded stochastic Runge-Kutta (SRK) methods with strong order 1.5 which have a natural embedding of strong order 1.0 methods. This allows for the derivation of an error estimate which requires no additional function evaluations. Next we derive a general method to reject the time steps without losing information about the future Brownian path termed Rejection Sampling with Memory (RSwM). This method utilizes a stack data structure to do rejection sampling, costing only a few floating point calculations. We show numerically that the methods generate statistically-correct and tolerance-controlled solutions. Lastly, we show that this form of adaptivity can be applied to systems of equations, and demonstrate that it solves a stiff biological model 12.28x faster than common fixed timestep algorithms. Our approach only requires the solution to a bridging problem and thus lends itself to natural generalizations beyond SDEs.

The Dynamic Programming Method of Stochastic Differential Game for Functional Forward-Backward Stochastic System

Directory of Open Access Journals (Sweden)

Shaolin Ji

2013-01-01

Full Text Available This paper is devoted to a stochastic differential game (SDG of decoupled functional forward-backward stochastic differential equation (FBSDE. For our SDG, the associated upper and lower value functions of the SDG are defined through the solution of controlled functional backward stochastic differential equations (BSDEs. Applying the Girsanov transformation method introduced by Buckdahn and Li (2008, the upper and the lower value functions are shown to be deterministic. We also generalize the Hamilton-Jacobi-Bellman-Isaacs (HJBI equations to the path-dependent ones. By establishing the dynamic programming principal (DPP, we derive that the upper and the lower value functions are the viscosity solutions of the corresponding upper and the lower path-dependent HJBI equations, respectively.

Variational and potential formulation for stochastic partial differential equations

International Nuclear Information System (INIS)

Munoz S, A G; Ojeda, J; Sierra D, P; Soldovieri, T

2006-01-01

Recently there has been interest in finding a potential formulation for stochastic partial differential equations (SPDEs). The rationale behind this idea lies in obtaining all the dynamical information of the system under study from one single expression. In this letter we formally provide a general Lagrangian formalism for SPDEs using the Hojman et al method. We show that it is possible to write the corresponding effective potential starting from an s-equivalent Lagrangian, and that this potential is able to reproduce all the dynamics of the system once a special differential operator has been applied. This procedure can be used to study the complete time evolution and spatial inhomogeneities of the system under consideration, and is also suitable for the statistical mechanics description of the problem. (letter to the editor)

Invariant measures for stochastic nonlinear beam and wave equations

Czech Academy of Sciences Publication Activity Database

Brzezniak, Z.; Ondreját, Martin; Seidler, Jan

2016-01-01

Roč. 260, č. 5 (2016), s. 4157-4179 ISSN 0022-0396 R&D Projects: GA ČR GAP201/10/0752 Institutional support: RVO:67985556 Keywords : stochastic partial differential equation * stochastic beam equation * stochastic wave equation * invariant measure Subject RIV: BA - General Mathematics Impact factor: 1.988, year: 2016 http://library.utia.cas.cz/separaty/2016/SI/ondrejat-0453412.pdf

A stochastic differential equation analysis of cerebrospinal fluid dynamics

Directory of Open Access Journals (Sweden)

Raman Kalyan

2011-01-01

Full Text Available Abstract Background Clinical measurements of intracranial pressure (ICP over time show fluctuations around the deterministic time path predicted by a classic mathematical model in hydrocephalus research. Thus an important issue in mathematical research on hydrocephalus remains unaddressed--modeling the effect of noise on CSF dynamics. Our objective is to mathematically model the noise in the data. Methods The classic model relating the temporal evolution of ICP in pressure-volume studies to infusions is a nonlinear differential equation based on natural physical analogies between CSF dynamics and an electrical circuit. Brownian motion was incorporated into the differential equation describing CSF dynamics to obtain a nonlinear stochastic differential equation (SDE that accommodates the fluctuations in ICP. Results The SDE is explicitly solved and the dynamic probabilities of exceeding critical levels of ICP under different clinical conditions are computed. A key finding is that the probabilities display strong threshold effects with respect to noise. Above the noise threshold, the probabilities are significantly influenced by the resistance to CSF outflow and the intensity of the noise. Conclusions Fluctuations in the CSF formation rate increase fluctuations in the ICP and they should be minimized to lower the patient's risk. The nonlinear SDE provides a scientific methodology for dynamic risk management of patients. The dynamic output of the SDE matches the noisy ICP data generated by the actual intracranial dynamics of patients better than the classic model used in prior research.

A Stochastic Differential Equation Model for the Spread of HIV amongst People Who Inject Drugs

Directory of Open Access Journals (Sweden)

Yanfeng Liang

2016-01-01

Full Text Available We introduce stochasticity into the deterministic differential equation model for the spread of HIV amongst people who inject drugs (PWIDs studied by Greenhalgh and Hay (1997. This was based on the original model constructed by Kaplan (1989 which analyses the behaviour of HIV/AIDS amongst a population of PWIDs. We derive a stochastic differential equation (SDE for the fraction of PWIDs who are infected with HIV at time. The stochasticity is introduced using the well-known standard technique of parameter perturbation. We first prove that the resulting SDE for the fraction of infected PWIDs has a unique solution in (0, 1 provided that some infected PWIDs are initially present and next construct the conditions required for extinction and persistence. Furthermore, we show that there exists a stationary distribution for the persistence case. Simulations using realistic parameter values are then constructed to illustrate and support our theoretical results. Our results provide new insight into the spread of HIV amongst PWIDs. The results show that the introduction of stochastic noise into a model for the spread of HIV amongst PWIDs can cause the disease to die out in scenarios where deterministic models predict disease persistence.

A Stochastic Differential Equation Model for the Spread of HIV amongst People Who Inject Drugs.

Science.gov (United States)

Liang, Yanfeng; Greenhalgh, David; Mao, Xuerong

2016-01-01

We introduce stochasticity into the deterministic differential equation model for the spread of HIV amongst people who inject drugs (PWIDs) studied by Greenhalgh and Hay (1997). This was based on the original model constructed by Kaplan (1989) which analyses the behaviour of HIV/AIDS amongst a population of PWIDs. We derive a stochastic differential equation (SDE) for the fraction of PWIDs who are infected with HIV at time. The stochasticity is introduced using the well-known standard technique of parameter perturbation. We first prove that the resulting SDE for the fraction of infected PWIDs has a unique solution in (0, 1) provided that some infected PWIDs are initially present and next construct the conditions required for extinction and persistence. Furthermore, we show that there exists a stationary distribution for the persistence case. Simulations using realistic parameter values are then constructed to illustrate and support our theoretical results. Our results provide new insight into the spread of HIV amongst PWIDs. The results show that the introduction of stochastic noise into a model for the spread of HIV amongst PWIDs can cause the disease to die out in scenarios where deterministic models predict disease persistence.

A micro-macro acceleration method for the Monte Carlo simulation of stochastic differential equations

DEFF Research Database (Denmark)

Debrabant, Kristian; Samaey, Giovanni; Zieliński, Przemysław

2017-01-01

We present and analyse a micro-macro acceleration method for the Monte Carlo simulation of stochastic differential equations with separation between the (fast) time-scale of individual trajectories and the (slow) time-scale of the macroscopic function of interest. The algorithm combines short...

Yosida approximations of stochastic differential equations in infinite dimensions and applications

CERN Document Server

Govindan, T E

2016-01-01

This research monograph brings together, for the first time, the varied literature on Yosida approximations of stochastic differential equations (SDEs) in infinite dimensions and their applications into a single cohesive work. The author provides a clear and systematic introduction to the Yosida approximation method and justifies its power by presenting its applications in some practical topics such as stochastic stability and stochastic optimal control. The theory assimilated spans more than 35 years of mathematics, but is developed slowly and methodically in digestible pieces. The book begins with a motivational chapter that introduces the reader to several different models that play recurring roles throughout the book as the theory is unfolded, and invites readers from different disciplines to see immediately that the effort required to work through the theory that follows is worthwhile. From there, the author presents the necessary prerequisite material, and then launches the reader into the main discussi...

On Exponential Hedging and Related Quadratic Backward Stochastic Differential Equations

International Nuclear Information System (INIS)

Sekine, Jun

2006-01-01

The dual optimization problem for the exponential hedging problem is addressed with a cone constraint. Without boundedness conditions on the terminal payoff and the drift of the Ito-type controlled process, the backward stochastic differential equation, which has a quadratic growth term in the drift, is derived as a necessary and sufficient condition for optimality via a variational method and dynamic programming. Further, solvable situations are given, in which the value and the optimizer are expressed in closed forms with the help of the Clark-Haussmann-Ocone formula

Derivation of stochastic differential equations for scrape-off layer plasma fluctuations from experimentally measured statistics

Energy Technology Data Exchange (ETDEWEB)

Mekkaoui, Abdessamad [IEK-4 Forschungszentrum Juelich 52428 (Germany)

2013-07-01

A method to derive stochastic differential equations for intermittent plasma density dynamics in magnetic fusion edge plasma is presented. It uses a measured first four moments (mean, variance, Skewness and Kurtosis) and the correlation time of turbulence to write a Pearson equation for the probability distribution function of fluctuations. The Fokker-Planck equation is then used to derive a Langevin equation for the plasma density fluctuations. A theoretical expectations are used as a constraints to fix the nonlinearity structure of the stochastic differential equation. In particular when the quadratically nonlinear dynamics is assumed, then it is shown that the plasma density is driven by a multiplicative Wiener process and evolves on the turbulence correlation time scale, while the linear growth is quadratically damped by the fluctuation level. Strong criteria for statistical discrimination of experimental time series are proposed as an alternative to the Kurtosis-Skewness scaling. This scaling is broadly used in contemporary literature to characterize edge turbulence, but it is inappropriate because a large family of distributions could share this scaling. Strong criteria allow us to focus on the relevant candidate distribution and approach a nonlinear structure of edge turbulence model.

Stochastic parameterizing manifolds and non-Markovian reduced equations stochastic manifolds for nonlinear SPDEs II

CERN Document Server

Chekroun, Mickaël D; Wang, Shouhong

2015-01-01

In this second volume, a general approach is developed to provide approximate parameterizations of the "small" scales by the "large" ones for a broad class of stochastic partial differential equations (SPDEs). This is accomplished via the concept of parameterizing manifolds (PMs), which are stochastic manifolds that improve, for a given realization of the noise, in mean square error the partial knowledge of the full SPDE solution when compared to its projection onto some resolved modes. Backward-forward systems are designed to give access to such PMs in practice. The key idea consists of representing the modes with high wave numbers as a pullback limit depending on the time-history of the modes with low wave numbers. Non-Markovian stochastic reduced systems are then derived based on such a PM approach. The reduced systems take the form of stochastic differential equations involving random coefficients that convey memory effects. The theory is illustrated on a stochastic Burgers-type equation.

Picard Approximation of Stochastic Differential Equations and Application to LIBOR Models

DEFF Research Database (Denmark)

Papapantoleon, Antonis; Skovmand, David

The aim of this work is to provide fast and accurate approximation schemes for the Monte Carlo pricing of derivatives in LIBOR market models. Standard methods can be applied to solve the stochastic differential equations of the successive LIBOR rates but the methods are generally slow. Our...... exponential to quadratic using truncated expansions of the product terms. We include numerical illustrations of the accuracy and speed of our method pricing caplets, swaptions and forward rate agreements....

Degenerate parabolic stochastic partial differential equations

Czech Academy of Sciences Publication Activity Database

span class="emphasis">Hofmanová, Martinaspan>

2013-01-01

Roč. 123, č. 12 (2013), s. 4294-4336 ISSN 0304-4149 R&D Projects: GA ČR GAP201/10/0752 Institutional support: RVO:67985556 Keywords : kinetic solutions * degenerate stochastic parabolic equations Subject RIV: BA - General Mathematics Impact factor: 1.046, year: 2013 http://library.utia.cas.cz/separaty/2013/SI/hofmanova-0397241.pdf

Nonintrusive Polynomial Chaos Expansions for Sensitivity Analysis in Stochastic Differential Equations

KAUST Repository

Jimenez, M. Navarro; Le Maî tre, O. P.; Knio, Omar

2017-01-01

A Galerkin polynomial chaos (PC) method was recently proposed to perform variance decomposition and sensitivity analysis in stochastic differential equations (SDEs), driven by Wiener noise and involving uncertain parameters. The present paper extends the PC method to nonintrusive approaches enabling its application to more complex systems hardly amenable to stochastic Galerkin projection methods. We also discuss parallel implementations and the variance decomposition of the derived quantity of interest within the framework of nonintrusive approaches. In particular, a novel hybrid PC-sampling-based strategy is proposed in the case of nonsmooth quantities of interest (QoIs) but smooth SDE solution. Numerical examples are provided that illustrate the decomposition of the variance of QoIs into contributions arising from the uncertain parameters, the inherent stochastic forcing, and joint effects. The simulations are also used to support a brief analysis of the computational complexity of the method, providing insight on the types of problems that would benefit from the present developments.

Nonintrusive Polynomial Chaos Expansions for Sensitivity Analysis in Stochastic Differential Equations

KAUST Repository

Jimenez, M. Navarro

2017-04-18

A Galerkin polynomial chaos (PC) method was recently proposed to perform variance decomposition and sensitivity analysis in stochastic differential equations (SDEs), driven by Wiener noise and involving uncertain parameters. The present paper extends the PC method to nonintrusive approaches enabling its application to more complex systems hardly amenable to stochastic Galerkin projection methods. We also discuss parallel implementations and the variance decomposition of the derived quantity of interest within the framework of nonintrusive approaches. In particular, a novel hybrid PC-sampling-based strategy is proposed in the case of nonsmooth quantities of interest (QoIs) but smooth SDE solution. Numerical examples are provided that illustrate the decomposition of the variance of QoIs into contributions arising from the uncertain parameters, the inherent stochastic forcing, and joint effects. The simulations are also used to support a brief analysis of the computational complexity of the method, providing insight on the types of problems that would benefit from the present developments.

Perturbation Solutions for Random Linear Structural Systems subject to Random Excitation using Stochastic Differential Equations

DEFF Research Database (Denmark)

Köyluoglu, H.U.; Nielsen, Søren R.K.; Cakmak, A.S.

1994-01-01

perturbation method using stochastic differential equations. The joint statistical moments entering the perturbation solution are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vector and their first and second derivatives with respect......The paper deals with the first and second order statistical moments of the response of linear systems with random parameters subject to random excitation modelled as white-noise multiplied by an envelope function with random parameters. The method of analysis is basically a second order...... to the random parameters of the problem. Equations for partial derivatives are obtained from the partial differentiation of the equations of motion. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. General formulation is given...

A Stochastic Maximum Principle for a Stochastic Differential Game of a Mean-Field Type

Energy Technology Data Exchange (ETDEWEB)

Hosking, John Joseph Absalom, E-mail: j.j.a.hosking@cma.uio.no [University of Oslo, Centre of Mathematics for Applications (CMA) (Norway)

2012-12-15

We construct a stochastic maximum principle (SMP) which provides necessary conditions for the existence of Nash equilibria in a certain form of N-agent stochastic differential game (SDG) of a mean-field type. The information structure considered for the SDG is of a possible asymmetric and partial type. To prove our SMP we take an approach based on spike-variations and adjoint representation techniques, analogous to that of S. Peng (SIAM J. Control Optim. 28(4):966-979, 1990) in the optimal stochastic control context. In our proof we apply adjoint representation procedures at three points. The first-order adjoint processes are defined as solutions to certain mean-field backward stochastic differential equations, and second-order adjoint processes of a first type are defined as solutions to certain backward stochastic differential equations. Second-order adjoint processes of a second type are defined as solutions of certain backward stochastic equations of a type that we introduce in this paper, and which we term conditional mean-field backward stochastic differential equations. From the resulting representations, we show that the terms relating to these second-order adjoint processes of the second type are of an order such that they do not appear in our final SMP equations. A comparable situation exists in an article by R. Buckdahn, B. Djehiche, and J. Li (Appl. Math. Optim. 64(2):197-216, 2011) that constructs a SMP for a mean-field type optimal stochastic control problem; however, the approach we take of using these second-order adjoint processes of a second type to deal with the type of terms that we refer to as the second form of quadratic-type terms represents an alternative to a development, to our setting, of the approach used in their article for their analogous type of term.

A Stochastic Maximum Principle for a Stochastic Differential Game of a Mean-Field Type

International Nuclear Information System (INIS)

Hosking, John Joseph Absalom

2012-01-01

We construct a stochastic maximum principle (SMP) which provides necessary conditions for the existence of Nash equilibria in a certain form of N-agent stochastic differential game (SDG) of a mean-field type. The information structure considered for the SDG is of a possible asymmetric and partial type. To prove our SMP we take an approach based on spike-variations and adjoint representation techniques, analogous to that of S. Peng (SIAM J. Control Optim. 28(4):966–979, 1990) in the optimal stochastic control context. In our proof we apply adjoint representation procedures at three points. The first-order adjoint processes are defined as solutions to certain mean-field backward stochastic differential equations, and second-order adjoint processes of a first type are defined as solutions to certain backward stochastic differential equations. Second-order adjoint processes of a second type are defined as solutions of certain backward stochastic equations of a type that we introduce in this paper, and which we term conditional mean-field backward stochastic differential equations. From the resulting representations, we show that the terms relating to these second-order adjoint processes of the second type are of an order such that they do not appear in our final SMP equations. A comparable situation exists in an article by R. Buckdahn, B. Djehiche, and J. Li (Appl. Math. Optim. 64(2):197–216, 2011) that constructs a SMP for a mean-field type optimal stochastic control problem; however, the approach we take of using these second-order adjoint processes of a second type to deal with the type of terms that we refer to as the second form of quadratic-type terms represents an alternative to a development, to our setting, of the approach used in their article for their analogous type of term.

New Exact Solutions for the Wick-Type Stochastic Kudryashov–Sinelshchikov Equation

International Nuclear Information System (INIS)

Ray, S. Saha; Singh, S.

2017-01-01

In this article, exact solutions of Wick-type stochastic Kudryashov–Sinelshchikov equation have been obtained by using improved Sub-equation method. We have used Hermite transform for transforming the Wick-type stochastic Kudryashov–Sinelshchikov equation to deterministic partial differential equation. Also we have applied inverse Hermite transform for obtaining a set of stochastic solutions in the white noise space. (paper)

Differential equations driven by rough paths with jumps

Science.gov (United States)

Friz, Peter K.; Zhang, Huilin

2018-05-01

We develop the rough path counterpart of Itô stochastic integration and differential equations driven by general semimartingales. This significantly enlarges the classes of (Itô/forward) stochastic differential equations treatable with pathwise methods. A number of applications are discussed.

On Some Fractional Stochastic Integrodifferential Equations in Hilbert Space

Directory of Open Access Journals (Sweden)

Hamdy M. Ahmed

2009-01-01

Full Text Available We study a class of fractional stochastic integrodifferential equations considered in a real Hilbert space. The existence and uniqueness of the Mild solutions of the considered problem is also studied. We also give an application for stochastic integropartial differential equations of fractional order.

Stationary distributions of stochastic processes described by a linear neutral delay differential equation

International Nuclear Information System (INIS)

Frank, T D

2005-01-01

Stationary distributions of processes are derived that involve a time delay and are defined by a linear stochastic neutral delay differential equation. The distributions are Gaussian distributions. The variances of the Gaussian distributions are either monotonically increasing or decreasing functions of the time delays. The variances become infinite when fixed points of corresponding deterministic processes become unstable. (letter to the editor)

Lp Theory for Super-Parabolic Backward Stochastic Partial Differential Equations in the Whole Space

International Nuclear Information System (INIS)

Du Kai; Qiu, Jinniao; Tang Shanjian

2012-01-01

This paper is concerned with semi-linear backward stochastic partial differential equations (BSPDEs for short) of super-parabolic type. An L p -theory is given for the Cauchy problem of BSPDEs, separately for the case of p∈(1,2] and for the case of p∈(2,∞). A comparison theorem is also addressed.

Analytical determination of the bifurcation thresholds in stochastic differential equations with delayed feedback.

Science.gov (United States)

Gaudreault, Mathieu; Drolet, François; Viñals, Jorge

2010-11-01

Analytical expressions for pitchfork and Hopf bifurcation thresholds are given for a nonlinear stochastic differential delay equation with feedback. Our results assume that the delay time τ is small compared to other characteristic time scales, not a significant limitation close to the bifurcation line. A pitchfork bifurcation line is found, the location of which depends on the conditional average , where x(t) is the dynamical variable. This conditional probability incorporates the combined effect of fluctuation correlations and delayed feedback. We also find a Hopf bifurcation line which is obtained by a multiple scale expansion around the oscillatory solution near threshold. We solve the Fokker-Planck equation associated with the slowly varying amplitudes and use it to determine the threshold location. In both cases, the predicted bifurcation lines are in excellent agreement with a direct numerical integration of the governing equations. Contrary to the known case involving no delayed feedback, we show that the stochastic bifurcation lines are shifted relative to the deterministic limit and hence that the interaction between fluctuation correlations and delay affect the stability of the solutions of the model equation studied.

Workshop on quantum stochastic differential equations for the quantum simulation of physical systems

Science.gov (United States)

2016-09-22

that would be complimentary to the efforts at ARL. One the other hand, topological quantum field theories have a dual application to topological...Witten provided a path-integral definition of the Jones polynomial using a three-dimensional Chern-Simons quantum field theory (QFT) based on a non...topology, quantum field theory , quantum stochastic differential equations, quantum computing REPORT DOCUMENTATION PAGE 11. SPONSOR/MONITOR’S REPORT

Infinite time interval backward stochastic differential equations with continuous coefficients.

Science.gov (United States)

Zong, Zhaojun; Hu, Feng

2016-01-01

In this paper, we study the existence theorem for [Formula: see text] [Formula: see text] solutions to a class of 1-dimensional infinite time interval backward stochastic differential equations (BSDEs) under the conditions that the coefficients are continuous and have linear growths. We also obtain the existence of a minimal solution. Furthermore, we study the existence and uniqueness theorem for [Formula: see text] [Formula: see text] solutions of infinite time interval BSDEs with non-uniformly Lipschitz coefficients. It should be pointed out that the assumptions of this result is weaker than that of Theorem 3.1 in Zong (Turkish J Math 37:704-718, 2013).

A Posteriori Error Analysis of Stochastic Differential Equations Using Polynomial Chaos Expansions

KAUST Repository

Butler, T.; Dawson, C.; Wildey, T.

2011-01-01

We develop computable a posteriori error estimates for linear functionals of a solution to a general nonlinear stochastic differential equation with random model/source parameters. These error estimates are based on a variational analysis applied to stochastic Galerkin methods for forward and adjoint problems. The result is a representation for the error estimate as a polynomial in the random model/source parameter. The advantage of this method is that we use polynomial chaos representations for the forward and adjoint systems to cheaply produce error estimates by simple evaluation of a polynomial. By comparison, the typical method of producing such estimates requires repeated forward/adjoint solves for each new choice of random parameter. We present numerical examples showing that there is excellent agreement between these methods. © 2011 Society for Industrial and Applied Mathematics.

Development of a restricted state space stochastic differential equation model for bacterial growth in rich media

DEFF Research Database (Denmark)

Møller, Jan Kloppenborg; Philipsen, Kirsten Riber; Christiansen, Lasse Engbo

2012-01-01

In the present study, bacterial growth in a rich media is analysed in a Stochastic Differential Equation (SDE) framework. It is demonstrated that the SDE formulation and smoothened state estimates provide a systematic framework for data driven model improvements, using random walk hidden states...

Emergent user behavior on Twitter modelled by a stochastic differential equation.

Science.gov (United States)

Mollgaard, Anders; Mathiesen, Joachim

2015-01-01

Data from the social-media site, Twitter, is used to study the fluctuations in tweet rates of brand names. The tweet rates are the result of a strongly correlated user behavior, which leads to bursty collective dynamics with a characteristic 1/f noise. Here we use the aggregated "user interest" in a brand name to model collective human dynamics by a stochastic differential equation with multiplicative noise. The model is supported by a detailed analysis of the tweet rate fluctuations and it reproduces both the exact bursty dynamics found in the data and the 1/f noise.

Individualism in plant populations: using stochastic differential equations to model individual neighbourhood-dependent plant growth.

Science.gov (United States)

Lv, Qiming; Schneider, Manuel K; Pitchford, Jonathan W

2008-08-01

We study individual plant growth and size hierarchy formation in an experimental population of Arabidopsis thaliana, within an integrated analysis that explicitly accounts for size-dependent growth, size- and space-dependent competition, and environmental stochasticity. It is shown that a Gompertz-type stochastic differential equation (SDE) model, involving asymmetric competition kernels and a stochastic term which decreases with the logarithm of plant weight, efficiently describes individual plant growth, competition, and variability in the studied population. The model is evaluated within a Bayesian framework and compared to its deterministic counterpart, and to several simplified stochastic models, using distributional validation. We show that stochasticity is an important determinant of size hierarchy and that SDE models outperform the deterministic model if and only if structural components of competition (asymmetry; size- and space-dependence) are accounted for. Implications of these results are discussed in the context of plant ecology and in more general modelling situations.

Nonparametric estimation of stochastic differential equations with sparse Gaussian processes.

Science.gov (United States)

García, Constantino A; Otero, Abraham; Félix, Paulo; Presedo, Jesús; Márquez, David G

2017-08-01

The application of stochastic differential equations (SDEs) to the analysis of temporal data has attracted increasing attention, due to their ability to describe complex dynamics with physically interpretable equations. In this paper, we introduce a nonparametric method for estimating the drift and diffusion terms of SDEs from a densely observed discrete time series. The use of Gaussian processes as priors permits working directly in a function-space view and thus the inference takes place directly in this space. To cope with the computational complexity that requires the use of Gaussian processes, a sparse Gaussian process approximation is provided. This approximation permits the efficient computation of predictions for the drift and diffusion terms by using a distribution over a small subset of pseudosamples. The proposed method has been validated using both simulated data and real data from economy and paleoclimatology. The application of the method to real data demonstrates its ability to capture the behavior of complex systems.

Kramers-Moyal expansion for stochastic differential equations with single and multiple delays: Applications to financial physics and neurophysics

International Nuclear Information System (INIS)

Frank, T.D.

2007-01-01

We present a generalized Kramers-Moyal expansion for stochastic differential equations with single and multiple delays. In particular, we show that the delay Fokker-Planck equation derived earlier in the literature is a special case of the proposed Kramers-Moyal expansion. Applications for bond pricing and a self-inhibitory neuron model are discussed

Stochastic partial differential equations

CERN Document Server

Lototsky, Sergey V

2017-01-01

Taking readers with a basic knowledge of probability and real analysis to the frontiers of a very active research discipline, this textbook provides all the necessary background from functional analysis and the theory of PDEs. It covers the main types of equations (elliptic, hyperbolic and parabolic) and discusses different types of random forcing. The objective is to give the reader the necessary tools to understand the proofs of existing theorems about SPDEs (from other sources) and perhaps even to formulate and prove a few new ones. Most of the material could be covered in about 40 hours of lectures, as long as not too much time is spent on the general discussion of stochastic analysis in infinite dimensions. As the subject of SPDEs is currently making the transition from the research level to that of a graduate or even undergraduate course, the book attempts to present enough exercise material to fill potential exams and homework assignments. Exercises appear throughout and are usually directly connected ...

Stochastic equations for complex systems theoretical and computational topics

CERN Document Server

Bessaih, Hakima

2015-01-01

Mathematical analyses and computational predictions of the behavior of complex systems are needed to effectively deal with weather and climate predictions, for example, and the optimal design of technical processes. Given the random nature of such systems and the recognized relevance of randomness, the equations used to describe such systems usually need to involve stochastics.  The basic goal of this book is to introduce the mathematics and application of stochastic equations used for the modeling of complex systems. A first focus is on the introduction to different topics in mathematical analysis. A second focus is on the application of mathematical tools to the analysis of stochastic equations. A third focus is on the development and application of stochastic methods to simulate turbulent flows as seen in reality.  This book is primarily oriented towards mathematics and engineering PhD students, young and experienced researchers, and professionals working in the area of stochastic differential equations ...

PKPD model of interleukin-21 effects on thermoregulation in monkeys - Application and evaluation of stochastic differential equations

DEFF Research Database (Denmark)

Overgaard, Rune Viig; Holford, Nick; Rytved, K. A.

2007-01-01

Purpose To describe the pharmacodynamic effects of recombinant human interleukin-21 (IL-21) on core body temperature in cynomolgus monkeys using basic mechanisms of heat regulation. A major effort was devoted to compare the use of ordinary differential equations (ODEs) with stochastic differentia...

A Convergence Result for the Euler-Maruyama Method for a Simple Stochastic Differential Equation with Discontinuous Drift

DEFF Research Database (Denmark)

Simonsen, Maria; Schiøler, Henrik; Leth, John-Josef

2014-01-01

The Euler-Maruyama method is applied to a simple stochastic differential equation (SDE) with discontinuous drift. Convergence aspects are investigated in the case, where the Euler-Maruyama method is simulated in dyadic points. A strong rate of convergence is presented for the numerical simulations...

Evaluation of stochastic differential equation approximation of ion channel gating models.

Science.gov (United States)

Bruce, Ian C

2009-04-01

Fox and Lu derived an algorithm based on stochastic differential equations for approximating the kinetics of ion channel gating that is simpler and faster than "exact" algorithms for simulating Markov process models of channel gating. However, the approximation may not be sufficiently accurate to predict statistics of action potential generation in some cases. The objective of this study was to develop a framework for analyzing the inaccuracies and determining their origin. Simulations of a patch of membrane with voltage-gated sodium and potassium channels were performed using an exact algorithm for the kinetics of channel gating and the approximate algorithm of Fox & Lu. The Fox & Lu algorithm assumes that channel gating particle dynamics have a stochastic term that is uncorrelated, zero-mean Gaussian noise, whereas the results of this study demonstrate that in many cases the stochastic term in the Fox & Lu algorithm should be correlated and non-Gaussian noise with a non-zero mean. The results indicate that: (i) the source of the inaccuracy is that the Fox & Lu algorithm does not adequately describe the combined behavior of the multiple activation particles in each sodium and potassium channel, and (ii) the accuracy does not improve with increasing numbers of channels.

Incorporating prior knowledge induced from stochastic differential equations in the classification of stochastic observations.

Science.gov (United States)

Zollanvari, Amin; Dougherty, Edward R

2016-12-01

In classification, prior knowledge is incorporated in a Bayesian framework by assuming that the feature-label distribution belongs to an uncertainty class of feature-label distributions governed by a prior distribution. A posterior distribution is then derived from the prior and the sample data. An optimal Bayesian classifier (OBC) minimizes the expected misclassification error relative to the posterior distribution. From an application perspective, prior construction is critical. The prior distribution is formed by mapping a set of mathematical relations among the features and labels, the prior knowledge, into a distribution governing the probability mass across the uncertainty class. In this paper, we consider prior knowledge in the form of stochastic differential equations (SDEs). We consider a vector SDE in integral form involving a drift vector and dispersion matrix. Having constructed the prior, we develop the optimal Bayesian classifier between two models and examine, via synthetic experiments, the effects of uncertainty in the drift vector and dispersion matrix. We apply the theory to a set of SDEs for the purpose of differentiating the evolutionary history between two species.

Structure and properties of Hughston's stochastic extension of the Schroedinger equation

International Nuclear Information System (INIS)

Adler, Stephen L.; Horwitz, Lawrence P.

2000-01-01

Hughston has recently proposed a stochastic extension of the Schroedinger equation, expressed as a stochastic differential equation on projective Hilbert space. We derive new projective Hilbert space identities, which we use to give a general proof that Hughston's equation leads to state vector collapse to energy eigenstates, with collapse probabilities given by the quantum mechanical probabilities computed from the initial state. We discuss the relation of Hughston's equation to earlier work on norm-preserving stochastic equations, and show that Hughston's equation can be written as a manifestly unitary stochastic evolution equation for the pure state density matrix. We discuss the behavior of systems constructed as direct products of independent subsystems, and briefly address the question of whether an energy-based approach, such as Hughston's, suffices to give an objective interpretation of the measurement process in quantum mechanics. (c) 2000 American Institute of Physics

Stochastic Calculus and Differential Equations for Physics and Finance

Science.gov (United States)

McCauley, Joseph L.

2013-02-01

1. Random variables and probability distributions; 2. Martingales, Markov, and nonstationarity; 3. Stochastic calculus; 4. Ito processes and Fokker-Planck equations; 5. Selfsimilar Ito processes; 6. Fractional Brownian motion; 7. Kolmogorov's PDEs and Chapman-Kolmogorov; 8. Non Markov Ito processes; 9. Black-Scholes, martingales, and Feynman-Katz; 10. Stochastic calculus with martingales; 11. Statistical physics and finance, a brief history of both; 12. Introduction to new financial economics; 13. Statistical ensembles and time series analysis; 14. Econometrics; 15. Semimartingales; References; Index.

Tuning of Controller for Type 1 Diabetes Treatment with Stochastic Differential Equations

DEFF Research Database (Denmark)

Duun-Henriksen, Anne Katrine; Boiroux, Dimitri; Schmidt, Signe

2012-01-01

due to the noise corrupted observations from the CGM. In this paper we present a method to estimate the optimal Kalman gain in the controller based on stochastic differential equation modeling. With this model type we could estimate the process noise and observation noise separately based on data from......People with type 1 diabetes need several insulin injections every day to keep their blood glucose level in the normal range and thereby avoiding the acute and long term complications of diabetes. One of the recent treatments consists of a pump injecting insulin into the subcutaneous layer combined...

Existence, uniqueness, and stability of stochastic neutral functional differential equations of Sobolev-type

Energy Technology Data Exchange (ETDEWEB)

Yang, Xuetao; Zhu, Quanxin, E-mail: zqx22@126.com [School of Mathematical Sciences and Institute of Mathematics, Nanjing Normal University, Nanjing 210023, Jiangsu (China)

2015-12-15

In this paper, we are mainly concerned with a class of stochastic neutral functional differential equations of Sobolev-type with Poisson jumps. Under two different sets of conditions, we establish the existence of the mild solution by applying the Leray-Schauder alternative theory and the Sadakovskii’s fixed point theorem, respectively. Furthermore, we use the Bihari’s inequality to prove the Osgood type uniqueness. Also, the mean square exponential stability is investigated by applying the Gronwall inequality. Finally, two examples are given to illustrate the theory results.

Analyzing a stochastic time series obeying a second-order differential equation.

Science.gov (United States)

Lehle, B; Peinke, J

2015-06-01

The stochastic properties of a Langevin-type Markov process can be extracted from a given time series by a Markov analysis. Also processes that obey a stochastically forced second-order differential equation can be analyzed this way by employing a particular embedding approach: To obtain a Markovian process in 2N dimensions from a non-Markovian signal in N dimensions, the system is described in a phase space that is extended by the temporal derivative of the signal. For a discrete time series, however, this derivative can only be calculated by a differencing scheme, which introduces an error. If the effects of this error are not accounted for, this leads to systematic errors in the estimation of the drift and diffusion functions of the process. In this paper we will analyze these errors and we will propose an approach that correctly accounts for them. This approach allows an accurate parameter estimation and, additionally, is able to cope with weak measurement noise, which may be superimposed to a given time series.

T-Stability of the Heun Method and Balanced Method for Solving Stochastic Differential Delay Equations

Directory of Open Access Journals (Sweden)

Xiaolin Zhu

2014-01-01

Full Text Available This paper studies the T-stability of the Heun method and balanced method for solving stochastic differential delay equations (SDDEs. Two T-stable conditions of the Heun method are obtained for two kinds of linear SDDEs. Moreover, two conditions under which the balanced method is T-stable are obtained for two kinds of linear SDDEs. Some numerical examples verify the theoretical results proposed.

Continuity in a pathwise sense with respect to the coefficients of solutions of stochastic differential equations

DEFF Research Database (Denmark)

Knudsen, Thomas Skov

1997-01-01

For stochastic differential equations (SDEs) of the form dX(t) = b(X)(t)) dt + sigma(X(t))dW(t) where b and sigma are Lipschitz continuous, it is shown that if we consider a fixed sigma is an element of C-5, bounded and with bounded derivatives, the random field of solutions is pathwise locally...

Efficient Galerkin solution of stochastic fractional differential equations using second kind Chebyshev wavelets

Directory of Open Access Journals (Sweden)

Fakhrodin Mohammadi

2017-10-01

Full Text Available ‎Stochastic fractional differential equations (SFDEs have been used for modeling many physical problems in the fields of turbulance‎, ‎heterogeneous‎, ‎flows and matrials‎, ‎viscoelasticity and electromagnetic theory‎. ‎In this paper‎, ‎an‎ efficient wavelet Galerkin method based on the second kind Chebyshev wavelets are proposed for approximate solution of SFDEs‎. ‎In ‎this ‎app‎roach‎‎, ‎o‎perational matrices of the second kind Chebyshev wavelets ‎are used ‎for reducing SFDEs to a linear system of algebraic equations that can be solved easily‎. ‎C‎onvergence and error analysis of the proposed method is ‎considered‎.‎ ‎Some numerical examples are performed to confirm the applicability and efficiency of the proposed method‎.

Approaches for modeling within subject variability in pharmacometric count data analysis: dynamic inter-occasion variability and stochastic differential equations.

Science.gov (United States)

Deng, Chenhui; Plan, Elodie L; Karlsson, Mats O

2016-06-01

Parameter variation in pharmacometric analysis studies can be characterized as within subject parameter variability (WSV) in pharmacometric models. WSV has previously been successfully modeled using inter-occasion variability (IOV), but also stochastic differential equations (SDEs). In this study, two approaches, dynamic inter-occasion variability (dIOV) and adapted stochastic differential equations, were proposed to investigate WSV in pharmacometric count data analysis. These approaches were applied to published count models for seizure counts and Likert pain scores. Both approaches improved the model fits significantly. In addition, stochastic simulation and estimation were used to explore further the capability of the two approaches to diagnose and improve models where existing WSV is not recognized. The results of simulations confirmed the gain in introducing WSV as dIOV and SDEs when parameters vary randomly over time. Further, the approaches were also informative as diagnostics of model misspecification, when parameters changed systematically over time but this was not recognized in the structural model. The proposed approaches in this study offer strategies to characterize WSV and are not restricted to count data.

Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation

KAUST Repository

Mélykúti, Bence; Burrage, Kevin; Zygalakis, Konstantinos C.

2010-01-01

The Chemical Langevin Equation (CLE), which is a stochastic differential equation driven by a multidimensional Wiener process, acts as a bridge between the discrete stochastic simulation algorithm and the deterministic reaction rate equation when

The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models

International Nuclear Information System (INIS)

Rupšys, P.

2015-01-01

A system of stochastic differential equations (SDE) with mixed-effects parameters and multivariate normal copula density function were used to develop tree height model for Scots pine trees in Lithuania. A two-step maximum likelihood parameter estimation method is used and computational guidelines are given. After fitting the conditional probability density functions to outside bark diameter at breast height, and total tree height, a bivariate normal copula distribution model was constructed. Predictions from the mixed-effects parameters SDE tree height model calculated during this research were compared to the regression tree height equations. The results are implemented in the symbolic computational language MAPLE

The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models

Energy Technology Data Exchange (ETDEWEB)

Rupšys, P. [Aleksandras Stulginskis University, Studenų g. 11, Akademija, Kaunas district, LT – 53361 Lithuania (Lithuania)

2015-10-28

A system of stochastic differential equations (SDE) with mixed-effects parameters and multivariate normal copula density function were used to develop tree height model for Scots pine trees in Lithuania. A two-step maximum likelihood parameter estimation method is used and computational guidelines are given. After fitting the conditional probability density functions to outside bark diameter at breast height, and total tree height, a bivariate normal copula distribution model was constructed. Predictions from the mixed-effects parameters SDE tree height model calculated during this research were compared to the regression tree height equations. The results are implemented in the symbolic computational language MAPLE.

Numerical simulation of stochastic point kinetic equation in the dynamical system of nuclear reactor

International Nuclear Information System (INIS)

Saha Ray, S.

2012-01-01

Highlights: ► In this paper stochastic neutron point kinetic equations have been analyzed. ► Euler–Maruyama method and Strong Taylor 1.5 order method have been discussed. ► These methods are applied for the solution of stochastic point kinetic equations. ► Comparison between the results of these methods and others are presented in tables. ► Graphs for neutron and precursor sample paths are also presented. -- Abstract: In the present paper, the numerical approximation methods, applied to efficiently calculate the solution for stochastic point kinetic equations () in nuclear reactor dynamics, are investigated. A system of Itô stochastic differential equations has been analyzed to model the neutron density and the delayed neutron precursors in a point nuclear reactor. The resulting system of Itô stochastic differential equations are solved over each time-step size. The methods are verified by considering different initial conditions, experimental data and over constant reactivities. The computational results indicate that the methods are simple and suitable for solving stochastic point kinetic equations. In this article, a numerical investigation is made in order to observe the random oscillations in neutron and precursor population dynamics in subcritical and critical reactors.

Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations.

Science.gov (United States)

Tornøe, Christoffer W; Overgaard, Rune V; Agersø, Henrik; Nielsen, Henrik A; Madsen, Henrik; Jonsson, E Niclas

2005-08-01

The objective of the present analysis was to explore the use of stochastic differential equations (SDEs) in population pharmacokinetic/pharmacodynamic (PK/PD) modeling. The intra-individual variability in nonlinear mixed-effects models based on SDEs is decomposed into two types of noise: a measurement and a system noise term. The measurement noise represents uncorrelated error due to, for example, assay error while the system noise accounts for structural misspecifications, approximations of the dynamical model, and true random physiological fluctuations. Since the system noise accounts for model misspecifications, the SDEs provide a diagnostic tool for model appropriateness. The focus of the article is on the implementation of the Extended Kalman Filter (EKF) in NONMEM for parameter estimation in SDE models. Various applications of SDEs in population PK/PD modeling are illustrated through a systematic model development example using clinical PK data of the gonadotropin releasing hormone (GnRH) antagonist degarelix. The dynamic noise estimates were used to track variations in model parameters and systematically build an absorption model for subcutaneously administered degarelix. The EKF-based algorithm was successfully implemented in NONMEM for parameter estimation in population PK/PD models described by systems of SDEs. The example indicated that it was possible to pinpoint structural model deficiencies, and that valuable information may be obtained by tracking unexplained variations in parameters.

On parametric domain for asymptotic stability with probability one of zero solution of linear Ito stochastic differential equations

International Nuclear Information System (INIS)

Phan Thanh An; Phan Le Na; Ngo Quoc Chung

2004-05-01

We describe a practical implementation for finding parametric domain for asymptotic stability with probability one of zero solution of linear Ito stochastic differential equations based on Korenevskij and Mitropolskij's sufficient condition and our sufficient conditions. Numerical results show that all of these sufficient conditions are crucial in the implementation. (author)

Stochastic Optimal Prediction with Application to Averaged Euler Equations

Energy Technology Data Exchange (ETDEWEB)

Bell, John [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Chorin, Alexandre J. [Univ. of California, Berkeley, CA (United States); Crutchfield, William [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

2017-04-24

Optimal prediction (OP) methods compensate for a lack of resolution in the numerical solution of complex problems through the use of an invariant measure as a prior measure in the Bayesian sense. In first-order OP, unresolved information is approximated by its conditional expectation with respect to the invariant measure. In higher-order OP, unresolved information is approximated by a stochastic estimator, leading to a system of random or stochastic differential equations. We explain the ideas through a simple example, and then apply them to the solution of Averaged Euler equations in two space dimensions.

Effective computation of stochastic protein kinetic equation by reducing stiffness via variable transformation

Energy Technology Data Exchange (ETDEWEB)

Wang, Lijin, E-mail: ljwang@ucas.ac.cn [School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049 (China)

2016-06-08

The stochastic protein kinetic equations can be stiff for certain parameters, which makes their numerical simulation rely on very small time step sizes, resulting in large computational cost and accumulated round-off errors. For such situation, we provide a method of reducing stiffness of the stochastic protein kinetic equation by means of a kind of variable transformation. Theoretical and numerical analysis show effectiveness of this method. Its generalization to a more general class of stochastic differential equation models is also discussed.

Probabilistic Forecasts of Wind Power Generation by Stochastic Differential Equation Models

DEFF Research Database (Denmark)

Møller, Jan Kloppenborg; Zugno, Marco; Madsen, Henrik

2016-01-01

The increasing penetration of wind power has resulted in larger shares of volatile sources of supply in power systems worldwide. In order to operate such systems efficiently, methods for reliable probabilistic forecasts of future wind power production are essential. It is well known...... that the conditional density of wind power production is highly dependent on the level of predicted wind power and prediction horizon. This paper describes a new approach for wind power forecasting based on logistic-type stochastic differential equations (SDEs). The SDE formulation allows us to calculate both state......-dependent conditional uncertainties as well as correlation structures. Model estimation is performed by maximizing the likelihood of a multidimensional random vector while accounting for the correlation structure defined by the SDE formulation. We use non-parametric modelling to explore conditional correlation...

Derivation of exact master equation with stochastic description: dissipative harmonic oscillator.

Science.gov (United States)

Li, Haifeng; Shao, Jiushu; Wang, Shikuan

2011-11-01

A systematic procedure for deriving the master equation of a dissipative system is reported in the framework of stochastic description. For the Caldeira-Leggett model of the harmonic-oscillator bath, a detailed and elementary derivation of the bath-induced stochastic field is presented. The dynamics of the system is thereby fully described by a stochastic differential equation, and the desired master equation would be acquired with statistical averaging. It is shown that the existence of a closed-form master equation depends on the specificity of the system as well as the feature of the dissipation characterized by the spectral density function. For a dissipative harmonic oscillator it is observed that the correlation between the stochastic field due to the bath and the system can be decoupled, and the master equation naturally results. Such an equation possesses the Lindblad form in which time-dependent coefficients are determined by a set of integral equations. It is proved that the obtained master equation is equivalent to the well-known Hu-Paz-Zhang equation based on the path-integral technique. The procedure is also used to obtain the master equation of a dissipative harmonic oscillator in time-dependent fields.

Modeling delay in genetic networks: from delay birth-death processes to delay stochastic differential equations.

Science.gov (United States)

Gupta, Chinmaya; López, José Manuel; Azencott, Robert; Bennett, Matthew R; Josić, Krešimir; Ott, William

2014-05-28

Delay is an important and ubiquitous aspect of many biochemical processes. For example, delay plays a central role in the dynamics of genetic regulatory networks as it stems from the sequential assembly of first mRNA and then protein. Genetic regulatory networks are therefore frequently modeled as stochastic birth-death processes with delay. Here, we examine the relationship between delay birth-death processes and their appropriate approximating delay chemical Langevin equations. We prove a quantitative bound on the error between the pathwise realizations of these two processes. Our results hold for both fixed delay and distributed delay. Simulations demonstrate that the delay chemical Langevin approximation is accurate even at moderate system sizes. It captures dynamical features such as the oscillatory behavior in negative feedback circuits, cross-correlations between nodes in a network, and spatial and temporal information in two commonly studied motifs of metastability in biochemical systems. Overall, these results provide a foundation for using delay stochastic differential equations to approximate the dynamics of birth-death processes with delay.

Modeling delay in genetic networks: From delay birth-death processes to delay stochastic differential equations

Energy Technology Data Exchange (ETDEWEB)

Gupta, Chinmaya; López, José Manuel; Azencott, Robert; Ott, William [Department of Mathematics, University of Houston, Houston, Texas 77004 (United States); Bennett, Matthew R. [Department of Biochemistry and Cell Biology, Rice University, Houston, Texas 77204, USA and Institute of Biosciences and Bioengineering, Rice University, Houston, Texas 77005 (United States); Josić, Krešimir [Department of Mathematics, University of Houston, Houston, Texas 77004 (United States); Department of Biology and Biochemistry, University of Houston, Houston, Texas 77204 (United States)

2014-05-28

Delay is an important and ubiquitous aspect of many biochemical processes. For example, delay plays a central role in the dynamics of genetic regulatory networks as it stems from the sequential assembly of first mRNA and then protein. Genetic regulatory networks are therefore frequently modeled as stochastic birth-death processes with delay. Here, we examine the relationship between delay birth-death processes and their appropriate approximating delay chemical Langevin equations. We prove a quantitative bound on the error between the pathwise realizations of these two processes. Our results hold for both fixed delay and distributed delay. Simulations demonstrate that the delay chemical Langevin approximation is accurate even at moderate system sizes. It captures dynamical features such as the oscillatory behavior in negative feedback circuits, cross-correlations between nodes in a network, and spatial and temporal information in two commonly studied motifs of metastability in biochemical systems. Overall, these results provide a foundation for using delay stochastic differential equations to approximate the dynamics of birth-death processes with delay.

Modeling delay in genetic networks: From delay birth-death processes to delay stochastic differential equations

International Nuclear Information System (INIS)

Gupta, Chinmaya; López, José Manuel; Azencott, Robert; Ott, William; Bennett, Matthew R.; Josić, Krešimir

2014-01-01

Delay is an important and ubiquitous aspect of many biochemical processes. For example, delay plays a central role in the dynamics of genetic regulatory networks as it stems from the sequential assembly of first mRNA and then protein. Genetic regulatory networks are therefore frequently modeled as stochastic birth-death processes with delay. Here, we examine the relationship between delay birth-death processes and their appropriate approximating delay chemical Langevin equations. We prove a quantitative bound on the error between the pathwise realizations of these two processes. Our results hold for both fixed delay and distributed delay. Simulations demonstrate that the delay chemical Langevin approximation is accurate even at moderate system sizes. It captures dynamical features such as the oscillatory behavior in negative feedback circuits, cross-correlations between nodes in a network, and spatial and temporal information in two commonly studied motifs of metastability in biochemical systems. Overall, these results provide a foundation for using delay stochastic differential equations to approximate the dynamics of birth-death processes with delay

A model and variance reduction method for computing statistical outputs of stochastic elliptic partial differential equations

International Nuclear Information System (INIS)

Vidal-Codina, F.; Nguyen, N.C.; Giles, M.B.; Peraire, J.

2015-01-01

We present a model and variance reduction method for the fast and reliable computation of statistical outputs of stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the hybridizable discontinuous Galerkin (HDG) discretization of elliptic partial differential equations (PDEs), which allows us to obtain high-order accurate solutions of the governing PDE; (2) the reduced basis method for a new HDG discretization of the underlying PDE to enable real-time solution of the parameterized PDE in the presence of stochastic parameters; and (3) a multilevel variance reduction method that exploits the statistical correlation among the different reduced basis approximations and the high-fidelity HDG discretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the reduced basis approximations. Furthermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the reduced basis approximations and the sizes of Monte Carlo samples to achieve a given error tolerance. We provide numerical examples to demonstrate the performance of the proposed method

Weak-noise limit of a piecewise-smooth stochastic differential equation.

Science.gov (United States)

Chen, Yaming; Baule, Adrian; Touchette, Hugo; Just, Wolfram

2013-11-01

We investigate the validity and accuracy of weak-noise (saddle-point or instanton) approximations for piecewise-smooth stochastic differential equations (SDEs), taking as an illustrative example a piecewise-constant SDE, which serves as a simple model of Brownian motion with solid friction. For this model, we show that the weak-noise approximation of the path integral correctly reproduces the known propagator of the SDE at lowest order in the noise power, as well as the main features of the exact propagator with higher-order corrections, provided the singularity of the path integral associated with the nonsmooth SDE is treated with some heuristics. We also show that, as in the case of smooth SDEs, the deterministic paths of the noiseless system correctly describe the behavior of the nonsmooth SDE in the low-noise limit. Finally, we consider a smooth regularization of the piecewise-constant SDE and study to what extent this regularization can rectify some of the problems encountered when dealing with discontinuous drifts and singularities in SDEs.

Modelling ventricular fibrillation coarseness during cardiopulmonary resuscitation by mixed effects stochastic differential equations.

Science.gov (United States)

Gundersen, Kenneth; Kvaløy, Jan Terje; Eftestøl, Trygve; Kramer-Johansen, Jo

2015-10-15

For patients undergoing cardiopulmonary resuscitation (CPR) and being in a shockable rhythm, the coarseness of the electrocardiogram (ECG) signal is an indicator of the state of the patient. In the current work, we show how mixed effects stochastic differential equations (SDE) models, commonly used in pharmacokinetic and pharmacodynamic modelling, can be used to model the relationship between CPR quality measurements and ECG coarseness. This is a novel application of mixed effects SDE models to a setting quite different from previous applications of such models and where using such models nicely solves many of the challenges involved in analysing the available data. Copyright © 2015 John Wiley & Sons, Ltd.

Differential form representation of stochastic electromagnetic fields

Science.gov (United States)

Haider, Michael; Russer, Johannes A.

2017-09-01

In this work, we revisit the theory of stochastic electromagnetic fields using exterior differential forms. We present a short overview as well as a brief introduction to the application of differential forms in electromagnetic theory. Within the framework of exterior calculus we derive equations for the second order moments, describing stochastic electromagnetic fields. Since the resulting objects are continuous quantities in space, a discretization scheme based on the Method of Moments (MoM) is introduced for numerical treatment. The MoM is applied in such a way, that the notation of exterior calculus is maintained while we still arrive at the same set of algebraic equations as obtained for the case of formulating the theory using the traditional notation of vector calculus. We conclude with an analytic calculation of the radiated electric field of two Hertzian dipole, excited by uncorrelated random currents.

Modeling real-time balancing power demands in wind power systems using stochastic differential equations

International Nuclear Information System (INIS)

Olsson, Magnus; Perninge, Magnus; Soeder, Lennart

2010-01-01

The inclusion of wind power into power systems has a significant impact on the demand for real-time balancing power due to the stochastic nature of wind power production. The overall aim of this paper is to present probabilistic models of the impact of large-scale integration of wind power on the continuous demand in MW for real-time balancing power. This is important not only for system operators, but also for producers and consumers since they in most systems through various market solutions provide balancing power. Since there can occur situations where the wind power variations cancel out other types of deviations in the system, models on an hourly basis are not sufficient. Therefore the developed model is in continuous time and is based on stochastic differential equations (SDE). The model can be used within an analytical framework or in Monte Carlo simulations. (author)

Jacobian elliptic function expansion solutions of nonlinear stochastic equations

International Nuclear Information System (INIS)

Wei Caimin; Xia Zunquan; Tian Naishuo

2005-01-01

Jacobian elliptic function expansion method is extended and applied to construct the exact solutions of the nonlinear Wick-type stochastic partial differential equations (SPDEs) and some new exact solutions are obtained via this method and Hermite transformation

A New Control Paradigm for Stochastic Differential Equations

Science.gov (United States)

Schmid, Matthias J. A.

This study presents a novel comprehensive approach to the control of dynamic systems under uncertainty governed by stochastic differential equations (SDEs). Large Deviations (LD) techniques are employed to arrive at a control law for a large class of nonlinear systems minimizing sample path deviations. Thereby, a paradigm shift is suggested from point-in-time to sample path statistics on function spaces. A suitable formal control framework which leverages embedded Freidlin-Wentzell theory is proposed and described in detail. This includes the precise definition of the control objective and comprises an accurate discussion of the adaptation of the Freidlin-Wentzell theorem to the particular situation. The new control design is enabled by the transformation of an ill-posed control objective into a well-conditioned sequential optimization problem. A direct numerical solution process is presented using quadratic programming, but the emphasis is on the development of a closed-form expression reflecting the asymptotic deviation probability of a particular nominal path. This is identified as the key factor in the success of the new paradigm. An approach employing the second variation and the differential curvature of the effective action is suggested for small deviation channels leading to the Jacobi field of the rate function and the subsequently introduced Jacobi field performance measure. This closed-form solution is utilized in combination with the supplied parametrization of the objective space. For the first time, this allows for an LD based control design applicable to a large class of nonlinear systems. Thus, Minimum Large Deviations (MLD) control is effectively established in a comprehensive structured framework. The construction of the new paradigm is completed by an optimality proof for the Jacobi field performance measure, an interpretive discussion, and a suggestion for efficient implementation. The potential of the new approach is exhibited by its extension

The solution of the neutron point kinetics equation with stochastic extension: an analysis of two moments

Energy Technology Data Exchange (ETDEWEB)

Silva, Milena Wollmann da; Vilhena, Marco Tullio M.B.; Bodmann, Bardo Ernst J.; Vasques, Richard, E-mail: milena.wollmann@ufrgs.br, E-mail: vilhena@mat.ufrgs.br, E-mail: bardobodmann@ufrgs.br, E-mail: richard.vasques@fulbrightmail.org [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica

2015-07-01

The neutron point kinetics equation, which models the time-dependent behavior of nuclear reactors, is often used to understand the dynamics of nuclear reactor operations. It consists of a system of coupled differential equations that models the interaction between (i) the neutron population; and (II) the concentration of the delayed neutron precursors, which are radioactive isotopes formed in the fission process that decay through neutron emission. These equations are deterministic in nature, and therefore can provide only average values of the modeled populations. However, the actual dynamical process is stochastic: the neutron density and the delayed neutron precursor concentrations vary randomly with time. To address this stochastic behavior, Hayes and Allen have generalized the standard deterministic point kinetics equation. They derived a system of stochastic differential equations that can accurately model the random behavior of the neutron density and the precursor concentrations in a point reactor. Due to the stiffness of these equations, this system was numerically implemented using a stochastic piecewise constant approximation method (Stochastic PCA). Here, we present a study of the influence of stochastic fluctuations on the results of the neutron point kinetics equation. We reproduce the stochastic formulation introduced by Hayes and Allen and compute Monte Carlo numerical results for examples with constant and time-dependent reactivity, comparing these results with stochastic and deterministic methods found in the literature. Moreover, we introduce a modified version of the stochastic method to obtain a non-stiff solution, analogue to a previously derived deterministic approach. (author)

The solution of the neutron point kinetics equation with stochastic extension: an analysis of two moments

International Nuclear Information System (INIS)

Silva, Milena Wollmann da; Vilhena, Marco Tullio M.B.; Bodmann, Bardo Ernst J.; Vasques, Richard

2015-01-01

The neutron point kinetics equation, which models the time-dependent behavior of nuclear reactors, is often used to understand the dynamics of nuclear reactor operations. It consists of a system of coupled differential equations that models the interaction between (i) the neutron population; and (II) the concentration of the delayed neutron precursors, which are radioactive isotopes formed in the fission process that decay through neutron emission. These equations are deterministic in nature, and therefore can provide only average values of the modeled populations. However, the actual dynamical process is stochastic: the neutron density and the delayed neutron precursor concentrations vary randomly with time. To address this stochastic behavior, Hayes and Allen have generalized the standard deterministic point kinetics equation. They derived a system of stochastic differential equations that can accurately model the random behavior of the neutron density and the precursor concentrations in a point reactor. Due to the stiffness of these equations, this system was numerically implemented using a stochastic piecewise constant approximation method (Stochastic PCA). Here, we present a study of the influence of stochastic fluctuations on the results of the neutron point kinetics equation. We reproduce the stochastic formulation introduced by Hayes and Allen and compute Monte Carlo numerical results for examples with constant and time-dependent reactivity, comparing these results with stochastic and deterministic methods found in the literature. Moreover, we introduce a modified version of the stochastic method to obtain a non-stiff solution, analogue to a previously derived deterministic approach. (author)

Stochastic Landau equation with time-dependent drift

International Nuclear Information System (INIS)

Swift, J.B.; Hohenberg, P.C.; Ahlers, G.

1991-01-01

The stochastic differential equation τ 0 ∂ tA =ε(t)A-g 3 A 3 +bar f(t), where bar f(t) is Gaussian white noise, is studied for arbitrary time dependence of ε(t). In particular, cases are considered where ε(t) goes through the bifurcation of the deterministic system, which occurs at ε=0. In the limit of weak noise an approximate analytic expression generalizing earlier work of Suzuki [Phys. Lett. A 67, 339 (1978); Prog. Theor. Phys. (Kyoto) Suppl. 64, 402 (1978)] is obtained for the time-dependent distribution function P(A,t). The results compare favorably with a numerical simulation of the stochastic equation for the case of a linear ramp (both increasing and decreasing) and for a periodic time dependence of ε(t). The procedure can be generalized to an arbitrary deterministic part ∂ tA =D(A,t)+bar f(t), but the deterministic equation may then have to be solved numerically

Differential form representation of stochastic electromagnetic fields

Directory of Open Access Journals (Sweden)

M. Haider

2017-09-01

Full Text Available In this work, we revisit the theory of stochastic electromagnetic fields using exterior differential forms. We present a short overview as well as a brief introduction to the application of differential forms in electromagnetic theory. Within the framework of exterior calculus we derive equations for the second order moments, describing stochastic electromagnetic fields. Since the resulting objects are continuous quantities in space, a discretization scheme based on the Method of Moments (MoM is introduced for numerical treatment. The MoM is applied in such a way, that the notation of exterior calculus is maintained while we still arrive at the same set of algebraic equations as obtained for the case of formulating the theory using the traditional notation of vector calculus. We conclude with an analytic calculation of the radiated electric field of two Hertzian dipole, excited by uncorrelated random currents.

A stochastic differential equation model of diurnal cortisol patterns

Science.gov (United States)

Brown, E. N.; Meehan, P. M.; Dempster, A. P.

2001-01-01

Circadian modulation of episodic bursts is recognized as the normal physiological pattern of diurnal variation in plasma cortisol levels. The primary physiological factors underlying these diurnal patterns are the ultradian timing of secretory events, circadian modulation of the amplitude of secretory events, infusion of the hormone from the adrenal gland into the plasma, and clearance of the hormone from the plasma by the liver. Each measured plasma cortisol level has an error arising from the cortisol immunoassay. We demonstrate that all of these three physiological principles can be succinctly summarized in a single stochastic differential equation plus measurement error model and show that physiologically consistent ranges of the model parameters can be determined from published reports. We summarize the model parameters in terms of the multivariate Gaussian probability density and establish the plausibility of the model with a series of simulation studies. Our framework makes possible a sensitivity analysis in which all model parameters are allowed to vary simultaneously. The model offers an approach for simultaneously representing cortisol's ultradian, circadian, and kinetic properties. Our modeling paradigm provides a framework for simulation studies and data analysis that should be readily adaptable to the analysis of other endocrine hormone systems.

Using random walk in models specified by stochastic differential equations to determine the best expression for the bacterial growth rate

DEFF Research Database (Denmark)

method allows us to develop a new expression for the growth rate. The method is based on the stochastic continuous-discrete time state-space model, with a continuous-time state equation (a stochastic differential equation, SDE) combined with a discrete-time measurement equation. In our study the SDE...... described by Kristensen et. al [2]. The resulting time series allows us graphically to inspect the functional dependence of the growth rate on the substrate content. From the method described above we find three new plausible expressions for μ(S). Therefore we apply the likelihood-ratio test to compare...... for the Monod expression. Thus, the method was applied to successfully determine a significant better expression for the substrate dependent growth expression, and we find the method generally applicable for model development. References [1] Kristensen NR, Madsen H, Jørgensen, SB (2004) A method for systematic...

Model Identification Using Stochastic Differential Equation Grey-Box Models in Diabetes

DEFF Research Database (Denmark)

Duun-Henriksen, Anne Katrine; Schmidt, Signe; Røge, Rikke Meldgaard

2013-01-01

are uncorrelated and provides the possibility to pinpoint model deficiencies. METHODS: An identifiable model of the glucoregulatory system in a type 1 diabetes mellitus (T1DM) patient is used as the basis for development of a stochastic-differential-equation-based grey-box model (SDE-GB). The parameters...... in a significant improvement in the prediction and uncorrelated errors. Tracking of the "peak time of meal absorption" parameter showed that the absorption rate varied according to meal type. CONCLUSION: This study shows the potential of using SDE-GBs in diabetes modeling. Improved model predictions were obtained...... are estimated on clinical data from four T1DM patients. The optimal SDE-GB is determined from likelihood-ratio tests. Finally, parameter tracking is used to track the variation in the "time to peak of meal response" parameter. RESULTS: We found that the transformation of the ODE model into an SDE-GB resulted...

Mean-Square Convergence of Drift-Implicit One-Step Methods for Neutral Stochastic Delay Differential Equations with Jump Diffusion

Directory of Open Access Journals (Sweden)

Lin Hu

2011-01-01

Full Text Available A class of drift-implicit one-step schemes are proposed for the neutral stochastic delay differential equations (NSDDEs driven by Poisson processes. A general framework for mean-square convergence of the methods is provided. It is shown that under certain conditions global error estimates for a method can be inferred from estimates on its local error. The applicability of the mean-square convergence theory is illustrated by the stochastic θ-methods and the balanced implicit methods. It is derived from Theorem 3.1 that the order of the mean-square convergence of both of them for NSDDEs with jumps is 1/2. Numerical experiments illustrate the theoretical results. It is worth noting that the results of mean-square convergence of the stochastic θ-methods and the balanced implicit methods are also new.

A stochastic solution of the advective transport equation with uncertainty

International Nuclear Information System (INIS)

Williams, M.M.R.

1991-01-01

A model has been developed for calculating the transport of water-borne radionuclides through layers of porous materials, such as rock or clay. The model is based upon a purely advective transport equation, in which the fluid velocity is a random variable, thereby simulating dispersion in a more realistic manner than the ad hoc introduction of a dispersivity. In addition to a random velocity field, which is an observable physical phenomenon, allowance is made for uncertainty in our knowledge of the parameters which enter the equation, e.g. the retardation coefficient. This too, is assumed to be a random variable and contributes to the stochasticity of the resulting partial differential equation of transport. The stochastic differential equation can be solved analytically and then ensemble averages taken over the associated probability distribution of velocity and retardation coefficient. A method based upon a novel form of the central limit theorem of statistics is employed to obtain tractable solutions of a system consisting of many serial legs of varying properties. One interesting conclusion is that the total flux out of a medium is significantly underestimated by using the deterministic solution with an average transit time compared with that from the stochastically averaged solution. The theory is illustrated numerically for a number of physically relevant cases. (author) 8 figs., 4 tabs., 7 refs

CRPropa 3.1—a low energy extension based on stochastic differential equations

Energy Technology Data Exchange (ETDEWEB)

Merten, Lukas; Tjus, Julia Becker; Eichmann, Björn [Theoretische Physik IV: Plasma-Astroteilchenphysik, Ruhr-Universität Bochum, Universitätsstrasse 150, 44801 Bochum (Germany); Fichtner, Horst [Theoretische Physik IV: Weltraum- und Astrophysik, Ruhr-Universität Bochum, Universitätsstrasse 150, 44801 Bochum (Germany); Sigl, Günter, E-mail: lukas.merten@rub.de, E-mail: julia.tjus@rub.de, E-mail: hf@tp4.rub.de, E-mail: eiche@tp4.rub.de, E-mail: guenter.sigl@desy.de [II Institut für Theoretische Physik, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg (Germany)

2017-06-01

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ{sub ∥} and perpendicular κ{sub ⊥} diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

CRPropa 3.1—a low energy extension based on stochastic differential equations

International Nuclear Information System (INIS)

Merten, Lukas; Tjus, Julia Becker; Eichmann, Björn; Fichtner, Horst; Sigl, Günter

2017-01-01

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ ∥ and perpendicular κ ⊥ diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

CRPropa 3.1—a low energy extension based on stochastic differential equations

Science.gov (United States)

Merten, Lukas; Becker Tjus, Julia; Fichtner, Horst; Eichmann, Björn; Sigl, Günter

2017-06-01

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ∥ and perpendicular κ⊥ diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

On stochastic differential equations with arbitrarily slow convergence rates for strong approximation in two space dimensions.

Science.gov (United States)

Gerencsér, Máté; Jentzen, Arnulf; Salimova, Diyora

2017-11-01

In a recent article (Jentzen et al. 2016 Commun. Math. Sci. 14 , 1477-1500 (doi:10.4310/CMS.2016.v14.n6.a1)), it has been established that, for every arbitrarily slow convergence speed and every natural number d ∈{4,5,…}, there exist d -dimensional stochastic differential equations with infinitely often differentiable and globally bounded coefficients such that no approximation method based on finitely many observations of the driving Brownian motion can converge in absolute mean to the solution faster than the given speed of convergence. In this paper, we strengthen the above result by proving that this slow convergence phenomenon also arises in two ( d =2) and three ( d =3) space dimensions.

Price dynamics of the financial markets using the stochastic differential equation for a potential double well

Science.gov (United States)

Lima, L. S.; Miranda, L. L. B.

2018-01-01

We have used the Itô's stochastic differential equation for the double well with additive white noise as a mathematical model for price dynamics of the financial market. We have presented a model which allows us to test within the same framework the comparative explanatory power of rational agents versus irrational agents, with respect to the facts of financial markets. We have obtained the mean price in terms of the β parameter that represents the force of the randomness term of the model.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics.

Science.gov (United States)

Cotter, C J; Gottwald, G A; Holm, D D

2017-09-01

In Holm (Holm 2015 Proc. R. Soc. A 471 , 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow.

Exponentially Stable Stationary Solutions for Stochastic Evolution Equations and Their Perturbation

International Nuclear Information System (INIS)

Caraballo, Tomas; Kloeden, Peter E.; Schmalfuss, Bjoern

2004-01-01

We consider the exponential stability of stochastic evolution equations with Lipschitz continuous non-linearities when zero is not a solution for these equations. We prove the existence of anon-trivial stationary solution which is exponentially stable, where the stationary solution is generated by the composition of a random variable and the Wiener shift. We also construct stationary solutions with the stronger property of attracting bounded sets uniformly. The existence of these stationary solutions follows from the theory of random dynamical systems and their attractors. In addition, we prove some perturbation results and formulate conditions for the existence of stationary solutions for semilinear stochastic partial differential equations with Lipschitz continuous non-linearities

The Pathwise Numerical Approximation of Stationary Solutions of Semilinear Stochastic Evolution Equations

International Nuclear Information System (INIS)

Caraballo, T.; Kloeden, P.E.

2006-01-01

Under a one-sided dissipative Lipschitz condition on its drift, a stochastic evolution equation with additive noise of the reaction-diffusion type is shown to have a unique stochastic stationary solution which pathwise attracts all other solutions. A similar situation holds for each Galerkin approximation and each implicit Euler scheme applied to these Galerkin approximations. Moreover, the stationary solution of the Euler scheme converges pathwise to that of the Galerkin system as the stepsize tends to zero and the stationary solutions of the Galerkin systems converge pathwise to that of the evolution equation as the dimension increases. The analysis is carried out on random partial and ordinary differential equations obtained from their stochastic counterparts by subtraction of appropriate Ornstein-Uhlenbeck stationary solutions

Indirect Inference for Stochastic Differential Equations Based on Moment Expansions

KAUST Repository

Ballesio, Marco

2016-01-06

We provide an indirect inference method to estimate the parameters of timehomogeneous scalar diffusion and jump diffusion processes. We obtain a system of ODEs that approximate the time evolution of the first two moments of the process by the approximation of the stochastic model applying a second order Taylor expansion of the SDE s infinitesimal generator in the Dynkin s formula. This method allows a simple and efficient procedure to infer the parameters of such stochastic processes given the data by the maximization of the likelihood of an approximating Gaussian process described by the two moments equations. Finally, we perform numerical experiments for two datasets arising from organic and inorganic fouling deposition phenomena.

Application of Monte Carlo method to solving boundary value problem of differential equations

International Nuclear Information System (INIS)

Zuo Yinghong; Wang Jianguo

2012-01-01

This paper introduces the foundation of the Monte Carlo method and the way how to generate the random numbers. Based on the basic thought of the Monte Carlo method and finite differential method, the stochastic model for solving the boundary value problem of differential equations is built. To investigate the application of the Monte Carlo method to solving the boundary value problem of differential equations, the model is used to solve Laplace's equations with the first boundary condition and the unsteady heat transfer equation with initial values and boundary conditions. The results show that the boundary value problem of differential equations can be effectively solved with the Monte Carlo method, and the differential equations with initial condition can also be calculated by using a stochastic probability model which is based on the time-domain finite differential equations. Both the simulation results and theoretical analyses show that the errors of numerical results are lowered as the number of simulation particles is increased. (authors)

Perturbation theory for continuous stochastic equations

International Nuclear Information System (INIS)

Chechetkin, V.R.; Lutovinov, V.S.

1987-01-01

The various general perturbational schemes for continuous stochastic equations are considered. These schemes have many analogous features with the iterational solution of Schwinger equation for S-matrix. The following problems are discussed: continuous stochastic evolution equations for probability distribution functionals, evolution equations for equal time correlators, perturbation theory for Gaussian and Poissonian additive noise, perturbation theory for birth and death processes, stochastic properties of systems with multiplicative noise. The general results are illustrated by diffusion-controlled reactions, fluctuations in closed systems with chemical processes, propagation of waves in random media in parabolic equation approximation, and non-equilibrium phase transitions in systems with Poissonian breeding centers. The rate of irreversible reaction X + X → A (Smoluchowski process) is calculated with the use of general theory based on continuous stochastic equations for birth and death processes. The threshold criterion and range of fluctuational region for synergetic phase transition in system with Poissonian breeding centers are also considered. (author)

Stochastic differential equations as a tool to regularize the parameter estimation problem for continuous time dynamical systems given discrete time measurements.

Science.gov (United States)

Leander, Jacob; Lundh, Torbjörn; Jirstrand, Mats

2014-05-01

In this paper we consider the problem of estimating parameters in ordinary differential equations given discrete time experimental data. The impact of going from an ordinary to a stochastic differential equation setting is investigated as a tool to overcome the problem of local minima in the objective function. Using two different models, it is demonstrated that by allowing noise in the underlying model itself, the objective functions to be minimized in the parameter estimation procedures are regularized in the sense that the number of local minima is reduced and better convergence is achieved. The advantage of using stochastic differential equations is that the actual states in the model are predicted from data and this will allow the prediction to stay close to data even when the parameters in the model is incorrect. The extended Kalman filter is used as a state estimator and sensitivity equations are provided to give an accurate calculation of the gradient of the objective function. The method is illustrated using in silico data from the FitzHugh-Nagumo model for excitable media and the Lotka-Volterra predator-prey system. The proposed method performs well on the models considered, and is able to regularize the objective function in both models. This leads to parameter estimation problems with fewer local minima which can be solved by efficient gradient-based methods. Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.

Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances

Directory of Open Access Journals (Sweden)

Wugan Cai

2017-02-01

Full Text Available Understanding the stochastic nature of emissions allowances is crucial for risk management in emissions trading markets. In this study, we discuss the emissions allowances spot price within the European Union Emissions Trading Scheme: Powernext and European Climate Exchange. To compare the fitness of five stochastic differential equations (SDEs to the European Union allowances spot price, we apply regression theory to obtain the point and interval estimations for the parameters of the SDEs. An empirical evaluation demonstrates that the mean reverting square root process (MRSRP has the best fitness of five SDEs to forecast the spot price. To reduce the degree of smog, we develop a new trading scheme in which firms have to hand many more allowances to the government when they emit one unit of air pollution on heavy pollution days, versus one allowance on clean days. Thus, we set up the SDE MRSRP model with Markovian switching to analyse the evolution of the spot price in such a scheme. The analysis shows that the allowances spot price will not jump too much in the new scheme. The findings of this study could contribute to developing a new type of emissions trading.

A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data

KAUST Repository

Babuška, Ivo; Nobile, Fabio; Tempone, Raul

2010-01-01

This work proposes and analyzes a stochastic collocation method for solving elliptic partial differential equations with random coefficients and forcing terms. These input data are assumed to depend on a finite number of random variables. The method consists of a Galerkin approximation in space and a collocation in the zeros of suitable tensor product orthogonal polynomials (Gauss points) in the probability space, and naturally leads to the solution of uncoupled deterministic problems as in the Monte Carlo approach. It treats easily a wide range of situations, such as input data that depend nonlinearly on the random variables, diffusivity coefficients with unbounded second moments, and random variables that are correlated or even unbounded. We provide a rigorous convergence analysis and demonstrate exponential convergence of the “probability error” with respect to the number of Gauss points in each direction of the probability space, under some regularity assumptions on the random input data. Numerical examples show the effectiveness of the method. Finally, we include a section with developments posterior to the original publication of this work. There we review sparse grid stochastic collocation methods, which are effective collocation strategies for problems that depend on a moderately large number of random variables.

Stochastic Analysis 2010

CERN Document Server

Crisan, Dan

2011-01-01

"Stochastic Analysis" aims to provide mathematical tools to describe and model high dimensional random systems. Such tools arise in the study of Stochastic Differential Equations and Stochastic Partial Differential Equations, Infinite Dimensional Stochastic Geometry, Random Media and Interacting Particle Systems, Super-processes, Stochastic Filtering, Mathematical Finance, etc. Stochastic Analysis has emerged as a core area of late 20th century Mathematics and is currently undergoing a rapid scientific development. The special volume "Stochastic Analysis 2010" provides a sa

Backward Stochastic Riccati Equations and Infinite Horizon L-Q Optimal Control with Infinite Dimensional State Space and Random Coefficients

International Nuclear Information System (INIS)

Guatteri, Giuseppina; Tessitore, Gianmario

2008-01-01

We study the Riccati equation arising in a class of quadratic optimal control problems with infinite dimensional stochastic differential state equation and infinite horizon cost functional. We allow the coefficients, both in the state equation and in the cost, to be random.In such a context backward stochastic Riccati equations are backward stochastic differential equations in the whole positive real axis that involve quadratic non-linearities and take values in a non-Hilbertian space. We prove existence of a minimal non-negative solution and, under additional assumptions, its uniqueness. We show that such a solution allows to perform the synthesis of the optimal control and investigate its attractivity properties. Finally the case where the coefficients are stationary is addressed and an example concerning a controlled wave equation in random media is proposed

Introduction to fractional and pseudo-differential equations with singular symbols

CERN Document Server

Umarov, Sabir

2015-01-01

The book systematically presents the theories of pseudo-differential operators with symbols singular in dual variables, fractional order derivatives, distributed and variable order fractional derivatives, random walk approximants, and applications of these theories to various initial and multi-point boundary value problems for pseudo-differential equations. Fractional Fokker-Planck-Kolmogorov equations associated with a large class of stochastic processes are presented. A complex version of the theory of pseudo-differential operators with meromorphic symbols based on the recently introduced complex Fourier transform is developed and applied for initial and boundary value problems for systems of complex differential and pseudo-differential equations.

Illness-death model: statistical perspective and differential equations.

Science.gov (United States)

Brinks, Ralph; Hoyer, Annika

2018-01-27

The aim of this work is to relate the theory of stochastic processes with the differential equations associated with multistate (compartment) models. We show that the Kolmogorov Forward Differential Equations can be used to derive a relation between the prevalence and the transition rates in the illness-death model. Then, we prove mathematical well-definedness and epidemiological meaningfulness of the prevalence of the disease. As an application, we derive the incidence of diabetes from a series of cross-sections.

A stochastic differential equation model for the foraging behavior of fish schools.

Science.gov (United States)

Tạ, Tôn Việt; Nguyen, Linh Thi Hoai

2018-03-15

Constructing models of living organisms locating food sources has important implications for understanding animal behavior and for the development of distribution technologies. This paper presents a novel simple model of stochastic differential equations for the foraging behavior of fish schools in a space including obstacles. The model is studied numerically. Three configurations of space with various food locations are considered. In the first configuration, fish swim in free but limited space. All individuals can find food with large probability while keeping their school structure. In the second and third configurations, they move in limited space with one and two obstacles, respectively. Our results reveal that the probability of foraging success is highest in the first configuration, and smallest in the third one. Furthermore, when school size increases up to an optimal value, the probability of foraging success tends to increase. When it exceeds an optimal value, the probability tends to decrease. The results agree with experimental observations.

A stochastic differential equation model for the foraging behavior of fish schools

Science.gov (United States)

Tạ, Tôn ệt, Vi; Hoai Nguyen, Linh Thi

2018-05-01

Constructing models of living organisms locating food sources has important implications for understanding animal behavior and for the development of distribution technologies. This paper presents a novel simple model of stochastic differential equations for the foraging behavior of fish schools in a space including obstacles. The model is studied numerically. Three configurations of space with various food locations are considered. In the first configuration, fish swim in free but limited space. All individuals can find food with large probability while keeping their school structure. In the second and third configurations, they move in limited space with one and two obstacles, respectively. Our results reveal that the probability of foraging success is highest in the first configuration, and smallest in the third one. Furthermore, when school size increases up to an optimal value, the probability of foraging success tends to increase. When it exceeds an optimal value, the probability tends to decrease. The results agree with experimental observations.

Stochastic porous media equations

CERN Document Server

Barbu, Viorel; Röckner, Michael

2016-01-01

Focusing on stochastic porous media equations, this book places an emphasis on existence theorems, asymptotic behavior and ergodic properties of the associated transition semigroup. Stochastic perturbations of the porous media equation have reviously been considered by physicists, but rigorous mathematical existence results have only recently been found. The porous media equation models a number of different physical phenomena, including the flow of an ideal gas and the diffusion of a compressible fluid through porous media, and also thermal propagation in plasma and plasma radiation. Another important application is to a model of the standard self-organized criticality process, called the "sand-pile model" or the "Bak-Tang-Wiesenfeld model". The book will be of interest to PhD students and researchers in mathematics, physics and biology.

Predictive performance for population models using stochastic differential equations applied on data from an oral glucose tolerance test

DEFF Research Database (Denmark)

Møller, Jonas Bech; Overgaard, R.V.; Madsen, Henrik

2010-01-01

Several articles have investigated stochastic differential equations (SDEs) in PK/PD models, but few have quantitatively investigated the benefits to predictive performance of models based on real data. Estimation of first phase insulin secretion which reflects beta-cell function using models of ...... obtained from the glucose tolerance tests. Since, the estimation time of extended models was not heavily increased compared to basic models, the applied method is concluded to have high relevance not only in theory but also in practice....

Exponential mean-square stability of two classes of theta Milstein methods for stochastic delay differential equations

Science.gov (United States)

Rouz, Omid Farkhondeh; Ahmadian, Davood; Milev, Mariyan

2017-12-01

This paper establishes exponential mean square stability of two classes of theta Milstein methods, namely split-step theta Milstein (SSTM) method and stochastic theta Milstein (STM) method, for stochastic differential delay equations (SDDEs). We consider the SDDEs problem under a coupled monotone condition on drift and diffusion coefficients, as well as a necessary linear growth condition on the last term of theta Milstein method. It is proved that the SSTM method with θ ∈ [0, ½] can recover the exponential mean square stability of the exact solution with some restrictive conditions on stepsize, but for θ ∈ (½, 1], we proved that the stability results hold for any stepsize. Then, based on the stability results of SSTM method, we examine the exponential mean square stability of the STM method and obtain the similar stability results to that of the SSTM method. In the numerical section the figures show thevalidity of our claims.

Applied analysis and differential equations

CERN Document Server

Cârj, Ovidiu

2007-01-01

This volume contains refereed research articles written by experts in the field of applied analysis, differential equations and related topics. Well-known leading mathematicians worldwide and prominent young scientists cover a diverse range of topics, including the most exciting recent developments. A broad range of topics of recent interest are treated: existence, uniqueness, viability, asymptotic stability, viscosity solutions, controllability and numerical analysis for ODE, PDE and stochastic equations. The scope of the book is wide, ranging from pure mathematics to various applied fields such as classical mechanics, biomedicine, and population dynamics.

On the pth moment estimates of solutions to stochastic functional differential equations in the G-framework.

Science.gov (United States)

Faizullah, Faiz

2016-01-01

The aim of the current paper is to present the path-wise and moment estimates for solutions to stochastic functional differential equations with non-linear growth condition in the framework of G-expectation and G-Brownian motion. Under the nonlinear growth condition, the pth moment estimates for solutions to SFDEs driven by G-Brownian motion are proved. The properties of G-expectations, Hölder's inequality, Bihari's inequality, Gronwall's inequality and Burkholder-Davis-Gundy inequalities are used to develop the above mentioned theory. In addition, the path-wise asymptotic estimates and continuity of pth moment for the solutions to SFDEs in the G-framework, with non-linear growth condition are shown.

Functional stochastic differential equations: mathematical theory of nonlinear parabolic systems with applications in field theory and statistical mechanics

International Nuclear Information System (INIS)

Doering, C.R.

1985-01-01

Applications of nonlinear parabolic stochastic differential equations with additive colored noise in equilibrium and nonequilibrium statistical mechanics and quantum field theory are developed in detail, providing a new unified mathematical approach to many problems. The existence and uniqueness of solutions to these equations is established, and some of the properties of the solutions are investigated. In particular, asymptotic expansions for the correlation functions of the solutions are introduced and compared to rigorous nonperturbative bounds on the moments. It is found that the perturbative analysis is in qualitative disagreement with the exact result in models corresponding to cut-off self-interacting nonperturbatively renormalizable scalar quantum field theories. For these theories the nonlinearities cannot be considered as perturbations of the linearized theory

System Entropy Measurement of Stochastic Partial Differential Systems

Directory of Open Access Journals (Sweden)

Bor-Sen Chen

2016-03-01

Full Text Available System entropy describes the dispersal of a system’s energy and is an indication of the disorder of a physical system. Several system entropy measurement methods have been developed for dynamic systems. However, most real physical systems are always modeled using stochastic partial differential dynamic equations in the spatio-temporal domain. No efficient method currently exists that can calculate the system entropy of stochastic partial differential systems (SPDSs in consideration of the effects of intrinsic random fluctuation and compartment diffusion. In this study, a novel indirect measurement method is proposed for calculating of system entropy of SPDSs using a Hamilton–Jacobi integral inequality (HJII-constrained optimization method. In other words, we solve a nonlinear HJII-constrained optimization problem for measuring the system entropy of nonlinear stochastic partial differential systems (NSPDSs. To simplify the system entropy measurement of NSPDSs, the global linearization technique and finite difference scheme were employed to approximate the nonlinear stochastic spatial state space system. This allows the nonlinear HJII-constrained optimization problem for the system entropy measurement to be transformed to an equivalent linear matrix inequalities (LMIs-constrained optimization problem, which can be easily solved using the MATLAB LMI-toolbox (MATLAB R2014a, version 8.3. Finally, several examples are presented to illustrate the system entropy measurement of SPDSs.

Boosting Bayesian parameter inference of nonlinear stochastic differential equation models by Hamiltonian scale separation.

Science.gov (United States)

Albert, Carlo; Ulzega, Simone; Stoop, Ruedi

2016-04-01

Parameter inference is a fundamental problem in data-driven modeling. Given observed data that is believed to be a realization of some parameterized model, the aim is to find parameter values that are able to explain the observed data. In many situations, the dominant sources of uncertainty must be included into the model for making reliable predictions. This naturally leads to stochastic models. Stochastic models render parameter inference much harder, as the aim then is to find a distribution of likely parameter values. In Bayesian statistics, which is a consistent framework for data-driven learning, this so-called posterior distribution can be used to make probabilistic predictions. We propose a novel, exact, and very efficient approach for generating posterior parameter distributions for stochastic differential equation models calibrated to measured time series. The algorithm is inspired by reinterpreting the posterior distribution as a statistical mechanics partition function of an object akin to a polymer, where the measurements are mapped on heavier beads compared to those of the simulated data. To arrive at distribution samples, we employ a Hamiltonian Monte Carlo approach combined with a multiple time-scale integration. A separation of time scales naturally arises if either the number of measurement points or the number of simulation points becomes large. Furthermore, at least for one-dimensional problems, we can decouple the harmonic modes between measurement points and solve the fastest part of their dynamics analytically. Our approach is applicable to a wide range of inference problems and is highly parallelizable.

A stochastic differential equations approach for the description of helium bubble size distributions in irradiated metals

Energy Technology Data Exchange (ETDEWEB)

Seif, Dariush, E-mail: dariush.seif@iwm-extern.fraunhofer.de [Fraunhofer Institut für Werkstoffmechanik, Freiburg 79108 (Germany); Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095-1597 (United States); Ghoniem, Nasr M. [Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095-1597 (United States)

2014-12-15

A rate theory model based on the theory of nonlinear stochastic differential equations (SDEs) is developed to estimate the time-dependent size distribution of helium bubbles in metals under irradiation. Using approaches derived from Itô’s calculus, rate equations for the first five moments of the size distribution in helium–vacancy space are derived, accounting for the stochastic nature of the atomic processes involved. In the first iteration of the model, the distribution is represented as a bivariate Gaussian distribution. The spread of the distribution about the mean is obtained by white-noise terms in the second-order moments, driven by fluctuations in the general absorption and emission of point defects by bubbles, and fluctuations stemming from collision cascades. This statistical model for the reconstruction of the distribution by its moments is coupled to a previously developed reduced-set, mean-field, rate theory model. As an illustrative case study, the model is applied to a tungsten plasma facing component under irradiation. Our findings highlight the important role of stochastic atomic fluctuations on the evolution of helium–vacancy cluster size distributions. It is found that when the average bubble size is small (at low dpa levels), the relative spread of the distribution is large and average bubble pressures may be very large. As bubbles begin to grow in size, average bubble pressures decrease, and stochastic fluctuations have a lessened effect. The distribution becomes tighter as it evolves in time, corresponding to a more uniform bubble population. The model is formulated in a general way, capable of including point defect drift due to internal temperature and/or stress gradients. These arise during pulsed irradiation, and also during steady irradiation as a result of externally applied or internally generated non-homogeneous stress fields. Discussion is given into how the model can be extended to include full spatial resolution and how the

A stochastic differential equations approach for the description of helium bubble size distributions in irradiated metals

Science.gov (United States)

Seif, Dariush; Ghoniem, Nasr M.

2014-12-01

A rate theory model based on the theory of nonlinear stochastic differential equations (SDEs) is developed to estimate the time-dependent size distribution of helium bubbles in metals under irradiation. Using approaches derived from Itô's calculus, rate equations for the first five moments of the size distribution in helium-vacancy space are derived, accounting for the stochastic nature of the atomic processes involved. In the first iteration of the model, the distribution is represented as a bivariate Gaussian distribution. The spread of the distribution about the mean is obtained by white-noise terms in the second-order moments, driven by fluctuations in the general absorption and emission of point defects by bubbles, and fluctuations stemming from collision cascades. This statistical model for the reconstruction of the distribution by its moments is coupled to a previously developed reduced-set, mean-field, rate theory model. As an illustrative case study, the model is applied to a tungsten plasma facing component under irradiation. Our findings highlight the important role of stochastic atomic fluctuations on the evolution of helium-vacancy cluster size distributions. It is found that when the average bubble size is small (at low dpa levels), the relative spread of the distribution is large and average bubble pressures may be very large. As bubbles begin to grow in size, average bubble pressures decrease, and stochastic fluctuations have a lessened effect. The distribution becomes tighter as it evolves in time, corresponding to a more uniform bubble population. The model is formulated in a general way, capable of including point defect drift due to internal temperature and/or stress gradients. These arise during pulsed irradiation, and also during steady irradiation as a result of externally applied or internally generated non-homogeneous stress fields. Discussion is given into how the model can be extended to include full spatial resolution and how the

A stochastic differential equations approach for the description of helium bubble size distributions in irradiated metals

International Nuclear Information System (INIS)

Seif, Dariush; Ghoniem, Nasr M.

2014-01-01

A rate theory model based on the theory of nonlinear stochastic differential equations (SDEs) is developed to estimate the time-dependent size distribution of helium bubbles in metals under irradiation. Using approaches derived from Itô’s calculus, rate equations for the first five moments of the size distribution in helium–vacancy space are derived, accounting for the stochastic nature of the atomic processes involved. In the first iteration of the model, the distribution is represented as a bivariate Gaussian distribution. The spread of the distribution about the mean is obtained by white-noise terms in the second-order moments, driven by fluctuations in the general absorption and emission of point defects by bubbles, and fluctuations stemming from collision cascades. This statistical model for the reconstruction of the distribution by its moments is coupled to a previously developed reduced-set, mean-field, rate theory model. As an illustrative case study, the model is applied to a tungsten plasma facing component under irradiation. Our findings highlight the important role of stochastic atomic fluctuations on the evolution of helium–vacancy cluster size distributions. It is found that when the average bubble size is small (at low dpa levels), the relative spread of the distribution is large and average bubble pressures may be very large. As bubbles begin to grow in size, average bubble pressures decrease, and stochastic fluctuations have a lessened effect. The distribution becomes tighter as it evolves in time, corresponding to a more uniform bubble population. The model is formulated in a general way, capable of including point defect drift due to internal temperature and/or stress gradients. These arise during pulsed irradiation, and also during steady irradiation as a result of externally applied or internally generated non-homogeneous stress fields. Discussion is given into how the model can be extended to include full spatial resolution and how the

Differential Equations Models to Study Quorum Sensing.

Science.gov (United States)

Pérez-Velázquez, Judith; Hense, Burkhard A

2018-01-01

Mathematical models to study quorum sensing (QS) have become an important tool to explore all aspects of this type of bacterial communication. A wide spectrum of mathematical tools and methods such as dynamical systems, stochastics, and spatial models can be employed. In this chapter, we focus on giving an overview of models consisting of differential equations (DE), which can be used to describe changing quantities, for example, the dynamics of one or more signaling molecule in time and space, often in conjunction with bacterial growth dynamics. The chapter is divided into two sections: ordinary differential equations (ODE) and partial differential equations (PDE) models of QS. Rates of change are represented mathematically by derivatives, i.e., in terms of DE. ODE models allow describing changes in one independent variable, for example, time. PDE models can be used to follow changes in more than one independent variable, for example, time and space. Both types of models often consist of systems (i.e., more than one equation) of equations, such as equations for bacterial growth and autoinducer concentration dynamics. Almost from the onset, mathematical modeling of QS using differential equations has been an interdisciplinary endeavor and many of the works we revised here will be placed into their biological context.

High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr; Cohen, David; Vilmart, Gilles; Zygalakis, Konstantinos C.

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration

xSPDE: Extensible software for stochastic equations

Directory of Open Access Journals (Sweden)

Simon Kiesewetter

2016-01-01

Full Text Available We introduce an extensible software toolbox, xSPDE, for solving ordinary and partial stochastic differential equations. The toolbox makes extensive use of vector and parallel methods. Inputs are exceptionally simple, to reduce the learning curve, with default options for all of the many input parameters. The code calculates functional means, correlations and spectra, checks for errors in both time-step and sampling, and provides several choices of algorithm. Most aspects of the code, including the numerical algorithm, have a modular functional design to allow user modifications.

Computable Error Estimates for Finite Element Approximations of Elliptic Partial Differential Equations with Rough Stochastic Data

KAUST Repository

Hall, Eric Joseph

2016-12-08

We derive computable error estimates for finite element approximations of linear elliptic partial differential equations with rough stochastic coefficients. In this setting, the exact solutions contain high frequency content that standard a posteriori error estimates fail to capture. We propose goal-oriented estimates, based on local error indicators, for the pathwise Galerkin and expected quadrature errors committed in standard, continuous, piecewise linear finite element approximations. Derived using easily validated assumptions, these novel estimates can be computed at a relatively low cost and have applications to subsurface flow problems in geophysics where the conductivities are assumed to have lognormal distributions with low regularity. Our theory is supported by numerical experiments on test problems in one and two dimensions.

Optimal control of stochastic difference Volterra equations an introduction

CERN Document Server

Shaikhet, Leonid

2015-01-01

This book showcases a subclass of hereditary systems, that is, systems with behaviour depending not only on their current state but also on their past history; it is an introduction to the mathematical theory of optimal control for stochastic difference Volterra equations of neutral type. As such, it will be of much interest to researchers interested in modelling processes in physics, mechanics, automatic regulation, economics and finance, biology, sociology and medicine for all of which such equations are very popular tools. The text deals with problems of optimal control such as meeting given performance criteria, and stabilization, extending them to neutral stochastic difference Volterra equations. In particular, it contrasts the difference analogues of solutions to optimal control and optimal estimation problems for stochastic integral Volterra equations with optimal solutions for corresponding problems in stochastic difference Volterra equations. Optimal Control of Stochastic Difference Volterra Equation...

On the existence of weak solutions of quantum stochastic differential ...

African Journals Online (AJOL)

We establish further results concerning the existence, uniqueness and stability of weak solutions of quantum stochastic differential equations (QSDEs). Our results are achieved by considering a more general Lipschit condition on the coefficients than our previous considerations in [1]. We exhibit a class of Lipschitzian ...

Advances in nonlinear partial differential equations and stochastics

CERN Document Server

Kawashima, S

1998-01-01

In the past two decades, there has been great progress in the theory of nonlinear partial differential equations. This book describes the progress, focusing on interesting topics in gas dynamics, fluid dynamics, elastodynamics etc. It contains ten articles, each of which discusses a very recent result obtained by the author. Some of these articles review related results.

Control of Stochastic Master Equation Models of Genetic Regulatory Networks by Approximating Their Average Behavior

Science.gov (United States)

Umut Caglar, Mehmet; Pal, Ranadip

2010-10-01

The central dogma of molecular biology states that ``information cannot be transferred back from protein to either protein or nucleic acid.'' However, this assumption is not exactly correct in most of the cases. There are a lot of feedback loops and interactions between different levels of systems. These types of interactions are hard to analyze due to the lack of data in the cellular level and probabilistic nature of interactions. Probabilistic models like Stochastic Master Equation (SME) or deterministic models like differential equations (DE) can be used to analyze these types of interactions. SME models based on chemical master equation (CME) can provide detailed representation of genetic regulatory system, but their use is restricted by the large data requirements and computational costs of calculations. The differential equations models on the other hand, have low calculation costs and much more adequate to generate control procedures on the system; but they are not adequate to investigate the probabilistic nature of interactions. In this work the success of the mapping between SME and DE is analyzed, and the success of a control policy generated by DE model with respect to SME model is examined. Index Terms--- Stochastic Master Equation models, Differential Equation Models, Control Policy Design, Systems biology

PKPD model of interleukin-21 effects on thermoregulation in monkeys--application and evaluation of stochastic differential equations.

Science.gov (United States)

Overgaard, Rune Viig; Holford, Nick; Rytved, Klaus A; Madsen, Henrik

2007-02-01

To describe the pharmacodynamic effects of recombinant human interleukin-21 (IL-21) on core body temperature in cynomolgus monkeys using basic mechanisms of heat regulation. A major effort was devoted to compare the use of ordinary differential equations (ODEs) with stochastic differential equations (SDEs) in pharmacokinetic pharmacodynamic (PKPD) modelling. A temperature model was formulated including circadian rhythm, metabolism, heat loss, and a thermoregulatory set-point. This model was formulated as a mixed-effects model based on SDEs using NONMEM. The effects of IL-21 were on the set-point and the circadian rhythm of metabolism. The model was able to describe a complex set of IL-21 induced phenomena, including 1) disappearance of the circadian rhythm, 2) no effect after first dose, and 3) high variability after second dose. SDEs provided a more realistic description with improved simulation properties, and further changed the model into one that could not be falsified by the autocorrelation function. The IL-21 induced effects on thermoregulation in cynomolgus monkeys are explained by a biologically plausible model. The quality of the model was improved by the use of SDEs.

Robustness of Quadratic Hedging Strategies in Finance via Backward Stochastic Differential Equations with Jumps

International Nuclear Information System (INIS)

Di Nunno, Giulia; Khedher, Asma; Vanmaele, Michèle

2015-01-01

We consider a backward stochastic differential equation with jumps (BSDEJ) which is driven by a Brownian motion and a Poisson random measure. We present two candidate-approximations to this BSDEJ and we prove that the solution of each candidate-approximation converges to the solution of the original BSDEJ in a space which we specify. We use this result to investigate in further detail the consequences of the choice of the model to (partial) hedging in incomplete markets in finance. As an application, we consider models in which the small variations in the price dynamics are modeled with a Poisson random measure with infinite activity and models in which these small variations are modeled with a Brownian motion or are cut off. Using the convergence results on BSDEJs, we show that quadratic hedging strategies are robust towards the approximation of the market prices and we derive an estimation of the model risk

Robustness of Quadratic Hedging Strategies in Finance via Backward Stochastic Differential Equations with Jumps

Energy Technology Data Exchange (ETDEWEB)

Di Nunno, Giulia, E-mail: giulian@math.uio.no [University of Oslo, Center of Mathematics for Applications (Norway); Khedher, Asma, E-mail: asma.khedher@tum.de [Technische Universität München, Chair of Mathematical Finance (Germany); Vanmaele, Michèle, E-mail: michele.vanmaele@ugent.be [Ghent University, Department of Applied Mathematics, Computer Science and Statistics (Belgium)

2015-12-15

We consider a backward stochastic differential equation with jumps (BSDEJ) which is driven by a Brownian motion and a Poisson random measure. We present two candidate-approximations to this BSDEJ and we prove that the solution of each candidate-approximation converges to the solution of the original BSDEJ in a space which we specify. We use this result to investigate in further detail the consequences of the choice of the model to (partial) hedging in incomplete markets in finance. As an application, we consider models in which the small variations in the price dynamics are modeled with a Poisson random measure with infinite activity and models in which these small variations are modeled with a Brownian motion or are cut off. Using the convergence results on BSDEJs, we show that quadratic hedging strategies are robust towards the approximation of the market prices and we derive an estimation of the model risk.

Stochastic Computational Approach for Complex Nonlinear Ordinary Differential Equations

International Nuclear Information System (INIS)

Khan, Junaid Ali; Raja, Muhammad Asif Zahoor; Qureshi, Ijaz Mansoor

2011-01-01

We present an evolutionary computational approach for the solution of nonlinear ordinary differential equations (NLODEs). The mathematical modeling is performed by a feed-forward artificial neural network that defines an unsupervised error. The training of these networks is achieved by a hybrid intelligent algorithm, a combination of global search with genetic algorithm and local search by pattern search technique. The applicability of this approach ranges from single order NLODEs, to systems of coupled differential equations. We illustrate the method by solving a variety of model problems and present comparisons with solutions obtained by exact methods and classical numerical methods. The solution is provided on a continuous finite time interval unlike the other numerical techniques with comparable accuracy. With the advent of neuroprocessors and digital signal processors the method becomes particularly interesting due to the expected essential gains in the execution speed. (general)

Efficient Multilevel and Multi-index Sampling Methods in Stochastic Differential Equations

KAUST Repository

Haji-Ali, Abdul Lateef

2016-05-22

of this thesis is the novel Multi-index Monte Carlo (MIMC) method which is an extension of MLMC in high dimensional problems with significant computational savings. Under reasonable assumptions on the weak and variance convergence, which are related to the mixed regularity of the underlying problem and the discretization method, the order of the computational complexity of MIMC is, at worst up to a logarithmic factor, independent of the dimensionality of the underlying parametric equation. We also apply the same multi-index methodology to another sampling method, namely the Stochastic Collocation method. Hence, the novel Multi-index Stochastic Collocation method is proposed and is shown to be more efficient in problems with sufficient mixed regularity than our novel MIMC method and other standard methods. Finally, MIMC is applied to approximate quantities of interest of stochastic particle systems in the mean-field when the number of particles tends to infinity. To approximate these quantities of interest up to an error tolerance, TOL, MIMC has a computational complexity of O(TOL-2log(TOL)2). This complexity is achieved by building a hierarchy based on two discretization parameters: the number of time steps in an Milstein scheme and the number of particles in the particle system. Moreover, we use a partitioning estimator to increase the correlation between two stochastic particle systems with different sizes. In comparison, the optimal computational complexity of MLMC in this case is O(TOL-3) and the computational complexity of Monte Carlo is O(TOL-4).

On Lipschitzian quantum stochastic differential inclusions

International Nuclear Information System (INIS)

Ekhaguere, G.O.S.

1990-12-01

Quantum stochastic differential inclusions are introduced and studied within the framework of the Hudson-Parthasarathy formulation of quantum stochastic calculus. Results concerning the existence of solutions of a Lipschitzian quantum stochastic differential inclusion and the relationship between the solutions of such an inclusion and those of its convexification are presented. These generalize the Filippov existence theorem and the Filippov-Wazewski Relaxation Theorem for classical differential inclusions to the present noncommutative setting. (author). 9 refs

Delayed Stochastic Linear-Quadratic Control Problem and Related Applications

Directory of Open Access Journals (Sweden)

Li Chen

2012-01-01

stochastic differential equations (FBSDEs with Itô’s stochastic delay equations as forward equations and anticipated backward stochastic differential equations as backward equations. Especially, we present the optimal feedback regulator for the time delay system via a new type of Riccati equations and also apply to a population optimal control problem.

Inviscid limit of stochastic damped 2D Navier–Stokes equations

International Nuclear Information System (INIS)

Bessaih, Hakima; Ferrario, Benedetta

2014-01-01

We consider the inviscid limit of the stochastic damped 2D Navier–Stokes equations. We prove that, when the viscosity vanishes, the stationary solution of the stochastic damped Navier–Stokes equations converges to a stationary solution of the stochastic damped Euler equation and that the rate of dissipation of enstrophy converges to zero. In particular, this limit obeys an enstrophy balance. The rates are computed with respect to a limit measure of the unique invariant measure of the stochastic damped Navier–Stokes equations. (paper)

Approximate reduction of linear population models governed by stochastic differential equations: application to multiregional models.

Science.gov (United States)

Sanz, Luis; Alonso, Juan Antonio

2017-12-01

In this work we develop approximate aggregation techniques in the context of slow-fast linear population models governed by stochastic differential equations and apply the results to the treatment of populations with spatial heterogeneity. Approximate aggregation techniques allow one to transform a complex system involving many coupled variables and in which there are processes with different time scales, by a simpler reduced model with a fewer number of 'global' variables, in such a way that the dynamics of the former can be approximated by that of the latter. In our model we contemplate a linear fast deterministic process together with a linear slow process in which the parameters are affected by additive noise, and give conditions for the solutions corresponding to positive initial conditions to remain positive for all times. By letting the fast process reach equilibrium we build a reduced system with a lesser number of variables, and provide results relating the asymptotic behaviour of the first- and second-order moments of the population vector for the original and the reduced system. The general technique is illustrated by analysing a multiregional stochastic system in which dispersal is deterministic and the rate growth of the populations in each patch is affected by additive noise.

RBSDE's with jumps and the related obstacle problems for integral-partial differential equations

Institute of Scientific and Technical Information of China (English)

FAN; Yulian

2006-01-01

The author proves, when the noise is driven by a Brownian motion and an independent Poisson random measure, the one-dimensional reflected backward stochastic differential equation with a stopping time terminal has a unique solution. And in a Markovian framework, the solution can provide a probabilistic interpretation for the obstacle problem for the integral-partial differential equation.

Stochastic differential equation (SDE) model of opening gold share price of bursa saham malaysia

Science.gov (United States)

Hussin, F. N.; Rahman, H. A.; Bahar, A.

2017-09-01

Black and Scholes option pricing model is one of the most recognized stochastic differential equation model in mathematical finance. Two parameter estimation methods have been utilized for the Geometric Brownian model (GBM); historical and discrete method. The historical method is a statistical method which uses the property of independence and normality logarithmic return, giving out the simplest parameter estimation. Meanwhile, discrete method considers the function of density of transition from the process of diffusion normal log which has been derived from maximum likelihood method. These two methods are used to find the parameter estimates samples of Malaysians Gold Share Price data such as: Financial Times and Stock Exchange (FTSE) Bursa Malaysia Emas, and Financial Times and Stock Exchange (FTSE) Bursa Malaysia Emas Shariah. Modelling of gold share price is essential since fluctuation of gold affects worldwide economy nowadays, including Malaysia. It is found that discrete method gives the best parameter estimates than historical method due to the smallest Root Mean Square Error (RMSE) value.

Density-based Monte Carlo filter and its applications in nonlinear stochastic differential equation models.

Science.gov (United States)

Huang, Guanghui; Wan, Jianping; Chen, Hui

2013-02-01

Nonlinear stochastic differential equation models with unobservable state variables are now widely used in analysis of PK/PD data. Unobservable state variables are usually estimated with extended Kalman filter (EKF), and the unknown pharmacokinetic parameters are usually estimated by maximum likelihood estimator. However, EKF is inadequate for nonlinear PK/PD models, and MLE is known to be biased downwards. A density-based Monte Carlo filter (DMF) is proposed to estimate the unobservable state variables, and a simulation-based M estimator is proposed to estimate the unknown parameters in this paper, where a genetic algorithm is designed to search the optimal values of pharmacokinetic parameters. The performances of EKF and DMF are compared through simulations for discrete time and continuous time systems respectively, and it is found that the results based on DMF are more accurate than those given by EKF with respect to mean absolute error. Copyright © 2012 Elsevier Ltd. All rights reserved.

Homogenization of the evolution Stokes equation in a perforated domain with a stochastic Fourier boundary condition

KAUST Repository

Bessaih, Hakima

2015-04-01

The evolution Stokes equation in a domain containing periodically distributed obstacles subject to Fourier boundary condition on the boundaries is considered. We assume that the dynamic is driven by a stochastic perturbation on the interior of the domain and another stochastic perturbation on the boundaries of the obstacles. We represent the solid obstacles by holes in the fluid domain. The macroscopic (homogenized) equation is derived as another stochastic partial differential equation, defined in the whole non perforated domain. Here, the initial stochastic perturbation on the boundary becomes part of the homogenized equation as another stochastic force. We use the twoscale convergence method after extending the solution with 0 in the holes to pass to the limit. By Itô stochastic calculus, we get uniform estimates on the solution in appropriate spaces. In order to pass to the limit on the boundary integrals, we rewrite them in terms of integrals in the whole domain. In particular, for the stochastic integral on the boundary, we combine the previous idea of rewriting it on the whole domain with the assumption that the Brownian motion is of trace class. Due to the particular boundary condition dealt with, we get that the solution of the stochastic homogenized equation is not divergence free. However, it is coupled with the cell problem that has a divergence free solution. This paper represents an extension of the results of Duan and Wang (Comm. Math. Phys. 275:1508-1527, 2007), where a reaction diffusion equation with a dynamical boundary condition with a noise source term on both the interior of the domain and on the boundary was studied, and through a tightness argument and a pointwise two scale convergence method the homogenized equation was derived. © American Institute of Mathematical Sciences.

Stochastic Levy Divergence and Maxwell's Equations

Directory of Open Access Journals (Sweden)

B. O. Volkov

2015-01-01

Full Text Available One of the main reasons for interest in the Levy Laplacian and its analogues such as Levy d'Alembertian is a connection of these operators with gauge fields. The theorem proved by Accardi, Gibillisco and Volovich stated that a connection in a bundle over a Euclidean space or over a Minkowski space is a solution of the Yang-Mills equations if and only if the corresponding parallel transport to the connection is a solution of the Laplace equation for the Levy Laplacian or of the d'Alembert equation for the Levy d'Alembertian respectively (see [5, 6]. There are two approaches to define Levy type operators, both of which date back to the original works of Levy [7]. The first is that the Levy Laplacian (or Levy d'Alembertian is defined as an integral functional generated by a special form of the second derivative. This approach is used in the works [5, 6], as well as in the paper [8] of Leandre and Volovich, where stochastic Levy-Laplacian is discussed. Another approach to the Levy Laplacian is defining it as the Cesaro mean of second order derivatives along the family of vectors, which is an orthonormal basis in the Hilbert space. This definition of the Levy Laplacian is used for the description of solutions of the Yang-Mills equations in the paper [10].The present work shows that the definitions of the Levy Laplacian and the Levy d'Alembertian based on Cesaro averaging of the second order directional derivatives can be transferred to the stochastic case. In the article the values of these operators on a stochastic parallel transport associated with a connection (vector potential are found. In this case, unlike the deterministic case and the stochastic case of Levy Laplacian from [8], these values are not equal to zero if the vector potential corresponding to the stochastic parallel transport is a solution of the Maxwell's equations. As a result, two approaches to definition of the Levy Laplacian in the stochastic case give different operators. This

Parameter estimation and change-point detection from Dynamic Contrast Enhanced MRI data using stochastic differential equations.

Science.gov (United States)

Cuenod, Charles-André; Favetto, Benjamin; Genon-Catalot, Valentine; Rozenholc, Yves; Samson, Adeline

2011-09-01

Dynamic Contrast Enhanced imaging (DCE-imaging) following a contrast agent bolus allows the extraction of information on tissue micro-vascularization. The dynamic signals obtained from DCE-imaging are modeled by pharmacokinetic compartmental models which integrate the Arterial Input Function. These models use ordinary differential equations (ODEs) to describe the exchanges between the arterial and capillary plasma and the extravascular-extracellular space. Their least squares fitting takes into account measurement noises but fails to deal with unpredictable fluctuations due to external/internal sources of variations (patients' anxiety, time-varying parameters, measurement errors in the input function, etc.). Adding Brownian components to the ODEs leads to stochastic differential equations (SDEs). In DCE-imaging, SDEs are discretely observed with an additional measurement noise. We propose to estimate the parameters of these noisy SDEs by maximum likelihood, using the Kalman filter. In DCE-imaging, the contrast agent injected in vein arrives in plasma with an unknown time delay. The delay parameter induces a change-point in the drift of the SDE and ODE models, which is estimated also. Estimations based on the SDE and ODE pharmacokinetic models are compared to real DCE-MRI data. They show that the use of SDE provides robustness in the estimation results. A simulation study confirms these results. Copyright © 2011 Elsevier Inc. All rights reserved.

Stochastic Perron's method and elementary strategies for zero-sum differential games

OpenAIRE

Sîrbu, Mihai

2013-01-01

We develop here the Stochastic Perron Method in the framework of two-player zero-sum differential games. We consider the formulation of the game where both players play, symmetrically, feed-back strategies (as in [CR09] or [PZ12]) as opposed to the Elliott-Kalton formulation prevalent in the literature. The class of feed-back strategies we use is carefully chosen so that the state equation admits strong solutions and the technicalities involved in the Stochastic Perron Method carry through in...

International Conference on Delay Differential and Difference Equations and Applications

CERN Document Server

Pituk, Mihály; Recent Advances in Delay Differential and Difference Equations

2014-01-01

Delay differential and difference equations serve as models for a range of processes in biology, physics, engineering, and control theory. In this volume, the participants of the International Conference on Delay Differential and Difference Equations and Applications, Balatonfüred, Hungary, July 15-19, 2013 present recent research in this quickly-evolving field. The papers relate to the existence, asymptotic, and oscillatory properties of the solutions; stability theory; numerical approximations; and applications to real world phenomena using deterministic and stochastic discrete and continuous dynamical systems.

Stochastic 2-D Navier-Stokes Equation

International Nuclear Information System (INIS)

Menaldi, J.L.; Sritharan, S.S.

2002-01-01

In this paper we prove the existence and uniqueness of strong solutions for the stochastic Navier-Stokes equation in bounded and unbounded domains. These solutions are stochastic analogs of the classical Lions-Prodi solutions to the deterministic Navier-Stokes equation. Local monotonicity of the nonlinearity is exploited to obtain the solutions in a given probability space and this significantly improves the earlier techniques for obtaining strong solutions, which depended on pathwise solutions to the Navier-Stokes martingale problem where the probability space is also obtained as a part of the solution

Quantum theory of open systems based on stochastic differential equations of generalized Langevin (non-Wiener) type

International Nuclear Information System (INIS)

Basharov, A. M.

2012-01-01

It is shown that the effective Hamiltonian representation, as it is formulated in author’s papers, serves as a basis for distinguishing, in a broadband environment of an open quantum system, independent noise sources that determine, in terms of the stationary quantum Wiener and Poisson processes in the Markov approximation, the effective Hamiltonian and the equation for the evolution operator of the open system and its environment. General stochastic differential equations of generalized Langevin (non-Wiener) type for the evolution operator and the kinetic equation for the density matrix of an open system are obtained, which allow one to analyze the dynamics of a wide class of localized open systems in the Markov approximation. The main distinctive features of the dynamics of open quantum systems described in this way are the stabilization of excited states with respect to collective processes and an additional frequency shift of the spectrum of the open system. As an illustration of the general approach developed, the photon dynamics in a single-mode cavity without losses on the mirrors is considered, which contains identical intracavity atoms coupled to the external vacuum electromagnetic field. For some atomic densities, the photons of the cavity mode are “locked” inside the cavity, thus exhibiting a new phenomenon of radiation trapping and non-Wiener dynamics.

Stochastic Analysis : A Series of Lectures

CERN Document Server

Dozzi, Marco; Flandoli, Franco; Russo, Francesco

2015-01-01

This book presents in thirteen refereed survey articles an overview of modern activity in stochastic analysis, written by leading international experts. The topics addressed include stochastic fluid dynamics and regularization by noise of deterministic dynamical systems; stochastic partial differential equations driven by Gaussian or Lévy noise, including the relationship between parabolic equations and particle systems, and wave equations in a geometric framework; Malliavin calculus and applications to stochastic numerics; stochastic integration in Banach spaces; porous media-type equations; stochastic deformations of classical mechanics and Feynman integrals and stochastic differential equations with reflection. The articles are based on short courses given at the Centre Interfacultaire Bernoulli of the Ecole Polytechnique Fédérale de Lausanne, Switzerland, from January to June 2012. They offer a valuable resource not only for specialists, but also for other researchers and Ph.D. students in the fields o...

Nonzero-Sum Stochastic Differential Portfolio Games under a Markovian Regime Switching Model

Directory of Open Access Journals (Sweden)

Chaoqun Ma

2015-01-01

Full Text Available We consider a nonzero-sum stochastic differential portfolio game problem in a continuous-time Markov regime switching environment when the price dynamics of the risky assets are governed by a Markov-modulated geometric Brownian motion (GBM. The market parameters, including the bank interest rate and the appreciation and volatility rates of the risky assets, switch over time according to a continuous-time Markov chain. We formulate the nonzero-sum stochastic differential portfolio game problem as two utility maximization problems of the sum process between two investors’ terminal wealth. We derive a pair of regime switching Hamilton-Jacobi-Bellman (HJB equations and two systems of coupled HJB equations at different regimes. We obtain explicit optimal portfolio strategies and Feynman-Kac representations of the two value functions. Furthermore, we solve the system of coupled HJB equations explicitly in a special case where there are only two states in the Markov chain. Finally we provide comparative statics and numerical simulation analysis of optimal portfolio strategies and investigate the impact of regime switching on optimal portfolio strategies.

Effective action for stochastic partial differential equations

International Nuclear Information System (INIS)

Hochberg, David; Molina-Paris, Carmen; Perez-Mercader, Juan; Visser, Matt

1999-01-01

Stochastic partial differential equations (SPDEs) are the basic tool for modeling systems where noise is important. SPDEs are used for models of turbulence, pattern formation, and the structural development of the universe itself. It is reasonably well known that certain SPDEs can be manipulated to be equivalent to (nonquantum) field theories that nevertheless exhibit deep and important relationships with quantum field theory. In this paper we systematically extend these ideas: We set up a functional integral formalism and demonstrate how to extract all the one-loop physics for an arbitrary SPDE subject to arbitrary Gaussian noise. It is extremely important to realize that Gaussian noise does not imply that the field variables undergo Gaussian fluctuations, and that these nonquantum field theories are fully interacting. The limitation to one loop is not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught us that one-loop physics is often quite adequate to give a good description of the salient issues. The limitation to one loop does, however, offer marked technical advantages: Because at one loop almost any field theory can be rendered finite using zeta function technology, we can sidestep the complications inherent in the Martin-Siggia-Rose formalism (the SPDE analog of the Becchi-Rouet-Stora-Tyutin formalism used in QFT) and instead focus attention on a minimalist approach that uses only the physical fields (this ''direct approach'' is the SPDE analog of canonical quantization using physical fields). After setting up the general formalism for the characteristic functional (partition function), we show how to define the effective action to all loops, and then focus on the one-loop effective action and its specialization to constant fields: the effective potential. The physical interpretation of the effective action and effective potential for SPDEs is addressed and we show that key features carry over from QFT to the case of

Effective action for stochastic partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Hochberg, David [Laboratorio de Astrofisica Espacial y Fisica Fundamental, Apartado 50727, 28080 Madrid, (Spain); Centro de Astrobiologia, INTA, Carratera Ajalvir, Km. 4, 28850 Torrejon, Madrid, (Spain); Molina-Paris, Carmen [Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (United States); Perez-Mercader, Juan [Laboratorio de Astrofisica Espacial y Fisica Fundamental, Apartado 50727, 28080 Madrid, (Spain); Visser, Matt [Physics Department, Washington University, Saint Louis, Missouri 63130-4899 (United States)

1999-12-01

Stochastic partial differential equations (SPDEs) are the basic tool for modeling systems where noise is important. SPDEs are used for models of turbulence, pattern formation, and the structural development of the universe itself. It is reasonably well known that certain SPDEs can be manipulated to be equivalent to (nonquantum) field theories that nevertheless exhibit deep and important relationships with quantum field theory. In this paper we systematically extend these ideas: We set up a functional integral formalism and demonstrate how to extract all the one-loop physics for an arbitrary SPDE subject to arbitrary Gaussian noise. It is extremely important to realize that Gaussian noise does not imply that the field variables undergo Gaussian fluctuations, and that these nonquantum field theories are fully interacting. The limitation to one loop is not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught us that one-loop physics is often quite adequate to give a good description of the salient issues. The limitation to one loop does, however, offer marked technical advantages: Because at one loop almost any field theory can be rendered finite using zeta function technology, we can sidestep the complications inherent in the Martin-Siggia-Rose formalism (the SPDE analog of the Becchi-Rouet-Stora-Tyutin formalism used in QFT) and instead focus attention on a minimalist approach that uses only the physical fields (this ''direct approach'' is the SPDE analog of canonical quantization using physical fields). After setting up the general formalism for the characteristic functional (partition function), we show how to define the effective action to all loops, and then focus on the one-loop effective action and its specialization to constant fields: the effective potential. The physical interpretation of the effective action and effective potential for SPDEs is addressed and we show that key features carry over from

Stochastic nonlinear beam equations

Czech Academy of Sciences Publication Activity Database

Brzezniak, Z.; Maslowski, Bohdan; Seidler, Jan

2005-01-01

Roč. 132, č. 1 (2005), s. 119-149 ISSN 0178-8051 R&D Projects: GA ČR(CZ) GA201/01/1197 Institutional research plan: CEZ:AV0Z10190503 Keywords : stochastic beam equation * stability Subject RIV: BA - General Mathematics Impact factor: 0.896, year: 2005

Comparison of Control Approaches in Genetic Regulatory Networks by Using Stochastic Master Equation Models, Probabilistic Boolean Network Models and Differential Equation Models and Estimated Error Analyzes

Science.gov (United States)

Caglar, Mehmet Umut; Pal, Ranadip

2011-03-01

Central dogma of molecular biology states that ``information cannot be transferred back from protein to either protein or nucleic acid''. However, this assumption is not exactly correct in most of the cases. There are a lot of feedback loops and interactions between different levels of systems. These types of interactions are hard to analyze due to the lack of cell level data and probabilistic - nonlinear nature of interactions. Several models widely used to analyze and simulate these types of nonlinear interactions. Stochastic Master Equation (SME) models give probabilistic nature of the interactions in a detailed manner, with a high calculation cost. On the other hand Probabilistic Boolean Network (PBN) models give a coarse scale picture of the stochastic processes, with a less calculation cost. Differential Equation (DE) models give the time evolution of mean values of processes in a highly cost effective way. The understanding of the relations between the predictions of these models is important to understand the reliability of the simulations of genetic regulatory networks. In this work the success of the mapping between SME, PBN and DE models is analyzed and the accuracy and affectivity of the control policies generated by using PBN and DE models is compared.

Modeling of long-range memory processes with inverse cubic distributions by the nonlinear stochastic differential equations

Science.gov (United States)

Kaulakys, B.; Alaburda, M.; Ruseckas, J.

2016-05-01

A well-known fact in the financial markets is the so-called ‘inverse cubic law’ of the cumulative distributions of the long-range memory fluctuations of market indicators such as a number of events of trades, trading volume and the logarithmic price change. We propose the nonlinear stochastic differential equation (SDE) giving both the power-law behavior of the power spectral density and the long-range dependent inverse cubic law of the cumulative distribution. This is achieved using the suggestion that when the market evolves from calm to violent behavior there is a decrease of the delay time of multiplicative feedback of the system in comparison to the driving noise correlation time. This results in a transition from the Itô to the Stratonovich sense of the SDE and yields a long-range memory process.

Fitting Nonlinear Ordinary Differential Equation Models with Random Effects and Unknown Initial Conditions Using the Stochastic Approximation Expectation-Maximization (SAEM) Algorithm.

Science.gov (United States)

Chow, Sy-Miin; Lu, Zhaohua; Sherwood, Andrew; Zhu, Hongtu

2016-03-01

The past decade has evidenced the increased prevalence of irregularly spaced longitudinal data in social sciences. Clearly lacking, however, are modeling tools that allow researchers to fit dynamic models to irregularly spaced data, particularly data that show nonlinearity and heterogeneity in dynamical structures. We consider the issue of fitting multivariate nonlinear differential equation models with random effects and unknown initial conditions to irregularly spaced data. A stochastic approximation expectation-maximization algorithm is proposed and its performance is evaluated using a benchmark nonlinear dynamical systems model, namely, the Van der Pol oscillator equations. The empirical utility of the proposed technique is illustrated using a set of 24-h ambulatory cardiovascular data from 168 men and women. Pertinent methodological challenges and unresolved issues are discussed.

Convergence and stability of the exponential Euler method for semi-linear stochastic delay differential equations.

Science.gov (United States)

Zhang, Ling

2017-01-01

The main purpose of this paper is to investigate the strong convergence and exponential stability in mean square of the exponential Euler method to semi-linear stochastic delay differential equations (SLSDDEs). It is proved that the exponential Euler approximation solution converges to the analytic solution with the strong order [Formula: see text] to SLSDDEs. On the one hand, the classical stability theorem to SLSDDEs is given by the Lyapunov functions. However, in this paper we study the exponential stability in mean square of the exact solution to SLSDDEs by using the definition of logarithmic norm. On the other hand, the implicit Euler scheme to SLSDDEs is known to be exponentially stable in mean square for any step size. However, in this article we propose an explicit method to show that the exponential Euler method to SLSDDEs is proved to share the same stability for any step size by the property of logarithmic norm.

Convergence and stability of the exponential Euler method for semi-linear stochastic delay differential equations

Directory of Open Access Journals (Sweden)

Ling Zhang

2017-10-01

Full Text Available Abstract The main purpose of this paper is to investigate the strong convergence and exponential stability in mean square of the exponential Euler method to semi-linear stochastic delay differential equations (SLSDDEs. It is proved that the exponential Euler approximation solution converges to the analytic solution with the strong order 1 2 $\\frac{1}{2}$ to SLSDDEs. On the one hand, the classical stability theorem to SLSDDEs is given by the Lyapunov functions. However, in this paper we study the exponential stability in mean square of the exact solution to SLSDDEs by using the definition of logarithmic norm. On the other hand, the implicit Euler scheme to SLSDDEs is known to be exponentially stable in mean square for any step size. However, in this article we propose an explicit method to show that the exponential Euler method to SLSDDEs is proved to share the same stability for any step size by the property of logarithmic norm.

Backward stochastic differential equations with two distinct reflecting barriers and quadratic growth generator

Directory of Open Access Journals (Sweden)

2006-01-01

Full Text Available We show the existence of a solution for the double-barrier reflected BSDE when the barriers are completely separate and the generator is continuous with quadratic growth. As an application, we solve the risk-sensitive mixed zero-sum stochastic differential game. In addition we deal with recallable options under Knightian uncertainty.

On the Langevin equation for stochastic quantization of gravity

International Nuclear Information System (INIS)

Nakazawa, Naohito.

1989-10-01

We study the Langevin equation for stochastic quantization of gravity. By introducing two independent variables with a second-class constraint for the gravitational field, we formulate a pair of the Langevin equations for gravity which couples with white noises. After eliminating the multiplier field for the second-class constraint, we show that the equations leads to stochastic quantization of gravity including an unique superspace metric. (author)

Quantum theory of open systems based on stochastic differential equations of generalized Langevin (non-Wiener) type

Energy Technology Data Exchange (ETDEWEB)

Basharov, A. M., E-mail: basharov@gmail.com [National Research Centre ' Kurchatov Institute,' (Russian Federation)

2012-09-15

It is shown that the effective Hamiltonian representation, as it is formulated in author's papers, serves as a basis for distinguishing, in a broadband environment of an open quantum system, independent noise sources that determine, in terms of the stationary quantum Wiener and Poisson processes in the Markov approximation, the effective Hamiltonian and the equation for the evolution operator of the open system and its environment. General stochastic differential equations of generalized Langevin (non-Wiener) type for the evolution operator and the kinetic equation for the density matrix of an open system are obtained, which allow one to analyze the dynamics of a wide class of localized open systems in the Markov approximation. The main distinctive features of the dynamics of open quantum systems described in this way are the stabilization of excited states with respect to collective processes and an additional frequency shift of the spectrum of the open system. As an illustration of the general approach developed, the photon dynamics in a single-mode cavity without losses on the mirrors is considered, which contains identical intracavity atoms coupled to the external vacuum electromagnetic field. For some atomic densities, the photons of the cavity mode are 'locked' inside the cavity, thus exhibiting a new phenomenon of radiation trapping and non-Wiener dynamics.

Non-cooperative stochastic differential game theory of generalized Markov jump linear systems

CERN Document Server

Zhang, Cheng-ke; Zhou, Hai-ying; Bin, Ning

2017-01-01

This book systematically studies the stochastic non-cooperative differential game theory of generalized linear Markov jump systems and its application in the field of finance and insurance. The book is an in-depth research book of the continuous time and discrete time linear quadratic stochastic differential game, in order to establish a relatively complete framework of dynamic non-cooperative differential game theory. It uses the method of dynamic programming principle and Riccati equation, and derives it into all kinds of existence conditions and calculating method of the equilibrium strategies of dynamic non-cooperative differential game. Based on the game theory method, this book studies the corresponding robust control problem, especially the existence condition and design method of the optimal robust control strategy. The book discusses the theoretical results and its applications in the risk control, option pricing, and the optimal investment problem in the field of finance and insurance, enriching the...

Set-Valued Stochastic Equation with Set-Valued Square Integrable Martingale

Directory of Open Access Journals (Sweden)

Li Jun-Gang

2017-01-01

Full Text Available In this paper, we shall introduce the stochastic integral of a stochastic process with respect to set-valued square integrable martingale. Then we shall give the Aumann integral measurable theorem, and give the set-valued stochastic Lebesgue integral and set-valued square integrable martingale integral equation. The existence and uniqueness of solution to set-valued stochastic integral equation are proved. The discussion will be useful in optimal control and mathematical finance in psychological factors.

Legendre Wavelet Operational Matrix Method for Solution of Riccati Differential Equation

Directory of Open Access Journals (Sweden)

S. Balaji

2014-01-01

Full Text Available A Legendre wavelet operational matrix method (LWM is presented for the solution of nonlinear fractional-order Riccati differential equations, having variety of applications in quantum chemistry and quantum mechanics. The fractional-order Riccati differential equations converted into a system of algebraic equations using Legendre wavelet operational matrix. Solutions given by the proposed scheme are more accurate and reliable and they are compared with recently developed numerical, analytical, and stochastic approaches. Comparison shows that the proposed LWM approach has a greater performance and less computational effort for getting accurate solutions. Further existence and uniqueness of the proposed problem are given and moreover the condition of convergence is verified.

Maximal stochastic transport in the Lorenz equations

Energy Technology Data Exchange (ETDEWEB)

Agarwal, Sahil, E-mail: sahil.agarwal@yale.edu [Program in Applied Mathematics, Yale University, New Haven (United States); Wettlaufer, J.S., E-mail: john.wettlaufer@yale.edu [Program in Applied Mathematics, Yale University, New Haven (United States); Departments of Geology & Geophysics, Mathematics and Physics, Yale University, New Haven (United States); Mathematical Institute, University of Oxford, Oxford (United Kingdom); Nordita, Royal Institute of Technology and Stockholm University, Stockholm (Sweden)

2016-01-08

We calculate the stochastic upper bounds for the Lorenz equations using an extension of the background method. In analogy with Rayleigh–Bénard convection the upper bounds are for heat transport versus Rayleigh number. As might be expected, the stochastic upper bounds are larger than the deterministic counterpart of Souza and Doering [1], but their variation with noise amplitude exhibits interesting behavior. Below the transition to chaotic dynamics the upper bounds increase monotonically with noise amplitude. However, in the chaotic regime this monotonicity depends on the number of realizations in the ensemble; at a particular Rayleigh number the bound may increase or decrease with noise amplitude. The origin of this behavior is the coupling between the noise and unstable periodic orbits, the degree of which depends on the degree to which the ensemble represents the ergodic set. This is confirmed by examining the close returns plots of the full solutions to the stochastic equations and the numerical convergence of the noise correlations. The numerical convergence of both the ensemble and time averages of the noise correlations is sufficiently slow that it is the limiting aspect of the realization of these bounds. Finally, we note that the full solutions of the stochastic equations demonstrate that the effect of noise is equivalent to the effect of chaos.

Stochastic Galerkin methods for the steady-state Navier–Stokes equations

Energy Technology Data Exchange (ETDEWEB)

Sousedík, Bedřich, E-mail: sousedik@umbc.edu [Department of Mathematics and Statistics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250 (United States); Elman, Howard C., E-mail: elman@cs.umd.edu [Department of Computer Science and Institute for Advanced Computer Studies, University of Maryland, College Park, MD 20742 (United States)

2016-07-01

We study the steady-state Navier–Stokes equations in the context of stochastic finite element discretizations. Specifically, we assume that the viscosity is a random field given in the form of a generalized polynomial chaos expansion. For the resulting stochastic problem, we formulate the model and linearization schemes using Picard and Newton iterations in the framework of the stochastic Galerkin method, and we explore properties of the resulting stochastic solutions. We also propose a preconditioner for solving the linear systems of equations arising at each step of the stochastic (Galerkin) nonlinear iteration and demonstrate its effectiveness for solving a set of benchmark problems.

Differential Equations Compatible with KZ Equations

International Nuclear Information System (INIS)

Felder, G.; Markov, Y.; Tarasov, V.; Varchenko, A.

2000-01-01

We define a system of 'dynamical' differential equations compatible with the KZ differential equations. The KZ differential equations are associated to a complex simple Lie algebra g. These are equations on a function of n complex variables z i taking values in the tensor product of n finite dimensional g-modules. The KZ equations depend on the 'dual' variable in the Cartan subalgebra of g. The dynamical differential equations are differential equations with respect to the dual variable. We prove that the standard hypergeometric solutions of the KZ equations also satisfy the dynamical equations. As an application we give a new determinant formula for the coordinates of a basis of hypergeometric solutions

STOCHASTIC FLOWS OF MAPPINGS

Institute of Scientific and Technical Information of China (English)

无

2007-01-01

In this paper, the stochastic flow of mappings generated by a Feller convolution semigroup on a compact metric space is studied. This kind of flow is the generalization of superprocesses of stochastic flows and stochastic diffeomorphism induced by the strong solutions of stochastic differential equations.

Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE

Science.gov (United States)

Ansmann, Gerrit

2018-04-01

We present a family of Python modules for the numerical integration of ordinary, delay, or stochastic differential equations. The key features are that the user enters the derivative symbolically and it is just-in-time-compiled, allowing the user to efficiently integrate differential equations from a higher-level interpreted language. The presented modules are particularly suited for large systems of differential equations such as those used to describe dynamics on complex networks. Through the selected method of input, the presented modules also allow almost complete automatization of the process of estimating regular as well as transversal Lyapunov exponents for ordinary and delay differential equations. We conceptually discuss the modules' design, analyze their performance, and demonstrate their capabilities by application to timely problems.

Differential equations

CERN Document Server

Barbu, Viorel

2016-01-01

This textbook is a comprehensive treatment of ordinary differential equations, concisely presenting basic and essential results in a rigorous manner. Including various examples from physics, mechanics, natural sciences, engineering and automatic theory, Differential Equations is a bridge between the abstract theory of differential equations and applied systems theory. Particular attention is given to the existence and uniqueness of the Cauchy problem, linear differential systems, stability theory and applications to first-order partial differential equations. Upper undergraduate students and researchers in applied mathematics and systems theory with a background in advanced calculus will find this book particularly useful. Supplementary topics are covered in an appendix enabling the book to be completely self-contained.

Derivation of the Schroedinger equation from stochastic mechanics

International Nuclear Information System (INIS)

Wallstrom, T.C.

1988-01-01

The thesis is divided into four largely independent chapters. The first three chapters treat mathematical problems in the theory of stochastic mechanics. The fourth chapter deals with stochastic mechanisms as a physical theory and shows that the Schroedinger equation cannot be derived from existing formulations of stochastic mechanics, as had previously been believed. Since the drift coefficients of stochastic mechanical diffusions are undefined on the nodes, or zeros of the density, an important problem has been to show that the sample paths stay away from the nodes. In Chapter 1, it is shown that for a smooth wavefunction, the closest approach to the nodes can be bounded solely in terms of the time-integrated energy. The ergodic properties of stochastic mechanical diffusions are greatly complicated by the tendency of the particles to avoid the nodes. In Chapter 2, it is shown that a sufficient condition for a stationary process to be ergodic is that there exist positive t and c such that for all x and y, p t (x,y) > cp(y), and this result is applied to show that the set of spin-1/2 diffusions is uniformly ergodic. Nelson has conjectured that in the limit as the particle's moment of inertia I goes to zero, the projections of the Bopp-Haag-Dankel diffusions onto IR 3 converge to a Markovian limit process. This conjecture is proved for the spin-1/2 case in Chapter 3, and the limit process identified as the diffusion naturally associated with the solution to the regular Pauli equation. In Chapter 4 it is shown that the general solution of the stochastic Newton equation does not correspond to a solution of the Schroedinger equation

Stochastic processes in cell biology

CERN Document Server

Bressloff, Paul C

2014-01-01

This book develops the theory of continuous and discrete stochastic processes within the context of cell biology.  A wide range of biological topics are covered including normal and anomalous diffusion in complex cellular environments, stochastic ion channels and excitable systems, stochastic calcium signaling, molecular motors, intracellular transport, signal transduction, bacterial chemotaxis, robustness in gene networks, genetic switches and oscillators, cell polarization, polymerization, cellular length control, and branching processes. The book also provides a pedagogical introduction to the theory of stochastic process – Fokker Planck equations, stochastic differential equations, master equations and jump Markov processes, diffusion approximations and the system size expansion, first passage time problems, stochastic hybrid systems, reaction-diffusion equations, exclusion processes, WKB methods, martingales and branching processes, stochastic calculus, and numerical methods.   This text is primarily...

Almost Periodic Solutions for Impulsive Fractional Stochastic Evolution Equations

Directory of Open Access Journals (Sweden)

Toufik Guendouzi

2014-08-01

Full Text Available In this paper, we consider the existence of square-mean piecewise almost periodic solutions for impulsive fractional stochastic evolution equations involving Caputo fractional derivative. The main results are obtained by means of the theory of operators semi-group, fractional calculus, fixed point technique and stochastic analysis theory and methods adopted directly from deterministic fractional equations. Some known results are improved and generalized.

Differential equations a dynamical systems approach ordinary differential equations

CERN Document Server

Hubbard, John H

1991-01-01

This is a corrected third printing of the first part of the text Differential Equations: A Dynamical Systems Approach written by John Hubbard and Beverly West. The authors' main emphasis in this book is on ordinary differential equations. The book is most appropriate for upper level undergraduate and graduate students in the fields of mathematics, engineering, and applied mathematics, as well as the life sciences, physics and economics. Traditional courses on differential equations focus on techniques leading to solutions. Yet most differential equations do not admit solutions which can be written in elementary terms. The authors have taken the view that a differential equations defines functions; the object of the theory is to understand the behavior of these functions. The tools the authors use include qualitative and numerical methods besides the traditional analytic methods. The companion software, MacMath, is designed to bring these notions to life.

Partial Differential Equations

CERN Document Server

1988-01-01

The volume contains a selection of papers presented at the 7th Symposium on differential geometry and differential equations (DD7) held at the Nankai Institute of Mathematics, Tianjin, China, in 1986. Most of the contributions are original research papers on topics including elliptic equations, hyperbolic equations, evolution equations, non-linear equations from differential geometry and mechanics, micro-local analysis.

Quantum stochastic calculus associated with quadratic quantum noises

International Nuclear Information System (INIS)

Ji, Un Cig; Sinha, Kalyan B.

2016-01-01

We first study a class of fundamental quantum stochastic processes induced by the generators of a six dimensional non-solvable Lie †-algebra consisting of all linear combinations of the generalized Gross Laplacian and its adjoint, annihilation operator, creation operator, conservation, and time, and then we study the quantum stochastic integrals associated with the class of fundamental quantum stochastic processes, and the quantum Itô formula is revisited. The existence and uniqueness of solution of a quantum stochastic differential equation is proved. The unitarity conditions of solutions of quantum stochastic differential equations associated with the fundamental processes are examined. The quantum stochastic calculus extends the Hudson-Parthasarathy quantum stochastic calculus

Quantum stochastic calculus associated with quadratic quantum noises

Energy Technology Data Exchange (ETDEWEB)

Ji, Un Cig, E-mail: uncigji@chungbuk.ac.kr [Department of Mathematics, Research Institute of Mathematical Finance, Chungbuk National University, Cheongju, Chungbuk 28644 (Korea, Republic of); Sinha, Kalyan B., E-mail: kbs-jaya@yahoo.co.in [Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore-64, India and Department of Mathematics, Indian Institute of Science, Bangalore-12 (India)

2016-02-15

We first study a class of fundamental quantum stochastic processes induced by the generators of a six dimensional non-solvable Lie †-algebra consisting of all linear combinations of the generalized Gross Laplacian and its adjoint, annihilation operator, creation operator, conservation, and time, and then we study the quantum stochastic integrals associated with the class of fundamental quantum stochastic processes, and the quantum Itô formula is revisited. The existence and uniqueness of solution of a quantum stochastic differential equation is proved. The unitarity conditions of solutions of quantum stochastic differential equations associated with the fundamental processes are examined. The quantum stochastic calculus extends the Hudson-Parthasarathy quantum stochastic calculus.

Partial differential equation methods for stochastic dynamic optimization: an application to wind power generation with energy storage.

Science.gov (United States)

Johnson, Paul; Howell, Sydney; Duck, Peter

2017-08-13

A mixed financial/physical partial differential equation (PDE) can optimize the joint earnings of a single wind power generator (WPG) and a generic energy storage device (ESD). Physically, the PDE includes constraints on the ESD's capacity, efficiency and maximum speeds of charge and discharge. There is a mean-reverting daily stochastic cycle for WPG power output. Physically, energy can only be produced or delivered at finite rates. All suppliers must commit hourly to a finite rate of delivery C , which is a continuous control variable that is changed hourly. Financially, we assume heavy 'system balancing' penalties in continuous time, for deviations of output rate from the commitment C Also, the electricity spot price follows a mean-reverting stochastic cycle with a strong evening peak, when system balancing penalties also peak. Hence the economic goal of the WPG plus ESD, at each decision point, is to maximize expected net present value (NPV) of all earnings (arbitrage) minus the NPV of all expected system balancing penalties, along all financially/physically feasible future paths through state space. Given the capital costs for the various combinations of the physical parameters, the design and operating rules for a WPG plus ESD in a finite market may be jointly optimizable.This article is part of the themed issue 'Energy management: flexibility, risk and optimization'. © 2017 The Author(s).

Perturbations of linear delay differential equations at the verge of instability.

Science.gov (United States)

Lingala, N; Namachchivaya, N Sri

2016-06-01

The characteristic equation for a linear delay differential equation (DDE) has countably infinite roots on the complex plane. This paper considers linear DDEs that are on the verge of instability, i.e., a pair of roots of the characteristic equation lies on the imaginary axis of the complex plane and all other roots have negative real parts. It is shown that when small noise perturbations are present, the probability distribution of the dynamics can be approximated by the probability distribution of a certain one-dimensional stochastic differential equation (SDE) without delay. This is advantageous because equations without delay are easier to simulate and one-dimensional SDEs are analytically tractable. When the perturbations are also linear, it is shown that the stability depends on a specific complex number. The theory is applied to study oscillators with delayed feedback. Some errors in other articles that use multiscale approach are pointed out.

Stochastic analysis of complex reaction networks using binomial moment equations.

Science.gov (United States)

Barzel, Baruch; Biham, Ofer

2012-09-01

The stochastic analysis of complex reaction networks is a difficult problem because the number of microscopic states in such systems increases exponentially with the number of reactive species. Direct integration of the master equation is thus infeasible and is most often replaced by Monte Carlo simulations. While Monte Carlo simulations are a highly effective tool, equation-based formulations are more amenable to analytical treatment and may provide deeper insight into the dynamics of the network. Here, we present a highly efficient equation-based method for the analysis of stochastic reaction networks. The method is based on the recently introduced binomial moment equations [Barzel and Biham, Phys. Rev. Lett. 106, 150602 (2011)]. The binomial moments are linear combinations of the ordinary moments of the probability distribution function of the population sizes of the interacting species. They capture the essential combinatorics of the reaction processes reflecting their stoichiometric structure. This leads to a simple and transparent form of the equations, and allows a highly efficient and surprisingly simple truncation scheme. Unlike ordinary moment equations, in which the inclusion of high order moments is prohibitively complicated, the binomial moment equations can be easily constructed up to any desired order. The result is a set of equations that enables the stochastic analysis of complex reaction networks under a broad range of conditions. The number of equations is dramatically reduced from the exponential proliferation of the master equation to a polynomial (and often quadratic) dependence on the number of reactive species in the binomial moment equations. The aim of this paper is twofold: to present a complete derivation of the binomial moment equations; to demonstrate the applicability of the moment equations for a representative set of example networks, in which stochastic effects play an important role.

Quantum Ito's formula and stochastic evolutions

International Nuclear Information System (INIS)

Hudson, R.L.; Parthasarathy, K.R.

1984-01-01

Using only the Boson canonical commutation relations and the Riemann-Lebesgue integral we construct a simple theory of stochastic integrals and differentials with respect to the basic field operator processes. This leads to a noncommutative Ito product formula, a realisation of the classical Poisson process in Fock space which gives a noncommutative central limit theorem, the construction of solutions of certain noncommutative stochastic differential equations, and finally to the integration of certain irreversible equations of motion governed by semigroups of completely positive maps. The classical Ito product formula for stochastic differentials with respect to Brownian motion and the Poisson process is a special case. (orig.)

Renormalization group equations in the stochastic quantization scheme

International Nuclear Information System (INIS)

Pugnetti, S.

1987-01-01

We show that there exists a remarkable link between the stochastic quantization and the theory of critical phenomena and dynamical statistical systems. In the stochastic quantization of a field theory, the stochastic Green functions coverge to the quantum ones when the frictious time goes to infinity. We therefore use the typical techniques of the Renormalization Group equations developed in the framework of critical phenomena to discuss some features of the convergence of the stochastic theory. We are also able, in this way, to compute some dynamical critical exponents and give new numerical valuations for them. (orig.)

Quantum cybernetics: a new perspective for Nelson's stochastic theory, nonlocality, and the Klein-Gordon equation

Science.gov (United States)

Grössing, Gerhard

2002-04-01

The Klein-Gordon equation is shown to be equivalent to coupled partial differential equations for a sub-quantum Brownian movement of a “particle”, which is both passively affected by, and actively affecting, a diffusion process of its generally nonlocal environment. This indicates circularly causal, or “cybernetic”, relationships between “particles” and their surroundings. Moreover, in the relativistic domain, the original stochastic theory of Nelson is shown to hold as a limiting case only, i.e., for a vanishing quantum potential.

Modelling the heat dynamics of a monitored Test Reference Environment for Building Integrated Photovoltaic systems using stochastic differential equations

DEFF Research Database (Denmark)

Lodi, C.; Bacher, Peder; Cipriano, J.

2012-01-01

reduce the ventilation thermal losses of the building by pre-heating the fresh air. Furthermore, by decreasing PV module temperature, the ventilation air heat extraction can simultaneously increase electrical and thermal energy production of the building. A correct prediction of the PV module temperature...... and heat transfer coefficients is fundamental in order to improve the thermo-electrical production.The considered grey-box models are composed of a set of continuous time stochastic differential equations, holding the physical description of the system, combined with a set of discrete time measurement......This paper deals with grey-box modelling of the energy transfer of a double skin Building Integrated Photovoltaic (BIPV) system. Grey-box models are based on a combination of prior physical knowledge and statistics, which enable identification of the unknown parameters in the system and accurate...

A stochastic model of multiple scattering of charged particles: process, transport equation and solutions

International Nuclear Information System (INIS)

Papiez, L.; Moskvin, V.; Tulovsky, V.

2001-01-01

The process of angular-spatial evolution of multiple scattering of charged particles can be described by a special case of Boltzmann integro-differential equation called Lewis equation. The underlying stochastic process for this evolution is the compound Poisson process on the surface of the unit sphere. The significant portion of events that constitute compound Poisson process that describes multiple scattering have diffusional character. This property allows to analyze the process of angular-spatial evolution of multiple scattering of charged particles as combination of soft and hard collision processes and compute appropriately its transition densities. These computations provide a method of the approximate solution to the Lewis equation. (orig.)

Polynomial chaos methods for hyperbolic partial differential equations numerical techniques for fluid dynamics problems in the presence of uncertainties

CERN Document Server

Pettersson, Mass Per; Nordström, Jan

2015-01-01

This monograph presents computational techniques and numerical analysis to study conservation laws under uncertainty using the stochastic Galerkin formulation. With the continual growth of computer power, these methods are becoming increasingly popular as an alternative to more classical sampling-based techniques. The approach described in the text takes advantage of stochastic Galerkin projections applied to the original conservation laws to produce a large system of modified partial differential equations, the solutions to which directly provide a full statistical characterization of the effect of uncertainties. Polynomial Chaos Methods of Hyperbolic Partial Differential Equations focuses on the analysis of stochastic Galerkin systems obtained for linear and non-linear convection-diffusion equations and for a systems of conservation laws; a detailed well-posedness and accuracy analysis is presented to enable the design of robust and stable numerical methods. The exposition is restricted to one spatial dime...

Stochastic optimal control in infinite dimension dynamic programming and HJB equations

CERN Document Server

Fabbri, Giorgio; Święch, Andrzej

2017-01-01

Providing an introduction to stochastic optimal control in infinite dimension, this book gives a complete account of the theory of second-order HJB equations in infinite-dimensional Hilbert spaces, focusing on its applicability to associated stochastic optimal control problems. It features a general introduction to optimal stochastic control, including basic results (e.g. the dynamic programming principle) with proofs, and provides examples of applications. A complete and up-to-date exposition of the existing theory of viscosity solutions and regular solutions of second-order HJB equations in Hilbert spaces is given, together with an extensive survey of other methods, with a full bibliography. In particular, Chapter 6, written by M. Fuhrman and G. Tessitore, surveys the theory of regular solutions of HJB equations arising in infinite-dimensional stochastic control, via BSDEs. The book is of interest to both pure and applied researchers working in the control theory of stochastic PDEs, and in PDEs in infinite ...

Different seeds to solve the equations of stochastic point kinetics using the Euler-Maruyama method

International Nuclear Information System (INIS)

Suescun D, D.; Oviedo T, M.

2017-09-01

In this paper, a numerical study of stochastic differential equations that describe the kinetics in a nuclear reactor is presented. These equations, known as the stochastic equations of punctual kinetics they model temporal variations in neutron population density and concentrations of deferred neutron precursors. Because these equations are probabilistic in nature (since random oscillations in the neutrons and population of precursors were considered to be approximately normally distributed, and these equations also possess strong coupling and stiffness properties) the proposed method for the numerical simulations is the Euler-Maruyama scheme that provides very good approximations for calculating the neutron population and concentrations of deferred neutron precursors. The method proposed for this work was computationally tested for different seeds, initial conditions, experimental data and forms of reactivity for a group of precursors and then for six groups of deferred neutron precursors at each time step with 5000 Brownian movements per seed. In a paper reported in the literature, the Euler-Maruyama method was proposed, but there are many doubts about the reported values, in addition to not reporting the seed used, so in this work is expected to rectify the reported values. After taking the average of the different seeds used to generate the pseudo-random numbers the results provided by the Euler-Maruyama scheme will be compared in mean and standard deviation with other methods reported in the literature and results of the deterministic model of the equations of the punctual kinetics. This comparison confirms in particular that the Euler-Maruyama scheme is an efficient method to solve the equations of stochastic point kinetics but different from the values found and reported by another author. The Euler-Maruyama method is simple and easy to implement, provides acceptable results for neutron population density and concentration of deferred neutron precursors and

Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation

KAUST Repository

Mélykúti, Bence

2010-01-01

The Chemical Langevin Equation (CLE), which is a stochastic differential equation driven by a multidimensional Wiener process, acts as a bridge between the discrete stochastic simulation algorithm and the deterministic reaction rate equation when simulating (bio)chemical kinetics. The CLE model is valid in the regime where molecular populations are abundant enough to assume their concentrations change continuously, but stochastic fluctuations still play a major role. The contribution of this work is that we observe and explore that the CLE is not a single equation, but a parametric family of equations, all of which give the same finite-dimensional distribution of the variables. On the theoretical side, we prove that as many Wiener processes are sufficient to formulate the CLE as there are independent variables in the equation, which is just the rank of the stoichiometric matrix. On the practical side, we show that in the case where there are m1 pairs of reversible reactions and m2 irreversible reactions there is another, simple formulation of the CLE with only m1 + m2 Wiener processes, whereas the standard approach uses 2 m1 + m2. We demonstrate that there are considerable computational savings when using this latter formulation. Such transformations of the CLE do not cause a loss of accuracy and are therefore distinct from model reduction techniques. We illustrate our findings by considering alternative formulations of the CLE for a human ether a-go-go related gene ion channel model and the Goldbeter-Koshland switch. © 2010 American Institute of Physics.

A Proposed Stochastic Finite Difference Approach Based on Homogenous Chaos Expansion

Directory of Open Access Journals (Sweden)

O. H. Galal

2013-01-01

Full Text Available This paper proposes a stochastic finite difference approach, based on homogenous chaos expansion (SFDHC. The said approach can handle time dependent nonlinear as well as linear systems with deterministic or stochastic initial and boundary conditions. In this approach, included stochastic parameters are modeled as second-order stochastic processes and are expanded using Karhunen-Loève expansion, while the response function is approximated using homogenous chaos expansion. Galerkin projection is used in converting the original stochastic partial differential equation (PDE into a set of coupled deterministic partial differential equations and then solved using finite difference method. Two well-known equations were used for efficiency validation of the method proposed. First one being the linear diffusion equation with stochastic parameter and the second is the nonlinear Burger's equation with stochastic parameter and stochastic initial and boundary conditions. In both of these examples, the probability distribution function of the response manifested close conformity to the results obtained from Monte Carlo simulation with optimized computational cost.

Stability of Equilibrium Points of Fractional Difference Equations with Stochastic Perturbations

Directory of Open Access Journals (Sweden)

Shaikhet Leonid

2008-01-01

Full Text Available It is supposed that the fractional difference equation , has an equilibrium point and is exposed to additive stochastic perturbations type of that are directly proportional to the deviation of the system state from the equilibrium point . It is shown that known results in the theory of stability of stochastic difference equations that were obtained via V. Kolmanovskii and L. Shaikhet general method of Lyapunov functionals construction can be successfully used for getting of sufficient conditions for stability in probability of equilibrium points of the considered stochastic fractional difference equation. Numerous graphical illustrations of stability regions and trajectories of solutions are plotted.

Fractional Stochastic Field Theory

Science.gov (United States)

Honkonen, Juha

2018-02-01

Models describing evolution of physical, chemical, biological, social and financial processes are often formulated as differential equations with the understanding that they are large-scale equations for averages of quantities describing intrinsically random processes. Explicit account of randomness may lead to significant changes in the asymptotic behaviour (anomalous scaling) in such models especially in low spatial dimensions, which in many cases may be captured with the use of the renormalization group. Anomalous scaling and memory effects may also be introduced with the use of fractional derivatives and fractional noise. Construction of renormalized stochastic field theory with fractional derivatives and fractional noise in the underlying stochastic differential equations and master equations and the interplay between fluctuation-induced and built-in anomalous scaling behaviour is reviewed and discussed.

Random ordinary differential equations and their numerical solution

CERN Document Server

Han, Xiaoying

2017-01-01

This book is intended to make recent results on the derivation of higher order numerical schemes for random ordinary differential equations (RODEs) available to a broader readership, and to familiarize readers with RODEs themselves as well as the closely associated theory of random dynamical systems. In addition, it demonstrates how RODEs are being used in the biological sciences, where non-Gaussian and bounded noise are often more realistic than the Gaussian white noise in stochastic differential equations (SODEs).   RODEs are used in many important applications and play a fundamental role in the theory of random dynamical systems.  They can be analyzed pathwise with deterministic calculus, but require further treatment beyond that of classical ODE theory due to the lack of smoothness in their time variable. Although classical numerical schemes for ODEs can be used pathwise for RODEs, they rarely attain their traditional order since the solutions of RODEs do not have sufficient smoothness to have Taylor ...

Stochastic effects on the nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Flessas, G P; Leach, P G L; Yannacopoulos, A N

2004-01-01

The aim of this article is to provide a brief review of recent advances in the field of stochastic effects on the nonlinear Schroedinger equation. The article reviews rigorous and perturbative results. (review article)

Stochastic Analysis with Financial Applications

CERN Document Server

Kohatsu-Higa, Arturo; Sheu, Shuenn-Jyi

2011-01-01

Stochastic analysis has a variety of applications to biological systems as well as physical and engineering problems, and its applications to finance and insurance have bloomed exponentially in recent times. The goal of this book is to present a broad overview of the range of applications of stochastic analysis and some of its recent theoretical developments. This includes numerical simulation, error analysis, parameter estimation, as well as control and robustness properties for stochastic equations. This book also covers the areas of backward stochastic differential equations via the (non-li

Momentum Maps and Stochastic Clebsch Action Principles

Science.gov (United States)

Cruzeiro, Ana Bela; Holm, Darryl D.; Ratiu, Tudor S.

2018-01-01

We derive stochastic differential equations whose solutions follow the flow of a stochastic nonlinear Lie algebra operation on a configuration manifold. For this purpose, we develop a stochastic Clebsch action principle, in which the noise couples to the phase space variables through a momentum map. This special coupling simplifies the structure of the resulting stochastic Hamilton equations for the momentum map. In particular, these stochastic Hamilton equations collectivize for Hamiltonians that depend only on the momentum map variable. The Stratonovich equations are derived from the Clebsch variational principle and then converted into Itô form. In comparing the Stratonovich and Itô forms of the stochastic dynamical equations governing the components of the momentum map, we find that the Itô contraction term turns out to be a double Poisson bracket. Finally, we present the stochastic Hamiltonian formulation of the collectivized momentum map dynamics and derive the corresponding Kolmogorov forward and backward equations.

On the Stochastic Wave Equation with Nonlinear Damping

International Nuclear Information System (INIS)

Kim, Jong Uhn

2008-01-01

We discuss an initial boundary value problem for the stochastic wave equation with nonlinear damping. We establish the existence and uniqueness of a solution. Our method for the existence of pathwise solutions consists of regularization of the equation and data, the Galerkin approximation and an elementary measure-theoretic argument. We also prove the existence of an invariant measure when the equation has pure nonlinear damping

On solutions of neutral stochastic delay Volterra equations with singular kernels

Directory of Open Access Journals (Sweden)

Xiaotai Wu

2012-08-01

Full Text Available In this paper, existence, uniqueness and continuity of the adapted solutions for neutral stochastic delay Volterra equations with singular kernels are discussed. In addition, continuous dependence on the initial date is also investigated. Finally, stochastic Volterra equation with the kernel of fractional Brownian motion is studied to illustrate the effectiveness of our results.

Brownian motion, martingales, and stochastic calculus

CERN Document Server

Le Gall, Jean-François

2016-01-01

This book offers a rigorous and self-contained presentation of stochastic integration and stochastic calculus within the general framework of continuous semimartingales. The main tools of stochastic calculus, including Itô’s formula, the optional stopping theorem and Girsanov’s theorem, are treated in detail alongside many illustrative examples. The book also contains an introduction to Markov processes, with applications to solutions of stochastic differential equations and to connections between Brownian motion and partial differential equations. The theory of local times of semimartingales is discussed in the last chapter. Since its invention by Itô, stochastic calculus has proven to be one of the most important techniques of modern probability theory, and has been used in the most recent theoretical advances as well as in applications to other fields such as mathematical finance. Brownian Motion, Martingales, and Stochastic Calculus provides a strong theoretical background to the reader interested i...

Nonlinear differential equations

Energy Technology Data Exchange (ETDEWEB)

Dresner, L.

1988-01-01

This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics.

Nonlinear differential equations

International Nuclear Information System (INIS)

Dresner, L.

1988-01-01

This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics

Stochastic differential calculus for Gaussian and non-Gaussian noises: A critical review

Science.gov (United States)

Falsone, G.

2018-03-01

In this paper a review of the literature works devoted to the study of stochastic differential equations (SDEs) subjected to Gaussian and non-Gaussian white noises and to fractional Brownian noises is given. In these cases, particular attention must be paid in treating the SDEs because the classical rules of the differential calculus, as the Newton-Leibnitz one, cannot be applied or are applicable with many difficulties. Here all the principal approaches solving the SDEs are reported for any kind of noise, highlighting the negative and positive properties of each one and making the comparisons, where it is possible.

Differential equations for dummies

CERN Document Server

Holzner, Steven

2008-01-01

The fun and easy way to understand and solve complex equations Many of the fundamental laws of physics, chemistry, biology, and economics can be formulated as differential equations. This plain-English guide explores the many applications of this mathematical tool and shows how differential equations can help us understand the world around us. Differential Equations For Dummies is the perfect companion for a college differential equations course and is an ideal supplemental resource for other calculus classes as well as science and engineering courses. It offers step-by-step techniques, practical tips, numerous exercises, and clear, concise examples to help readers improve their differential equation-solving skills and boost their test scores.

Wong-Zakai approximations and attractors for stochastic reaction-diffusion equations on unbounded domains

Science.gov (United States)

Wang, Xiaohu; Lu, Kening; Wang, Bixiang

2018-01-01

In this paper, we study the Wong-Zakai approximations given by a stationary process via the Wiener shift and their associated long term behavior of the stochastic reaction-diffusion equation driven by a white noise. We first prove the existence and uniqueness of tempered pullback attractors for the Wong-Zakai approximations of stochastic reaction-diffusion equation. Then, we show that the attractors of Wong-Zakai approximations converges to the attractor of the stochastic reaction-diffusion equation for both additive and multiplicative noise.

EXPONENTIAL ERGODICITY FOR STOCHASTIC BURGERS AND 2D NAVIER-STOKES EQUATIONS

CERN Document Server

Goldys, B

2004-01-01

It is shown that transition measures of the stochastic Navier-Stokes equation in dimension 2 converge exponentially fast to the corresponding invariant measures in the distance of total variation. As a corollary we obtain the existence of spectral gap for a related semigroup obtained by a sort of ground state trasformation. Analogous results are proved for the stochastic Burgers equation.

Numerical studies of the stochastic Korteweg-de Vries equation

International Nuclear Information System (INIS)

Lin Guang; Grinberg, Leopold; Karniadakis, George Em

2006-01-01

We present numerical solutions of the stochastic Korteweg-de Vries equation for three cases corresponding to additive time-dependent noise, multiplicative space-dependent noise and a combination of the two. We employ polynomial chaos for discretization in random space, and discontinuous Galerkin and finite difference for discretization in physical space. The accuracy of the stochastic solutions is investigated by comparing the first two moments against analytical and Monte Carlo simulation results. Of particular interest is the interplay of spatial discretization error with the stochastic approximation error, which is examined for different orders of spatial and stochastic approximation

Inhomogeneous Linear Random Differential Equations with Mutual Correlations between Multiplicative, Additive and Initial-Value Terms

NARCIS (Netherlands)

Roerdink, J.B.T.M.

1981-01-01

The cumulant expansion for linear stochastic differential equations is extended to the general case in which the coefficient matrix, the inhomogeneous part and the initial condition are all random and, moreover, statistically interdependent. The expansion now involves not only the autocorrelation

Maximum principle for a stochastic delayed system involving terminal state constraints.

Science.gov (United States)

Wen, Jiaqiang; Shi, Yufeng

2017-01-01

We investigate a stochastic optimal control problem where the controlled system is depicted as a stochastic differential delayed equation; however, at the terminal time, the state is constrained in a convex set. We firstly introduce an equivalent backward delayed system depicted as a time-delayed backward stochastic differential equation. Then a stochastic maximum principle is obtained by virtue of Ekeland's variational principle. Finally, applications to a state constrained stochastic delayed linear-quadratic control model and a production-consumption choice problem are studied to illustrate the main obtained result.

Stochastic Linear Quadratic Optimal Control Problems

International Nuclear Information System (INIS)

Chen, S.; Yong, J.

2001-01-01

This paper is concerned with the stochastic linear quadratic optimal control problem (LQ problem, for short) for which the coefficients are allowed to be random and the cost functional is allowed to have a negative weight on the square of the control variable. Some intrinsic relations among the LQ problem, the stochastic maximum principle, and the (linear) forward-backward stochastic differential equations are established. Some results involving Riccati equation are discussed as well

Large deviations for solutions to stochastic recurrence equations under Kesten's condition

DEFF Research Database (Denmark)

Buraczewski, Dariusz; Damek, Ewa; Mikosch, Thomas Valentin

2013-01-01

In this paper we prove large deviations results for partial sums constructed from the solution to a stochastic recurrence equation. We assume Kesten’s condition [17] under which the solution of the stochastic recurrence equation has a marginal distribution with power law tails, while the noise...... sequence of the equations can have light tails. The results of the paper are analogs of those obtained by A.V. and S.V. Nagaev [21, 22] in the case of partial sums of iid random variables. In the latter case, the large deviation probabilities of the partial sums are essentially determined by the largest...... step size of the partial sum. For the solution to a stochastic recurrence equation, the magnitude of the large deviation probabilities is again given by the tail of the maximum summand, but the exact asymptotic tail behavior is also influenced by clusters of extreme values, due to dependencies...

Partial differential equations

CERN Document Server

Evans, Lawrence C

2010-01-01

This text gives a comprehensive survey of modern techniques in the theoretical study of partial differential equations (PDEs) with particular emphasis on nonlinear equations. The exposition is divided into three parts: representation formulas for solutions; theory for linear partial differential equations; and theory for nonlinear partial differential equations. Included are complete treatments of the method of characteristics; energy methods within Sobolev spaces; regularity for second-order elliptic, parabolic, and hyperbolic equations; maximum principles; the multidimensional calculus of variations; viscosity solutions of Hamilton-Jacobi equations; shock waves and entropy criteria for conservation laws; and, much more.The author summarizes the relevant mathematics required to understand current research in PDEs, especially nonlinear PDEs. While he has reworked and simplified much of the classical theory (particularly the method of characteristics), he primarily emphasizes the modern interplay between funct...

Feynman path integral related to stochastic schroedinger equation

International Nuclear Information System (INIS)

Belavkin, V.P.; Smolyanov, O.G.

1998-01-01

The derivation of the Schroedinger equation describing the continuous measurement process is presented. The representation of the solution of the stochastic Schroedinger equation for continuous measurements is obtained by means of the Feynman path integral. The connection with the heuristic approach to the description of continuous measurements is considered. The connection with the Senon paradox is established [ru

Gaussian approximations for stochastic systems with delay: Chemical Langevin equation and application to a Brusselator system

International Nuclear Information System (INIS)

Brett, Tobias; Galla, Tobias

2014-01-01

We present a heuristic derivation of Gaussian approximations for stochastic chemical reaction systems with distributed delay. In particular, we derive the corresponding chemical Langevin equation. Due to the non-Markovian character of the underlying dynamics, these equations are integro-differential equations, and the noise in the Gaussian approximation is coloured. Following on from the chemical Langevin equation, a further reduction leads to the linear-noise approximation. We apply the formalism to a delay variant of the celebrated Brusselator model, and show how it can be used to characterise noise-driven quasi-cycles, as well as noise-triggered spiking. We find surprisingly intricate dependence of the typical frequency of quasi-cycles on the delay period

Gaussian approximations for stochastic systems with delay: chemical Langevin equation and application to a Brusselator system.

Science.gov (United States)

Brett, Tobias; Galla, Tobias

2014-03-28

We present a heuristic derivation of Gaussian approximations for stochastic chemical reaction systems with distributed delay. In particular, we derive the corresponding chemical Langevin equation. Due to the non-Markovian character of the underlying dynamics, these equations are integro-differential equations, and the noise in the Gaussian approximation is coloured. Following on from the chemical Langevin equation, a further reduction leads to the linear-noise approximation. We apply the formalism to a delay variant of the celebrated Brusselator model, and show how it can be used to characterise noise-driven quasi-cycles, as well as noise-triggered spiking. We find surprisingly intricate dependence of the typical frequency of quasi-cycles on the delay period.

Treatment of constraints in the stochastic quantization method and covariantized Langevin equation

International Nuclear Information System (INIS)

Ikegami, Kenji; Kimura, Tadahiko; Mochizuki, Riuji

1993-01-01

We study the treatment of the constraints in the stochastic quantization method. We improve the treatment of the stochastic consistency condition proposed by Namiki et al. by suitably taking into account the Ito calculus. Then we obtain an improved Langevin equation and the Fokker-Planck equation which naturally leads to the correct path integral quantization of the constrained system as the stochastic equilibrium state. This treatment is applied to an O(N) non-linear σ model and it is shown that singular terms appearing in the improved Langevin equation cancel out the δ n (0) divergences in one loop order. We also ascertain that the above Langevin equation, rewritten in terms of independent variables, is actually equivalent to the one in the general-coordinate transformation covariant and vielbein-rotation invariant formalism. (orig.)

Time-ordered product expansions for computational stochastic system biology

International Nuclear Information System (INIS)

Mjolsness, Eric

2013-01-01

The time-ordered product framework of quantum field theory can also be used to understand salient phenomena in stochastic biochemical networks. It is used here to derive Gillespie’s stochastic simulation algorithm (SSA) for chemical reaction networks; consequently, the SSA can be interpreted in terms of Feynman diagrams. It is also used here to derive other, more general simulation and parameter-learning algorithms including simulation algorithms for networks of stochastic reaction-like processes operating on parameterized objects, and also hybrid stochastic reaction/differential equation models in which systems of ordinary differential equations evolve the parameters of objects that can also undergo stochastic reactions. Thus, the time-ordered product expansion can be used systematically to derive simulation and parameter-fitting algorithms for stochastic systems. (paper)

Differential Equation over Banach Algebra

OpenAIRE

Kleyn, Aleks

2018-01-01

In the book, I considered differential equations of order $1$ over Banach $D$-algebra: differential equation solved with respect to the derivative; exact differential equation; linear homogeneous equation. In noncommutative Banach algebra, initial value problem for linear homogeneous equation has infinitely many solutions.

Ordinary differential equations

CERN Document Server

Greenberg, Michael D

2014-01-01

Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps

Higher-order Solution of Stochastic Diffusion equation with Nonlinear Losses Using WHEP technique

KAUST Repository

El-Beltagy, Mohamed A.

2014-01-06

Using Wiener-Hermite expansion with perturbation (WHEP) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. The Wiener-Hermite expansion is the only known expansion that handles the white/colored noise exactly. The main statistics, such as the mean, covariance, and higher order statistical moments, can be calculated by simple formulae involving only the deterministic Wiener-Hermite coefficients. In this poster, the WHEP technique is used to solve the 2D diffusion equation with nonlinear losses and excited with white noise. The solution will be obtained numerically and will be validated and compared with the analytical solution that can be obtained from any symbolic mathematics package such as Mathematica.

Higher-order Solution of Stochastic Diffusion equation with Nonlinear Losses Using WHEP technique

KAUST Repository

El-Beltagy, Mohamed A.; Al-Mulla, Noah

2014-01-01

Using Wiener-Hermite expansion with perturbation (WHEP) technique in the solution of the stochastic partial differential equations (SPDEs) has the advantage of converting the problem to a system of deterministic equations that can be solved efficiently using the standard deterministic numerical methods [1]. The Wiener-Hermite expansion is the only known expansion that handles the white/colored noise exactly. The main statistics, such as the mean, covariance, and higher order statistical moments, can be calculated by simple formulae involving only the deterministic Wiener-Hermite coefficients. In this poster, the WHEP technique is used to solve the 2D diffusion equation with nonlinear losses and excited with white noise. The solution will be obtained numerically and will be validated and compared with the analytical solution that can be obtained from any symbolic mathematics package such as Mathematica.

Ordinary differential equations

CERN Document Server

Miller, Richard K

1982-01-01

Ordinary Differential Equations is an outgrowth of courses taught for a number of years at Iowa State University in the mathematics and the electrical engineering departments. It is intended as a text for a first graduate course in differential equations for students in mathematics, engineering, and the sciences. Although differential equations is an old, traditional, and well-established subject, the diverse backgrounds and interests of the students in a typical modern-day course cause problems in the selection and method of presentation of material. In order to compensate for this diversity,

Uncertain differential equations

CERN Document Server

Yao, Kai

2016-01-01

This book introduces readers to the basic concepts of and latest findings in the area of differential equations with uncertain factors. It covers the analytic method and numerical method for solving uncertain differential equations, as well as their applications in the field of finance. Furthermore, the book provides a number of new potential research directions for uncertain differential equation. It will be of interest to researchers, engineers and students in the fields of mathematics, information science, operations research, industrial engineering, computer science, artificial intelligence, automation, economics, and management science.

Diffusion approximations to the chemical master equation only have a consistent stochastic thermodynamics at chemical equilibrium.

Science.gov (United States)

Horowitz, Jordan M

2015-07-28

The stochastic thermodynamics of a dilute, well-stirred mixture of chemically reacting species is built on the stochastic trajectories of reaction events obtained from the chemical master equation. However, when the molecular populations are large, the discrete chemical master equation can be approximated with a continuous diffusion process, like the chemical Langevin equation or low noise approximation. In this paper, we investigate to what extent these diffusion approximations inherit the stochastic thermodynamics of the chemical master equation. We find that a stochastic-thermodynamic description is only valid at a detailed-balanced, equilibrium steady state. Away from equilibrium, where there is no consistent stochastic thermodynamics, we show that one can still use the diffusive solutions to approximate the underlying thermodynamics of the chemical master equation.

An introduction to probability and stochastic processes

CERN Document Server

Melsa, James L

2013-01-01

Geared toward college seniors and first-year graduate students, this text is designed for a one-semester course in probability and stochastic processes. Topics covered in detail include probability theory, random variables and their functions, stochastic processes, linear system response to stochastic processes, Gaussian and Markov processes, and stochastic differential equations. 1973 edition.

Environmental vs Demographic Stochasticity in Population Growth

OpenAIRE

Braumann, C. A.

2010-01-01

Compares the effect on population growth of envinonmental stochasticity (random environmental variations described by stochastic differential equations) with demographic stochasticity (random variations in births and deaths described by branching processes and birth-and-death processes), in the density-independent and the density-dependent cases.

Set-valued and fuzzy stochastic integral equations driven by semimartingales under Osgood condition

Directory of Open Access Journals (Sweden)

Malinowski Marek T.

2015-01-01

Full Text Available We analyze the set-valued stochastic integral equations driven by continuous semimartingales and prove the existence and uniqueness of solutions to such equations in the framework of the hyperspace of nonempty, bounded, convex and closed subsets of the Hilbert space L2 (consisting of square integrable random vectors. The coefficients of the equations are assumed to satisfy the Osgood type condition that is a generalization of the Lipschitz condition. Continuous dependence of solutions with respect to data of the equation is also presented. We consider equations driven by semimartingale Z and equations driven by processes A;M from decomposition of Z, where A is a process of finite variation and M is a local martingale. These equations are not equivalent. Finally, we show that the analysis of the set-valued stochastic integral equations can be extended to a case of fuzzy stochastic integral equations driven by semimartingales under Osgood type condition. To obtain our results we use the set-valued and fuzzy Maruyama type approximations and Bihari’s inequality.

QUANTUM STOCHASTIC PROCESSES: BOSON AND FERMION BROWNIAN MOTION

Directory of Open Access Journals (Sweden)

A.E.Kobryn

2003-01-01

Full Text Available Dynamics of quantum systems which are stochastically perturbed by linear coupling to the reservoir can be studied in terms of quantum stochastic differential equations (for example, quantum stochastic Liouville equation and quantum Langevin equation. In order to work it out one needs to define the quantum Brownian motion. As far as only its boson version has been known until recently, in the present paper we present the definition which makes it possible to consider the fermion Brownian motion as well.

A fixed point approach towards stability of delay differential equations with applications to neural networks

NARCIS (Netherlands)

Chen, Guiling

2013-01-01

This thesis studies asymptotic behavior and stability of determinsitic and stochastic delay differential equations. The approach used in this thesis is based on fixed point theory, which does not resort to any Liapunov function or Liapunov functional. The main contribution of this thesis is to study

Generalized differential transform method to differential-difference equation

International Nuclear Information System (INIS)

Zou Li; Wang Zhen; Zong Zhi

2009-01-01

In this Letter, we generalize the differential transform method to solve differential-difference equation for the first time. Two simple but typical examples are applied to illustrate the validity and the great potential of the generalized differential transform method in solving differential-difference equation. A Pade technique is also introduced and combined with GDTM in aim of extending the convergence area of presented series solutions. Comparisons are made between the results of the proposed method and exact solutions. Then we apply the differential transform method to the discrete KdV equation and the discrete mKdV equation, and successfully obtain solitary wave solutions. The results reveal that the proposed method is very effective and simple. We should point out that generalized differential transform method is also easy to be applied to other nonlinear differential-difference equation.

Introduction to partial differential equations

CERN Document Server

Greenspan, Donald

2000-01-01

Designed for use in a one-semester course by seniors and beginning graduate students, this rigorous presentation explores practical methods of solving differential equations, plus the unifying theory underlying the mathematical superstructure. Topics include basic concepts, Fourier series, second-order partial differential equations, wave equation, potential equation, heat equation, approximate solution of partial differential equations, and more. Exercises appear at the ends of most chapters. 1961 edition.

New analytic solutions of stochastic coupled KdV equations

International Nuclear Information System (INIS)

Dai Chaoqing; Chen Junlang

2009-01-01

In this paper, firstly, we use the exp-function method to seek new exact solutions of the Riccati equation. Then, with the help of Hermit transformation, we employ the Riccati equation and its new exact solutions to find new analytic solutions of the stochastic coupled KdV equation in the white noise environment. As some special examples, some analytic solutions can degenerate into these solutions reported in open literatures.

A Stochastic Delay Model For Pricing Debt And Loan Guarantees: Theoretical Results

OpenAIRE

Kemajou, Elisabeth; Mohammed, Salah-Eldin; Tambue, Antoine

2012-01-01

We consider that the price of a firm follows a non linear stochastic delay differential equation. We also assume that any claim value whose value depends on firm value and time follows a non linear stochastic delay differential equation. Using self-financed strategy and replication we are able to derive a Random Partial Differential Equation (RPDE) satisfied by any corporate claim whose value is a function of firm value and time. Under specific final and boundary conditions, we solve the RPDE...

On the well-posedness of the stochastic Allen–Cahn equation in two dimensions

International Nuclear Information System (INIS)

Ryser, Marc D.; Nigam, Nilima; Tupper, Paul F.

2012-01-01

White noise-driven nonlinear stochastic partial differential equations (SPDEs) of parabolic type are frequently used to model physical systems in space dimensions d = 1, 2, 3. Whereas existence and uniqueness of weak solutions to these equations are well established in one dimension, the situation is different for d ⩾ 2. Despite their popularity in the applied sciences, higher dimensional versions of these SPDE models are generally assumed to be ill-posed by the mathematics community. We study this discrepancy on the specific example of the two dimensional Allen–Cahn equation driven by additive white noise. Since it is unclear how to define the notion of a weak solution to this equation, we regularize the noise and introduce a family of approximations. Based on heuristic arguments and numerical experiments, we conjecture that these approximations exhibit divergent behavior in the continuum limit. The results strongly suggest that shrinking the mesh size in simulations of the two-dimensional white noise-driven Allen–Cahn equation does not lead to the recovery of a physically meaningful limit.

Representing Sudden Shifts in Intensive Dyadic Interaction Data Using Differential Equation Models with Regime Switching.

Science.gov (United States)

Chow, Sy-Miin; Ou, Lu; Ciptadi, Arridhana; Prince, Emily B; You, Dongjun; Hunter, Michael D; Rehg, James M; Rozga, Agata; Messinger, Daniel S

2018-06-01

A growing number of social scientists have turned to differential equations as a tool for capturing the dynamic interdependence among a system of variables. Current tools for fitting differential equation models do not provide a straightforward mechanism for diagnosing evidence for qualitative shifts in dynamics, nor do they provide ways of identifying the timing and possible determinants of such shifts. In this paper, we discuss regime-switching differential equation models, a novel modeling framework for representing abrupt changes in a system of differential equation models. Estimation was performed by combining the Kim filter (Kim and Nelson State-space models with regime switching: classical and Gibbs-sampling approaches with applications, MIT Press, Cambridge, 1999) and a numerical differential equation solver that can handle both ordinary and stochastic differential equations. The proposed approach was motivated by the need to represent discrete shifts in the movement dynamics of [Formula: see text] mother-infant dyads during the Strange Situation Procedure (SSP), a behavioral assessment where the infant is separated from and reunited with the mother twice. We illustrate the utility of a novel regime-switching differential equation model in representing children's tendency to exhibit shifts between the goal of staying close to their mothers and intermittent interest in moving away from their mothers to explore the room during the SSP. Results from empirical model fitting were supplemented with a Monte Carlo simulation study to evaluate the use of information criterion measures to diagnose sudden shifts in dynamics.

Basis adaptation and domain decomposition for steady-state partial differential equations with random coefficients

Energy Technology Data Exchange (ETDEWEB)

Tipireddy, R.; Stinis, P.; Tartakovsky, A. M.

2017-12-01

We present a novel approach for solving steady-state stochastic partial differential equations (PDEs) with high-dimensional random parameter space. The proposed approach combines spatial domain decomposition with basis adaptation for each subdomain. The basis adaptation is used to address the curse of dimensionality by constructing an accurate low-dimensional representation of the stochastic PDE solution (probability density function and/or its leading statistical moments) in each subdomain. Restricting the basis adaptation to a specific subdomain affords finding a locally accurate solution. Then, the solutions from all of the subdomains are stitched together to provide a global solution. We support our construction with numerical experiments for a steady-state diffusion equation with a random spatially dependent coefficient. Our results show that highly accurate global solutions can be obtained with significantly reduced computational costs.

Elements of partial differential equations

CERN Document Server

Sneddon, Ian Naismith

1957-01-01

Geared toward students of applied rather than pure mathematics, this volume introduces elements of partial differential equations. Its focus is primarily upon finding solutions to particular equations rather than general theory.Topics include ordinary differential equations in more than two variables, partial differential equations of the first and second orders, Laplace's equation, the wave equation, and the diffusion equation. A helpful Appendix offers information on systems of surfaces, and solutions to the odd-numbered problems appear at the end of the book. Readers pursuing independent st

Stochastic B-series and order conditions for exponential integrators

DEFF Research Database (Denmark)

Arara, Alemayehu Adugna; Debrabant, Kristian; Kværnø, Anne

2018-01-01

We discuss stochastic differential equations with a stiff linear part and their approximation by stochastic exponential integrators. Representing the exact and approximate solutions using B-series and rooted trees, we derive the order conditions for stochastic exponential integrators. The resulting...

Solving polynomial differential equations by transforming them to linear functional-differential equations

OpenAIRE

Nahay, John Michael

2008-01-01

We present a new approach to solving polynomial ordinary differential equations by transforming them to linear functional equations and then solving the linear functional equations. We will focus most of our attention upon the first-order Abel differential equation with two nonlinear terms in order to demonstrate in as much detail as possible the computations necessary for a complete solution. We mention in our section on further developments that the basic transformation idea can be generali...

Stochastic models for atmospheric dispersion

DEFF Research Database (Denmark)

Ditlevsen, Ove Dalager

2003-01-01

Simple stochastic differential equation models have been applied by several researchers to describe the dispersion of tracer particles in the planetary atmospheric boundary layer and to form the basis for computer simulations of particle paths. To obtain the drift coefficient, empirical vertical...... positions close to the boundaries. Different rules have been suggested in the literature with justifications based on simulation studies. Herein the relevant stochastic differential equation model is formulated in a particular way. The formulation is based on the marginal transformation of the position...... velocity distributions that depend on height above the ground both with respect to standard deviation and skewness are substituted into the stationary Fokker/Planck equation. The particle position distribution is taken to be uniform *the well/mixed condition( and also a given dispersion coefficient...

Introduction to differential equations

CERN Document Server

Taylor, Michael E

2011-01-01

The mathematical formulations of problems in physics, economics, biology, and other sciences are usually embodied in differential equations. The analysis of the resulting equations then provides new insight into the original problems. This book describes the tools for performing that analysis. The first chapter treats single differential equations, emphasizing linear and nonlinear first order equations, linear second order equations, and a class of nonlinear second order equations arising from Newton's laws. The first order linear theory starts with a self-contained presentation of the exponen

A primer on stochastic epidemic models: Formulation, numerical simulation, and analysis

Directory of Open Access Journals (Sweden)

Linda J.S. Allen

2017-05-01

Full Text Available Some mathematical methods for formulation and numerical simulation of stochastic epidemic models are presented. Specifically, models are formulated for continuous-time Markov chains and stochastic differential equations. Some well-known examples are used for illustration such as an SIR epidemic model and a host-vector malaria model. Analytical methods for approximating the probability of a disease outbreak are also discussed. Keywords: Branching process, Continuous-time Markov chain, Minor outbreak, Stochastic differential equation, 2000 MSC: 60H10, 60J28, 92D30

Non-Markovian stochastic Schroedinger equations: Generalization to real-valued noise using quantum-measurement theory

International Nuclear Information System (INIS)

Gambetta, Jay; Wiseman, H.M.

2002-01-01

Do stochastic Schroedinger equations, also known as unravelings, have a physical interpretation? In the Markovian limit, where the system on average obeys a master equation, the answer is yes. Markovian stochastic Schroedinger equations generate quantum trajectories for the system state conditioned on continuously monitoring the bath. For a given master equation, there are many different unravelings, corresponding to different sorts of measurement on the bath. In this paper we address the non-Markovian case, and in particular the sort of stochastic Schroedinger equation introduced by Strunz, Diosi, and Gisin [Phys. Rev. Lett. 82, 1801 (1999)]. Using a quantum-measurement theory approach, we rederive their unraveling that involves complex-valued Gaussian noise. We also derive an unraveling involving real-valued Gaussian noise. We show that in the Markovian limit, these two unravelings correspond to heterodyne and homodyne detection, respectively. Although we use quantum-measurement theory to define these unravelings, we conclude that the stochastic evolution of the system state is not a true quantum trajectory, as the identity of the state through time is a fiction

Intermittency for stochastic partial differential equations driven by strongly inhomogeneous space-time white noises

Science.gov (United States)

Xie, Bin

2018-01-01

In this paper, the main topic is to investigate the intermittent property of the one-dimensional stochastic heat equation driven by an inhomogeneous Brownian sheet, which is a noise deduced from the study of the catalytic super-Brownian motion. Under some proper conditions on the catalytic measure of the inhomogeneous Brownian sheet, we show that the solution is weakly full intermittent based on the estimates of moments of the solution. In particular, it is proved that the second moment of the solution grows at the exponential rate. The novelty is that the catalytic measure relative to the inhomogeneous noise is not required to be absolutely continuous with respect to the Lebesgue measure on R.

Role of statistical linearization in the solution of nonlinear stochastic equations

International Nuclear Information System (INIS)

Budgor, A.B.

1977-01-01

The solution of a generalized Langevin equation is referred to as a stochastic process. If the external forcing function is Gaussian white noise, the forward Kolmogarov equation yields the transition probability density function. Nonlinear problems must be handled by approximation procedures e.g., perturbation theories, eigenfunction expansions, and nonlinear optimization procedures. After some comments on the first two of these, attention is directed to the third, and the method of statistical linearization is used to demonstrate a relation to the former two. Nonlinear stochastic systems exhibiting sustained or forced oscillations and the centered nonlinear Schroedinger equation in the presence of Gaussian white noise excitation are considered as examples. 5 figures, 2 tables

Arithmetic differential equations on $GL_n$, I: differential cocycles

OpenAIRE

Buium, Alexandru; Dupuy, Taylor

2013-01-01

The theory of differential equations has an arithmetic analogue in which derivatives are replaced by Fermat quotients. One can then ask what is the arithmetic analogue of a linear differential equation. The study of usual linear differential equations is the same as the study of the differential cocycle from $GL_n$ into its Lie algebra given by the logarithmic derivative. However we prove here that there are no such cocycles in the context of arithmetic differential equations. In sequels of t...

Pathwise Strategies for Stochastic Differential Games with an Erratum to “Stochastic Differential Games with Asymmetric Information”

International Nuclear Information System (INIS)

Cardaliaguet, P.; Rainer, C.

2013-01-01

We introduce a new notion of pathwise strategies for stochastic differential games. This allows us to give a correct meaning to some statement asserted in Cardaliaguet and Rainer (Appl. Math. Optim. 59: 1–36, 2009)

Beginning partial differential equations

CERN Document Server

O'Neil, Peter V

2011-01-01

A rigorous, yet accessible, introduction to partial differential equations-updated in a valuable new edition Beginning Partial Differential Equations, Second Edition provides a comprehensive introduction to partial differential equations (PDEs) with a special focus on the significance of characteristics, solutions by Fourier series, integrals and transforms, properties and physical interpretations of solutions, and a transition to the modern function space approach to PDEs. With its breadth of coverage, this new edition continues to present a broad introduction to the field, while also addres

White noise solutions to the stochastic mKdV equation

International Nuclear Information System (INIS)

Zhang Zhongjun; Wei Caimin

2009-01-01

In this paper, we present the white noise solutions of the stochastic mKdV equation via the Hermite transformation and variable-coefficient generalized projected Ricatti equation expansion method. These solutions include white noise solitary wave solutions, white noise soliton-like solutions and white noise trigonometric function solutions.

Exact Finite-Difference Schemes for d-Dimensional Linear Stochastic Systems with Constant Coefficients

Directory of Open Access Journals (Sweden)

Peng Jiang

2013-01-01

Full Text Available The authors attempt to construct the exact finite-difference schemes for linear stochastic differential equations with constant coefficients. The explicit solutions to Itô and Stratonovich linear stochastic differential equations with constant coefficients are adopted with the view of providing exact finite-difference schemes to solve them. In particular, the authors utilize the exact finite-difference schemes of Stratonovich type linear stochastic differential equations to solve the Kubo oscillator that is widely used in physics. Further, the authors prove that the exact finite-difference schemes can preserve the symplectic structure and first integral of the Kubo oscillator. The authors also use numerical examples to prove the validity of the numerical methods proposed in this paper.

Stochastic calculus an introduction through theory and exercises

CERN Document Server

Baldi, Paolo

2017-01-01

This book provides a comprehensive introduction to the theory of stochastic calculus and some of its applications. It is the only textbook on the subject to include more than two hundred exercises with complete solutions. After explaining the basic elements of probability, the author introduces more advanced topics such as Brownian motion, martingales and Markov processes. The core of the book covers stochastic calculus, including stochastic differential equations, the relationship to partial differential equations, numerical methods and simulation, as well as applications of stochastic processes to finance. The final chapter provides detailed solutions to all exercises, in some cases presenting various solution techniques together with a discussion of advantages and drawbacks of the methods used. Stochastic Calculus will be particularly useful to advanced undergraduate and graduate students wishing to acquire a solid understanding of the subject through the theory and exercises. Including full mathematical ...

Global existence and regularity for the 3D stochastic primitive equations of the ocean and atmosphere with multiplicative white noise

Science.gov (United States)

Debussche, A.; Glatt-Holtz, N.; Temam, R.; Ziane, M.

2012-07-01

The primitive equations (PEs) are a basic model in the study of large scale oceanic and atmospheric dynamics. These systems form the analytical core of the most advanced general circulation models. For this reason and due to their challenging nonlinear and anisotropic structure, the PEs have recently received considerable attention from the mathematical community. On the other hand, in view of the complex multi-scale nature of the earth's climate system, many uncertainties appear that should be accounted for in the basic dynamical models of atmospheric and oceanic processes. In the climate community stochastic methods have come into extensive use in this connection. For this reason there has appeared a need to further develop the foundations of nonlinear stochastic partial differential equations in connection with the PEs and more generally. In this work we study a stochastic version of the PEs. We establish the global existence and uniqueness of strong, pathwise solutions for these equations in dimension 3 for the case of a nonlinear multiplicative noise. The proof makes use of anisotropic estimates, L^{p}_{t}L^{q}_{x} estimates on the pressure and stopping time arguments.

Global existence and regularity for the 3D stochastic primitive equations of the ocean and atmosphere with multiplicative white noise

International Nuclear Information System (INIS)

Debussche, A; Glatt-Holtz, N; Temam, R; Ziane, M

2012-01-01

The primitive equations (PEs) are a basic model in the study of large scale oceanic and atmospheric dynamics. These systems form the analytical core of the most advanced general circulation models. For this reason and due to their challenging nonlinear and anisotropic structure, the PEs have recently received considerable attention from the mathematical community. On the other hand, in view of the complex multi-scale nature of the earth's climate system, many uncertainties appear that should be accounted for in the basic dynamical models of atmospheric and oceanic processes. In the climate community stochastic methods have come into extensive use in this connection. For this reason there has appeared a need to further develop the foundations of nonlinear stochastic partial differential equations in connection with the PEs and more generally. In this work we study a stochastic version of the PEs. We establish the global existence and uniqueness of strong, pathwise solutions for these equations in dimension 3 for the case of a nonlinear multiplicative noise. The proof makes use of anisotropic estimates, L p t L q x estimates on the pressure and stopping time arguments

Numerical schemes for dynamically orthogonal equations of stochastic fluid and ocean flows

International Nuclear Information System (INIS)

Ueckermann, M.P.; Lermusiaux, P.F.J.; Sapsis, T.P.

2013-01-01

The quantification of uncertainties is critical when systems are nonlinear and have uncertain terms in their governing equations or are constrained by limited knowledge of initial and boundary conditions. Such situations are common in multiscale, intermittent and non-homogeneous fluid and ocean flows. The dynamically orthogonal (DO) field equations provide an adaptive methodology to predict the probability density functions of such flows. The present work derives efficient computational schemes for the DO methodology applied to unsteady stochastic Navier–Stokes and Boussinesq equations, and illustrates and studies the numerical aspects of these schemes. Semi-implicit projection methods are developed for the mean and for the DO modes, and time-marching schemes of first to fourth order are used for the stochastic coefficients. Conservative second-order finite-volumes are employed in physical space with new advection schemes based on total variation diminishing methods. Other results include: (i) the definition of pseudo-stochastic pressures to obtain a number of pressure equations that is linear in the subspace size instead of quadratic; (ii) symmetric advection schemes for the stochastic velocities; (iii) the use of generalized inversion to deal with singular subspace covariances or deterministic modes; and (iv) schemes to maintain orthonormal modes at the numerical level. To verify our implementation and study the properties of our schemes and their variations, a set of stochastic flow benchmarks are defined including asymmetric Dirac and symmetric lock-exchange flows, lid-driven cavity flows, and flows past objects in a confined channel. Different Reynolds number and Grashof number regimes are employed to illustrate robustness. Optimal convergence under both time and space refinements is shown as well as the convergence of the probability density functions with the number of stochastic realizations.

Hyperbolic partial differential equations

CERN Document Server

Witten, Matthew

1986-01-01

Hyperbolic Partial Differential Equations III is a refereed journal issue that explores the applications, theory, and/or applied methods related to hyperbolic partial differential equations, or problems arising out of hyperbolic partial differential equations, in any area of research. This journal issue is interested in all types of articles in terms of review, mini-monograph, standard study, or short communication. Some studies presented in this journal include discretization of ideal fluid dynamics in the Eulerian representation; a Riemann problem in gas dynamics with bifurcation; periodic M

On the interpretations of Langevin stochastic equation in different coordinate systems

International Nuclear Information System (INIS)

Martinez, E.; Lopez-Diaz, L.; Torres, L.; Alejos, O.

2004-01-01

The stochastic Langevin Landau-Lifshitz equation is usually utilized in micromagnetics formalism to account for thermal effects. Commonly, two different interpretations of the stochastic integrals can be made: Ito and Stratonovich. In this work, the Langevin-Landau-Lifshitz (LLL) equation is written in both Cartesian and Spherical coordinates. If Spherical coordinates are employed, the noise is additive, and therefore, Ito and Stratonovich solutions are equal. This is not the case when (LLL) equation is written in Cartesian coordinates. In this case, the Langevin equation must be interpreted in the Stratonovich sense in order to reproduce correct statistical results. Nevertheless, the statistics of the numerical results obtained from Euler-Ito and Euler-Stratonovich schemes are equivalent due to the additional numerical constraint imposed in Cartesian system after each time step, which itself assures that the magnitude of the magnetization is preserved

A stochastic version of the Price equation reveals the interplay of deterministic and stochastic processes in evolution

Directory of Open Access Journals (Sweden)

Rice Sean H

2008-09-01

Full Text Available Abstract Background Evolution involves both deterministic and random processes, both of which are known to contribute to directional evolutionary change. A number of studies have shown that when fitness is treated as a random variable, meaning that each individual has a distribution of possible fitness values, then both the mean and variance of individual fitness distributions contribute to directional evolution. Unfortunately the most general mathematical description of evolution that we have, the Price equation, is derived under the assumption that both fitness and offspring phenotype are fixed values that are known exactly. The Price equation is thus poorly equipped to study an important class of evolutionary processes. Results I present a general equation for directional evolutionary change that incorporates both deterministic and stochastic processes and applies to any evolving system. This is essentially a stochastic version of the Price equation, but it is derived independently and contains terms with no analog in Price's formulation. This equation shows that the effects of selection are actually amplified by random variation in fitness. It also generalizes the known tendency of populations to be pulled towards phenotypes with minimum variance in fitness, and shows that this is matched by a tendency to be pulled towards phenotypes with maximum positive asymmetry in fitness. This equation also contains a term, having no analog in the Price equation, that captures cases in which the fitness of parents has a direct effect on the phenotype of their offspring. Conclusion Directional evolution is influenced by the entire distribution of individual fitness, not just the mean and variance. Though all moments of individuals' fitness distributions contribute to evolutionary change, the ways that they do so follow some general rules. These rules are invisible to the Price equation because it describes evolution retrospectively. An equally general

On matrix fractional differential equations

Directory of Open Access Journals (Sweden)

Adem Kılıçman

2017-01-01

Full Text Available The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objective of this article is to discuss the Laplace transform method based on operational matrices of fractional derivatives for solving several kinds of linear fractional differential equations. Moreover, we present the operational matrices of fractional derivatives with Laplace transform in many applications of various engineering systems as control system. We present the analytical technique for solving fractional-order, multi-term fractional differential equation. In other words, we propose an efficient algorithm for solving fractional matrix equation.

Stochastic integration in Banach spaces theory and applications

CERN Document Server

Mandrekar, Vidyadhar

2015-01-01

Considering Poisson random measures as the driving sources for stochastic (partial) differential equations allows us to incorporate jumps and to model sudden, unexpected phenomena. By using such equations the present book introduces a new method for modeling the states of complex systems perturbed by random sources over time, such as interest rates in financial markets or temperature distributions in a specific region. It studies properties of the solutions of the stochastic equations, observing the long-term behavior and the sensitivity of the solutions to changes in the initial data. The authors consider an integration theory of measurable and adapted processes in appropriate Banach spaces as well as the non-Gaussian case, whereas most of the literature only focuses on predictable settings in Hilbert spaces. The book is intended for graduate students and researchers in stochastic (partial) differential equations, mathematical finance and non-linear filtering and assumes a knowledge of the required integrati...

Differential equations extended to superspace

International Nuclear Information System (INIS)

Torres, J.; Rosu, H.C.

2003-01-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

Differential equations extended to superspace

Energy Technology Data Exchange (ETDEWEB)

Torres, J. [Instituto de Fisica, Universidad de Guanajuato, A.P. E-143, Leon, Guanajuato (Mexico); Rosu, H.C. [Instituto Potosino de Investigacion Cientifica y Tecnologica, A.P. 3-74, Tangamanga, San Luis Potosi (Mexico)

2003-07-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

Beginning partial differential equations

CERN Document Server

O'Neil, Peter V

2014-01-01

A broad introduction to PDEs with an emphasis on specialized topics and applications occurring in a variety of fields Featuring a thoroughly revised presentation of topics, Beginning Partial Differential Equations, Third Edition provides a challenging, yet accessible,combination of techniques, applications, and introductory theory on the subjectof partial differential equations. The new edition offers nonstandard coverageon material including Burger's equation, the telegraph equation, damped wavemotion, and the use of characteristics to solve nonhomogeneous problems. The Third Edition is or

Differential equations problem solver

CERN Document Server

Arterburn, David R

2012-01-01

REA's Problem Solvers is a series of useful, practical, and informative study guides. Each title in the series is complete step-by-step solution guide. The Differential Equations Problem Solver enables students to solve difficult problems by showing them step-by-step solutions to Differential Equations problems. The Problem Solvers cover material ranging from the elementary to the advanced and make excellent review books and textbook companions. They're perfect for undergraduate and graduate studies.The Differential Equations Problem Solver is the perfect resource for any class, any exam, and

Sparse learning of stochastic dynamical equations

Science.gov (United States)

Boninsegna, Lorenzo; Nüske, Feliks; Clementi, Cecilia

2018-06-01

With the rapid increase of available data for complex systems, there is great interest in the extraction of physically relevant information from massive datasets. Recently, a framework called Sparse Identification of Nonlinear Dynamics (SINDy) has been introduced to identify the governing equations of dynamical systems from simulation data. In this study, we extend SINDy to stochastic dynamical systems which are frequently used to model biophysical processes. We prove the asymptotic correctness of stochastic SINDy in the infinite data limit, both in the original and projected variables. We discuss algorithms to solve the sparse regression problem arising from the practical implementation of SINDy and show that cross validation is an essential tool to determine the right level of sparsity. We demonstrate the proposed methodology on two test systems, namely, the diffusion in a one-dimensional potential and the projected dynamics of a two-dimensional diffusion process.

Stochastic integration by parts and functional Itô calculus

CERN Document Server

Vives, Josep

2016-01-01

This volume contains lecture notes from the courses given by Vlad Bally and Rama Cont at the Barcelona Summer School on Stochastic Analysis (July 2012). The notes of the course by Vlad Bally, co-authored with Lucia Caramellino, develop integration by parts formulas in an abstract setting, extending Malliavin's work on abstract Wiener spaces. The results are applied to prove absolute continuity and regularity results of the density for a broad class of random processes. Rama Cont's notes provide an introduction to the Functional Itô Calculus, a non-anticipative functional calculus that extends the classical Itô calculus to path-dependent functionals of stochastic processes. This calculus leads to a new class of path-dependent partial differential equations, termed Functional Kolmogorov Equations, which arise in the study of martingales and forward-backward stochastic differential equations. This book will appeal to both young and senior researchers in probability and stochastic processes, as well as to pract...

Hyperbolic white noise functional solutions of Wick-type stochastic compound KdV-Burgers equations

International Nuclear Information System (INIS)

Han Xiu; Xie Yingchao

2009-01-01

Variable coefficient and Wick-type stochastic compound KdV-Burgers equations are investigated. By using white noise analysis, Hermite transform and the hyperbolic function method, we obtain a number of Wick versions of hyperbolic white noise functional solutions and hyperbolic function solutions for Wick-type stochastic and variable coefficient compound KdV-Burgers equations, respectively.

Pathwise Strategies for Stochastic Differential Games with an Erratum to 'Stochastic Differential Games with Asymmetric Information'

Energy Technology Data Exchange (ETDEWEB)

Cardaliaguet, P., E-mail: cardaliaguet@ceremade.dauphine.fr [Universite Paris-Dauphine, Ceremade (France); Rainer, C., E-mail: Catherine.Rainer@univ-brest.fr [Universite de Bretagne Occidentale (France)

2013-08-01

We introduce a new notion of pathwise strategies for stochastic differential games. This allows us to give a correct meaning to some statement asserted in Cardaliaguet and Rainer (Appl. Math. Optim. 59: 1-36, 2009)

Study of the stochastic point reactor kinetic equation

International Nuclear Information System (INIS)

Gotoh, Yorio

1980-01-01

Diagrammatic technique is used to solve the stochastic point reactor kinetic equation. The method gives exact results which are derived from Fokker-Plank theory. A Green's function dressed with the clouds of noise is defined, which is a transfer function of point reactor with fluctuating reactivity. An integral equation for the correlation function of neutron power is derived using the following assumptions: 1) Green's funntion should be dressed with noise, 2) The ladder type diagrams only contributes to the correlation function. For a white noise and the one delayed neutron group approximation, the norm of the integral equation and the variance to mean-squared ratio are analytically obtained. (author)

4th International Conference on Particle Systems and Partial Differential Equations

CERN Document Server

Soares, Ana

2017-01-01

'This book addresses mathematical problems motivated by various applications in physics, engineering, chemistry and biology. It gathers the lecture notes from the mini-course presented by Jean-Christophe Mourrat on the construction of the various stochastic “basic” terms involved in the formulation of the dynamic Ö4  theory in three space dimensions, as well as selected contributions presented at the fourth meeting on Particle Systems and PDEs, which was held at the University of Minho’s Centre of Mathematics in December 2015. The purpose of the conference was to bring together prominent researchers working in the fields of particle systems and partial differential equations, offering them a forum to present their recent results and discuss their topics of expertise. The meeting was also intended to present to a vast and varied public, including young researchers, the area of interacting particle systems, its underlying motivation, and its relation to partial differential equations.  The book w...

Linear determining equations for differential constraints

International Nuclear Information System (INIS)

Kaptsov, O V

1998-01-01

A construction of differential constraints compatible with partial differential equations is considered. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the classical determining equations used in the search for admissible Lie operators. As applications of this approach equations of an ideal incompressible fluid and non-linear heat equations are discussed

Homogenization of the stochastic Navier–Stokes equation with a stochastic slip boundary condition

KAUST Repository

Bessaih, Hakima

2015-11-02

The two-dimensional Navier–Stokes equation in a perforated domain with a dynamical slip boundary condition is considered. We assume that the dynamic is driven by a stochastic perturbation on the interior of the domain and another stochastic perturbation on the boundaries of the holes. We consider a scaling (ᵋ for the viscosity and 1 for the density) that will lead to a time-dependent limit problem. However, the noncritical scaling (ᵋ, β > 1) is considered in front of the nonlinear term. The homogenized system in the limit is obtained as a Darcy’s law with memory with two permeabilities and an extra term that is due to the stochastic perturbation on the boundary of the holes. The nonhomogeneity on the boundary contains a stochastic part that yields in the limit an additional term in the Darcy’s law. We use the two-scale convergence method after extending the solution with 0 inside the holes to pass to the limit. By Itô stochastic calculus, we get uniform estimates on the solution in appropriate spaces. Due to the stochastic integral, the pressure that appears in the variational formulation does not have enough regularity in time. This fact made us rely only on the variational formulation for the passage to the limit on the solution. We obtain a variational formulation for the limit that is solution of a Stokes system with two pressures. This two-scale limit gives rise to three cell problems, two of them give the permeabilities while the third one gives an extra term in the Darcy’s law due to the stochastic perturbation on the boundary of the holes.

Solving Ordinary Differential Equations

Science.gov (United States)

Krogh, F. T.

1987-01-01

Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Algebraic entropy for differential-delay equations

OpenAIRE

Viallet, Claude M.

2014-01-01

We extend the definition of algebraic entropy to a class of differential-delay equations. The vanishing of the entropy, as a structural property of an equation, signals its integrability. We suggest a simple way to produce differential-delay equations with vanishing entropy from known integrable differential-difference equations.

Stochastic differential equtions with non-lipschitz coefficients:II. Dependence with respect to initial values

OpenAIRE

Fang, Shizan; Zhang, Tusheng

2003-01-01

14 pages; The existence of the unique strong solution for a class of stochastic differential equations with non-Lipschitz coefficients was established recently. In this paper, we shall investigate the dependence with respect to the initial values. We shall prove that the non confluence of solutions holds under our general conditions. To obtain a continuous version, the modulus of continuity of coefficients is assumed to be less than $\\dis |x-y|\\log{1\\over|x-y|}$. In this case, it will give ri...

Application of Stochastic Partial Differential Equations to Reservoir Property Modelling

KAUST Repository

Potsepaev, R.; Farmer, C.L.

2010-01-01

in parametric space. In order to sample in physical space we introduce a stochastic elliptic PDE with tensor coefficients, where the tensor is related to correlation anisotropy and its variation is physical space.

Research on nonlinear stochastic dynamical price model

International Nuclear Information System (INIS)

Li Jiaorui; Xu Wei; Xie Wenxian; Ren Zhengzheng

2008-01-01

In consideration of many uncertain factors existing in economic system, nonlinear stochastic dynamical price model which is subjected to Gaussian white noise excitation is proposed based on deterministic model. One-dimensional averaged Ito stochastic differential equation for the model is derived by using the stochastic averaging method, and applied to investigate the stability of the trivial solution and the first-passage failure of the stochastic price model. The stochastic price model and the methods presented in this paper are verified by numerical studies

Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves

DEFF Research Database (Denmark)

Eldeberky, Y.; Madsen, Per A.

1999-01-01

and stochastic formulations are solved numerically for the case of cross shore motion of unidirectional waves and the results are verified against laboratory data for wave propagation over submerged bars and over a plane slope. Outside the surf zone the two model predictions are generally in good agreement......This paper presents a new and more accurate set of deterministic evolution equations for the propagation of fully dispersive, weakly nonlinear, irregular, multidirectional waves. The equations are derived directly from the Laplace equation with leading order nonlinearity in the surface boundary...... is significantly underestimated for larger wave numbers. In the present work we correct this inconsistency. In addition to the improved deterministic formulation, we present improved stochastic evolution equations in terms of the energy spectrum and the bispectrum for multidirectional waves. The deterministic...

Pesin’s entropy formula for stochastic flows of diffeomorphisms

Institute of Scientific and Technical Information of China (English)

刘培东

1996-01-01

Pesin’s entropy formula relating entropy and Lyapunov exponents within the context of random dynamical systems generated by (discrete or continuous) stochastic flows of diffeomorphisms (including solution flows of stochastic differential equations on manifolds) is proved.

Iterative Splitting Methods for Differential Equations

CERN Document Server

Geiser, Juergen

2011-01-01

Iterative Splitting Methods for Differential Equations explains how to solve evolution equations via novel iterative-based splitting methods that efficiently use computational and memory resources. It focuses on systems of parabolic and hyperbolic equations, including convection-diffusion-reaction equations, heat equations, and wave equations. In the theoretical part of the book, the author discusses the main theorems and results of the stability and consistency analysis for ordinary differential equations. He then presents extensions of the iterative splitting methods to partial differential

Semi-groups of operators and some of their applications to partial differential equations

International Nuclear Information System (INIS)

Kisynski, J.

1976-01-01

Basic notions and theorems of the theory of one-parameter semi-groups of linear operators are given, illustrated by some examples concerned with linear partial differential operators. For brevity, some important and widely developed parts of the semi-group theory such as the general theory of holomorphic semi-groups or the theory of temporally inhomogeneous evolution equations are omitted. This omission includes also the very important application of semi-groups to investigating stochastic processes. (author)

Stochastic reliability analysis using Fokker Planck equations

International Nuclear Information System (INIS)

Hari Prasad, M.; Rami Reddy, G.; Srividya, A.; Verma, A.K.

2011-01-01

The Fokker-Planck equation describes the time evolution of the probability density function of the velocity of a particle, and can be generalized to other observables as well. It is also known as the Kolmogorov forward equation (diffusion). Hence, for any process, which evolves with time, the probability density function as a function of time can be represented with Fokker-Planck equation. In stochastic reliability analysis one is more interested in finding out the reliability or failure probability of the components or structures as a function of time rather than instantaneous failure probabilities. In this analysis the variables are represented with random processes instead of random variables. A random processes can be either stationary or non stationary. If the random process is stationary then the failure probability doesn't change with time where as in the case of non stationary processes the failure probability changes with time. In the present paper Fokker Planck equations have been used to find out the probability density function of the non stationary random processes. In this paper a flow chart has been provided which describes step by step process for carrying out stochastic reliability analysis using Fokker-Planck equations. As a first step one has to identify the failure function as a function of random processes. Then one has to solve the Fokker-Planck equation for each random process. In this paper the Fokker-Planck equation has been solved by using Finite difference method. As a result one gets the probability density values of the random process in the sample space as well as time space. Later at each time step appropriate probability distribution has to be identified based on the available probability density values. For checking the better fitness of the data Kolmogorov-Smirnov Goodness of fit test has been performed. In this way one can find out the distribution of the random process at each time step. Once one has the probability distribution

New derivation of quantum equations from classical stochastic arguments

OpenAIRE

Bergeron, H.

2003-01-01

In a previous article [H. Bergeron, J. Math. Phys. 42, 3983 (2001)], we presented a method to obtain a continuous transition from classical to quantum mechanics starting from the usual phase space formulation of classical mechanics. This procedure was based on a Koopman-von Neumann approach where classical equations are reformulated into a quantumlike form. In this article, we develop a different derivation of quantum equations, based on purely classical stochastic arguments, taking some elem...

Differential equation analysis in biomedical science and engineering ordinary differential equation applications with R

CERN Document Server

Schiesser, William E

2014-01-01

Features a solid foundation of mathematical and computational tools to formulate and solve real-world ODE problems across various fields With a step-by-step approach to solving ordinary differential equations (ODEs), Differential Equation Analysis in Biomedical Science and Engineering: Ordinary Differential Equation Applications with R successfully applies computational techniques for solving real-worldODE problems that are found in a variety of fields, including chemistry, physics, biology,and physiology. The book provides readers with the necessary knowledge to reproduce andextend the comp