Dilute Suspensions Research Articles

Statin-treated RBC dynamics in a microfluidic porous-like network

The impact of therapeutic interventions on red blood cell (RBC) deformability and microscale transport is investigated, using statins as an exemplar. Human RBCs were treated in vitro with two commonly prescribed statins, atorvastatin and rosuvastatin, at clinically relevant concentrations. Changes in RBC deformability were quantified using a microfluidic-based ektacytometer and expressed in terms of the elongation index. Dilute suspensions of the statin-treated RBCs were then perfused through a microfluidic pillar array, at a constant flow rate and negligible inertia, and imaged. Particle Tracking Velocimetry (PTV) was applied to track RBCs, identify preferential paths and estimate their velocities, whereas image processing was used to estimate cell dynamics, perfusion metrics and distributions. The findings were compared against those of healthy, untreated cells. Statins enhanced RBC deformability in agreement with literature. The extent of enhancement was found to be statin-dependent. The softer statin-treated cells were found to flow in straight, less tortuous paths, spend more time inside the pillar array and exhibit lower velocities compared to healthy RBCs, attributed to their enhanced deformation and longer shape recovery time upon impact with the array posts. The in vitro microfluidic approach demonstrated here may serve as a monitoring tool to personalise and maximise the outcome of a therapeutic treatment.

Anomalous Long-Ranged Influence of an Inclusion in Momentum-Conserving Active Fluids

We show that an inclusion placed inside a dilute Stokesian suspension of microswimmers induces power-law number-density modulations and flows. These take a different form depending on whether the inclusion is held fixed by an external force—for example, an optical tweezer—or if it is free. When the inclusion is held in place, the far-field fluid flow is a Stokeslet, while the microswimmer density decays as 1/r2+ε, with r the distance from the inclusion and ε an anomalous exponent which depends on the symmetry of the inclusion and varies continuously as a function of a dimensionless number characterizing the relative amplitudes of the convective and diffusive effects. The angular dependence takes a nontrivial form which depends on the same dimensionless number. When the inclusion is free to move, the far-field fluid flow is a stresslet, and the microswimmer density decays as 1/r2 with a simple angular dependence. These long-range modulations mediate long-range interactions between inclusions that we characterize. Published by the American Physical Society 2024

Tailoring enzyme modified starch through high-moisture extrusion: Unraveling structure-property relationships

The conventional method of producing enzyme modified starch (ES) involves using a dilute starch suspension, which necessitated substantial heat expenditure for drying the product and resulted in significant energy waste. Improved extrusion cooking technology (IECT) could extrude materials under high moisture, and it is a new physical modification technology. A comprehensive investigation was undertaken to produce enzyme modified starch (ES) with varying dextrose equivalent (DE) using IECT. The study systematically explored the effects of different screw speeds on the resulting IECT-ES properties. And principal component analysis and artificial neural network analysis were used to explore the correlation between data and the synaptic weight relationships. Results indicated that Maltese cross of IECT-ES disappeared, the R1047/1022 and R995/1022 decreased to 0.7305 and 0.9012 respectively, the full width at half maximum at 480 cm−1 increased to 26.81 %, light transmittance increased to 21.20 % and DE value rose to 14.19 % at 300 rpm. The properties evolution of IECT-ES was closely related to changes in structure and conditions. This study aimed to evaluate the effectiveness of IECT as an alternative method for producing enzyme modified starch. The findings provide valuable insights into more energy-efficient starch modification processes.

Two-Dimensional MXene Flakes with Large Second Harmonic Generation and Unique Surface Responses.

MXenes are a family of two-dimensional (2D) materials with broad and varied applications in biology, materials science, photonics, and environmental remediation owing to their layered structure and high surface area-to-volume ratio. MXenes have exhibited significant nonlinear optical characteristics, which have been primarily explored in the context of photonics applications, yet the second-harmonic generation (SHG) behavior of MXenes remains an unexplored aspect of their optical properties. Herein, we demonstrate and quantify large second-order responses of 2D Ti3C2Tx MXenes both in aqueous solutions and on a silicon substrate for the first time. MXene flakes showed strong second-harmonic scattering (SHS) in a dilute suspension with a sensitivity of less than 0.1 μg/mL. Angle-dependent SHS experiments further found that the second-order responses originate from coherent 2D dipole radiation. Through confocal and atomic force microscopies, we found that the intense SHG signal from free-standing MXene flakes increases exponentially with decreasing thickness, while two-photon fluorescence increases linearly with thickness. The second-order susceptibility of the MXenes was determined to be 3.6 pm V-1 with a thickness of 10 nm, almost twice of that for an often-used SHG crystal, beta barium borate. We further explored surface properties of the MXene sheets by investigating the SHS responses upon addition of organic dye molecules to the system. It was found that the adsorption of crystal violet (CV) obeys a Langmuir adsorption model while the addition of malachite green (MG) resulted in almost no change in SHG intensity, even though the adsorption capacities for both CV (61.3 ± 1.7 mg/g) and MG (54.8 ± 2.8 mg/g) are similar. Such a stark difference in adsorption characteristics between cationic organic CV and MG dyes is likely due to their distinct orientational orderings on the MXene surfaces. This work opens many possibilities for the further employment of the family of 2D materials in photonics, optics, and surface catalysis applications.

An evaluation of the analytical and biological robustness of a method for quantifying iron in individual red blood cells via single-cell tandem ICP-mass spectrometry

This work evaluated the analytical and biological robustness of iron (Fe) determination in individual red blood cells (RBCs) via single-cell ICP-MS (SC-ICP-MS). RBCs were separated from other whole blood constituents using Ficoll-Paque™ density gradient centrifugation. While fixation with paraformaldehyde (PFA) led to RBC lysis, the use of glutaraldehyde (GA) left the RBCs intact and permitted storage of RBC suspensions in ultra-pure water for up to 16 months at 4 °C. GA-fixation also rendered the RBCs sufficiently robust to maintain integrity during their introduction into the ICP by means of nebulization of dilute suspensions. Obtaining quantitative data on a cell-per-cell basis required determination of the transport efficiency using the particle size method and external calibration against aqueous Fe standard solutions. The average Fe content per RBC obtained using SC-ICP-MS (103 fg/cell) agreed well with the value obtained using solution-based ICP-MS obtained after cell pellet digestion and with values obtained from literature. Variation of the cell number density in the suspensions analyzed between 1.5 × 105 and 6.0 × 105 cells per mL did not affect the result. Identical results from one-week interval blood drawings from healthy individuals demonstrate biological consistency. Compared to bulk analysis, the SC-ICP-MS approach offers the added value of providing information on the cell-to-cell heterogeneity in Fe content.

Hybrid particle-phase field model and renormalized surface tension in dilute suspensions of nanoparticles.

We present a two-phase field model and a hybrid particle-phase field model to simulate dilute colloidal sedimentation and flotation near a liquid-gas interface (or fluid-fluid interface in general). Both models are coupled to the incompressible Stokes equation, which is solved numerically using a combination of sine and regular Fourier transforms to account for the no-slip boundary conditions at the boundaries. The continuum two-phase field model allows us to analytically solve the equilibrium interfacial profile using a perturbative approach, demonstrating excellent agreement with numerical simulations. Notably, we show that strong coupling to particle dynamics can significantly alter the liquid-gas interface, thereby modifying the liquid-gas interfacial tension. In particular, we show that the renormalized surface tension is monotonically decreasing with increasing colloidal particle concentration and decreasing buoyant mass.

Prediction of Transport and Deposition of Porous Particles in the Respiratory System Using Eulerian-Lagrangian Approach.

Deep lung delivery is crucial for respiratory disease treatment. Although nano and submicron particles exhibited a good deposition efficiency in deep regions of the lung, powder nonuniformity and particle agglomeration reduce their efficiency. Inhalation of porous particles (PPs) can overcome the mentioned challenges due to their larger size and low-density. The present study numerically investigates the deposition and penetration efficiency of orally inhaled PPs. A revised drag coefficient was implemented for PP transport. A realistic mouth-throat to the fifth generation of the lung was reconstructed from CT-scan images. A dilute suspension of uniformly distributed particles was considered at three inhalation flow rates (15, 30, and 45 L/min). Governing equations of the flow field and particle transport are solved using an Eulerian-Lagrangian approach. The results demonstrate that inhaling PPs significantly reduces the total and regional deposition of particles. There was also a critical porosity value under moderate and high inhalation flow rates for large PPs. Below this critical value, PP deposition efficiency substantially decreases. Additionally, it was also found that under low inhalation flow rates, the impact of porosity value is negligible. Almost 95% of the PPs penetrate the lower branches. These findings provide particle engineers and pharmaceutics with profound insight into developing novel inhalation techniques and drug delivery methods for deep lung delivery.

Enhanced Diffusion and Non-Gaussian Displacements of Colloids in Quasi-2D Suspensions of Motile Bacteria

In the real world, active agents interact with surrounding passive objects, thus introducing additional degrees of complexity. The relative contributions of far-field hydrodynamic and near-field contact interactions to the anomalous diffusion of passive particles in suspensions of active swimmers remain a subject of ongoing debate. We constructed a quasi-two-dimensional microswimmer–colloid mixed system by taking advantage of Serratia marcescens’ tendency to become trapped at the air–water interface to investigate the origins of the enhanced diffusion and non-Gaussianity of the displacement distributions of passive colloidal tracers. Our findings reveal that the diffusion behavior of colloidal particles exhibits a strong dependence on bacterial density. At moderate densities, the collective dynamics of bacteria dominate the diffusion of tracer particles. In dilute bacterial suspensions, although there are multiple dynamic types present, near-field contact interactions such as collisions play a major role in the enhancement of colloidal transport and the emergence of non-Gaussian displacement distributions characterized by heavy exponential tails in short times. Despite the distinct types of microorganisms and their diverse self-propulsion mechanisms, a generality in the diffusion behavior of passive colloids and their underlying dynamics is observed.

Energy costs of the sponge pump

Sponges pump large volumes of water to obtain nutrition from a dilute suspension of food particles and dissolved organic matter. Our review of the literature shows pumping power ranging from 0.85% up to 28% of the total metabolic respiration rate Rtot, where the pumping power (Pp = pumping rate × pump pressure) is the energy per unit time transferred to the water flow. Some published values of pumping costs are high while others are low due to incorrect estimates of some pressure losses in the aquiferous system of sponges. However, as improved descriptions of detailed structures and function, especially of the filter apparatus of the choanocytes, have become available, our corrected values for pressure losses show the respiration-specific pumping costs of Pp/Rtot = 1.3 to 22%. The stated pumping cost is based on the minimal (theoretical) pump power delivered to the water, which excludes a mechanical efficiency of pump action and a metabolic efficiency that expresses the fraction of total metabolism that is devoted to driving the pump action, considerations that should be included in future studies.

Detailed comparison of recent—and dissimilar—effective-medium models incorporating dependent scattering

We compare the predictions of two recently derived effective-medium models for the effective refractive index of a turbid suspension of particles. The two formulas are notoriously dissimilar; both are based on the quasi-crystalline approximation, but the approximations used beyond this point are entirely different. Nevertheless, for dilute suspensions both reduce to the well-established van de Hulst formula. The dissimilarities between the formulas are evident for dense suspensions, where dependent-scattering effects are important. When they might coincide is, therefore, not clear. The purpose of this work is to explore the range of particle parameters and volume fractions for which both models are applicable. Our results show that, rather surprisingly, the models produce very similar curves of the real and imaginary parts of the effective refractive index for volume fractions up to 0.4 and for particles comparable to, and larger than, the wavelength, as well as for a fairly large range of refractive-index contrasts between the particles and the surrounding medium. These results significantly increase our confidence in the validity of both models.

Ideal Collision Rate of Symmetric Hexapods in a Simple Shear Flow.

An ideal collision rate (ICR) is defined as the average rate of contact between two particles that translate and rotate with the imposed fluid flow in the absence of interparticle interactions. ICR in a simple shear flow provides an estimate of the collision rate in a flowing dilute particle suspension, and its value is known only for traditional convex shapes such as spheres and cylinders. In this work, we compute the ICR for a family of symmetric hexapods (shown in Figure 1) that are particles composed of three orthogonal cylinders with coinciding centers (forming six arms) with at least two cylinders having the same length. We employed the finite-element method to obtain the rotational dynamics of hexapods and Monte Carlo simulations to calculate their ICR. Our results indicate that the ICR for hexapods is not directly proportional to the volume of the particle, in contrast to what is observed for convex shapes such as spheres and cylinders. For hexapods of the same size, the ICR could be similar for particles with order-of-magnitude differences in their volumes and it could also vary by an order of magnitude for particles with similar volumes. Specifically for hexapods termed branched fibers, which have one longer cylinder, the ICR changes by an order of magnitude even when the shorter cylinders are significantly smaller than the longer cylinder. This change is attributed to the increase in the tumbling frequency of the particle due to hydrodynamic forces acting on the shorter arms, in addition to the increased probability of collision afforded by them. Using asymptotic theory for high-aspect ratio cylinders, we showed that the tumbling period of hexapods was proportional to an algebraic power of the ratio of its arm lengths and has a weaker logarithm relationship with the aspect ratio of its arms. The collision cross-section provided the relative cross-streamline particle separations for the most likely binary collisions in the suspensions, and its value was sensitive to the arm lengths, as well. The collision rate of hexapods also could not be estimated within an order of magnitude from existing geometric models, such as the ICR of a circumscribing convex shape, ICR for one of the individual cylinders of the hexapod, or using the particle volume times the shear rate. Our work indicates that collision rate for non-convex particles in shear flows critically depends on the shape of the particle due to nontrivial changes in the particle's orientational dynamics, and the ICR calculation serves as a more reliable method for estimating their true collision rates in the suspension.

Examination of flow birefringence induced by the shear components along the optical axis using a parallel-plate-type rheometer

In this study, the concept of rheo-optics was applied to explore the flow birefringence induced by the stress components along the camera’s optical axis since it is often overlooked in the traditional theories of photoelastic flow measurement. A novel aspect of this research is that it involved conducting polarization measurements on simple shear flows, specifically from a perspective in which a shear-velocity gradient exists along the camera’s optical axis. A parallel-plate-type rheometer and a polarization camera are employed for these systematic measurements. The experimental findings for dilute aqueous cellulose nanocrystal suspensions demonstrate that the flow birefringence can be expressed as a power law based on the power of the second invariant of the deformation-rate tensor. This suggests that flow birefringence can be universally characterized by the coordinate-independent invariants and a pre-factor determined by the direction of polarization measurement. By adjusting the nonlinear term in the stress-optic law, its applicability could be expanded to include three-dimensional fluid stress fields in which the stress is distributed along the camera’s optical axis.

Instabilities and particle-induced patterns in co-rotating suspension Taylor–Couette flow

The first experimental results on pattern transitions in the co-rotation regime (i.e. the rotation ratio $\varOmega = \omega _o/\omega _i > 0$ , where $\omega _i$ and $\omega _o$ are the angular speeds of the inner and outer cylinders, respectively) of the Taylor–Couette flow (TCF) are reported for a neutrally buoyant suspension of non-colloidal particles, up to a particle volume fraction of $\phi = 0.3$ . While the stationary Taylor vortex flow (TVF) is the primary bifurcating state in dilute suspensions ( $\phi \leq ~0.05$ ), the non-axisymmetric oscillatory states, such as the spiral vortex flow (SVF) and the ribbon (RIB), appear as primary bifurcations with increasing particle loading, with an overall de-stabilization of the primary bifurcating states (TVF/SVF/RIB) being found with increasing $\phi$ for all $\varOmega \geq ~0$ . At small co-rotations ( $\varOmega \sim 0$ ), the particles play the dual role of stabilization ( $\phi < 0.1$ ) and destabilization ( $\phi \geq ~0.1$ ) on the secondary/tertiary oscillatory states. The distinctive features of the ‘particle-induced’ spiral vortices are identified and contrasted with those of the ‘fluid-induced’ spirals that operate in the counter-rotation regime.

Collective gradient sensing by dilute swimming bacteria without clustering

We characterize the taxis enhancement of swimming bacteria by collective migration without apparent clustering. We confine a dilute suspension in a shallow channel and evaluate the thermotaxis response to local heating and diffusion. By combining cell tracking analysis and numerical simulation based on simple modeling, we show that the alignment interaction suppresses orientation fluctuation, strengthens migration bias, and also prevents the dispersion of accumulated population. The results show a prominent example of how a collective motion of active matter implements a biological function. Published by the American Physical Society 2024

Dynamic behavior and mechanistic analysis in filtration of dilute particulate suspensions through nanofibrous membrane

Nanofibrous membranes have garnered significant interest in membrane filtration due to their high porosity, small pore sizes, and complex interconnected pore structures. Although many studies focus on size exclusion in surface filtration, the filtration behavior and mechanisms for particles smaller than the pore size have been relatively underexplored. In this study, membrane filtration using nanofibrous membranes was employed to elucidate the filtration behavior and mechanisms for polystyrene latex suspensions with particles smaller than the pore size. The results indicated that reducing fiber diameter decreases the pore size and enhances particle rejection, with particle retention occurring within the intricate three-dimensional nanofiber pore network. In contrast, the effect of the membrane’s basis weight was less significant compared to fiber diameter. This observation can be attributed to the concept of effective thickness; once this thickness is exceeded, particles are prevented from penetrating deeper into the membrane. The filtration model was analyzed using Hermans-Bredée blockage filtration theory and cake filtration theory, considering the mass of trapped particles within the membranes. The model characterizes the filtration process as initially following standard blocking, progressing to cake filtration once a sufficient number of pores become clogged.

Origin of Basaltic Subplinian Eruption at Shishaldin Volcano (Alaska): A Vigorously Degassing Magma Reservoir Hosting Small Bubbles

AbstractThe 1999 basaltic eruption of Shishaldin volcano (Alaska) displayed a transition between Subplinian and Strombolian activity. Strombolian bubbles indicate the presence of a periodically unstable foam at the top of magma reservoir. In contrast, a long foam, whose rupture led to the eruptive column, was also able to collect in the conduit. Laboratory experiments show that long foams can be produced in a conduit by the spreading of a stable foam accumulated at the top of the reservoir. The existence of a Taylor bubble at the onset of the Subplinian phase, also reproduced by my laboratory experiments, suggests that the foam in the reservoir was just at the transition between stable and unstable. This constrains the bubble diameter prior to the Subplinian phase to be 0.034–0.038 mm when using the foam dimensionless analysis and the underlying gas flux (0.52–0.80 m3/s). The increase in bubble diameter and potentially gas flux prior to the Strombolian activity, 0.81–1.4 m3/s, is sufficient to explain the foam in transition to be unstable. The radius of the magma reservoir is small, 200–210 m, as expected. The bubble diameter is the smallest of those estimated for classical basaltic eruptions (Etna, Kı̄lauea, Erta 'Ale), while the gas flux is among the largest. A dilute suspension of small and isolated bubbles cannot explain the large gas flux at Shishaldin. This implies numerous bubbles with a gas volume fraction ≥0.63−2%, a regime for which the bubbles form bubble clusters. The diameter of these bubble clusters, 3.0–5.4 mm, is sufficient to explain large gas fluxes.

Two regimes of dilute turbulent settling suspensions under shear

When turbulent flow is laden with negatively buoyant particles, their mean distribution in the direction of gravity can induce stable density gradients that penalize turbulent fluctuations. This effect is studied numerically for shear-driven flow with dilute noninertial sediment. The turbulent dynamics and sediment transport depend critically on particle settling velocity vs, splitting into two regimes: homogeneous weakly stratified turbulence and flow with developed turbulence atop an intermittent boundary layer. At intermediate vs, neither state can be sustained and the flow laminarizes. Published by the American Physical Society 2024

High PZT-power actuation for the separation of polydisperse glycine crystals at high flow rates and solid content using acoustophoresis

Acoustofluidics is a technique that can be used for particle separation purposes in both chemical and biological fields, with the technique appreciated for its capability of handling delicate material with no need for functionalization of particles/cells. Due to the generation of heat during the actuation of a ceramic piezoelectric element (PZT) and due to clogging, acoustofluidic separation is currently limited to dilute suspensions at a low flow rate. PZT actuation at a frequency that induces an acoustic standing wave drives the migration of particles toward the center of the channel. To overcome heating due to actuation at high power, a Peltier-based cooling system maintains the temperature of the PZT to a set point, keeping stable performance for up to a tested time of 24 min. Separation of suspensions of 10 wt.% and 5 wt.% of polydisperse glycine crystals in 5 v% methanol – ethanol buffer was achieved at a flow rate of 0.2 L/h. The cut-off found for the separation was 29.24±1.87μm and 19.51±1.19μm for 5 wt.% and 10 wt.%, respectively.

Porous structures impact on particle dynamics of non-Brownian and noncolloidal suspensions

This study aims to provide valuable insights into the impact of porous structures on particle dynamics in non-Brownian, non-colloidal suspension flows at very low Reynolds numbers. Two experimental approaches, Particle Image Velocimetry (PIV) with refractive index matching and Optical Flow Tracking Velocimetry (OFTV) were employed to analyze very dilute suspensions over various porous media models. The study considered three different porous structures with permeabilities ranging from 0.7 to 0.9 and three different thicknesses ranging from 0.2 cm to 0.5 cm, while the suspension bulk volume fraction was maintained at 3 %. In the PIV analysis, we observed that decreasing the porous permeability resulted in the maximum velocity location within the free flow region moving closer towards the interface between the flow and the porous media. We further quantified the effect of the porous structure on the suspension by characterizing interface properties, such as dimensionless slip velocity, shear rate, and slip length. These interface properties were found to be influenced by both the thickness and permeability of the porous media. Next, we analyzed particle migration due to the presence of porous structures using OFTV for very dilute suspensions of 1 %, 2 %, and 3 %, considering a porous medium with known physical properties and thickness. The study revealed two local concentration maxima: one within the free flow region on top of the rod arrays used to create the porous structure and a second along the rods' centerline inside the porous media model.

Sedimentation dynamics of passive particles in dilute bacterial suspensions: emergence of bioconvection

Microorganisms are ubiquitous in nature and technology. They inhabit diverse environments, ranging from small river tributaries and lakes, to oceans, as well as wastewater treatment plants and food manufacturing. In many of these environments, microorganisms coexist with settling particles. Here, we investigate the effects of microbial activity (swimming E. coli) on the settling dynamics of passive colloidal particles using particle tracking methods. Our results reveal the existence of two distinct regimes in the correlation length scale ( $L_u$ ) and the effective diffusivity of the colloidal particles ( $D_{eff}$ ), with increasing bacterial concentration ( $\phi _b$ ). At low $\phi _b$ , the parameters $L_u$ and $D_{eff}$ increase monotonically with increasing $\phi _b$ . Beyond critical $\phi _b$ , a second regime is found where both $D_{eff}$ and $L_u$ are independent of $\phi _b$ . We demonstrate that the transition between these regimes is characterized by the emergence of bioconvection. We use experimentally measured particle-scale quantities $L_u$ and $D_{eff}$ to predict the critical bacterial concentration for the diffusion–bioconvection transition.