BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Alison L. Marsden

Alison L. Marsden contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phphysics.flu-dynmath.NANumerical Analysisphysics.med-phTissues and OrgansApplicationsComputational Engineering, Finance, and ScienceQuantitative Methods

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2026arXiv

Simulations Predict Improved Valve Performance Without Direct Leaflet Intervention After Neonatal Truncus Arteriosus Repair

Truncus arteriosus (TA) is a rare and severe congenital heart disease. Quadricuspid valve morphology occurs in 25% of all TA patients and is linked to regurgitation and increased risk of re-operation. It remains unclear how hemodynamic changes after TA repair alter valve performance. This study simulated pre- and postoperative conditions in a neonatal TA patient to investigate valve performance without direct intervention. We hypothesize that valve performance before and after truncal repair can be predicted in-silico, matching in-vivo imaging and identifying mechanisms how hemodynamic changes after repair will reduce valve regurgitation without direct intervention. Pre- and postoperative CT images of a neonatal patient with quadricuspid valve were segmented. Free edge length and geometric height from the patient's echocardiogram were used to model the valve. For the preoperative condition, ventricular pressures were set equal modeling an unrestricted ventricular septal defect. Systemic and pulmonary resistances were tuned based on the patient's Qp:Qs ratio. For the postoperative condition, boundary conditions were modified to mimic patient-specific hemodynamics after TA repair. The preoperative simulation confirmed mild valve regurgitation seen in-vivo. Interaction between asymmetric flow and surrounding vessel resulted in asymmetric opening and closing. Poor central coaptation led to a central regurgitant jet toward the septum. Altered postoperative hemodynamics improved coaptation and eliminated regurgitation, as seen in-vivo. This modeling approach reproduced in-vivo pre- and postoperative valve performance and identified mechanisms improving coaptation after TA repair. TA repair led to elimination of regurgitation due to enhanced central coaptation. Thus, altered postoperative hemodynamic conditions after TA repair may improve valve performance without direct leaflet intervention.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Automated generation of 0D and 1D reduced-order models of patient-specific blood flow

Three-dimensional (3D) cardiovascular fluid dynamics simulations typically require hours to days of computing time on a high-performance computing cluster. One-dimensional (1D) and lumped-parameter zero-dimensional (0D) models show great promise for accurately predicting blood bulk flow and pressure waveforms with only a fraction of the cost. They can also accelerate uncertainty quantification, optimization, and design parameterization studies. Despite several prior studies generating 1D and 0D models and comparing them to 3D solutions, these were typically limited to either 1D or 0D and a singular category of vascular anatomies. This work proposes a fully automated and openly available framework to generate and simulate 1D and 0D models from 3D patient-specific geometries, automatically detecting vessel junctions and stenosis segments. Our only input is the 3D geometry; we do not use any prior knowledge from 3D simulations. All computational tools presented in this work are implemented in the open-source software platform SimVascular. We demonstrate the reduced-order approximation quality against rigid-wall 3D solutions in a comprehensive comparison with N=72 publicly available models from various anatomies, vessel types, and disease conditions. Relative average approximation errors of flows and pressures typically ranged from 1% to 10% for both 1D and 0D models, measured at the outlets of terminal vessel branches. In general, 0D model errors were only slightly higher than 1D model errors despite requiring only a third of the 1D runtime. Automatically generated ROMs can significantly speed up model development and shift the computational load from high-performance machines to personal computers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Controlled Comparison of Simulated Hemodynamics across Tricuspid and Bicuspid Aortic Valves

Bicuspid aortic valve is the most common congenital heart defect, affecting 1-2% of the global population. Patients with bicuspid valves frequently develop dilation and aneurysms of the ascending aorta. Both hemodynamic and genetic factors are believed to contribute to dilation, yet the precise mechanism underlying this progression remains under debate. Controlled comparisons of hemodynamics in patients with different forms of bicuspid valve disease are challenging because of confounding factors, and simulations offer the opportunity for direct and systematic comparisons. Using fluid-structure interaction simulations, we simulate flows through multiple aortic valve models in a patient-specific geometry. The aortic geometry is based on a healthy patient with no known aortic or valvular disease, which allows us to isolate the hemodynamic consequences of changes to the valve alone. Four fully-passive, elastic model valves are studied: a tricuspid valve and bicuspid valves with fusion of the left- and right-, right- and non-, and non- and left-coronary cusps. The resulting tricuspid flow is relatively uniform, with little secondary or reverse flow, and little to no pressure gradient across the valve. The bicuspid cases show localized jets of forward flow, excess streamwise momentum, elevated secondary and reverse flow, and clinically significant levels of stenosis. Localized high flow rates correspond to locations of dilation observed in patients, with the location related to which valve cusps are fused. Thus, the simulations support the hypothesis that chronic exposure to high local flow contributes to localized dilation and aneurysm formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical investigation of abdominal aortic aneurysm hemodynamics using the reduced unified continuum formulation for vascular fluid-structure interaction

We recently demonstrated the reduction of the unified continuum and variational multiscale formulation to a computationally efficient fluid-structure interaction (FSI) formulation via three sound modeling assumptions pertaining to the vascular wall. Similar to the coupled momentum method introduced by Figueroa et al., the resulting semi-discrete formulation yields a monolithically coupled FSI system posed in an Eulerian frame of reference with only a minor modification of the fluid boundary integral. To achieve uniform second-order temporal accuracy and user-controlled high-frequency algorithmic damping, we adopt the generalized-$α$ method for uniform temporal discretization of the entire coupled system. In conjunction with a fully consistent, segregated predictor multi-corrector algorithm preserving the block structure of the incompressible Navier-Stokes equations in the implicit solver's associated linear system, a three-level nested block preconditioner is adopted for improved representation of the Schur complement. In this work, we apply our reduced unified continuum formulation to an appropriately prestressed patient-specific abdominal aortic aneurysm and investigate the effects of varying spatial distributions of wall properties on hemodynamic and vascular wall quantities of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Validation of the Reduced Unified Continuum Formulation Against In Vitro 4D-Flow MRI

In our recent work, we introduced the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) and demonstrated enhanced solver accuracy, scalability, and performance compared to conventional approaches. We further verified the formulation against Womersley's deformable wall theory. In this study, we assessed its performance in a compliant patient-specific aortic model by leveraging 3D printing, 2D magnetic resonance imaging (MRI), and 4D-flow MRI to extract high-resolution anatomical and hemodynamic information from an in vitro flow circuit. To accurately reflect experimental conditions, we additionally enabled in-plane vascular motion at each inlet and outlet, and implemented viscoelastic external tissue support and vascular tissue prestressing. Validation of our formulation is achieved through close quantitative agreement in pressures, lumen area changes, pulse wave velocity, and early systolic velocities, as well as qualitative agreement in late systolic flow structures. Our validated suite of FSI techniques can be used to investigate vascular disease initiation, progression, and treatment at a computational cost on the same order as that of rigid-walled simulations. This study is the first to validate a cardiovascular FSI formulation against an in vitro flow circuit involving a compliant vascular phantom of complex patient-specific anatomy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A reduced unified continuum formulation for vascular fluid-structure interaction

We recently derived the unified continuum and variational multiscale formulation for fluid-structure interaction (FSI) using the Gibbs free energy. Restricting our attention to vascular FSI, we now reduce this arbitrary Lagrangian-Eulerian (ALE) formulation by adopting three assumptions for the vascular wall. The resulting reduced unified continuum formulation achieves monolithic FSI coupling in the Eulerian frame through a simple modification of the fluid boundary integral. While ostensibly similar to the semi-discrete formulation of the coupled momentum method, its underlying derivation does not rely on an assumption of a fictitious body force in the elastodynamics sub-problem and therefore represents a direct simplification of the ALE method. Uniform temporal discretization is performed via the generalized-$α$ scheme. In contrast to the predominant approach yielding only first-order accuracy for pressure, we collocate both pressure and velocity at the intermediate time step to achieve uniform second-order temporal accuracy. In conjunction with quadratic tetrahedral elements, our methodology offers higher-order temporal and spatial accuracy for quantities of clinical interest. Furthermore, without loss of consistency, a segregated predictor multi-corrector algorithm is developed to preserve the same block structure as for the incompressible Navier-Stokes equations in the implicit solver's associated linear system. Block preconditioning of a monolithically coupled FSI system is therefore made possible for the first time. Compared to alternative preconditioners, our three-level nested block preconditioner, which improves representation of the Schur complement, demonstrates robust performance over a wide range of physical parameters. We present verification against Womersley's deformable wall theory and additionally develop practical modeling techniques for clinical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unified continuum and variational multiscale formulation for fluids, solids, and fluid-structure interaction

We develop a unified continuum modeling framework for viscous fluids and hyperelastic solids using the Gibbs free energy as the thermodynamic potential. This framework naturally leads to a pressure primitive variable formulation for the continuum body, which is well-behaved in both compressible and incompressible regimes. Our derivation also provides a rational justification of the isochoric-volumetric additive split of free energies in nonlinear continuum mechanics. The variational multiscale analysis is performed for the continuum model to construct a foundation for numerical discretization. We first consider the continuum body instantiated as a hyperelastic material and develop a variational multiscale formulation for the hyper-elastodynamic problem. The generalized-alpha method is applied for temporal discretization. A segregated algorithm for the nonlinear solver is designed and carefully analyzed. Second, we apply the new formulation to construct a novel unified formulation for fluid-solid coupled problems. The variational multiscale formulation is utilized for spatial discretization in both fluid and solid subdomains. The generalized-alpha method is applied for the whole continuum body, and optimal high-frequency dissipation is achieved in both fluid and solid subproblems. A new predictor multi-corrector algorithm is developed based on the segregated algorithm to attain a good balance between robustness and efficiency. The efficacy of the new formulations is examined in several benchmark problems. The results indicate that the proposed modeling and numerical methodologies constitute a promising technology for biomedical and engineering applications, particularly those necessitating incompressible models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluid-structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation

In this work, we present a computational fluid-structure interaction (FSI) study for a healthy patient-specific pulmonary arterial tree using the unified continuum and variational multiscale (VMS) formulation we previously developed. The unified framework is particularly well-suited for FSI, as the fluid and solid sub-problems are addressed in essentially the same manner and can thus be uniformly integrated in time with the generalized-$α$ method. In addition, the VMS formulation provides a mechanism for large-eddy simulation in the fluid sub-problem and pressure stabilization in the solid sub-problem. The FSI problem is solved in a quasi-direct approach, in which the pressure and velocity in the unified continuum body are first solved, and the solid displacement is then obtained via a segregated algorithm and prescribed as a boundary condition for the mesh motion. Results of the pulmonary arterial FSI simulation are presented and compared against those of a rigid wall simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-fidelity estimators for coronary circulation models under clinically-informed data uncertainty

Numerical models are increasingly used for non-invasive diagnosis and treatment planning in coronary artery disease, where service-based technologies have proven successful in identifying hemodynamically significant and hence potentially dangerous vascular anomalies. Despite recent progress towards clinical adoption, many results in the field are still based on a deterministic characterization of blood flow, with no quantitative assessment of the variability of simulation outputs due to uncertainty from multiple sources. In this study, we focus on parameters that are essential to construct accurate patient-specific representations of the coronary circulation, such as aortic pressure waveform, intramyocardial pressure and quantify how their uncertainty affects clinically relevant model outputs. We construct a deformable model of the left coronary artery subject to a prescribed inlet pressure and with open-loop outlet boundary conditions, treating fluid-structure interaction through an Arbitrary-Lagrangian-Eulerian frame of reference. Random input uncertainty is estimated directly from repeated clinical measurements from intra-coronary catheterization and complemented by literature data. We also achieve significant computational cost reductions in uncertainty propagation thanks to multifidelity Monte Carlo estimators of the outputs of interest, leveraging the ability to generate, at practically no cost, one- and zero-dimensional low-fidelity representations of left coronary artery flow, with appropriate boundary conditions. The results demonstrate how the use of multi-fidelity control variate estimators leads to significant reductions in variance and accuracy improvements with respect to traditional Monte-Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics

Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. We propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these low-fidelity and high-fidelity models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. We seek to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to Monte Carlo and multilevel Monte Carlo estimators based only on three-dimensional simulations. We demonstrate significant reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities and healthy models show larger reductions than local quantities and diseased model, as the latter rely more heavily on the highest fidelity model evaluations. In all cases, our workflow coupling Dakota's MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effects of clinically-derived parametric data uncertainty in patient-specific coronary simulations with deformable walls

Cardiovascular simulations are increasingly used for non-invasive diagnosis of cardiovascular disease, to guide treatment decisions, and in the design of medical devices. Quantitative assessment of the variability of simulation outputs due to input uncertainty is a key step toward further integration of cardiovascular simulations in the clinical workflow. In this study, we present uncertainty quantification in computational models of the coronary circulation to investigate the effect of uncertain parameters, including coronary pressure waveform, intramyocardial pressure, morphometry exponent, and the vascular wall Young's modulus. We employ a left coronary artery model with deformable vessel walls, simulated via an ALE framework for FSI, with a prescribed inlet pressure and open-loop lumped parameter network outlet boundary conditions. Stochastic modeling of the uncertain inputs is determined from intra-coronary catheterization data or gathered from the literature. Uncertainty propagation is performed using several approaches including Monte Carlo, Quasi MC, stochastic collocation, and multiwavelet stochastic expansion. Variabilities in QoI, including branch pressure, flow, wall shear stress, and wall deformation are assessed. We find that uncertainty in inlet pressures and intramyocardial pressures significantly affect all resulting QoIs, while uncertainty in elastic modulus only affects the mechanical response of the vascular wall. Variability in the morphometry exponent has little effect on coronary hemodynamics or wall mechanics. Finally, we compare convergence behaviors of statistics of QoIs using several uncertainty propagation methods. From the simulation results, we conclude that the multi-wavelet stochastic expansion shows superior accuracy and performance against Quasi Monte Carlo and stochastic collocation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The nested block preconditioning technique for the incompressible Navier-Stokes equations with emphasis on hemodynamic simulations

We develop a novel iterative solution method for the incompressible Navier-Stokes equations with boundary conditions coupled with reduced models. The iterative algorithm is designed based on the variational multiscale formulation and the generalized-$α$ scheme. The spatiotemporal discretization leads to a block structure of the resulting consistent tangent matrix in the Newton-Raphson procedure. As a generalization of the conventional block preconditioners, a three-level nested block preconditioner is introduced to attain a better representation of the Schur complement, which plays a key role in the overall algorithm robustness and efficiency. This approach provides a flexible, algorithmic way to handle the Schur complement for problems involving multiscale and multiphysics coupling. The solution method is implemented and benchmarked against experimental data from the nozzle challenge problem issued by the US Food and Drug Administration. The robustness, efficiency, and parallel scalability of the proposed technique are then examined in several settings, including moderately high Reynolds number flows and physiological flows with strong resistance effect due to coupled downstream vasculature models. Two patient-specific hemodynamic simulations, covering systemic and pulmonary flows, are performed to further corroborate the efficacy of the proposed methodology.

0 citations0 reviewsNo code linkedOpen access