1 / 30

* Institute of Thermal Technology, Silesian University of Technology, Poland

Application of the Proper Orthogonal Decomposition for Bayesian Estimation of flow parameters in porous medium. Zbigniew Buliński * , Helcio R.B. Orlande ** , Tomasz Krysiński * , Łukasz Ziółkowski * , Sebastian Werle **.

rleak
Download Presentation

* Institute of Thermal Technology, Silesian University of Technology, Poland

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Application of the Proper Orthogonal Decomposition for Bayesian Estimation of flow parameters in porous medium Zbigniew Buliński*, Helcio R.B. Orlande**, Tomasz Krysiński*, Łukasz Ziółkowski*, Sebastian Werle** * Institute of Thermal Technology, Silesian University of Technology, Poland ** Department of Mechanical Engineering, Politecnica/COPPE Federal University of Rio de Janeiro – UFRJ, Rio de Janeiro, Brasil New Trends in Parameter Identification for Mathematical Models 2017 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  2. Motivation • According to EU regulations biomass is treated as a renewable source of Energy which might be used directly (burning) or indirectly (gasification), • The European Environment Agency estimates that around 235 Mtoe/year of biomass could be made available in EU by 2020 without harming the environment (toe ≈ 42GJ), • Effective utilisation of this source of energy needs reliable mathematical model of heat and fluid flow, • In some applications mathematical model of biomass burning or gasifiaction may be based on the porous medium formulation (fixed bed), • Porous medium (material) is a solid material (solid matrix) containing interconnected voids which filled with a fluid – two phase system, • Usually model parameters are unknown. 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  3. Mathematical modelling of porous material Direct modelling • Detailed geometrical model of the porous material • Huge computational cost • Detailed solution field Velocity field Mesh size ~2.2 M cells Domain size: 3x3x3 mm 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  4. Mathematical modeling of porous material Volume (sometimes also time) averaged flow equations • Darcy Law • Forchheimer’s equation (for isotropicporous medium) • Hsu and Cheng, 1990, Vafai and Kim, 1990 refinement 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  5. Mathematical modeling of porous material Volume (sometimes also time) averaged flow equations • Darcy Law • Forchheimer’s equation (for isotropicporous medium) • Hsu and Cheng, 1990, Vafai and Kim, 1990 refinement 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  6. Problem formulation Vector of unknown parameters Measured data 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  7. Statistical inverse problem formulation Main features of the Statistical inverse problem formulation • All unknown quantities (estimated parameters) and input data (measurements) are treated as random variables, • Knowledge on the realisation of all these quantities is given in the form of the probability distribution, • Solution of the considered problem is given in a form of probability distribution of unknown quantities, 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  8. Bayesian inverse problem formulation Bayes formula Likelihood distribution 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  9. Bayesian inverse problem formulation Prior distribution Gaussian prior Uniform prior 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  10. Bayesian inverse problem formulation Posterior distribution of unknown parameters No explicit analytical formulation of the distribution is known, obtaining distribution moments needs numerical integration usually by sampling. 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  11. Metropolis-Hastings algorithm 1. Draw a Candidate PointK+ from a known proposal distribution p(K+,Ki-1), 2. Calculate the acceptance factor: 3. Generate a random number R that is uniformly distributed on (0,1), 4. Accept candidate if R < a, set Ki= K+,otherwise set Ki = Ki-1, 5. Return to step 1. 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  12. Metropolis-Hastings algorithm Proposal distribution – random walk 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  13. Posterior distribution of unknown parameters Exact values 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  14. Posterior distribution of unknown parameters Porosity Pore diameter Exact values Reduced Order Model Computation time ~72h 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  15. Proper Orthogonal Decomposition (POD) Single snapshot yi is discretized response (image) of the field under investigationfor given spatial coordinates and time System Snapshot Matrix of snapshots 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  16. Proper Orthogonal Decomposition (POD) Steady state flow problem – a snapshot represents discretised pressure field inside computational domain Matrix of snapshots Snapshot System 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  17. Proper Orthogonal Decomposition (POD) Steady state flow problem – a snapshot represents discretised pressure field inside computational domain System Snapshot 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  18. Proper Orthogonal Decomposition (POD) Retrieve structure (anatomy) of the snapshot matrix by SVD Matrix of amplitudes Singular vectors resolving the response space Product describing influence of input space May be interpreted as a variables separation 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  19. Proper Orthogonal Decomposition (POD) Truncate the input space basis 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  20. Proper Orthogonal Decomposition (POD) Approximation of the amplitude matrix with a set of basis functions Finally, response of the system for any set of input parameters K can be calculated with trained POD-RBF network 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  21. Proper Orthogonal Decomposition (POD) Number of snapshots – full factorization Snapshots number 11x11 Snapshots number 21x21 Snapshots number 31x31 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  22. Proper Orthogonal Decomposition (POD) Number of snapshots – full factorization 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  23. Proper Orthogonal Decomposition (POD) Type of approximation function MQD IMQD Tspline Gauss 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  24. Reduced Order Model – POD Pressure field - exact model Pressure field – reduced POD model 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  25. Reduced Order Model – POD 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  26. Approximation Error Model (AEM) Posterior distribution Modified likelihood distribution Prior distribution 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  27. Posterior distribution of unknown parameters (AEM) Exact values 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  28. Posterior distribution of unknown parameters (AEM) Porosity Pore diameter Exact values Computation time 5 min 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  29. Summary • Bayesian inverse methodology was developed to estimate porosity and pore diameter of the porous material based on the pressure drop measurements, • Reduced Order Model based on the Proper Order Decomposition was developed to reduce computational time • Approximation Error Model (AEM) was used to incorporate approximation errors into the inference procedure giving promising results, • Developed procedure is going to be applied to real-life problem of estimation of flow parameters for the biomass. 30.10 – 03.11.2017, Rio de Janeiro, Brasil

  30. The financial support of The National Science Centre, Poland, within the grant UMO-2014/15/D/ST8/02767  is gratefully acknowledged herewith. Investigation of heat and mass transfer in porous media with chemical reactions withuseof statistical inverse methodology zbigniew.bulinski@polsl.pl 30.10 – 03.11.2017, Rio de Janeiro, Brasil

More Related