Nested and parallel sparse algorithms for arterial fluid mechanics computations with boundary layer mesh refinement

Manguoğlu, Murat
Sameh, Ahmed H.
Tezduyar, Tayfun E.
Arterial fluid-structure interaction (FSI) computations involve a number of numerical challenges. Because blood flow is incompressible, iterative solution of the fluid mechanics part of the linear equation system at every nonlinear iteration of each time step is one of those challenges, especially for computations over slender domains and in the presence of boundary layer mesh refinement. In this paper we address that challenge. As test cases, we use equation systems from stabilized finite element computation of a bifurcating middle cerebral artery segment with aneurysm, with thin layers of elements near the arterial wall. We show how the preconditioning techniques, we propose for solving these large sparse nonsymmetric systems, perform at different time steps of the computation over a cardiac cycle. We also present a new hybrid parallel sparse linear system solver 'DD-Spike' and demonstrate its scalability. Copyright (C) 2010 John Wiley & Sons, Ltd.


Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme
Bozkaya, Canan; Tezer, Münevver (Wiley, 2006-06-20)
A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydro-dynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time-dependent convection-diffusion-type equations. These equations are solved by using DRBEM which treats the time and the space deriva...
The DRBEM solution of incompressible MHD flow equations
Bozkaya, Nuray; Tezer, Münevver (Wiley, 2011-12-10)
This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier-Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function-vorticity-magnetic induction-current density formulation of the full MHD equations in 2D. The stream f...
Time-domain BEM solution of convection-diffusion-type MHD equations
Bozkaya, N.; Tezer, Münevver (Wiley, 2008-04-20)
The two-dimensional convection-diffusion-type equations are solved by using the boundary element method (BEM) based on the time-dependent fundamental solution. The emphasis is given on the solution of magnetohydrodynamic (MHD) duct flow problems with arbitrary wall conductivity. The boundary and time integrals in the BEM formulation are computed numerically assuming constant variations of the unknowns on both the boundary elements and the time intervals. Then, the solution is advanced to the steady-state it...
Improvements to compressible Euler methods for low-Mach number flows
Sabanca, M; Brenner, G; Alemdaroglu, N (Wiley, 2000-09-30)
In the present study improvements to numerical algorithms for the solution of the compressible Euler equations at low Mach numbers are investigated. To solve flow problems for a wide range of Mach numbers, from the incompressible limit to supersonic speeds, preconditioning techniques are frequently employed. On the other hand, one can achieve the same aim by using a suitably modified acoustic damping method. The solution algorithm presently under consideration is based on Roe's approximate Riemann solver [R...
Two-level finite element method with a stabilizing subgrid for the incompressible Navier-Stokes equations
NESLİTÜRK, ALİ İHSAN; Aydın Bayram, Selma; Tezer, Münevver (Wiley, 2008-10-20)
We consider the Galerkin finite element method for the incompressible Navier-Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite-dimensional spaces emploved consist of piecewise continuous linear interpolants enriched with the residual-free bubble functions. To find the bubble part of the Solution, a two-level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier-Stokes equation is ...
Citation Formats
M. Manguoğlu, A. H. Sameh, and T. E. Tezduyar, “Nested and parallel sparse algorithms for arterial fluid mechanics computations with boundary layer mesh refinement,” INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, pp. 135–149, 2011, Accessed: 00, 2020. [Online]. Available: