Abstract
Real life convection-diffusion problems are characterized by their inherent or externally induced uncertainties in the design parameters. This paper presents a spectral stochastic finite element semi-Lagrangian method for numerical solution of convection-diffusion equations with uncertainty. Using the spectral decomposition, the stochastic variational problem is reformulated to a set of deterministic variational problems to be solved for each Wiener polynomial chaos. To obtain the chaos coefficients in the corresponding deterministic convection-diffusion equations, we implement a semi-Lagrangian method in the finite element framework. Once this representation is computed, statistics of the numerical solution can be easily evaluated. These numerical techniques associate the geometrical flexibility of the finite element method with the ability offered by the semi-Lagrangian method to solve convection-dominated problems using time steps larger than its Eulerian counterpart. Numerical results are shown for a convection-diffusion problem driven with stochastic velocity and for an incompressible viscous flow problem with a random force. In both examples, the proposed method demonstrates its ability to better maintain the shape of the solution in the presence of uncertainties and steep gradients.
Similar content being viewed by others
References
Adler, R.J., Taylor, J.E.: Random Fields and Geometry. Birkhäuser, Basel (2005)
Armaly, B.F., Durst, F., Periera, J.C.F., Schönung, B.: Experimental and theoretical investigation of backward-facing step flow. J. Fluid Mech. 127, 473–496 (1983)
Asokan, V., Narayanan, B., Zabaras, N.: Variational multiscale stabilized FEM formulations for transport equations: stochastic advection-diffusion and incompressible stochastic Navier-Stokes equations. J. Comput. Physics. 202, 94–133 (2005)
Babuska, I., Tempone, R., Zouraris, G.E.: Galerkin finite element approximations of stochastic elliptic partial differential equations. SIAM J. Numer. Anal. 42, 800–825 (2004)
Cameron, R.H., Martin, W.T.: The orthogonal development of nonlinear functionals in series of Fourier-Hermite functionals. Ann. Math. 48, 385–392 (1947)
Chorin, A.J.: Gaussian fields and random flow. J. Fluid Mech. 63, 21–32 (1947)
Dagan, G., Neuman, S.P.: Subsurface Flow and Transport: A Stochastic Approach. International Hydrology Series. Cambridge University Press, Cambridge (1997)
Douglas, J., Russell, T.F.: Numerical methods for convection dominated diffusion problems based on combining the method of characteristics with finite elements or finite differences. SIAM J. Numer. Anal. 19, 871–885 (1982)
El-Amrani, M., Seaïd, M.: A finite element modified method of characteristics for convective heat transport. Numer. Methods Partial Differ. Equ. 24, 776–798 (2008)
El-Amrani, M., Seaïd, M.: Eulerian-Lagrangian time-stepping methods for convection-dominated problems. Int. J. Comput. Math. 85, 421–439 (2008)
El-Amrani, M., Seaïd, M.: Convergence and stability of finite element modified method of characteristics for the incompressible Navier-Stokes equations. J. Numer. Math. 15, 101–135 (2007)
El-Amrani, M., Seaïd, M.: Numerical simulation of natural and mixed convection flows by Galerkin-characteristic method. Int. J. Numer. Meth. Fluids. 53, 1819–1845 (2007)
Ghanem, R.: Probabilistic characterization of transport in heterogeneous media. Comput. Methods Appl. Mech. Eng. 158, 199–220 (1998)
Gottlieb, S., Shu, C.W., Tadmor, E.: Strong stability preserving high order time integration methods. SIAM Rev. 43, 89–112 (2001)
Gottlieb, D., Orszag, S.A.: Numerical Analysis of Spectral Methods: Theory and Applications. SIAM, Philadelphia (1977)
Kim, J.U.: On the stochastic Euler equations in a two-dimensional domain. SIAM J. Math. Anal. 33, 1211–1227 (2002)
Kim, J., Moin, P.: Application of a fractional-step method to incompressible Navier-Stokes equations. J. Comput. Phys. 59, 308–323 (1985)
Le Maître, O.P., Knio, O.M., Debusschere, B.J., Najm, H.N., Ghanem, R.G.: A multigrid solver for two-dimensional stochastic diffusion equations. Comput. Methods Appl. Mech. Eng. 192, 41–42 (2003)
Manouzi, H., Seaïd, M., Zahri, M.: Wick-stochastic finite element solution of reaction-diffusion problems. J. Comput. Appl. Math. 203, 516–532 (2007)
Matthies, H.G., Brenner, C.E., Bucher, C.G., Soares, C.G.: Uncertainties in probabilistic numerical analysis of structures and solids-Stochastic finite elements. Struct. Saf. 19, 283–336 (1997)
Morgan, K., Periaux, J., Thomasset, F.: Analysis of Laminar Flow over a Backward Facing Step. Notes on Numerical Fluid Mechanics, vol. 9 (1984)
Phenix, B.D., Dinaro, J.L., Tatang, M.A., Tester, J.W., Howard, J.B., McRae, G.J.: Incorporation of parametric uncertainty into complex kinetic mechanisms: application to hydrogen oxidation in supercritical water. Combust. Flame. 112, 132–146 (1998)
Pironneau, O.: On the transport-diffusion algorithm and its applications to the Navier-Stokes equations. Numer. Math. 38, 309–332 (1982)
Robert, A.: A stable numerical integration scheme for the primitive meteorological equations. Atmos. Ocean 19, 35–46 (1981)
Seaïd, M.: On the quasi-monotone modified method of characteristics for transport-diffusion problems with reactive sources. Comput. Methods Appl. Math. 2, 186–210 (2002)
Seaïd, M.: Semi-Lagrangian integration schemes for viscous incompressible flows. Comput. Methods Appl. Math. 4, 1–18 (2002)
Shu, C.W.: Total variation diminishing time discretizations. SIAM J. Sci. Stat. Comput. 9, 1073–1084 (1988)
Temperton, C., Staniforth, A.: An efficient two-time-level semi-Lagrangian semi-implicit integration scheme. Q. J. R. Meteorol. Soc. 113, 1025–1039 (1987)
Wan, X., Xiu, D., Karniadakis, G.E.: Stochastic solutions for the two-dimensional advection-diffusion equation. SIAM J. Sci. Comput. 26, 578–590 (2004)
Wiener, N.: The homogeneous chaos. Am. J. Math. 60, 897–936 (1938)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
El-Amrani, M., Seaïd, M. A Spectral Stochastic Semi-Lagrangian Method for Convection-Diffusion Equations with Uncertainty. J Sci Comput 39, 371–393 (2009). https://doi.org/10.1007/s10915-009-9273-5
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-009-9273-5