Researcher profile

Ju Liu

Ju Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
7topics
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

7 published item(s)

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.

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.

preprint2022arXiv

Vascular fluid-structure interaction: unified continuum formulation, image-based mesh generation pipeline, and scalable fully implicit solver technology

We propose a computational framework for vascular fluid-structure interaction (FSI), focusing on biomechanical modeling, geometric modeling, and solver technology. The biomechanical model is constructed based on the unified continuum formulation. We highlight that the chosen time integration scheme differs from existing implicit FSI integration methods in that it is indeed second-order accurate, does not suffer from the overshoot phenomenon, and optimally dissipates high-frequency modes in both subproblems. We propose a pipeline for generating subject-specific meshes for FSI analysis for anatomically realistic geometric modeling. Unlike most existing methodologies that operate directly on the wall surface mesh, our pipeline starts from the image segmentation stage. With high-quality surface meshes obtained, the volumetric meshes are then generated, guaranteeing a boundary-layered mesh in the fluid subdomain and a matching mesh across the fluid-solid interface. In the last, we propose a combined suite of nonlinear and linear solver technologies. Invoking a segregated algorithm within the Newton-Raphson iteration, the problem reduces to solving two linear systems in the multi-corrector stage. The first linear system can be addressed by the algebraic multigrid (AMG) method. The matrix related to the balance equations presents a two-by-two block structure in both subproblems. Using the Schur complement reduction (SCR) technique reduces the problem to solving matrices of smaller sizes of the elliptic type, and the AMG method again becomes a natural candidate. The benefit of the unified formulation is demonstrated in parallelizing the solution algorithms as the number of unknowns matches in both subdomains. We use the Greenshields-Weller benchmark as well as a patient-specific vascular model to demonstrate the robustness, efficiency, and scalability of the overall FSI solver technology.

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.

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.

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.

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.