Search Results

Now showing 1 - 2 of 2
  • Item
    Reconstruction of flow domain boundaries from velocity data via multi-step optimization of distributed resistance
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2022) Pártl, Ondřej; Wilbrandt, Ulrich; Mura, Joaquín; Caiazzo, Alfonso
    We reconstruct the unknown shape of a flow domain using partially available internal velocity measurements. This inverse problem is motivated by applications in cardiovascular imaging where motion-sensitive protocols, such as phase-contrast MRI, can be used to recover three-dimensional velocity fields inside blood vessels. In this context, the information about the domain shape serves to quantify the severity of pathological conditions, such as vessel obstructions. We consider a flow modeled by a linear Brinkman problem with a fictitious resistance accounting for the presence of additional boundaries. To reconstruct these boundaries, we employ a multi-step gradient-based variational method to compute a resistance that minimizes the difference between the computed flow velocity and the available data. Afterward, we apply different post-processing steps to reconstruct the shape of the internal boundaries. To limit the overall computational cost, we use a stabilized equal-order finite element method. We prove the stability and the well-posedness of the considered optimization problem. We validate our method on three-dimensional examples based on synthetic velocity data and using realistic geometries obtained from cardiovascular imaging.
  • Item
    An assessment of solvers for algebraically stabilized discretizations of convection-diffusion-reaction equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021) Jha, Abhinav; Pártl, Ondřej; Ahmed, Naveed; Kuzmin, Dmitri
    We consider flux-corrected finite element discretizations of 3D convection-dominated transport problems and assess the computational efficiency of algorithms based on such approximations. The methods under investigation include flux-corrected transport schemes and monolithic limiters. We discretize in space using a continuous Galerkin method and P1 or Q1 finite elements. Time integration is performed using the Crank-Nicolson method or an explicit strong stability preserving Runge-Kutta method. Nonlinear systems are solved using a fixed-point iteration method, which requires solution of large linear systems at each iteration or time step. The great variety of options in the choice of discretization methods and solver components calls for a dedicated comparative study of existing approaches. To perform such a study, we define new 3D test problems for time dependent and stationary convection-diffusion-reaction equations. The results of our numerical experiments illustrate how the limiting technique, time discretization and solver impact on the overall performance.