Adjoint-accelerated Bayesian Inference for joint reconstruction and segmentation of Flow-MRI images

We formulate and solve a generalized inverse Navier–Stokes boundary value problem for velocity field reconstruction and simultaneous boundary segmentation of noisy Flow-MRI velocity images. We use a Bayesian framework that combines CFD, Gaussian processes, adjoint methods, and shape optimization in a unified and rigorous manner.
With this framework, we find the velocity field and flow boundaries (i.e. the digital twin) that are most likely to have produced a given noisy image. We also calculate the posterior covariances of the unknown parameters and thereby deduce the uncertainty in the reconstructed flow. First, we verify this method on synthetic noisy images of flows. Then we apply it to experimental phase contrast magnetic resonance (PC-MRI) images of an axisymmetric flow at low and high SNRs. We show that this method successfully reconstructs and segments the low SNR images, producing noiseless velocity fields that match the high SNR images, using 30 times less data.
This framework also provides additional flow information, such as the pressure field and wall shear stress, accurately and with known error bounds. We demonstrate this further on a 3-D in-vitro flow through a 3D-printed aorta and 3-D in-vivo flow through a carotid artery.