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
    Analysis of a stabilized penalty-free Nitsche method for the Brinkman, Stokes, and Darcy problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Blank, Laura; Caiazzo, Alfonso; Chouly, Franz; Lozinski, Alexei; Mura, Joaquin
    In this paper we study the Brinkman model as a unified framework to allow the transition between the Darcy and the Stokes problems. We propose an unconditionally stable low-order finite element approach, which is robust with respect to the whole range of physical parameters, and is based on the combination of stabilized equal-order finite elements with a non-symmetric penalty free Nitsche method for the weak imposition of essential boundary conditions. In particular, we study the properties of the penalty-free Nitsche formulation for the Brinkman setting, extending a recently reported analysis for the case of incompressible elasticity (T. Boiveau & E. Burman. IMA J. Numer. Anal. 36 (2016), no.2, 770-795). Focusing on the two-dimensional case, we obtain optimal a priori error estimates in a mesh-dependent norm, which, converging to natural norms in the cases of Stokes or Darcy flows, allows to extend the results also to these limits. Moreover, we show that, in order to obtain robust estimates also in the Darcy limit, the formulation shall be equipped with a Grad-Div stabilization and an additional stabilization to control the discontinuities of the normal velocity along the boundary. The conclusions of the analysis are supported by numerical simulations.