libMoM: a library for stochastic simulations in engineering using statistical moments

Stochastic simulations are becoming increasingly important in numerous engineering applications. The solution to the governing equations are complicated due to the high-dimensional spaces and the presence of randomness. In this paper we present libMoM (http://libmom.sourceforge.net), a software library to solve various types of Stochastic Differential Equations (SDE) as well as estimate statistical distributions from the moments. The library provides a suite of tools to solve various SDEs using the method of moments (MoM) as well as estimate statistical distributions from the moments using moment matching algorithms. For a large class of problems, MoM provide efficient solutions compared with other stochastic simulation techniques such as Monte Carlo (MC). In the physical sciences, the moments of the distribution are usually the primary quantities of interest. The library enables the solution of moment equations derived from a variety of SDEs, with closure using non-standard Gaussian quadrature. In engineering risk assessment and decision making, statistical distributions are required. The library implements tools for fitting the Generalized Lambda Distribution (GLD) with the given moments. The objectives of this paper are (1) to briefly outline the theory behind moment methods for solving SDEs/estimation of statistical distributions; (2) describe the organization of the software and user interfaces; (3) discuss use of standard software engineering tools for regression testing, aid collaboration, distribution and further development. A number of representative examples of the use of libMoM in various engineering applications are presented and future areas of research are discussed.

[1]  Marshall Freimer,et al.  a study of the generalized tukey lambda family , 1988 .

[2]  Yunjun Xu,et al.  A quadrature-based method of moments for nonlinear filtering , 2009, Autom..

[3]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[4]  Marc Massot,et al.  Numerical simulation of spray coalescence in an Eulerian framework: Direct quadrature method of moments and multi-fluid method , 2007, J. Comput. Phys..

[5]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[6]  Robert McGraw,et al.  Numerical advection of correlated tracers: preserving particle size/composition moment sequences during transport of aerosol mixtures , 2007 .

[7]  Ofodike A. Ezekoye,et al.  Smoke buildup and light scattering in a cylindrical cavity above a uniform flow , 2005 .

[8]  Jasper V. Stokman,et al.  Orthogonal Polynomials of Several Variables , 2001, J. Approx. Theory.

[9]  S. Pratsinis.  Simultaneous nucleation, condensation, and coagulation in aerosol reactors , 1988 .

[10]  Frances Y. Kuo,et al.  Remark on algorithm 659: Implementing Sobol's quasirandom sequence generator , 2003, TOMS.

[11]  Jean-Pierre Gaspard,et al.  Continued-fraction technique for tight-binding systems. A generalized-moments method , 1982 .

[12]  M. Frenklach.  Method of moments with interpolative closure , 2002 .

[13]  A. Lakhany,et al.  Estimating the parameters of the generalized lambda distribution , 2000 .

[14]  I. M. Sobol,et al.  Random and quasirandom sequences: Numerical estimates of uniformity of distribution , 1993 .

[15]  Robert McGraw,et al.  Description of Aerosol Dynamics by the Quadrature Method of Moments , 1997 .

[16]  Doraiswami Ramkrishna,et al.  Population Balances: Theory and Applications to Particulate Systems in Engineering , 2000 .

[17]  Yuan Xu,et al.  Orthogonal Polynomials of Several Variables: Subject index , 2001 .

[18]  Yuan Xu,et al.  Orthogonal Polynomials of Several Variables , 2014, 1701.02709.

[19]  Prakash Vedula,et al.  Direct Quadrature Method of Moments Solution of the Fokker-Planck Equation for Stochastic Processes in Aeroelasticity , 2008 .

[20]  D. E. Rosner,et al.  Bivariate Extension of the Quadrature Method of Moments for Modeling Simultaneous Coagulation and Sintering of Particle Populations. , 2001, Journal of colloid and interface science.

[21]  Ofodike A. Ezekoye,et al.  Evaluation of the 1-point quadrature approximation in QMOM for combined aerosol growth laws , 2003 .

[22]  Carlos A. Dorao,et al.  The quadrature method of moments and its relationship with the method of weighted residuals , 2006 .

[23]  Frances Y. Kuo,et al.  Constructing Sobol Sequences with Better Two-Dimensional Projections , 2008, SIAM J. Sci. Comput..

[24]  Ofodike A. Ezekoye,et al.  Treatment of size-dependent aerosol transport processes using quadrature based moment methods , 2006 .

[25]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[26]  William H. Press,et al.  Numerical recipes in Fortran 77 : the art of scientificcomputing. , 1992 .

[27]  Ofodike A. Ezekoye,et al.  Treatment of design fire uncertainty using Quadrature Method of Moments , 2008 .

[28]  Rodney O. Fox,et al.  Optimal Moment Sets for Multivariate Direct Quadrature Method of Moments , 2009 .

[29]  Rochan Raj Upadhyay,et al.  Simulation of population balance equations using quadrature based moment methods , 2006 .

[30]  D. Sherrington.  Stochastic Processes in Physics and Chemistry , 1983 .

[31]  Prakash Vedula,et al.  Direct quadrature method of moments solution of the Fokker–Planck equation , 2008 .

[32]  R. McGraw,et al.  Method of Moments Solutions for a Laminar Flow Aerosol Reactor Model , 2001 .

[33]  Roy G. Gordon,et al.  Error Bounds in Equilibrium Statistical Mechanics , 1968 .

[34]  William H. Press,et al.  Numerical Recipes in Fortran 77: The Art of Scientific Computing 2nd Editionn - Volume 1 of Fortran Numerical Recipes , 1992 .

[35]  S. Friedlander.  Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics , 2000 .

[36]  Thomas E. Wehrly.  Fitting Statistical Distributions: The Generalized Lambda Distribution and Generalized Bootstrap Methods , 2002, Technometrics.

[37]  Robert McGraw,et al.  Representation of generally mixed multivariate aerosols by the quadrature method of moments: I. Statistical foundation , 2004 .

[38]  Daniele Marchisio,et al.  Solution of population balance equations using the direct quadrature method of moments , 2005 .

[39]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[40]  S. Katz,et al.  Some problems in particle technology: A statistical mechanical formulation , 1964 .

[41]  W. J. Thron,et al.  Encyclopedia of Mathematics and its Applications. , 1982 .

[42]  Robert McGraw,et al.  Chemically resolved aerosol dynamics for internal mixtures by the quadrature method of moments , 2003 .