Fractional Advection Dispersion Equation and Martingale Problem: Difference between pages

From nonlocal pde
(Difference between pages)
Jump to navigation Jump to search
imported>Nestor
No edit summary
 
imported>Nestor
No edit summary
 
Line 1: Line 1:
{{stub}}
{{stub}}


The Fractional Advection Dispersion (fADE) was introduced <ref name="MBB1999"/> (and further developed and justified <ref name="BMW"/> <ref name="MBB2001"/> <ref name ="BMW2001"/>) to better account for super-diffusive spreading of tracer particles in aquifers.  It has quite a few variations at this point, but a couple of characteristic examples are
(Classical Martingale Problem) Given a (local) linear operator $\mathcal{L}: C^2(\mathbb{R}^n) \to C(\mathbb{R}^n)$ Martingale Problem for $\mathcal{L}$ consists in finding for each $x_0 \in \mathbb{R}^d$ a probability measure $\mathbb{P}^{x_0}$ over the space of all continuous functions $X: [0,+\infty) \to \mathbb{R}^d$  such that


\[
\begin{equation*}
u_t = c\frac{\partial u}{\partial x} + p\frac{\partial^\alpha u}{\partial x^\alpha} + q\frac{\partial^\alpha u}{\partial (-x)^\alpha} ,
\mathbb{P}^{x_0}\left ( X(0)=x_0 \right ) = 1
\]
\end{equation*}
where the fractional derivatives are defined as
\begin{equation}
\frac{\partial^{\alpha}}{\partial e^{\alpha}} u(x) = c_{n,\alpha}\int_0^\infty \frac{u(x)-u(x-ye)}{y}y^{-\alpha}dy\ \ \text{and}\ \ \frac{\partial^{\alpha}}{\partial (-e)^{\alpha}} u(x) = c_{n,\alpha}\int_0^\infty \frac{u(x)-u(x+ye)}{y}y^{-\alpha}dy .
\end{equation}
It is important to point out that these fractional derivatives are examples of one-dimensional [[Linear integro-differential operator | linear integro-differential operators]] with non-symmetric kernels (in this case, e.g. $K(y)=\mathbb{1}_{\{y\geq 0\}}|y|^{-1-\alpha}$).


and whenever $f \in C^2(\mathbb{R}^d)$ we have that


== See also ==
\begin{equation*}
f(X(t))-f(X(0))-\int_0^t \mathcal{L} f(X(s)) \;ds
\end{equation*}


* [[Drift-diffusion equations]]
is a [[Martingale| Local Martingale]] under $\mathbb{P}^{x_0}$. If for any $x \in \mathbb{R}^d$ there exists a unique $\mathbb{P}^x$ satisfying the above conditions we say that the Martingale Problem for $\mathcal{L}$ is well posed.


== References ==
Note that existence for the Martingale Problem implies existence of solutions to the Cauchy problem for the operator $\mathcal{L}$ with initial data in $C^2(\mathbb{R}^d)$, indeed, if $f \in C^2(\mathbb{R}^d)$ then we can define
{{reflist|refs=
 
<ref name="MBB1999">{{Citation | last1=Meerschaert | first1=M.M. | last2=Benson | first2=D. J. | last3=Bäumer | first3=B. | title=Multidimensional advection and fractional dispersion | publisher=APS | year=1999 | journal=[[Physical Review|Physical Review E]] | issn=1539-3755 | volume=59 | issue=5 | pages=5026}}</ref>
\begin{equation*}
<ref name="BMW">{{Citation | last1=Benson | first1=D. J. | last2=Wheatcraft | first2=S.W. | last3=Meerschaert | first3=M.M. | title=Application of a fractional advection-dispersion equation | year=2000 | journal=Water Resources Research | volume=36 | issue=6 | pages=1403–1412}}</ref>
u(x,t) = \mathbb{E}_{\mathbb{P}^x} [ f(X_t) ]
<ref name="MBB2001">{{Citation | last1=Meerschaert | first1=M.M. | last2=Benson | first2=D. J. | last3=Baeumer | first3=B. | title=Operator Lévy motion and multiscaling anomalous diffusion | publisher=APS | year=2001 | journal=[[Physical Review|Physical Review E]] | issn=1539-3755 | volume=63 | issue=2 | pages=021112}}</ref>
\end{equation*}
<ref name="BMW2001">{{Citation | last1=Benson | first1=D. J. | last2=Schumer | first2=R. | last3=Meerschaert | first3=M.M. | last4=Wheatcraft | first4=S.W. | title=Fractional dispersion, Lévy motion, and the MADE tracer tests | publisher=[[Springer-Verlag]] | location=Berlin, New York | year=2001 | journal=Transport in Porous Media | volume=42 | issue=1 | pages=211–240}}</ref>
 
}}
Then it can be shown that $u(x,t)$ is a classical solution to the problem
 
\begin{equation*}
\left \{ \begin{array}{ll}
u_t  = \mathcal{L}u & \text{ in } \mathbb{R}^d\times \mathbb{R}_+\\
u\;\;  = f & \text{ for } t=0
\end{array}\right.
\end{equation*}
 
Thus at first glance it would seem that  well-posedness for the Martingale Problem for $\mathcal{L}$ is a stronger fact than that the well-posedness for the corresponding Cauchy problem, but often one can turn this around and show the well-posedness for the Martingale Problem after understanding the Cauchy problem for $\mathcal{L}$ well enough.
 
The oldest example of a Martingale Problem is when $\mathcal{L}= \Delta$, the classical Laplacian, in this case the unique solution $\mathbb{P}^{x}$ corresponds to the standard Brownian motion or Wiener process. Using this as a base case and using basic Ito Calculus one can solve the corresponding problem for a generic second order operator
 
\begin{equation*}
\mathcal{L} u = \text{Tr}(A(x)D^2u(x))+b(x)\cdot \nabla u(x)+c(x)u(x)
\end{equation*}
 
where $A(x)\geq 0$, $c(x)\leq 0$ and $A,b,c$ are all Lipschitz functions of $x$. The case where the coefficients are merely continuous is much harder and it is an important result of Stroock and Varadhan that says still in this case the Martingale Problem is well posed.

Latest revision as of 19:48, 19 November 2012

This article is a stub. You can help this nonlocal wiki by expanding it.

(Classical Martingale Problem) Given a (local) linear operator $\mathcal{L}: C^2(\mathbb{R}^n) \to C(\mathbb{R}^n)$ Martingale Problem for $\mathcal{L}$ consists in finding for each $x_0 \in \mathbb{R}^d$ a probability measure $\mathbb{P}^{x_0}$ over the space of all continuous functions $X: [0,+\infty) \to \mathbb{R}^d$ such that

\begin{equation*} \mathbb{P}^{x_0}\left ( X(0)=x_0 \right ) = 1 \end{equation*}

and whenever $f \in C^2(\mathbb{R}^d)$ we have that

\begin{equation*} f(X(t))-f(X(0))-\int_0^t \mathcal{L} f(X(s)) \;ds \end{equation*}

is a Local Martingale under $\mathbb{P}^{x_0}$. If for any $x \in \mathbb{R}^d$ there exists a unique $\mathbb{P}^x$ satisfying the above conditions we say that the Martingale Problem for $\mathcal{L}$ is well posed.

Note that existence for the Martingale Problem implies existence of solutions to the Cauchy problem for the operator $\mathcal{L}$ with initial data in $C^2(\mathbb{R}^d)$, indeed, if $f \in C^2(\mathbb{R}^d)$ then we can define

\begin{equation*} u(x,t) = \mathbb{E}_{\mathbb{P}^x} [ f(X_t) ] \end{equation*}

Then it can be shown that $u(x,t)$ is a classical solution to the problem

\begin{equation*} \left \{ \begin{array}{ll} u_t = \mathcal{L}u & \text{ in } \mathbb{R}^d\times \mathbb{R}_+\\ u\;\; = f & \text{ for } t=0 \end{array}\right. \end{equation*}

Thus at first glance it would seem that well-posedness for the Martingale Problem for $\mathcal{L}$ is a stronger fact than that the well-posedness for the corresponding Cauchy problem, but often one can turn this around and show the well-posedness for the Martingale Problem after understanding the Cauchy problem for $\mathcal{L}$ well enough.

The oldest example of a Martingale Problem is when $\mathcal{L}= \Delta$, the classical Laplacian, in this case the unique solution $\mathbb{P}^{x}$ corresponds to the standard Brownian motion or Wiener process. Using this as a base case and using basic Ito Calculus one can solve the corresponding problem for a generic second order operator

\begin{equation*} \mathcal{L} u = \text{Tr}(A(x)D^2u(x))+b(x)\cdot \nabla u(x)+c(x)u(x) \end{equation*}

where $A(x)\geq 0$, $c(x)\leq 0$ and $A,b,c$ are all Lipschitz functions of $x$. The case where the coefficients are merely continuous is much harder and it is an important result of Stroock and Varadhan that says still in this case the Martingale Problem is well posed.