Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Considerations for Advection-Diffusion Problems in Cardiovascular Hemodynamics
174
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Numerical simulations of cardiovascular mass transport pose significant challenges due to the wide range of Peclet numbers and backflow at Neumann boundaries. In this paper we present and discuss several numerical tools to address these challenges in the context of a stabilized finite element computational framework. To overcome numerical instabilities when backflow occurs at Neumann boundaries, we propose an approach based on the prescription of the total flux. In addition, we introduce a consistent flux outflow boundary condition and demonstrate its superior performance over the traditional zero diffusive flux boundary condition. Lastly, we discuss discontinuity capturing (DC) stabilization techniques to address the well-known oscillatory behavior of the solution near the concentration front in advection-dominated flows.We present numerical examples in both idealized and patient-specific geometries to demonstrate the efficacy of the proposed procedures. The three contributions dis-cussed in this paper enable to successfully address commonly found challenges when simulating mass transport processes in cardiovascular flows.
rate research
Read More
Computing the solution of linear systems of equations is invariably the most time consuming task in the numerical solutions of PDEs in many fields of computational science. In this study, we focus on the numerical simulation of cardiovascular hemodynamics with rigid and deformable walls, discretized in space and time through the variational multi-scale finite element method. We focus on three approaches: the problem agnostic generalized minimum residual (GMRES) and stabilized bi-conjugate gradient (BICGS) methods, and a recently proposed, problem specific, bi-partitioned (BIPN) method. We also perform a comparative analysis of several preconditioners, including diagonal, block-diagonal, incomplete factorization, multi-grid, and resistance based methods. Solver performance and matrix characteristics (diagonal dominance, symmetry, sparsity, bandwidth and spectral properties) are first examined for an idealized cylindrical geometry with physiologic boundary conditions and then successively tested on several patient-specific anatomies representative of realistic cardiovascular simulation problems. Incomplete factorization pre-conditioners provide the best performance and results in terms of both strong and weak scalability. The BIPN method was found to outperform other methods in patient-specific models with rigid walls. In models with deformable walls, BIPN was outperformed by BICG with diagonal and Incomplete LU preconditioners.
A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid using a finite difference scheme. A volume-of-fluid formulation (Hirt and Nichols (1981, J. Comput. Phys., 39, 201)), which has been widely used for multiphase flow simulations, is applied to describing the multi-component geometry. The temporal change in the solid deformation is described in the Eulerian frame by updating a left Cauchy-Green deformation tensor, which is used to express constitutive equations for nonlinear Mooney-Rivlin materials. In this paper, various verifications and validations of the present full Eulerian method, which solves the fluid and solid motions on a fixed grid, are demonstrated, and the numerical accuracy involved in the fluid-structure coupling problems is examined.
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 Dakotas MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.
This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection-diffusion equations. However, to the best of the authors knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. The algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.
In this paper, a multiple-distribution-function lattice Boltzmann method (MDF-LBM) with multiple-relaxation-time model is proposed for incompressible Navier-Stokes equations (NSEs) which are considered as the coupled convection-diffusion equations (CDEs). Through direct Taylor expansion analysis, we show that the Navier-Stokes equations can be recovered correctly from the present MDF-LBM, and additionally, it is also found that the velocity and pressure can be directly computed through the zero and first-order moments of distribution function. Then in the framework of present MDF-LBM, we develop a locally computational scheme for the velocity gradient where the first-order moment of the non-equilibrium distribution is used, this scheme is also extended to calculate the velocity divergence, strain rate tensor, shear stress and vorticity. Finally, we also conduct some simulations to test the MDF-LBM, and find that the numerical results not only agree with some available analytical and numerical solutions, but also have a second-order convergence rate in space.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
