Search

Although coronary artery bypass graft (CABG) surgery is a well-established intervention, graft failure can occur, and the underlying mechanisms remain incompletely understood. The purpose of this prospective study is to utilize computational fluid dynamics (CFD) to investigate how graft hemodynamics one month post surgery may vary among graft types, which have different long-term patency rates. Twenty-four grafts from 10 participants (64.6 ± 8.5 years, 9 men) were scanned with coronary CT angiography and 4D flow MRI one month after CABG surgery. Grafts included 10 left internal mammary arteries (LIMA), 3 radial arteries (RA), and 11 saphenous vein grafts (SVG). Image-guided CFD was used to quantify blood flow rate and wall area exposed to abnormal wall shear stress (WSS). Arterial grafts had a lower abnormal WSS area than venous grafts (17.9% vs. 70.1%; p = 0.001), and a similar trend was observed for LIMA vs. SVG (13.8% vs. 70.1%; p = 0.001). Abnormal WSS area correlated positively to lumen diameter (p

In this paper, a generalized macromodeling approach is presented to simulate complex electromagnetic (EM) surfaces consisting of unit cells with connected conductors. Macromodels of each unit cell are produced by applying the equivalence principle on fictitious surfaces encapsulating them. Unit cells often consist of multiple dielectric layers and conductor traces, featuring multiscale structures. Challenges arise when a current-carrying conductor trace traverses the fictitious surface. Hence, a new method based on half Rao-Wilton-Glisson basis functions is proposed to accurately ensure the continuity of the surface currents and avoid singularities at the intersections. The accuracy of the proposed approach is validated by comparing the results with commercial solvers for different EM surfaces.

Surface integral equation (SIE) methods are of great interest for the efficient electromagnetic modeling of various devices, from integrated circuits to antenna arrays. Existing acceleration algorithms for SIEs, such as the adaptive integral method (AIM), enable the fast approximation of interactions between well-separated mesh elements. Nearby interactions involve the singularity of the kernel, and must instead be computed accurately with direct integration at each frequency of interest, which can be computationally expensive. We propose a novel algorithm for reducing the cost-per-frequency of near-region computations for both homogeneous and layered background media. In the proposed extended AIM (AIMx), the SIE operators are decomposed into a frequency-independent term containing the singularity of the kernel, and a nonsingular frequency-dependent term. Direct integration is only required for the frequency-independent term, and can be reused at each frequency, leading to significantly faster frequency sweeps. The frequency-dependent term is captured with good accuracy via fast Fourier transform (FFT)-based acceleration even in the near region, as confirmed with an error analysis. The accuracy and efficiency of the proposed method are demonstrated through numerical examples drawn from several applications, and CPU times are significantly reduced by factors ranging from three to 16., 15 pages, 18 figures, accepted for publication in the IEEE Transactions on Antennas and Propagation

Coronary artery bypass graft surgery is a common intervention for coronary artery disease; however, it suffers from graft failure, and the underlying mechanisms are not fully understood. To better understand the relation between graft hemodynamics and surgical outcomes, we performed computational fluid dynamics simulations with deformable vessel walls in 10 study participants (24 bypass grafts) based on CT and 4D flow MRI one month after surgery to quantify lumen diameter, wall shear stress (WSS), and related hemodynamic measures. A second CT acquisition was performed one year after surgery to quantify lumen remodeling. Compared to venous grafts, left internal mammary artery grafts experienced lower abnormal WSS (

The electromagnetic modeling of dielectrics is necessary in several radiation and scattering problems, such as those involving dielectric antennas and metasurfaces. The surface integral method, accelerated with algorithms such as the adaptive integral method (AIM), require computing interactions between nearby mesh elements with direct integration at each frequency point, which can be computationally expensive. In this work, we propose an extended AIM (AIMx) for the fast analysis of dielectric media, which is easy to implement within existing AIM-based code. The surface integral operators are decomposed into frequency-independent and frequency-dependent parts such that direct integration is only necessary for the frequency-independent part, which can be performed once and reused for each frequency point. Several realistic numerical examples demonstrate significant CPU time savings with no loss of accuracy.

The choice of appropriate boundary conditions is a fundamental step in computational fluid dynamics (CFD) simulations of the cardiovascular system. Boundary conditions, in fact, highly affect the computed pressure and flow rates, and consequently haemodynamic indicators such as wall shear stress, which are of clinical interest. Devising automated procedures for the selection of boundary conditions is vital to achieve repeatable simulations. However, the most common techniques do not automatically assimilate patient-specific data, relying instead on expensive and time-consuming manual tuning procedures. In this work, we propose a technique for the automated estimation of outlet boundary conditions based on optimal control. The values of resistive boundary conditions are set as control variables and optimized to match available patient-specific data. Experimental results on four aortic arches demonstrate that the proposed framework can assimilate 4D-Flow MRI data more accurately than two other common techniques based on Murray's law and Ohm's law. This article is protected by copyright. All rights reserved.

A new boundary element formulation is presented for the electromagnetic analysis of interconnects in layered media, which accurately captures skin effect in conductors of arbitrary shape and size. Results show superior CPU-time and memory performance compared to several state-of-the-art techniques.

A novel parallel solver based on the adaptive integral method (AIM) is proposed for the electromagnetic analysis of electrical interconnects in layered media. We show that graph partitioning techniques can be used to optimally distribute, across thousands of processes, the computations related to both matrix filling and system solution. The proposed workload distribution strategy is compared to existing techniques through a scalability study on a large realistic interposer model in layered media.

A novel single-source surface integral formulation is proposed for efficient modeling of penetrable media, without the need for dual basis functions. The differential current density concept is leveraged to yield a well conditioned system matrix. The new formulation leads to significantly reduced GMRES iterations compared to existing single-source formulations, over a broad range of frequency and material parameters.