/pde-and-uq

Application examples of UQ in PDEs

Primary LanguageJupyter Notebook

Partial differential equation with stochastic inputs

I am working on reproducing the results from this paper with Proteus code. One objective is to learn Proteus which is DG-FEM code, another is to understand the paper better.

Structure

The repository is split onto several different directories

  • gpc - contains example generalized polynomial chaos setting for the Burgers 1D problegpc - contains example generalized polynomial chaos setting for the Burgers 1D problem.
  • graphs - sample code that I used to visualise some of the results
  • high_dimension - example setup for an exponential correlation function. The point is here to prove sensitivity of high dimensional IC on the forcing term.
  • notebooks - iPython notebooks which serve as documentation
  • riemann - Riemann problem examples

Current issues

  • The multi-dimensional case with forcing doesn't seem to work. I can only see small effect of the forcing.
  • I was unable to rerun the kinetic case in 3D. The third dimension would correspond the the forcing.