No abstract

Abstract. We investigate several robust preconditioners for solving the saddle-point linear systems that arise from spatial discretization of unsteady and steady variable-coefficient Stokes equations on a uniform staggered grid. Building on the success of using the classical projection method as a preconditioner for the coupled velocity-pressure system [B. E. Griffith, J. Comp. Phys., 228 (2009), pp. 7565-7595 ], as well as established techniques for steady and unsteady Stokes flow in the finite-element literature, we construct preconditioners that employ independent generalized Helmholtz and Poisson solvers for the velocity and pressure subproblems. We demonstrate that only a single cycle of a standard geometric multigrid algorithm serves as an effective inexact solver for each of these subproblems. Contrary to traditional wisdom, we find that the Stokes problem can be solved nearly as efficiently as the independent pressure and velocity subproblems, making the overall cost of solving the Stokes system comparable to the cost of classical projection or fractional step methods for incompressible flow, even for steady flow and in the presence of large density and viscosity contrasts. Two of the five preconditioners considered here are found to be robust to GMRES restarts and to increasing problem size, making them suitable for large-scale problems. Our work opens many possibilities for constructing novel unsplit temporal integrators for finite-volume spatial discretizations of the equations of low Mach and incompressible flow dynamics.

Computational models have the potential to provide precise estimates of stresses and strains associated with sites of coronary plaque rupture. However, lack of adequate mathematical description of diseased human vessel wall mechanical properties is hindering computational accuracy. The goal of this study is to characterize the behavior of diseased human coronary and carotid arteries using planar biaxial testing. Diseased coronary specimens exhibit relatively high stiffness (50–210 kPa) and low extensibility (1–10%) at maximum equibiaxial stress (250 kPa) compared to human carotid specimens and values commonly reported for porcine coronary arteries. A thick neointimal layer observed histologically appears to be associated with heightened stiffness and the direction of anisotropy of the specimens. Fung, Choi-Vito and modified Mooney-Rivlin constitutive equations fit the multiaxial data from multiple stress protocols well, and parameters from representative coronary specimens were utilized in a finite element model with fluid-solid interactions. Computed locations of maximal stress and strain are substantially altered, and magnitudes of maximum principal stress (48–65 kPa) and strain (6.5–8%) in the vessel wall are lower than previously predicted using parameters from uniaxial tests. Taken together, the results demonstrate the importance of utilizing disease-matched multiaxial constitutive relationships within patient-specific computational models to accurately predict stress and strain within diseased coronary arteries.

No abstract