Bayesian inverse Navier-Stokes problems: joint flow field reconstruction and parameter learning
- URL: http://arxiv.org/abs/2406.18464v3
- Date: Tue, 10 Dec 2024 23:33:51 GMT
- Title: Bayesian inverse Navier-Stokes problems: joint flow field reconstruction and parameter learning
- Authors: Alexandros Kontogiannis, Scott V. Elgersma, Andrew J. Sederman, Matthew P. Juniper,
- Abstract summary: We formulate and solve a Bayesian inverse Navier-Stokes (N-S) problem that assimilates velocimetry data.
We learn the unknown N-S parameters, including the boundary position.
We then use this method to reconstruct magnetic resonance velocimetry data of a 3D steady laminar flow.
- Score: 44.62264781248436
- License:
- Abstract: We formulate and solve a Bayesian inverse Navier-Stokes (N-S) problem that assimilates velocimetry data in order to jointly reconstruct a 3D flow field and learn the unknown N-S parameters, including the boundary position. By hardwiring a generalised N-S problem, and regularising its unknown parameters using Gaussian prior distributions, we learn the most likely parameters in a collapsed search space. The most likely flow field reconstruction is then the N-S solution that corresponds to the learned parameters. We develop the method in the variational setting and use a stabilised Nitsche weak form of the N-S problem that permits the control of all N-S parameters. To regularise the inferred the geometry, we use a viscous signed distance field (vSDF) as an auxiliary variable, which is given as the solution of a viscous Eikonal boundary value problem. We devise an algorithm that solves this inverse problem, and numerically implement it using an adjoint-consistent stabilised cut-cell finite element method. We then use this method to reconstruct magnetic resonance velocimetry (flow-MRI) data of a 3D steady laminar flow through a physical model of an aortic arch for two different Reynolds numbers and signal-to-noise ratio (SNR) levels (low/high). We find that the method can accurately i) reconstruct the low SNR data by filtering out the noise/artefacts and recovering flow features that are obscured by noise, and ii) reproduce the high SNR data without overfitting. Although the framework that we develop applies to 3D steady laminar flows in complex geometries, it readily extends to time-dependent laminar and Reynolds-averaged turbulent flows, as well as non-Newtonian (e.g. viscoelastic) fluids.
Related papers
- Bayesian inference of mean velocity fields and turbulence models from flow MRI [0.0]
We devise an algorithm that learns the most likely parameters of an effective viscosity model, and estimates their uncertainties, from mean flow data of a turbulent flow.
The algorithm successfully reconstructs the mean flow field and learns the most likely turbulence model parameters without overfitting.
arXiv Detail & Related papers (2024-12-15T18:07:36Z) - Geometry-aware PINNs for Turbulent Flow Prediction [3.193631734016982]
A novel geometry aware parametric PINN surrogate model has been created, which can predict flow fields for NACA 4 digit airfoils in turbulent conditions.
The models have been validated with unknown NACA airfoils at unseen Reynolds numbers.
arXiv Detail & Related papers (2024-12-02T20:30:53Z) - Learning rheological parameters of non-Newtonian fluids from velocimetry data [46.2482873419289]
We devise an algorithm that learns the most likely Carreau parameters of a shear-thinning fluid.
We show that the algorithm can successfully reconstruct the flow field by learning the most likely Carreau parameters.
arXiv Detail & Related papers (2024-08-05T16:27:38Z) - Decomposed Diffusion Sampler for Accelerating Large-Scale Inverse
Problems [64.29491112653905]
We propose a novel and efficient diffusion sampling strategy that synergistically combines the diffusion sampling and Krylov subspace methods.
Specifically, we prove that if tangent space at a denoised sample by Tweedie's formula forms a Krylov subspace, then the CG with the denoised data ensures the data consistency update to remain in the tangent space.
Our proposed method achieves more than 80 times faster inference time than the previous state-of-the-art method.
arXiv Detail & Related papers (2023-03-10T07:42:49Z) - Orthogonal Matrix Retrieval with Spatial Consensus for 3D Unknown-View
Tomography [58.60249163402822]
Unknown-view tomography (UVT) reconstructs a 3D density map from its 2D projections at unknown, random orientations.
The proposed OMR is more robust and performs significantly better than the previous state-of-the-art OMR approach.
arXiv Detail & Related papers (2022-07-06T21:40:59Z) - Physics-informed compressed sensing for PC-MRI: an inverse Navier-Stokes
problem [78.20667552233989]
We formulate a physics-informed compressed sensing (PICS) method for the reconstruction of velocity fields from noisy and sparse magnetic resonance signals.
We find that the method is capable of reconstructing and segmenting the velocity fields from sparsely-sampled signals.
arXiv Detail & Related papers (2022-07-04T14:51:59Z) - Simultaneous boundary shape estimation and velocity field de-noising in
Magnetic Resonance Velocimetry using Physics-informed Neural Networks [70.7321040534471]
Magnetic resonance velocimetry (MRV) is a non-invasive technique widely used in medicine and engineering to measure the velocity field of a fluid.
Previous studies have required the shape of the boundary (for example, a blood vessel) to be known a priori.
We present a physics-informed neural network that instead uses the noisy MRV data alone to infer the most likely boundary shape and de-noised velocity field.
arXiv Detail & Related papers (2021-07-16T12:56:09Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.