@misc{11953,
  author = {Harish Narayanan},
  title = {A Goal-Oriented Error-Controlled Solver for Biomedical Flows},
  abstract = {Accurate modelling of biomedical flows is of substantial interest because of its immediate applicability to a variety of fields, including the design of medical devices, the planning of surgical procedures and most fundamentally, the scientific understanding of healthy and diseased biological function. The complexity in this problem arises from the fact that physically-interesting flows occur in complex, patient-specific geometries and often occur in conjunction with other physics such as the finite deformation of the vascular walls and biochemical reaction-diffusion processes. The primary aim of this study is the development of a robust computational framework for biomedical flow simulation. Central to this development are adaptive mesh techniques based on a posteriori error estimates which rely on estimates of the errors of physical quantities of interest (goal quantities) to determine regions where the residual needs to be reduced. These estimates are calculated from the computed solution and the solution of an auxiliary (dual) problem containing information about the stability of the flow equations being solved and the goal quantities. Highlighting our current application of interest, this exposition will present ideas in the context of a study of blood flow in an artery containing an aneurysm. We begin with an overview of the biology, which motivates the need for accurate computation of quantities which are implicated in the growth of aneurysms (such as the shear forces exerted by the blood flow on their vessel walls). Then, we proceed to present the finite-element scheme for solving the governing Navier-Stokes equations for blood flow and provide insight into the numerical analysis to determine the required error indicators for the goal quantities. Finally, we present an implementation of this adaptive scheme in FEniCS and apply it to our problem, resulting in the construction of optimal computational meshes to compute our goal quantities related to aneurysm growth to a specified accuracy requirement with minimal work.},
  year = {2009},
  journal = {Talk at the Fifth M.I.T. Conference on Computational Fluid and Solid Mechanics, Cambridge, Massachusetts},
}