Skip to main content
Springer Nature Link
Log in

Ensemble Time-Stepping Algorithm for the Convection-Diffusion Equation with Random Diffusivity

  • Published:
Journal of Scientific Computing Aims and scope Submit manuscript

Abstract

In this paper, we develop two ensemble time-stepping algorithms to solve the convection-diffusion equation with random diffusion coefficients, forcing terms and initial conditions based on the pseudo-spectral stochastic collocation method. The key step of the pseudo-spectral stochastic collocation method is to solve a number of deterministic problems derived from the original stochastic convection-diffusion equation. In general, a common way to solve the set of deterministic problems is by using the backward differentiation formula, which requires us to store the coefficient matrix and right-hand-side vector multiple times, and solve them one by one. However, the proposed algorithm only need to solve a single linear system with one shared coefficient matrix and multiple right-hand-side vectors, reducing both storage required and computational cost of the solution process. The stability and error analysis of the first- and second-order ensemble time-stepping algorithms are provided. Several numerical experiments are presented to confirm the theoretical analyses and verify the feasibility and effectiveness of the proposed method.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4

Similar content being viewed by others

References

  1. Xiu, D., Hesthaven, J.S.: High-order collocation methods for differential equations with random inputs. SIAM J. Sci. Comput. 27, 1118–1139 (2005)

    Article  MathSciNet  MATH  Google Scholar 

  2. Xiu, D.: Efficient collocational approach for parametric uncertainty analysis. Commun. Comput. Phys. 2, 293–309 (2007)

    MathSciNet  MATH  Google Scholar 

  3. Xiu, D.: Fast numerical methods for stochastic computations: a review. Commun. Comput. Phys. 5, 242–272 (2009)

    MathSciNet  MATH  Google Scholar 

  4. Xiu, D.: Stochastic collocation methods: a survey. In: Ghanem, R., Higdon, D., Owhadi, H. (eds.) Handbook of Uncertainty Quantification. Springer, Cham (2017)

    Google Scholar 

  5. 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)

    Article  MathSciNet  MATH  Google Scholar 

  6. Babuska, I., Nobile, F., Tempone, R.: A stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Numer. Anal. 45, 1005–1034 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  7. Babuska, I., Nobile, F., Tempone, R.: A stochastic collocation method for elliptic partial differential equations with random input data. SIAM Rev. 52, 317–355 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  8. Nobile, F., Tempone, R., Webster, C.G.: A sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J. Numer. Anal. 46, 2309–2345 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  9. Zhang, G., Gunzburger, M.: Error analysis of a stochastic collocation method for parabolic partial differential equations with random input data. SIAM J. Numer. Anal. 50, 1922–1940 (2012)

    Article  MathSciNet  MATH  Google Scholar 

  10. Zhang, Q., Li, Z., Zhang, Z.: A sparse grid stochastic collocation method for elliptic interface problems with random input. J. Sci. Comput. 67, 262–280 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  11. Tang, T., Zhou, T.: Recent developments in high order numerical methods for uncertainty quantification. Sci. Sin. Math. 45, 891–928 (2015)

    Google Scholar 

  12. Tang, T., Zhou, T.: On discrete least-squares projection in unbounded domain with random evaluations and its application to parametric uncertainty quantification. SIAM J. Sci. Comput. 36, A2272–A2295 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  13. Jakeman, J.D., Narayan, A., Zhou, T.: A generalized sampling and preconditioning scheme for sparse approximation of polynomial chaos expansions. SIAM J. Sci. Comput. 39, A1114–A1144 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  14. Narayan, A., Zhou, T.: Stochastic collocation on unstructured multivariate meshes. Commun. Comput. Phys. 18, 1–36 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  15. Guo, L., Narayan, A., Zhou, T., Chen, Y.: Stochastic collocation methods via \(\ell _1\) minimization using randomized quadratures. SIAM J. Sci. Comput. 39, A333–A359 (2017)

    Article  MATH  Google Scholar 

  16. Chen, L., Zheng, B., Lin, G., Voulgarakis, N.: A two-level stochastic collocation method for semilinear elliptic equations with random coefficients. J. Comput. Appl. Math. 315, 195–207 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  17. Fishman, G.: Monte Carlo: Concepts, Algorithms, and Applications. Springer, New York (1996)

    Book  MATH  Google Scholar 

  18. Jiang, N., Layton, W.: An algorithm for fast calculation of flow ensembles. Int. J. Uncertain. Quantif. 4, 273–301 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  19. Jiang, N., Kaya, S., Layton, W.: Analysis of model variance for ensemble based turbulence modeling. Comput. Meth. Appl. Math. 15, 173–188 (2015)

    MathSciNet  MATH  Google Scholar 

  20. Jiang, N.: A higher order ensemble simulation algorithm for fluid flows. J. Sci. Comput. 64, 264–288 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  21. Gunzburger, M., Jiang, N., Wang, Z.: An efficient algorithm for simulating ensembles of parameterized flow problems, IMA J. Numer. Anal. (2018). https://doi.org/10.1093/imanum/dry029

  22. Gunzburger, M., Jiang, N., Wang, Z.: A second-order time-stepping scheme for simulating ensembles of parameterized flow problems. Comput. Methods Appl. Math. (2018). https://doi.org/10.1515/cmam-2017-0051

  23. Fiordilino, J.A.: A second order ensemble timestepping algorithm for natural convection. SIAM J. Numer. Anal. 56, 816–837 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  24. Fiordilino, J.A., Khankan, S.: Ensemble timestepping algorithms for natural convection. Int. J. Numer. Anal. Model. 15, 524–551 (2018)

    MathSciNet  MATH  Google Scholar 

  25. Mohebujjaman, M., Rebholz, L.: An efficient algorithm for computation of MHD flow ensembles. Comput. Methods Appl. Math. 17, 121–137 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  26. Fiordilino, J.A.: Ensemble timestepping algorithms for the heat equation with uncertain conductivity. Numer. Meth. Partial. Differ. Eqs. 34, 1901–1916 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  27. Luo, Y., Wang, Z.: An ensemble algorithm for numerical solutions to deterministic and random parabolic PDEs. SIAM J. Numer. Anal. 56, 859–876 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  28. Phoon, K.K., Huang, S.P., Quek, S.T.: Simulation of second-order processes using Karhunen-Loève expansion. Comput. Struct. 80, 1049–1060 (2002)

    Article  Google Scholar 

  29. Zhang, D., Lu, Z.: An efficient, high-order perturbation approach for flow in random porous media via Karhunen-Loève and polynomial expansions. J. Comput. Phys. 194, 773–794 (2004)

    Article  MATH  Google Scholar 

  30. Øksendal, B.: Stochastic Differential Equations: an Introduction with Applications. Springer, Berlin (2003)

    Book  MATH  Google Scholar 

  31. Brenner, S., Scott, R.: The Mathematical Theory of Finite Element Methods, vol. 15. Springer, New York (2007)

    Google Scholar 

  32. Hecht, F.: New development in freefem++. J. Numer. Math. 20, 251–265 (2012)

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors thank Professor William Layton for his helpful comments and suggestions.

Funding

Ning Li is supported by China Scholarship Council Grant 201707010004, Joseph Fiordilino is supported by the DoD SMART Scholarship and partially supported by NSF Grants CBET 1609120 and DMS 1522267, Xinlong Feng is supported by the NSF of China (No. 11671345, No. 11362021).

Author information

Authors and Affiliations

  1. College of Information and Management of Science, Henan Agricultural University, Zhengzhou, 450002, People’s Republic of China

    Ning Li

  2. Department of Mathematics, University of Pittsburgh, Pittsburgh, PA, 15260, USA

    Joseph Fiordilino

  3. College of Mathematics and Systems Science, Xinjiang University, Ürümqi, 830046, People’s Republic of China

    Xinlong Feng

Authors
  1. Ning Li
    View author publications

    You can also search for this author inPubMed Google Scholar

  2. Joseph Fiordilino
    View author publications

    You can also search for this author inPubMed Google Scholar

  3. Xinlong Feng
    View author publications

    You can also search for this author inPubMed Google Scholar

Corresponding author

Correspondence to Xinlong Feng.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, N., Fiordilino, J. & Feng, X. Ensemble Time-Stepping Algorithm for the Convection-Diffusion Equation with Random Diffusivity. J Sci Comput 79, 1271–1293 (2019). https://doi.org/10.1007/s10915-018-0890-8

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10915-018-0890-8

Keywords

Mathematics Subject Classification