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...

Statistical Methods for Stochastic Differential Equations

CERN Document Server

Kessler, Mathieu; Sorensen, Michael

2012-01-01

The seventh volume in the SemStat series, Statistical Methods for Stochastic Differential Equations presents current research trends and recent developments in statistical methods for stochastic differential equations. Written to be accessible to both new students and seasoned researchers, each self-contained chapter starts with introductions to the topic at hand and builds gradually towards discussing recent research. The book covers Wiener-driven equations as well as stochastic differential equations with jumps, including continuous-time ARMA processes and COGARCH processes. It presents a sp

Stochastic differential equation model to Prendiville processes

International Nuclear Information System (INIS)

Granita; Bahar, Arifah

2015-01-01

The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution

Stochastic differential equation model to Prendiville processes

Energy Technology Data Exchange (ETDEWEB)

Granita, E-mail: granitafc@gmail.com [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); Bahar, Arifah [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); UTM Center for Industrial & Applied Mathematics (UTM-CIAM) (Malaysia)

2015-10-22

The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution.

New evolution equations for the joint response-excitation probability density function of stochastic solutions to first-order nonlinear PDEs

Science.gov (United States)

Venturi, D.; Karniadakis, G. E.

2012-08-01

By using functional integral methods we determine new evolution equations satisfied by the joint response-excitation probability density function (PDF) associated with the stochastic solution to first-order nonlinear partial differential equations (PDEs). The theory is presented for both fully nonlinear and for quasilinear scalar PDEs subject to random boundary conditions, random initial conditions or random forcing terms. Particular applications are discussed for the classical linear and nonlinear advection equations and for the advection-reaction equation. By using a Fourier-Galerkin spectral method we obtain numerical solutions of the proposed response-excitation PDF equations. These numerical solutions are compared against those obtained by using more conventional statistical approaches such as probabilistic collocation and multi-element probabilistic collocation methods. It is found that the response-excitation approach yields accurate predictions of the statistical properties of the system. In addition, it allows to directly ascertain the tails of probabilistic distributions, thus facilitating the assessment of rare events and associated risks. The computational cost of the response-excitation method is order magnitudes smaller than the one of more conventional statistical approaches if the PDE is subject to high-dimensional random boundary or initial conditions. The question of high-dimensionality for evolution equations involving multidimensional joint response-excitation PDFs is also addressed.

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

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 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.

Evolution and mass extinctions as lognormal stochastic processes

Science.gov (United States)

Maccone, Claudio

2014-10-01

In a series of recent papers and in a book, this author put forward a mathematical model capable of embracing the search for extra-terrestrial intelligence (SETI), Darwinian Evolution and Human History into a single, unified statistical picture, concisely called Evo-SETI. The relevant mathematical tools are: (1) Geometric Brownian motion (GBM), the stochastic process representing evolution as the stochastic increase of the number of species living on Earth over the last 3.5 billion years. This GBM is well known in the mathematics of finances (Black-Sholes models). Its main features are that its probability density function (pdf) is a lognormal pdf, and its mean value is either an increasing or, more rarely, decreasing exponential function of the time. (2) The probability distributions known as b-lognormals, i.e. lognormals starting at a certain positive instant b>0 rather than at the origin. These b-lognormals were then forced by us to have their peak value located on the exponential mean-value curve of the GBM (Peak-Locus theorem). In the framework of Darwinian Evolution, the resulting mathematical construction was shown to be what evolutionary biologists call Cladistics. (3) The (Shannon) entropy of such b-lognormals is then seen to represent the `degree of progress' reached by each living organism or by each big set of living organisms, like historic human civilizations. Having understood this fact, human history may then be cast into the language of b-lognormals that are more and more organized in time (i.e. having smaller and smaller entropy, or smaller and smaller `chaos'), and have their peaks on the increasing GBM exponential. This exponential is thus the `trend of progress' in human history. (4) All these results also match with SETI in that the statistical Drake equation (generalization of the ordinary Drake equation to encompass statistics) leads just to the lognormal distribution as the probability distribution for the number of extra

Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation

Directory of Open Access Journals (Sweden)

Hamidreza Rezazadeh

2014-05-01

Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.

Stochastic evolutions and hadronization of highly excited hadronic matter

International Nuclear Information System (INIS)

Carruthers, P.

1984-01-01

Stochastic ingredients of high energy hadronic collisions are analyzed, with emphasis on multiplicity distributions. The conceptual simplicity of the k-cell negative binomial distribution is related to the evolution of probability distributions via the Fokker-Planck and related equations. The connection to underlying field theory ideas is sketched. 17 references

Exact many-body dynamics with stochastic one-body density matrix evolution

International Nuclear Information System (INIS)

Lacroix, D.

2004-05-01

In this article, we discuss some properties of the exact treatment of the many-body problem with stochastic Schroedinger equation (SSE). Starting from the SSE theory, an equivalent reformulation is proposed in terms of quantum jumps in the density matrix space. The technical details of the derivation a stochastic version of the Liouville von Neumann equation are given. It is shown that the exact Many-Body problem could be replaced by an ensemble of one-body density evolution, where each density matrix evolves according to its own mean-field augmented by a one-body noise. (author)

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.)

Lyapunov functionals and stability of stochastic functional differential equations

CERN Document Server

Shaikhet, Leonid

2013-01-01

Stability conditions for functional differential equations can be obtained using Lyapunov functionals. Lyapunov Functionals and Stability of Stochastic Functional Differential Equations describes the general method of construction of Lyapunov functionals to investigate the stability of differential equations with delays. This work continues and complements the author’s previous book Lyapunov Functionals and Stability of Stochastic Difference Equations, where this method is described for discrete- and continuous-time difference equations. The text begins with a description of the peculiarities of deterministic and stochastic functional differential equations. There follow basic definitions for stability theory of stochastic hereditary systems, and a formal procedure of Lyapunov functionals construction is presented. Stability investigation is conducted for stochastic linear and nonlinear differential equations with constant and distributed delays. The proposed method is used for stability investigation of di...

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.

Quantum dynamical time evolutions as stochastic flows on phase space

International Nuclear Information System (INIS)

Combe, P.; Rodriguez, R.; Guerra, F.; Sirigue, M.; Sirigue-Collin, M.

1984-01-01

We are mainly interested in describing the time development of the Wigner functions by means of stochastic processes. In the second section we recall the main properties of the Wigner functions as well as those of their Fourier transform. In the next one we derive the evolution equation of these functions for a class of Hamiltonians and we give a probabilistic expression for the solution of these equations by means of a stochastic flow in phase space which reminds of the classical flows. In the last section we remark that the previously defined flow can be extended to the bounded continuous functions on phase space and that this flow conserves the cone generated by the Wigner functions. (orig./HSI)

A microscopic derivation of stochastic differential equations

International Nuclear Information System (INIS)

Arimitsu, Toshihico

1996-01-01

With the help of the formulation of Non-Equilibrium Thermo Field Dynamics, a unified canonical operator formalism is constructed for the quantum stochastic differential equations. In the course of its construction, it is found that there are at least two formulations, i.e. one is non-hermitian and the other is hermitian. Having settled which framework should be satisfied by the quantum stochastic differential equations, a microscopic derivation is performed for these stochastic differential equations by extending the projector methods. This investigation may open a new field for quantum systems in order to understand the deeper meaning of dissipation

On the fundamental equation of nonequilibrium statistical physics—Nonequilibrium entropy evolution equation and the formula for entropy production rate

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

In this paper the author presents an overview on his own research works. More than ten years ago, we proposed a new fundamental equation of nonequilibrium statistical physics in place of the present Liouville equation. That is the stochastic velocity type’s Langevin equation in 6N dimensional phase space or its equivalent Liouville diffusion equation. This equation is time-reversed asymmetrical. It shows that the form of motion of particles in statistical thermodynamic systems has the drift-diffusion duality, and the law of motion of statistical thermodynamics is expressed by a superposition of both the law of dynamics and the stochastic velocity and possesses both determinism and probability. Hence it is different from the law of motion of particles in dynamical systems. The stochastic diffusion motion of the particles is the microscopic origin of macroscopic irreversibility. Starting from this fundamental equation the BBGKY diffusion equation hierarchy, the Boltzmann collision diffusion equation, the hydrodynamic equations such as the mass drift-diffusion equation, the Navier-Stokes equation and the thermal conductivity equation have been derived and presented here. What is more important, we first constructed a nonlinear evolution equation of nonequilibrium entropy density in 6N, 6 and 3 dimensional phase space, predicted the existence of entropy diffusion. This entropy evolution equation plays a leading role in nonequilibrium entropy theory, it reveals that the time rate of change of nonequilibrium entropy density originates together from its drift, diffusion and production in space. From this evolution equation, we presented a formula for entropy production rate (i.e. the law of entropy increase) in 6N and 6 dimensional phase space, proved that internal attractive force in nonequilibrium system can result in entropy decrease while internal repulsive force leads to another entropy increase, and derived a common expression for this entropy decrease rate or

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...

Stochastic solution of population balance equations for reactor networks

International Nuclear Information System (INIS)

Menz, William J.; Akroyd, Jethro; Kraft, Markus

2014-01-01

This work presents a sequential modular approach to solve a generic network of reactors with a population balance model using a stochastic numerical method. Full-coupling to the gas-phase is achieved through operator-splitting. The convergence of the stochastic particle algorithm in test networks is evaluated as a function of network size, recycle fraction and numerical parameters. These test cases are used to identify methods through which systematic and statistical error may be reduced, including by use of stochastic weighted algorithms. The optimal algorithm was subsequently used to solve a one-dimensional example of silicon nanoparticle synthesis using a multivariate particle model. This example demonstrated the power of stochastic methods in resolving particle structure by investigating the transient and spatial evolution of primary polydispersity, degree of sintering and TEM-style images. Highlights: •An algorithm is presented to solve reactor networks with a population balance model. •A stochastic method is used to solve the population balance equations. •The convergence and efficiency of the reported algorithms are evaluated. •The algorithm is applied to simulate silicon nanoparticle synthesis in a 1D reactor. •Particle structure is reported as a function of reactor length and time

Stochastic Differential Equations and Kondratiev Spaces

Energy Technology Data Exchange (ETDEWEB)

Vaage, G.

1995-05-01

The purpose of this mathematical thesis was to improve the understanding of physical processes such as fluid flow in porous media. An example is oil flowing in a reservoir. In the first of five included papers, Hilbert space methods for elliptic boundary value problems are used to prove the existence and uniqueness of a large family of elliptic differential equations with additive noise without using the Hermite transform. The ideas are then extended to the multidimensional case and used to prove existence and uniqueness of solution of the Stokes equations with additive noise. The second paper uses functional analytic methods for partial differential equations and presents a general framework for proving existence and uniqueness of solutions to stochastic partial differential equations with multiplicative noise, for a large family of noises. The methods are applied to equations of elliptic, parabolic as well as hyperbolic type. The framework presented can be extended to the multidimensional case. The third paper shows how the ideas from the second paper can be extended to study the moving boundary value problem associated with the stochastic pressure equation. The fourth paper discusses a set of stochastic differential equations. The fifth paper studies the relationship between the two families of Kondratiev spaces used in the thesis. 102 refs.

Bipartite Fuzzy Stochastic Differential Equations with Global Lipschitz Condition

Directory of Open Access Journals (Sweden)

Marek T. Malinowski

2016-01-01

Full Text Available We introduce and analyze a new type of fuzzy stochastic differential equations. We consider equations with drift and diffusion terms occurring at both sides of equations. Therefore we call them the bipartite fuzzy stochastic differential equations. Under the Lipschitz and boundedness conditions imposed on drifts and diffusions coefficients we prove existence of a unique solution. Then, insensitivity of the solution under small changes of data of equation is examined. Finally, we mention that all results can be repeated for solutions to bipartite set-valued stochastic differential equations.

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

On Volatility Induced Stationarity for Stochastic Differential Equations

DEFF Research Database (Denmark)

Albin, J.M.P.; Astrup Jensen, Bjarne; Muszta, Anders

2006-01-01

This article deals with stochastic differential equations with volatility induced stationarity. We study of theoretical properties of such equations, as well as numerical aspects, together with a detailed study of three examples.......This article deals with stochastic differential equations with volatility induced stationarity. We study of theoretical properties of such equations, as well as numerical aspects, together with a detailed study of three examples....

Stochastic modeling of mode interactions via linear parabolized stability equations

Science.gov (United States)

Ran, Wei; Zare, Armin; Hack, M. J. Philipp; Jovanovic, Mihailo

2017-11-01

Low-complexity approximations of the Navier-Stokes equations have been widely used in the analysis of wall-bounded shear flows. In particular, the parabolized stability equations (PSE) and Floquet theory have been employed to capture the evolution of primary and secondary instabilities in spatially-evolving flows. We augment linear PSE with Floquet analysis to formally treat modal interactions and the evolution of secondary instabilities in the transitional boundary layer via a linear progression. To this end, we leverage Floquet theory by incorporating the primary instability into the base flow and accounting for different harmonics in the flow state. A stochastic forcing is introduced into the resulting linear dynamics to model the effect of nonlinear interactions on the evolution of modes. We examine the H-type transition scenario to demonstrate how our approach can be used to model nonlinear effects and capture the growth of the fundamental and subharmonic modes observed in direct numerical simulations and experiments.

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.

Diffusion equations and the time evolution of foreign exchange rates

Energy Technology Data Exchange (ETDEWEB)

Figueiredo, Annibal; Castro, Marcio T. de [Institute of Physics, Universidade de Brasília, Brasília DF 70910-900 (Brazil); Fonseca, Regina C.B. da [Department of Mathematics, Instituto Federal de Goiás, Goiânia GO 74055-110 (Brazil); Gleria, Iram, E-mail: iram@fis.ufal.br [Institute of Physics, Federal University of Alagoas, Brazil, Maceió AL 57072-900 (Brazil)

2013-10-01

We investigate which type of diffusion equation is most appropriate to describe the time evolution of foreign exchange rates. We modify the geometric diffusion model assuming a non-exponential time evolution and the stochastic term is the sum of a Wiener noise and a jump process. We find the resulting diffusion equation to obey the Kramers–Moyal equation. Analytical solutions are obtained using the characteristic function formalism and compared with empirical data. The analysis focus on the first four central moments considering the returns of foreign exchange rate. It is shown that the proposed model offers a good improvement over the classical geometric diffusion model.

Diffusion equations and the time evolution of foreign exchange rates

Science.gov (United States)

Figueiredo, Annibal; de Castro, Marcio T.; da Fonseca, Regina C. B.; Gleria, Iram

2013-10-01

We investigate which type of diffusion equation is most appropriate to describe the time evolution of foreign exchange rates. We modify the geometric diffusion model assuming a non-exponential time evolution and the stochastic term is the sum of a Wiener noise and a jump process. We find the resulting diffusion equation to obey the Kramers-Moyal equation. Analytical solutions are obtained using the characteristic function formalism and compared with empirical data. The analysis focus on the first four central moments considering the returns of foreign exchange rate. It is shown that the proposed model offers a good improvement over the classical geometric diffusion model.

Diffusion equations and the time evolution of foreign exchange rates

International Nuclear Information System (INIS)

Figueiredo, Annibal; Castro, Marcio T. de; Fonseca, Regina C.B. da; Gleria, Iram

2013-01-01

We investigate which type of diffusion equation is most appropriate to describe the time evolution of foreign exchange rates. We modify the geometric diffusion model assuming a non-exponential time evolution and the stochastic term is the sum of a Wiener noise and a jump process. We find the resulting diffusion equation to obey the Kramers–Moyal equation. Analytical solutions are obtained using the characteristic function formalism and compared with empirical data. The analysis focus on the first four central moments considering the returns of foreign exchange rate. It is shown that the proposed model offers a good improvement over the classical geometric diffusion model.

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)

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

Stochastic partial differential equations an introduction

CERN Document Server

Liu, Wei

2015-01-01

This book provides an introduction to the theory of stochastic partial differential equations (SPDEs) of evolutionary type. SPDEs are one of the main research directions in probability theory with several wide ranging applications. Many types of dynamics with stochastic influence in nature or man-made complex systems can be modelled by such equations. The theory of SPDEs is based both on the theory of deterministic partial differential equations, as well as on modern stochastic analysis. Whilst this volume mainly follows the ‘variational approach’, it also contains a short account on the ‘semigroup (or mild solution) approach’. In particular, the volume contains a complete presentation of the main existence and uniqueness results in the case of locally monotone coefficients. Various types of generalized coercivity conditions are shown to guarantee non-explosion, but also a systematic approach to treat SPDEs with explosion in finite time is developed. It is, so far, the only book where the latter and t...

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.

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

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.

Symmetries of stochastic differential equations: A geometric approach

Energy Technology Data Exchange (ETDEWEB)

De Vecchi, Francesco C., E-mail: francesco.devecchi@unimi.it; Ugolini, Stefania, E-mail: stefania.ugolini@unimi.it [Dipartimento di Matematica, Università degli Studi di Milano, via Saldini 50, Milano (Italy); Morando, Paola, E-mail: paola.morando@unimi.it [DISAA, Università degli Studi di Milano, via Celoria 2, Milano (Italy)

2016-06-15

A new notion of stochastic transformation is proposed and applied to the study of both weak and strong symmetries of stochastic differential equations (SDEs). The correspondence between an algebra of weak symmetries for a given SDE and an algebra of strong symmetries for a modified SDE is proved under suitable regularity assumptions. This general approach is applied to a stochastic version of a two dimensional symmetric ordinary differential equation and to the case of two dimensional Brownian motion.

An introduction to stochastic differential equations

CERN Document Server

Evans, Lawrence C

2014-01-01

These notes provide a concise introduction to stochastic differential equations and their application to the study of financial markets and as a basis for modeling diverse physical phenomena. They are accessible to non-specialists and make a valuable addition to the collection of texts on the topic. -Srinivasa Varadhan, New York University This is a handy and very useful text for studying stochastic differential equations. There is enough mathematical detail so that the reader can benefit from this introduction with only a basic background in mathematical analysis and probability. -George Papa

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

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 ...

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.

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...

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 equations from linear to fully nonlinear theory

CERN Document Server

Zhang, Jianfeng

2017-01-01

This book provides a systematic and accessible approach to stochastic differential equations, backward stochastic differential equations, and their connection with partial differential equations, as well as the recent development of the fully nonlinear theory, including nonlinear expectation, second order backward stochastic differential equations, and path dependent partial differential equations. Their main applications and numerical algorithms, as well as many exercises, are included. The book focuses on ideas and clarity, with most results having been solved from scratch and most theories being motivated from applications. It can be considered a starting point for junior researchers in the field, and can serve as a textbook for a two-semester graduate course in probability theory and stochastic analysis. It is also accessible for graduate students majoring in financial engineering.

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)

Analysis of stability for stochastic delay integro-differential equations.

Science.gov (United States)

Zhang, Yu; Li, Longsuo

2018-01-01

In this paper, we concern stability of numerical methods applied to stochastic delay integro-differential equations. For linear stochastic delay integro-differential equations, it is shown that the mean-square stability is derived by the split-step backward Euler method without any restriction on step-size, while the Euler-Maruyama method could reproduce the mean-square stability under a step-size constraint. We also confirm the mean-square stability of the split-step backward Euler method for nonlinear stochastic delay integro-differential equations. The numerical experiments further verify the theoretical results.

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

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

Asymptotic analysis for functional stochastic differential equations

CERN Document Server

Bao, Jianhai; Yuan, Chenggui

2016-01-01

This brief treats dynamical systems that involve delays and random disturbances. The study is motivated by a wide variety of systems in real life in which random noise has to be taken into consideration and the effect of delays cannot be ignored. Concentrating on such systems that are described by functional stochastic differential equations, this work focuses on the study of large time behavior, in particular, ergodicity. This brief is written for probabilists, applied mathematicians, engineers, and scientists who need to use delay systems and functional stochastic differential equations in their work. Selected topics from the brief can also be used in a graduate level topics course in probability and stochastic processes.

Numerical resolution of N-dimensional Fokker-Plank stochastic equations

International Nuclear Information System (INIS)

Garcia-Olivares, A.; Muoz, A.

1992-01-01

This document describes the use of a library of programs able to solve stochastic Fokker-Planck equations in a N-dimensional space. the input data are essentially: (i) the initial distribution of the stochastic variable, (ii) the drift and fluctuation coefficients as a function of the state (which can be obtained from the transition probabilities between neighboring states) and (iii) some parameters controlling the run. A last version of the library accepts sources and sinks defined in the states space. The output is the temporal evolution of the probability distribution in the space defined by a N-dimensional grid. Some applications and readings in Synergetics, Self-Organization, transport phenomena, Ecology and other fields are suggested. If the probability distribution is interpreted as a distribution of particles then the codes can be used to solve the N-dimensional problem of advection-diffusion. (author) 21 fig. 16 ref

Numerical Resolution of N-dimensional Fokker-Planck stochastic equations

International Nuclear Information System (INIS)

Garcia-Olivares, R. A.; Munoz Roldan, A.

1992-01-01

This document describes the use of a library of programs able to solve stochastic Fokker-Planck equations in a N-dimensional space. The input data are essentially: (i) the initial distribution of the stochastic variable, (ii) the drift and fluctuation coefficients as a function of the state (which can be obtained from the transition probabilities between neighboring states) and (iii) some parameters controlling the run. A last version of the library accepts sources and sinks defined in the states space. The output is the temporal evolution of the probability distribution in the space defined by a N-dimensional grid. Some applications and readings in Synergetic, Self-Organization, transport phenomena, Ecology and other fields are suggested. If the probability distribution is interpreted as a distribution of particles then the codes can be used to solve the N-dimensional problem of advection-diffusion. (Author) 16 refs

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.

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.

Neutral Backward Stochastic Functional Differential Equations and Their Application

OpenAIRE

Wei, Wenning

2013-01-01

In this paper we are concerned with a new type of backward equations with anticipation which we call neutral backward stochastic functional differential equations. We obtain the existence and uniqueness and prove a comparison theorem. As an application, we discuss the optimal control of neutral stochastic functional differential equations, establish a Pontryagin maximum principle, and give an explicit optimal value for the linear optimal control.

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)

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

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.)

Expressing stochastic unravellings using random evolution operators

International Nuclear Information System (INIS)

Salgado, D; Sanchez-Gomez, J L

2002-01-01

We prove how the form of the most general invariant stochastic unravelling for Markovian (recently given in the literature by Wiseman and Diosi) and non-Markovian but Lindblad-type open quantum systems can be attained by imposing a single mathematical condition upon the random evolution operator of the system, namely a.s. trace preservation (a.s. stands for almost surely). The use of random operators ensures the complete positivity of the density operator evolution and characterizes the linear/non-linear character of the evolution in a straightforward way. It is also shown how three quantum stochastic evolution models - continuous spontaneous localization, quantum state diffusion and quantum mechanics with universal position localization - appear as concrete choices for the noise term of the evolution random operators are assumed. We finally conjecture how these operators may in the future be used in two different directions: both to connect quantum stochastic evolution models with random properties of space-time and to handle noisy quantum logical gates

Operations involving momentum variables in non-Hamiltonian evolution equations

International Nuclear Information System (INIS)

Benatti, F.; Ghirardi, G.C.; Rimini, A.; Weber, T.

1988-02-01

Non-Hamiltonian evolution equations have been recently considered for the description of various physical processes. Among this type of equations the class which has been more extensively studied is the one usually referred to as Quantum Dynamical Semigroup equations (QDS). In particular an equation of the QDS type has been considered as the basis for a model, called Quantum Mechanics with Spontaneous Localization (QMSL), which has been shown to exhibit some very interesting features allowing to overcome most of the conceptual difficulties of standard quantum theory, QMSL assumes a modification of the pure Schroedinger evolution by assuming the occurrence, at random times, of stochastic processes for the wave function corresponding formally to approximate position measurements. In this paper, we investigate the consequences of modifying and/or enlarging the class of the considered stochastic processes, by considering the spontaeous occurrence of approximate momentum and of simultaneous position and momentum measurements. It is shown that the considered changes in the elementary processes have unacceptable consequences. In particular they either lead to drastic modifications in the dynamics of microsystems or are completely useless from the point of view of the conceptual advantages that one was trying to get from QMSL. The present work supports therefore the idea that QMSL, as originally formulated, can be taken as the basic scheme for the generalizations which are still necessary in order to make it appropriate for the description of systems of identical particles and to meet relativistic requirements. (author). 14 refs

Operations involving momentum variables in non-Hamiltonian evolution equation

International Nuclear Information System (INIS)

Benatti, F.; Ghirardi, G.C.; Weber, T.; Rimini, A.

1988-01-01

Non-Hamiltonian evolution equations have been recently considered for the description of various physical processes. Among these types of equations the class which has been more extensively studied is the one usually referred to as quantum-dynamical semi-group equations (QDS). In particular an equation of the QDS type has been considered as the basis for a model, called quantum mechanics with spontaneous localization (QMSL), which has been shown to exhibit some very interesting features allowing us to overcome most of the conceptual difficulties of standard quantum theory. QMSL assumes a modification of the pure Schroedinger evolution by assuming the occurrence, at random times, of stochastic processes for the wave function corresponding formally to approximate position measurements. In this paper the consequences of modifying and/or enlarging the class of the considered stochastic processes, by considering the spontaneous occurrence of approximate momentum and of simultaneous position and momentum measurements, are investigated. It is shown that the considered changes in the elementary processes have unacceptable consequences. In particular they either lead to drastic modification in the dynamics of microsystems or are completely useless from the point of view of the conceptual advantages that one was trying to get from QMSL. The present work supports therefore the idea that QMSL, as originally formulated, can be taken as the basic scheme for the generalizations which are still necessary in order to make it appropriate for the description of systems of identical particles and to meet relativistic requirements

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.

Evolution in fluctuating environments: decomposing selection into additive components of the Robertson-Price equation.

Science.gov (United States)

Engen, Steinar; Saether, Bernt-Erik

2014-03-01

We analyze the stochastic components of the Robertson-Price equation for the evolution of quantitative characters that enables decomposition of the selection differential into components due to demographic and environmental stochasticity. We show how these two types of stochasticity affect the evolution of multivariate quantitative characters by defining demographic and environmental variances as components of individual fitness. The exact covariance formula for selection is decomposed into three components, the deterministic mean value, as well as stochastic demographic and environmental components. We show that demographic and environmental stochasticity generate random genetic drift and fluctuating selection, respectively. This provides a common theoretical framework for linking ecological and evolutionary processes. Demographic stochasticity can cause random variation in selection differentials independent of fluctuating selection caused by environmental variation. We use this model of selection to illustrate that the effect on the expected selection differential of random variation in individual fitness is dependent on population size, and that the strength of fluctuating selection is affected by how environmental variation affects the covariance in Malthusian fitness between individuals with different phenotypes. Thus, our approach enables us to partition out the effects of fluctuating selection from the effects of selection due to random variation in individual fitness caused by demographic stochasticity. © 2013 The Author(s). Evolution © 2013 The Society for the Study of Evolution.

Potential in stochastic differential equations: novel construction

International Nuclear Information System (INIS)

Ao, P

2004-01-01

There is a whole range of emergent phenomena in a complex network such as robustness, adaptiveness, multiple-equilibrium, hysteresis, oscillation and feedback. Those non-equilibrium behaviours can often be described by a set of stochastic differential equations. One persistent important question is the existence of a potential function. Here we demonstrate that a dynamical structure built into stochastic differential equation allows us to construct such a global optimization potential function. We present an explicit construction procedure to obtain the potential and relevant quantities. In the procedure no reference to the Fokker-Planck equation is needed. The availability of the potential suggests that powerful statistical mechanics tools can be used in nonequilibrium situations. (letter to the editor)

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.

A general comparison theorem for backward stochastic differential equations

OpenAIRE

Cohen, Samuel N.; Elliott, Robert J.; Pearce, Charles E. M.

2010-01-01

A useful result when dealing with backward stochastic differential equations is the comparison theorem of Peng (1992). When the equations are not based on Brownian motion, the comparison theorem no longer holds in general. In this paper we present a condition for a comparison theorem to hold for backward stochastic differential equations based on arbitrary martingales. This theorem applies to both vector and scalar situations. Applications to the theory of nonlinear expectat...

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.

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....

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

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.

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.

Stochastic wave-function unravelling of the generalized Lindblad equation using correlated states

International Nuclear Information System (INIS)

Moodley, Mervlyn; Nsio Nzundu, T; Paul, S

2012-01-01

We perform a stochastic wave-function unravelling of the generalized Lindblad master equation using correlated states, a combination of the system state vectors and the environment population. The time-convolutionless projection operator method using correlated projection superoperators is applied to a two-state system, a qubit, that is coupled to an environment consisting of two energy bands which are both populated. These results are compared to the data obtained from Monte Carlo wave-function simulations based on the unravelling of the master equation. We also show a typical quantum trajectory and the average time evolution of the state vector on the Bloch sphere. (paper)

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.

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.

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.

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

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...

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...

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...

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.

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.

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.

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.)

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.

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 ...

Stochastic evolution of cosmological parameters in the early universe

Indian Academy of Sciences (India)

We develop a stochastic formulation of cosmology in the early universe, after considering the scatter in the redshift-apparent magnitude diagram in the early epochs as an observational evidence for the non-deterministic evolution of early universe. We consider the stochastic evolution of density parameter in the early ...

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

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.

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.

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

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.

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...

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)

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

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.

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

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.

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.

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.

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.

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.

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)

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

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....

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...

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....

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...

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

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.

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...

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

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.

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...

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.)

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...

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

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.

Antisymmetrized molecular dynamics of wave packets with stochastic incorporation of Vlasov equation

International Nuclear Information System (INIS)

Ono, Akira; Horiuchi, Hisashi.

1996-01-01

The first purpose of this report is to present an extended AMD model which can generally describe such minor branching processes by removing the restriction on the one-body distribution function. This is done not by generalizing the wave packets to arbitrary single-particle wave functions but by representing the diffused and/or deformed wave packet as an ensemble of Gaussian wave packets. In other words, stochastic displacements are given to the wave packets in phase space so that the ensemble-average of the time evolution of the one-body distribution function is essentially equivalent to the solution of Vlasov equation which does not have any restriction on the shape of wave packets. This new model is called AMD-V. Although AMD-V is equivalent to Vlasov equation in the instantaneous time evolution of the one-body distribution function for an AMD wave function, AMD-V describes the branching into channels and the fluctuation of the mean field which are caused by the spreading or the splitting of the single-particle wave function. The second purpose of this report is to show the drastic effect of this new stochastic process of wave packet splitting on the dynamics of heavy ion collisions, especially in the fragmentation mechanism. We take the 40 Ca + 40 Ca system at the incident energy 35 MeV/nucleon. It will be shown that the reproduction of data by the AMD-V calculation is surprisingly good. We will see that the effect of the wave packet diffusion is crucially important to remove the spurious binary feature of the AMD calculation and to enable the multi-fragment final state. (J.P.N.)

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)

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).

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.

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.

Environmental Stochasticity and the Speed of Evolution

Science.gov (United States)

Danino, Matan; Kessler, David A.; Shnerb, Nadav M.

2018-03-01

Biological populations are subject to two types of noise: demographic stochasticity due to fluctuations in the reproductive success of individuals, and environmental variations that affect coherently the relative fitness of entire populations. The rate in which the average fitness of a community increases has been considered so far using models with pure demographic stochasticity; here we present some theoretical considerations and numerical results for the general case where environmental variations are taken into account. When the competition is pairwise, fitness fluctuations are shown to reduce the speed of evolution, while under global competition the speed increases due to environmental stochasticity.

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.

Stochastic interpretation of the advection-diffusion equation and its relevance to bed load transport

Science.gov (United States)

Ancey, C.; Bohorquez, P.; Heyman, J.

2015-12-01

The advection-diffusion equation is one of the most widespread equations in physics. It arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Phenomenological laws are usually sufficient to derive this equation and interpret its terms. Stochastic models can also be used to derive it, with the significant advantage that they provide information on the statistical properties of particle activity. These models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. Among these stochastic models, the most common approach consists of random walk models. For instance, they have been used to model the random displacement of tracers in rivers. Here we explore an alternative approach, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. Birth-death Markov processes are well suited to this objective. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received no attention. We therefore look into the possibility of deriving the advection-diffusion equation (with a source term) within the framework of birth-death Markov processes. We show that in the continuum limit (when the cell size becomes vanishingly small), we can derive an advection-diffusion equation for particle activity. Yet while this derivation is formally valid in the continuum limit, it runs into difficulty in practical applications involving cells or meshes of finite length. Indeed, within our stochastic framework, particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due

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.

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.

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...

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

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

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.

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.

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.

Boussinesq evolution equations

DEFF Research Database (Denmark)

Bredmose, Henrik; Schaffer, H.; Madsen, Per A.

2004-01-01

This paper deals with the possibility of using methods and ideas from time domain Boussinesq formulations in the corresponding frequency domain formulations. We term such frequency domain models "evolution equations". First, we demonstrate that the numerical efficiency of the deterministic...... Boussinesq evolution equations of Madsen and Sorensen [Madsen, P.A., Sorensen, O.R., 1993. Bound waves and triad interactions in shallow water. Ocean Eng. 20 359-388] can be improved by using Fast Fourier Transforms to evaluate the nonlinear terms. For a practical example of irregular waves propagating over...... a submerged bar, it is demonstrated that evolution equations utilising FFT can be solved around 100 times faster than the corresponding time domain model. Use of FFT provides an efficient bridge between the frequency domain and the time domain. We utilise this by adapting the surface roller model for wave...

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.

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)....

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...

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 approach to chemical evolution

International Nuclear Information System (INIS)

Copi, C.J.

1997-01-01

Observations of elemental abundances in the Galaxy have repeatedly shown an intrinsic scatter as a function of time and metallicity. The standard approach to chemical evolution does not attempt to address this scatter in abundances since only the mean evolution is followed. In this work, the scatter is addressed via a stochastic approach to solving chemical evolution models. Three simple chemical evolution scenarios are studied using this stochastic approach: a closed box model, an infall model, and an outflow model. These models are solved for the solar neighborhood in a Monte Carlo fashion. The evolutionary history of one particular region is determined randomly based on the star formation rate and the initial mass function. Following the evolution in an ensemble of such regions leads to the predicted spread in abundances expected, based solely on different evolutionary histories of otherwise identical regions. In this work, 13 isotopes are followed, including the light elements, the CNO elements, a few α-elements, and iron. It is found that the predicted spread in abundances for a 10 5 M circle-dot region is in good agreement with observations for the α-elements. For CN, the agreement is not as good, perhaps indicating the need for more physics input for low-mass stellar evolution. Similarly for the light elements, the predicted scatter is quite small, which is in contradiction to the observations of 3 He in HII regions. The models are tuned for the solar neighborhood so that good agreement with HII regions is not expected. This has important implications for low-mass stellar evolution and on using chemical evolution to determine the primordial light-element abundances in order to test big bang nucleosynthesis. copyright 1997 The American Astronomical Society

Master equations and the theory of stochastic path integrals

Science.gov (United States)

Weber, Markus F.; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a ‘generating functional’, which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a ‘forward’ and a ‘backward’ path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from

Master equations and the theory of stochastic path integrals.

Science.gov (United States)

Weber, Markus F; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a 'generating functional', which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a 'forward' and a 'backward' path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from them. Upon

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.

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.

Kinetic theory of age-structured stochastic birth-death processes

Science.gov (United States)

Greenman, Chris D.; Chou, Tom

2016-01-01

Classical age-structured mass-action models such as the McKendrick-von Foerster equation have been extensively studied but are unable to describe stochastic fluctuations or population-size-dependent birth and death rates. Stochastic theories that treat semi-Markov age-dependent processes using, e.g., the Bellman-Harris equation do not resolve a population's age structure and are unable to quantify population-size dependencies. Conversely, current theories that include size-dependent population dynamics (e.g., mathematical models that include carrying capacity such as the logistic equation) cannot be easily extended to take into account age-dependent birth and death rates. In this paper, we present a systematic derivation of a new, fully stochastic kinetic theory for interacting age-structured populations. By defining multiparticle probability density functions, we derive a hierarchy of kinetic equations for the stochastic evolution of an aging population undergoing birth and death. We show that the fully stochastic age-dependent birth-death process precludes factorization of the corresponding probability densities, which then must be solved by using a Bogoliubov--Born--Green--Kirkwood--Yvon-like hierarchy. Explicit solutions are derived in three limits: no birth, no death, and steady state. These are then compared with their corresponding mean-field results. Our results generalize both deterministic models and existing master equation approaches by providing an intuitive and efficient way to simultaneously model age- and population-dependent stochastic dynamics applicable to the study of demography, stem cell dynamics, and disease evolution.

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

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

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).

Stochastic Averaging Principle for Spatial Birth-and-Death Evolutions in the Continuum

Science.gov (United States)

Friesen, Martin; Kondratiev, Yuri

2018-06-01

We study a spatial birth-and-death process on the phase space of locally finite configurations Γ^+ × Γ^- over R}^d. Dynamics is described by an non-equilibrium evolution of states obtained from the Fokker-Planck equation and associated with the Markov operator L^+(γ ^-) + 1/ɛ L^-, ɛ > 0. Here L^- describes the environment process on Γ^- and L^+(γ ^-) describes the system process on Γ^+, where γ ^- indicates that the corresponding birth-and-death rates depend on another locally finite configuration γ ^- \\in Γ^-. We prove that, for a certain class of birth-and-death rates, the corresponding Fokker-Planck equation is well-posed, i.e. there exists a unique evolution of states μ _t^{ɛ } on Γ^+ × Γ^-. Moreover, we give a sufficient condition such that the environment is ergodic with exponential rate. Let μ _{inv} be the invariant measure for the environment process on Γ^-. In the main part of this work we establish the stochastic averaging principle, i.e. we prove that the marginal of μ _t^{ɛ } onto Γ^+ converges weakly to an evolution of states on {Γ}^+ associated with the averaged Markov birth-and-death operator {\\overline{L}} = \\int _{Γ}^- L^+(γ ^-)d μ _{inv}(γ ^-).

Stochastic Averaging Principle for Spatial Birth-and-Death Evolutions in the Continuum

Science.gov (United States)

Friesen, Martin; Kondratiev, Yuri

2018-04-01

We study a spatial birth-and-death process on the phase space of locally finite configurations Γ^+ × Γ^- over R^d . Dynamics is described by an non-equilibrium evolution of states obtained from the Fokker-Planck equation and associated with the Markov operator L^+(γ ^-) + 1/ɛ L^- , ɛ > 0 . Here L^- describes the environment process on Γ^- and L^+(γ ^-) describes the system process on Γ^+ , where γ ^- indicates that the corresponding birth-and-death rates depend on another locally finite configuration γ ^- \\in Γ^- . We prove that, for a certain class of birth-and-death rates, the corresponding Fokker-Planck equation is well-posed, i.e. there exists a unique evolution of states μ _t^{ɛ } on Γ^+ × Γ^- . Moreover, we give a sufficient condition such that the environment is ergodic with exponential rate. Let μ _{inv} be the invariant measure for the environment process on Γ^- . In the main part of this work we establish the stochastic averaging principle, i.e. we prove that the marginal of μ _t^{ɛ } onto Γ^+ converges weakly to an evolution of states on Γ^+ associated with the averaged Markov birth-and-death operator \\overline{L} = \\int _{Γ}^-}L^+(γ ^-)d μ _{inv}(γ ^-).

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.

Nonlinear evolution equations

CERN Document Server

Uraltseva, N N

1995-01-01

This collection focuses on nonlinear problems in partial differential equations. Most of the papers are based on lectures presented at the seminar on partial differential equations and mathematical physics at St. Petersburg University. Among the topics explored are the existence and properties of solutions of various classes of nonlinear evolution equations, nonlinear imbedding theorems, bifurcations of solutions, and equations of mathematical physics (Navier-Stokes type equations and the nonlinear Schrödinger equation). The book will be useful to researchers and graduate students working in p

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.

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.

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.

A stochastic Friedman Universe with dissipation

International Nuclear Information System (INIS)

Gruszczak, J.

1985-01-01

A probabilistic measure is constructed for the radiation-filled Friedman Universe with bulk viscosity and the equation of state perturbed by a ''white noise''. The corresponding Fokker-Planck equation is solved. In the stochastic evolution singularities turn out to be irrelevant. 3 refs., 1 fig. (author)

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

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.

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.

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.

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

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.

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)

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.

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 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.

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

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.

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...

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...

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.

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.

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.

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,

Space-time-modulated stochastic processes

Science.gov (United States)

Giona, Massimiliano

2017-10-01

Starting from the physical problem associated with the Lorentzian transformation of a Poisson-Kac process in inertial frames, the concept of space-time-modulated stochastic processes is introduced for processes possessing finite propagation velocity. This class of stochastic processes provides a two-way coupling between the stochastic perturbation acting on a physical observable and the evolution of the physical observable itself, which in turn influences the statistical properties of the stochastic perturbation during its evolution. The definition of space-time-modulated processes requires the introduction of two functions: a nonlinear amplitude modulation, controlling the intensity of the stochastic perturbation, and a time-horizon function, which modulates its statistical properties, providing irreducible feedback between the stochastic perturbation and the physical observable influenced by it. The latter property is the peculiar fingerprint of this class of models that makes them suitable for extension to generic curved-space times. Considering Poisson-Kac processes as prototypical examples of stochastic processes possessing finite propagation velocity, the balance equations for the probability density functions associated with their space-time modulations are derived. Several examples highlighting the peculiarities of space-time-modulated processes are thoroughly analyzed.

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.

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.

Stochastic Effects for the Reaction-Duffing Equation with Wick-Type Product

Directory of Open Access Journals (Sweden)

Jin Hyuk Choi

2016-01-01

Full Text Available We construct new explicit solutions of the Wick-type stochastic reaction-Duffing equation arising from mathematical physics with the help of the white noise theory and the system technique. Based on these exact solutions, we also discuss the influences of stochastic effects for dynamical behaviors according to functions h1(t, h2(t, and Brownian motion B(t which are the solitary wave group velocities.

Stochastic equations theory and applications in acoustics, hydrodynamics, magnetohydrodynamics, and radiophysics

CERN Document Server

Klyatskin, Valery I

2015-01-01

This monograph set presents a consistent and self-contained framework of stochastic dynamic systems with maximal possible completeness. Volume 1 presents the basic concepts, exact results, and asymptotic approximations of the theory of stochastic equations on the basis of the developed functional approach. This approach offers a possibility of both obtaining exact solutions to stochastic problems for a number of models of fluctuating parameters and constructing various asymptotic buildings. Ideas of statistical topography are used to discuss general issues of generating coherent structures from chaos with probability one, i.e., almost in every individual realization of random parameters. The general theory is illustrated with certain problems and applications of stochastic mathematical physics in various fields such as mechanics, hydrodynamics, magnetohydrodynamics, acoustics, optics, and radiophysics.  

Computational stochastic model of ions implantation

Energy Technology Data Exchange (ETDEWEB)

Zmievskaya, Galina I., E-mail: zmi@gmail.ru; Bondareva, Anna L., E-mail: bal310775@yandex.ru [M.V. Keldysh Institute of Applied Mathematics RAS, 4,Miusskaya sq., 125047 Moscow (Russian Federation); Levchenko, Tatiana V., E-mail: tatlevchenko@mail.ru [VNII Geosystem Russian Federal Center, Varshavskoye roadway, 8, Moscow (Russian Federation); Maino, Giuseppe, E-mail: giuseppe.maino@enea.it [Scuola di Lettere e BeniCulturali, University di Bologna, sede di Ravenna, via Mariani 5, 48100 Ravenna (Italy)

2015-03-10

Implantation flux ions into crystal leads to phase transition /PT/ 1-st kind. Damaging lattice is associated with processes clustering vacancies and gaseous bubbles as well their brownian motion. System of stochastic differential equations /SDEs/ Ito for evolution stochastic dynamical variables corresponds to the superposition Wiener processes. The kinetic equations in partial derivatives /KE/, Kolmogorov-Feller and Einstein-Smolukhovskii, were formulated for nucleation into lattice of weakly soluble gases. According theory, coefficients of stochastic and kinetic equations uniquely related. Radiation stimulated phase transition are characterized by kinetic distribution functions /DFs/ of implanted clusters versus their sizes and depth of gas penetration into lattice. Macroscopic parameters of kinetics such as the porosity and stress calculated in thin layers metal/dielectric due to Xe{sup ++} irradiation are attracted as example. Predictions of porosity, important for validation accumulation stresses in surfaces, can be applied at restoring of objects the cultural heritage.

A Dynamic BI–Orthogonal Field Equation Approach to Efficient Bayesian Inversion

Directory of Open Access Journals (Sweden)

Tagade Piyush M.

2017-06-01

Full Text Available This paper proposes a novel computationally efficient stochastic spectral projection based approach to Bayesian inversion of a computer simulator with high dimensional parametric and model structure uncertainty. The proposed method is based on the decomposition of the solution into its mean and a random field using a generic Karhunen-Loève expansion. The random field is represented as a convolution of separable Hilbert spaces in stochastic and spatial dimensions that are spectrally represented using respective orthogonal bases. In particular, the present paper investigates generalized polynomial chaos bases for the stochastic dimension and eigenfunction bases for the spatial dimension. Dynamic orthogonality is used to derive closed-form equations for the time evolution of mean, spatial and the stochastic fields. The resultant system of equations consists of a partial differential equation (PDE that defines the dynamic evolution of the mean, a set of PDEs to define the time evolution of eigenfunction bases, while a set of ordinary differential equations (ODEs define dynamics of the stochastic field. This system of dynamic evolution equations efficiently propagates the prior parametric uncertainty to the system response. The resulting bi-orthogonal expansion of the system response is used to reformulate the Bayesian inference for efficient exploration of the posterior distribution. The efficacy of the proposed method is investigated for calibration of a 2D transient diffusion simulator with an uncertain source location and diffusivity. The computational efficiency of the method is demonstrated against a Monte Carlo method and a generalized polynomial chaos approach.

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 models with power-law tails the equation X = AX + B

CERN Document Server

Buraczewski, Dariusz; Mikosch, Thomas

2016-01-01

In this monograph the authors give a systematic approach to the probabilistic properties of the fixed point equation X=AX+B. A probabilistic study of the stochastic recurrence equation X_t=A_tX_{t-1}+B_t for real- and matrix-valued random variables A_t, where (A_t,B_t) constitute an iid sequence, is provided. The classical theory for these equations, including the existence and uniqueness of a stationary solution, the tail behavior with special emphasis on power law behavior, moments and support, is presented. The authors collect recent asymptotic results on extremes, point processes, partial sums (central limit theory with special emphasis on infinite variance stable limit theory), large deviations, in the univariate and multivariate cases, and they further touch on the related topics of smoothing transforms, regularly varying sequences and random iterative systems. The text gives an introduction to the Kesten-Goldie theory for stochastic recurrence equations of the type X_t=A_tX_{t-1}+B_t. It provides the c...

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 ...

Probability density function evolution of power systems subject to stochastic variation of renewable energy

Science.gov (United States)

Wei, J. Q.; Cong, Y. C.; Xiao, M. Q.

2018-05-01

As renewable energies are increasingly integrated into power systems, there is increasing interest in stochastic analysis of power systems.Better techniques should be developed to account for the uncertainty caused by penetration of renewables and consequently analyse its impacts on stochastic stability of power systems. In this paper, the Stochastic Differential Equations (SDEs) are used to represent the evolutionary behaviour of the power systems. The stationary Probability Density Function (PDF) solution to SDEs modelling power systems excited by Gaussian white noise is analysed. Subjected to such random excitation, the Joint Probability Density Function (JPDF) solution to the phase angle and angular velocity is governed by the generalized Fokker-Planck-Kolmogorov (FPK) equation. To solve this equation, the numerical method is adopted. Special measure is taken such that the generalized FPK equation is satisfied in the average sense of integration with the assumed PDF. Both weak and strong intensities of the stochastic excitations are considered in a single machine infinite bus power system. The numerical analysis has the same result as the one given by the Monte Carlo simulation. Potential studies on stochastic behaviour of multi-machine power systems with random excitations are discussed at the end.

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 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 ...

Evolution equations for Killing fields

International Nuclear Information System (INIS)

Coll, B.

1977-01-01

The problem of finding necessary and sufficient conditions on the Cauchy data for Einstein equations which insure the existence of Killing fields in a neighborhood of an initial hypersurface has been considered recently by Berezdivin, Coll, and Moncrief. Nevertheless, it can be shown that the evolution equations obtained in all these cases are of nonstrictly hyperbolic type, and, thus, the Cauchy data must belong to a special class of functions. We prove here that, for the vacuum and Einstein--Maxwell space--times and in a coordinate independent way, one can always choose, as evolution equations for the Killing fields, a strictly hyperbolic system: The above theorems can be thus extended to all Cauchy data for which the Einstein evolution problem has been proved to be well set

Conditional stability in determination of initial data for stochastic parabolic equations

International Nuclear Information System (INIS)

Yuan, Ganghua

2017-01-01

In this paper, we solve two kinds of inverse problems in determination of the initial data for stochastic parabolic equations. One is determination of the initial data by lateral boundary observation on arbitrary portion of the boundary, the second one is determination of the initial data by internal observation in a subregion inside the domain. We obtain conditional stability for the two kinds of inverse problems. To prove the results, we estimate the initial data by a terminal observation near the initial time, then we estimate this terminal observation by lateral boundary observation on arbitrary portion of the boundary or internal observation in a subregion inside the domain. To achieve those goals, we derive several new Carleman estimates for stochastic parabolic equations in this paper. (paper)

Conditional stability in determination of initial data for stochastic parabolic equations

Science.gov (United States)

Yuan, Ganghua

2017-03-01

In this paper, we solve two kinds of inverse problems in determination of the initial data for stochastic parabolic equations. One is determination of the initial data by lateral boundary observation on arbitrary portion of the boundary, the second one is determination of the initial data by internal observation in a subregion inside the domain. We obtain conditional stability for the two kinds of inverse problems. To prove the results, we estimate the initial data by a terminal observation near the initial time, then we estimate this terminal observation by lateral boundary observation on arbitrary portion of the boundary or internal observation in a subregion inside the domain. To achieve those goals, we derive several new Carleman estimates for stochastic parabolic equations in this paper.

Stochastic methods for the description of multiparticle production

International Nuclear Information System (INIS)

Carruthers, P.

1984-01-01

Dynamical questions in the evolution of excited hadronic matter are reviewed, with emphasis on KNO scaling and its possible violation. It is suggested that the KNO distributions is described by a stochastic evolution of the Fokker-Planck type related to underlying field theory by coupled rate equations approximated by Langevin equations with noise. Refined correlation analysis of data, especially the use of intensity interferometry techniques, is recommended for data analysis. 26 references

Variance estimates for transport in stochastic media by means of the master equation

International Nuclear Information System (INIS)

Pautz, S. D.; Franke, B. C.; Prinja, A. K.

2013-01-01

The master equation has been used to examine properties of transport in stochastic media. It has been shown previously that not only may the Levermore-Pomraning (LP) model be derived from the master equation for a description of ensemble-averaged transport quantities, but also that equations describing higher-order statistical moments may be obtained. We examine in greater detail the equations governing the second moments of the distribution of the angular fluxes, from which variances may be computed. We introduce a simple closure for these equations, as well as several models for estimating the variances of derived transport quantities. We revisit previous benchmarks for transport in stochastic media in order to examine the error of these new variance models. We find, not surprisingly, that the errors in these variance estimates are at least as large as the corresponding estimates of the average, and sometimes much larger. We also identify patterns in these variance estimates that may help guide the construction of more accurate models. (authors)

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.

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...

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.

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.

Stochastic Effects; Application in Nuclear Physics

International Nuclear Information System (INIS)

Mazonka, O.

2000-04-01

Stochastic effects in nuclear physics refer to the study of the dynamics of nuclear systems evolving under stochastic equations of motion. In this dissertation we restrict our attention to classical scattering models. We begin with introduction of the model of nuclear dynamics and deterministic equations of evolution. We apply a Langevin approach - an additional property of the model, which reflect the statistical nature of low energy nuclear behaviour. We than concentrate our attention on the problem of calculating tails of distribution functions, which actually is the problem of calculating probabilities of rare outcomes. Two general strategies are proposed. Result and discussion follow. Finally in the appendix we consider stochastic effects in nonequilibrium systems. A few exactly solvable models are presented. For one model we show explicitly that stochastic behaviour in a microscopic description can lead to ordered collective effects on the macroscopic scale. Two others are solved to confirm the predictions of the fluctuation theorem. (author)

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)

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.

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.

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 ...

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...

Fractional diffusion equation with distributed-order material derivative. Stochastic foundations

International Nuclear Information System (INIS)

Magdziarz, M; Teuerle, M

2017-01-01

In this paper, we present the stochastic foundations of fractional dynamics driven by the fractional material derivative of distributed-order type. Before stating our main result, we present the stochastic scenario which underlies the dynamics given by the fractional material derivative. Then we introduce the Lévy walk process of distributed-order type to establish our main result, which is the scaling limit of the considered process. It appears that the probability density function of the scaling limit process fulfills, in a weak sense, the fractional diffusion equation with the material derivative of distributed-order type. (paper)

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.

Time evolution of one-dimensional gapless models from a domain wall initial state: stochastic Loewner evolution continued?

International Nuclear Information System (INIS)

Calabrese, Pasquale; Hagendorf, Christian; Doussal, Pierre Le

2008-01-01

We study the time evolution of quantum one-dimensional gapless systems evolving from initial states with a domain wall. We generalize the path integral imaginary time approach that together with boundary conformal field theory allows us to derive the time and space dependence of general correlation functions. The latter are explicitly obtained for the Ising universality class, and the typical behavior of one- and two-point functions is derived for the general case. Possible connections with the stochastic Loewner evolution are discussed and explicit results for one-point time dependent averages are obtained for generic κ for boundary conditions corresponding to stochastic Loewner evolution. We use this set of results to predict the time evolution of the entanglement entropy and obtain the universal constant shift due to the presence of a domain wall in the initial state

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.

Stochastic quantization of field theories on the lattice and supersymmetrical models

International Nuclear Information System (INIS)

Aldazabal, Gerardo.

1984-01-01

Several aspects of the stochastic quantization method are considered. Specifically, field theories on the lattice and supersymmetrical models are studied. A non-linear sigma model is studied firstly, and it is shown that it is possible to obtain evolution equations written directly for invariant quantities. These ideas are generalized to obtain Langevin equations for the Wilson loops of non-abelian lattice gauge theories U (N) and SU (N). In order to write these equations, some different ways of introducing the constraints which the fields must satisfy are discussed. It is natural to have a strong coupling expansion in these equations. The correspondence with quantum field theory is established, and it is noticed that at all orders in the perturbation theory, Langevin equations reduce to Schwinger-Dyson equations. From another point of view, stochastic quantization is applied to large N matrix models on the lattice. As a result, a simple and systematic way of building reduced models is found. Referring to stochastic quantization in supersymmetric theories, a simple supersymmetric model is studied. It is shown that it is possible to write an evolution equation for the superfield wich leads to quantum field theory results in equilibrium. As the Langevin equation preserves supersymmetry, the property of dimensional reduction known for the quantum model is shown to be valid at all times. (M.E.L.) [es

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.

Large Deviations for Stochastic Tamed 3D Navier-Stokes Equations

International Nuclear Information System (INIS)

Roeckner, Michael; Zhang, Tusheng; Zhang Xicheng

2010-01-01

In this paper, using weak convergence method, we prove a large deviation principle of Freidlin-Wentzell type for the stochastic tamed 3D Navier-Stokes equations driven by multiplicative noise, which was investigated in (Roeckner and Zhang in Probab. Theory Relat. Fields 145(1-2), 211-267, 2009).

A stochastic Galerkin method for the Euler equations with Roe variable transformation

KAUST Repository

Pettersson, Per; Iaccarino, Gianluca; Nordströ m, Jan

2014-01-01

The Euler equations subject to uncertainty in the initial and boundary conditions are investigated via the stochastic Galerkin approach. We present a new fully intrusive method based on a variable transformation of the continuous equations. Roe variables are employed to get quadratic dependence in the flux function and a well-defined Roe average matrix that can be determined without matrix inversion.In previous formulations based on generalized polynomial chaos expansion of the physical variables, the need to introduce stochastic expansions of inverse quantities, or square roots of stochastic quantities of interest, adds to the number of possible different ways to approximate the original stochastic problem. We present a method where the square roots occur in the choice of variables, resulting in an unambiguous problem formulation.The Roe formulation saves computational cost compared to the formulation based on expansion of conservative variables. Moreover, the Roe formulation is more robust and can handle cases of supersonic flow, for which the conservative variable formulation fails to produce a bounded solution. For certain stochastic basis functions, the proposed method can be made more effective and well-conditioned. This leads to increased robustness for both choices of variables. We use a multi-wavelet basis that can be chosen to include a large number of resolution levels to handle more extreme cases (e.g. strong discontinuities) in a robust way. For smooth cases, the order of the polynomial representation can be increased for increased accuracy. © 2013 Elsevier Inc.

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

Continuum modeling of ion-beam eroded surfaces under normal incidence: Impact of stochastic fluctuations

International Nuclear Information System (INIS)

Dreimann, Karsten; Linz, Stefan J.

2010-01-01

Graphical abstract: Deterministic surface pattern (left) and its stochastic counterpart (right) arising in a stochastic damped Kuramoto-Sivashinsky equation that serves as a model equation for ion-beam eroded surfaces and is systematically investigated. - Abstract: Using a recently proposed field equation for the surface evolution of ion-beam eroded semiconductor target materials under normal incidence, we systematically explore the impact of additive stochastic fluctuations that are permanently present during the erosion process. Specifically, we investigate the dependence of the surface roughness, the underlying pattern forming properties and the bifurcation behavior on the strength of the fluctuations.

Stochastic weighted particle methods for population balance equations with coagulation, fragmentation and spatial inhomogeneity

International Nuclear Information System (INIS)

Lee, Kok Foong; Patterson, Robert I.A.; Wagner, Wolfgang; Kraft, Markus

2015-01-01

Graphical abstract: -- Highlights: •Problems concerning multi-compartment population balance equations are studied. •A class of fragmentation weight transfer functions is presented. •Three stochastic weighted algorithms are compared against the direct simulation algorithm. •The numerical errors of the stochastic solutions are assessed as a function of fragmentation rate. •The algorithms are applied to a multi-dimensional granulation model. -- Abstract: This paper introduces stochastic weighted particle algorithms for the solution of multi-compartment population balance equations. In particular, it presents a class of fragmentation weight transfer functions which are constructed such that the number of computational particles stays constant during fragmentation events. The weight transfer functions are constructed based on systems of weighted computational particles and each of it leads to a stochastic particle algorithm for the numerical treatment of population balance equations. Besides fragmentation, the algorithms also consider physical processes such as coagulation and the exchange of mass with the surroundings. The numerical properties of the algorithms are compared to the direct simulation algorithm and an existing method for the fragmentation of weighted particles. It is found that the new algorithms show better numerical performance over the two existing methods especially for systems with significant amount of large particles and high fragmentation rates.

Stochastic weighted particle methods for population balance equations with coagulation, fragmentation and spatial inhomogeneity

Energy Technology Data Exchange (ETDEWEB)

Lee, Kok Foong [Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA (United Kingdom); Patterson, Robert I.A.; Wagner, Wolfgang [Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstraße 39, 10117 Berlin (Germany); Kraft, Markus, E-mail: mk306@cam.ac.uk [Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA (United Kingdom); School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, Singapore, 637459 (Singapore)

2015-12-15

Graphical abstract: -- Highlights: •Problems concerning multi-compartment population balance equations are studied. •A class of fragmentation weight transfer functions is presented. •Three stochastic weighted algorithms are compared against the direct simulation algorithm. •The numerical errors of the stochastic solutions are assessed as a function of fragmentation rate. •The algorithms are applied to a multi-dimensional granulation model. -- Abstract: This paper introduces stochastic weighted particle algorithms for the solution of multi-compartment population balance equations. In particular, it presents a class of fragmentation weight transfer functions which are constructed such that the number of computational particles stays constant during fragmentation events. The weight transfer functions are constructed based on systems of weighted computational particles and each of it leads to a stochastic particle algorithm for the numerical treatment of population balance equations. Besides fragmentation, the algorithms also consider physical processes such as coagulation and the exchange of mass with the surroundings. The numerical properties of the algorithms are compared to the direct simulation algorithm and an existing method for the fragmentation of weighted particles. It is found that the new algorithms show better numerical performance over the two existing methods especially for systems with significant amount of large particles and high fragmentation rates.

A new evolution equation

International Nuclear Information System (INIS)

Laenen, E.

1995-01-01

We propose a new evolution equation for the gluon density relevant for the region of small x B . It generalizes the GLR equation and allows deeper penetration in dense parton systems than the GLR equation does. This generalization consists of taking shadowing effects more comprehensively into account by including multigluon correlations, and allowing for an arbitrary initial gluon distribution in a hadron. We solve the new equation for fixed α s . We find that the effects of multigluon correlations on the deep-inelastic structure function are small. (orig.)

Analysis and reconstruction of stochastic coupled map lattice models

International Nuclear Information System (INIS)

Coca, Daniel; Billings, Stephen A.

2003-01-01

The Letter introduces a general stochastic coupled lattice map model together with an algorithm to estimate the nodal equations involved based only on a small set of observable variables and in the presence of stochastic perturbations. More general forms of the Frobenius-Perron and the transfer operators, which describe the evolution of densities under the action of the CML transformation, are derived

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...

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.

Stochastic coalescence in finite systems: an algorithm for the numerical solution of the multivariate master equation.

Science.gov (United States)

Alfonso, Lester; Zamora, Jose; Cruz, Pedro

2015-04-01

The stochastic approach to coagulation considers the coalescence process going in a system of a finite number of particles enclosed in a finite volume. Within this approach, the full description of the system can be obtained from the solution of the multivariate master equation, which models the evolution of the probability distribution of the state vector for the number of particles of a given mass. Unfortunately, due to its complexity, only limited results were obtained for certain type of kernels and monodisperse initial conditions. In this work, a novel numerical algorithm for the solution of the multivariate master equation for stochastic coalescence that works for any type of kernels and initial conditions is introduced. The performance of the method was checked by comparing the numerically calculated particle mass spectrum with analytical solutions obtained for the constant and sum kernels, with an excellent correspondence between the analytical and numerical solutions. In order to increase the speedup of the algorithm, software parallelization techniques with OpenMP standard were used, along with an implementation in order to take advantage of new accelerator technologies. Simulations results show an important speedup of the parallelized algorithms. This study was funded by a grant from Consejo Nacional de Ciencia y Tecnologia de Mexico SEP-CONACYT CB-131879. The authors also thanks LUFAC® Computacion SA de CV for CPU time and all the support provided.

Characterisation of Exponential Convergence to Nonequilibrium Limits for Stochastic Volterra Equations

Directory of Open Access Journals (Sweden)

John A. D. Appleby

2008-01-01

Full Text Available This paper considers necessary and sufficient conditions for the solution of a stochastically and deterministically perturbed Volterra equation to converge exponentially to a nonequilibrium and nontrivial limit. Convergence in an almost sure and pth mean sense is obtained.

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...

Projection scheme for a reflected stochastic heat equation with additive noise

Science.gov (United States)

Higa, Arturo Kohatsu; Pettersson, Roger

2005-02-01

We consider a projection scheme as a numerical solution of a reflected stochastic heat equation driven by a space-time white noise. Convergence is obtained via a discrete contraction principle and known convergence results for numerical solutions of parabolic variational inequalities.

Lie symmetries for systems of evolution equations

Science.gov (United States)

Paliathanasis, Andronikos; Tsamparlis, Michael

2018-01-01

The Lie symmetries for a class of systems of evolution equations are studied. The evolution equations are defined in a bimetric space with two Riemannian metrics corresponding to the space of the independent and dependent variables of the differential equations. The exact relation of the Lie symmetries with the collineations of the bimetric space is determined.

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

Stochastic cooling for beginners

International Nuclear Information System (INIS)

Moehl, D.

1984-01-01

These two lectures have been prepared to give a simple introduction to the principles. In Part I we try to explain stochastic cooling using the time-domain picture which starts from the pulse response of the system. In Part II the discussion is repeated, looking more closely at the frequency-domain response. An attempt is made to familiarize the beginners with some of the elementary cooling equations, from the 'single particle case' up to equations which describe the evolution of the particle distribution. (orig.)

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.

Stochastic substitute for coupled rate equations in the modeling of highly ionized transient plasmas

International Nuclear Information System (INIS)

Eliezer, S.; Falquina, R.; Minguez, E.

1994-01-01

Plasmas produced by intense laser pulses incident on solid targets often do not satisfy the conditions for local thermodynamic equilibrium, and so cannot be modeled by transport equations relying on equations of state. A proper description involves an excessively large number of coupled rate equations connecting many quantum states of numerous species having different degrees of ionization. Here we pursue a recent suggestion to model the plasma by a few dominant states perturbed by a stochastic driving force. The driving force is taken to be a Poisson impulse process, giving a Langevin equation which is equivalent to a Fokker-Planck equation for the probability density governing the distribution of electron density. An approximate solution to the Langevin equation permits calculation of the characteristic relaxation rate. An exact stationary solution to the Fokker-Planck equation is given as a function of the strength of the stochastic driving force. This stationary solution is used, along with a Laplace transform, to convert the Fokker-Planck equation to one of Schroedinger type. We consider using the classical Hamiltonian formalism and the WKB method to obtain the time-dependent solution

On the small time asymptotics of 3D stochastic primitive equations

OpenAIRE

Dong, Zhao; Zhang, Rangrang

2017-01-01

In this paper, we establish a small time large deviation principle for the strong solution of 3D stochastic primitive equations driven by multiplicative noise. Both the small noise and the small, but highly nonlinear, unbounded nonlinear terms should be taken into consideration.

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

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

Continuous stochastic approach to birth and death processes and co-operative behaviour of systems far from equilibrium

Energy Technology Data Exchange (ETDEWEB)

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

1986-09-11

The continuous stochastic formalism for the description of systems with birth and death processes randomly distributed in space is developed with the use of local birth and death operators and local generalization of the corresponding Chapman-Kolmogorov equation. The functional stochastic equation for the evolution of the probability functional is derived and its modifications for evolution of the characteristic functional and the first passage time problem are given. The corresponding evolution equations for equal-time correlators are also derived. The results are generalized then on the exothermic and endothermic chemical reactions. As examples of the particular applications of the results the small fluctuations near stable equilibrium state and fluctuations in mono-molecular reactions, Lotka-Volterra model, Schloegl reaction and brusselator are considered. It is shown that the two-dimensional Lotka-Volterra model may exhibit synergetic phase transition analogous to the topological transition of the Kosterlitz-Thouless-Berezinskii type. At the end of the paper some general consequences from stochastic evolution of the birth and death processes are discussed and the arguments on their importance in evolution of populations, cellular dynamics and in applications to various chemical and biological problems are presented.

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

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...

Decomposition of a hierarchy of nonlinear evolution equations

International Nuclear Information System (INIS)

Geng Xianguo

2003-01-01

The generalized Hamiltonian structures for a hierarchy of nonlinear evolution equations are established with the aid of the trace identity. Using the nonlinearization approach, the hierarchy of nonlinear evolution equations is decomposed into a class of new finite-dimensional Hamiltonian systems. The generating function of integrals and their generator are presented, based on which the finite-dimensional Hamiltonian systems are proved to be completely integrable in the Liouville sense. As an application, solutions for the hierarchy of nonlinear evolution equations are reduced to solving the compatible Hamiltonian systems of ordinary differential equations

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.).

Advanced functional evolution equations and inclusions

CERN Document Server

Benchohra, Mouffak

2015-01-01

This book presents up-to-date results on abstract evolution equations and differential inclusions in infinite dimensional spaces. It covers equations with time delay and with impulses, and complements the existing literature in functional differential equations and inclusions. The exposition is devoted to both local and global mild solutions for some classes of functional differential evolution equations and inclusions, and other densely and non-densely defined functional differential equations and inclusions in separable Banach spaces or in Fréchet spaces. The tools used include classical fixed points theorems and the measure-of non-compactness, and each chapter concludes with a section devoted to notes and bibliographical remarks. This monograph is particularly useful for researchers and graduate students studying pure and applied mathematics, engineering, biology and all other applied sciences.

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.

Mean Field Games for Stochastic Growth with Relative Utility

Energy Technology Data Exchange (ETDEWEB)

Huang, Minyi, E-mail: mhuang@math.carleton.ca [Carleton University, School of Mathematics and Statistics (Canada); Nguyen, Son Luu, E-mail: sonluu.nguyen@upr.edu [University of Puerto Rico, Department of Mathematics (United States)

2016-12-15

This paper considers continuous time stochastic growth-consumption optimization in a mean field game setting. The individual capital stock evolution is determined by a Cobb–Douglas production function, consumption and stochastic depreciation. The individual utility functional combines an own utility and a relative utility with respect to the population. The use of the relative utility reflects human psychology, leading to a natural pattern of mean field interaction. The fixed point equation of the mean field game is derived with the aid of some ordinary differential equations. Due to the relative utility interaction, our performance analysis depends on some ratio based approximation error estimate.

Mean Field Games for Stochastic Growth with Relative Utility

International Nuclear Information System (INIS)

Huang, Minyi; Nguyen, Son Luu

2016-01-01

This paper considers continuous time stochastic growth-consumption optimization in a mean field game setting. The individual capital stock evolution is determined by a Cobb–Douglas production function, consumption and stochastic depreciation. The individual utility functional combines an own utility and a relative utility with respect to the population. The use of the relative utility reflects human psychology, leading to a natural pattern of mean field interaction. The fixed point equation of the mean field game is derived with the aid of some ordinary differential equations. Due to the relative utility interaction, our performance analysis depends on some ratio based approximation error estimate.

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.

Population density equations for stochastic processes with memory kernels

Science.gov (United States)

Lai, Yi Ming; de Kamps, Marc

2017-06-01

We present a method for solving population density equations (PDEs)-a mean-field technique describing homogeneous populations of uncoupled neurons—where the populations can be subject to non-Markov noise for arbitrary distributions of jump sizes. The method combines recent developments in two different disciplines that traditionally have had limited interaction: computational neuroscience and the theory of random networks. The method uses a geometric binning scheme, based on the method of characteristics, to capture the deterministic neurodynamics of the population, separating the deterministic and stochastic process cleanly. We can independently vary the choice of the deterministic model and the model for the stochastic process, leading to a highly modular numerical solution strategy. We demonstrate this by replacing the master equation implicit in many formulations of the PDE formalism by a generalization called the generalized Montroll-Weiss equation—a recent result from random network theory—describing a random walker subject to transitions realized by a non-Markovian process. We demonstrate the method for leaky- and quadratic-integrate and fire neurons subject to spike trains with Poisson and gamma-distributed interspike intervals. We are able to model jump responses for both models accurately to both excitatory and inhibitory input under the assumption that all inputs are generated by one renewal process.

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

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

Beam life-time with intrabeam scattering and stochastic cooling

International Nuclear Information System (INIS)

Wei, J.; Ruggiero, A.G.

1991-01-01

A transport equation has been derived in terms of the longitudinal action variable to describe the time evolution of the longitudinal density distribution of a bunched hadron beam in the presence of intrabeam scattering and stochastic cooling. A computer program has been developed to numerically solve this equation. Both beam loss and bunch-shape evolution have been investigated for the 197 Au 79+ beams during the 10-hour storage in the Relativistic Heavy Ion Collider currently under construction at the Brookhaven National Laboratory. 9 refs., 1 fig

On the maximal noise for stochastic and QCD travelling waves

International Nuclear Information System (INIS)

Peschanski, Robi

2008-01-01

Using the relation of a set of nonlinear Langevin equations to reaction-diffusion processes, we note the existence of a maximal strength of the noise for the stochastic travelling wave solutions of these equations. Its determination is obtained using the field-theoretical analysis of branching-annihilation random walks near the directed percolation transition. We study its consequence for the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation. For the related Langevin equation modeling the quantum chromodynamic nonlinear evolution of gluon density with rapidity, the physical maximal-noise limit may appear before the directed percolation transition, due to a shift in the travelling-wave speed. In this regime, an exact solution is known from a coalescence process. Universality and other open problems and applications are discussed in the outlook

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

Parameter-free resolution of the superposition of stochastic signals

Energy Technology Data Exchange (ETDEWEB)

Scholz, Teresa, E-mail: tascholz@fc.ul.pt [Center for Theoretical and Computational Physics, University of Lisbon (Portugal); Raischel, Frank [Center for Geophysics, IDL, University of Lisbon (Portugal); Closer Consulting, Av. Eng. Duarte Pacheco Torre 1 15" 0, 1070-101 Lisboa (Portugal); Lopes, Vitor V. [DEIO-CIO, University of Lisbon (Portugal); UTEC–Universidad de Ingeniería y Tecnología, Lima (Peru); Lehle, Bernd; Wächter, Matthias; Peinke, Joachim [Institute of Physics and ForWind, Carl-von-Ossietzky University of Oldenburg, Oldenburg (Germany); Lind, Pedro G. [Institute of Physics and ForWind, Carl-von-Ossietzky University of Oldenburg, Oldenburg (Germany); Institute of Physics, University of Osnabrück, Osnabrück (Germany)

2017-01-30

This paper presents a direct method to obtain the deterministic and stochastic contribution of the sum of two independent stochastic processes, one of which is an Ornstein–Uhlenbeck process and the other a general (non-linear) Langevin process. The method is able to distinguish between the stochastic processes, retrieving their corresponding stochastic evolution equations. This framework is based on a recent approach for the analysis of multidimensional Langevin-type stochastic processes in the presence of strong measurement (or observational) noise, which is here extended to impose neither constraints nor parameters and extract all coefficients directly from the empirical data sets. Using synthetic data, it is shown that the method yields satisfactory results.

Perturbative approach to non-Markovian stochastic Schroedinger equations

International Nuclear Information System (INIS)

Gambetta, Jay; Wiseman, H.M.

2002-01-01

In this paper we present a perturbative procedure that allows one to numerically solve diffusive non-Markovian stochastic Schroedinger equations, for a wide range of memory functions. To illustrate this procedure numerical results are presented for a classically driven two-level atom immersed in an environment with a simple memory function. It is observed that as the order of the perturbation is increased the numerical results for the ensemble average state ρ red (t) approach the exact reduced state found via Imamog-barlu ' s enlarged system method [Phys. Rev. A 50, 3650 (1994)

Stochastic processes, slaves and supersymmetry

International Nuclear Information System (INIS)

Drummond, I T; Horgan, R R

2012-01-01

We extend the work of Tănase-Nicola and Kurchan on the structure of diffusion processes and the associated supersymmetry algebra by examining the responses of a simple statistical system to external disturbances of various kinds. We consider both the stochastic differential equations (SDEs) for the process and the associated diffusion equation. The influence of the disturbances can be understood by augmenting the original SDE with an equation for slave variables. The evolution of the slave variables describes the behaviour of line elements carried along in the stochastic flow. These line elements, together with the associated surface and volume elements constructed from them, provide the basis of the supersymmetry properties of the theory. For ease of visualization, and in order to emphasize a helpful electromagnetic analogy, we work in three dimensions. The results are all generalizable to higher dimensions and can be specialized to one and two dimensions. The electromagnetic analogy is a useful starting point for calculating asymptotic results at low temperature that can be compared with direct numerical evaluations. We also examine the problems that arise in a direct numerical simulation of the stochastic equation together with the slave equations. We pay special attention to the dependence of the slave variable statistics on temperature. We identify in specific models the critical temperature below which the slave variable distribution ceases to have a variance and consider the effect on estimates of susceptibilities. (paper)

Front propagation and clustering in the stochastic nonlocal Fisher equation

Science.gov (United States)

Ganan, Yehuda A.; Kessler, David A.

2018-04-01

In this work, we study the problem of front propagation and pattern formation in the stochastic nonlocal Fisher equation. We find a crossover between two regimes: a steadily propagating regime for not too large interaction range and a stochastic punctuated spreading regime for larger ranges. We show that the former regime is well described by the heuristic approximation of the system by a deterministic system where the linear growth term is cut off below some critical density. This deterministic system is seen not only to give the right front velocity, but also predicts the onset of clustering for interaction kernels which give rise to stable uniform states, such as the Gaussian kernel, for sufficiently large cutoff. Above the critical cutoff, distinct clusters emerge behind the front. These same features are present in the stochastic model for sufficiently small carrying capacity. In the latter, punctuated spreading, regime, the population is concentrated on clusters, as in the infinite range case, which divide and separate as a result of the stochastic noise. Due to the finite interaction range, if a fragment at the edge of the population separates sufficiently far, it stabilizes as a new cluster, and the processes begins anew. The deterministic cutoff model does not have this spreading for large interaction ranges, attesting to its purely stochastic origins. We show that this mode of spreading has an exponentially small mean spreading velocity, decaying with the range of the interaction kernel.

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...

On Continuous Selection Sets of Non-Lipschitzian Quantum Stochastic Evolution Inclusions

Directory of Open Access Journals (Sweden)

Sheila Bishop

2015-01-01

Full Text Available We establish existence of a continuous selection of multifunctions associated with quantum stochastic evolution inclusions under a general Lipschitz condition. The coefficients here are multifunctions but not necessarily Lipschitz.

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.

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.

Stochastic evolution of the Universe: A possible dynamical process ...

Indian Academy of Sciences (India)

C Sivakumar

2017-12-11

Dec 11, 2017 ... https://doi.org/10.1007/s12043-017-1491-z. Stochastic evolution of the Universe: A possible dynamical process leading to fractal structures. C SIVAKUMAR. Department of Physics, Maharaja's College, Ernakulam 682 011, India. E-mail: thrisivc@yahoo.com. MS received 6 July 2016; revised 26 June 2017; ...

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...

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...

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.

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.

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.

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.

Reduced equations of motion for quantum systems driven by diffusive Markov processes.

Science.gov (United States)

Sarovar, Mohan; Grace, Matthew D

2012-09-28

The expansion of a stochastic Liouville equation for the coupled evolution of a quantum system and an Ornstein-Uhlenbeck process into a hierarchy of coupled differential equations is a useful technique that simplifies the simulation of stochastically driven quantum systems. We expand the applicability of this technique by completely characterizing the class of diffusive Markov processes for which a useful hierarchy of equations can be derived. The expansion of this technique enables the examination of quantum systems driven by non-Gaussian stochastic processes with bounded range. We present an application of this extended technique by simulating Stark-tuned Förster resonance transfer in Rydberg atoms with nonperturbative position fluctuations.

Systems of evolution equations and the singular perturbation method

International Nuclear Information System (INIS)

Mika, J.

Several fundamental theorems are presented important for the solution of linear evolution equations in the Banach space. The algorithm is deduced extending the solution of the system of singularly perturbed evolution equations into an asymptotic series with respect to a small positive parameter. The asymptotic convergence is shown of an approximate solution to the accurate solution. Singularly perturbed evolution equations of the resonance type were analysed. The special role is considered of the asymptotic equivalence of P1 equations obtained as the first order approximation if the spherical harmonics method is applied to the linear Boltzmann equation, and the diffusion equations of the linear transport theory where the small parameter approaches zero. (J.B.)

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.

Stochastic quantization of topological field theory: generalized Langevin equation with memory kernel

International Nuclear Information System (INIS)

Menezes, G.; Svaiter, N.F.

2006-04-01

We use the method of stochastic quantization in a topological field theory defined in an Euclidean space, assuming a Langevin equation with a memory kernel. We show that our procedure for the Abelian Chern-Simons theory converges regardless of the nature of the Chern-Simons coefficient. (author)

The Asymptotic Behaviour of a Stochastic 3D LANS-α Model

International Nuclear Information System (INIS)

Caraballo, Tomas; Marquez-Duran, Antonio M.; Real, Jose

2006-01-01

The long-time behaviour of a stochastic 3D LANS-α model on a bounded domain is analysed. First, we reformulate the model as an abstract problem. Next, we establish sufficient conditions ensuring the existence of stationary (steady state) solutions of this abstract nonlinear stochastic evolution equation, and study the stability properties of the model. Finally, we analyse the effects produced by stochastic perturbations in the deterministic version of the system (persistence of exponential stability as well as possible stabilisation effects produced by the noise). The general results are applied to our stochastic LANS-α system throughout the paper

Quantitative study of quasi-one-dimensional Bose gas experiments via the stochastic Gross-Pitaevskii equation

International Nuclear Information System (INIS)

Cockburn, S. P.; Gallucci, D.; Proukakis, N. P.

2011-01-01

The stochastic Gross-Pitaevskii equation is shown to be an excellent model for quasi-one-dimensional Bose gas experiments, accurately reproducing the in situ density profiles recently obtained in the experiments of Trebbia et al.[Phys. Rev. Lett. 97, 250403 (2006)] and van Amerongen et al.[Phys. Rev. Lett. 100, 090402 (2008)] and the density fluctuation data reported by Armijo et al.[Phys. Rev. Lett. 105, 230402 (2010)]. To facilitate such agreement, we propose and implement a quasi-one-dimensional extension to the one-dimensional stochastic Gross-Pitaevskii equation for the low-energy, axial modes, while atoms in excited transverse modes are treated as independent ideal Bose gases.

Calibration of a stochastic health evolution model using NHIS data

Science.gov (United States)

Gupta, Aparna; Li, Zhisheng

2011-10-01

This paper presents and calibrates an individual's stochastic health evolution model. In this health evolution model, the uncertainty of health incidents is described by a stochastic process with a finite number of possible outcomes. We construct a comprehensive health status index (HSI) to describe an individual's health status, as well as a health risk factor system (RFS) to classify individuals into different risk groups. Based on the maximum likelihood estimation (MLE) method and the method of nonlinear least squares fitting, model calibration is formulated in terms of two mixed-integer nonlinear optimization problems. Using the National Health Interview Survey (NHIS) data, the model is calibrated for specific risk groups. Longitudinal data from the Health and Retirement Study (HRS) is used to validate the calibrated model, which displays good validation properties. The end goal of this paper is to provide a model and methodology, whose output can serve as a crucial component of decision support for strategic planning of health related financing and risk management.

Stochastic bosonization for a d ≥ 3 Fermi system

International Nuclear Information System (INIS)

Accardi, L.; Lu, Y.G.; Mastropietro, V.

1997-01-01

We consider a system of fermions interacting via an external field and we prove, in d ≥ 3, that a suitable collective operator, bilinear in the fermionic fields, in the stochastic limit becomes a boson quantum brownian motion. The evolution operator after the limit satisfies a quantum stochastic differential equation, in which the imaginary part of the Ito correction is the ground state shift while its real part is the lifetime of the ground state. (orig.)

Numerical Resolution of N-dimensional Fokker-Planck stochastic equations; Resolucion Numerica de Ecuaciones Estocasticas de tipo Fokker-Planck en Varias Dimensiones

Energy Technology Data Exchange (ETDEWEB)

Garcia-Olivares, R A; Munoz Roldan, A

1992-07-01

This document describes the use of a library of programs able to solve stochastic Fokker-Planck equations in a N-dimensional space. The input data are essentially: (i) the initial distribution of the stochastic variable, (ii) the drift and fluctuation coefficients as a function of the state (which can be obtained from the transition probabilities between neighboring states) and (iii) some parameters controlling the run. A last version of the library accepts sources and sinks defined in the states space. The output is the temporal evolution of the probability distribution in the space defined by a N-dimensional grid. Some applications and readings in Synergetic, Self-Organization, transport phenomena, Ecology and other fields are suggested. If the probability distribution is interpreted as a distribution of particles then the codes can be used to solve the N-dimensional problem of advection-diffusion. (Author) 16 refs.

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.

Physical entropy, information entropy and their evolution equations

Institute of Scientific and Technical Information of China (English)

无

2001-01-01

Inspired by the evolution equation of nonequilibrium statistical physics entropy and the concise statistical formula of the entropy production rate, we develop a theory of the dynamic information entropy and build a nonlinear evolution equation of the information entropy density changing in time and state variable space. Its mathematical form and physical meaning are similar to the evolution equation of the physical entropy: The time rate of change of information entropy density originates together from drift, diffusion and production. The concise statistical formula of information entropy production rate is similar to that of physical entropy also. Furthermore, we study the similarity and difference between physical entropy and information entropy and the possible unification of the two statistical entropies, and discuss the relationship among the principle of entropy increase, the principle of equilibrium maximum entropy and the principle of maximum information entropy as well as the connection between them and the entropy evolution equation.

Stochastic star formation and the evolution of galaxies

International Nuclear Information System (INIS)

Seiden, P.E.; Schulman, L.S.; Gerola, H.

1979-01-01

The mechanism of stochastic self-propagating star formation has previously been invoked to explain the origin of spiral arms in galaxies. In this paper we extend the application of this mechanism to account for the diversity of morphological types and the evolution of galaxies. The new property that arises from consideration of this mechanism is that the rate of star formation exhibits the critical behavior of a phase transition. This is a general property of the system and is not strongly dependent on the details of the star--interstellar gas interaction. Examination of the properties of this phase transition provides a general scenario for the evolution of galaxies and the origin of the various morphological types

Stochastic beam dynamics in storage rings

International Nuclear Information System (INIS)

Pauluhn, A.

1993-12-01

In this thesis several approaches to stochastic dynamics in storage rings are investigated. In the first part the theory of stochastic differential equations and Fokker-Planck equations is used to describe the processes which have been assumed to be Markov processes. The mathematical theory of Markov processes is well known. Nevertheless, analytical solutions can be found only in special cases and numerical algorithms are required. Several numerical integration schemes for stochastic differential equations will therefore be tested in analytical solvable examples and then applied to examples from accelerator physics. In particular the stochastically perturbed synchrotron motion is treated. For the special case of a double rf system several perturbation theoretical methods for deriving the Fokker-Planck equation in the action variable are used and compared with numerical results. The second part is concerned with the dynamics of electron storage rings. Due to the synchrotron radiation the electron motion is influenced by damping and exciting forces. An algorithm for the computation of the density function in the phase space of such a dissipative stochastically excited system is introduced. The density function contains all information of a process, e.g. it determines the beam dimensions and the lifetime of a stored electron beam. The new algorithm consists in calculating a time propagator for the density function. By means of this propagator the time evolution of the density is modelled very computing time efficient. The method is applied to simple models of the beam-beam interaction (one-dimensional, round beams) and the results of the density calculations are compared with results obtained from multiparticle tracking. Furthermore some modifications of the algorithm are introduced to improve its efficiency concerning computing time and storage requirements. Finally, extensions to two-dimensional beam-beam models are described. (orig.)

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...

Aggregation patterns from nonlocal interactions: Discrete stochastic and continuum modeling

KAUST Repository

Hackett-Jones, Emily J.

2012-04-17

Conservation equations governed by a nonlocal interaction potential generate aggregates from an initial uniform distribution of particles. We address the evolution and formation of these aggregating steady states when the interaction potential has both attractive and repulsive singularities. Currently, no existence theory for such potentials is available. We develop and compare two complementary solution methods, a continuous pseudoinverse method and a discrete stochastic lattice approach, and formally show a connection between the two. Interesting aggregation patterns involving multiple peaks for a simple doubly singular attractive-repulsive potential are determined. For a swarming Morse potential, characteristic slow-fast dynamics in the scaled inverse energy is observed in the evolution to steady state in both the continuous and discrete approaches. The discrete approach is found to be remarkably robust to modifications in movement rules, related to the potential function. The comparable evolution dynamics and steady states of the discrete model with the continuum model suggest that the discrete stochastic approach is a promising way of probing aggregation patterns arising from two- and three-dimensional nonlocal interaction conservation equations. © 2012 American Physical Society.

Nonequilibrium steady state in open quantum systems: Influence action, stochastic equation and power balance

International Nuclear Information System (INIS)

Hsiang, J.-T.; Hu, B.L.

2015-01-01

The existence and uniqueness of a steady state for nonequilibrium systems (NESS) is a fundamental subject and a main theme of research in statistical mechanics for decades. For Gaussian systems, such as a chain of classical harmonic oscillators connected at each end to a heat bath, and for classical anharmonic oscillators under specified conditions, definitive answers exist in the form of proven theorems. Answering this question for quantum many-body systems poses a challenge for the present. In this work we address this issue by deriving the stochastic equations for the reduced system with self-consistent backaction from the two baths, calculating the energy flow from one bath to the chain to the other bath, and exhibiting a power balance relation in the total (chain + baths) system which testifies to the existence of a NESS in this system at late times. Its insensitivity to the initial conditions of the chain corroborates to its uniqueness. The functional method we adopt here entails the use of the influence functional, the coarse-grained and stochastic effective actions, from which one can derive the stochastic equations and calculate the average values of physical variables in open quantum systems. This involves both taking the expectation values of quantum operators of the system and the distributional averages of stochastic variables stemming from the coarse-grained environment. This method though formal in appearance is compact and complete. It can also easily accommodate perturbative techniques and diagrammatic methods from field theory. Taken all together it provides a solid platform for carrying out systematic investigations into the nonequilibrium dynamics of open quantum systems and quantum thermodynamics. -- Highlights: •Nonequilibrium steady state (NESS) for interacting quantum many-body systems. •Derivation of stochastic equations for quantum oscillator chain with two heat baths. •Explicit calculation of the energy flow from one bath to the

Simulation of Stochastic Processes by Coupled ODE-PDE

Science.gov (United States)

Zak, Michail

2008-01-01

A document discusses the emergence of randomness in solutions of coupled, fully deterministic ODE-PDE (ordinary differential equations-partial differential equations) due to failure of the Lipschitz condition as a new phenomenon. It is possible to exploit the special properties of ordinary differential equations (represented by an arbitrarily chosen, dynamical system) coupled with the corresponding Liouville equations (used to describe the evolution of initial uncertainties in terms of joint probability distribution) in order to simulate stochastic processes with the proscribed probability distributions. The important advantage of the proposed approach is that the simulation does not require a random-number generator.

Solution of the chemical master equation by radial basis functions approximation with interface tracking

NARCIS (Netherlands)

Kryven, I.; Röblitz, S; Schütte, C.

2015-01-01

Background: The chemical master equation is the fundamental equation of stochastic chemical kinetics. This differential-difference equation describes temporal evolution of the probability density function for states of a chemical system. A state of the system, usually encoded as a vector, represents

Stochastic description of supersymmetric fields with values in a manifold

International Nuclear Information System (INIS)

Hoba, Z.

1986-01-01

This paper discusses the mathematical problem of the imaginary time quantum mechanics of a particle moving in Euclidean space as considered from the theory of diffusion processes. The diffusion process is defined by a stochastic equation; the equation describes the diffusion process as a time evolution of a Brownian particle in a force field. The paper considers a Brownian particle on a Riemannian manifold

Subordination principle for fractional evolution equations

NARCIS (Netherlands)

Bazhlekova, E.G.

2000-01-01

The abstract Cauchy problem for the fractional evolution equation Daa = Au, a > 0, (1) where A is a closed densely de??ned operator in a Banach space, is investigated. The subordination principle, presented earlier in [J. P r ??u s s, Evolutionary In- tegral Equations and Applications. Birkh??auser,

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.

On the history of a stochastic ansatz for solving the transport equation

International Nuclear Information System (INIS)

Williams, M.M.R.

2010-01-01

A very useful approximate tool for understanding the role of random material properties on solutions of the transport equation is described and its historical derivation given. The development of this stochastic tool, from its introduction by Randall, to its use in describing current problems involving dichotomic or pseudo-dichotomic Markov processes is discussed.

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 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...

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

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.

Stochastic modeling of stock price process induced from the conjugate heat equation

Science.gov (United States)

Paeng, Seong-Hun

2015-02-01

Currency can be considered as a ruler for values of commodities. Then the price is the measured value by the ruler. We can suppose that inflation and variation of exchange rate are caused by variation of the scale of the ruler. In geometry, variation of the scale means that the metric is time-dependent. The conjugate heat equation is the modified heat equation which satisfies the heat conservation law for the time-dependent metric space. We propose a new model of stock prices by using the stochastic process whose transition probability is determined by the kernel of the conjugate heat equation. Our model of stock prices shows how the volatility term is affected by inflation and exchange rate. This model modifies the Black-Scholes equation in light of inflation and exchange rate.

Electroweak evolution equations

International Nuclear Information System (INIS)

Ciafaloni, Paolo; Comelli, Denis

2005-01-01

Enlarging a previous analysis, where only fermions and transverse gauge bosons were taken into account, we write down infrared-collinear evolution equations for the Standard Model of electroweak interactions computing the full set of splitting functions. Due to the presence of double logs which are characteristic of electroweak interactions (Bloch-Nordsieck violation), new infrared singular splitting functions have to be introduced. We also include corrections related to the third generation Yukawa couplings

Generalized Langevin equation with colored noise description of the stochastic oscillations of accretion disks

International Nuclear Information System (INIS)

Harko, Tiberiu; Leung, Chun Sing; Mocanu, Gabriela

2014-01-01

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects interacting with their external medium based on a generalized Langevin equation with colored noise and on the fluctuation-dissipation theorems. The former accounts for the general memory and retarded effects of the frictional force. The presence of the memory effects influences the response of the disk to external random interactions, and it modifies the dynamical behavior of the disk, as well as the energy dissipation processes. The generalized Langevin equation of the motion of the disk in the vertical direction is studied numerically, and the vertical displacements, velocities, and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The power spectral distribution of the disk luminosity is also obtained. As a possible astrophysical application of the formalism we investigate the possibility that the intra-day variability of the active galactic nuclei may be due to the stochastic disk instabilities. The perturbations due to colored/nontrivially correlated noise induce a complicated disk dynamics, which could explain some astrophysical observational features related to disk variability. (orig.)

Generalized Langevin equation with colored noise description of the stochastic oscillations of accretion disks

Energy Technology Data Exchange (ETDEWEB)

Harko, Tiberiu [University College London, Department of Mathematics, London (United Kingdom); Leung, Chun Sing [Polytechnic University, Department of Applied Mathematics, Hong Kong (China); Mocanu, Gabriela [Babes-Bolyai University, Faculty of Physics, Cluj-Napoca (Romania)

2014-05-15

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects interacting with their external medium based on a generalized Langevin equation with colored noise and on the fluctuation-dissipation theorems. The former accounts for the general memory and retarded effects of the frictional force. The presence of the memory effects influences the response of the disk to external random interactions, and it modifies the dynamical behavior of the disk, as well as the energy dissipation processes. The generalized Langevin equation of the motion of the disk in the vertical direction is studied numerically, and the vertical displacements, velocities, and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The power spectral distribution of the disk luminosity is also obtained. As a possible astrophysical application of the formalism we investigate the possibility that the intra-day variability of the active galactic nuclei may be due to the stochastic disk instabilities. The perturbations due to colored/nontrivially correlated noise induce a complicated disk dynamics, which could explain some astrophysical observational features related to disk variability. (orig.)

Generalized Langevin equation with colored noise description of the stochastic oscillations of accretion disks

Science.gov (United States)

Harko, Tiberiu; Leung, Chun Sing; Mocanu, Gabriela

2014-05-01

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects interacting with their external medium based on a generalized Langevin equation with colored noise and on the fluctuation-dissipation theorems. The former accounts for the general memory and retarded effects of the frictional force. The presence of the memory effects influences the response of the disk to external random interactions, and it modifies the dynamical behavior of the disk, as well as the energy dissipation processes. The generalized Langevin equation of the motion of the disk in the vertical direction is studied numerically, and the vertical displacements, velocities, and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The power spectral distribution of the disk luminosity is also obtained. As a possible astrophysical application of the formalism we investigate the possibility that the intra-day variability of the active galactic nuclei may be due to the stochastic disk instabilities. The perturbations due to colored/nontrivially correlated noise induce a complicated disk dynamics, which could explain some astrophysical observational features related to disk variability.

EFT Beyond the Horizon: Stochastic Inflation and How Primordial Quantum Fluctuations Go Classical

CERN Document Server

Burgess, C P; Tasinato, G; Williams, M

2015-01-01

We identify the effective theory describing inflationary super-Hubble scales and show it to be a special case of effective field theories appropriate to open systems. Open systems allow information to be exchanged between the degrees of freedom of interest and those that are integrated out, such as for particles moving through a fluid. Strictly speaking they cannot in general be described by an effective lagrangian; rather the appropriate `low-energy' limit is instead a Lindblad equation describing the evolution of the density matrix of the slow degrees of freedom. We derive the equation relevant to super-Hubble modes of quantum fields in near-de Sitter spacetimes and derive two implications. We show the evolution of the diagonal density-matrix elements quickly approaches the Fokker-Planck equation of Starobinsky's stochastic inflationary picture. This provides an alternative first-principles derivation of this picture's stochastic noise and drift, as well as its leading corrections. (An application computes ...

Stochastic heat and Burgers equations and their singularities II - Analytical Properties and Limiting Distributions

CERN Document Server

Davies, I M; Zhao, H

2004-01-01

We study the inviscid limit, $\\mu\\to 0$, of the stochastic viscous Burgers equation, for the velocity field $v^{\\mu}(x,t)$, $t>0$, $x\\in\\mathbb R^d$,\\frac{\\partial{v^{\\mu}}}{\\partial{t}} + (v^{\\mu}\\cdot\

Environmental Noise Could Promote Stochastic Local Stability of Behavioral Diversity Evolution

Science.gov (United States)

Zheng, Xiu-Deng; Li, Cong; Lessard, Sabin; Tao, Yi

2018-05-01

In this Letter, we investigate stochastic stability in a two-phenotype evolutionary game model for an infinite, well-mixed population undergoing discrete, nonoverlapping generations. We assume that the fitness of a phenotype is an exponential function of its expected payoff following random pairwise interactions whose outcomes randomly fluctuate with time. We show that the stochastic local stability of a constant interior equilibrium can be promoted by the random environmental noise even if the system may display a complicated nonlinear dynamics. This result provides a new perspective for a better understanding of how environmental fluctuations may contribute to the evolution of behavioral diversity.

Using Equation-Free Computation to Accelerate Network-Free Stochastic Simulation of Chemical Kinetics.

Science.gov (United States)

Lin, Yen Ting; Chylek, Lily A; Lemons, Nathan W; Hlavacek, William S

2018-06-21

The chemical kinetics of many complex systems can be concisely represented by reaction rules, which can be used to generate reaction events via a kinetic Monte Carlo method that has been termed network-free simulation. Here, we demonstrate accelerated network-free simulation through a novel approach to equation-free computation. In this process, variables are introduced that approximately capture system state. Derivatives of these variables are estimated using short bursts of exact stochastic simulation and finite differencing. The variables are then projected forward in time via a numerical integration scheme, after which a new exact stochastic simulation is initialized and the whole process repeats. The projection step increases efficiency by bypassing the firing of numerous individual reaction events. As we show, the projected variables may be defined as populations of building blocks of chemical species. The maximal number of connected molecules included in these building blocks determines the degree of approximation. Equation-free acceleration of network-free simulation is found to be both accurate and efficient.

Stochastic symplectic and multi-symplectic methods for nonlinear Schrödinger equation with white noise dispersion

Energy Technology Data Exchange (ETDEWEB)

Cui, Jianbo, E-mail: jianbocui@lsec.cc.ac.cn [Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, 100190 (China); Hong, Jialin, E-mail: hjl@lsec.cc.ac.cn [Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, 100190 (China); Liu, Zhihui, E-mail: liuzhihui@lsec.cc.ac.cn [Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, 100190 (China); Zhou, Weien, E-mail: weienzhou@nudt.edu.cn [College of Science, National University of Defense Technology, Changsha 410073 (China)

2017-08-01

We indicate that the nonlinear Schrödinger equation with white noise dispersion possesses stochastic symplectic and multi-symplectic structures. Based on these structures, we propose the stochastic symplectic and multi-symplectic methods, which preserve the continuous and discrete charge conservation laws, respectively. Moreover, we show that the proposed methods are convergent with temporal order one in probability. Numerical experiments are presented to verify our theoretical results.

Stochastic symplectic and multi-symplectic methods for nonlinear Schrödinger equation with white noise dispersion

International Nuclear Information System (INIS)

Cui, Jianbo; Hong, Jialin; Liu, Zhihui; Zhou, Weien

2017-01-01

We indicate that the nonlinear Schrödinger equation with white noise dispersion possesses stochastic symplectic and multi-symplectic structures. Based on these structures, we propose the stochastic symplectic and multi-symplectic methods, which preserve the continuous and discrete charge conservation laws, respectively. Moreover, we show that the proposed methods are convergent with temporal order one in probability. Numerical experiments are presented to verify our theoretical results.

On a new series of integrable nonlinear evolution equations

International Nuclear Information System (INIS)

Ichikawa, Y.H.; Wadati, Miki; Konno, Kimiaki; Shimizu, Tohru.

1980-10-01

Recent results of our research are surveyed in this report. The derivative nonlinear Schroedinger equation for the circular polarized Alfven wave admits the spiky soliton solutions for the plane wave boundary condition. The nonlinear equation for complex amplitude associated with the carrier wave is shown to be a generalized nonlinear Schroedinger equation, having the ordinary cubic nonlinear term and the derivative of cubic nonlinear term. A generalized scheme of the inverse scattering transformation has confirmed that superposition of the A-K-N-S scheme and the K-N scheme for the component equations valids for the generalized nonlinear Schroedinger equation. Then, two types of new integrable nonlinear evolution equation have been derived from our scheme of the inverse scattering transformation. One is the type of nonlinear Schroedinger equation, while the other is the type of Korteweg-de Vries equation. Brief discussions are presented for physical phenomena, which could be accounted by the second type of the new integrable nonlinear evolution equation. Lastly, the stationary solitary wave solutions have been constructed for the integrable nonlinear evolution equation of the second type. These solutions have peculiar structure that they are singular and discrete. It is a new challenge to construct singular potentials by the inverse scattering transformation. (author)

Recognition of Equations Using a Two-Dimensional Stochastic Context-Free Grammar

Science.gov (United States)

Chou, Philip A.

1989-11-01

We propose using two-dimensional stochastic context-free grammars for image recognition, in a manner analogous to using hidden Markov models for speech recognition. The value of the approach is demonstrated in a system that recognizes printed, noisy equations. The system uses a two-dimensional probabilistic version of the Cocke-Younger-Kasami parsing algorithm to find the most likely parse of the observed image, and then traverses the corresponding parse tree in accordance with translation formats associated with each production rule, to produce eqn I troff commands for the imaged equation. In addition, it uses two-dimensional versions of the Inside/Outside and Baum re-estimation algorithms for learning the parameters of the grammar from a training set of examples. Parsing the image of a simple noisy equation currently takes about one second of cpu time on an Alliant FX/80.

Stochastic weighted particle methods for population balance equations

International Nuclear Information System (INIS)

Patterson, Robert I.A.; Wagner, Wolfgang; Kraft, Markus

2011-01-01

Highlights: → Weight transfer functions for Monte Carlo simulation of coagulation. → Efficient support for single-particle growth processes. → Comparisons to analytic solutions and soot formation problems. → Better numerical accuracy for less common particles. - Abstract: A class of coagulation weight transfer functions is constructed, each member of which leads to a stochastic particle algorithm for the numerical treatment of population balance equations. These algorithms are based on systems of weighted computational particles and the weight transfer functions are constructed such that the number of computational particles does not change during coagulation events. The algorithms also facilitate the simulation of physical processes that change single particles, such as growth, or other surface reactions. Four members of the algorithm family have been numerically validated by comparison to analytic solutions to simple problems. Numerical experiments have been performed for complex laminar premixed flame systems in which members of the class of stochastic weighted particle methods were compared to each other and to a direct simulation algorithm. Two of the weighted algorithms have been shown to offer performance advantages over the direct simulation algorithm in situations where interest is focused on the larger particles in a system. The extent of this advantage depends on the particular system and on the quantities of interest.

Stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable duffing oscillator and bifurcation of moment equation

International Nuclear Information System (INIS)

Zhang Guangjun; Xu Jianxue; Wang Jue; Yue Zhifeng; Zou Hailin

2009-01-01

In this paper stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable Duffing oscillator is analyzed by moment method. This kind of novel transition refers to the one among three potential well on two sides of bifurcation point of original system at the presence of internal noise. Several conclusions are drawn. First, the semi-analytical result of stochastic resonance induced by the novel random transitions of two-dimensional weak damping bistable Duffing oscillator can be obtained, and the semi-analytical result is qualitatively compatible with the one of Monte Carlo simulation. Second, a bifurcation of double-branch fixed point curves occurs in the moment equations with noise intensity as their bifurcation parameter. Third, the bifurcation of moment equations corresponds to stochastic resonance of original system. Finally, the mechanism of stochastic resonance is presented from another viewpoint through analyzing the energy transfer induced by the bifurcation of moment equation.

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...

Bound states of quarks calculated with stochastic integration of the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Salomon, M.

1992-07-01

We have computed the masses, wave functions and sea quark content of mesons in their ground state by integrating the Bethe-Salpeter equation with a stochastic algorithm. This method allows the inclusion of a large set of diagrams. Inspection of the kernel of the equation shows that q-q-bar pairs with similar constituent masses in a singlet spin state exhibit a high bound state which is not present in other pairs. The pion, kaon and eta belongs to this category. 19 refs., 2 figs., 2 tabs

Asymptotic behavior of non-autonomous stochastic parabolic equations with nonlinear Laplacian principal part

Directory of Open Access Journals (Sweden)

Bixiang Wang

2013-08-01

Full Text Available We prove the existence and uniqueness of random attractors for the p-Laplace equation driven simultaneously by non-autonomous deterministic and stochastic forcing. The nonlinearity of the equation is allowed to have a polynomial growth rate of any order which may be greater than p. We further establish the upper semicontinuity of random attractors as the intensity of noise approaches zero. In addition, we show the pathwise periodicity of random attractors when all non-autonomous deterministic forcing terms are time periodic.

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...

Existence and Uniqueness of Solutions to the Stochastic Porous Media Equations of Saturated Flows

International Nuclear Information System (INIS)

Ciotir, Ioana

2010-01-01

This paper proves the existence and uniqueness of nonnegative solutions for the stochastic porous media equations with multiplicative noise, infinite jump and discontinuous diffusivity function relevant in description of saturation processes in underground water infiltration in a bounded domain of R 3 .

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.

Moving interfaces and quasilinear parabolic evolution equations

CERN Document Server

Prüss, Jan

2016-01-01

In this monograph, the authors develop a comprehensive approach for the mathematical analysis of a wide array of problems involving moving interfaces. It includes an in-depth study of abstract quasilinear parabolic evolution equations, elliptic and parabolic boundary value problems, transmission problems, one- and two-phase Stokes problems, and the equations of incompressible viscous one- and two-phase fluid flows. The theory of maximal regularity, an essential element, is also fully developed. The authors present a modern approach based on powerful tools in classical analysis, functional analysis, and vector-valued harmonic analysis. The theory is applied to problems in two-phase fluid dynamics and phase transitions, one-phase generalized Newtonian fluids, nematic liquid crystal flows, Maxwell-Stefan diffusion, and a variety of geometric evolution equations. The book also includes a discussion of the underlying physical and thermodynamic principles governing the equations of fluid flows and phase transitions...

Spatial evolution equation of wind wave growth

Institute of Scientific and Technical Information of China (English)

WANG; Wei; (王; 伟); SUN; Fu; (孙; 孚); DAI; Dejun; (戴德君)

2003-01-01

Based on the dynamic essence of air-sea interactions, a feedback type of spatial evolution equation is suggested to match reasonably the growing process of wind waves. This simple equation involving the dominant factors of wind wave growth is able to explain the transfer of energy from high to low frequencies without introducing the concept of nonlinear wave-wave interactions, and the results agree well with observations. The rate of wave height growth derived in this dissertation is applicable to both laboratory and open sea, which solidifies the physical basis of using laboratory experiments to investigate the generation of wind waves. Thus the proposed spatial evolution equation provides a new approach for the research on dynamic mechanism of air-sea interactions and wind wave prediction.

Ito's formula in UMD Banach spaces and regularity of solution of the Zakai equation

NARCIS (Netherlands)

Brzezniak, Z.; Van Neerven, J.M.A.M.; Veraar, M.C.; Weis, L.

2008-01-01

Using the theory of stochastic integration for processes with values in a UMD Banach space developed recently by the authors, an Itô formula is proved which is applied to prove the existence of strong solutions for a class of stochastic evolution equations in UMD Banach spaces. The abstract results

Stochastic Ocean Predictions with Dynamically-Orthogonal Primitive Equations

Science.gov (United States)

Subramani, D. N.; Haley, P., Jr.; Lermusiaux, P. F. J.

2017-12-01

The coastal ocean is a prime example of multiscale nonlinear fluid dynamics. Ocean fields in such regions are complex and intermittent with unstationary heterogeneous statistics. Due to the limited measurements, there are multiple sources of uncertainties, including the initial conditions, boundary conditions, forcing, parameters, and even the model parameterizations and equations themselves. For efficient and rigorous quantification and prediction of these uncertainities, the stochastic Dynamically Orthogonal (DO) PDEs for a primitive equation ocean modeling system with a nonlinear free-surface are derived and numerical schemes for their space-time integration are obtained. Detailed numerical studies with idealized-to-realistic regional ocean dynamics are completed. These include consistency checks for the numerical schemes and comparisons with ensemble realizations. As an illustrative example, we simulate the 4-d multiscale uncertainty in the Middle Atlantic/New York Bight region during the months of Jan to Mar 2017. To provide intitial conditions for the uncertainty subspace, uncertainties in the region were objectively analyzed using historical data. The DO primitive equations were subsequently integrated in space and time. The probability distribution function (pdf) of the ocean fields is compared to in-situ, remote sensing, and opportunity data collected during the coincident POSYDON experiment. Results show that our probabilistic predictions had skill and are 3- to 4- orders of magnitude faster than classic ensemble schemes.

An interpolation between the wave and diffusion equations through the fractional evolution equations Dirac like

International Nuclear Information System (INIS)

Pierantozzi, T.; Vazquez, L.

2005-01-01

Through fractional calculus and following the method used by Dirac to obtain his well-known equation from the Klein-Gordon equation, we analyze a possible interpolation between the Dirac and the diffusion equations in one space dimension. We study the transition between the hyperbolic and parabolic behaviors by means of the generalization of the D'Alembert formula for the classical wave equation and the invariance under space and time inversions of the interpolating fractional evolution equations Dirac like. Such invariance depends on the values of the fractional index and is related to the nonlocal property of the time fractional differential operator. For this system of fractional evolution equations, we also find an associated conserved quantity analogous to the Hamiltonian for the classical Dirac case

Stochastic control theory dynamic programming principle

CERN Document Server

Nisio, Makiko

2015-01-01

This book offers a systematic introduction to the optimal stochastic control theory via the dynamic programming principle, which is a powerful tool to analyze control problems. First we consider completely observable control problems with finite horizons. Using a time discretization we construct a nonlinear semigroup related to the dynamic programming principle (DPP), whose generator provides the Hamilton–Jacobi–Bellman (HJB) equation, and we characterize the value function via the nonlinear semigroup, besides the viscosity solution theory. When we control not only the dynamics of a system but also the terminal time of its evolution, control-stopping problems arise. This problem is treated in the same frameworks, via the nonlinear semigroup. Its results are applicable to the American option price problem. Zero-sum two-player time-homogeneous stochastic differential games and viscosity solutions of the Isaacs equations arising from such games are studied via a nonlinear semigroup related to DPP (the min-ma...

Poisson Stochastic Process and Basic Schauder and Sobolev Estimates in the Theory of Parabolic Equations

Science.gov (United States)

Krylov, N. V.; Priola, E.

2017-09-01

We show, among other things, how knowing Schauder or Sobolev-space estimates for the one-dimensional heat equation allows one to derive their multidimensional analogs for equations with coefficients depending only on the time variable with the same constants as in the case of the one-dimensional heat equation. The method is quite general and is based on using the Poisson stochastic process. It also applies to equations involving non-local operators. It looks like no other methods are available at this time and it is a very challenging problem to find a purely analytical approach to proving such results.

Switching of bound vector solitons for the coupled nonlinear Schrödinger equations with nonhomogenously stochastic perturbations

International Nuclear Information System (INIS)

Sun Zhiyuan; Yu Xin; Liu Ying; Gao Yitian

2012-01-01

We investigate the dynamics of the bound vector solitons (BVSs) for the coupled nonlinear Schrödinger equations with the nonhomogenously stochastic perturbations added on their dispersion terms. Soliton switching (besides soliton breakup) can be observed between the two components of the BVSs. Rate of the maximum switched energy (absolute values) within the fixed propagation distance (about 10 periods of the BVSs) enhances in the sense of statistics when the amplitudes of stochastic perturbations increase. Additionally, it is revealed that the BVSs with enhanced coherence are more robust against the perturbations with nonhomogenous stochasticity. Diagram describing the approximate borders of the splitting and non-splitting areas is also given. Our results might be helpful in dynamics of the BVSs with stochastic noises in nonlinear optical fibers or with stochastic quantum fluctuations in Bose-Einstein condensates.

Switching of bound vector solitons for the coupled nonlinear Schrödinger equations with nonhomogenously stochastic perturbations

Science.gov (United States)

Sun, Zhi-Yuan; Gao, Yi-Tian; Yu, Xin; Liu, Ying

2012-12-01

We investigate the dynamics of the bound vector solitons (BVSs) for the coupled nonlinear Schrödinger equations with the nonhomogenously stochastic perturbations added on their dispersion terms. Soliton switching (besides soliton breakup) can be observed between the two components of the BVSs. Rate of the maximum switched energy (absolute values) within the fixed propagation distance (about 10 periods of the BVSs) enhances in the sense of statistics when the amplitudes of stochastic perturbations increase. Additionally, it is revealed that the BVSs with enhanced coherence are more robust against the perturbations with nonhomogenous stochasticity. Diagram describing the approximate borders of the splitting and non-splitting areas is also given. Our results might be helpful in dynamics of the BVSs with stochastic noises in nonlinear optical fibers or with stochastic quantum fluctuations in Bose-Einstein condensates.

Weierstrass Elliptic Function Solutions to Nonlinear Evolution Equations

International Nuclear Information System (INIS)

Yu Jianping; Sun Yongli

2008-01-01

This paper is based on the relations between projection Riccati equations and Weierstrass elliptic equation, combined with the Groebner bases in the symbolic computation. Then the novel method for constructing the Weierstrass elliptic solutions to the nonlinear evolution equations is given by using the above relations

Extended Jacobi Elliptic Function Rational Expansion Method and Its Application to (2+1)-Dimensional Stochastic Dispersive Long Wave System

International Nuclear Information System (INIS)

Song Lina; Zhang Hongqing

2007-01-01

In this work, by means of a generalized method and symbolic computation, we extend the Jacobi elliptic function rational expansion method to uniformly construct a series of stochastic wave solutions for stochastic evolution equations. To illustrate the effectiveness of our method, we take the (2+1)-dimensional stochastic dispersive long wave system as an example. We not only have obtained some known solutions, but also have constructed some new rational formal stochastic Jacobi elliptic function solutions.

Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations

International Nuclear Information System (INIS)

Wu, Fuke; Tian, Tianhai; Rawlings, James B.; Yin, George

2016-01-01

The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766–1793 (1996); ibid. 56, 1794–1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.

Emmy Noether and Linear Evolution Equations

Directory of Open Access Journals (Sweden)

P. G. L. Leach

2013-01-01

Full Text Available Noether’s Theorem relates the Action Integral of a Lagrangian with symmetries which leave it invariant and the first integrals consequent upon the variational principle and the existence of the symmetries. These each have an equivalent in the Schrödinger Equation corresponding to the Lagrangian and by extension to linear evolution equations in general. The implications of these connections are investigated.

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

Effective evolution equations from quantum mechanics

OpenAIRE

Leopold, Nikolai

2018-01-01

The goal of this thesis is to provide a mathematical rigorous derivation of the Schrödinger-Klein-Gordon equations, the Maxwell-Schrödinger equations and the defocusing cubic nonlinear Schrödinger equation in two dimensions. We study the time evolution of the Nelson model (with ultraviolet cutoff) in a limit where the number N of charged particles gets large while the coupling of each particle to the radiation field is of order N^{−1/2}. At time zero it is assumed that almost all charges a...

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....

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.

Effects of stochastic noise on dynamical decoupling procedures

Energy Technology Data Exchange (ETDEWEB)

Bernad, Jozsef Zsolt; Frydrych, Holger; Alber, Gernot [Institut fuer Angewandte Physik, Technische Universitaet Darmstadt, D-64289 Darmstadt (Germany)

2013-07-01

Dynamical decoupling is a well-established technique to protect quantum systems from unwanted influences of their environment by exercising active control. It has been used experimentally to drastically increase the lifetime of qubit states in various implementations. The efficiency of different dynamical decoupling schemes defines the lifetime. However, errors in control operations always limit this efficiency. We propose a stochastic model as a possible description of imperfect control pulses and discuss the impact of this kind of error on different decoupling schemes. In the limit of continuous control, i.e. if the number of pulses N → ∞, we derive a stochastic differential equation for the evolution of the density operator of the controlled system and its environment. In the context of this modified time evolution we discuss possibilities of protecting qubit states against environmental noise.

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).

A stochastic collocation method for the second order wave equation with a discontinuous random speed

KAUST Repository

Motamed, Mohammad; Nobile, Fabio; Tempone, Raul

2012-01-01

In this paper we propose and analyze a stochastic collocation method for solving the second order wave equation with a random wave speed and subjected to deterministic boundary and initial conditions. The speed is piecewise smooth in the physical

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.

Interpretation of non-Markovian stochastic Schroedinger equations as a hidden-variable theory

International Nuclear Information System (INIS)

Gambetta, Jay; Wiseman, H.M.

2003-01-01

Do diffusive non-Markovian stochastic Schroedinger equations (SSEs) for open quantum systems have a physical interpretation? In a recent paper [Phys. Rev. A 66, 012108 (2002)] we investigated this question using the orthodox interpretation of quantum mechanics. We found that the solution of a non-Markovian SSE represents the state the system would be in at that time if a measurement was performed on the environment at that time, and yielded a particular result. However, the linking of solutions at different times to make a trajectory is, we concluded, a fiction. In this paper we investigate this question using the modal (hidden variable) interpretation of quantum mechanics. We find that the noise function z(t) appearing in the non-Markovian SSE can be interpreted as a hidden variable for the environment. That is, some chosen property (beable) of the environment has a definite value z(t) even in the absence of measurement on the environment. The non-Markovian SSE gives the evolution of the state of the system 'conditioned' on this environment hidden variable. We present the theory for diffusive non-Markovian SSEs that have as their Markovian limit SSEs corresponding to homodyne and heterodyne detection, as well as one which has no Markovian limit

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.

From stochastic phase-space evolution to brownian motion in collective space

Energy Technology Data Exchange (ETDEWEB)

Benhassine, B. (Lab. de Physique Nucleaire/ CNRS et Univ. de Nantes, 44 Nantes (France)); Farine, M. (Lab. de Physique Nucleaire/ CNRS et Univ. de Nantes, 44 Nantes (France) Ecole Navale, Lamveoc-Loulmic, 29 Brest-Naval (France)); Hernandez, E.S. (Dept. de Fisica - Facultad de Ciencias Exactas y Naturales, Univ. de Buenos Aires (Argentina)); Idier, D. (Lab. de Physique Nucleaire/ CNRS et Univ. de Nantes, 44 Nantes (France)); Remaud, B. (Lab. de Physique Nucleaire/ CNRS et Univ. de Nantes, 44 Nantes (France)); Sebille, F. (Lab. de Physique Nucleaire/ CNRS et Univ. de Nantes, 44 Nantes (France))

1994-01-24

Within the framework of stochastic transport equations in phase space, we study the dynamics of fluctuations on collective variables in homogeneous fermion systems. The transport coefficients are formally deduced in the relaxation-time approximation and a general method to compute dynamically the dispersions of collective observables is proposed as a set of coupled equations: respectively, the BUU/Landau-Vlasov equation for the average phase-space trajectories and the equations for the averages and dispersions of the observables. Independently, we derive the general covariance matrix of phase-space fluctuations and then by projection, the dispersion on collective variables at equilibrium. Detailed numerical applications of the formalism are given; they show that the dynamics of fluctuations can be extracted from noisy numerical simulations and that the leading parameter for collective fluctuations is the excitation energy, whatever is its degree of thermalization. (orig.)

From stochastic phase-space evolution to brownian motion in collective space

International Nuclear Information System (INIS)

Benhassine, B.; Farine, M.; Hernandez, E.S.; Idier, D.; Remaud, B.; Sebille, F.

1994-01-01

Within the framework of stochastic transport equations in phase space, we study the dynamics of fluctuations on collective variables in homogeneous fermion systems. The transport coefficients are formally deduced in the relaxation-time approximation and a general method to compute dynamically the dispersions of collective observables is proposed as a set of coupled equations: respectively, the BUU/Landau-Vlasov equation for the average phase-space trajectories and the equations for the averages and dispersions of the observables. Independently, we derive the general covariance matrix of phase-space fluctuations and then by projection, the dispersion on collective variables at equilibrium. Detailed numerical applications of the formalism are given; they show that the dynamics of fluctuations can be extracted from noisy numerical simulations and that the leading parameter for collective fluctuations is the excitation energy, whatever is its degree of thermalization. (orig.)

Nonlocal higher order evolution equations

KAUST Repository

Rossi, Julio D.; Schö nlieb, Carola-Bibiane

2010-01-01

In this article, we study the asymptotic behaviour of solutions to the nonlocal operator ut(x, t)1/4(-1)n-1 (J*Id -1)n (u(x, t)), x ∈ ℝN, which is the nonlocal analogous to the higher order local evolution equation vt(-1)n-1(Δ)nv. We prove

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)

Linear stochastic neutron transport theory

International Nuclear Information System (INIS)

Lewins, J.

1978-01-01

A new and direct derivation of the Bell-Pal fundamental equation for (low power) neutron stochastic behaviour in the Boltzmann continuum model is given. The development includes correlation of particle emission direction in induced and spontaneous fission. This leads to generalizations of the backward and forward equations for the mean and variance of neutron behaviour. The stochastic importance for neutron transport theory is introduced and related to the conventional deterministic importance. Defining equations and moment equations are derived and shown to be related to the backward fundamental equation with the detector distribution of the operational definition of stochastic importance playing the role of an adjoint source. (author)

Effective evolution equations from many-body quantum mechanics

International Nuclear Information System (INIS)

Benedikter, Niels Patriz

2014-01-01

Systems of interest in physics often consist of a very large number of interacting particles. In certain physical regimes, effective non-linear evolution equations are commonly used as an approximation for making predictions about the time-evolution of such systems. Important examples are Bose-Einstein condensates of dilute Bose gases and degenerate Fermi gases. While the effective equations are well-known in physics, a rigorous justification is very difficult. However, a rigorous derivation is essential to precisely understand the range and the limits of validity and the quality of the approximation. In this thesis, we prove that the time evolution of Bose-Einstein condensates in the Gross-Pitaevskii regime can be approximated by the time-dependent Gross-Pitaevskii equation, a cubic non-linear Schroedinger equation. We then turn to fermionic systems and prove that the evolution of a degenerate Fermi gas can be approximated by the time-dependent Hartree-Fock equation (TDHF) under certain assumptions on the semiclassical structure of the initial data. Finally, we extend the latter result to fermions with relativistic kinetic energy. All our results provide explicit bounds on the error as the number of particles becomes large. A crucial methodical insight on bosonic systems is that correlations can be modeled by Bogolyubov transformations. We construct initial data appropriate for the Gross-Pitaevskii regime using a Bogolyubov transformation acting on a coherent state, which amounts to studying squeezed coherent states. As a crucial insight for fermionic systems, we point out a semiclassical structure in states close to the ground state of fermions in a trap. As a convenient language for studying the dynamics of fermionic systems, we use particle-hole transformations.

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.

Almost Surely Asymptotic Stability of Exact and Numerical Solutions for Neutral Stochastic Pantograph Equations

Directory of Open Access Journals (Sweden)

Zhanhua Yu

2011-01-01

Full Text Available We study the almost surely asymptotic stability of exact solutions to neutral stochastic pantograph equations (NSPEs, and sufficient conditions are obtained. Based on these sufficient conditions, we show that the backward Euler method (BEM with variable stepsize can preserve the almost surely asymptotic stability. Numerical examples are demonstrated for illustration.

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...

Physical interpretation of Monte Carlo wave-function and stochastic Schroedinger equation methods for cavity quantum electrodynamics

International Nuclear Information System (INIS)

Kist, Tarso B.L.; Orszag, M.; Davidovich, L.

1997-01-01

The dynamics of open system is frequently modeled in terms of a small system S coupled to a reservoir R, the last having an infinitely larger number of degree of freedom than S. Usually the dynamics of the S variables may be of interest, which can be studied using either Langevin equations, or master equations, or yet the path integral formulation. Useful alternatives for the master equation method are the Monte Carlo Wave-function method (MCWF), and Stochastic Schroedinger Equations (SSE's). The methods MCWF and SSE's recently experienced a fast development both in their theoretical background and applications to the study of the dissipative quantum systems dynamics in quantum optics. Even though these alternatives can be shown to be formally equivalent to the master equation approach, they are often regarded as mathematical tricks, with no relation to a concrete physical evolution of the system. The advantage of using them is that one has to deal with state vectors, instead of density matrices, thus reducing the total amount of matrix elements to be calculated. In this work, we consider the possibility of giving a physical interpretation to these methods, in terms of continuous measurements made on the evolving system. We show that physical realizations of the two methods are indeed possible, for a mode of the electromagnetic field in a cavity interacting with a continuum of modes corresponding to the field outside the cavity. Two schemes are proposed, consisting of a mode of the electromagnetic field interacting with a beam of Rydberg two-level atoms. In these schemes, the field mode plays the role of a small system and the atomic beam plays the role of a reservoir (infinitely larger number of degrees of freedom at finite temperature, the interaction between them being given by the Jaynes-Cummings model

Existence of solutions of abstract fractional impulsive semilinear evolution equations

Directory of Open Access Journals (Sweden)

K. Balachandran

2010-01-01

Full Text Available In this paper we prove the existence of solutions of fractional impulsive semilinear evolution equations in Banach spaces. A nonlocal Cauchy problem is discussed for the evolution equations. The results are obtained using fractional calculus and fixed point theorems. An example is provided to illustrate the theory.

Stochastic dynamics modeling solute transport in porous media modeling solute transport in porous media

CERN Document Server

Kulasiri, Don

2002-01-01

Most of the natural and biological phenomena such as solute transport in porous media exhibit variability which can not be modeled by using deterministic approaches. There is evidence in natural phenomena to suggest that some of the observations can not be explained by using the models which give deterministic solutions. Stochastic processes have a rich repository of objects which can be used to express the randomness inherent in the system and the evolution of the system over time. The attractiveness of the stochastic differential equations (SDE) and stochastic partial differential equations (SPDE) come from the fact that we can integrate the variability of the system along with the scientific knowledge pertaining to the system. One of the aims of this book is to explaim some useufl concepts in stochastic dynamics so that the scientists and engineers with a background in undergraduate differential calculus could appreciate the applicability and appropriateness of these developments in mathematics. The ideas ...

Hamiltonian structures of some non-linear evolution equations

International Nuclear Information System (INIS)

Tu, G.Z.

1983-06-01

The Hamiltonian structure of the O(2,1) non-linear sigma model, generalized AKNS equations, are discussed. By reducing the O(2,1) non-linear sigma model to its Hamiltonian form some new conservation laws are derived. A new hierarchy of non-linear evolution equations is proposed and shown to be generalized Hamiltonian equations with an infinite number of conservation laws. (author)

GPELab, a Matlab toolbox to solve Gross-Pitaevskii equations II: Dynamics and stochastic simulations

Science.gov (United States)

Antoine, Xavier; Duboscq, Romain

2015-08-01

GPELab is a free Matlab toolbox for modeling and numerically solving large classes of systems of Gross-Pitaevskii equations that arise in the physics of Bose-Einstein condensates. The aim of this second paper, which follows (Antoine and Duboscq, 2014), is to first present the various pseudospectral schemes available in GPELab for computing the deterministic and stochastic nonlinear dynamics of Gross-Pitaevskii equations (Antoine, et al., 2013). Next, the corresponding GPELab functions are explained in detail. Finally, some numerical examples are provided to show how the code works for the complex dynamics of BEC problems.

Stochastic dynamic modeling of regular and slow earthquakes

Science.gov (United States)

Aso, N.; Ando, R.; Ide, S.

2017-12-01

Both regular and slow earthquakes are slip phenomena on plate boundaries and are simulated by a (quasi-)dynamic modeling [Liu and Rice, 2005]. In these numerical simulations, spatial heterogeneity is usually considered not only for explaining real physical properties but also for evaluating the stability of the calculations or the sensitivity of the results on the condition. However, even though we discretize the model space with small grids, heterogeneity at smaller scales than the grid size is not considered in the models with deterministic governing equations. To evaluate the effect of heterogeneity at the smaller scales we need to consider stochastic interactions between slip and stress in a dynamic modeling. Tidal stress is known to trigger or affect both regular and slow earthquakes [Yabe et al., 2015; Ide et al., 2016], and such an external force with fluctuation can also be considered as a stochastic external force. A healing process of faults may also be stochastic, so we introduce stochastic friction law. In the present study, we propose a stochastic dynamic model to explain both regular and slow earthquakes. We solve mode III problem, which corresponds to the rupture propagation along the strike direction. We use BIEM (boundary integral equation method) scheme to simulate slip evolution, but we add stochastic perturbations in the governing equations, which is usually written in a deterministic manner. As the simplest type of perturbations, we adopt Gaussian deviations in the formulation of the slip-stress kernel, external force, and friction. By increasing the amplitude of perturbations of the slip-stress kernel, we reproduce complicated rupture process of regular earthquakes including unilateral and bilateral ruptures. By perturbing external force, we reproduce slow rupture propagation at a scale of km/day. The slow propagation generated by a combination of fast interaction at S-wave velocity is analogous to the kinetic theory of gasses: thermal

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.

Weak and Strong Order of Convergence of a Semidiscrete Scheme for the Stochastic Nonlinear Schrodinger Equation

International Nuclear Information System (INIS)

Bouard, Anne de; Debussche, Arnaud

2006-01-01

In this article we analyze the error of a semidiscrete scheme for the stochastic nonlinear Schrodinger equation with power nonlinearity. We consider supercritical or subcritical nonlinearity and the equation can be either focusing or defocusing. Allowing sufficient spatial regularity we prove that the numerical scheme has strong order 1/2 in general and order 1 if the noise is additive. Furthermore, we also prove that the weak order is always 1

Stochastic representation of a class of non-Markovian completely positive evolutions

International Nuclear Information System (INIS)

Budini, Adrian A.

2004-01-01

By modeling the interaction of an open quantum system with its environment through a natural generalization of the classical concept of continuous time random walk, we derive and characterize a class of non-Markovian master equations whose solution is a completely positive map. The structure of these master equations is associated with a random renewal process where each event consist in the application of a superoperator over a density matrix. Strong nonexponential decay arise by choosing different statistics of the renewal process. As examples we analyze the stochastic and averaged dynamics of simple systems that admit an analytical solution. The problem of positivity in quantum master equations induced by memory effects [S. M. Barnett and S. Stenholm, Phys. Rev. A 64, 033808 (2001)] is clarified in this context

Diffusion equations and hard collisions in multiple scattering of charged particles

International Nuclear Information System (INIS)

Papiez, Lech; Tulovsky, Vladimir

1998-01-01

The processes of angular-spatial evolution of multiple scattering of charged particles are described by the Lewis (special case of Boltzmann) integro-differential equation. The underlying stochastic process for this evolution is the compound Poisson process with transition densities satisfying the Lewis equation. In this paper we derive the Lewis equation from the compound Poisson process and show that the effective method of the solution of this equation can be based on the idea of decomposition of the compound Poisson process into processes of soft and hard collisions. Formulas for transition densities of soft and hard collision processes are provided in this paper together with the formula expressing the general solution of the Lewis equation in terms of those transition densities

Diffusion equations and hard collisions in multiple scattering of charged particles

Energy Technology Data Exchange (ETDEWEB)

Papiez, Lech [Department of Radiation Oncology, Indiana University, Indianapolis, IN (United States); Tulovsky, Vladimir [Department of Mathematics, St. John' s College, Staten Island, New York, NY (United States)

1998-09-01

The processes of angular-spatial evolution of multiple scattering of charged particles are described by the Lewis (special case of Boltzmann) integro-differential equation. The underlying stochastic process for this evolution is the compound Poisson process with transition densities satisfying the Lewis equation. In this paper we derive the Lewis equation from the compound Poisson process and show that the effective method of the solution of this equation can be based on the idea of decomposition of the compound Poisson process into processes of soft and hard collisions. Formulas for transition densities of soft and hard collision processes are provided in this paper together with the formula expressing the general solution of the Lewis equation in terms of those transition densities.

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

Invariant measure for the stochastic Navier-Stokes equations in unbounded 2D domains

Czech Academy of Sciences Publication Activity Database

Brzezniak, Z.; Motyl, E.; Ondreját, Martin

2017-01-01

Roč. 45, č. 5 (2017), s. 3145-3201 ISSN 0091-1798 R&D Projects: GA ČR(CZ) GA15-08819S Institutional support: RVO:67985556 Keywords : Invariant measure * bw-Feller semigroup * stochastic Navier–Stokes equation Subject RIV: BA - General Mathematics OBOR OECD: Statistics and probability Impact factor: 1.940, year: 2016 http://library.utia.cas.cz/separaty/2017/SI/ondrejat-0478383.pdf

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.

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.

QCD evolution equations for high energy partons in nuclear matter

CERN Document Server

Kinder-Geiger, Klaus; Geiger, Klaus; Mueller, Berndt

1994-01-01

We derive a generalized form of Altarelli-Parisi equations to decribe the time evolution of parton distributions in a nuclear medium. In the framework of the leading logarithmic approximation, we obtain a set of coupled integro- differential equations for the parton distribution functions and equations for the virtuality (``age'') distribution of partons. In addition to parton branching processes, we take into account fusion and scattering processes that are specific to QCD in medium. Detailed balance between gain and loss terms in the resulting evolution equations correctly accounts for both real and virtual contributions which yields a natural cancellation of infrared divergences.

Stochastic solutions to the Schrodinger equation for fermions

International Nuclear Information System (INIS)

Arnow, D.M.

1981-01-01

An exact stochastic method has been developed for generating the antisymmetric eigensolution of lowest index and its associated eigenvalue for the Schrodinger wave equation in 3N dimensions. The method is called the Green's function Monte Carlo method for fermions (FGFMC) because it is based on a Monte Carlo solution to the integral form of the Schrodinger equation (using Green's function) and because it is the fermion class of particles in physics which require antisymmetric solutions. The solution consists of two sets of 3N-dimensional points, [R/sub j/ + ] and [R/sub j/ - ], distributed by density functions psi + and psi - , whose difference, psi + -psi - , is proportional to the eigensolution, psi/sub F/. The FGFMC method is successfully applied to a one dimensional problem and a nine dimensional problem, the results of which are presented here. These results demonstrate that this method can be successfully applied to small physical problems on medium-scale computing machines. The key to this success was the transformation of the problem from exponential to linear cost as a function of accuracy. The strong dependence on dimensionality, however, currently results in an exponential cost as a function of problem size, and this, until overcome, imposes a severe barrier to calculations on large systems

From stochastic phase space evolution to Brownian motion in collective space

International Nuclear Information System (INIS)

Benhassine, B.; Farine, M.; Hernandez, E.S.; Idier, D.; Remaud, B.; Sebille, F.

1993-01-01

Within the framework of stochastic transport equations in phase space, the dynamics of fluctuations on collective variables in homogeneous fermion systems is studied. The transport coefficients are formally deduced in the relaxation time approximation and a general method to compute dynamically the dispersions of collective observables is proposed as a set of coupled equations. Independently, the general covariance matrix of phase space fluctuations and the dispersion on collective variables at equilibrium are derived. Detailed numerical applications show that dynamics of fluctuations can be extracted from noisy numerical simulations and that the leading parameter for collective fluctuations is the excitation energy whatever is its degree of thermalization. (authors). 16 refs., 12 figs

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 Systems Uncertainty Quantification and Propagation

CERN Document Server

Grigoriu, Mircea

2012-01-01

Uncertainty is an inherent feature of both properties of physical systems and the inputs to these systems that needs to be quantified for cost effective and reliable designs. The states of these systems satisfy equations with random entries, referred to as stochastic equations, so that they are random functions of time and/or space. The solution of stochastic equations poses notable technical difficulties that are frequently circumvented by heuristic assumptions at the expense of accuracy and rigor. The main objective of Stochastic Systems is to promoting the development of accurate and efficient methods for solving stochastic equations and to foster interactions between engineers, scientists, and mathematicians. To achieve these objectives Stochastic Systems presents: ·         A clear and brief review of essential concepts on probability theory, random functions, stochastic calculus, Monte Carlo simulation, and functional analysis   ·          Probabilistic models for random variables an...

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.

Sign reversals of the output autocorrelation function for the stochastic Bernoulli-Verhulst equation

Energy Technology Data Exchange (ETDEWEB)

Lumi, N., E-mail: Neeme.Lumi@tlu.ee; Mankin, R., E-mail: Romi.Mankin@tlu.ee [Institute of Mathematics and Natural Sciences, Tallinn University, 29 Narva Road, 10120 Tallinn (Estonia)

2015-10-28

We consider a stochastic Bernoulli-Verhulst equation as a model for population growth processes. The effect of fluctuating environment on the carrying capacity of a population is modeled as colored dichotomous noise. Relying on the composite master equation an explicit expression for the stationary autocorrelation function (ACF) of population sizes is found. On the basis of this expression a nonmonotonic decay of the ACF by increasing lag-time is shown. Moreover, in a certain regime of the noise parameters the ACF demonstrates anticorrelation as well as related sign reversals at some values of the lag-time. The conditions for the appearance of this highly unexpected effect are also discussed.

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.

Prolongation Structure of Semi-discrete Nonlinear Evolution Equations

International Nuclear Information System (INIS)

Bai Yongqiang; Wu Ke; Zhao Weizhong; Guo Hanying

2007-01-01

Based on noncommutative differential calculus, we present a theory of prolongation structure for semi-discrete nonlinear evolution equations. As an illustrative example, a semi-discrete model of the nonlinear Schroedinger equation is discussed in terms of this theory and the corresponding Lax pairs are also given.

Controlled Nonlinear Stochastic Delay Equations: Part II: Approximations and Pipe-Flow Representations

International Nuclear Information System (INIS)

Kushner, Harold J.

2012-01-01

This is the second part of a work dealing with key issues that have not been addressed in the modeling and numerical optimization of nonlinear stochastic delay systems. We consider new classes of models, such as those with nonlinear functions of several controls (such as products), each with is own delay, controlled random Poisson measure driving terms, admissions control with delayed retrials, and others. Part I was concerned with issues concerning the class of admissible controls and their approximations, since the classical definitions are inadequate for our models. This part is concerned with transportation equation representations and their approximations. Such representations of nonlinear stochastic delay models have been crucial in the development of numerical algorithms with much reduced memory and computational requirements. The representations for the new models are not obvious and are developed. They also provide a template for the adaptation of the Markov chain approximation numerical methods.

Nonlinear evolution equations having a physical meaning

International Nuclear Information System (INIS)

Nakach, R.

1976-06-01

The non stationary self-similar solutions of the nonlinear evolution equations which can be solved by the inverse scattering method are studied. It turns out, as shown by means of several examples, that when the L linear operator associated with these equations, is of second order and only then, the self-similar solutions can be expressed in terms of the various Painleve's transcendents [fr

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).

Critical spaces for quasilinear parabolic evolution equations and applications

Science.gov (United States)

Prüss, Jan; Simonett, Gieri; Wilke, Mathias

2018-02-01

We present a comprehensive theory of critical spaces for the broad class of quasilinear parabolic evolution equations. The approach is based on maximal Lp-regularity in time-weighted function spaces. It is shown that our notion of critical spaces coincides with the concept of scaling invariant spaces in case that the underlying partial differential equation enjoys a scaling invariance. Applications to the vorticity equations for the Navier-Stokes problem, convection-diffusion equations, the Nernst-Planck-Poisson equations in electro-chemistry, chemotaxis equations, the MHD equations, and some other well-known parabolic equations are given.

Soliton evolution and radiation loss for the Korteweg--de Vries equation

International Nuclear Information System (INIS)

Kath, W.L.; Smyth, N.F.

1995-01-01

The time-dependent behavior of solutions of the Korteweg--de Vries (KdV) equation for nonsoliton initial conditions is considered. While the exact solution of the KdV equation can in principle be obtained using the inverse scattering transform, in practice it can be extremely difficult to obtain information about a solution's transient evolution by this method. As an alternative, we present here an approximate method for investigating this transient evolution which is based upon the conservation laws associated with the KdV equation. Initial conditions which form one or two solitons are considered, and the resulting approximate evolution is found to be in good agreement with the numerical solution of the KdV equation. Justification for the approximations employed is also given by way of the linearized inverse scattering solution of the KdV equation. In addition, the final soliton state determined from the approximate equations agrees very well with the final state determined from the exact inverse scattering transform solution

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.

Hill functions for stochastic gene regulatory networks from master equations with split nodes and time-scale separation

Science.gov (United States)

Lipan, Ovidiu; Ferwerda, Cameron

2018-02-01

The deterministic Hill function depends only on the average values of molecule numbers. To account for the fluctuations in the molecule numbers, the argument of the Hill function needs to contain the means, the standard deviations, and the correlations. Here we present a method that allows for stochastic Hill functions to be constructed from the dynamical evolution of stochastic biocircuits with specific topologies. These stochastic Hill functions are presented in a closed analytical form so that they can be easily incorporated in models for large genetic regulatory networks. Using a repressive biocircuit as an example, we show by Monte Carlo simulations that the traditional deterministic Hill function inaccurately predicts time of repression by an order of two magnitudes. However, the stochastic Hill function was able to capture the fluctuations and thus accurately predicted the time of repression.

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.

A general nonlinear evolution equation for irreversible conservative approach to stable equilibrium

International Nuclear Information System (INIS)

Beretta, G.P.

1986-01-01

This paper addresses a mathematical problem relevant to the question of nonequilibrium and irreversibility, namely, that of ''designing'' a general evolution equation capable of describing irreversible but conservative relaxtion towards equilibrium. The objective is to present an interesting mathematical solution to this design problem, namely, a new nonlinear evolution equation that satisfies a set of very stringent relevant requirements. Three different frameworks are defined from which the new equation could be adopted, with entirely different interpretations. Some useful well-known mathematics involving Gram determinants are presented and a nonlinear evolution equation is given which meets the stringent design specifications

Semigroup methods for evolution equations on networks

CERN Document Server

Mugnolo, Delio

2014-01-01

This concise text is based on a series of lectures held only a few years ago and originally intended as an introduction to known results on linear hyperbolic and parabolic equations.  Yet the topic of differential equations on graphs, ramified spaces, and more general network-like objects has recently gained significant momentum and, well beyond the confines of mathematics, there is a lively interdisciplinary discourse on all aspects of so-called complex networks. Such network-like structures can be found in virtually all branches of science, engineering and the humanities, and future research thus calls for solid theoretical foundations.      This book is specifically devoted to the study of evolution equations – i.e., of time-dependent differential equations such as the heat equation, the wave equation, or the Schrödinger equation (quantum graphs) – bearing in mind that the majority of the literature in the last ten years on the subject of differential equations of graphs has been devoted to ellip...

Spatial stochasticity and non-continuum effects in gas flows

Energy Technology Data Exchange (ETDEWEB)

Dadzie, S. Kokou, E-mail: k.dadzie@glyndwr.ac.uk [Mechanical and Aeronautical Engineering, Glyndwr University, Mold Road, Wrexham LL11 2AW (United Kingdom); Reese, Jason M., E-mail: jason.reese@strath.ac.uk [Department of Mechanical and Aerospace Engineering, University of Strathclyde, Glasgow G1 1XJ (United Kingdom)

2012-02-06

We investigate the relationship between spatial stochasticity and non-continuum effects in gas flows. A kinetic model for a dilute gas is developed using strictly a stochastic molecular model reasoning, without primarily referring to either the Liouville or the Boltzmann equations for dilute gases. The kinetic equation, a stochastic version of the well-known deterministic Boltzmann equation for dilute gas, is then associated with a set of macroscopic equations for the case of a monatomic gas. Tests based on a heat conduction configuration and sound wave dispersion show that spatial stochasticity can explain some non-continuum effects seen in gases. -- Highlights: ► We investigate effects of molecular spatial stochasticity in non-continuum regime. ► Present a simplify spatial stochastic kinetic equation. ► Present a spatial stochastic macroscopic flow equations. ► Show effects of the new model on sound wave dispersion prediction. ► Show effects of the new approach in density profiles in a heat conduction.

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.

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.

Nonlocal higher order evolution equations

KAUST Repository

Rossi, Julio D.

2010-06-01

In this article, we study the asymptotic behaviour of solutions to the nonlocal operator ut(x, t)1/4(-1)n-1 (J*Id -1)n (u(x, t)), x ∈ ℝN, which is the nonlocal analogous to the higher order local evolution equation vt(-1)n-1(Δ)nv. We prove that the solutions of the nonlocal problem converge to the solution of the higher order problem with the right-hand side given by powers of the Laplacian when the kernel J is rescaled in an appropriate way. Moreover, we prove that solutions to both equations have the same asymptotic decay rate as t goes to infinity. © 2010 Taylor & Francis.

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...

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.

Symmetry Reduction and Cauchy Problems for a Class of Fourth-Order Evolution Equations

International Nuclear Information System (INIS)

Li Jina; Zhang Shunli

2008-01-01

We exploit higher-order conditional symmetry to reduce initial-value problems for evolution equations to Cauchy problems for systems of ordinary differential equations (ODEs). We classify a class of fourth-order evolution equations which admit certain higher-order generalized conditional symmetries (GCSs) and give some examples to show the main reduction procedure. These reductions cannot be derived within the framework of the standard Lie approach, which hints that the technique presented here is something essential for the dimensional reduction of evolution equations

The Liouville equation for flavour evolution of neutrinos and neutrino wave packets

Energy Technology Data Exchange (ETDEWEB)

Hansen, Rasmus Sloth Lundkvist; Smirnov, Alexei Yu., E-mail: rasmus@mpi-hd.mpg.de, E-mail: smirnov@mpi-hd.mpg.de [Max-Planck-Institut für Kernphysik, Saupfercheckweg 1, 69117 Heidelberg (Germany)

2016-12-01

We consider several aspects related to the form, derivation and applications of the Liouville equation (LE) for flavour evolution of neutrinos. To take into account the quantum nature of neutrinos we derive the evolution equation for the matrix of densities using wave packets instead of Wigner functions. The obtained equation differs from the standard LE by an additional term which is proportional to the difference of group velocities. We show that this term describes loss of the propagation coherence in the system. In absence of momentum changing collisions, the LE can be reduced to a single derivative equation over a trajectory coordinate. Additional time and spatial dependence may stem from initial (production) conditions. The transition from single neutrino evolution to the evolution of a neutrino gas is considered.

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

Finite difference evolution equations and quantum dynamical semigroups

International Nuclear Information System (INIS)

Ghirardi, G.C.; Weber, T.

1983-12-01

We consider the recently proposed [Bonifacio, Lett. Nuovo Cimento, 37, 481 (1983)] coarse grained description of time evolution for the density operator rho(t) through a finite difference equation with steps tau, and we prove that there exists a generator of the quantum dynamical semigroup type yielding an equation giving a continuous evolution coinciding at all time steps with the one induced by the coarse grained description. The map rho(0)→rho(t) derived in this way takes the standard form originally proposed by Lindblad [Comm. Math. Phys., 48, 119 (1976)], even when the map itself (and, therefore, the corresponding generator) is not bounded. (author)

Periodic feedback stabilization for linear periodic evolution equations

CERN Document Server

Wang, Gengsheng

2016-01-01

This book introduces a number of recent advances regarding periodic feedback stabilization for linear and time periodic evolution equations. First, it presents selected connections between linear quadratic optimal control theory and feedback stabilization theory for linear periodic evolution equations. Secondly, it identifies several criteria for the periodic feedback stabilization from the perspective of geometry, algebra and analyses respectively. Next, it describes several ways to design periodic feedback laws. Lastly, the book introduces readers to key methods for designing the control machines. Given its coverage and scope, it offers a helpful guide for graduate students and researchers in the areas of control theory and applied mathematics.

Spectral transform and solvability of nonlinear evolution equations

International Nuclear Information System (INIS)

Degasperis, A.

1979-01-01

These lectures deal with an exciting development of the last decade, namely the resolving method based on the spectral transform which can be considered as an extension of the Fourier analysis to nonlinear evolution equations. Since many important physical phenomena are modeled by nonlinear partial wave equations this method is certainly a major breakthrough in mathematical physics. We follow the approach, introduced by Calogero, which generalizes the usual Wronskian relations for solutions of a Sturm-Liouville problem. Its application to the multichannel Schroedinger problem will be the subject of these lectures. We will focus upon dynamical systems described at time t by a multicomponent field depending on one space coordinate only. After recalling the Fourier technique for linear evolution equations we introduce the spectral transform method taking the integral equations of potential scattering as an example. The second part contains all the basic functional relationships between the fields and their spectral transforms as derived from the Wronskian approach. In the third part we discuss a particular class of solutions of nonlinear evolution equations, solitons, which are considered by many physicists as a first step towards an elementary particle theory, because of their particle-like behaviour. The effect of the polarization time-dependence on the motion of the soliton is studied by means of the corresponding spectral transform, leading to new concepts such as the 'boomeron' and the 'trappon'. The rich dynamic structure is illustrated by a brief report on the main results of boomeron-boomeron and boomeron-trappon collisions. In the final section we discuss further results concerning important properties of the solutions of basic nonlinear equations. We introduce the Baecklund transform for the special case of scalar fields and demonstrate how it can be used to generate multisoliton solutions and how the conservation laws are obtained. (HJ)

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.

Stochastic Evolution Dynamic of the Rock-Scissors-Paper Game Based on a Quasi Birth and Death Process.

Science.gov (United States)

Yu, Qian; Fang, Debin; Zhang, Xiaoling; Jin, Chen; Ren, Qiyu

2016-06-27

Stochasticity plays an important role in the evolutionary dynamic of cyclic dominance within a finite population. To investigate the stochastic evolution process of the behaviour of bounded rational individuals, we model the Rock-Scissors-Paper (RSP) game as a finite, state dependent Quasi Birth and Death (QBD) process. We assume that bounded rational players can adjust their strategies by imitating the successful strategy according to the payoffs of the last round of the game, and then analyse the limiting distribution of the QBD process for the game stochastic evolutionary dynamic. The numerical experiments results are exhibited as pseudo colour ternary heat maps. Comparisons of these diagrams shows that the convergence property of long run equilibrium of the RSP game in populations depends on population size and the parameter of the payoff matrix and noise factor. The long run equilibrium is asymptotically stable, neutrally stable and unstable respectively according to the normalised parameters in the payoff matrix. Moreover, the results show that the distribution probability becomes more concentrated with a larger population size. This indicates that increasing the population size also increases the convergence speed of the stochastic evolution process while simultaneously reducing the influence of the noise factor.

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.

Stochastic Field evolution of disoriented chiral condensates

International Nuclear Information System (INIS)

Bettencourt, Luis M.A.

2003-01-01

I present a summary of recent work [1] where we describe the time-evolution of a region of disoriented chiral condensate via Langevin field equations for the linear σ model. We analyze the model in equilibrium, paying attention to subtracting ultraviolet divergent classical terms and replacing them by their finite quantum counter-parts. We use results from lattice gauge theory and chiral perturbation theory to fix nonuniversal constants. The result is a ultraviolet cutoff independent theory that reproduces quantitatively the expected equilibrium behavior of pion and σ quantum fields. We also estimate the viscosity η(T), which controls the dynamical timescale in the Langevin equation, so that the near equilibrium dynamical response agrees with theoretical expectations

Parametrization of open systems with effective quadratic hamiltonians plus stochastic force

International Nuclear Information System (INIS)

Hernandez, E.S.; Mizrahi, S.S.

1981-12-01

The evolution generated by general dissipative Hamiltonians is analyzed when a stochastic force is included. A mapping technique allows to easily write the equations of motion for the observables of interest. A general dissipativity condition is extracted, whose fullfilment guarantees that thermal equilibrium is reached as the final stage of the evolution. Several existing frictional Hamiltonians are examined and it is seen that the correlation of the fluctuating force is essential to the destruction of a constant of motion inherent to pure quantal behaviour. (Author) [pt