Maximum-Entropy Density Estimation for MRI Stochastic Surrogate Models

Stochastic spectral methods can generate accurate compact stochastic models for electromagnetic problems with material and geometric uncertainties. This letter presents an improved implementation of the maximum-entropy algorithm to compute the density function of an obtained generalized polynomial chaos expansion in magnetic resonance imaging (MRI) applications. Instead of using statistical moments, we show that the expectations of some orthonormal polynomials can be better constraints for the optimization flow. The proposed algorithm is coupled with a finite element-boundary element method (FEM-BEM) domain decomposition field solver to obtain a robust computational prototyping for MRI problems with low- and high-dimensional uncertainties.