A parallel sparse algorithm targeting arterial fluid mechanics computations

2011-09-01
Manguoğlu, Murat
Sameh, Ahmed H.
Tezduyar, Tayfun E.
Iterative solution of large sparse nonsymmetric linear equation systems is one of the numerical challenges in arterial fluid-structure interaction computations. This is because the fluid mechanics parts of the fluid + structure block of the equation system that needs to be solved at every nonlinear iteration of each time step corresponds to incompressible flow, the computational domains include slender parts, and accurate wall shear stress calculations require boundary layer mesh refinement near the arterial walls. We propose a hybrid parallel sparse algorithm, domain-decomposing parallel solver (DDPS), to address this challenge. As the test case, we use a fluid mechanics equation system generated by starting with an arterial shape and flow field coming from an FSI computation and performing two time steps of fluid mechanics computation with a prescribed arterial shape change, also coming from the FSI computation. We show how the DDPS algorithm performs in solving the equation system and demonstrate the scalability of the algorithm.
COMPUTATIONAL MECHANICS

Suggestions

An Asymptotic-Numerical Hybrid Method for Solving Singularly Perturbed Linear Delay Differential Equations
Cengizci, Süleyman (Hindawi Limited, 2017)
In thiswork, approximations to the solutions of singularly perturbed second-order linear delay differential equations are studied. We firstly use two-term Taylor series expansion for the delayed convection term and obtain a singularly perturbed ordinary differential equation (ODE). Later, an efficient and simple asymptotic method so called Successive Complementary Expansion Method (SCEM) is employed to obtain a uniformly valid approximation to this corresponding singularly perturbed ODE. As the final step, ...
A NOVEL PARTITIONING METHOD FOR ACCELERATING THE BLOCK CIMMINO ALGORITHM
Torun, F. Sukru; Manguoğlu, Murat; Aykanat, Cevdet (2018-01-01)
We propose a novel block-row partitioning method in order to improve the convergence rate of the block Cimmino algorithm for solving general sparse linear systems of equations. The convergence rate of the block Cimmino algorithm depends on the orthogonality among the block rows obtained by the partitioning method. The proposed method takes numerical orthogonality among block rows into account by proposing a row inner-product graph model of the coefficient matrix. In the graph partitioning formulation define...
A Fully Implicit Finite Volume Lattice Boltzmann Method for Turbulent Flow
Cevik, Fatih; Albayrak, Kahraman (2017-08-01)
Almost all schemes existed in the literature to solve the Lattice Boltzmann Equation like stream & collide, finite difference, finite element, finite volume schemes are explicit. However, it is known fact that implicit methods utilizes better stability and faster convergence compared to the explicit methods. In this paper, a method named herein as Implicit Finite Volume Lattice BoltzmannMethod (IFVLBM) for incompressible laminar and turbulent flows is proposed and it is applied to some 2D benchmark test cas...
A modal superposition method for non-linear structures
Kuran, B; Özgüven, Hasan Nevzat (Elsevier BV, 1996-01-25)
The dynamic response of multi-degree of freedom (MDOF) non-linear structures is usually determined by the numerical integration of equations of motion. This is computationally very costly for steady state response analysis. In this study, a powerful and economical method is developed for the harmonic response analysis of non-linear structures. In this method, the equations of motion are first converted into a set of non-linear algebraic equations, and then the number of equations to be solved is reduced by ...
A Multithreaded Recursive and Nonrecursive Parallel Sparse Direct Solver
Bölükbaşı, Ercan Selçuk (2016-01-01)
Sparse linear system of equations often arises after discretization of the partial differential equations (PDEs) such as computational fluid dynamics, material science, and structural engineering. There are, however, sparse linear systems that are not governed by PDEs, some examples of such applications are circuit simulations, power network analysis, and social network analysis. For solution of sparse linear systems one can choose using either a direct or an iterative method. Direct solvers are based on so...
Citation Formats
M. Manguoğlu, A. H. Sameh, and T. E. Tezduyar, “A parallel sparse algorithm targeting arterial fluid mechanics computations,” COMPUTATIONAL MECHANICS, pp. 377–384, 2011, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/32846.