A Novel Size-Based Centrifugal Microfluidic Design to Enrich and Magnetically Isolate Circulating Tumor Cells from Blood Cells through Biocompatible Magnetite-Arginine Nanoparticles.

This paper presents a novel centrifugal microfluidic approach (so-called lab-on-a-CD) for magnetic circulating tumor cell (CTC) separation from the other healthy cells according to their physical and acquired chemical properties. This study enhances the efficiency of CTC isolation, crucial for cancer diagnosis, prognosis, and therapy. CTCs are cells that break away from primary tumors and travel through the bloodstream; however, isolating CTCs from blood cells is difficult due to their low numbers and diverse characteristics. The proposed microfluidic device consists of two sections: a passive section that uses inertial force and bifurcation law to sort CTCs into different streamlines based on size and shape and an active section that uses magnetic forces along with Dean drag, inertial, and centrifugal forces to capture magnetized CTCs at the downstream of the microchannel. The authors designed, simulated, fabricated, and tested the device with cultured cancer cells and human cells. We also proposed a cost-effective method to mitigate the surface roughness and smooth surfaces created by micromachines and a unique pulsatile technique for flow control to improve separation efficiency. The possibility of a device with fewer layers to improve the leaks and alignment concerns was also demonstrated. The fabricated device could quickly handle a large volume of samples and achieve a high separation efficiency (93%) of CTCs at an optimal angular velocity. The paper shows the feasibility and potential of the proposed centrifugal microfluidic approach to satisfy the pumping, cell sorting, and separating functions for CTC separation.

Spacing effects on flows around two square cylinders in staggered arrangement via LBM.

This study presents a computational analysis of fluid flow characteristics around two staggered arranged square cylinders using the Lattice Boltzmann Method (LBM). With Reynolds number (Re) fixed at 200, numerical simulations explore the influence of varying gap ratios (G) ranging from 0 to 10 times the cylinder size. Emphasis is placed on understanding the impact of cylinders spacing on flow structure mechanisms and induced forces. Investigation of fluid flow parameters includes vorticity behavior, pressure streamlines, and variations in drag and lift coefficients alongside the Strouhal number under different values of G. From the results, four distinct flow patterns emerge: single bluff body flow, flip flopping flow, modulated synchronized flow, and synchronized flow, each exhibiting unique characteristics. This study reveals the strong dependence of fluid forces on G, with low spacing values leading to complex vortex structures and fluctuating forces influenced by jet flow effects. At higher spacing values, proximity effects between cylinders diminish, resulting in a smoother periodic flow. The Strouhal number, average drag force and the rms values of drag and lift force coefficients vary abruptly at narrow gaps and become smooth at higher gap ratios. Unlike the tandem and side-by-side arrangements the staggered cylinders arrangement is found to have significant impact on the pressure variations around both cylinders. Overall, this research could contribute to a comprehensive understanding of staggered cylinder arrangements and their implications for engineering applications.

The effects of drag models and operating condition on the simulation of photocatalytic degradation of pesticides in a fluidized bed reactor.

Fluidized bed reactor can enhance mass transfer and increase reaction rate. Numerical simulation helps to optimize fluidized bed reactors. The present paper models the photocatalytic oxidation of 2,4-dichlorophenoxy acetic acid in a fluidized bed reactor using the Eulerian-Eulerian model. The drag models have influences on the distribution of catalysts particles. The bed expands under the fluid flow and reaches a quasi-steady height at approximately 3s. The asymmetric distribution of catalysts with respect to the axis plane is predicted. The Gidaspow model predicts the nearly same bed expansion with the experimental data, whereas the Syamlal-O ' Brien model overestimates it. The simulation results at two pH values are in accordance with the experimental data. The removal in the continuous stirred tank reactor decreases with the increase in flow rate.

Towards an accurate rolling resistance: Estimating intra-cycle load distribution between front and rear wheels during wheelchair propulsion from inertial sensors.

Accurate assessment of rolling resistance is important for wheelchair propulsion analyses. However, the commonly used drag and deceleration tests are reported to underestimate rolling resistance up to 6% due to the (neglected) influence of trunk motion. The first aim of this study was to investigate the accuracy of using trunk and wheelchair kinematics to predict the intra-cyclical load distribution, more particularly front wheel loading, during hand-rim wheelchair propulsion. Secondly, the study compared the accuracy of rolling resistance determined from the predicted load distribution with the accuracy of drag test-based rolling resistance. Twenty-five able-bodied participants performed hand-rim wheelchair propulsion on a large motor-driven treadmill. During the treadmill sessions, front wheel load was assessed with load pins to determine the load distribution between the front and rear wheels. Accordingly, a machine learning model was trained to predict front wheel load from kinematic data. Based on two inertial sensors (attached to the trunk and wheelchair) and the machine learning model, front wheel load was predicted with a mean absolute error (MAE) of 3.8% (or 1.8 kg). Rolling resistance determined from the predicted load distribution (MAE: 0.9%, mean error (ME): 0.1%) was more accurate than drag test-based rolling resistance (MAE: 2.5%, ME: -1.3%).

How Magnetic Erosion Affects the Drag-Based Kinematics of Fast Coronal Mass Ejections.

In order to advance our understanding of the dynamic interactions between coronal mass ejections (CMEs) and the magnetized solar wind, we investigate the impact of magnetic erosion on the well-known aerodynamic drag force acting on CMEs traveling faster than the ambient solar wind. In particular, we start by generating empirical relationships for the basic physical parameters of CMEs that conserve their mass and magnetic flux. Furthermore, we examine the impact of the virtual mass on the equation of motion by studying a variable-mass system. We next implement magnetic reconnection into CME propagation, which erodes part of the CME magnetic flux and outer-shell mass, on the drag acting on CMEs, and we determine its impact on their time and speed of arrival at 1 AU. Depending on the strength of the magnetic erosion, the leading edge of the magnetic structure can reach near-Earth space up to ≈ three hours later, compared to the non-eroded case. Therefore, magnetic erosion may have a significant impact on the propagation of fast CMEs and on predictions of their arrivals at 1 AU. Finally, the modeling indicates that eroded CMEs may experience a significant mass decrease. Since such a decrease is not observed in the corona, the initiation distance of erosion may lie beyond the field-of-view of coronagraphs (i.e. 30R⊙).

Vortex structure and fluid force changed by altering whole-body kinematic parameters during underwater undulatory swimming.

Swimmers generate vortices around their bodies during underwater undulatory swimming (UUS). Alteration of UUS movement would induce changes in vortex structure and fluid force. This study investigated whether a skilled swimmer's movement generated an effective vortex and fluid force for increasing the UUS velocity. A three-dimensional digital model and kinematic data yielded during UUS with maximum effort were collected for one skilled and one unskilled swimmer. The skilled swimmer's UUS kinematics were input into the skilled swimmer's model (SK-SM) and unskilled swimmer's model (SK-USM), followed by the kinematics of the unskilled swimmer (USK-USM and USK-SM, respectively). The vortex area, circulation, and peak drag force were determined using computational fluid dynamics. A larger vortex with greater circulation at the ventral side of the trunk and a greater circulation vortex behind the swimmer were observed in SK-USM compared to USK-USM. USK-SM generated a smaller vortex on the ventral side of the trunk and behind the swimmer, with a weaker circulation behind the swimmer compared to SK-SM. The peak drag force was larger for SK-USM than for USK-USM. Our results indicate that an effective vortex for propulsion was generated when a skilled swimmer's UUS kinematics was input in the other swimmer's model.

A novel approach to directly measuring wheel and caster rolling resistance accurately predicts user-wheelchair system-level rolling resistance.

Introduction: Clinical practice guidelines for preservation of upper extremity recommend minimizing wheelchair propulsion forces. Our ability to make quantitative recommendations about the effects of wheelchair configuration changes is limited by system-level tests to measure rolling resistance (RR). We developed a method that directly measures caster and propulsion wheel RR at a component-level. The study purpose is to assess accuracy and consistency of component-level estimates of system-level RR. Methods: The RR of N = 144 simulated unique wheelchair-user systems were estimated using our novel component-level method and compared to system-level RR measured by treadmill drag tests, representing combinations of caster types/diameters, rear wheel types/diameters, loads, and front-rear load distributions. Accuracy was assessed by Bland-Altman limits of agreement (LOA) and consistency by intraclass correlation (ICC). Results: Overall ICC was 0.94, 95% CI [0.91-0.95]. Component-level estimates were systematically lower than system-level (-1.1 N), with LOA +/-1.3 N. RR force differences between methods were constant over the range of test conditions. Conclusion: Component-level estimates of wheelchair-user system RR are accurate and consistent when compared to a system-level test method, evidenced by small absolute LOA and high ICC. Combined with a prior study on precision, this study helps to establish validity for this RR test method.

Dynamics and hydrodynamic efficiency of diving beetle while swimming.

Diving beetle, an excellent biological prototype for bionic underwater vehicles, can achieve forward swimming, backward swimming, and flexible cornering by swinging its two powerful hind legs. An in-depth study of the propulsion performance of them will contribute to the micro underwater vehicles. In this paper, the kinematic and dynamic parameters, and the hydrodynamic efficiency of the diving beetle are studied by analysis of swimming videos using Motion Capture Technology, combined with CFD simulations. The results show that the hind legs of diving beetle can achieve high propulsion force and low return resistance during one propulsion cycle at both forward and backward swimming modes. The propulsion efficiencies of forward and backward swimming are 0.47 and 0.30, respectively. Although the efficiency of backward swimming is lower, the diving beetle can reach a higher speed in a short time at this mode, which can help it avoid natural enemies. At backward swimming mode, there is a long period of passive swing of hind legs, larger drag exists at higher speed during the recovery stroke, which reduces the propulsion efficiency to a certain extent. Reasonable planning of the swing speed of the hind legs during the power stroke and the recovery stroke can obtain the highest propulsion efficiency of this propulsion method. This work will be useful for the development of a bionic propulsion system of micro underwater vehicle.

A Simple Approach to Estimate the Drag Coefficients of a Submerged Floater.

The calculation of the drag force is a fundamental requirement in the design of any submerged system intended for marine exploration. The calculation can be performed by analytic analysis, numerical modeling, or by a direct calculation performed in a designated testing facility. However, for complex structures and especially those with a non-rigid design, the analytic and numerical analyses are not sufficiently accurate, while the direct calculation is a costly operation. In this paper, we propose a simple approach for how to calculate the drag coefficient in-situ. Aimed specifically at the complex case of elastic objects whose modeling via Computer-Aided Design (CAD) is challenging, our approach evaluates the relation between the object's speed at steady-state and its mass to extract the drag coefficient in any desired direction, the hydro-static force, and, when relevant, also the thruster's force. We demonstrate our approach for the special case of a highly complex elastic-shaped floater that profiles the water column. The analysis of two such floaters in two different sea environments shows accurate evaluation results and supports our claim for robustness. In particular, the simplicity of the approach makes it appealing for any arbitrary shaped object.

Three-dimensional modeling and experimentation of microfluidic devices driven by surface acoustic wave.

Surface acoustic wave (SAW) technology is proving to be an effective tool for manipulating micro-nano particles. In this paper, we present a fully-coupled 3D model of standing SAW acoustofluidic devices for obtaining particle motion. The "improved limiting velocity method" (ILVM) was used to investigate the distribution of acoustic pressure and acoustic streaming in microchannel. The results show that the distribution of acoustic pressure and acoustic streaming on the piezoelectric substrate surface perpendicular to the acoustic wave propagation direction is inhomogeneous. The motion of micro-particles with diameters of 0.5-, 5-, and 10 µm is then simulated to investigate the interaction of acoustic radiation force and drag force caused by pressure and acoustic streaming. We demonstrate that micro and nanoparticles can move in three dimensions when acoustic radiation force and acoustic streaming interact. This result and method are critical for designing SAW microfluidic chips and controlling particle motion precisely.

Membrane-organic solute interactions in asymmetric flow field flow fractionation: Interplay of hydrodynamic and electrostatic forces.

The structure and size characterization of organic matter (OM) using flow field-flow fractionation (FFFF) is interesting due to the numerous interactions of OM in aquatic systems and water treatment processes. The estimation of hydrodynamic and electrostatic forces involved in the fractionation of OM over different molecular weight cut-off (MWCO) membranes is vital for a better understanding of the FFFF process. This work aims to understand the membrane-OM interactive forces with respect to membrane MWCO, solute molecular weight, flow rates, solution pH and ionic strength. Polystyrene sulfonate sodium salt (PSS) of molecular weights 10, 30 and 65 kDa were used as model organic solutes for fractionation over ultrafiltration (UF) membranes of MWCO 1-30 kDa. Maximum fractionation of PSS was achieved by using a tight membrane of 1 kDa MWCO at the conditions of high permeate flow rate (1.5-2.0 mL·min-1), low concentrate flow rate (0.2-0.3 mL·min-1) and low ionic strength (10 mM). The better fractionation corresponds to high permeate drag force and low concentrate drag force. A low membrane-solute DLVO interaction is favourable for the retention of a small solute. This study illustrated that FFFF characteristics can be analyzed based on membrane-solute interactive forces controlled by selected flow, size and charge parameters.

Force and Velocity Analysis of Particles Manipulated by Toroidal Vortex on Optoelectrokinetic Microfluidic Platform.

The rapid electrokinetic patterning (REP) technique has been demonstrated to enable dynamic particle manipulation in biomedical applications. Previous studies on REP have generally considered particles with a size less than 5 µm. In this study, a REP platform was used to manipulate polystyrene particles with a size of 3~11 µm in a microfluidic channel sandwiched between two ITO conductive glass plates. The effects of the synergy force produced by the REP electrothermal vortex on the particle motion were investigated numerically for fixed values of the laser power, AC driving voltage, and AC driving frequency, respectively. The simulation results showed that the particles were subject to a competition effect between the drag force produced by the toroidal vortex, which prompted the particles to recirculate in the bulk flow adjacent to the laser illumination spot on the lower electrode, and the trapping force produced by the particle and electrode interactions, which prompted the particles to aggregate in clusters on the surface of the illuminated spot. The experimental results showed that as the laser power increased, the toroidal flow range over which the particles circulated in the bulk flow increased, while the cluster range over which the particles were trapped on the electrode surface reduced. The results additionally showed that the particle velocity increased with an increasing laser power, particularly for particles with a smaller size. The excitation frequency at which the particles were trapped on the illuminated hot-spot reduced as the particle size increased. The force and velocity of polystyrene particles by the REP toroidal vortex has implications for further investigating the motion behavior at the biological cell level.

Aging-related structural change in 3D extracellular matrix affects its mechanics.

Aging induces marked alterations in the structural, mechanical, and transport properties in the extracellular matrix (ECM). To provide computational data on the impact of aging-related changes on ECM mechanical quantities and transport properties, we developed a computational model for the aging-related ECM fibrous network. A finite volume method was utilized to calculate the velocity field, pressure loss, hydraulic conductivity and drag force. Our results quantitatively demonstrated that the hydraulic conductivity in most of the aging ECM-mimetic fibrous networks tends to be significantly lower than young ECM-mimetic fibrous networks, while pressure loss and drag force show the opposite trend. All these findings highlight that such altered mechanical quantities and transport properties during aging can be important biological cues to assess the aging process and eventually provide insights in treating aging-related diseases.

Analysis of Movements and Behavior of Bighead Carps (Hypophthalmichthys nobilis) Considering Fish Passage Energetics in an Experimental Vertical Slot Fishway.

An understanding of fish movement behavior in response to flow field variables is important for exploring the hydrodynamic strategies of fish in fish passages. In this paper, bighead carps were taken as an example. The fish movement behavior response to water flow field information by means of estimating the energetic expenditure using an IBM approach in an experimental fishway was investigated. Fish swimming velocity, drag force, and energy expenditure were analyzed in varied flow conditions related to hydraulic variables, including velocity (V), turbulent kinetic energy (TKE), and strain rate (SR). The result indicated that the fish will require more energy in high TKE zones. This study provides a reference for optimizing the design of fish passages and fisheries management. This method can be applied to assess the efficiency of fish bypass structures and conduct fish survival studies.

Quantitative Analysis of Drag Force for Task-Specific Micromachine at Low Reynolds Numbers.

Micromotors have spread widely in order to meet the needs of new applications, including cell operation, drug delivery, biosensing, precise surgery and environmental decontamination, due to their small size, low energy consumption and large propelling power, especially the newly designed multifunctional micromotors that combine many extra shape features in one device. Features such as rod-like receptors, dendritic biosensors and ball-like catalyzing enzymes are added to the outer surface of the tubular micromotor during fabrication to perform their special mission. However, the structural optimization of motion performance is still unclear. The main factor restricting the motion performance of the micromotors is the drag forces. The complex geometry of a micromotor makes its dynamic behavior more complicated in a fluid environment. This study aimed to design the optimum structure of tubular micromotors with minimum drag forces and obtain the magnitude of drag forces considering both the internal and external fluids of the micromotors. By using the computational fluid dynamics software Fluent 18.0 (ANSYS), the drag force and the drag coefficient of different conical micromotors were calculated. Moreover, the influence of the Reynolds numbers Re, the semi-cone angle δ and the ratios ξ and η on the drag coefficient was analyzed. The results show the drag force monotonically increased with Reynolds numbers Re and the ratio η. The extreme point of the drag curve is reached when the semi-cone angle δ is 8° and the ratio ξ is 3.846. This work provides theoretical support and guidance for optimizing the design and development of conical micromotors.

A real-time and modular approach for quick detection and mechanism exploration of DPIs with different carrier particle sizes.

Dry powder inhalers (DPIs) had been widely used in lung diseases on account of direct pulmonary delivery, good drug stability and satisfactory patient compliance. However, an indistinct understanding of pulmonary delivery processes (PDPs) hindered the development of DPIs. Most current evaluation methods explored the PDPs with over-simplified models, leading to uncompleted investigations of the whole or partial PDPs. In the present research, an innovative modular process analysis platform (MPAP) was applied to investigate the detailed mechanisms of each PDP of DPIs with different carrier particle sizes (CPS). The MPAP was composed of a laser particle size analyzer, an inhaler device, an artificial throat and a pre-separator, to investigate the fluidization and dispersion, transportation, detachment and deposition process of DPIs. The release profiles of drug, drug aggregation and carrier were monitored in real-time. The influence of CPS on PDPs and corresponding mechanisms were explored. The powder properties of the carriers were investigated by the optical profiler and Freeman Technology four powder rheometer. The next generation impactor was employed to explore the aerosolization performance of DPIs. The novel MPAP was successfully applied in exploring the comprehensive mechanism of PDPs, which had enormous potential to be used to investigate and develop DPIs.

Aerodynamic interaction between in-line runners: new insights on the drafting strategy in running.

This study aims at modelling the aerodynamic interaction between a world-class runner and several pacers running in line, the objective being to determine the best drafting position in terms of potential speed gain and running time. Computational Fluid Dynamics calculations were performed to determine the aerodynamic drag forces exerted on the runners. Secondly, we estimated the metabolic savings for each of the runner's positions in the pack of pacers. Finally, we estimated a potential equivalent running speed and a corresponding running time gain for each of the runner's positions. Running second in a group of 5 runners would be the most effective drafting strategy, with a reduction of aerodynamic drag reaching 63.3%, corresponding to an improvement of 4.4% of the running economy. Furthermore, decreasing the drag forces acting on the runner would results in greater running speed. For example, a 63.3% reduction in the aerodynamic forces yields a 2.9% faster running speed (20.7 to 21.31 km/h) for an equal metabolic cost. Running in the wake of a leading runner (in a pack of five runners) over a marathon is estimated to provide a significant time saving of 3 min 28 s.

Kinetic Monte Carlo method applied to micrometric particle detachment mechanisms by aerodynamic forces.

The formulation of a kinetic Monte Carlo simulation to account for the different possible mechanisms present in the problem of resuspension of aerosol particles is addressed as an extension of a former model Benitoet al(2016J. Aerosol Sci.10026-37). The re-entrainment of micrometer particles to airflow when detached from a surface by aerodynamic forces is modeled using the similitude of the problem with the desorption process from heterogeneous surfaces. Depending on the relative role of the intervening forces, three main mechanisms for movement initiation can be present:rolling,slidingandlifting-off. Three different transition probabilities are defined for each mechanism and the corresponding transition rates calculated for the kinetic process to be simulated. The decisive factor for the development of the model is to set an appropriate dynamical hierarchy to simulate correctly the evolution of the transition rates as the airflow velocity increases, reflecting the stochastic nature of the process, not always fully captured by other Monte Carlo approaches. The model is applied to spherical and elongated particles on a flat surface, reproducing qualitatively well the experimental trends found by other authors for the case of particles with different shapes. It is also demonstrated that, for elongated particles, the main mechanism assisting the detachment is not rolling but sliding, underscoring the need for an adequate choice of the particles shape and detachment mechanism when looking for the critical conditions for particle removal from surfaces.

Flying Squirrel-Inspired Motion Control of a Light-Deformed Pt-PAzoMA Micromotor through Drag Force Manipulation.

Micromotors require stable and precise motion control for complex tasks such as microsurgery, drug delivery in vivo, or environmental monitoring ex vivo. However, a continuous control signal is needed for micromotors to achieve motion control during their whole journey, which hinders their application in areas where external control stimuli are limited or unavailable. Fortunately, nature suggests an excellent solution that flying squirrel exhibits motion tuning capability by deforming itself when jumping off a tall tree. Inspired by this, we propose a Pt-PAzoMA Janus micromotor that precisely changes its shape (from a spherical shape to an elliptical shape) under a brief light signal (450 nm) and maintains this deformation until next signal reception. The deformed elliptical micromotor performs relatively low-speed motion compared to the spherical one, which is further confirmed by massive simulation results. In addition, by investigating motion behavior experimentally and theoretically, it is proved that the motion modulation is caused by the drag force changing brought from the deformation. This method represents a different route to regulate the motion of micromotors without a continuous signal, which is useful in application scenarios where the environmental control signal is inaccessible/limited or long-time operation with minimum energy input is required to maintain motion manipulation. With further function modification, this kind of shape-changing micromotor has potential in optimizing drug diffusion efficiency by speed altering and long-term monitoring at the diseased area by confining the active range of the micromotor in the targeted area through deformation.

Numerical study of acoustophoretic manipulation of particles in microfluidic channels.

The microfluidic technology based on surface acoustic waves (SAW) has been developing rapidly, as it can precisely manipulate fluid flow and particle motion at microscales. We hereby present a numerical study of the transient motion of suspended particles in a microchannel. In conventional studies, only the microchannel's bottom surface generates SAW and only the final positions of the particles are analyzed. In our study, the microchannel is sandwiched by two identical SAW transducers at both the bottom and top surfaces while the channel's sidewalls are made of poly-dimethylsiloxane (PDMS). Based on the perturbation theory, the suspended particles are subject to two types of forces, namely the Acoustic Radiation Force (ARF) and the Stokes Drag Force (SDF), which correspond to the first-order acoustic field and the second-order streaming field, respectively. We use the Finite Element Method (FEM) to compute the fluid responses and particle trajectories. Our numerical model is shown to be accurate by verifying against previous experimental and numerical results. We have determined the threshold particle size that divides the SDF-dominated regime and the ARF-dominated regime. By examining the time scale of the particle movement, we provide guidelines on the device design and operation.