Variation Of Wall Shear Stress Research Articles

Research on heat transfer characteristics and prediction methods of supercritical liquified natural gas flowing horizontally under the influence of buoyancy

An exhaustive elucidation of the heat transfer performances and mechanisms of supercritical LNG under water-bath heating conditions is paramount for optimizing heat transfer and designing related heat exchangers. This study employs experiments and numerical simulations to explore the heat transfer mechanisms and prediction methods of supercritical LNG flowing horizontally under buoyancy. The experimental investigation was executed within a small SCV using the N2-Ar mixture as a substitute medium for safety considerations. Concurrently, the accuracy of the numerical model employed is rigorously authenticated through a comparative analysis with the empirical findings obtained from the experimental endeavor. The investigation delves into the impacts of operating conditions, structural parameters, and medium composition on heat transfer through numerical simulations. The influence mechanism of buoyancy on heat transfer is meticulously explicated by examining the distribution of flow fields and thermophysical properties within the boundary layer. The introduction of a novel criterion, Bufull, derived from the analysis of wall shear stress variation caused by buoyancy, facilitates a quantitative assessment of the impact of buoyancy on heat transfer. Additionally, a simplified buoyancy criterion, FBu, is introduced, dependent solely on fundamental operating parameters and tailored for engineering applications, with a determined critical value of 1/ FBu = 60.3. Based on the dramatically variable thermophysical and transport properties, that occurs near the critical point, as well as buoyancy correction, a correlation is formulated for predicting mixed convection heat transfer, exhibiting an MAE of merely 3.04% and accurately predicting 81.28% of the data within a ± 5% margin of error, surpassing existing correlations. Finally, the proposed heat transfer calculation method is validated using five groups of operational data from two in-service SCVs. The maximum absolute discrepancy between the computed outlet temperature and the monitored value is 1.95K. The findings of this paper can support practical applications.

Numerical study of Casson dusty hybrid nanofluid flow with volume fraction over a Darcy–Forchheimer porous medium

AbstractThis study develops a mathematical model for dusty Casson hybrid nanofluid flow over a vertical stretching sheet, considering the influence of volume fractions for the dust particles and nanoparticles, magnetic field and Darcy–Forchheimer penetrating medium. The nanoparticles used are TC4 (Ti‐6Al‐4V) and Nichrome (80% Ni and 20% Cr), suspended in unused engine oil (EO). The governing partial differential equations are transformed into ordinary ones using suitable transforms. Solutions are obtained numerically via the Matlab built‐in program Bvp4c technique and represent graphically the velocity, temperature, and concentration profiles for both the fluid and dust phases. Additionally, tabular forms present wall shear stress and heat transfer rate variations under different parameter values. It is found that the enhancing volume fraction parameter helps to boost the magnitude of shear stress, heat, and mass diffusivity.

Relative Residence Time Can Account for Half of the Anatomical Variation in Fatty Streak Prevalence Within the Right Coronary Artery.

The patchy anatomical distribution of atherosclerosis has been attributed to variation in haemodynamic wall shear stress (WSS). The consensus is that low WSS and a high Oscillatory Shear Index (OSI) trigger the disease. We found that atherosclerosis at aortic branch sites correlates threefold better with transverse WSS (transWSS), a metric which quantifies multidirectional near-wall flow. Coronary artery disease has greater clinical significance than aortic disease but computation of WSS metrics is complicated by the substantial vessel motion occurring during each cardiac cycle. Here we present the first comparison of the distribution of atherosclerosis with WSS metrics computed for moving coronary arteries. Maps of WSS metrics were computed using dynamic geometries reconstructed from angiograms of ten non-stenosed human right coronary arteries (RCAs). They were compared with maps of fatty streak prevalence derived from a previous study of 1852 RCAs. Time average WSS (TAWSS), OSI, transWSS and the cross-flow index (CFI), a non-dimensional form of the transWSS, gave non-significant or significant but low spatial correlations with lesion prevalence. The highest correlation coefficient (0.71) was for the relative residence time (RRT), a metric that decreases with TAWSS and increases with OSI. The coefficient was not changed if RRT was calculated using CFI, which captures multidirectional WSS only, rather than OSI, which encompasses both multidirectional and oscillatory WSS. Contrary to our earlier findings in the aorta, low WSS in combination with highly multidirectional flow correlates best with lesion location in the RCA, explaining approximately half of its anatomical variation.

Numerical prediction of passive speed and performance for multistage pump without power drive in natural flow process

With the development of engineering applications and the increase in system complexity, some particular fields, such as liquid rocket engine turbopumps, aircraft engine fuel systems, and marine natural flow cooling systems, are increasingly focusing on the performance characteristics of pumps under natural flow conditions. The pump is in the form of resistance components under natural flow conditions without a power drive. The impeller undergoes passive rotation by the impact of inlet flow. Due to the specificity of its operating conditions and performance indicators, the pump's natural flow performance cannot be evaluated by regular methods. Therefore, this paper proposed a numerical prediction method for pump natural flow performance based on a coupled computational fluid dynamics coupled with six-degrees-of freedom model. The performance of a multistage pump with guide vanes was evaluated under different natural flow conditions, and the accuracy was verified by experimental measurements. The transient variation mode of pump performance parameters with time at the initial stage of natural flow impact was analyzed. The flow field's transient evolution characteristics and the wall shear stress variation during natural flow were investigated. It was found that the impeller's passive rotational speed increases linearly with the natural flow rate, while the hydraulic loss shows an exponentially increasing trend. Meanwhile, the natural flow loss coefficient shows an exponentially decreasing trend and gradually tends to a stable value. The high turbulent kinetic energy inside the impeller is mainly distributed in the flow separation region and large velocity gradients. The distribution of shear stresses is closely related to the flow behavior inside the pump and the geometrical features of the hydraulic components.

FSGe: A fast and strongly-coupled 3D fluid–solid-growth interaction method

Equilibrated fluid–solid-growth (FSGe) is a fast, open source, three-dimensional (3D) computational platform for simulating interactions between instantaneous hemodynamics and long-term vessel wall adaptation through mechanobiologically equilibrated growth and remodeling (G&R). Such models can capture evolving geometry, composition, and material properties in health and disease and following clinical interventions. In traditional G&R models, this feedback is modeled through highly simplified fluid solutions, neglecting local variations in blood pressure and wall shear stress (WSS). FSGe overcomes these inherent limitations by strongly coupling the 3D Navier–Stokes equations for blood flow with a 3D equilibrated constrained mixture model (CMMe) for vascular tissue G&R. CMMe allows one to predict long-term evolved mechanobiological equilibria from an original homeostatic state at a computational cost equivalent to that of a standard hyperelastic material model. In illustrative computational examples, we focus on the development of a stable aortic aneurysm in a mouse model to highlight key differences in growth patterns between FSGe and solid-only G&R models. We show that FSGe is especially important in blood vessels with asymmetric stimuli. Simulation results reveal greater local variation in fluid-derived WSS than in intramural stress (IMS). Thus, differences between FSGe and G&R models became more pronounced with the growing influence of WSS relative to pressure. Future applications in highly localized disease processes, such as for lesion formation in atherosclerosis, can now include spatial and temporal variations of WSS.

A Pseudo-Spectral Method for Wall Shear Stress Estimation from Doppler Ultrasound Imaging in Coronary Arteries.

The Wall Shear Stress (WSS) is the component tangential to the boundary of the normal stress tensor in an incompressible fluid, and it has been recognized as a quantity of primary importance in predicting possible adverse events in cardiovascular diseases, in general, and in coronary diseases, in particular. The quantification of the WSS in patient-specific settings can be achieved by performing a Computational Fluid Dynamics (CFD) analysis based on patient geometry, or it can be retrieved by a numerical approximation based on blood flow velocity data, e.g., ultrasound (US) Doppler measurements. This paper presents a novel method for WSS quantification from 2D vector Doppler measurements. Images were obtained through unfocused plane waves and transverse oscillation to acquire both in-plane velocity components. These velocity components were processed using pseudo-spectral differentiation techniques based on Fourier approximations of the derivatives to compute the WSS. Our Pseudo-Spectral Method (PSM) is tested in two vessel phantoms, straight and stenotic, where a steady flow of 15mL/min is applied. The method is successfully validated against CFD simulations and compared against current techniques based on the assumption of a parabolic velocity profile. The PSM accurately detected Wall Shear Stress (WSS) variations in geometries differing from straight cylinders, and is less sensitive to measurement noise. In particular, when using synthetic data (noise free, e.g., generated by CFD) on cylindrical geometries, the Poiseuille-based methods and PSM have comparable accuracy; on the contrary, when using the data retrieved from US measures, the average error of the WSS obtained with the PSM turned out to be 3 to 9 times smaller than that obtained by state-of-the-art methods. The pseudo-spectral approach allows controlling the approximation errors in the presence of noisy data. This gives a more accurate alternative to the present standard and a less computationally expensive choice compared to CFD, which also requires high-quality data to reconstruct the vessel geometry.

Numerical investigation on flow-induced wall shear stress variation of metastatic cancer cells in lymphatics with elastic valves

Hematogenous metastasis occurs when cancer cells detach from the extracellular matrix in the primary tumor into the bloodstream or lymphatic system. Elucidating the response of metastatic tumor cells in suspension to the flow conditions in lymphatics with valves from a mechanical/fluidic perspective is necessary. A physiologically relevant computational model of a lymphatic vessel with valves was constructed using fully coupled fluid–cell–vessel interactions to investigate the effects of lymphatic vessel contractility, valve properties, and cell size and stiffness on the variations in magnitude and gradient of the flow-induced wall shear stress (WSS) experienced by suspended tumor cells. Results indicated that the maximum WSSmax increased with the increments in cell diameter, vessel contraction amplitude, and valve stiffness. The decrease in vessel contraction period and valve aspect ratio also increased the maximum WSSmax. The influence of the properties of the valve on the WSS was more significant among the factors mentioned above. The maximum WSSmax acting on the cancer cell when the cell reversed the direction of its motion in the valve region increased by 0.5–1.4 times that before the cell entered the valve region. The maximum change in WSS was in the range of 0.004–0.028 Pa/µm depending on the factors studied. They slightly exceeded the values associated with breast cancer cell apoptosis. The results of this study provide biofluid mechanics-based support for mechanobiological research on the metastasis of metastatic cancer cells in suspension within the lymphatics.

Unsteady electroosmotic flow of Carreau–Newtonian fluids through a cylindrical tube

The present study aims to investigate the fully developed electroosmotic flow of Carreau–Newtonian fluids within a circular tube, considering two different boundary conditions of zeta potential. The unsteady behavior of Carreau–Newtonian fluids, is attributed to the pulsatile pressure-driven flow. Studying two distinct formulations, slip and no-slip zeta potentials, illustrates the relative movement of fluid particles and the charged surface. The proposed study’s framework is divided into two regions: a central region containing a non-Newtonian Carreau fluid exhibiting both shear-thinning and thickening behavior, and a peripheral region containing a Newtonian fluid with the effect of an electrical double layer at the solid wall of the cylindrical tube. The Poisson–Boltzmann equation under the assumption of Debye–Hückel approximation governs the potential distribution due to the electric double layer, facilitating the examination of electroosmotic flow through tube. Obtaining analytical solutions for the momentum equations in both regions becomes challenging due to the inclusion of pulsatile pressure-driven flow and a non-Newtonian Carreau fluid. Asymptotic series expansions are employed, utilizing small parameters such as the Weissenberg number (We≪1) and a pulsatile Reynolds number (α≪1), to derive velocity expressions for both regions. Wall shear stress fluctuations, as indicated by time-averaged wall shear stress and oscillatory shear index are assessed to gauge the temporal variation of wall shear stress from the dominant blood flow direction. By utilizing potential and hydrodynamic variables, an asymptotic formulation for the streaming potential is developed to ascertain the asymptotic expression of electrokinetic energy conversion efficiency and scrutinize the impacts of pertinent physical parameters on it. The proposed study prominently shows that the Debye–Hückel parameter, Weissenberg number, pulsatile Reynolds number, and time-dependent pressure gradient significantly influence the fluid velocity, flow rate, and flow resistance. The notable determination of the present study is that the velocity amplitude and electrokinetic energy conversion efficiency are enlarged in the slip-dependent case relative to the no-slip zeta potential. It is noteworthy that as inertial forces become more prominent compared to viscous forces, both time-average wall shear stress and oscillatory shear index experience a slight decrease, although the reduction in oscillatory shear index is minimal owing to the absence of obstruction within the tube wall. The NDSolve numerical scheme in Mathematica software is employed to calculate numerical solutions for the governing equations, which are subsequently verified against results obtained through asymptotic analysis. The findings of this study are expected to deepen our understanding of unsteady electroosmotic flows of Carreau fluid and contribute to the design and development of advanced microfluidic devices with enhanced performance for various applications.

Mathematical modeling of creeping electromagnetohydrodynamic peristaltic propulsion in an annular gap between sinusoidally deforming permeable and impermeable curved tubes

The present work investigates the creeping peristaltic propulsion of viscid fluid in an annular gap between sinusoidally deforming permeable and impermeable curved tubes of similar shape under the influence of an externally imposed electric and magnetic field. In this model, the outer tube with a permeable wall surface is supposed to satisfy the Saffman slip condition. The flow equations are simplified by the estimation of a large wavelength in comparison with the radius of the external tube. An analytical solution for the axial velocity is obtained in the computational software MATHEMATICA. Graphical analyses are conducted to explore the variations in wall shear stress, velocity, pressure rise, frictional force, and stream function with respect to different emergent parameters, providing insight into the underlying physics of the flow phenomena. An investigation of the effects of the Hartmann number and electric field strength on the flow through a gap between deformable tubes with curved structures has important implications for a variety of engineering applications, including mechanical and biomedical engineering. The streamlines are plotted to discuss fluid trapping and visualize the flow pattern of the viscid fluid inside the curved annular domain. A comparative analysis of fluid transport induced by sinusoidal, triangular, trapezoidal, and square wave shapes is encountered with the help of streamlined contour diagrams. The comparison of pressure gradients in three different models is also discussed to gain insight due to fluid–structure interaction. A gap in the body of recently published literature is filled by the results discussed in this paper.

The impact of choledochal cysts on bile fluid dynamics: A perspective using computational fluid dynamics and surface mapping technique

Choledochal cysts (CCs) are an important risk factor for cholangiocarcinoma, though their etiology remains debated. Given the vital role of bile fluid in digestive processes within the biliary system, examining such mechanisms from the perspective of bile fluid dynamics may offer additional insights for clinical use. This study utilized magnetic resonance imaging (MRI)-based patient-specific scans for detailed reconstruction and further employed the computational fluid dynamic method to assess the physiological functions of each system, including refilling and emptying processes. The impact of bile rheological property was also examined. Key biomechanical parameters—pressure and wall shear stress (WSS)—were displayed on a two-dimensional plane via surface mapping for enhanced visualization and comparative analysis. Outcomes demonstrated a significant reduction in bile flow velocity in CCs patients due to common bile duct's anatomical features and bile's shear-thinning, non-Newtonian nature, with a notable increase in pressure drop observed. In healthy biliary systems, WSS variations were minimal; however, in CCs patients, extreme WSS differences were found, with the highest WSS in the segmental bile duct and the lowest in the dilatation area, presenting a magnitude difference of approximately 1000. CCs one showed WSS levels 100–250 times higher than healthy ones in the common bile duct. Bile rheological properties substantially affect pressure and WSS patterns, particularly WSS, where pathological bile caused a tenfold increase in WSS compared to healthy bile. The findings aimed to enhance the understanding of biliary fluid mechanics in CCs and offer insights into selected fluidic variables for future microfluidic chip experiments.

Evaluation of the shrinkage process of a neck remnant after stent-coil treatment of a cerebral aneurysm using silent magnetic resonance angiography and computational fluid dynamics analysis: illustrative case.

Silent magnetic resonance angiography (MRA) mitigates metal artifacts, facilitating clear visualization of neck remnants after stent and coil embolization of cerebral aneurysms. This study aims to scrutinize hemodynamics at the neck remnant by employing silent MRA and computational fluid dynamics. The authors longitudinally tracked images of a partially thrombosed anterior communicating artery aneurysm's neck remnant, which had been treated with stent-assisted coil embolization, using silent MRA over a decade. Computational fluid dynamics delineated the neck remnant's reduction process, evaluating hemodynamic parameters such as flow rate, wall shear stress magnitude and vector, and streamlines. The neck remnant exhibited diminishing surface area, volume, neck size, dome depth, and aspect ratio. Its reduction correlated with a decline in the flow rate ratio of the remnant dome to the inflow parent artery. Analysis delineated regions within the contracting neck remnant characterized by consistently low average wall shear stress magnitude and variation, accompanied by notable variations in wall shear stress vector directionality. Evaluation of neck remnants after stent-coil embolization is possible through silent MRA and computational fluid dynamics. Predicting the neck remnant reduction may be achievable through hemodynamic parameter analysis.

Wall shear stress measurement of turbulent bubbly flows using laser Doppler displacement sensor

A novel wall shear stress meter is developed to promote studies of frictional drag reduction by injecting air bubbles into a turbulent boundary layer. Such studies require a tool having a time resolution sufficiently high to capture the unsteady effect of turbulent characteristics due to various bubble passage patterns. The present study adopts the principle of the heterodyne interferometer to capture of time variation of wall shear stress in turbulent flow up to 500 Hz in sampling frequency. We demonstrate the performance of the tool by applying it to bubbly flows beneath a flat-bottom model ship running at 3 m/s in a towing water reservoir, where bubbles experience wavy passage around a few Hz in frequency. The data obtained by the proposed shear stress meter agree well with Dean's formula for single-phase turbulent channel flow. The scattering of the data across three trials is sufficiently small, at approximately 3 %. This high reproducibility is a great advantage over conventional direct measurements, which have low reproducibility and a large scatter of the data of approximately 15 %. The obtained results represent reasonable drag reductions by injecting bubbles with different void fractions. We also evaluate histogram of the wall shear stress fluctuation to characterize the effect of drag reduction by unsteady passage of bubbles.

A Parametric 3D Model of Human Airways for Particle Drug Delivery and Deposition

The treatment for asthma and chronic obstructive pulmonary disease relies on forced inhalation of drug particles. Their distribution is essential for maximizing the outcomes. Patient-specific computational fluid dynamics (CFD) simulations can be used to optimize these therapies. In this regard, this study focuses on creating a parametric model of the human respiratory tract from which synthetic anatomies for particle deposition analysis through CFD simulation could be derived. A baseline geometry up to the fourth generation of bronchioles was extracted from a CT dataset. Radial basis function (RBF) mesh morphing acting on a dedicated tree structure was used to modify this baseline mesh, extracting 1000 synthetic anatomies. A total of 26 geometrical parameters affecting branch lengths, angles, and diameters were controlled. Morphed models underwent CFD simulations to analyze airflow and particle dynamics. Mesh morphing was crucial in generating high-quality computational grids, with 96% of the synthetic database being immediately suitable for accurate CFD simulations. Variations in wall shear stress, particle accretion rate, and turbulent kinetic energy across different anatomies highlighted the impact of the anatomical shape on drug delivery and deposition. The study successfully demonstrates the potential of tree-structure-based RBF mesh morphing in generating parametric airways for drug delivery studies.

Computational fluid dynamics study on three‐dimensional polymeric scaffolds to predict wall shear stress using machine learning models for bone tissue engineering applications

AbstractGeometrical patterns and dimensions of the polymeric scaffold play a major role in controlling the degradation and mechanical stimuli for osteogenic differentiation. Wall shear stress (WSS) analysis of scaffold provides a better understanding of the body fluid flow dynamics. A computational fluid dynamics (CFD) study was carried out to understand velocity profile and WSS distribution when the strands are arranged in rectangular and triangular pitch for the different strand diameters and spacing. The number of scaffold surfaces with less than 30 mPa and maximum and average WSS was estimated to check the suitability of the scaffold for loading stem cells. This situation is favorable to induce osteogenic activity and cell viability. Higher spacing/pitch between the strands increases the chances of scaffold surface having WSS less than 30 mPa. When the spacing and diameter are smaller, there is no significant variation in WSS and pressure drop between rectangular and triangular pitch arrangement is observed. Machine learning (ML) models were developed to predict WSS distribution and to reduce the computational cost involved in solving the Navier–Stokes equation. XG Boost and support vector machine (SVM) models outperform the other models in predicting the WSS with high R2 and five‐fold cross‐validation accuracy and are helpful in predicting the optimal design parameters of a three‐dimensional scaffold.

Abstract 18625: Mechano-Transcriptomic Interaction Mediates Spatial Variations in Early Vascular Fibrosis

Background: The spatiotemporal variations of wall shear stress (WSS) modulate vascular genomics and function. Despite the advances in imaging modalities and computational modeling, the link between local hemodynamics and genetic pathways underlying atherosclerosis development remains elusive. We hypothesized that differential WSS patterns in the aorta modulate fibroblast gene expression in response to high fat diet (HFD) leading to vascular stiffness. Methods: We performed ultrasound-based 4-D+time moving boundary computational fluid dynamics (CFD) on the aorta of Ldlr -/- mice fed a normal diet (ND) or HFD during 4 weeks. B-mode images (Vevo 3100) of the aortic arch were used to reconstruct a 3-D model of the lumen structure by revolving the vessel wall around the centerline in the aorta and the branches. Computational domain reconstruction was followed by image registration to generate a moving boundary model. WSS contours and velocity profiles were acquired by solving Navier-Stokes equations using svFSI code. The longitudinal strains along the outer curvature of the aortic arch over several cardiac cycles were measured, and the changes were compared with the effects of HFD. Mouse aortas were processed for single cell RNA sequencing and analyzed using the Scanpy package. Results and Conclusion: Our analysis of transcriptomics profiles revealed increased Acta2 and Postn expression in the myofibroblast subclusters, suggesting myofibroblast differentiation in response to HFD in the aortic arch. This was confirmed and localized to the lesser curvature of the aortic arch through Periostin immunostaining. We found that the increase in Postn and Acta2 is associated with reduced arterial strain or increased arterial stiffness measured in the aortic arch. In parallel, CFD analysis demonstrated that an elevated WSS along the outer curvature vs. low WSS and disturbed flow along the inner curvature in early diastole results in differential gene expression to regulate myofibroblast differentiation and modulates vascular stiffness.

Detection of Abnormal Wall Shear Stress and Oscillatory Shear Index via Ultrasound Vector Flow Imaging as Possible Indicators for Arteriovenous Fistula Stenosis in Hemodialysis.

The arteriovenous fistula (AVF) is an essential vascular access for hemodialysis patients. AVF stenosis may occur at sites with abnormal wall shear stress (WSS) and oscillatory shear index (OSI), which are caused by the complex flow in the AVF. At present, an effective method for rapid determination of the WSS and OSI of the AVF is lacking. The objective of this study was to apply an ultrasound-based method for determination of the WSS and OSI to explore the risk sites of the AVF. In this study, the ultrasound vector flow imaging technique V Flow was applied to measure the WSS and OSI at four different regions of the AVF to detect and analyze the risk sites: (i) anastomosis region, (ii) curved region, (iii) proximal vein and (iv) distal vein. Twenty-one patients were included in this study. The relative residence time was calculated based on the measured WSS and OSI. The curved region had the lowest WSS; the anastomosis region had a significantly higher OSI (p < 0.05) compared with the venous regions, and the curved region had a significantly higher RRT (p < 0.05) compared with the proximal vein region. V Flow is a feasible tool for studying WSS variations in AVF. The possible risk site in the AVF may be located in the anastomosis and curved regions, where the latter could present a higher risk for AVF stenosis.

The variation in wall shear stress induced by a linear train of soft particles in channel flow

We investigate the wall shear stress variation induced by soft particles modeled as capsules migrating in a channel. Interestingly, the wall shear stress exhibits a roughly linear increase in both its global maximum value and variation magnitude with an increase in the normalized overall lateral position of the capsules when they stabilize as a linear train. Furthermore, when a single capsule stabilizes in the channel centerline, the streamlines show an M-like shape in the upper part and a W-like shape in the lower part of the channel. Meanwhile, we use the vertical velocity to explain the typical peak-valley-peak structure in the wall shear stress profile. For multiple capsules, the contours of the vertical velocity can also be employed to determine the locations where the peaks or valleys in the wall shear stress occur. These findings enhance our comprehension of the variation in wall shear stress caused by soft particles.

Is spontaneous coronary artery dissection (SCAD) related to local anatomy and hemodynamics? An exploratory study

AimsSpontaneous coronary artery dissection (SCAD) is an increasingly diagnosed cause of myocardial infarction with unclear pathophysiology. The aim of the study was to test if vascular segments site of SCAD present distinctive local anatomy and hemodynamic profiles. MethodsCoronary arteries with spontaneously healed SCAD (confirmed by follow-up angiography) underwent three-dimensional reconstruction, morphometric analysis with definition of vessel local curvature and torsion, and computational fluid dynamics (CFD) simulations with derivation of time-averaged wall shear stress (TAWSS) and topological shear variation index (TSVI). The (reconstructed) healed proximal SCAD segment was visually inspected for co-localization with curvature, torsion, and CFD-derived quantities hot spots. ResultsThirteen vessels with healed SCAD underwent the morpho-functional analysis. Median time between baseline and follow-up coronary angiograms was 57 (interquartile range [IQR] 45–95) days. In seven cases (53.8%), SCAD was classified as type 2b and occurred in the left anterior descending artery or near a bifurcation. In all cases (100%), at least one hot spot co-localized within the healed proximal SCAD segment, in 9 cases (69.2%) ≥ 3 hot spots were identified. Healed SCAD in proximity of a coronary bifurcation presented lower TAWSS peak values (6.65 [IQR 6.20–13.20] vs. 3.81 [2.53–5.17] Pa, p = 0.008) and hosted less frequently TSVI hot spots (100% vs. 57.1%, p = 0.034). ConclusionVascular segments of healed SCAD were characterized by high curvature/torsion and WSS profiles reflecting increased local flow disturbances. Hence, a pathophysiological role of the interaction between vessel anatomy and shear forces in SCAD is hypothesized.

Detection of Early Endothelial Dysfunction by Optoacoustic Tomography.

Variations in vascular wall shear stress are often presumed to result in the formation of atherosclerotic lesions at specific arterial regions, where continuous laminar flow is disturbed. The influences of altered blood flow dynamics and oscillations on the integrity of endothelial cells and the endothelial layer have been extensively studied in vitro and in vivo. Under pathological conditions, the Arg-Gly-Asp (RGD) motif binding integrin αvβ3 has been identified as a relevant target, as it induces endothelial cell activation. Animal models for in vivo imaging of endothelial dysfunction (ED) mainly rely on genetically modified knockout models that develop endothelial damage and atherosclerotic plaques upon hypercholesterolemia (ApoE-/- and LDLR-/-), thereby depicting late-stage pathophysiology. The visualization of early ED, however, remains a challenge. Therefore, a carotid artery cuff model of low and oscillating shear stress was applied in CD-1 wild-type mice, which should be able to show the effects of altered shear stress on a healthy endothelium, thus revealing alterations in early ED. Multispectral optoacoustic tomography (MSOT) was assessed as a non-invasive and highly sensitive imaging technique for the detection of an intravenously injected RGD-mimetic fluorescent probe in a longitudinal (2-12 weeks) study after surgical cuff intervention of the right common carotid artery (RCCA). Images were analyzed concerning the signal distribution upstream and downstream of the implanted cuff, as well as on the contralateral side as a control. Subsequent histological analysis was applied to delineate the distribution of relevant factors within the carotid vessel walls. Analysis revealed a significantly enhanced fluorescent signal intensity in the RCCA upstream of the cuff compared to the contralateral healthy side and the downstream region at all time points post-surgery. The most obvious differences were recorded at 6 and 8 weeks after implantation. Immunohistochemistry revealed a high degree of αv-positivity in this region of the RCCA, but not in the left common carotid artery (LCCA) or downstream of the cuff. In addition, macrophages could be detected by CD68 immunohistochemistry in the RCCA, showing ongoing inflammatory processes. In conclusion, MSOT is capable of delineating alterations in endothelial cell integrity in vivo in the applied model of early ED, where an elevated expression of integrin αvβ3 was detected within vascular structures.

Oscillatory behaviors of multiple shock waves to upstream disturbances

The oscillatory response of multiple shock waves to upstream disturbances in a supersonic flow is studied numerically in a constant area rectangular duct. The flow is accelerated through a nozzle with an exit Mach number of 1.75 and continues in the constant area duct, where multiple shock waves are formed. To investigate the effect of upstream disturbance on shock oscillations, three parameters are varied systematically: upstream turbulent intensity, frequency of upstream pressure fluctuation, and amplitude of upstream pressure fluctuation. The wall shear stress variation along the duct length provides the location of separation and reattachment points in the flow field. The wall pressure frequency spectra were used to investigate the low-frequency unsteadiness in shock oscillations. The power spectral density of the wall static pressure and the probability density function (PDF) of shock location are analyzed, and the results suggest that as the upstream turbulent intensity is increased, the dominant frequency of oscillation is increased and the shock oscillations become more symmetrical. As the upstream disturbance frequency is increased, the shock oscillations become more symmetrical and follow the Gaussian curve closely. The shock wave oscillates with the same upstream excitation frequency when the upstream disturbance amplitude is increased. At large values of upstream disturbance amplitude, the PDF shows a large deviation from the Gaussian, and the rms amplitude of shock oscillation increases monotonously. At higher amplitudes of upstream disturbance excitation, the traces of shock train leading-edge location display path-dependence characteristics.