Persistent corotation of the large-scale flow of thermal convection and an immersed free body.

Inspired by the superrotation of the Earth's solid core, we investigate the dynamics of a free-rotating body as it interacts with the large-scale circulation (LSC) of the Rayleigh-Bénard thermal convection in a cylindrical container. A surprising and persistent corotation of both the free body and the LSC emerges, breaking the axial symmetry of the system. The corotational speed increases monotonically with the intensity of thermal convection, measured by the Rayleigh number Ra, which is proportional to the temperature difference between the heated bottom and cooled top. The rotational direction occasionally and spontaneously reverses, occurring more frequently at higher Ra. The reversal events follow a Poisson process; it is feasible that flow fluctuations randomly interrupt and reestablish the rotation-sustaining mechanism. This corotation is powered by thermal convection alone and promoted by the addition of a free body, enriching the classical dynamical system.

Dissolution-driven propulsion of floating solids.

We show that unconstrained asymmetric dissolving solids floating in a fluid can move rectilinearly as a result of attached density currents which occur along their inclined surfaces. Solids in the form of boats composed of centimeter-scale sugar and salt slabs attached to a buoy are observed to move rapidly in water with speeds up to 5 mm/s determined by the inclination angle and orientation of the dissolving surfaces. While symmetric boats drift slowly, asymmetric boats are observed to accelerate rapidly along a line before reaching a terminal velocity when their drag matches the thrust generated by dissolution. By visualizing the flow around the body, we show that the boat velocity is always directed opposite to the horizontal component of the density current. We derive the thrust acting on the body from its measured kinematics and show that the propulsion mechanism is consistent with the unbalanced momentum generated by the attached density current. We obtain an analytical formula for the body speed depending on geometry and material properties and show that it captures the observed trends reasonably. Our analysis shows that the gravity current sets the scale of the body speed consistent with our observations, and we estimate that speeds can grow slowly as the cube root of the length of the inclined dissolving surface. The dynamics of dissolving solids demonstrated here applies equally well to solids undergoing phase change and may enhance the drift of melting icebergs, besides unraveling a primal strategy by which to achieve locomotion in active matter.

Physical mechanisms of red blood cell splenic filtration.

The splenic interendothelial slits fulfill the essential function of continuously filtering red blood cells (RBCs) from the bloodstream to eliminate abnormal and aged cells. To date, the process by which 8 [Formula: see text]m RBCs pass through 0.3 [Formula: see text]m-wide slits remains enigmatic. Does the slit caliber increase during RBC passage as sometimes suggested? Here, we elucidated the mechanisms that govern the RBC retention or passage dynamics in slits by combining multiscale modeling, live imaging, and microfluidic experiments on an original device with submicron-wide physiologically calibrated slits. We observed that healthy RBCs pass through 0.28 [Formula: see text]m-wide rigid slits at 37 °C. To achieve this feat, they must meet two requirements. Geometrically, their surface area-to-volume ratio must be compatible with a shape in two tether-connected equal spheres. Mechanically, the cells with a low surface area-to-volume ratio (28% of RBCs in a 0.4 [Formula: see text]m-wide slit) must locally unfold their spectrin cytoskeleton inside the slit. In contrast, activation of the mechanosensitive PIEZO1 channel is not required. The RBC transit time through the slits follows a [Formula: see text]1 and [Formula: see text]3 power law with in-slit pressure drop and slip width, respectively. This law is similar to that of a Newtonian fluid in a two-dimensional Poiseuille flow, showing that the dynamics of RBCs is controlled by their cytoplasmic viscosity. Altogether, our results show that filtration through submicron-wide slits is possible without further slit opening. Furthermore, our approach addresses the critical need for in vitro evaluation of splenic clearance of diseased or engineered RBCs for transfusion and drug delivery.

Magnetohydrodynamic levitation for high-performance flexible pumps.

We use magnetohydrodynamic levitation as a means to create a soft, elastomeric, solenoid-driven pump (ESP). We present a theoretical framework and fabrication of a pump designed to address the unique challenges of soft robotics, maintaining pumping performance under deformation. Using a permanent magnet as a piston and ferrofluid as a liquid seal, we model and construct a deformable displacement pump. The magnet is driven back and forth along the length of a flexible core tube by a series of solenoids made of thin conductive wire. The magnet piston is kept concentric within the tube by Maxwell stresses within the ferrofluid and magnetohydrodynamic levitation, as viscous lift pressure is created due to its forward velocity. The centering of the magnet reduces shear stresses during pumping and improves efficiency. We provide a predictive model and capture the transient nonlinear dynamics of the magnet during operation, leading to a parametric performance curve characterizing the ESP, enabling goal-driven design. In our experimental validation, we report a shut-off pressure of 2 to 8 kPa and run-out flow rate of 50 to 320 mL⋅min-1, while subject to deformation of its own length scale, drawing a total of 0.17 W. This performance leads to the highest reported duty point (i.e., pressure and flow rate provided under load) for a pump that operates under deformation of its own length scale. We then integrate the pump into an elastomeric chassis and squeeze it through a tortuous pathway while providing continuous fluid pressure and flow rate; the vehicle then emerges at the other end and propels itself swimming.

Eulerian simulation of complex suspensions and biolocomotion in three dimensions.

We present a numerical method specifically designed for simulating three-dimensional fluid-structure interaction (FSI) problems based on the reference map technique (RMT). The RMT is a fully Eulerian FSI numerical method that allows fluids and large-deformation elastic solids to be represented on a single fixed computational grid. This eliminates the need for meshing complex geometries typical in other FSI approaches and greatly simplifies the coupling between fluid and solids. We develop a three-dimensional implementation of the RMT, parallelized using the distributed memory paradigm, to simulate incompressible FSI with neo-Hookean solids. As part of our method, we develop a field extrapolation scheme that works efficiently in parallel. Through representative examples, we demonstrate the method's suitability in investigating many-body and active systems, as well as its accuracy and convergence. The examples include settling of a mixture of heavy and buoyant soft ellipsoids, lid-driven cavity flow containing a soft sphere, and swimmers actuated via active stress.

An Integrated Experimental-Computational Study of Vocal Fold Vibration in Type I Thyroplasty.

Subject-specific computational modeling of vocal fold (VF) vibration was integrated with an ex vivo animal experiment of type 1 thyroplasty to study the effect of the implant on the vocal fold vibration. In the experiment, a rabbit larynx was used to simulate type 1 thyroplasty, where one side of the vocal fold was medialized with a trans-muscular suture while the other side was medialized with a silastic implant. Vocal fold vibration was then achieved by flowing air through the larynx and was filmed with a high-speed camera. The three-dimensional computational model was built upon the pre-operative scan of the laryngeal anatomy. This subject-specific model was used to simulate the vocal fold medialization and then the fluid-structure interaction (FSI) of the vocal fold. Model validation was done by comparing the vocal fold displacement with postoperative scan (for medialization), and by comparing the vibratory characteristics with the high-speed images (for vibration). These comparisons showed the computational model successfully captured the effect of the implant and thus has the potential for presurgical planning.

Essential role of papillae flexibility in nectar capture by bees.

Many bees possess a tongue resembling a brush composed of a central rod (glossa) covered by elongated papillae, which is dipped periodically into nectar to collect this primary source of energy. In vivo measurements show that the amount of nectar collected per lap remains essentially constant for sugar concentrations lower than 50% but drops significantly for a concentration around 70%. To understand this variation of the ingestion rate with the sugar content of nectar, we investigate the dynamics of fluid capture by Bombus terrestris as a model system. During the dipping process, the papillae, which initially adhere to the glossa, unfold when immersed in the nectar. Combining in vivo investigations, macroscopic experiments with flexible rods, and an elastoviscous theoretical model, we show that the capture mechanism is governed by the relaxation dynamics of the bent papillae, driven by their elastic recoil slowed down through viscous dissipation. At low sugar concentrations, the papillae completely open before the tongue retracts out of nectar and thus, fully contribute to the fluid capture. In contrast, at larger concentrations corresponding to the drop of the ingestion rate, the viscous dissipation strongly hinders the papillae opening, reducing considerably the amount of nectar captured. This study shows the crucial role of flexible papillae, whose aspect ratio determines the optimal nectar concentration, to understand quantitatively the capture of nectar by bees and how physics can shed some light on the degree of adaptation of a specific morphological trait.

Neuromechanical wave resonance in jellyfish swimming.

For organisms to have robust locomotion, their neuromuscular organization must adapt to constantly changing environments. In jellyfish, swimming robustness emerges when marginal pacemakers fire action potentials throughout the bell's motor nerve net, which signals the musculature to contract. The speed of the muscle activation wave is dictated by the passage times of the action potentials. However, passive elastic material properties also influence the emergent kinematics, with time scales independent of neuromuscular organization. In this multimodal study, we examine the interplay between these two time scales during turning. A three-dimensional computational fluid-structure interaction model of a jellyfish was developed to determine the resulting emergent kinematics, using bidirectional muscular activation waves to actuate the bell rim. Activation wave speeds near the material wave speed yielded successful turns, with a 76-fold difference in turning rate between the best and worst performers. Hyperextension of the margin occurred only at activation wave speeds near the material wave speed, suggesting resonance. This hyperextension resulted in a 34-fold asymmetry in the circulation of the vortex ring between the inside and outside of the turn. Experimental recording of the activation speed confirmed that jellyfish actuate within this range, and flow visualization using particle image velocimetry validated the corresponding fluid dynamics of the numerical model. This suggests that neuromechanical wave resonance plays an important role in the robustness of an organism's locomotory system and presents an undiscovered constraint on the evolution of flexible organisms. Understanding these dynamics is essential for developing actuators in soft body robotics and bioengineered pumps.

Fluid-Solid Interaction Analysis for Developing In-Situ Strain and Flow Sensors for Prosthetic Valve Monitoring.

Transcatheter aortic valve implantation (TAVI) was initially developed for adult patients, but there is a growing interest to expand this procedure to younger individuals with longer life expectancies. However, the gradual degradation of biological valve leaflets in transcatheter heart valves (THV) presents significant challenges for this extension. This study aimed to establish a multiphysics computational framework to analyze structural and flow measurements of TAVI and evaluate the integration of optical fiber and photoplethysmography (PPG) sensors for monitoring valve function. A two-way fluid-solid interaction (FSI) analysis was performed on an idealized aortic vessel before and after the virtual deployment of the SAPIEN 3 Ultra (S3) THV. Subsequently, an analytical analysis was conducted to estimate the PPG signal using computational flow predictions and to analyze the effect of different pressure gradients and distances between PPG sensors. Circumferential strain estimates from the embedded optical fiber in the FSI model were highest in the sinus of Valsalva; however, the optimal fiber positioning was found to be distal to the sino-tubular junction to minimize bending effects. The findings also demonstrated that positioning PPG sensors both upstream and downstream of the bioprosthesis can be used to effectively assess the pressure gradient across the valve. We concluded that computational modeling allows sensor design to quantify vessel wall strain and pressure gradients across valve leaflets, with the ultimate goal of developing low-cost monitoring systems for detecting valve deterioration.

Research on the Deformation Law of Foundation Excavation and Support Based on Fluid-Solid Coupling Theory.

To investigate the impact of underground water seepage and soil stress fields on the deformation of excavation and support structures, this study initially identified the key influencing factors on excavation deformation. Subsequently, through a finite element simulation analysis using Plaxis, this study explored the effects of critical factors, such as the excavation support form, groundwater lowering depth, permeability coefficient, excavation layer, and sequence on excavation deformation. Furthermore, a comprehensive consideration of various adverse factors was integrated to establish excavation support early warning thresholds, and optimal dewatering strategies. Finally, this study validated the simulation analysis through an on-site in situ testing with wireless sensors in the context of a physical construction site. The research results indicate that the internal support system within the excavation piles exhibited better stability compared to the external anchor support system, resulting in a 34.5% reduction in the overall deformation. Within the internal support system, the factors influencing the excavation deformation were ranked in the following order: water level (35.5%) > permeability coefficient (17.62%) > excavation layer (11.4%). High water levels, high permeability coefficients, and multi-layered soils were identified as the most unfavorable factors for excavation deformation. The maximum deformation under the coupled effect of these factors was established as the excavation support early warning threshold, and the optimal dewatering strategy involved lowering the water level at the excavation to 0.5 m below the excavation face. The on-site in situ monitoring data obtained through wireless sensors exhibited low discrepancies compared to the finite element simulation data, indicating the high precision of the finite element model for considering the fluid-structure interaction.

A Magneto-Electric Device for Fluid Pipelines with Vibration Damping and Vibration Energy Harvesting.

This study introduces an innovative energy harvesting system designed for industrial applications such as fluid pipelines, air conditioning ducts, sewer systems, and subsea oil pipelines. The system integrates magneto-electric flow coupling and utilizes a dynamic vibration absorber (DVA) to mitigate the vibrations induced by fluid flow while simultaneously harvesting energy through magnetic dipole-dipole interactions in a vibration energy harvester (VEH). The theoretical models, based on Hamilton's Principle and the Biot-Savart Law, were validated through comprehensive experiments. The results indicate the superior performance of the small-magnet system over the large-magnet system in both damping and power generation. The study analyzed the frequency response and energy conversion efficiency across different parameters, including the DVA mass, spring constant, and placement location. The experimental findings demonstrated significant vibration reduction and increased voltage output, validating the theoretical model. This research offers new avenues for energy harvesting systems in pipeline infrastructures, potentially enhancing energy efficiency and structural integrity.

A mechanically consistent unified formulation for fluid-porous-structure-contact interaction.

Fluid-structure interaction with contact poses profound mathematical and numerical challenges, particularly when considering realistic contact scenarios and the influence of surface roughness. Computationally, contact introduces challenges in altering the fluid domain topology and preserving stress balance. This work introduces a new mathematical framework for a unified continuum description of fluid-porous-structure-contact interaction (FPSCI), leveraging the Navier-Stokes-Brinkman (NSB) equations to incorporate porous effects within the surface asperities in the contact region. Our approach maintains mechanical consistency during contact, circumventing issues associated with contact models and complex interface coupling conditions, allowing for the modeling of tangential creeping flows due to surface roughness. The unified continuum and variational multiscale formulation ensure robustness by enabling stable and unified integration of fluid, porous, and solid sub-problems. Computational efficiency and ease of implementation - key advantages of our approach - are demonstrated by solving two benchmark problems of a falling ball and an idealized heart valve. This research has broad implications for fields reliant on accurate fluid-structure interactions and promising advancements in modeling and numerical simulation techniques.

Mesoscopic ring element growth and deformation induced biofilm streamer evolution in microfluidic channels.

In a fluid environment, biofilms usually form and grow into streamers attached to solid surfaces. Existing research on single streamers studied their formation and failure modes. In the experiment on biofilm growth in a microfluidic channel, we found that rings composed of bacteria and an extracellular matrix are important elements on a mesoscopic scale. In the fluid environment, the failure of these ring elements causes damage to streamers. We simulated the growth and deformation of the ring structure in the micro-channel using multi-agent simulation and fluid-structure coupling of a porous elastic body. Based on this, we simulated the biofilm evolution involving multi-ring deformation, which provides a new length scale to study the biofilm streamer dynamics in fluid environments.

Integrated particle image velocimetry and fluid-structure interaction analysis for patient-specific abdominal aortic aneurysm studies.

BACKGROUND: Understanding the hemodynamics of an abdominal aortic aneurysm (AAA) is crucial for risk assessment and treatment planning. This study introduces a low-cost, patient-specific in vitro AAA model to investigate hemodynamics using particle image velocimetry (PIV) and flow-simulating circuit, validated through fluid-structure interaction (FSI) simulations. METHODS: In this study, 3D printing was employed to manufacture a flexible patient-specific AAA phantom using a lost-core casting technique. A pulsatile flow circuit was constructed using off-the-shelf components. A particle image velocimetry (PIV) setup was built using an affordable laser source and global shutter camera, and finally, the flow field inside the AAA was analyzed using open-source software. Fluid-structure interaction (FSI) simulations were performed to enhance our understanding of the flow field, and the results were validated by PIV analysis. Both steady-state and transient flow conditions were investigated. RESULTS: Our experimental setup replicated physiological conditions, analyzing arterial wall deformations and flow characteristics within the aneurysm. Under constant flow, peak wall deformations and flow velocities showed deviations within - 12% to + 27% and - 7% to + 5%, respectively, compared to FSI simulations. Pulsatile flow conditions further demonstrated a strong correlation (Pearson coefficient 0.85) in flow velocities and vectors throughout the cardiac cycle. Transient phenomena, particularly the formation and progression of vortex structures during systole, were consistently depicted between experimental and numerical models. CONCLUSIONS: By bridging high-fidelity experimental observations with comprehensive computational analyses, this study underscores the potential of integrated methodologies in enhancing our understanding of AAA pathophysiology. The convergence of realistic AAA phantoms, precise PIV measurements at affordable cost point, and validated FSI models heralds a new paradigm in vascular research, with significant implications for personalized medicine and bioengineering innovations.

Simulation of the female pelvic mobility and vesical pressure changes employing fluid-structure interaction method.

INTRODUCTION AND HYPOTHESIS: This study aims to develop a fluid-structural interaction (FSI) method to pinpoint the effects of pressure changes inside the bladder and their impact on the supporting structure and the urethra mobility. METHODS: A physiological model of the nulliparous female pelvis, including the organs, supportive structures, and urine, was developed based on magnetic resonance images. Soft tissues with nonlinear hyperelastic material characteristics were modeled. The Navier-Stokes equations governing the fluid flow within the computational domain (urine) were solved. The urine and soft tissue interactions were simulated by the FSI method. The vesical pressure and its impact on the urethral mobility and supportive structures were investigated during the Valsalva maneuver. Moreover, the simulation results were validated by comparing with a urodynamic test and other research. RESULTS: The results demonstrated that the vesical pressure simulated by the FSI method could predict the nonlinear behavior of the urodynamic test pressure. The urethra retropubic bladder neck and the bladder neck-pubic bone angle changed 58.92% and -55.76%, respectively. The retropubic urethral length distance changed by -48.74%. The error compared to the statistical results of other research is < 5%. CONCLUSIONS: The total deformation and mobility of the urethra predicted by the FSI model were consistent with clinical observations in a subject. The urethra supports dependence on the tissues' mechanical properties, interaction between the tissues, and effect of urine fluid inside the bladder. This simulation effectively depicts the patterns of urethra mobility, which provides a better understanding of the behavior of the pelvic floor.

The reliability and utility of on-site CT-derived fractional flow reserve (FFR) based on fluid structure interactions: comparison with FFRCT based on computational fluid dynamics, invasive FFR, and resting full-cycle ratio.

Fractional flow reserve (FFR) derived off-site by coronary computed tomography angiography (CCTA) (FFRCT) is obtained by applying the principles of computational fluid dynamics. This study aimed to validate the overall reliability of on-site CCTA-derived FFR based on fluid structure interactions (CT-FFR) and assess its clinical utility compared with FFRCT, invasive FFR, and resting full-cycle ratio (RFR). We calculated the CT-FFR for 924 coronary vessels in 308 patients who underwent CCTA for clinically suspected coronary artery disease. Of these patients, 35 patients with at least one obstructive stenosis (> 50%) detected on CCTA underwent both CT-FFR and FFRCT for further investigation. Furthermore, 24 and 20 patients underwent invasive FFR and RFR in addition to CT-FFR, respectively. The inter-observer correlation (r) of CT-FFR was 0.93 (95% confidence interval [CI] 0.85-0.97, P < 0.0001) with a mean absolute difference of - 0.0042 (limits of agreement - 0.073, 0.064); 97.3% of coronary arteries without obstructive lesions on CCTA had negative results for ischemia on CT-FFR (> 0.80). The correlation coefficient between CT-FFR and FFRCT for 105 coronary vessels was 0.87 (95% CI 0.82-0.91, P < 0.0001) with a mean absolute difference of - 0.012 (limits of agreement - 0.12, 0.10). CT-FFR correlated well with both invasive FFR (r = 0.66, 95% CI 0.36-0.84, P = 0.0003) and RFR (r = 0.78, 95% CI 0.51-0.91, P < 0.0001). These data suggest that CT-FFR can potentially substitute for FFRCT and correlates closely with invasive FFR and RFR with high reproducibility. Our findings should be proven by further clinical investigation in a larger cohort.

Framework Development for Patient-Specific Compliant Aortic Dissection Phantom Model Fabrication: Magnetic Resonance Imaging Validation and Deep-Learning Segmentation.

Type B aortic dissection is a life-threatening medical emergency that can result in rupture of the aorta. Due to the complexity of patient-specific characteristics, only limited information on flow patterns in dissected aortas has been reported in the literature. Leveraging the medical imaging data for patient-specific in vitro modeling can complement the hemodynamic understanding of aortic dissections. We propose a new approach toward fully automated patient-specific type B aortic dissection model fabrication. Our framework uses a novel deep-learning-based segmentation for negative mold manufacturing. Deep-learning architectures were trained on a dataset of 15 unique computed tomography scans of dissection subjects and were blind-tested on 4 sets of scans, which were targeted for fabrication. Following segmentation, the three-dimensional models were created and printed using polyvinyl alcohol. These models were then coated with latex to create compliant patient-specific phantom models. The magnetic resonance imaging (MRI) structural images demonstrate the ability of the introduced manufacturing technique for creating intimal septum walls and tears based on patient-specific anatomy. The in vitro experiments show the fabricated phantoms generate physiologically-accurate pressure results. The deep-learning models also show high similarity metrics between manual segmentation and autosegmentation where Dice metric is as high as 0.86. The proposed deep-learning-based negative mold manufacturing method facilitates an inexpensive, reproducible, and physiologically-accurate patient-specific phantom model fabrication suitable for aortic dissection flow modeling.

Ultra-sharp pinnacles sculpted by natural convective dissolution.

The evolution of landscapes, landforms, and other natural structures involves highly interactive physical and chemical processes that often lead to intriguing shapes and recurring motifs. Particularly intricate and fine-scale features characterize the so-called karst morphologies formed by mineral dissolution into water. An archetypal form is the tall, slender, and sharply tipped karst pinnacle or rock spire that appears in multitudes in striking landforms called stone forests, but whose formative mechanisms remain unclear due to complex, fluctuating, and incompletely understood developmental conditions. Here, we demonstrate that exceedingly sharp spires also form under the far-simpler conditions of a solid dissolving into a surrounding liquid. Laboratory experiments on solidified sugars in water show that needlelike pinnacles, as well as bed-of-nails-like arrays of pinnacles, emerge robustly from the dissolution of solids with smooth initial shapes. Although the liquid is initially quiescent and no external flow is imposed, persistent flows are generated along the solid boundary as dense, solute-laden fluid descends under gravity. We use these observations to motivate a mathematical model that links such boundary-layer flows to the shape evolution of the solid. Dissolution induces these natural convective flows that, in turn, enhance dissolution rates, and simulations show that this feedback drives the shape toward a finite-time singularity or blow-up of apex curvature that is cut off once the pinnacle tip reaches microscales. This autogenic mechanism produces ultra-fine structures as an attracting state or natural consequence of the coupled processes at work in the closed solid-fluid system.

Thinner biological tissues induce leaflet flutter in aortic heart valve replacements.

Valvular heart disease has recently become an increasing public health concern due to the high prevalence of valve degeneration in aging populations. For patients with severely impacted aortic valves that require replacement, catheter-based bioprosthetic valve deployment offers a minimally invasive treatment option that eliminates many of the risks associated with surgical valve replacement. Although recent percutaneous device advancements have incorporated thinner, more flexible biological tissues to streamline safer deployment through catheters, the impact of such tissues in the complex, mechanically demanding, and highly dynamic valvular system remains poorly understood. The present work utilized a validated computational fluid-structure interaction approach to isolate the behavior of thinner, more compliant aortic valve tissues in a physiologically realistic system. This computational study identified and quantified significant leaflet flutter induced by the use of thinner tissues that initiated blood flow disturbances and oscillatory leaflet strains. The aortic flow and valvular dynamics associated with these thinner valvular tissues have not been previously identified and provide essential information that can significantly advance fundamental knowledge about the cardiac system and support future medical device innovation. Considering the risks associated with such observed flutter phenomena, including blood damage and accelerated leaflet deterioration, this study demonstrates the potentially serious impact of introducing thinner, more flexible tissues into the cardiac system.

A computer model of the effect of venous-arterial-venous extracorporeal membrane oxygenation venous return of oxygenated blood on right atrial recirculation and tricuspid valve blood oxygen saturation.

OBJECTIVE: Venous-arterial venous extracorporeal membrane oxygenation (V-AV ECMO), as a new clinical application of ECMO, showed great clinical application potential in the treatment of patients with combined cardiopulmonary failure. Given the complicated cannulation strategy of V-AV ECMO, its influence on the hemodynamics of the human circulatory system remained unclear. METHODS: In this paper, a fluid-structure interaction was used to study the effect of V-AV ECMO oxygenated blood shunt ratio on right atrial recirculation and tricuspid valve (TV) blood oxygen saturation. In this study, the right atrium, superior vena cava supplying cannulae and inferior vena cava draining cannulae model of a specific patient was constructed. Seven cases with shunt ratio of 12.50%, 18.75%, 25.00%, 31.25%, 37.50%, 43.75% and 50.00% were designed. RESULTS: The streamline diagram and velocity contour of oxygenated blood, recirculation fraction (RF), correlation of three variables (shunt ratio, RF, and oxygen saturation), and the oxygen saturation of blood at the TV were extracted for the study. Study results showed that, first, as the shunt ratio increased, the RF of the seven cases was 14.64%, 29.87%, 33.85%, 40.12%, 40.40%, 40.02%, and 38.09%. Second, with the increase of the shunt ratio, oxygen saturation of blood at the TV in seven cases was 82.1%, 82.5%, 83.3%, 83.3%, 84.0%, 84.6%, and 85.3%. CONCLUSIONS: In this study, the shunt ratio had a strong correlation with the RF and oxygen saturation of blood at the TV. As the shunt ratio increased, the RF initially increased and then stabilized. However, oxygen saturation of blood at the TV would increase with the increase of the shunt ratio, but the degree of increase was small. This research provided useful information for surgeons and operators using V-AV ECMO.