Adaptive Numerical Solutions of Fokker-Planck Equations in Computational Uncertainty Quantification

Author(s):  
Mani Razi ◽  
Peter Attar ◽  
Prakash Vedula
Entropy ◽  
2018 ◽  
Vol 20 (10) ◽  
pp. 760 ◽  
Author(s):  
Johan Anderson ◽  
Sara Moradi ◽  
Tariq Rafiq

The numerical solutions to a non-linear Fractional Fokker–Planck (FFP) equation are studied estimating the generalized diffusion coefficients. The aim is to model anomalous diffusion using an FFP description with fractional velocity derivatives and Langevin dynamics where Lévy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space. Distribution functions are found using numerical means for varying degrees of fractionality of the stable Lévy distribution as solutions to the FFP equation. The statistical properties of the distribution functions are assessed by a generalized normalized expectation measure and entropy and modified transport coefficient. The transport coefficient significantly increases with decreasing fractality which is corroborated by analysis of experimental data.


2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Limei Yan

A relatively new iterative Laplace transform method, which combines two methods; the iterative method and the Laplace transform method, is applied to obtain the numerical solutions of fractional Fokker-Planck equations. The method gives numerical solutions in the form of convergent series with easily computable components, requiring no linearization or small perturbation. The numerical results show that the approach is easy to implement and straightforward when applied to space-time fractional Fokker-Planck equations. The method provides a promising tool for solving space-time fractional partial differential equations.


2020 ◽  
Author(s):  
Christian Behnken ◽  
Matthias Wächter ◽  
Joachim Peinke

Abstract. The most intermittent behavior of atmospheric turbulence is found for very short time scales. Based on a concatenation of conditional probability density functions (cpdfs) of nested wind speeds increments, inspired by a Markov process in scale, we derive a short-time predictor for wind speed fluctuations around a non-stationary mean value and with a corresponding non-stationary variance. As a new quality this short time predictor enables a multipoint reconstruction of wind data. The used cpdfs are (1) directly estimated from historical data from the offshore research platform FINO1 and (2) obtained from numerical solutions of a family of Fokker-Planck equations in the scale domain. The explicit forms of the Fokker-Planck equations are estimated from the given wind data. A good agreement between the statistics of the generated synthetic wind speed fluctuations and the measured is found even on time scales below 1 s. This shows that our approach captures the short-time dynamics of real wind speed fluctuations very well. Our method is extended by taking the non-stationarity of the mean wind speed and its non-stationary variance into account.


Author(s):  
Manoj Kumar ◽  
Sapna Pandit

Purpose – The purpose of this paper is to discuss the application of the Haar wavelets for solving linear and nonlinear Fokker-Planck equations with appropriate initial and boundary conditions. Design/methodology/approach – Haar wavelet approach converts the problems into a system of linear algebraic equations and the obtained system is solved by Gauss-elimination method. Findings – The accuracy of the proposed scheme is demonstrated on three test examples. The numerical solutions prove that the proposed method is reliable and yields compatible results with the exact solutions. The scheme provides better results than the schemes [9, 14]. Originality/value – The developed scheme is a new scheme for Fokker-Planck equations. The scheme based on Haar wavelets is expended for nonlinear partial differential equations with variable coefficients.


Entropy ◽  
2020 ◽  
Vol 22 (8) ◽  
pp. 802
Author(s):  
Dimitra Maoutsa ◽  
Sebastian Reich ◽  
Manfred Opper

Fokker–Planck equations are extensively employed in various scientific fields as they characterise the behaviour of stochastic systems at the level of probability density functions. Although broadly used, they allow for analytical treatment only in limited settings, and often it is inevitable to resort to numerical solutions. Here, we develop a computational approach for simulating the time evolution of Fokker–Planck solutions in terms of a mean field limit of an interacting particle system. The interactions between particles are determined by the gradient of the logarithm of the particle density, approximated here by a novel statistical estimator. The performance of our method shows promising results, with more accurate and less fluctuating statistics compared to direct stochastic simulations of comparable particle number. Taken together, our framework allows for effortless and reliable particle-based simulations of Fokker–Planck equations in low and moderate dimensions. The proposed gradient–log–density estimator is also of independent interest, for example, in the context of optimal control.


2020 ◽  
Vol 5 (3) ◽  
pp. 1211-1223
Author(s):  
Christian Behnken ◽  
Matthias Wächter ◽  
Joachim Peinke

Abstract. The most intermittent behaviour of atmospheric turbulence is found for very short timescales. Based on a concatenation of conditional probability density functions (cpdf's) of nested wind speed increments, inspired by a Markov process in scale, we derive a short-time predictor for wind speed fluctuations around a non-stationary mean value and with a corresponding non-stationary variance. As a new quality this short-time predictor enables a multipoint reconstruction of wind data. The used cpdf's are (1) directly estimated from historical data from the offshore research platform FINO1 and (2) obtained from numerical solutions of a family of Fokker–Planck equations in the scale domain. The explicit forms of the Fokker–Planck equations are estimated from the given wind data. A good agreement between the statistics of the generated and measured synthetic wind speed fluctuations is found even on timescales below 1 s. This shows that our approach captures the short-time dynamics of real wind speed fluctuations very well. Our method is extended by taking the non-stationarity of the mean wind speed and its non-stationary variance into account.


Author(s):  
Luca Giuggioli ◽  
Zohar Neu

Noise and time delays, or history-dependent processes, play an integral part in many natural and man-made systems. The resulting interplay between random fluctuations and time non-locality are essential features of the emerging complex dynamics in non-Markov systems. While stochastic differential equations in the form of Langevin equations with additive noise for such systems exist, the corresponding probabilistic formalism is yet to be developed. Here we introduce such a framework via an infinite hierarchy of coupled Fokker–Planck equations for the n -time probability distribution. When the non-Markov Langevin equation is linear, we show how the hierarchy can be truncated at n  = 2 by converting the time non-local Langevin equation to a time-local one with additive coloured noise. We compare the resulting Fokker–Planck equations to an earlier version, solve them analytically and analyse the temporal features of the probability distributions that would allow to distinguish between Markov and non-Markov features. This article is part of the theme issue ‘Nonlinear dynamics of delay systems’.


Sign in / Sign up

Export Citation Format

Share Document