Aerodynamic Simulation of Small Airway Resistance: A New Imaging Biomarker for Chronic Obstructive Pulmonary Disease.

Purpose: To develop a novel method for calculating small airway resistance using computational fluid dynamics (CFD) based on CT data and evaluate its value to identify COPD. Patients and Methods: 24 subjects who underwent chest CT scans and pulmonary function tests between August 2020 and December 2020 were enrolled retrospectively. Subjects were divided into three groups: normal (10), high-risk (6), and COPD (8). The airway from the trachea down to the sixth generation of bronchioles was reconstructed by a 3D slicer. The small airway resistance (RSA) and RSA as a percentage of total airway resistance (RSA%) were calculated by CFD combined with airway resistance and FEV1 measured by pulmonary function test. A correlation analysis was conducted between RSA and pulmonary function parameters, including FEV1/FVC, FEV1% predicted, MEF50% predicted, MEF75% predicted and MMEF75/25% predicted. Results: The RSA and RSA% were significantly different among the three groups (p<0.05) and related to FEV1/FVC (r = -0.70, p < 0.001; r = -0.67, p < 0.001), FEV1% predicted (r = -0.60, p = 0.002; r = -0.57, p = 0.004), MEF50% predicted (r = -0.64, p = 0.001; r = -0.64, p = 0.001), MEF75% predicted (r = -0.71, p < 0.001; r = -0.60, p = 0.002) and MMEF 75/25% predicted (r = -0.64, p = 0.001; r = -0.64, p = 0.001). Conclusion: Airway CFD is a valuable method for estimating the small airway resistance, where the derived RSA will aid in the early diagnosis of COPD.

Urban flood modeling with a novel coupling method of surface and sewer hydrodynamic processes.

Drainage modeling that accurately captures urban storm inundation serves as the foundation for flood warning and drainage scheduling. In this paper, we proposed a novel coupling ideology that, by integrating 2D-1D and 1D-2D unidirectional processes, overcomes the drawback of the conventional unidirectional coupling approach that fails to properly represent the rainfall surface catchment dynamics, and provides more coherent hydrological implications compared to the bidirectional coupling concept. This paper first referred to a laboratory experimental case from the literature, applied and analyzed the coupling scheme proposed in this paper and the bidirectional coupling scheme that has been widely studied in recent years, compared the two coupling solutions in terms of the resulting accuracy and applicability, and discussed their respective strengths and weaknesses to validate the reliability of the proposed method. The verified proposed coupling scheme was then applied to the modeling of a real drainage system in a region of Nanjing, China, and the results proved that the coupling mechanism proposed in this study is of practical application value.

Hydrodynamic Assessment of Explanted Degenerated Transcatheter Aortic Valves: Novel Insights Into Noncalcific and Calcific Mechanisms.

BACKGROUND: The etiology of transcatheter aortic valve (TAV) degeneration is poorly understood, particularly noncalcific mechanisms. OBJECTIVES: The authors sought to investigate noncalcific and calcific mechanisms of TAV degeneration and evaluate their impact on leaflet function by bench testing, imaging, and histology. METHODS: TAV explants were obtained from the EXPLANT THV registry and clinical institutions. Hydrodynamic assessment was performed using a heart valve pulse duplicator system under physiological conditions. Micro-computed tomography, high-resolution photography, high speed video, and hematoxylin and eosin staining were used to evaluate the morphological appearance, leaflet kinematics, and calcium burden of TAVs. RESULTS: A total of 14 explants were evaluated: 10 self-expanding CoreValve/Evolut TAVs (Medtronic), 3 balloon-expandable SAPIEN 3 TAVs (Edwards Lifesciences), and 1 mechanically expandable Lotus TAV (Boston Scientific). The median patient age at explantation was 73.0 years (Q1-Q3: 64.5-80.0 years), with a time to explantation of 4 years 1 month (1 year 5 months to 4 years 11 months). Six TAV explants were found to have leaflet calcification (162.4 mm3; 58.8-603.0 mm3), and 8 had no calcification detectable by micro-computed tomography and histology. All samples had impaired leaflet kinematics. There was no significant difference in the hydrodynamic mean gradient between calcified (47.2 mm Hg; 26.6-74.1 mm Hg) and noncalcified (27.6 mm Hg; 15.2-36.7 mm Hg; P = 0.28) TAVs. Leaflet calcification had a weak but nonsignificant association with the hydrodynamic mean gradient (r = 0.42; P = 0.14). CONCLUSIONS: TAV function can be severely impacted by noncalcific and calcific mechanisms of tissue degeneration. Importantly, functional stenosis can occur in TAVs in the absence of obvious and significant calcification.

Enhanced microfluidic multi-target separation by positive and negative magnetophoresis.

We introduce magnetophoresis-based microfluidics for sorting biological targets using positive Magnetophoresis (pM) for magnetically labeled particles and negative Magnetophoresis (nM) for label-free particles. A single, externally magnetized ferromagnetic wire induces repulsive forces and is positioned across the focused sample flow near the main channel's closed end. We analyze magnetic attributes and separation performance under two transverse dual-mode magnetic configurations, examining magnetic fields, hydrodynamics, and forces on microparticles of varying sizes and properties. In pM, the dual-magnet arrangement (DMA) for sorting three distinct particles shows higher magnetic gradient generation and throughput than the single-magnet arrangement (SMA). In nM, the numerical results for SMA sorting of red blood cells (RBCs), white blood cells (WBCs), and prostate cancer cells (PC3-9) demonstrate superior magnetic properties and throughput compared to DMA. Magnetized wire linear movement is a key design parameter, allowing device customization. An automated device for handling more targets can be created by manipulating magnetophoretic repulsion forces. The transverse wire and magnet arrangement accommodate increased channel depth without sacrificing efficiency, yielding higher throughput than other devices. Experimental validation using soft lithography and 3D printing confirms successful sorting and separation, aligning well with numerical results. This demonstrates the successful sorting and separating of injected particles within a hydrodynamically focused sample in all systems. Both numerical and experimental findings indicate a separation accuracy of 100% across various Reynolds numbers. The primary channel dimensions measure 100 µm in height and 200 µm in width. N52 permanent magnets were employed in both numerical simulations and experiments. For numerical simulations, a remanent flux density of 1.48 T was utilized. In the experimental setup, magnets measuring 0.5 × 0.5 × 0.125 inches and 0.5 × 0.5 × 1 inch were employed. The experimental data confirm the device's capability to achieve 100% separation accuracy at a Reynolds number of 3. However, this study did not explore the potential impact of increased flow rates on separation accuracy.

Analysis of the urine flow characteristics inside catheters for intermittent catheter selection.

In this study, we conducted a numerical analysis on catheter sizes using computational fluid dynamics to assess urinary flow rates during intermittent catheterization (IC). The results revealed that the fluid (urine) movement within a catheter is driven by intravesical pressure, with friction against the catheter walls being the main hindrance to fluid movement. Higher-viscosity fluids experienced increased friction with increasing intravesical pressure, resulting in reduced fluid velocity, whereas lower-viscosity fluids experienced reduced friction under similar pressure, leading to increased fluid velocity. Regarding urine characteristics, the results indicated that bacteriuria, with lower viscosity, exhibited higher flow rates, whereas glucosuria exhibited the lowest flow rates. Additionally, velocity gradients decreased with increasing catheter diameters, reducing friction and enhancing fluid speed, while the friction increased with decreasing diameters, reducing fluid velocity. These findings confirm that flow rates increased with larger catheter sizes. Furthermore, in terms of specific gravity, the results showed that a 12Fr catheter did not meet the ISO-suggested average flow rate (50 cc/min). The significance of this study lies in its application of fluid dynamics to nursing, examining urinary flow characteristics in catheterization. It is expected to aid nurses in selecting appropriate catheters for intermittent catheterization based on urinary test results.

In silico intravascular optical coherence tomography (IVOCT) for quality assured imaging with reduced intervention.

In the clinical application of intravascular optical coherence tomography (IVOCT), it is necessary to flush opaque blood during image acquisition. However, there are no specific standards for how to perform low-dose but effective flushing. In this study, computational fluid dynamics (CFD) and optical models were integrated to numerically simulate the complete process of IVOCT, which includes blood flushing with normal saline followed by image acquisition. Moreover, an intermittent injection scheme was proposed, and its advantages over the conventionally adopted scheme of continuous injection were verified. The results show that intermittent injection can significantly reduce the dosage of normal saline (reduced by 44.4%) with only a slight sacrifice of image quality (reduced by 8.7%, but still acceptable). The developed model and key findings in this work can help surgeons practice optimized IVOCT operations and potentially lead to improved designs of the IVOCT equipment.

Aqueductal CSF stroke volume is associated with the burden of perivascular space enlargement in chronic adult hydrocephalus.

The inflow of CSF into perivascular spaces (PVS) in the brain is crucial for clearing waste molecules. Inefficiency in PVS flow leads to neurodegeneration. Failure of PVS flushing is associated with CSF flow impairment in the intracranial hydrodynamic condition of CSF hypo-pulsatility. However, enlarged PVS (ePVS), a finding indicative of PVS flow dysfunction, is also present in patients with derangement of CSF dynamics characterized by CSF hyper-pulsatility, which increases CSF flow. Intriguingly, two opposite intracranial hydrodynamic conditions would lead to the same result of impairing the PVS flushing. To investigate this issue, we assessed the subsistence of a dysfunctional interplay between CSF and PVS flows and, if the case, the mechanisms preventing a hyper-pulsatile brain from providing an effective PVS flushing. We analyzed the association between phase contrast MRI aqueductal CSF stroke volume (aqSV), a proxy of CSF pulsatility, and the burden of ePVS in chronic adult hydrocephalus, a disease involving a broad spectrum of intracranial hydrodynamics disturbances. In the 147 (85 males, 62 females) patients, the age at diagnosis ranged between 28 and 88 years (median 73 years). Ninety-seven patients had tri-ventriculomegaly and 50 tetra-ventriculomegaly. According to the extent of ePVS, 113 patients had a high ePVS burden, while 34 had a low ePVS burden. aqSV, which ranged between 0 and 562 µL (median 86 µL), was increased with respect to healthy subjects. Patients presenting with less ePVS burden had higher aqSV (p < 0.002, corrected for the multiple comparisons) than those with higher ePVS burden. The present study confirmed the association between CSF dynamics and PVS flow disturbances and demonstrated this association in intracranial hyper-pulsatility. Further studies should investigate the association between PVS flow failure and CSF hypo- and hyper-pulsatility as responsible/co-responsible for glymphatic failure in other neurodegenerative diseases, particularly in diseases in which CSF disturbances can be corrected, as in chronic adult hydrocephalus.

Collective movement of schooling fish reduces the costs of locomotion in turbulent conditions.

The ecological and evolutionary benefits of energy-saving in collective behaviors are rooted in the physical principles and physiological mechanisms underpinning animal locomotion. We propose a turbulence sheltering hypothesis that collective movements of fish schools in turbulent flow can reduce the total energetic cost of locomotion by shielding individuals from the perturbation of chaotic turbulent eddies. We test this hypothesis by quantifying energetics and kinematics in schools of giant danio (Devario aequipinnatus) and compared that to solitary individuals swimming under laminar and turbulent conditions over a wide speed range. We discovered that, when swimming at high speeds and high turbulence levels, fish schools reduced their total energy expenditure (TEE, both aerobic and anaerobic energy) by 63% to 79% compared to solitary fish (e.g., 228 versus 48 kj kg-1). Solitary individuals spend approximately 22% more kinematic effort (tail beat amplitude•frequency: 1.7 versus 1.4 BL s-1) to swim in turbulence at higher speeds than in laminar conditions. Fish schools swimming in turbulence reduced their three-dimensional group volume by 41% to 68% (at higher speeds, approximately 103 versus 33 cm3) and did not alter their kinematic effort compared to laminar conditions. This substantial energy saving highlights that schooling behaviors can mitigate turbulent disturbances by sheltering fish (within schools) from the eddies of sufficient kinetic energy that can disrupt locomotor gaits. Therefore, providing a more desirable internal hydrodynamic environment could be one of the ecological drivers underlying collective behaviors in a dense fluid environment.

Adaptive design of experiments to fit surrogate Gaussian process regression models allows fast sensitivity analysis of the input waveform for patient-specific 3D CFD models of liver radioembolization.

BACKGROUND AND OBJECTIVE: Patient-specific 3D computational fluid dynamics (CFD) models are increasingly being used to understand and predict transarterial radioembolization procedures used for hepatocellular carcinoma treatment. While sensitivity analyses of these CFD models can help to determine the most impactful input parameters, such analyses are computationally costly. Therefore, we aim to use surrogate modelling to allow relatively cheap sensitivity analysis. As an example, we compute Sobol's sensitivity indices for three input waveform shape parameters. METHODS: We extracted three characteristic shape parameters from our input mass flow rate waveform (peak systolic mass flow rate, heart rate, systolic duration) and defined our 3D input parameter space by varying these parameters within 75 %-125 % of their nominal values. To fit our surrogate model with a minimal number of costly CFD simulations, we developed an adaptive design of experiments (ADOE) algorithm. The ADOE uses 100 Latin hypercube sampled points in 3D input space to define the initial design of experiments (DOE). Subsequently, we re-sample input space with 10,000 Latin Hypercube sampled points and cheaply estimate the outputs using the surrogate model. In each of 27 equivolume bins which divide our input space, we determine the most uncertain prediction of the 10,000 points, compute the true outputs using CFD, and add these points to the DOE. For each ADOE iteration, we calculate Sobol's sensitivity indices, and we continue to add batches of 27 samples to the DOE until the Sobol indices have stabilized. RESULTS: We tested our ADOE algorithm on the Ishigami function and showed that we can reliably obtain Sobol's indices with an absolute error <0.1. Applying ADOE to our waveform sensitivity problem, we found that the first-order sensitivity indices were 0.0550, 0.0191 and 0.407 for the peak systolic mass flow rate, heart rate, and the systolic duration, respectively. CONCLUSIONS: Although the current study was an illustrative case, the ADOE allows reliable sensitivity analysis with a limited number of complex model evaluations, and performs well even when the optimal DOE size is a priori unknown. This enables us to identify the highest-impact input parameters of our model, and other novel, costly models in the future.

A computational fluid dynamics study to assess the impact of coughing on cerebrospinal fluid dynamics in Chiari type 1 malformation.

Chiari type 1 malformation is a neurological disorder characterized by an obstruction of the cerebrospinal fluid (CSF) circulation between the brain (intracranial) and spinal cord (spinal) compartments. Actions such as coughing might evoke spinal cord complications in patients with Chiari type 1 malformation, but the underlying mechanisms are not well understood. More insight into the impact of the obstruction on local and overall CSF dynamics can help reveal these mechanisms. Therefore, our previously developed computational fluid dynamics framework was used to establish a subject-specific model of the intracranial and upper spinal CSF space of a healthy control. In this model, we emulated a single cough and introduced porous zones to model a posterior (OBS-1), mild (OBS-2), and severe posterior-anterior (OBS-3) obstruction. OBS-1 and OBS-2 induced minor changes to the overall CSF pressures, while OBS-3 caused significantly larger changes with a decoupling between the intracranial and spinal compartment. Coughing led to a peak in overall CSF pressure. During this peak, pressure differences between the lateral ventricles and the spinal compartment were locally amplified for all degrees of obstruction. These results emphasize the effects of coughing and indicate that severe levels of obstruction lead to distinct changes in intracranial pressure.

Adapting to the Abyss: Passive Ventilation in the Deep-Sea Glass Sponge Euplectella aspergillum.

We analyze the flow physics inside the body cavity and downstream the deep-sea glass sponge Euplectella aspergillum. We provide evidence that the helical skeletal motifs of the sponge give rise to a rich fluid dynamic field, allowing the organism to scavenge flow from the bottom of the sea and promoting a spontaneous, organized vertical flow within its body cavity toward the osculum. Our analysis points at a functional adaptation of the organism, which can passively divert flow through the osculum in unfavorable, low ambient currents, with no need for active pumping, with potential repercussions in functional ecology, as well as the design of chemical reactors, air-treatment units, and civil and aeronaval structures.

Multiflagellate Swimming Controlled by Hydrodynamic Interactions.

Many eukaryotic microorganisms propelled by multiple flagella can swim very rapidly with distinct gaits. Here, we model a three-dimensional mutiflagellate swimmer, resembling the microalgae. When the flagella are actuated synchronously, the swimming efficiency can be enhanced or reduced by interflagella hydrodynamic interactions (HIs), determined by the intrinsic tilting angle of the flagella. The asynchronous gait with a phase difference between neighboring flagella can reduce oscillatory motion via the basal mechanical coupling. In the presence of a spherical body, simulations taking into account the flagella-body interactions reveal the advantage of anterior configuration compared with posterior configuration, where in the latter case an optimal flagella number arises. Apart from understanding the role of HIs in the multiflagellate microorganisms, this work could also guide laboratory fabrications of novel microswimmers.

Influence of bacterial swimming and hydrodynamics on attachment of phages.

Bacteriophages ("phages") are viruses that infect bacteria. Since they do not actively self-propel, phages rely on thermal diffusion to find target cells-but can also be advected by fluid flows, such as those generated by motile bacteria themselves in bulk fluids. How does the flow field generated by a swimming bacterium influence how it encounters phages? Here, we address this question using coupled molecular dynamics and lattice Boltzmann simulations of flagellated bacteria swimming through a bulk fluid containing uniformly-dispersed phages. We find that while swimming increases the rate at which phages attach to both the cell body and flagellar propeller, hydrodynamic interactions strongly suppress this increase at the cell body, but conversely enhance this increase at the flagellar bundle. Our results highlight the pivotal influence of hydrodynamics on the interactions between bacteria and phages, as well as other diffusible species, in microbial environments.

Are man-made hives valid thermal surrogates for natural honey bee nests (Apismellifera)?

Honey bees preferentially occupy thick walled tall narrow tree cavities and attach their combs directly to the nest wall, leaving periodic gaps. However, academic research and beekeeping are conducted in squat, thin walled man made hives, with a continuous gap between the combs and the walls and roof. Utilising a computational fluid dynamics (CFD) model of thermoregulating bees in complete nests in trees and thin walled man made hives, with the average size of tree comb gaps determined from honey bee occupied synthetic tree nests, this research compared the metabolic energy impacts of comb gaps and vertical movement of the thermoregulated brood area. This shows their heat transfer regimes are disparate, including: bee space above combs increases heat loss by up to ∼70%; hives, compared to tree nests, require at least 150% the density of honey bees to arrest convection across the brood area. Tree cavities have a larger vertical freedom, a greater thermal resistance and can make dense clustering redundant. With the thermal environment being critical to honey bees, the magnitude and scope of these differences suggest that some hive based behavioural research needs extra validation to be considered non-anthropogenic, and some bee keeping practices are sub-optimal.

Flow sensing on dragonfly wings.

One feature of animal wings is their embedded mechanosensory system that can support flight control. Insect wings are particularly interesting as they are highly deformable yet the actuation is limited to the wing base. It is established that strain sensors on insect wings can directly mediate reflexive control; however, little is known about airflow sensing by insect wings. What information can flow sensors capture and how can flow sensing benefit flight control? Here, we use the dragonfly (Sympetrum striolatum) as a model to explore the function of wing sensory bristles in the context of flight control. Combining our detailed anatomical reconstructions of both the sensor microstructures and wing architecture, we used computational fluid dynamics simulations to ask the following questions. (1) Are there strategic locations on wings that sample flow for estimating aerodynamically relevant parameters such as the local effective angle of attack? (2) Is the sensory bristle distribution on dragonfly wings optimal for flow sensing? (3) What is the aerodynamic effect of microstructures found near the sensory bristles on dragonfly wings? We discuss the benefits of flow sensing for flexible wings and how the evolved sensor placement affects information encoding.

The interplay of group size and flow velocity modulates fish exploratory behaviour.

Social facilitation is a well-known phenomenon where the presence of organisms belonging to the same species enhances an individual organism's performance in a specific task. As far as fishes are concerned, most studies on social facilitation have been conducted in standing-water conditions. However, for riverine species, fish are most commonly located in moving waters, and the effects of hydrodynamics on social facilitation remain largely unknown. To bridge this knowledge gap, we designed and performed flume experiments where the behaviour of wild juvenile Italian riffle dace (Telestes muticellus) in varying group sizes and at different mean flow velocities, was studied. An artificial intelligence (AI) deep learning algorithm was developed and employed to track fish positions in time and subsequently assess their exploration, swimming activity, and space use. Results indicate that energy-saving strategies dictated space use in flowing waters regardless of group size. Instead, exploration and swimming activity increased by increasing group size, but the magnitude of this enhancement (which quantifies social facilitation) was modulated by flow velocity. These results have implications for how future research efforts should be designed to understand the social dynamics of riverine fish populations, which can no longer ignore the contribution of hydrodynamics.

Advancements of the Fluidized Bed Fenton (FBF) Technology for wastewater treatment: Mechanism, mass and heat transfer.

Fluidized Bed Fenton (FBF) technology, a fusion of the Fenton method and fluidized bed reactor, has emerged as a superior alternative to conventional Fenton technology for treating organic industrial wastewater. This innovative approach has garnered significant attention from researchers in recent years. While earlier studies primarily focused on pollutant degradation in simulated wastewater and catalyst development, there has been a growing interest in examining the alterations in mass or heat transfer performance attributed to fluidized beds. This paper explores the factors that contribute to the effectiveness of Fluidized Bed Fenton technology in efficiently degrading various challenging organic pollutants, while also reducing iron sludge production and expanding the applicable pH range, through an analysis of reaction kinetics. Meanwhile, combined with the related work of fluid dynamics, the research related to mass and heat transfer inside the reactor of Fluidized Bed Fenton technology is summarized, and it is proposed that the use of computers to establish a suitable model of Fluidized Bed Fenton and solve it with the assistance of computational fluid dynamics (CFD) and other software will help to further explore the process of mass and heat transfer inside the fluidized bed, which will provide the basis for the future of the Fluidized Bed Fenton from the laboratory to the actual industrial application.

Computational study of effect of hybrid nanoparticles on hemodynamics and thermal transfer in ruptured arteries with pathological dilation.

The intended research aims to explore the convection phenomena of a hybrid nanofluid composed of gold and silver nanoparticles. This research is novel and significant because there is a lack of existing studies on the flow behavior of hybrid nanoparticles with important physical properties of blood base fluids, especially in the case of sidewall ruptured dilated arteries. The implementation of combined nanoparticles rather than unadulterated nanoparticles is one of the most crucial elements in boosting the thermal conduction of fluids. The research methodology encompasses the utilization of advanced bio-fluid dynamics software for simulating the flow of the nanofluid. The physical context elucidates the governing equations of momentum, mass, momentum, and energy in terms of partial differential equations. The results are displayed in both tabular and graphical forms to demonstrate the numerical and graphical solutions. The effect of physical parameters on velocity distribution is illustrated through graphs. Furthermore, the study's findings are unique and original, and these computational discoveries have not been published by any researcher before. The finding implies that utilizing hybrid nanoparticles as drug carriers holds great promise in mitigating the effects of blood flow, potentially enhancing drug delivery, and minimizing its impact on the body.

Effect of an accelerating metro cabin on the diffusion of cough droplets.

Coronaviruses being capable of spreading through droplet contamination have raised significant concerns regarding high-capacity public rail transport, such as the metro. Within a rapidly moving railcar cabin, the internal airflow lags behind the bulkhead, generating internally induced airflow that accelerates droplet dispersion within a non-inertial reference system. This study investigates the impact of acceleration on the diffusion of cough droplets of varying sizes using computational fluid dynamics. The modified k-ε equation in ANSYS® Fluent was utilized to simulate droplet diffusion under different body orientations by adjusting the inertial force correction source term. Results indicate that droplets in the middle size range (50-175 µm) are primarily influenced by inertial forces, whereas smaller droplets (3.5-20 µm) are predominantly controlled by air drag forces. Regardless of facial orientation, the outlet of high-capacity public rail transport poses the highest risk of infection.

Modeling cerebrospinal fluid dynamics across the entire intracranial space through integration of four-dimensional flow and intravoxel incoherent motion magnetic resonance imaging.

BACKGROUND: Bidirectional reciprocal motion of cerebrospinal fluid (CSF) was quantified using four-dimensional (4D) flow magnetic resonance imaging (MRI) and intravoxel incoherent motion (IVIM) MRI. To estimate various CSF motions in the entire intracranial region, we attempted to integrate the flow parameters calculated using the two MRI sequences. To elucidate how CSF dynamics deteriorate in Hakim's disease, an age-dependent chronic hydrocephalus, flow parameters were estimated from the two MRI sequences to assess CSF motion in the entire intracranial region. METHODS: This study included 127 healthy volunteers aged ≥ 20 years and 44 patients with Hakim's disease. On 4D flow MRI for measuring CSF motion, velocity encoding was set at 5 cm/s. For the IVIM MRI analysis, the diffusion-weighted sequence was set at six b-values (i.e., 0, 50, 100, 250, 500, and 1000 s/mm2), and the biexponential IVIM fitting method was adapted. The relationships between the fraction of incoherent perfusion (f) on IVIM MRI and 4D flow MRI parameters including velocity amplitude (VA), absolute maximum velocity, stroke volume, net flow volume, and reverse flow rate were comprehensively evaluated in seven locations in the ventricles and subarachnoid spaces. Furthermore, we developed a new parameter for fluid oscillation, the Fluid Oscillation Index (FOI), by integrating these two measurements. In addition, we investigated the relationship between the measurements and indices specific to Hakim's disease and the FOIs in the entire intracranial space. RESULTS: The VA on 4D flow MRI was significantly associated with the mean f-values on IVIM MRI. Therefore, we estimated VA that could not be directly measured on 4D flow MRI from the mean f-values on IVIM MRI in the intracranial CSF space, using the following formula; e0.2(f-85) + 0.25. To quantify fluid oscillation using one integrated parameter with weighting, FOI was calculated as VA × 10 + f × 0.02. In addition, the FOIs at the left foramen of Luschka had the strongest correlations with the Evans index (Pearson's correlation coefficient: 0.78). The other indices related with Hakim's disease were significantly associated with the FOIs at the cerebral aqueduct and bilateral foramina of Luschka. FOI at the cerebral aqueduct was also elevated in healthy controls aged ≥ 60 years. CONCLUSIONS: We estimated pulsatile CSF movements in the entire intracranial CSF space in healthy individuals and patients with Hakim's disease using FOI integrating VA from 4D flow MRI and f-values from IVIM MRI. FOI is useful for quantifying the CSF oscillation.