# Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format

- Handle URI:
- http://hdl.handle.net/10754/347255
- Title:
- Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format
- Authors:
- Abstract:
- We apply the Tensor Train (TT) decomposition to construct the tensor product Polynomial Chaos Expansion (PCE) of a random field, to solve the stochastic elliptic diffusion PDE with the stochastic Galerkin discretization, and to compute some quantities of interest (mean, variance, exceedance probabilities). We assume that the random diffusion coefficient is given as a smooth transformation of a Gaussian random field. In this case, the PCE is delivered by a complicated formula, which lacks an analytic TT representation. To construct its TT approximation numerically, we develop the new block TT cross algorithm, a method that computes the whole TT decomposition from a few evaluations of the PCE formula. The new method is conceptually similar to the adaptive cross approximation in the TT format, but is more efficient when several tensors must be stored in the same TT representation, which is the case for the PCE. Besides, we demonstrate how to assemble the stochastic Galerkin matrix and to compute the solution of the elliptic equation and its post-processing, staying in the TT format. We compare our technique with the traditional sparse polynomial chaos and the Monte Carlo approaches. In the tensor product polynomial chaos, the polynomial degree is bounded for each random variable independently. This provides higher accuracy than the sparse polynomial set or the Monte Carlo method, but the cardinality of the tensor product set grows exponentially with the number of random variables. However, when the PCE coefficients are implicitly approximated in the TT format, the computations with the full tensor product polynomial set become possible. In the numerical experiments, we confirm that the new methodology is competitive in a wide range of parameters, especially where high accuracy and high polynomial degrees are required.
- KAUST Department:
- Citation:
- @article{dolgov2015polynomial, title={Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format}, author={Dolgov, Sergey and Khoromskij, Boris N and Litvinenko, Alexander and Matthies, Hermann G}, journal={arXiv preprint arXiv:1503.03210}, year={2015} }
- Publisher:
- Journal:
- Issue Date:
- 11-Mar-2015
- Type:
- Article
- Description:
- Submitted to SIAM/ASA Journal on Uncertainty Quantification (JUQ)
- Sponsors:
- SRI Uncertainty Quantification and Extreme Computing Research Centers, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
- Additional Links:
- http://arxiv.org/abs/1503.03210

- Appears in Collections:
- Articles

# Full metadata record

DC Field | Value | Language |
---|---|---|

dc.contributor.author | Dolgov, S. | en |

dc.contributor.author | Khoromskij, B.N. | en |

dc.contributor.author | Litvinenko, Alexander | en |

dc.contributor.author | Matthies, H.G. | en |

dc.date.accessioned | 2015-03-29T05:56:27Z | en |

dc.date.available | 2015-03-29T05:56:27Z | en |

dc.date.issued | 2015-03-11 | en |

dc.identifier.citation | @article{dolgov2015polynomial, title={Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format}, author={Dolgov, Sergey and Khoromskij, Boris N and Litvinenko, Alexander and Matthies, Hermann G}, journal={arXiv preprint arXiv:1503.03210}, year={2015} } | en |

dc.identifier.uri | http://hdl.handle.net/10754/347255 | en |

dc.description | Submitted to SIAM/ASA Journal on Uncertainty Quantification (JUQ) | en |

dc.description.abstract | We apply the Tensor Train (TT) decomposition to construct the tensor product Polynomial Chaos Expansion (PCE) of a random field, to solve the stochastic elliptic diffusion PDE with the stochastic Galerkin discretization, and to compute some quantities of interest (mean, variance, exceedance probabilities). We assume that the random diffusion coefficient is given as a smooth transformation of a Gaussian random field. In this case, the PCE is delivered by a complicated formula, which lacks an analytic TT representation. To construct its TT approximation numerically, we develop the new block TT cross algorithm, a method that computes the whole TT decomposition from a few evaluations of the PCE formula. The new method is conceptually similar to the adaptive cross approximation in the TT format, but is more efficient when several tensors must be stored in the same TT representation, which is the case for the PCE. Besides, we demonstrate how to assemble the stochastic Galerkin matrix and to compute the solution of the elliptic equation and its post-processing, staying in the TT format. We compare our technique with the traditional sparse polynomial chaos and the Monte Carlo approaches. In the tensor product polynomial chaos, the polynomial degree is bounded for each random variable independently. This provides higher accuracy than the sparse polynomial set or the Monte Carlo method, but the cardinality of the tensor product set grows exponentially with the number of random variables. However, when the PCE coefficients are implicitly approximated in the TT format, the computations with the full tensor product polynomial set become possible. In the numerical experiments, we confirm that the new methodology is competitive in a wide range of parameters, especially where high accuracy and high polynomial degrees are required. | en |

dc.description.sponsorship | SRI Uncertainty Quantification and Extreme Computing Research Centers, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia | en |

dc.language.iso | en | en |

dc.publisher | arXiv, http://arxiv.org/abs/1503.03210 | en |

dc.relation.url | http://arxiv.org/abs/1503.03210 | en |

dc.subject | uncertainty quantification | en |

dc.subject | polynomial chaos expansion | en |

dc.subject | Karhunen-Lo`eve expansion | en |

dc.subject | stochastic Galerkin | en |

dc.subject | tensor product methods | en |

dc.subject | tensor train format | en |

dc.subject | adaptive cross approximation | en |

dc.subject | block cross | en |

dc.title | en | |

dc.type | Article | en |

dc.contributor.department | Max-Planck-Institut für Dynamik komplexer technischer Systeme, Sandtorstr. 1, 39106 Magdeburg, Germany | en |

dc.contributor.department | Max-Planck-Institut für Mathematik in den Naturwissenschaften, Inselstraße 22, 04103 Leipzig, Germany. | en |

dc.contributor.department | King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia. | en |

dc.contributor.department | Institute for Scientific Computing, Technische Universit ̈at Braunschweig, Hans- Sommerstr. 65, Brunswick, Germany. | en |

dc.identifier.journal | arXiv | en |

dc.eprint.version | Pre-print, submitted to SIAM/ASA Journal on Uncertainty Quantification (JUQ) | en |

All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.