TY - JOUR
T1 - Autonomous uncertainty quantification for discontinuous models using multivariate Padé approximations
AU - Ahlfeld, Richard
AU - Montomoli, Francesco
AU - Carnevale, Mauro
AU - Salvadori, Simone
PY - 2018/4/1
Y1 - 2018/4/1
N2 - Problems in turbomachinery computational fluid dynamics (CFD) are often characterized by nonlinear and discontinuous responses. Ensuring the reliability of uncertainty quantification (UQ) codes in such conditions, in an autonomous way, is challenging. In this work, we suggest a new approach that combines three state-of-the-art methods: multivariate Padé approximations, optimal quadrature subsampling (OQS), and statistical learning. Its main component is the generalized least-squares multivariate Padé- Legendre (PL) approximation. PL approximations are globally fitted rational functions that can accurately describe discontinuous nonlinear behavior. They need fewer model evaluations than local or adaptive methods and do not cause the Gibbs phenomenon like continuous polynomial chaos methods. A series of modifications of the Padé algorithm allows us to apply it to arbitrary input points instead of optimal quadrature locations. This property is particularly useful for industrial applications, where a database of CFD runs is already available, but not in optimal parameter locations. One drawback of the PL approximation is that it is nontrivial to ensure reliability. To improve stability, we suggest to couple it with OQS. Our reasoning is that least-squares errors, caused by an ill-conditioned design matrix, are the main source of error. Finally, we use statistical learning methods to check smoothness and convergence. The resulting method is shown to efficiently and correctly fit thousands of partly discontinuous response surfaces for an industrial film cooling and shock interaction problem using only nine CFD simulations.
AB - Problems in turbomachinery computational fluid dynamics (CFD) are often characterized by nonlinear and discontinuous responses. Ensuring the reliability of uncertainty quantification (UQ) codes in such conditions, in an autonomous way, is challenging. In this work, we suggest a new approach that combines three state-of-the-art methods: multivariate Padé approximations, optimal quadrature subsampling (OQS), and statistical learning. Its main component is the generalized least-squares multivariate Padé- Legendre (PL) approximation. PL approximations are globally fitted rational functions that can accurately describe discontinuous nonlinear behavior. They need fewer model evaluations than local or adaptive methods and do not cause the Gibbs phenomenon like continuous polynomial chaos methods. A series of modifications of the Padé algorithm allows us to apply it to arbitrary input points instead of optimal quadrature locations. This property is particularly useful for industrial applications, where a database of CFD runs is already available, but not in optimal parameter locations. One drawback of the PL approximation is that it is nontrivial to ensure reliability. To improve stability, we suggest to couple it with OQS. Our reasoning is that least-squares errors, caused by an ill-conditioned design matrix, are the main source of error. Finally, we use statistical learning methods to check smoothness and convergence. The resulting method is shown to efficiently and correctly fit thousands of partly discontinuous response surfaces for an industrial film cooling and shock interaction problem using only nine CFD simulations.
UR - http://www.scopus.com/inward/record.url?scp=85041295544&partnerID=8YFLogxK
U2 - 10.1115/1.4038826
DO - 10.1115/1.4038826
M3 - Article
AN - SCOPUS:85041295544
VL - 140
JO - Journal of Turbomachinery: Transactions of the ASME
JF - Journal of Turbomachinery: Transactions of the ASME
SN - 0889-504X
IS - 4
M1 - 041004
ER -