Flagellum Pumping Efficiency in Shear-Thinning Viscoelastic Fluids.

Microorganism motility often takes place within complex, viscoelastic fluid environments, e.g., sperm in cervicovaginal mucus and bacteria in biofilms.

Simulating cardiac fluid dynamics in the human heart.

Cardiac fluid dynamics fundamentally involves interactions between complex blood flows and the structural deformations of the muscular heart walls and the thin valve leaflets. There has been longstanding scientific, engineering, and medical interest in creating mathematical models of the heart that capture, explain, and predict these fluid-structure interactions (FSIs). However, existing computational models that account for interactions among the blood, the actively contracting myocardium, and the valves are limited in their abilities to predict valve performance, capture fine-scale flow features, or use realistic descriptions of tissue biomechanics.

Benchmarking the Immersed Boundary Method for Viscoelastic Flows.

We present and analyze a series of benchmark tests regarding the application of the immersed boundary (IB) method to viscoelastic flows through and around non-trivial, stationary geometries. The IB method is widely used to simulate biological fluid dynamics and other modeling scenarios in which a structure is immersed in a fluid. Although the IB method has been most commonly used to model systems involving viscous incompressible fluids, it also can be applied to visoelastic fluids, and has enabled the study of a wide variety of dynamical problems including the settling of vesicles and the swimming of elastic filaments in fluids modeled by the Oldroyd-B constitutive equation.

A stabilized linear finite element method for anisotropic poroelastodynamics with application to cardiac perfusion.

We propose a variational multiscale method stabilization of a linear finite element method for nonlinear poroelasticity. Our approach is suitable for the implicit time integration of poroelastic formulations in which the solid skeleton is anisotropic and incompressible. A detailed numerical methodology is presented for a monolithic formulation that includes both structural dynamics and Darcy flow.

Simulating Cardiac Fluid Dynamics in the Human Heart.

Cardiac fluid dynamics fundamentally involves interactions between complex blood flows and the structural deformations of the muscular heart walls and the thin, flexible valve leaflets. There has been longstanding scientific, engineering, and medical interest in creating mathematical models of the heart that capture, explain, and predict these fluid-structure interactions. However, existing computational models that account for interactions among the blood, the actively contracting myocardium, and the cardiac valves are limited in their abilities to predict valve performance, resolve fine-scale flow features, or use realistic descriptions of tissue biomechanics.

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction.

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface.

A model of fluid-structure and biochemical interactions for applications to subclinical leaflet thrombosis.

Subclinical leaflet thrombosis (SLT) is a potentially serious complication of aortic valve replacement with a bioprosthetic valve in which blood clots form on the replacement valve. SLT is associated with increased risk of transient ischemic attacks and strokes and can progress to clinical leaflet thrombosis. SLT following aortic valve replacement also may be related to subsequent structural valve deterioration, which can impair the durability of the valve replacement.

A Nodal Immersed Finite Element-Finite Difference Method.

The immersed finite element-finite difference (IFED) method is a computational approach to modeling interactions between a fluid and an immersed structure. The IFED method uses a finite element (FE) method to approximate the stresses, forces, and structural deformations on a and a finite difference (FD) method to approximate the momentum and enforce incompressibility of the entire fluid-structure system on a The fundamental approach used by this method follows the immersed boundary framework for modeling fluid-structure interaction (FSI), in which a force spreading operator prolongs structural forces to a Cartesian grid, and a velocity interpolation operator restricts a velocity field defined on that grid back onto the structural mesh. With an FE structural mechanics framework, force spreading first requires that the force itself be projected onto the finite element space.

Microphysiological model of PIK3CA-driven vascular malformations reveals a role of dysregulated Rac1 and mTORC1/2 in lesion formation.

Somatic activating mutations of are associated with development of vascular malformations (VMs). Here, we describe a microfluidic model of -driven VMs consisting of human umbilical vein endothelial cells expressing activating mutations embedded in three-dimensional hydrogels. We observed enlarged, irregular vessel phenotypes and the formation of cyst-like structures consistent with clinical signatures and not previously observed in cell culture models.

Rule-based definition of muscle bundles in patient-specific models of the left atrium.

Atrial fibrillation (AF) is the most common arrhythmia encountered clinically, and as the population ages, its prevalence is increasing. Although the CHADS VASc score is the most used risk-stratification system for stroke risk in AF, it lacks personalization. Patient-specific computer models of the atria can facilitate personalized risk assessment and treatment planning.

Patient-Specific Immersed Finite Element-Difference Model of Transcatheter Aortic Valve Replacement.

Transcatheter aortic valve replacement (TAVR) first received FDA approval for high-risk surgical patients in 2011 and has been approved for low-risk surgical patients since 2019. It is now the most common type of aortic valve replacement, and its use continues to accelerate. Computer modeling and simulation (CM&S) is a tool to aid in TAVR device design, regulatory approval, and indication in patient-specific care.

Semi-Automated Construction of Patient-Specific Aortic Valves from Computed Tomography Images.

This paper presents a semi-automatic method for the construction of volumetric models of the aortic valve using computed tomography angiography images. Although the aortic valve typically cannot be segmented directly from a computed tomography angiography image, the method described herein uses manually selected samples of an aortic segmentation derived from this image to inform the construction. These samples capture certain physiologic landmarks and are used to construct a volumetric valve model.

On the Lagrangian-Eulerian Coupling in the Immersed Finite Element/Difference Method.

The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method's performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions.

Bioprosthetic aortic valve diameter and thickness are directly related to leaflet fluttering: Results from a combined experimental and computational modeling study.

Objective: Bioprosthetic heart valves (BHVs) are commonly used in surgical and percutaneous valve replacement. The durability of percutaneous valve replacement is unknown, but surgical valves have been shown to require reintervention after 10 to 15 years. Further, smaller-diameter surgical BHVs generally experience higher rates of prosthesis-patient mismatch, which leads to higher rates of failure.

A Hybrid Semi-Lagrangian Cut Cell Method for Advection-Diffusion Problems with Robin Boundary Conditions in Moving Domains.

We present a new discretization approach to advection-diffusion problems with Robin boundary conditions on complex, time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. [8] to discretize the Laplace operator and Robin boundary condition.

Adherens junction engagement regulates functional patterning of the cardiac pacemaker cell lineage.

Cardiac pacemaker cells (CPCs) rhythmically initiate the electrical impulses that drive heart contraction. CPCs display the highest rate of spontaneous depolarization in the heart despite being subjected to inhibitory electrochemical conditions that should theoretically suppress their activity. While several models have been proposed to explain this apparent paradox, the actual molecular mechanisms that allow CPCs to overcome electrogenic barriers to their function remain poorly understood.

Models for plasma kinetics during simultaneous therapeutic plasma exchange and extracorporeal membrane oxygenation.

This paper focuses on the derivation and simulation of mathematical models describing new plasma fraction in blood for patients undergoing simultaneous extracorporeal membrane oxygenation and therapeutic plasma exchange. Models for plasma exchange with either veno-arterial or veno-venous extracorporeal membrane oxygenation are considered. Two classes of models are derived for each case, one in the form of an algebraic delay equation and another in the form of a system of delay differential equations.

A poroelastic immersed finite element framework for modelling cardiac perfusion and fluid-structure interaction.

Modern approaches to modelling cardiac perfusion now commonly describe the myocardium using the framework of poroelasticity. Cardiac tissue can be described as a saturated porous medium composed of the pore fluid (blood) and the skeleton (myocytes and collagen scaffold). In previous studies fluid-structure interaction in the heart has been treated in a variety of ways, but in most cases, the myocardium is assumed to be a hyperelastic fibre-reinforced material.

Mathematical modeling of the impact of recirculation on exchange kinetics in tandem extracorporeal membrane oxygenation and therapeutic plasma exchange.

Vascular access connection configurations during tandem extracorporeal membrane oxygenation (ECMO) and therapeutic plasma exchange (TPE) may impact exchange kinetics. In these tandem procedures, typically the TPE inlet line is proximal to the TPE return line with respect to blood flow in the ECMO device, maximizing the opportunity for replacement fluid homogenization within the ECMO circuit. However, if TPE inlet and return line connections are switched, recirculation-a phenomenon in which replacement fluid leaving the TPE return line is prematurely drawn into the TPE inlet line prior to satisfactory homogenization within the ECMO circuit-will occur.

Stabilization approaches for the hyperelastic immersed boundary method for problems of large-deformation incompressible elasticity.

The immersed boundary method is a mathematical framework for modeling fluid-structure interaction. This formulation describes the momentum, viscosity, and incompressibility of the fluid-structure system in Eulerian form, and it uses Lagrangian coordinates to describe the structural deformations, stresses, and resultant forces. Integral transforms with Dirac delta function kernels connect the Eulerian and Lagrangian frames.

Fluid-Structure Interaction Models of Bioprosthetic Heart Valve Dynamics in an Experimental Pulse Duplicator.

Computer modeling and simulation is a powerful tool for assessing the performance of medical devices such as bioprosthetic heart valves (BHVs) that promises to accelerate device design and regulation. This study describes work to develop dynamic computer models of BHVs in the aortic test section of an experimental pulse-duplicator platform that is used in academia, industry, and regulatory agencies to assess BHV performance. These computational models are based on a hyperelastic finite element extension of the immersed boundary method for fluid-structure interaction (FSI).

Analysis of multi-electron, multi-step homogeneous catalysis by rotating disc electrode voltammetry: theory, application, and obstacles.

Rotating disc electrode (RDE) voltammetry has been widely adopted for the study of heterogenized molecular electrocatalysts for multi-step fuel-forming reactions but this tool has never been comprehensively applied to their homogeneous analogues. Here, the utility and limitations of RDE techniques for mechanistic and kinetic analysis of homogeneous molecular catalysts that mediate multi-electron, multi-substrate redox transformations are explored. Using the ECEC' reaction mechanism as a case study, two theoretical models are derived based on the Nernst diffusion layer model and the Hale transformation.

An Immersed Interface Method for Discrete Surfaces.

Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-solid interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces.