Numerical solution of the Stratonovich- and Ito-Euler equations: Application to the stochastic piston problem

We consider a piston with a velocity perturbed by Brownian motion moving into a straight tube filled with a perfect gas at rest. The shock generated ahead of the piston can be located by solving the one-dimensional Euler equations driven by white noise using the Stratonovich or Ito formulations. We approximate the Brownian motion with its spectral truncation and subsequently apply stochastic collocation using either sparse grid or the quasi-Monte Carlo (QMC) method. In particular, we first transform the Euler equations with an unsteady stochastic boundary into stochastic Euler equations over a fixed domain with a time-dependent stochastic source term. We then solve the transformed equations by splitting them up into two parts, i.e., a 'deterministic part' and a 'stochastic part'. Numerical results verify the Stratonovich-Euler and Ito-Euler models against stochastic perturbation results, and demonstrate the efficiency of sparse grid and QMC for small and large random piston motions, respectively. The variance of shock location of the piston grows cubically in the case of white noise in contrast to colored noise reported in [1], where the variance of shock location grows quadratically with time for short times and linearly for longer times.

[1]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[2]  Henryk Wozniakowski,et al.  When Are Quasi-Monte Carlo Algorithms Efficient for High Dimensional Integrals? , 1998, J. Complex..

[3]  H. Holden,et al.  Conservation laws with a random source , 1997 .

[4]  E. Wong,et al.  On the Convergence of Ordinary Integrals to Stochastic Integrals , 1965 .

[5]  Bernt Øksendal,et al.  Stochastic differential equations (3rd ed.): an introduction with applications , 1992 .

[6]  O. L. Maître,et al.  Uncertainty propagation in CFD using polynomial chaos decomposition , 2006 .

[7]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[8]  Habib N. Najm,et al.  Uncertainty Quantification and Polynomial Chaos Techniques in Computational Fluid Dynamics , 2009 .

[9]  Frances Y. Kuo,et al.  Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications , 2011, J. Comput. Phys..

[10]  Zhongqiang Zhang,et al.  A Multistage Wiener Chaos Expansion Method for Stochastic Advection-Diffusion-Reaction Equations , 2012, SIAM J. Sci. Comput..

[11]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[12]  R. Lindsay,et al.  Elements of gasdynamics , 1957 .

[13]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[14]  J. M. Sek,et al.  On the L2-discrepancy for anchored boxes , 1998 .

[15]  W. J. Whiten,et al.  Computational investigations of low-discrepancy sequences , 1997, TOMS.

[16]  MatoušekJiří.  On the L2-discrepancy for anchored boxes , 1998 .

[17]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[18]  E. Wong,et al.  ON THE RELATION BETWEEN ORDINARY AND STOCHASTIC DIFFERENTIAL EQUATIONS , 1965 .

[19]  BabuskaIvo,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007 .

[20]  桜井 健郎,et al.  H. W.Liepmann and A.Roshko: Elements of Gasdynamics John Wiley and Sons, Inc., New York, 1957, 439頁, 15×23cm, $11.00. , 1957 .

[21]  H. Sussmann.  On the Gap Between Deterministic and Stochastic Ordinary Differential Equations , 1978 .

[22]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[23]  G E Karniadakis,et al.  The stochastic piston problem. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Fabio Nobile,et al.  A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[25]  Jirí Matousek,et al.  On the L2-Discrepancy for Anchored Boxes , 1998, J. Complex..

[26]  A. Genz,et al.  Fully symmetric interpolatory rules for multiple integrals over infinite regions with Gaussian weight , 1996 .