A Discontinuous Galerkin Method for Approximating the Stationary Distribution of Stochastic Fluid-Fluid Processes

The stochastic fluid-fluid model (SFFM) is a Markov process {(Xt,Yt,φt),t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{(X_t,Y_t,\varphi _t),t\ge 0\}$$\end{document}, where {φt,t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\varphi _t,{t\ge 0}\}$$\end{document} is a continuous-time Markov chain, the first fluid, {Xt,t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{X_t,t\ge 0\}$$\end{document}, is a classical stochastic fluid process driven by {φt,t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\varphi _t,t\ge 0\}$$\end{document}, and the second fluid, {Yt,t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_t,t\ge 0\}$$\end{document}, is driven by the pair {(Xt,φt),t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{(X_t,\varphi _t),t\ge 0\}$$\end{document}. Operator-analytic expressions for the stationary distribution of the SFFM, in terms of the infinitesimal generator of the process {(Xt,φt),t≥0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{(X_t,\varphi _t),t\ge 0\}$$\end{document}, are known. However, these operator-analytic expressions do not lend themselves to direct computation. In this paper the discontinuous Galerkin (DG) method is used to construct approximations to these operators, in the form of finite dimensional matrices, to enable computation. The DG approximations are used to construct approximations to the stationary distribution of the SFFM, and results are verified by simulation. The numerics demonstrate that the DG scheme can have a superior rate of convergence compared to other methods.