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.

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.

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.

Modeling and experimental study of the intervention forces between the guidewire and blood vessels.

A profound investigation of the interaction mechanics between blood vessels and guidewires is necessary to achieve safe intervention. An interactive force model between guidewires and blood vessels is established based on cardiovascular fluid dynamics theory and contact mechanics, considering two intervention phases (straight intervention and contact intervention at a corner named "J-vessel"). The contributing factors of the force model, including intervention conditions, guidewire characteristics, and intravascular environment, are analyzed. A series of experiments were performed to validate the availability of the interactive force model and explore the effects of influential factors on intervention force. The intervention force data were collected using a 2-DOF mechanical testing system instrumented with a force sensor. The guidewire diameter and material were found to significantly impact the intervention force. Additionally, the intervention force was influenced by factors such as blood viscosity, blood vessel wall thickness, blood flow velocity, as well as the interventional velocity and interventional mode. The experiment of the intervention in a coronary artery physical vascular model confirms the practicality validation of the predicted force model and can provide an optimized interventional strategy for vascular interventional surgery. The enhanced intervention strategy has resulted in a considerable reduction of approximately 21.97 % in the force exerted on blood vessels, effectively minimizing the potential for complications associated with the interventional surgery.

Ediacaran marine animal forests and the ventilation of the oceans.

The rise of animals across the Ediacaran-Cambrian transition marked a step-change in the history of life, from a microbially dominated world to the complex macroscopic biosphere we see today.1,2,3 While the importance of bioturbation and swimming in altering the structure and function of Earth systems is well established,4,5,6 the influence of epifaunal animals on the hydrodynamics of marine environments is not well understood. Of particular interest are the oldest "marine animal forests,"7 which comprise a diversity of sessile soft-bodied organisms dominated by the fractally branching rangeomorphs.8,9 Typified by fossil assemblages from the Ediacaran of Mistaken Point, Newfoundland,8,10,11 these ancient communities might have played a pivotal role in structuring marine environments, similar to modern ecosystems,7,12,13 but our understanding of how they impacted fluid flow in the water column is limited. Here, we use ecological modeling and computational flow simulations to explore how Ediacaran marine animal forests influenced their surrounding environment. Our results reveal how organism morphology and community structure and composition combined to impact vertical mixing of the surrounding water. We find that Mistaken Point communities were capable of generating high-mixing conditions, thereby likely promoting gas and nutrient transport within the "canopy." This mixing could have served to enhance local-scale oxygen concentrations and redistribute resources like dissolved organic carbon. Our work suggests that Ediacaran marine animal forests may have contributed to the ventilation of the oceans over 560 million years ago, well before the Cambrian explosion of animals.

Foraging mechanisms in excavate flagellates shed light on the functional ecology of early eukaryotes.

The phagotrophic flagellates described as "typical excavates" have been hypothesized to be morphologically similar to the Last Eukaryotic Common Ancestor and understanding the functional ecology of excavates may therefore help shed light on the ecology of these early eukaryotes. Typical excavates are characterized by a posterior flagellum equipped with a vane that beats in a ventral groove. Here, we combined flow visualization and observations of prey capture in representatives of the three clades of excavates with computational fluid dynamic modeling, to understand the functional significance of this cell architecture. We record substantial differences amongst species in the orientation of the vane and the beat plane of the posterior flagellum. Clearance rate magnitudes estimated from flow visualization and modeling are both like that of other similarly sized flagellates. The interaction between a vaned flagellum beating in a confinement is modeled to produce a very efficient feeding current at low energy costs, irrespective of the beat plane and vane orientation and of all other morphological variations. Given this predicted uniformity of function, we suggest that the foraging systems of typical excavates studied here may be good proxies to understand those potentially used by our distant ancestors more than 1 billion years ago.

A finite element analysis of periodontal ligament fluid mechanics response to occlusal loading based on hydro-mechanical coupling model.

OBJECTIVE: Considering fluid stimulation is one of the essential biomechanical signals for periodontal tissues, this study aims to characterizing fluid mechanics response during occlusal loading by a hydro-mechanical coupling model for periodontal ligament. DESIGN: Models simulating periodontium with normal bone height and with intraosseous defects were built with three mechanical modules: tooth, periodontal ligament and alveolar bone. Tooth was modeled as linear elastic, and periodontal ligament and alveolar bone as a hydro-mechanical coupling model. Transient analyses under dynamic occlusal loading were performed. Fluid dynamics within periodontal ligament space was simulated and visualized by post-processing module. RESULTS: Reciprocating oscillatory flow occurred within the periodontal ligament under occlusal loading. Higher pore pressure and fluid velocity were observed in furcation and apical regions compared to mid-root and cervical regions. Intraosseous defects increased pore pressure and fluid velocity within the periodontal ligament, most significantly near the defect. CONCLUSION: Based on the results of the hydro-mechanical coupling model, significant oscillatory fluid motion is observed within the periodontal ligament under occlusal loading. Particularly, higher fluid velocity is evident in the furcation and apical areas. Additionally, Intraosseous defects significantly enhance fluid motion within the periodontal ligament.

Transport and retention of sinking microplastics in a well-mixed estuary.

Estuaries have been shown to be potential hotspots of microplastic accumulation, but the hydrodynamic conditions and particle properties that control this process need further investigation. We have designed a series of numerical particle-tracking experiments to examine the sensitivity of retention in estuaries to particle size, particle density and varying tides and freshwater flow. At the end of the simulation, over 90 % of sinking particles are retained in the estuary, and the retention rate is further increased by high river runoff. In contrast, increased river discharge increases the number of marginally-buoyant (i.e. density close to estuarine water) particles that escape the estuary. Larger particle size tends to limit the downstream transport of sinking particles but can facilitate the transport of marginally-buoyant particles. Tidal asymmetry, vertical turbulent mixing and the vertical structure of the subtidal circulation are proposed as the underlying mechanisms controlling the fate of particles.

Engineered shapes using electrohydrodynamic atomization for an improved drug delivery.

The shapes of micro- and nano-products have profound influences on their functional performances, which has not received sufficient attention during the past several decades. Electrohydrodynamic atomization (EHDA) techniques, mainly include electrospinning and electrospraying, are facile in manipulate their products' shapes. In this review, the shapes generated using EHDA for modifying drug release profiles are reviewed. These shapes include linear nanofibers, round micro-/nano-particles, and beads-on-a-string hybrids. They can be further divided into different kinds of sub-shapes, and can be explored for providing the desired pulsatile release, sustained release, biphasic release, delayed release, and pH-sensitive release. Additionally, the shapes resulted from the organizations of electrospun nanofibers are discussed for drug delivery, and the shapes and inner structures can be considered together for developing novel drug delivery systems. In future, the shapes and the related shape-performance relationships at nanoscale, besides the size, inner structure and the related structure-performance relationships, would further play their important roles in promoting the further developments of drug delivery field. This article is categorized under: Therapeutic Approaches and Drug Discovery > Emerging Technologies.

Real-time position and pose prediction for a self-propelled undulatory swimmer in 3D space with artificial lateral line system.

This study aims to investigate the feasibility of using an artificial lateral line (ALL) system for predicting the real-time position and pose of an undulating swimmer with Carangiform swimming patterns. We established a 3D computational fluid dynamics simulation to replicate the swimming dynamics of a freely swimming mackerel under various motion parameters, calculating the corresponding pressure fields. Using the simulated lateral line data, we trained an artificial neural network to predict the centroid coordinates and orientation of the swimmer. A comprehensive analysis was further conducted to explore the impact of sensor quantity, distribution, noise amplitude and sampling intervals of the ALL array on predicting performance. Additionally, to quantitatively assess the reliability of the localization network, we trained another neural network to evaluate error magnitudes for different input signals. These findings provide valuable insights for guiding future research on mutual sensing and schooling in underwater robotic fish.

Critical hydrodynamic force levels for efficient removal of oral biofilms in simulated interdental spaces.

OBJECTIVES: Sonic toothbrushes generate hydrodynamic shear forces for oral biofilm removal on tooth surfaces, but the effective thresholds for biofilm removal remain unexplored. This in vitro study aimed to investigate various threshold values for hydrodynamic biofilm removal in vitro. MATERIALS AND METHODS: A specialized test bench was designed with a known water flow field within a gap, ensuring that hydrodynamic shear forces on the wall were solely dependent on the volume flow, which was quantifiable using an integrated flow meter and proven by a computational fluid dynamics simulation. A young 20 h supragingival six-species biofilm was developed on hydroxyapatite disks (∅ 5 mm) and applied into the test bench, subjecting them to ascending force levels ranging from 0 to 135 Pa. The remaining biofilms were quantified using colony forming units (CFU) and subjected to statistical analysis through one-way ANOVA. RESULTS: Volume flow measures < 0.1 l/s: Error 1% of reading were established with the test bench. Untreated biofilms (0 Pa, no hydrodynamic shear forces) reached 7.7E7 CFU/harvest and differed significantly from all treated biofilm groups. CFU reductions of up to 2.3E6 were detected using 20 Pa, and reductions of two orders of magnitude were reached above wall shear forces of 45 Pa (6.9E5). CONCLUSIONS: Critical hydrodynamic force levels of at least 20 Pa appear to be necessary to have a discernible impact on initial biofilm removal. CLINICAL RELEVANCE: Pure hydrodynamic forces alone are insufficient for adequate biofilm removal. The addition of antiseptics is essential to penetrate and disrupt hydrodynamically loosened biofilm structures effectively.

Hydrodynamic performance of Ordovician archaeostracan carapaces.

The diversification of macroscopic pelagic arthropods such as caryocaridid archaeostracans was a crucial aspect of the Great Ordovician Biodiversification Event, and the plankton revolution. A pelagic mode of life has been inferred for caryocaridids from their common presence in black graptolitic shales alongside carapace morphologies that appear streamlined. However, the hydrodynamic performance within the group and comparisons with other archaeostracans were lacking. Here we use a computational fluid dynamics approach to quantify the hydrodynamic performance of caryocaridids, and other early Palaeozoic archaeostracans including Arenosicaris inflata and Ordovician ceratiocaridids. We show that streamlining of the carapace was an important factor facilitating a pelagic mode of life in caryocaridids, in reducing the drag coefficient and facilitating a broader range of lift coefficients at different angles of attack. However, comparable hydrodynamic performance is also recovered for some ceratiocaridids. This suggests that alongside carapace streamlining, adaptations to appendages and thinning of the carapace were also important for a pelagic mode of life in Ordovician caryocaridids.

Effect of basket mesh size on the hydrodynamics of a partially filled (500 mL) USP rotating basket dissolution testing Apparatus 1.

The USP Rotating Basket Dissolution Testing Apparatus 1 is listed in the USP as one of the tools to assess dissolution of oral solid dosage forms. Baskets of different mesh sizes can be used to differentiate between dissolution profiles of different formulations. Here, we used Particle Image Velocimetry (PIV) to study the hydrodynamics of the USP Apparatus 1 using baskets with different mesh openings (10-, 20- and 40-mesh) revolving at 100 rpm, when the vessel was filled with 500 mL. The velocity profiles throughout the liquid were found to vary significantly using baskets of different mesh sizes, typically increasing with increased size of the opening of the basket mesh, especially for axial and radial velocities. This, in turn, resulted in a significantly different flow rate through the basket, which can be expected to significantly impact the dissolution rate of the drug product. A comparison between the results of this work with those of a previous study with a 900-mL fill volume (Sirasitthichoke et al., Intern. J. Pharmaceutics, 2021, 607: 120976), shows that although the hydrodynamics in the USP Apparatus 1 changed with fill level in the vessel, the flow rate through the basket was not significantly affected. This implies that tablets dissolving in the two systems would experience similar tablet-liquid medium mass transfer coefficients, and therefore similar initial dissolution rates, but different dissolution profiles because of the difference in volume.

Numerical simulation of the influence of nasal cycle on nasal airflow.

To study the characteristics of nasal airflow in the presence of nasal cycle by computational fluid dynamics. CT scan data of a healthy Chinese individual was used to construct a three-dimensional model of the nasal cavity to be used as simulation domain. A sinusoidal airflow velocity is set at the nasal cavity entrance to reproduce the breathing pattern of a healthy human. There was a significant difference in the cross-sectional area between the two sides of the nasal cavity. Particularly, the decongested side is characterized by a larger cross-section area, and consequently, by a larger volume with respect to the congested side. The airflow velocity, pressure, and nasal resistance were higher on the congested narrow side. The temperature regulation ability on the congested narrow side was stronger than that on the decongested wider side. During the nasal cycle, there are differences in the nasal cavity function between the congested and decongested sides. Therefore, when evaluating the impact of various factors on nasal cavity function, the nasal cycle should be considered.

An adjoint-based approach for the surgical correction of nasal septal deviations.

Deviations of the septal wall are widespread anatomic anomalies of the human nose; they vary significantly in shape and location, and often cause the obstruction of the nasal airways. When severe, septal deviations need to be surgically corrected by ear-nose-throat (ENT) specialists. Septoplasty, however, has a low success rate, owing to the lack of suitable standardized clinical tools for assessing type and severity of obstructions, and for surgery planning. Moreover, the restoration of a perfectly straight septal wall is often impossible and possibly unnecessary. This paper introduces a procedure, based on advanced patient-specific Computational Fluid Dynamics (CFD) simulations, to support ENT surgeons in septoplasty planning. The method hinges upon the theory of adjoint-based optimization, and minimizes a cost function that indirectly accounts for viscous losses. A sensitivity map is computed on the mucosal wall to provide the surgeon with a simple quantification of how much tissue removal at each location would contribute to easing the obstruction. The optimization procedure is applied to three representative nasal anatomies, reconstructed from CT scans of patients affected by complex septal deviations. The computed sensitivity consistently identifies all the anomalies correctly. Virtual surgery, i.e. morphing of the anatomies according to the computed sensitivity, confirms that the characteristics of the nasal airflow improve significantly after small anatomy changes derived from adjoint-based optimization.

An individualised 3D computational flow and particle model to predict the deposition of inhaled medicines - A case study using a nebuliser.

BACKGROUND AND OBJECTIVE: Drug inhalation is generally accepted as the preferred administration method for treating respiratory diseases. To achieve effective inhaled drug delivery for an individual, it is necessary to use an interdisciplinary approach that can cope with inter-individual differences. The paper aims to present an individualised pulmonary drug deposition model based on Computational Fluid and Particle Dynamics simulations within a time frame acceptable for clinical use. METHODS: We propose a model that can analyse the inhaled drug delivery efficiency based on the patient's airway geometry as well as breathing pattern, which has the potential to also serve as a tool for a sub-regional diagnosis of respiratory diseases. The particle properties and size distribution are taken for the case of drug inhalation by using nebulisers, as they are independent of the patient's breathing pattern. Finally, the inhaled drug doses that reach the deep airways of different lobe regions of the patient are studied. RESULTS: The numerical accuracy of the proposed model is verified by comparison with experimental results. The difference in total drug deposition fractions between the simulation and experimental results is smaller than 4.44% and 1.43% for flow rates of 60 l/min and 15 l/min, respectively. A case study involving a COVID-19 patient is conducted to illustrate the potential clinical use of the model. The study analyses the drug deposition fractions in relation to the breathing pattern, aerosol size distribution, and different lobe regions. CONCLUSIONS: The entire process of the proposed model can be completed within 48 h, allowing an evaluation of the deposition of the inhaled drug in an individual patient's lung within a time frame acceptable for clinical use. Achieving a 48-hour time window for a single evaluation of patient-specific drug delivery enables the physician to monitor the patient's changing conditions and potentially adjust the drug administration accordingly. Furthermore, we show that the proposed methodology also offers a possibility to be extended to a detection approach for some respiratory diseases.

Modelling lower-limb peripheral arterial disease using clinically available datasets: impact of inflow boundary conditions on hemodynamic indices for restenosis prediction.

BACKGROUND AND OBJECTIVES: The integration of hemodynamic markers as risk factors in restenosis prediction models for lower-limb peripheral arteries is hindered by fragmented clinical datasets. Computed tomography (CT) scans enable vessel geometry reconstruction and can be obtained at different times than the Doppler ultrasound (DUS) images, which provide information on blood flow velocity. Computational fluid dynamics (CFD) simulations allow the computation of near-wall hemodynamic indices, whose accuracy depends on the prescribed inlet boundary condition (BC), derived from the DUS images. This study aims to: (i) investigate the impact of different DUS-derived velocity waveforms on CFD results; (ii) test whether the same vessel areas, subjected to altered hemodynamics, can be detected independently of the applied inlet BC; (iii) suggest suitable DUS images to obtain reliable CFD results. METHODS: CFD simulations were conducted on three patients treated with bypass surgery, using patient-specific DUS-derived inlet BCs recorded at either the same or different time points than the CT scan. The impact of the chosen inflow condition on bypass hemodynamics was assessed in terms of wall shear stress (WSS)-derived quantities. Patient-specific critical thresholds for the hemodynamic indices were applied to identify critical luminal areas and compare the results with a reference obtained with a DUS image acquired in close temporal proximity to the CT scan. RESULTS: The main findings indicate that: (i) DUS-derived inlet velocity waveforms acquired at different time points than the CT scan led to statistically significantly different CFD results (p<0.001); (ii) the same luminal surface areas, exposed to low time-averaged WSS, could be identified independently of the applied inlet BCs; (iii) similar outcomes were observed for the other hemodynamic indices if the prescribed inlet velocity waveform had the same shape and comparable systolic acceleration time to the one recorded in close temporal proximity to the CT scan. CONCLUSIONS: Despite a lack of standardised data collection for diseased lower-limb peripheral arteries, an accurate estimation of luminal areas subjected to altered near-wall hemodynamics is possible independently of the applied inlet BC. This holds if the applied inlet waveform shares some characteristics - derivable from the DUS report - as one matching the acquisition time of the CT scan.