Publicou 6 edições por ano
ISSN Imprimir: 2152-5080
ISSN On-line: 2152-5099
Indexed in
VARIANCE REDUCTION METHODS AND MULTILEVEL MONTE CARLO STRATEGY FOR ESTIMATING DENSITIES OF SOLUTIONS TO RANDOM SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS
RESUMO
This paper concerns the estimation of the density function of the solution to a random nonautonomous second-order linear differential equation with analytic data processes. In a recent contribution, we proposed to express the density function as an expectation, and we used a standard Monte Carlo algorithm to approximate the expectation. Although the algorithms worked satisfactorily for most test problems, some numerical challenges emerged for others, due to large statistical errors. In these situations, the convergence of the Monte Carlo simulation slows down severely, and noisy features plague the estimates. In this paper, we focus on computational aspects and propose several variance reduction methods to remedy these issues and speed up the convergence. First, we introduce a pathwise selection of the approximating processes which aims at controlling the variance of the estimator. Second, we propose a hybrid method, combining Monte Carlo and deterministic quadrature rules, to estimate the expectation. Third, we exploit the series expansions of the solutions to design a multilevel Monte Carlo estimator. The proposed methods are implemented and tested on several numerical examples to highlight the theoretical discussions and demonstrate the significant improvements achieved.
-
Soong, T.T., Random Differential Equations in Science and Engineering, New York: Academic Press, 1973.
-
Neckel, T. andRupp, F., Random Differential Equations in Scientific Computing, London: Walter de Gruyter, 2013.
-
Strand, J.L., Random Ordinary Differential Equations, J. Differ. Equations, 7:538-553, 1970.
-
Smith, R.C., Uncertainty Quantification. Theory, Implementation, and Application, New York: SIAM Computational Science & Engineering, 2014.
-
Cortes, J.-C., Jodar, L., Camacho, J., and Villafuerte, L., Random Airy Type Differential Equations: Mean Square Exact and Numerical Solutions, Comput. Math. Appl., 60, pp. 1237-1244,2010.
-
Calbo, G., Cortes, J.-C., and Jodar, L., Random Hermite Differential Equations: Mean Square Power Series Solutions and Statistical Properties, Appl. Math. Comput., 218(7):3654-3666,2011.
-
Calatayud, J., Cortes, J.-C., and Jornet, M., Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation, Mediterr. J. Math., 16(68):1-14, 2019.
-
Calatayud, J., Cortes, J.-C., and Jornet, M., Computational Uncertainty Quantification for Random Nonautonomous Second Order Linear Differential Equations via Adapted gPC: A Comparative Case Study with Random Frobenius Method and Monte Carlo Simulation, Open Math, 16:1651-1666,2018.
-
Xiu, D., Numerical Methods for Stochastic Computations. A Spectral Method Approach, Cambridge Texts in Applied Mathematics, New York: Princeton University Press, 2010.
-
Golmankhaneh, A.A., Porghoveh,N.A., andBaleanu, D., Mean Square Solutions of Second-Order Random Differential Equations by Using Homotopy Analysis Method, Rom. Rep. Phys., 65:350-362, 2013.
-
Khudair, A.R., Ameen, A.A., and Khalaf, S.L., Mean Square Solutions of Second-Order Random Differential Equations by Using Variational Iteration Method, Appl. Math. Sci., 5:2505-2519, 2011.
-
Khudair, A.R., Ameen, A.A., and Khalaf, S.L., Mean Square Solutions of Second-Order Random Differential Equations by Using Adomian Decomposition Method, Appl. Math. Sci., 5:2521-2535, 2011.
-
Khudair, A.K., Haddad, S.A.M., and Khalaf, S.L., Mean Square Solutions of Second-Order Random Differential Equations by Using the Differential Transformation Method, Open J. Appl. Sci, 6:287-297, 2016.
-
Villafuerte, L. and Chen-Charpentier, B.M., A Random Differential Transform Method: Theory and Applications, Appl. Math. Lett., 25:1490-1494,2012.
-
Calatayud, J., Cortes, J.-C., Jornet, M., and Villafuerte, L., Random Nonautonomous Second Order Linear Differential Equations: Mean Square Analytic Solutions and Their Statistical Properties, Adv. Differ. Equations, 2018(392):1-29, 2018.
-
Calatayud, J., Cortes, J.-C., and Jornet, M., Some Notes to Extend the Study on Random Nonautonomous Second Order Linear Differential Equations Appearing in Mathematical Modeling, Math. Comput. Appl, 23(4):76-89, 2018.
-
Cortes, J.-C., Navarro-Quiles, A., Romero, J.-V., and Rosello, M.-D., Probabilistic Solution of Random Autonomous First-Order Linear Systems of Ordinary Differential Equations, Rom. Rep. Phys., 68(4):1397-1406, 2016.
-
Calatayud, J., Chen-Charpentier, B.M., Cortes, J.-C., and Jornet, M., Combining Polynomial Chaos Expansions and the Random Variable Transformation Technique to Approximate the Density Function of Stochastic Problems, Including Some Epidemiological Models, Symmetry, 11(1):43-60, 2019.
-
El-Tawil, M.A., The Approximate Solutions of Some Stochastic Differential Equations Using Transformations, Appl. Math. Comput., 164(1):167-178, 2005.
-
Dorini, F.A., Cecconello, M.S., and Dorini, M.B., On the Logistic Equation Subject to Uncertainties in the Environmental Carrying Capacity and Initial Population Density, Commun. Nonlinear Sci. Numer. Simul. 33:160-173, 2016.
-
Santos, L.T., Dorini, F.A., and Cunha, M.C.C., The Probability Density Function to the Random Linear Transport Equation, Appl. Math. Comput., 216(5):1524-1530, 2010.
-
Calatayud, J., Cortes, J.-C., and Jornet, M., The Damped Pendulum Random Differential Equation: A Comprehensive Stochastic Analysis via the Computation of the Probability Density Function, Physica A, 512:261-279, 2018.
-
Hussein, A. and Selim, M.M., Solution of the Stochastic Radiative Transfer Equation with Rayleigh Scattering Using RVT Technique, Appl. Math. Comput, 218(13):7193-7203,2012.
-
Casaban, M.C., Cortes, J.-C., Romero, J.V., and Rosello, M.D., Solving Random Homogeneous Linear Second-Order Differential Equations: A Full Probabilistic Description, Mediterr. J. Math, 13(6):3817-3836, 2016.
-
Jornet, M., Calatayud, J., LeMaitre, O.P., and Cortes, J.-C., Second Order Linear Differential Equations with Analytic Uncertainties: Stochastic Analysis via the Computation of the Probability Density Function, J. Comput. Appl. Math., 374(112770):1-20, 2020.
-
Wolfram Research, Inc., Mathematica, Version 11.2, Champaign, IL: Wolfram Research, Inc., 2017.
-
Giles, M.B., Multilevel Monte Carlo Path Simulation, Oper. Res., 56(3):607-617,2008.
-
Giles, M.B., Multilevel Monte Carlo Methods, Acta Numer, 24:259-328, 2015.
-
Cliffe, A., Giles, M.B., Scheichl, R., and Teckentrup, A., Multilevel Monte Carlo Methods and Applications to Elliptic PDEs with Random Coefficients, Comput. Visual. Sci., 14(1):3-15, 2011.
-
Pisaroni, M., Nobile, F.B., and Leyland, P., A Continuation Multi Level Monte Carlo Method for Uncertainty Quantification in Compressible Inviscid Aerodynamics, Comput. Methods Appl. Mech. Eng., 326(C):20-50,2017.
-
Beskos, A., Jasra, A., Law, K., Tempone, R., and Zhou, Y., Multilevel Sequential Monte Carlo Samplers, Stoch. Proc. Appl, 127(5):1417-1440, 2017.
-
Sobol, I.M., Global Sensitivity Indices for Nonlinear Mathematical Models and Their Monte Carlo Estimates, Math. Comput. Simul., 55(1-3):271-280,2001.
-
Saltelli, A., Ratto, M., Andres, T., Campolongo, F., Cariboni, J., Gatelli, D., Saisana, M., and Tarantola, S., Global Sensitivity Analysis, The Primer, London: John Wiley & Sons, 2008.
-
Saltelli, A., Annoni, P., Azzini, I., Campolongo, F., Ratto, M., and Tarantola, S., Variance based Sensitivity Analysis of Model Output. Design and Estimator for the Total Sensitivity Index, Comput. Phys. Commun., 181(2):259-270, 2010.
-
Sobol', I.M., Tarantola, S., Gatelli, D., Kucherenko, S., and Mauntz, W., Estimating the Approximation Error when Fixing Unessential Factors in Global Sensitivity Analysis, Reliab. Eng. Syst. Saf., 92(7):957-960,2007.
-
Efron, B. and Tibshirani, R., An Introduction to the Bootstrap, New York: CRC Press, 1994.
-
Davison, A.C. and Hinkley, D.V., Bootstrap Methods and Their Application, London: Cambridge University Press, 1997.
-
Chicone, C., Ordinary Differential Equations with Applications, New York: Springer Science+Business Media, 2006.
-
Calatayud Julia, Jornet Marc, Extending the applicability of the RVT technique for the randomized radioactive decay chain model, The European Physical Journal Plus, 137, 3, 2022. Crossref