Source author record

Charles Puelz

Charles Puelz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Tissues and Organs math.DS physics.flu-dyn math.SP Numerical Analysis physics.comp-ph physics.med-ph Quantitative Methods

Catalog footprint

What is connected

7works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A computational study of aortic reconstruction in single ventricle patients

Patients with hypoplastic left heart syndrome (HLHS) are born with an underdeveloped left heart. They typically receive a sequence of surgeries that result in a single ventricle physiology called the Fontan circulation. While these patients usually survive into early adulthood, they are at risk for medical complications, partially due to their lower than normal cardiac output, which leads to insufficient cerebral and gut perfusion. While clinical imaging data can provide detailed insight into cardiovascular function within the imaged region, it is difficult to use this data for assessing deficiencies in the rest of the body and for deriving blood pressure dynamics. Data from patients used in this paper include three dimensional, magnetic resonance angiograms (MRA), time-resolved phase contrast cardiac magnetic resonance images (4D-MRI) and sphygmomanometer blood pressure measurements. The 4D-MRI images provide detailed insight into velocity and flow in vessels within the imaged region, but they cannot predict flow in the rest of the body, nor do they provide values of blood pressure. To remedy these limitations, this study combines the MRA, 4D-MRI, and pressure data with 1D fluid dynamics models to predict hemodynamics in the major systemic arteries, including the cerebral and gut vasculature. To study effects post-surgery we compare simulations for an HLHS patient with simulations for a matched control patient that has double outlet right ventricle (DORV) physiology with a native aorta. Our results show that the HLHS patient has hypertensive pressures in the brain as well as reduced flow to the gut. Wave-intensity analysis suggests that the HLHS patient has irregular circulatory function during light upright exercise conditions and that predicted wall-shear stresses are lower than normal, suggesting the HLHS patient may have hypertension.

preprint2022arXiv

Optimal fenestration of the Fontan circulation

In this paper, we develop a pulsatile compartmental model of the Fontan circulation and use it to explore the effects of a fenestration added to this physiology. A fenestration is a shunt between the systemic and pulmonary veins that is added either at the time of Fontan conversion or at a later time for the treatment of complications. This shunt increases cardiac output and decreases systemic venous pressure. However, these hemodynamic benefits are achieved at the expense of a decrease in the arterial oxygen saturation. The model developed this paper incorporates fenestration size as a parameter and describes both blood flow and oxygen transport. It is calibrated to clinical data from Fontan patients, and we use it to study the impact of a fenestration on several hemodynamic variables. In certain scenarios corresponding to high-risk Fontan physiology, we demonstrate the existence of an optimal fenestration size that maximizes oxygen delivery to the systemic tissues.

preprint2020arXiv

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. In special cases, our models reduce to single compartment ones for plasma exchange that have been validated with experimental data. We also show that the algebraic delay equations are forward Euler discretizations of the delay differential equations, with timesteps equal to transit times through model compartments. Numerical simulations are performed to compare different model types, to investigate the impact of plasma device port switching on the efficiency of the exchange process, and to study the sensitivity of the models to their parameters.

preprint2019arXiv

A sharp interface method for an immersed viscoelastic solid

The immersed boundary-finite element method (IBFE) is an approach to describing the dynamics of an elastic structure immersed in an incompressible viscous fluid. In this formulation, there are discontinuities in the pressure and viscous stress at fluid-structure interfaces. The standard immersed boundary approach, which connects the Lagrangian and Eulerian variables via integral transforms with regularized Dirac delta function kernels, smooths out these discontinuities, which generally leads to low order accuracy. This paper describes an approach to accurately resolve pressure discontinuities for these types of formulations, in which the solid may undergo large deformations. Our strategy is to decompose the physical pressure field into a sum of two pressure-like fields, one defined on the entire computational domain, which includes both the fluid and solid subregions, and one defined only on the solid subregion. Each of these fields is continuous on its domain of definition, which enables high accuracy via standard discretization methods without sacrificing sharp resolution of the pressure discontinuity. Numerical tests demonstrate that this method improves rates of convergence for displacements, velocities, stresses, and pressures, as compared to the conventional IBFE method. Further, it produces much smaller errors at reasonable numbers of degrees of freedom. The performance of this method is tested on several cases with analytic solutions, a nontrivial benchmark problem of incompressible solid mechanics, and an example involving a thick, actively contracting torus.

preprint2016arXiv

Comparison of reduced models for blood flow using Runge-Kutta discontinuous Galerkin methods

One-dimensional blood flow models take the general form of nonlinear hyperbolic systems but differ greatly in their formulation. One class of models considers the physically conserved quantities of mass and momentum, while another class describes mass and velocity. Further, the averaging process employed in the model derivation requires the specification of the axial velocity profile; this choice differentiates models within each class. Discrepancies among differing models have yet to be investigated. In this paper, we systematically compare several reduced models of blood flow for physiologically relevant vessel parameters, network topology, and boundary data. The models are discretized by a class of Runge-Kutta discontinuous Galerkin methods.

preprint2015arXiv

Numerical method of characteristics for one-dimensional blood flow

Mathematical modeling at the level of the full cardiovascular system requires the numerical approximation of solutions to a one-dimensional nonlinear hyperbolic system describing flow in a single vessel. This model is often simulated by computationally intensive methods like finite elements and discontinuous Galerkin, while some recent applications require more efficient approaches (e.g. for real-time clinical decision support, phenomena occurring over multiple cardiac cycles, iterative solutions to optimization/inverse problems, and uncertainty quantification). Further, the high speed of pressure waves in blood vessels greatly restricts the time step needed for stability in explicit schemes. We address both cost and stability by presenting an efficient and unconditionally stable method for approximating solutions to diagonal nonlinear hyperbolic systems. Theoretical analysis of the algorithm is given along with a comparison of our method to a discontinuous Galerkin implementation. Lastly, we demonstrate the utility of the proposed method by implementing it on small and large arterial networks of vessels whose elastic and geometrical parameters are physiologically relevant.

preprint2014arXiv

Spectral Approximation for Quasiperiodic Jacobi Operators

Quasiperiodic Jacobi operators arise as mathematical models of quasicrystals and in more general studies of structures exhibiting aperiodic order. The spectra of these self-adjoint operators can be quite exotic, such as Cantor sets, and their fine properties yield insight into associated dynamical systems. Quasiperiodic operators can be approximated by periodic ones, the spectra of which can be computed via two finite dimensional eigenvalue problems. Since long periods are necessary to get detailed approximations, both computational efficiency and numerical accuracy become a concern. We describe a simple method for numerically computing the spectrum of a period-$K$ Jacobi operator in $O(K^2)$ operations, and use it to investigate the spectra of Schrödinger operators with Fibonacci, period doubling, and Thue-Morse potentials.