Numerical Simulation of Feller’s Diffusion Equation

This article is devoted to <span style="font-variant: small-caps;">Feller</span>&#8217;s diffusion equation, which arises naturally in probability and physics (e.g., wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties...

Full description

Bibliographic Details
Main Author: Denys Dutykh
Format: Article
Language:English
Published: MDPI AG 2019-11-01
Series:Mathematics
Subjects:
Online Access:https://www.mdpi.com/2227-7390/7/11/1067
id doaj-09c81470ba8d44419081fbfc85977b89
record_format Article
spelling doaj-09c81470ba8d44419081fbfc85977b892020-11-25T01:56:33ZengMDPI AGMathematics2227-73902019-11-01711106710.3390/math7111067math7111067Numerical Simulation of Feller’s Diffusion EquationDenys Dutykh0Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LAMA, 73000 Chambéry, FranceThis article is devoted to <span style="font-variant: small-caps;">Feller</span>&#8217;s diffusion equation, which arises naturally in probability and physics (e.g., wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties, we reformulate this equation using some ideas from the <span style="font-variant: small-caps;">Lagrangian</span> fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation, and the corresponding <span style="font-variant: small-caps;">Matlab</span> code is provided with this article under an open source license.https://www.mdpi.com/2227-7390/7/11/1067feller equationparabolic equationslagrangian schemefokker–planck equationprobability distribution
collection DOAJ
language English
format Article
sources DOAJ
author Denys Dutykh
spellingShingle Denys Dutykh
Numerical Simulation of Feller’s Diffusion Equation
Mathematics
feller equation
parabolic equations
lagrangian scheme
fokker–planck equation
probability distribution
author_facet Denys Dutykh
author_sort Denys Dutykh
title Numerical Simulation of Feller’s Diffusion Equation
title_short Numerical Simulation of Feller’s Diffusion Equation
title_full Numerical Simulation of Feller’s Diffusion Equation
title_fullStr Numerical Simulation of Feller’s Diffusion Equation
title_full_unstemmed Numerical Simulation of Feller’s Diffusion Equation
title_sort numerical simulation of feller’s diffusion equation
publisher MDPI AG
series Mathematics
issn 2227-7390
publishDate 2019-11-01
description This article is devoted to <span style="font-variant: small-caps;">Feller</span>&#8217;s diffusion equation, which arises naturally in probability and physics (e.g., wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties, we reformulate this equation using some ideas from the <span style="font-variant: small-caps;">Lagrangian</span> fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation, and the corresponding <span style="font-variant: small-caps;">Matlab</span> code is provided with this article under an open source license.
topic feller equation
parabolic equations
lagrangian scheme
fokker–planck equation
probability distribution
url https://www.mdpi.com/2227-7390/7/11/1067
work_keys_str_mv AT denysdutykh numericalsimulationoffellersdiffusionequation
_version_ 1724979342262403072