BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jinqiao Duan

Jinqiao Duan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.DSmath.PRmath.APmath-phmath.MPnlin.CDMachine Learningphysics.ao-phphysics.data-anmath.NAcond-mat.stat-mechmath.STphysics.flu-dynStatistics TheoryNumerical Analysisphysics.comp-phphysics.geo-phBiological Physicschao-dynInformation Theorymath.FAmath.ITmath.OCMethodologyMolecular Networksnlin.AOnlin.PSOther Quantitative BiologyPopulations and Evolution

Catalog footprint

What is connected

82works
29topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2024arXiv

Fourier neural operator based fluid-structure interaction for predicting the vesicle dynamics

Solving complex fluid-structure interaction (FSI) problems, characterized by nonlinear partial differential equations, is crucial in various scientific and engineering applications. Traditional computational fluid dynamics (CFD) solvers are insufficient to meet the growing requirements for large-scale and long-period simulations. Fortunately, the rapid advancement in neural networks, especially neural operator learning mappings between function spaces, has introduced novel approaches to tackle these challenges via data-driven modeling. In this paper, we propose a Fourier neural operator-based fluid-structure interaction solver (FNO-based FSI solver) for efficient simulation of FSI problems, where the solid solver based on the finite difference method is seamlessly integrated with the Fourier neural operator to predict incompressible flow using the immersed boundary method. We analyze the performance of the FNO-based FSI solver in the following three situations: training data with or without the steady state, training method with one-step label or multi-step labels, and prediction in interpolation or extrapolation. We find that the best performance for interpolation is achieved by training the operator with multi-step labels using steady-state data. Finally, we train the FNO-based FSI solver using this optimal training method and apply it to vesicle dynamics. The results show that the FNO-based FSI solver is capable of capturing the variations in the fluid and the vesicle.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Learning effective dynamics from data-driven stochastic systems

Multiscale stochastic dynamical systems have been widely adopted to a variety of scientific and engineering problems due to their capability of depicting complex phenomena in many real world applications. This work is devoted to investigating the effective dynamics for slow-fast stochastic dynamical systems. Given observation data on a short-term period satisfying some unknown slow-fast stochastic systems, we propose a novel algorithm including a neural network called Auto-SDE to learn invariant slow manifold. Our approach captures the evolutionary nature of a series of time-dependent autoencoder neural networks with the loss constructed from a discretized stochastic differential equation. Our algorithm is also validated to be accurate, stable and effective through numerical experiments under various evaluation metrics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An end-to-end deep learning approach for extracting stochastic dynamical systems with $α$-stable Lévy noise

Recently, extracting data-driven governing laws of dynamical systems through deep learning frameworks has gained a lot of attention in various fields. Moreover, a growing amount of research work tends to transfer deterministic dynamical systems to stochastic dynamical systems, especially those driven by non-Gaussian multiplicative noise. However, lots of log-likelihood based algorithms that work well for Gaussian cases cannot be directly extended to non-Gaussian scenarios which could have high error and low convergence issues. In this work, we overcome some of these challenges and identify stochastic dynamical systems driven by $α$-stable Lévy noise from only random pairwise data. Our innovations include: (1) designing a deep learning approach to learn both drift and diffusion coefficients for Lévy induced noise with $α$ across all values, (2) learning complex multiplicative noise without restrictions on small noise intensity, (3) proposing an end-to-end complete framework for stochastic systems identification under a general input data assumption, that is, $α$-stable random variable. Finally, numerical experiments and comparisons with the non-local Kramers-Moyal formulas with moment generating function confirm the effectiveness of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Onsager-Machlup approach to the most probable transition pathway for a genetic regulatory network

We investigate a quantitative network of gene expression dynamics describing the competence development in Bacillus subtilis. First, we introduce an Onsager-Machlup approach to quantify the most probable transition pathway for both excitable and bistable dynamics. Then, we apply a machine learning method to calculate the most probable transition pathway via the Euler-Lagrangian equation. Finally, we analyze how the noise intensity affects the transition phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data-driven method to learn the most probable transition pathway and stochastic differential equations

Transition phenomena between metastable states play an important role in complex systems due to noisy fluctuations. In this paper, the physics informed neural networks (PINNs) are presented to compute the most probable transition pathway. It is shown that the expected loss is bounded by the empirical loss. And the convergence result for the empirical loss is obtained. Then, a sampling method of rare events is presented to simulate the transition path by the Markovian bridge process. And we investigate the inverse problem to extract the stochastic differential equation from the most probable transition pathway data and the Markovian bridge process data, respectively. Finally, several numerical experiments are presented to verify the effectiveness of our methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker-Planck equations). We analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker-Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Detection of bistable structures via the Conley index and applications to biological systems

Bistability is a ubiquitous phenomenon in life sciences. In this paper, two kinds of bistable structures in dynamical systems are studied: One is two one-point attractors, another is a one-point attractor accompanied by a cycle attractor. By the Conley index theory, we prove that there exist other isolated invariant sets besides the two attractors, and also obtain the possible components and their configuration. Moreover, we find that there is always a separatrix or cycle separatrix, which separates the two attractors. Finally, the biological meanings and implications of these structures are given and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extracting Governing Laws from Sample Path Data of Non-Gaussian Stochastic Dynamical Systems

Advances in data science are leading to new progresses in the analysis and understanding of complex dynamics for systems with experimental and observational data. With numerous physical phenomena exhibiting bursting, flights, hopping, and intermittent features, stochastic differential equations with non-Gaussian Lévy noise are suitable to model these systems. Thus it is desirable and essential to infer such equations from available data to reasonably predict dynamical behaviors. In this work, we consider a data-driven method to extract stochastic dynamical systems with non-Gaussian asymmetric (rather than the symmetric) Lévy process, as well as Gaussian Brownian motion. We establish a theoretical framework and design a numerical algorithm to compute the asymmetric Lévy jump measure, drift and diffusion (i.e., nonlocal Kramers-Moyal formulas), hence obtaining the stochastic governing law, from noisy data. Numerical experiments on several prototypical examples confirm the efficacy and accuracy of this method. This method will become an effective tool in discovering the governing laws from available data sets and in understanding the mechanisms underlying complex random phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pointwise convergence problem of Ostrovsky equation with rough data and random data

In this paper, we consider the pointwise convergence problem of free Ostrovsky equation with rough data and random data. Firstly, we show the almost everywhere pointwise convergence of free Ostrovsky equation in $H^{s}(\mathbb{R})$ with $s\geq \frac{1}{4}$ with rough data. Secondly, we present counterexamples showing that the maximal function estimate related to the free Ostrovsky equation can fail if $s<\frac{1}{4}$. Finally, for every $x\in \mathbb{R}$, we show the almost surely pointwise convergence of free Ostrovsky equation in $L^{2}(\mathbb{R})$ with random data. The main tools are the density theorem, high-low frequency idea, Wiener decomposition and Lemmas 2.1-2.6 as well as the probabilistic estimates of some random series which are just Lemmas 3.2-3.4 in this paper. The main difficulty is that zero is the singular point of the phase functions of free Ostrovsky equation. We use high-low frequency idea to conquer the difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantifying model uncertainty for the observed non-Gaussian data by the Hellinger distance

Mathematical models for complex systems under random fluctuations often certain uncertain parameters. However, quantifying model uncertainty for a stochastic differential equation with an $α$-stable Lévy process is still lacking. Here, we propose an approach to infer all the uncertain non-Gaussian parameters and other system parameters by minimizing the Hellinger distance over the parameter space. The Hellinger distance measures the similarity between an empirical probability density of non-Gaussian observations and a solution (as a probability density) of the associated nonlocal Fokker-Planck equation. Numerical experiments verify that our method is feasible for estimating single and multiple parameters. Meanwhile, we find an optimal estimation interval of the estimated parameters. This method is beneficial for extracting governing dynamical system models under non-Gaussian fluctuations, as in the study of abrupt climate changes in the Dansgaard-Oeschger events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Data-Driven Approach for Discovering Stochastic Dynamical Systems with Non-Gaussian Levy Noise

With the rapid increase of valuable observational, experimental and simulating data for complex systems, great efforts are being devoted to discovering governing laws underlying the evolution of these systems. However, the existing techniques are limited to extract governing laws from data as either deterministic differential equations or stochastic differential equations with Gaussian noise. In the present work, we develop a new data-driven approach to extract stochastic dynamical systems with non-Gaussian symmetric Lévy noise, as well as Gaussian noise. First, we establish a feasible theoretical framework, by expressing the drift coefficient, diffusion coefficient and jump measure (i.e., anomalous diffusion) for the underlying stochastic dynamical system in terms of sample paths data. We then design a numerical algorithm to compute the drift, diffusion coefficient and jump measure, and thus extract a governing stochastic differential equation with Gaussian and non-Gaussian noise. Finally, we demonstrate the efficacy and accuracy of our approach by applying to several prototypical one-, two- and three-dimensional systems. This new approach will become a tool in discovering governing dynamical laws from noisy data sets, from observing or simulating complex phenomena, such as rare events triggered by random fluctuations with heavy as well as light tail statistical features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Logistic-Harvest Model with Allee Effect under Multiplicative Noise

This work is devoted to the study of a stochastic logistic growth model with and without the Allee effect. Such a model describes the evolution of a population under environmental stochastic fluctuations and is in the form of a stochastic differential equation driven by multiplicative Gaussian noise. With the help of the associated Fokker-Planck equation, we analyze the population extinction probability and the probability of reaching a large population size before reaching a small one. We further study the impact of the harvest rate, noise intensity, and the Allee effect on population evolution. The analysis and numerical experiments show that if the noise intensity and harvest rate are small, the population grows exponentially, and upon reaching the carrying capacity, the population size fluctuates around it. In the stochastic logistic-harvest model without the Allee effect, when noise intensity becomes small (or goes to zero), the stationary probability density becomes more acute and its maximum point approaches one. However, for large noise intensity and harvest rate, the population size fluctuates wildly and does not grow exponentially to the carrying capacity. So as far as biological meanings are concerned, we must catch at small values of noise intensity and harvest rate. Finally, we discuss the biological implications of our results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Machine Learning Framework for Computing the Most Probable Paths of Stochastic Dynamical Systems

The emergence of transition phenomena between metastable states induced by noise plays a fundamental role in a broad range of nonlinear systems. The computation of the most probable paths is a key issue to understand the mechanism of transition behaviors. Shooting method is a common technique for this purpose to solve the Euler-Lagrange equation for the associated action functional, while losing its efficacy in high-dimensional systems. In the present work, we develop a machine learning framework to compute the most probable paths in the sense of Onsager-Machlup action functional theory. Specifically, we reformulate the boundary value problem of Hamiltonian system and design a neural network to remedy the shortcomings of shooting method. The successful applications of our algorithms to several prototypical examples demonstrate its efficacy and accuracy for stochastic systems with both (Gaussian) Brownian noise and (non-Gaussian) Lévy noise. This novel approach is effective in exploring the internal mechanisms of rare events triggered by random fluctuations in various scientific fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An averaging principle for fractional stochastic differential equations with Lévy noise

This paper is devoted to the study of an averaging principle for fractional stochastic differential equations in Rnwith Lévy motion, using an integral transform method. We obtain a time-averaged equation under suitable assumptions. Furthermore, we show that the solutions of averaged equation approach the solutions of the original equation. Our results in this paper provide better understanding for effective approximation of fractional dynamical systems with non-Gaussian Lévy noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of multiscale methods for stochastic dynamical systems driven by $α$-stable processes

In this paper, we first analyze the strong and weak convergence of projective integration methods for multiscale stochastic dynamical systems driven by $α$-stable processes, which are used to estimate the effect that the fast components have on slow ones. Then we obtain the $p$th moment error bounds between the solution of slow component produced by projective integration method and the solution of effective system with $p \in \left(1, α\right)$. Finally, we corroborate our analytical results through a specific numerical example.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discovering transition phenomena from data of stochastic dynamical systems with Lévy noise

It is a challenging issue to analyze complex dynamics from observed and simulated data. An advantage of extracting dynamic behaviors from data is that this approach enables the investigation of nonlinear phenomena whose mathematical models are unavailable. The purpose of this present work is to extract information about transition phenomena (e.g., mean exit time and escape probability), from data of stochastic differential equations with non-Gaussian Lévy noise. As a tool in describing dynamical systems, the Koopman semigroup transforms a nonlinear system into a linear system, but at the cost of elevating a finite dimensional problem into an infinite dimensional one. In spite of this, using the relation between the stochastic Koopman semigroup and the infinitesimal generator of a stochastic differential equation, we learn the mean exit time and escape probability from data. Specifically, we first obtain a finite dimensional approximation of the infinitesimal generator by an extended dynamic mode decomposition algorithm. Then we identify the drift coefficient, diffusion coefficient and anomalous diffusion coefficient for the stochastic differential equation. Finally, we compute the mean exit time and escape probability by finite difference discretization of the associated nonlocal partial differential equations. This approach is applicable in extracting transition information from data of stochastic differential equations with either (Gaussian) Brownian motion or (non-Gaussian) Lévy motion. We present one - and two-dimensional examples to demonstrate the effectiveness of our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective reduction for a nonlocal Zakai stochastic partial differential equation in data assimilation

We study the effective reduction for a nonlocal stochastic partial differential equation with oscillating coefficients. The nonlocal operator in this stochastic partial differential equation is the generator of non-Gaussian Lévy processes, with either \textbf{integrable} or \textbf{non-integrable} jump kernels. We examine the limiting behavior of this equation as a scaling parameter tends to zero, and derive a reduced (local or nonlocal) effective equation. In particular, this work leads to an effective reduction for a data assimilation system with Lévy noise, by examining the corresponding nonlocal Zakai stochastic partial differential equation. We show that the probability density for the reduced data assimilation system approximates that for the original system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective Wave Factorization for a Stochastic Schrödinger Equation

We study the homogenization of a stochastic Schrödinger equation with a large periodic potential in solid state physics. Denoting by $\varepsilon$ the period, the potential is scaled as $\varepsilon^{-2}$. Under a generic assumption on the spectral properties of the associated cell problem, we prove that the solution can be approximately factorized as the product of a fast oscillating cell eigenfunction and of a slowly varying solution of an effective equation. Our method is based on two-scale convergence and Bloch waves theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identifying stochastic governing equations from data of the most probable transition trajectories

Extracting governing stochastic differential equation models from elusive data is crucial to understand and forecast dynamics for complex systems. We devise a method to extract the drift term and estimate the diffusion coefficient of a governing stochastic dynamical system, from its time-series data of the most probable transition trajectory. By the Onsager-Machlup theory, the most probable transition trajectory satisfies the corresponding Euler-Lagrange equation, which is a second order deterministic ordinary differential equation involving the drift term and diffusion coefficient. We first estimate the coefficients of the Euler-Lagrange equation based on the data of the most probable trajectory, and then we calculate the drift and diffusion coefficients of the governing stochastic dynamical system. These two steps involve sparse regression and optimization. Finally, we illustrate our method with an example and some discussions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linear Response Theory for Nonlinear Stochastic Differential Equations with $α$-stable Lévy Noises

We consider a nonlinear stochastic differential equation driven by an $α$-stable Lévy process ($1<α<2$). We first obtain some regularity results for the probability density of its invariant measure via establishing the a priori estimate of the corresponding stationary Fokker-Planck equation. Then by the a priori estimate of Kolmogorov backward equations and the perturbation property of Markov semigroup, we derive the response function and generalize the famous linear response theory in nonequilibrium statistical mechanics to non-Gaussian stochastic dynamic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximal Likely Phase Lines for a Reduced Ice Growth Model

We study the impact of Brownian noise on transitions between metastable equilibrium states in a stochastic ice sheet model. Two methods to accomplish different objectives are employed. The maximal likely trajectory by maximizing the probability density function and numerically solving the Fokker-Planck equation shows how the system will evolve over time. We have especially studied the maximal likely trajectories starting near the ice-free metastable state, and examined whether they evolve to or near the ice-covered metastable state for certain parameters, in order to gain insights into how the ice sheet formed. Furthermore, for the transition from ice-covered metastable state to the ice-free metastable state, we study the most probable path for various noise parameters via the Onsager-Machlup least action principle. This enables us to predict and visualize the melting process of the ice sheet if such a rare event ever does take place.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum Likelihood Estimation of Stochastic Differential Equations with Random Effects Driven by Fractional Brownian Motion

Stochastic differential equations and stochastic dynamics are good models to describe stochastic phenomena in real world. In this paper, we study N independent stochastic processes Xi(t) with real entries and the processes are determined by the stochastic differential equations with drift term relying on some random effects. We obtain the Girsanov-type formula of the stochastic differential equation driven by Fractional Brownian Motion through kernel transformation. Under some assumptions of the random effect, we estimate the parameter estimators by the maximum likelihood estimation and give some numerical simulations for the discrete observations. Results show that for the different H, the parameter estimator is closer to the true value as the amount of data increases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Most Probable Dynamics of Stochastic Dynamical Systems with Exponentially Light Jump Fluctuations

The emergence of the exit events from a bounded domain containing a stable fixed point induced by non-Gaussian Lévy fluctuations plays a pivotal role in practical physical systems. In the limit of weak noise, we develop a Hamiltonian formalism under the Lévy fluctuations with exponentially light jumps for one- and two-dimensional stochastic dynamical systems. This formalism is based on a recently proved large deviation principle for dynamical systems under non-Gaussian Lévy perturbations. We demonstrate how to compute the most probable exit path and the quasi-potential by several examples. Meanwhile, we explore the impacts of the jump measure on the quasi-potential quantitatively and on the most probable exit path qualitatively. Results show that the quasi-potential can be well estimated by an approximate analytical expression. Moreover, we discover that although the most probable exit paths are analogous to the Gaussian case for the isotropic noise, the anisotropic noise leads to significant changes of the structure of the exit paths. These findings shed light on the underlying qualitative mechanism and quantitative feature of the exit phenomenon induced by non-Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical analysis and applications of Fokker-Planck equations for stochastic dynamical systems with multiplicative $α$-stable noises

The Fokker-Planck equations (FPEs) for stochastic systems driven by additive symmetric $α$-stable noises may not adequately describe the time evolution for the probability densities of solution paths in some practical applications, such as hydrodynamical systems, porous media, and composite materials. As a continuation of previous works on additive case, the FPEs for stochastic dynamical systems with multiplicative symmetric $α$-stable noises are derived by the adjoint operator method, which satisfy the nonlocal partial differential equations. A finite difference method for solving the nonlocal Fokker-Planck equation (FPE) is constructed, which is shown to satisfy the discrete maximum principle and to be convergent. Moreover, an example is given to illustrate this method. For asymmetric case, general finite difference schemes are proposed, and some analyses of the corresponding numerical schemes are given. Furthermore, the corresponding result is successfully applied to the nonlinear filtering problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recurrent Solutions of a Nonautonomous Modified Swift-Hohenberg Equation

We consider recurrent solutions of the nonautonomous modified Swift-Hohenberg equation $$u_t+Δ^2u+2Δu+au+b|\nabla u|^2+u^3=g(t,x).$$ We employ Conley index theory to show that, if the forcing $g:\mathbb{R}\rightarrow L^2(Ω)$ is a recurrent function, then there are at least two recurrent solutions in $H_0^2(Ω)$ under appropriate assumptions on the parameters $a$, $b$ and $g$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker-Planck Equation and Physics-informed Neural Networks

The Fokker-Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants, hence requiring the development of a data-driven modeling approach. When the data available is directly on the PDF, then there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles' positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback-Leibler divergence to connect the stochastic samples with the FP equation, to simultaneously learn the equation and infer the multi-dimensional PDF at all times. In particular, we consider two types of inverse problems, type I where the FP equation is known but the initial PDF is unknown, and type II in which, in addition to unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Levy noise or a combination of both. We demonstrate the new PINN framework in detail in the one-dimensional case (1D) but we also provide results for up to 5D demonstrating that we can infer both the FP equation and} dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symplectic Euler scheme for Hamiltonian stochastic differential equations driven by Levy noise

This paper proposes a general symplectic Euler scheme for a class of Hamiltonian stochastic differential equations driven by L$\acute{e}$vy noise in the sense of Marcus form. The convergence of the symplectic Euler scheme for this Hamiltonian stochastic differential equations is investigated. Realizable numerical implementation of this scheme is also provided in details. Numerical experiments are presented to illustrate the effectiveness and superiority of the proposed method by the simulations of its orbits, symplectic structure and Hamlitonian.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The role of slow manifolds in parameter estimation for a multiscale stochastic system with $α$-stable Lévy noise

This work is about parameter estimation for a fast-slow stochastic system with non-Gaussian $α$-stable Lévy noise. When the observations are only available for slow components, a system parameter is estimated and the accuracy for this estimation is quantified by $p$-moment with $p\in(1, α)$, with the help of a reduced system through random slow manifold approximation. This method provides an advantage in computational complexity and cost, due to the dimension reduction in stochastic systems. To numerically illustrate this method, and to corroborate that the parameter estimator based on the reduced slow system is a good approximation for the true parameter value of the original system, a prototypical example is present.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The tipping times in an Arctic sea ice system under influence of extreme events

In light of the rapid recent retreat of Arctic sea ice, the extreme weather events triggering the variability in Arctic ice cover has drawn increasing attention. A non-Gaussian $α$-stable Lévy process is thought to be an appropriate model to describe such extreme event. The maximal likely trajectory, based on the nonlocal Fokker-Planck equation, is applied to a nonautonomous Arctic sea ice system under $α$-stable Lévy noise. Two types of tipping times, the early-warning tipping time and the disaster-happening tipping time, are used to predict the critical time for the maximal likely transition from a perennially ice-covered state to a seasonally ice-free one, and from a seasonally ice-free state to a perennially ice-free one, respectively. We find that the increased intensity of extreme events results in shorter warning time for sea ice melting, and that an enhanced greenhouse effect will intensify this influence, making the arrival of warning time significantly earlier. Meanwhile, for the enhanced greenhouse effect, we discover that increased intensity and frequency of extreme events will advance the disaster-happening tipping time, in which an ice-free state is maintained throughout the year in the Arctic Ocean. Finally, we identify values of Lévy index $α$ and noise intensity $ε$ in $αε$-space that can trigger a transition between the Arctic sea ice state. These results provide an effective theoretical framework for studying Arctic sea ice variations under the influence of extreme events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Onsager-Machlup Function as Lagrangian for the Most Probable Path of a Jump-diffusion Process

This work is devoted to deriving the Onsager-Machlup function for a class of stochastic dynamical systems under (non-Gaussian) Levy noise as well as (Gaussian) Brownian noise, and examining the corresponding most probable paths. This Onsager-Machlup function is the Lagrangian giving the most probable path connecting metastable states for jump-diffusion processes. This is done by applying the Girsanov transformation for measures induced by jump-diffusion processes. Moreover, we have found this Lagrangian function is consistent with the result in the special case of diffusion processes. Finally, we apply this new Onsager-Machlup function to investigate dynamical behaviors analytically and numerically in several examples. These include the transitions from one metastable state to another metastable state in a double-well system, with numerical experiments illustrating most probable transition paths for various noise parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Derivation of Fokker-Planck equations for stochastic dynamical systems under excitation of multiplicative non-Gaussian white noise

Fokker-Planck equations describe time evolution of probability densities of stochastic dynamical systems and play an important role in quantifying propagation and evolution of uncertainty. Although Fokker-Planck equations can be written explicitly for nonlinear dynamical systems excited by Gaussian white noise, they are not available in general for nonlinear dynamical systems excited by multiplicative non-Gaussian white noise. Marcus stochastic differential equations are often appropriate models in engineering and physics for stochastic dynamical systems excited by non-Gaussian white noise. In this paper, we derive explicit forms of Fokker-Planck equations for one dimensional systems modeled by Marcus stochastic differential equations under multiplicative non-Gaussian white noise. As examples to illustrate the theoretical results, the derived formula is used to obtain Fokker-Plank equations for nonlinear dynamical systems under excitation of (i) $α$-stable white noise; (ii) combined Gaussian and Poisson white noise, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A regularity result for the nonlocal Fokker-Planck equation with Ornstein-Uhlenbeck drift

Despite there are numerous theoretical studies of stochastic differential equations with a symmetric $α$-stable Lévy noise, very few regularity results exist in the case of $0<α\leq1$. In this paper, we study the fractional Fokker-Planck equation with Ornstein-Uhlenbeck drift, and prove that there exists a unique solution, which is $C^\infty$ in space for $t>0$ when $α\in (0, 2]$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical Inference for Transitions in Stochastic Systems with $α-$stable Lévy Noise

A goal of data assimilation is to infer stochastic dynamical behaviors with available observations. We consider transition phenomena between metastable states for a stochastic system with (non-Gaussian) $α-$stable Lévy noise. With either discrete time or continuous time observations, we infer such transitions by computing the corresponding nonlocal Zakai equation (and its discrete time counterpart) and examining the most probable orbits for the state system. Examples are presented to demonstrate this approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Approximation of the inertial manifold for a nonlocal dynamical system

We consider inertial manifolds and their approximation for a class of partial differential equations with a nonlocal Laplacian operator $-(-Δ)^{\fracα{2}}$, with $0<α<2$. The nonlocal or fractional Laplacian operator represents an anomalous diffusion effect. We first establish the existence of an inertial manifold and highlight the influence of the parameter $α$. Then we approximate the inertial manifold when a small normal diffusion $\varepsilon Δ$ (with $\varepsilon \in (0, 1)$) enters the system, and obtain the estimate for the Hausdorff semi-distance between the inertial manifolds with and without normal diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymptotic methods for stochastic dynamical systems with small non-Gaussian Lévy noise

The goal of the paper is to analytically examine escape probabilities for dynamical systems driven by symmetric $α$-stable Lévy motions. Since escape probabilities are solutions of a type of integro-differential equations (i.e., differential equations with nonlocal interactions), asymptotic methods are offered to solve these equations to obtain escape probabilities when noises are sufficiently small. Three examples are presented to illustrate the asymptotic methods, and asymptotic escape probability is compared with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Impacts of Noise on a Class of Partial Differential Equations

This paper is concerned with effects of noise on the solutions of partial differential equations. We first provide a sufficient condition to ensure the existence of a unique positive solution for a class of stochastic parabolic equations. Then, we prove that noise could induce singularities (finite time blow up of solutions). Finally, we show that a stochastic Allen-Cahn equation does not have finite time singularities and the unique solution exists globally.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Martingale and Weak Solutions for a Stochastic Nonlocal Burgers Equation on Bounded Intervals

This work is about the existence of martingale solutions and weak solutions for a stochastic nonlocal Burgers equation on bounded intervals. The existence of a martingale solution is shown by using a Galerkin approximation, Prokhorov's theorem and Skorokhod's embedding theorem. The same Galerkin approximation also leads to the existence of weak solution for the corresponding deterministic nonlocal Burgers equation on a bounded domain.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stationary Measures for Stochastic Differential Equations with Jumps

In the paper, stationary measures of stochastic differential equations with jumps are considered. Under some general conditions, existence of stationary measures is proved through Markov measures and Lyapunov functions. Moreover, for two special cases, stationary measures are given by solutions of Fokker-Planck equations and long time limits for the distributions of system states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A parameter estimation method based on random slow manifolds

A parameter estimation method is devised for a slow-fast stochastic dynamical system, where often only the slow component is observable. By using the observations only on the slow component, the system parameters are estimated by working on the slow system on the random slow manifold. This offers a benefit of dimension reduction in quantifying parameters in stochastic dynamical systems. An example is presented to illustrate this method, and verify that the parameter estimator based on the lower dimensional, reduced slow system is a good approximation of the parameter estimator for original slow-fast stochastic dynamical system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Approximation of Random Slow Manifolds and Settling of Inertial Particles under Uncertainty

A method is provided for approximating random slow manifolds of a class of slow-fast stochastic dynamical systems. Thus approximate, low dimensional, reduced slow systems are obtained analytically in the case of sufficiently large time scale separation. To illustrate this dimension reduction procedure, the impact of random environmental fluctuations on the settling motion of inertial particles in a cellular flow field is examined. It is found that noise delays settling for some particles but enhances settling for others. A deterministic stable manifold is an agent to facilitate this phenomenon. Overall, noise appears to delay the settling in an averaged sense.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Center manifolds for infinite dimensional random dynamical systems

Stochastic center manifolds theory are crucial in modelling the dynamical behavior of complex systems under stochastic influences. A multiplicative ergodic theorem on Hilbert space is proved to be satisfied to the exponential trichotomy condition. Then the existence of stochastic center manifolds for infinite dimensional random dynamical systems is shown under the assumption of exponential trichotomy. The theory provides a support for the discretisations of nonlinear stochastic partial differential equations with space-time white noise.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fokker-Planck Equations for Stochastic Dynamical Systems with Symmetric Lévy Motions

The Fokker-Planck equations for stochastic dynamical systems, with non-Gaussian $α-$stable symmetric Lévy motions, have a nonlocal or fractional Laplacian term. This nonlocality is the manifestation of the effect of non-Gaussian fluctuations. Taking advantage of the Toeplitz matrix structure of the time-space discretization, a fast and accurate numerical algorithm is proposed to simulate the nonlocal Fokker-Planck equations, under either absorbing or natural conditions. The scheme is shown to satisfy a discrete maximum principle and to be convergent. It is validated against a known exact solution and the numerical solutions obtained by using other methods. The numerical results for two prototypical stochastic systems, the Ornstein-Uhlenbeck system and the double-well system are shown.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global solutions for a nonlocal Ginzberg-Landau equation and a nonlocal Fokker-Plank equation

This work is devoted to the study of a nonlocal Ginzberg-Landau equation by the semigroup method and a nonlocal Fokker-Plank equation by the viscosity vanishing method. For the nonlocal Ginzberg-Landau equation, there exists a unique global solution in the set $C^0(\mathbb{R}^+,\,H_0^{\fracα{2}}(D))\cap L_{loc}(\mathbb{R}^+,\,H_0^α(D))$, for $α\in (0,\,2)$. For the nonlocal Fokker-Plank equation, the regularity of the solution is weaker than that of the nonlocal Ginzberg-Landau equation due to the drift term.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantifying model uncertainty in non-Gaussian dynamical systems with observations on mean exit time or escape probability

Complex systems are sometimes subject to non Gaussian alpha stable Levy fluctuations. A new method is devised to estimate this uncertain parameter and other system parameters, using observations on either mean exit time or escape probability for the system evolution. It is based on solving an inverse problem for a deterministic, non-local partial differential equation via numerical optimization. The existing methods for estimating parameters require observations on system state sample paths for long time periods or probability densities at large spatial ranges. The method proposed here, instead, requires observations on mean exit time or escape probability only for an arbitrarily small spatial domain. This new method is beneficial to systems for which mean exit time or escape probability is feasible to observe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Slow foliation of a slow-fast stochastic evolutionary system

This work is concerned with the dynamics of a slow-fast stochastic evolutionary system quantified with a scale parameter. An invariant foliation decomposes the state space into geometric regions of different dynamical regimes, and thus helps understand dynamics. A slow invariant foliation is established for this system. It is shown that the slow foliation converges to a critical foliation (i.e., the scale parameter is zero) in probability distribution, as the scale parameter tends to zero. The approximation of slow foliation is also constructed with error estimate in distribution. Furthermore, the geometric structure of the slow foliation is investigated: every fiber of the slow foliation parallels each other, with the slow manifold as a special fiber. In fact, when an arbitrarily chosen point of a fiber falls in the slow manifold, the fiber must be the slow manifold itself.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Slow Manifolds for Multi-Time-Scale Stochastic Evolutionary Systems

This article deals with invariant manifolds for infinite dimensional random dynamical systems with different time scales. Such a random system is generated by a coupled system of fast-slow stochastic evolutionary equations. Under suitable conditions, it is proved that an exponentially tracking random invariant manifold exists, eliminating the fast motion for this coupled system. It is further shown that if the scaling parameter tends to zero, the invariant manifold tends to a slow manifold which captures long time dynamics. As examples the results are applied to a few systems of coupled parabolic-hyperbolic partial differential equations, coupled parabolic partial differential-ordinary differential equations, and coupled hyperbolic-hyperbolic partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

State estimation under non-Gaussian Levy noise: A modified Kalman filtering method

The Kalman filter is extensively used for state estimation for linear systems under Gaussian noise. When non-Gaussian Lévy noise is present, the conventional Kalman filter may fail to be effective due to the fact that the non-Gaussian Lévy noise may have infinite variance. A modified Kalman filter for linear systems with non-Gaussian Lévy noise is devised. It works effectively with reasonable computational cost. Simulation results are presented to illustrate this non-Gaussian filtering method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic Averaging Principle for Dynamical Systems with Fractional Brownian Motion

Stochastic averaging for a class of stochastic differential equations (SDEs) with fractional Brownian motion, of the Hurst parameter H in the interval (1/2, 1), is investigated. An averaged SDE for the original SDE is proposed, and their solutions are quantitatively compared. It is shown that the solution of the averaged SDE converges to that of the original SDE in the sense of mean square and also in probability. It is further demonstrated that a similar averaging principle holds for SDEs under stochastic integral of pathwise backward and forward types. Two examples are presented and numerical simulations are carried out to illustrate the averaging principle.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An alternative expression of Di Paola and Falson's formula for stochastic dynamics

Di Paola and Falsone's formula is widely used in studying stochastic dynamics of nonlinear systems under Poisson white noise. In this short communication, an alternative expression is presented. Compared to Di Paola and Falsone's original expression, the alternative one is applicable under more general condition, and shows significantly improved performance in numerical implementation. The alternative expression turns out to be a special case of the Marcus integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approximating dynamics of a singularly perturbed stochastic wave equation with a random dynamical boundary condition

This work is concerned with a singularly perturbed stochastic nonlinear wave equation with a random dynamical boundary condition. A splitting skill is used to derive the approximating equation of the system in the sense of probability distribution, when the singular perturbation parameter is sufficiently small. The approximating equation is a stochastic parabolic equation when the power exponent of singular perturbation parameter is in $[1/2, 1)$, but a deterministic hyperbolic (wave) equation when the power exponent is in $(1, +\infty)$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymmetric non-Gaussian effects in a tumor growth model with immunization

The dynamical evolution of a tumor growth model, under immune surveillance and subject to asymmetric non-Gaussian $α$-stableLévy noise, is explored. The lifetime of a tumor staying in the range between the tumor-free state and the stable tumor state, and the likelihood of noise-inducing tumor extinction, are characterized by the mean exit time (also called mean residence time) and the escape probability, respectively. For various initial densities of tumor cells, the mean exit time and the escape probability are computed with different noise parameters. It is observed that unlike the Gaussian noise or symmetric non-Gaussian noise, the asymmetric non-Gaussian noise plays a constructive role in the tumor evolution in this simple model. By adjusting the noise parameters, the mean exit time can be shortened and the escape probability can be increased, simultaneously. This suggests that a tumor may be mitigated with higher probability in a shorter time, under certain external environmental stimuli.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Center manifolds for stochastic evolution equations

Stochastic invariant manifolds are crucial in modelling the dynamical behavior of dynamical systems under uncertainty. Under the assumption of exponential trichotomy, existence and smoothness of center manifolds for a class of stochastic evolution equations with linearly multiplicative noise are proved. The exponential attraction and approximation to center manifolds are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effective dynamics of stochastic wave equation with a random dynamical boundary condition

This work is devoted to the effective macroscopic dynamics of a weakly damped stochastic nonlinear wave equation with a random dynamical boundary condition. The white noises are taken into account not only in the model equation defined on a domain perforated with small holes, but also in the dynamical boundary condition on the boundaries of the small holes. An effective homogenized, macroscopic model is derived in the sense of probability distribution, which is a new stochastic wave equation on a unified domain, without small holes, with a usual static boundary condition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Elementary bifurcations for a simple dynamical system under non-Gaussian Levy noises

Nonlinear dynamical systems are sometimes under the influence of random fluctuations. It is desirable to examine possible bifurcations for stochastic dynamical systems when a parameter varies. A computational analysis is conducted to investigate bifurcations of a simple dynamical system under non-Gaussian α-stable Levy motions, by examining the changes in stationary probability density functions for the solution orbits of this stochastic system. The stationary probability density functions are obtained by numerically solving a non local Fokker-Planck equation. This allows numerically investigating phenomenological bifurcation, or P-bifurcation, for stochastic differential equations with non-Gaussian Levy noises.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ensemble Averaging for Dynamical Systems under Fast Oscillating Random Boundary Conditions

This paper is devoted to provide a theoretical underpinning for ensemble forecasting with rapid fluctuations in body forcing and in boundary conditions. Ensemble averaging principles are proved under suitable `mixing' conditions on random boundary conditions and on random body forcing. The ensemble averaged model is a nonlinear stochastic partial differential equation, with the deviation process (i.e., the approximation error process) quantified as the solution of a linear stochastic partial differential equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Escape Probability for Stochastic Dynamical Systems with Jumps

The escape probability is a deterministic concept that quantifies some aspects of stochastic dynamics. This issue has been investigated previously for dynamical systems driven by Gaussian Brownian motions. The present work considers escape probabilities for dynamical systems driven by non-Gaussian Lévy motions, especially symmetric $α$-stable Lévy motions. The escape probabilities are characterized as solutions of the Balayage-Dirichlet problems of certain partial differential-integral equations. Differences between escape probabilities for dynamical systems driven by Gaussian and non-Gaussian noises are highlighted. In certain special cases, analytic results for escape probabilities are given.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fokker-Planck equations for nonlinear dynamical systems driven by non-Gaussian Levy processes

The Fokker-Planck equations describe time evolution of probability densities of stochastic dynamical systems and are thus widely used to quantify random phenomena such as uncertainty propagation. For dynamical systems driven by non-Gaussian Lévy processes, however, it is difficult to obtain explicit forms of Fokker-Planck equations because the adjoint operators of the associated infinitesimal generators usually do not have exact formulation. In the present paper, Fokker- Planck equations are derived in terms of infinite series for nonlinear stochastic differential equations with non-Gaussian Lévy processes. A few examples are presented to illustrate the method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lévy Noise-Induced Stochastic Resonance in a Bistable System

Stochastic resonance phenomenon induced by non-Gaussian Lévy noise in a second-order bistable system is investigated. The signal-noise-ratio for different parameters is computed by an efficient numerical scheme. The influences of the noise intensity, stability index of Lévy noise and amplitude of external signal on the occurrence of stochastic resonance phenomenon are characterized. This implies that a high amplitude of signal not only enhances the output power spectrum of system but also promotes stochastic resonance, and a proper adjustment of Lévy noise intensity in a certain range enlarges the peak value of output power spectrum which is significant for stochastic resonance. Moreover, with the optimal damping parameter, lowering the stability index leads to larger fluctuations of Lévy noise, and further reduces the chance of the stochastic resonance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mean exit time and escape probability for dynamical systems driven by Levy noise

The mean first exit time and escape probability are utilized to quantify dynamical behaviors of stochastic differential equations with non-Gaussian alpha-stable type Levy motions. Both deterministic quantities are characterized by differential-integral equations(i.e.,differential equations with non local terms) but with different exterior conditions. The non-Gaussianity of noises manifests as nonlocality at the level of mean exit time and escape probability. An objective of this paper is to make mean exit time and escape probability as efficient computational tools, to the applied probability community, for quantifying stochastic dynamics. An accurate numerical scheme is developed and validated for computing the mean exit time and escape probability. Asymptotic solution for the mean exit time is given when the pure jump measure in the Levy motion is small. From both the analytical and numerical results, it is observed that the mean exit time depends strongly on the domain size and the value of alpha in the alpha-stable Levy jump measure. The mean exit time can measure which of the two competing factors in alpha-stable Levy motion, i.e. the jump frequency or the jump size, is dominant in helping a process exit a bounded domain. The escape probability is shown to vary with the underlying vector field(i.e.,drift). The mean exit time and escape probability could become discontinuous at the boundary of the domain, when the process is subject to certain deterministic potential and the value of alpha is in (0,1).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling nonlinear random vibration: Implication of the energy conservation law

Nonlinear random vibration under excitations of both Gaussian and Poisson white noises is considered. The model is based on stochastic differential equations, and the corresponding stochastic integrals are defined in such a way that the energy conservation law is satisfied. It is shown that Stratonovich integral and Di Paola-Falsone integral should be used for excitations of Gaussian and Poisson white noises, respectively, in order for the model to satisfy the underlining physical laws (e.g., energy conservation). Numerical examples are presented to illustrate the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-Gaussian dynamics of a tumor growth system with immunization

This paper is devoted to exploring the effects of non-Gaussian fluctuations on dynamical evolution of a tumor growth model with immunization, subject to non-Gaussian α-stable type Lévy noise. The corresponding deterministic model has two meaningful states which represent the state of tumor extinction and the state of stable tumor, respectively. To characterize the lifetime for different initial densities of tumor cells staying in the domain between these two states and the likelihood of crossing this domain, the mean exit time and the escape probability are quantified by numerically solving differential integral equations with appropriate exterior boundary conditions. The relationships between the dynamical properties and the noise parameters are examined. It is found that in the different stages of tumor, the noise parameters have different influence on the lifetime and the likelihood inducing tumor extinction. These results are relevant for determining efficient therapeutic regimes to induce the extinction of tumor cells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantifying Model Uncertainties in the Space of Probability Measures

Due to lack of scientific understanding, some mechanisms may be missing in mathematical modeling of complex phenomena in science and engineering. These mathematical models thus contain some uncertainties such as uncertain parameters. One method to estimate these parameters is based on pathwise observations, i.e., quantifying model uncertainty in the space of sample paths for system evolution. Another method is devised here to estimate uncertain parameters, or unknown system functions, based on experimental observations of probability distributions for system evolution. This is called the quantification of model uncertainties in the space of probability measures. A few examples are presented to demonstrate this method, analytically or numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulating Stochastic Inertial Manifolds by a Backward-Forward Approach

A numerical approach for the approximation of inertial manifolds of stochastic evolutionary equations with multiplicative noise is presented and illustrated. After splitting the stochastic evolutionary equations into a backward and a forward part, a numerical scheme is devised for solving this backward-forward stochastic system, and an ensemble of graphs representing the inertial manifold is consequently obtained. This numerical approach is tested in two illustrative examples: one is for a system of stochastic ordinary differential equations and the other is for a stochastic partial differential equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Synchronization of an evolving complex hyper-network

In this paper, the synchronization in a hyper-network of coupled dynamical systems is investigated for the first time. An evolving hyper-network model is proposed for better describing some complex systems. A concept of joint degree is introduced, and the evolving mechanism of hyper-network is given with respect to the joint degree. The hyper-degree distribution of the proposed evolving hyper-network is derived based on a rate equation method and obeys a power law distribution. Furthermore, the synchronization in a hyper-network of coupled dynamical systems is investigated for the first time. By calculating the joint degree matrix, several simple yet useful synchronization criteria are obtained and illustrated by several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Topological equivalence for discontinuous random dynamical systems and applications

After defining non-Gaussian Lévy processes for two-sided time, stochastic differential equations with such Lévy processes are considered. Solution paths for these stochastic differential equations have countable jump discontinuities in time. Topological equivalence (or conjugacy) for such an Itô stochastic differential equation and its transformed random differential equation is established. Consequently, a stochastic Hartman-Grobman theorem is proved for the linearization of the Itô stochastic differential equation. Furthermore, for Marcus stochastic differential equations,this topological equivalence is used to prove existence of global random attractors.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Approximation of invariant foliations for stochastic dynamical systems

Invariant foliations are geometric structures for describing and understanding the qualitative behaviors of nonlinear dynamical systems. For stochastic dynamical systems, however, these geometric structures themselves are complicated random sets. Thus it is desirable to have some techniques to approximate random invariant foliations. In this paper, invariant foliations are approximated for dynamical systems with small noisy perturbations, via asymptotic analysis. Namely, random invariant foliations are represented as a perturbation of the deterministic invariant foliations, with deviation errors estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global Well-posedness of the Stochastic Kuramoto-Sivashinsky Equation with Multiplicative Noise

Global well-posedness of the initial-boundary value problem for the stochastic Kuramoto-Sivashinsky equation in a bounded domain $D$ with a multiplicative noise is studied. It is shown that under suitable sufficient conditions, for any initial data $u_0\in L^2(D\times Ω)$ this problem has a unique global solution $u$ in the space $L^2(Ω,C([0,T],L^2({D})))$ for any $T>0$, and the solution map $u_0\mapsto u$ is Lipschitz continuous.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mean Exit Time and Escape Probability for a Tumor Growth System under Non-Gaussian Noise

Effects of non-Gaussian $α-$stable Lévy noise on the Gompertz tumor growth model are quantified by considering the mean exit time and escape probability of the cancer cell density from inside a safe or benign domain. The mean exit time and escape probability problems are formulated in a differential-integral equation with a fractional Laplacian operator. Numerical simulations are conducted to evaluate how the mean exit time and escape probability vary or bifurcates when $α$ changes. Some bifurcation phenomena are observed and their impacts are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolution systems of measures for stochastic flows

A new concept of {\em an evolution system of measures for stochastic flows} is considered. It corresponds to the notion of an invariant measure for random dynamical systems (or cocycles). The existence of evolution systems of measures for asymptotically compact stochastic flows is obtained. For a white noise stochastic flow, there exists a one to one correspondence between evolution systems of measures for a stochastic flow \emph{and} evolution systems of measures for the associated Markov transition semigroup. As an application, an alternative approach for evolution systems of measures of 2D stochastic Navier-Stokes equations with a time-periodic forcing term is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometric shape of invariant manifolds for a class of stochastic partial differential equations

Invariant manifolds play an important role in the study of the qualitative dynamical behaviors for nonlinear stochastic partial differential equations. However, the geometric shape of these manifolds is largely unclear. The purpose of the present paper is to try to describe the geometric shape of invariant manifolds for a class of stochastic partial differential equations with multiplicative white noises. The local geometric shape of invariant manifolds is approximated, which holds with significant likelihood. Furthermore, the result is compared with that for the corresponding deterministic partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large deviations for slow-fast stochastic partial differential equations

A large deviation principle is derived for stochastic partial differential equations with slow-fast components. The result shows that the rate function is exactly that of the averaged equation plus the fluctuating deviation which is a stochastic partial differential equation with small Gaussian perturbation. This also confirms the effectiveness of the approximation of the averaged equation plus the fluctuating deviation to the slow-fast stochastic partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantifying Uncertainties in Complex Systems

Uncertainties are abundant in complex systems. Mathematical models for these systems thus contain random effects or noises. The models are often in the form of stochastic differential equations, with some parameters to be determined by observations. The stochastic differential equations may be driven by Brownian motion, fractional Brownian motion or Lévy motion. After a brief overview of recent advances in estimating parameters in stochastic differential equations, various numerical algorithms for computing parameters are implemented. The numerical simulation results are shown to be consistent with theoretical analysis. Moreover, for fractional Brownian motion and $α-$stable Lévy motion, several algorithms are reviewed and implemented to numerically estimate the Hurst parameter $H$ and characteristic exponent $α$.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Impact of Boundary Conditions on Entrainment and Transport in Gravity Currents

Gravity currents have been studied numerically and experimentally both in the laboratory and in the ocean. The question of appropriate boundary conditions is still challenging for most complex flows. Gravity currents make no exception - appropriate, physically and mathematically sound boundary conditions are yet to be found. This task is further complicated by the technical limitations imposed by the current oceanographic techniques. In this paper, we make a first step toward a better understanding of the impact of boundary conditions on gravity currents. Specifically, we use direct numerical simulations to investigate the effect that the popular Neumann, and less popular Dirichlet boundary conditions on the bottom continental shelf have on the entrainment and transport of gravity currents. The finding is that gravity currents under these two different boundary conditions differ most in the way they transport heat from the top towards the bottom. This major difference occurs at medium temperature ranges. Entrainment and transport at high temperatures also show significant differences.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Determining functionals for random partial differential equations

Determining functionals are tools to describe the finite dimensional long-term dynamics of infinite dimensional dynamical systems. There also exist several applications to infinite dimensional {\em random} dynamical systems. In these applications the convergence condition of the trajectories of an infinite dimensional random dynamical system with respect to a finite set of linear functionals is assumed to be either in mean or {\em exponential} with respect to the convergence almost surely. In contrast to these ideas we introduce a convergence concept which is based on the convergence in probability. By this ansatz we get rid of the assumption of exponential convergence. In addition, setting the random terms to zero we obtain usual deterministic results.\\ We apply our results to the 2D Navier - Stokes equations forced by a white noise.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Dynamics of the Thermohaline Circulation under Wind forcing

The ocean thermohaline circulation, also called meridional overturning circulation, is caused by water density contrasts. This circulation has large capacity of carrying heat around the globe and it thus affects the energy budget and further affects the climate. We consider a thermohaline circulation model in the meridional plane under external wind forcing. We show that, when there is no wind forcing, the stream function and the density fluctuation (under appropriate metrics) tend to zero exponentially fast as time goes to infinity. With rapidly oscillating wind forcing, we obtain an averaging principle for the thermohaline circulation model. This averaging principle provides convergence results and comparison estimates between the original thermohaline circulation and the averaged thermohaline circulation, where the wind forcing is replaced by its time average. This establishes the validity for using the averaged thermohaline circulation model for numerical simulations at long time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Enstrophy Dynamics of Stochastically Forced Large-Scale Geophysical Flows

Enstrophy is an averaged measure of fluid vorticity. This quantity is particularly important in {\em rotating} geophysical flows. We investigate the dynamical evolution of enstrophy for large-scale quasi-geostrophic flows under random wind forcing. We obtain upper bounds on the enstrophy, as well as results establishing its Hölder continuity and describing the small-time asymptotics.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Ergodicity of Stochastically Forced Large Scale Geophysical Flows

We investigate the ergodicity of 2D large scale quasigeostrophic flows under random wind forcing. We show that the quasigeostrophic flows are ergodic under suitable conditions on the random forcing and on the fluid domain, and under no restrictions on viscosity, Ekman constant or Coriolis parameter. When these conditions are satisfied, then for any observable of the quasigeostrophic flows, its time average approximates the statistical ensemble average, as long as the time interval is sufficiently long.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Probabilistic Dynamics of Two-Layer Geophysical Flows

The two-layer quasigeostrophic flow model is an intermidiate system between the single-layer 2D barotropic flow model and the continuously stratified, 3D baroclinic flow model. This model is widely used to investigate basic mechanisms in geophysical flows, such as baroclinic effects, the Gulf Stream and subtropical gyres. The wind forcing acts only on the top layer. We consider the two-layer quasigeostrophic model under stochastic wind forcing. We first transformed this system into a coupled system of random partial differential equations and then show that the asymptotic probabilistic dynamics of this system depends only on the top fluid layer. Namely, in the probability sense and asymptotically, the dynamics of the two-layer quasigeostrophic fluid system is determinied by the top fluid layer, or, the bottom fluid layer is slaved by the top fluid layer. This conclusion is true provided that the Wiener process and the fluid parameters satisfy a certain condition. In particular, this latter condition is satisfied when the trace of the covariance operator of the Wiener process is controled by a certain upper bound, and the Ekman constant $r$ is sufficiently large. Note that the generalized time derivative of the Wiener process models the fluctuating part of the wind stress forcing on the top fluid layer, and the Ekman constant $r$ measures the rate for vorticity decay due to the friction in the bottom Ekman layer.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Averaging Principle for Quasi-Geostrophic Motions under Rapidly Oscillating Forcing

In this paper, the averaging principle for quasi-geostrophic motions with rapidly oscillating forcing is proved, both on finite but large time intervals and on the entire time axis. This includes comparison estimate, stability estimate, and convergence result between quasi-geostrophic motions and its averaged motions. Furthermore, the existence of almost periodic quasi-geostrophic motions and attractor convergence are also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The 3D Quasigeostrophic Fluid Dynamics under Random Forcing on Boundary

The three-dimensional baroclinic quasigeostrophic flow model has been widely used to study basic mechanisms in oceanic flows and climate dynamics. In this paper, we consider this flow model under random wind forcing and time-periodic fluctuations on fluid boundary (the interface between the oceans and the atmosphere). The time-periodic fluctuations are due to periodic rotation of the earth and thus periodic exposure of the earth to the solar radiation. After establishing the well-posedness of the baroclinic quasigeostrophic flow model in the state space, we demonstrate the existence of the random attractors, again in the state space. We also discuss the relevance of our result to climate modeling.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Chaotic Properties of Subshifts Generated by a Non-Periodic Recurrent Orbit

The chaotic properties of some subshift maps are investigated. These subshifts are the orbit closures of certain non-periodic recurrent points of a shift map. We first provide a review of basic concepts for dynamics of continuous maps in metric spaces. These concepts include nonwandering point, recurrent point, eventually periodic point, scrambled set, sensitive dependence on initial conditions, Robinson chaos, and topological entropy. Next we review the notion of shift maps and subshifts. Then we show that the one-sided subshifts generated by a non-periodic recurrent point are chaotic in the sense of Robinson. Moreover, we show that such a subshift has an infinite scrambled set if it has a periodic point. Finally, we give some examples and discuss the topological entropy of these subshifts, and present two open problems on the dynamics of subshifts.

0 citations0 reviewsNo code linkedOpen access