Microstructural diversity, nucleation paths, and phase behavior in binary mixtures of charged colloidal spheres

We study low-salt, binary aqueous suspensions of charged colloidal spheres of size ratio Γ = 0.57, number densities below the eutectic number density n_E, and number fractions of p = 1.00-0.40. The typical phase obtained by solidification from a homogeneous shear-melt is a substitutional alloy with a body centered cubic structure. In strictly gas-tight vials, the polycrystalline solid is stable against melting and further phase transformation for extended times. For comparison, we also prepare the same samples by slow, mechanically undisturbed deionization in commercial slit cells. These cells feature a complex but well reproducible sequence of global and local gradients in salt concentration, number density, and composition as induced by successive deionization, phoretic transport, and differential settling of the components, respectively. Moreover, they provide an extended bottom surface suitable for heterogeneous nucleation of the β-phase. We give a detailed qualitative characterization of the crystallization processes using imaging and optical microscopy. By contrast to the bulk samples, the initial alloy formation is not volume-filling, and we now observe also α- and β-phases with low solubility of the odd component. In addition to the initial homogeneous nucleation route, the interplay of gradients opens various further crystallization and transformation pathways leading to a great diversity of microstructures. Upon a subsequent increase in salt concentration, the crystals melt again. Wall-based, pebble-shaped β-phase crystals and facetted α-crystals melt last. Our observations suggest that the substitutional alloys formed in bulk experiments by homogeneous nucleation and subsequent growth are mechanically stable in the absence of solid-fluid interfaces but thermodynamically metastable.

Smart "Thrombus": Self-Localizing UCST-Type Microcage

Embolization is often used to block blood supply for controlling the growth of fibroids and malignant tumors, but limited by embolic agents lacking spontaneous targeting and post-treatment removal. So we first adopted nonionic poly(acrylamide-co-acrylonitrile) with an upper critical solution temperature (UCST) to build up self-localizing microcages by inverse emulsification. The results showed that these UCST-type microcages behaved with the appropriate phase-transition threshold value around 40 °C, and spontaneously underwent an expansion-fusion-fission cycle under the stimulus of mild temperature hyperthermia. Given the simultaneous local release of cargoes, this simple but smart microcage is expected to act as a multifunctional embolic agent for tumorous starving therapy, tumor chemotherapy, and imaging.

Unidirectional drying of a suspension of diffusiophoretic colloids under gravity

Recent experiments (K. Inoue and S. Inasawa, RSC Adv., 2020, 10, 15763-15768) and simulations (J.-B. Salmon and F. Doumenc, Phys. Rev. Fluids, 2020, 5, 024201) demonstrated the significant impact of gravity on unidirectional drying of a colloidal suspension. However, under gravity, the role of colloid transport induced by an electrolyte concentration gradient, a mechanism known as diffusiophoresis, is unexplored to date. In this work, we employ direct numerical simulations and develop a macrotransport theory to analyze the advective-diffusive transport of an electrolyte-colloid suspension in a unidirectional drying cell under the influence of gravity and diffusiophoresis. We report three key findings. First, drying a suspension of solute-attracted diffusiophoretic colloids causes the strongest phase separation and generates the thinnest colloidal layer compared to non-diffusiophoretic or solute-repelled colloids. Second, when colloids are strongly solute-repelled, diffusiophoresis prevents the formation of colloid concentration gradient and hence gravity has a negligible effect on colloidal layer formation. Third, our macrotransport theory predicts new scalings for the growth of the colloidal layer. The scalings match with direct numerical simulations and indicate that the colloidal layer produced by solute-repelled diffusiophoretic colloids could be an order of magnitude thicker compared to non-diffusiophoretic or solute-attracted colloids. Our results enable tailoring the separation of colloid-electrolyte suspensions by tuning the interactions between the solvent, electrolyte, and colloids under Earth's or microgravity, which is central to ground-based and in-space applications.

Development of a novel experimental technique for the measurement of residual wall layer thickness in water-oil displacement flows

The effective removal and displacement of fluids is important in many industrial and environmental applications, such as for operation and cleaning of process equipment, fluid injection in porous media for oil recovery or aquifer remediation, or for achieving subsurface zonal isolation in new or abandoned wells. The accurate measurement of the residual fluid wall film left behind after displacement by a cleaning fluid is a long-standing challenge, particularly so for very thin fluid films where the thickness can be of the order of micrometer. We focus on the characterization of oil films left on the wall of a horizontal pipe after the pipe has been displaced by water, and develop a novel, non-intrusive analytical technique that allows the use of relevant pipe materials. The oil that originally occupies the pipe is stained by a hydrophobic dye Nile red, and an intermediate organic solvent is used to collect the residual oil volume that remains after displacing the pipe with a known volume of water. Finally, ultraviolet-visible spectroscopy is used to measure the Nile red concentration in the collected fluid, which is proportional to the residual volume of oil in the pipe. We demonstrate the methodology by conducting experiments where the displacing fluid is injected at two different imposed velocities, and where the injected fluid volume is varied. As expected, we find a gradual thinning of the oil film with increasing injected fluid volume. We compare the measured film thicknesses to a displacement model based on the steady velocity profile in a pipe, and find that experiments consistently produce smaller film thicknesses. This developed technique allows quantification of displacement and cleaning mechanisms involved in immiscible displacements at laminar, transitional and turbulent regimes, for different non-Newtonian fluid pairs, and for different realistic pipe materials and surface roughnesses.

Bioconvection pattern of Euglena under periodical illumination

Microorganisms respond to environmental conditions and often spontaneously form highly ordered convection patterns. This mechanism has been well-studied from the viewpoint of self-organization. However, environmental conditions in nature are usually dynamic. Naturally, biological systems respond to temporal changes in environmental condition. To elucidate the response mechanisms in such a dynamic environment, we observed the bioconvection pattern of Euglena under periodical changes in illumination. It is known that Euglena shows localized bioconvection patterns under constant homogeneous illumination from the bottom. Periodical changes in light intensity induced two different types of spatiotemporal patterns: alternation of formation and decomposition over a long period and complicated transition of pattern over a short period. Our observations suggest that pattern formation in a periodically changing environment is of fundamental importance to the behavior of biological systems.

Drop splashing after impact onto immiscible pools of different viscosities

Droplet impact onto liquid pools is a canonical scenario relevant to numerous natural phenomena and industrial processes. However, despite their ubiquity, multi-fluid systems with the drop and pool consisting of different liquids are far less well understood. Our hypothesis is that the post-impact dynamics greatly depends on the pool-to-droplet viscosity ratioμ_p/μ_d, which we explore over a range of six orders of magnitude using a combination of experiments and theoretical approaches (mathematical modelling and direct numerical simulation). Our findings indicate that in this scenario the splashing threshold and the composition of the ejecta sheet are controlled by the viscosity ratio. We uncover that increasing the pool viscosity decreases the splashing threshold for high viscosity pools (μ_p/μ_d≳35) when the splash comes from the droplet. By contrast, for low viscosity pools, the splash sheet comes from the pool and increasing the pool viscosity increases the splashing threshold. Surprisingly, there are conditions for which no splashing is observed under the conditions attainable in our laboratory. Furthermore, considering the interface velocity together with asymptotic arguments underlying the generation of the ejecta has allowed us to understand meaningful variations in the pressure during impact and rationalise the observed changes in the splashing threshold.

Scale dependence in hydrodynamic regime for jumping on water

Momentum transfer from the water surface is strongly related to the dynamical scale and morphology of jumping animals. Here, we investigate the scale-dependent momentum transfer of various jumping organisms and engineered systems at an air-water interface. A simplified analytical model for calculating the maximum momentum transfer identifies an intermediate dynamical scale region highly disadvantageous for jumping on water. The Weber number of the systems should be designed far from 1 to achieve high jumping performance on water. We design a relatively large water-jumping robot in the drag-dominant scale range, having a high Weber number, for maximum jumping height and distance. The jumping robot, around 10 times larger than water striders, has a take-off speed of 3.6 m/s facilitated by drag-based propulsion, which is the highest value reported thus far. The scale-dependent hydrodynamics of water jumpers provides a useful framework for understanding nature and robotic system interacting with the water surface.

An Explicit-Correction-Force Scheme of IB-LBM Based on Interpolated Particle Distribution Function

In order to obtain a better numerical simulation method for fluid-structure interaction (FSI), the IB-LBM combining the lattice Boltzmann method (LBM) and immersed boundary method (IBM) has been studied more than a decade. For this purpose, an explicit correction force scheme of IB-LBM was proposed in this paper. Different from the current IB-LBMs, this paper introduced the particle distribution function to the interpolation process from the fluid grids to the immersed boundary at the mesoscopic level and directly applied the LBM force models to obtain the interface force with a simple form and explicit process. Then, in order to ensure the mass conservation in the local area of the interface, this paper corrected the obtained interface force with the correction matrix, forming the total explicit-correction-force (ECP) scheme of IB-LBM. The results of four numerical tests were used to verify the order of accuracy and effectiveness of the present method. The streamline penetration is limited and the numerical simulation with certain application significance is successful for complex boundary conditions such as the movable rigid bodies (free oscillation of the flapping foil) and flexible deformable bodies (free deformation of cylinders). In summary, we obtained a simple and alternative simulation method that can achieve good simulation results for engineering reference models with complex boundary problems.

Simulation of Multi-Phase Flow in Autoclaves Using a Coupled CFD-DPM Approach

In this work, a numerical simulation study on the mixing characteristics of multiphase flow in an autoclave was carried out using CFD technology. The Eulerian–Eulerian model and discrete phase model (DPM) were employed to investigate the solid holdup, critical suspension speed, nonuniformity of solid suspension, gas holdup distribution, bubble tracks, and residence time during stirring leaching in the autoclave. Experiments validate the accuracy of the numerical model, and the experimental values correspond well with the simulation results. The numerical simulation results show that the solid–liquid mixing is mainly affected by the axial flow, the best agitation speed is 400 rpm, and increasing the speed further cannot make the mixture more homogenous and buildup occurred above the autoclave. The calculated critical suspension speed is 406 rpm, which is slightly lower than that obtained from the empirical formula. The gas phase is mainly concentrated in the vortex area above the blade. When the gas phase is in a completely dispersed state (N = 300 rpm), the average residence time of the bubbles is 5.66 s.

New types of complex motion of a simple camphor boat

We discuss the motion of a rectangular camphor boat, considering the position of a camphor pill in relation to the boat's stern as the control parameter. The boat moves because the pill releases surface active molecules that decrease the surface tension and support the motion. We introduce a new experimental system in which the boat rotates on a long arm around the axis located at the centre of a Petri dish; thus, the motion is restricted to a circle and can be studied under stationary conditions for a long time. The experiments confirmed two previously reported modes of motion: continuous motion when the pill was located at the boat edge and pulsating (intermittent) motion if it was close to the boat centre (Suematsu et al., J. Phys. Chem. C, 2010, 114(21), 9876-9882). For intermediate pill locations, we observed a new, unreported type of motion characterised by oscillating speed (i.e. oscillating motion). Different modes of motion can be observed for the same pill location. The experimental results are qualitatively confirmed using a simple reaction-diffusion model of the boat evolution used in the above-mentioned paper.

A Nodal Immersed Finite Element-Finite Difference Method

The immersed finite element-finite difference (IFED) method is a computational approach to modeling interactions between a fluid and an immersed structure. The IFED method uses a finite element (FE) method to approximate the stresses, forces, and structural deformations on a structural mesh and a finite difference (FD) method to approximate the momentum and enforce incompressibility of the entire fluid-structure system on a Cartesian grid. The fundamental approach used by this method follows the immersed boundary framework for modeling fluid-structure interaction (FSI), in which a force spreading operator prolongs structural forces to a Cartesian grid, and a velocity interpolation operator restricts a velocity field defined on that grid back onto the structural mesh. With an FE structural mechanics framework, force spreading first requires that the force itself be projected onto the finite element space. Similarly, velocity interpolation requires projecting velocity data onto the FE basis functions. Consequently, evaluating either coupling operator requires solving a matrix equation at every time step. Mass lumping, in which the projection matrices are replaced by diagonal approximations, has the potential to accelerate this method considerably. This paper provides both numerical and computational analyses of the effects of this replacement for evaluating the force projection and for the IFED coupling operators. Constructing the coupling operators also requires determining the locations on the structure mesh where the forces and velocities are sampled. Here we show that sampling the forces and velocities at the nodes of the structural mesh is equivalent to using lumped mass matrices in the IFED coupling operators. A key theoretical result of our analysis is that if both of these approaches are used together, the IFED method permits the use of lumped mass matrices derived from nodal quadrature rules for any standard interpolatory element. This is different from standard FE methods, which require specialized treatments to accommodate mass lumping with higher-order shape functions. Our theoretical results are confirmed by numerical benchmarks, including standard solid mechanics tests and examination of a dynamic model of a bioprosthetic heart valve.

Modifying last layer in polyelectrolyte multilayer coatings for capillary electrophoresis of proteins

Protein adsorption on the inner wall of the fused silica capillary wall is an important concern for capillary electrophoresis (CE) analysis since it is mainly responsible for separation efficiency reduction. Successive Multiple Ionic-polymer Layers (SMIL) are used as capillary coatings to limit protein adsorption, but even low residual adsorption strongly impacts the separation efficiency, especially at high separation voltages. In this work, the influence of the chemical nature and the PEGylation of the polyelectrolyte deposited in the last layer of the SMIL coating was investigated on the separation performances of a mixture of four model intact proteins (myoglobin (Myo), trypsin inhibitor (TI), ribonuclease a (RNAse A) and lysozyme (Lyz)). Poly(allylamine hydrochloride) (PAH), polyethyleneimine (PEI), ε-poly(L-lysine) (εPLL) and α-poly(L-lysine) (αPLL) were compared before and after chemical modification with polyethyleneglycol (PEG) of different chain lengths. The experimental results obtained by performing electrophoretic separations at different separation voltages allowed determining the residual retention factor of the proteins onto the capillary wall via the determination of the plate height at different solute velocities and demonstrated a strong impact of the polycationic last layer on the electroosmotic mobility, the separation efficiency and the overall resolution. Properties of SMIL coatings were also characterized by quartz microbalance and atomic force microscopy, demonstrating a glassy structure of the films.

Acousto-dielectric tweezers for size-insensitive manipulation and biophysical characterization of single cells

The intrinsic biophysical properties of cells, such as mechanical, acoustic, and electrical properties, are valuable indicators of a cell's function and state. However, traditional single-cell biophysical characterization methods are hindered by limited measurable properties, time-consuming procedures, and complex system setups. This study presents acousto-dielectric tweezers that leverage the balance between controllable acoustophoretic and dielectrophoretic forces applied on cells through surface acoustic waves and alternating current electric fields, respectively. Particularly, the balanced acoustophoretic and dielectrophoretic forces can trap cells at equilibrium positions independent of the cell size to differentiate between various cell-intrinsic mechanical, acoustic, and electrical properties. Experimental results show our mechanism has the potential for applications in single-cell analysis, size-insensitive cell separation, and cell phenotyping, which are all primarily based on cells' intrinsic biophysical properties. Our results also show the measured equilibrium position of a cell can inversely determine multiple biophysical properties, including membrane capacitance, cytoplasm conductivity, and acoustic contrast factor. With these features, our acousto-dielectric tweezing mechanism is a valuable addition to the resources available for biophysical property-based biological and medical research.

The effects of aneurysmal subarachnoid hemorrhage on cerebral vessel diameter and flow velocity

BACKGROUND

Transcranial Doppler flow velocity is used to monitor for cerebral vasospasm following aneurysmal subarachnoid hemorrhage. Generally, blood flow velocities appear inversely related to the square of vessel diameter representing local fluid dynamics. However, studies of flow velocity-diameter relationships are few, and may identify vessels for which diameter changes are better correlated with Doppler velocity. We therefore studied a large retrospective cohort with concurrent transcranial Doppler velocities and angiographic vessel diameters.

METHODS

This is a single-site, retrospective, cohort study of adult patients with aneurysmal subarachnoid hemorrhage, approved by the UT Southwestern Medical Center Institutional Review Board. Study inclusion required transcranial Doppler measurements within </= 24 hours of vessel imaging. Vessels assessed were: bilateral anterior, middle, posterior cerebral arteries; internal carotid siphons; vertebral arteries; and basilar artery. Flow velocity-diameter relationships were constructed and fitted with a simple inverse power function. A greater influence of local fluid dynamics is suggested as power factors approach two.

RESULTS

98 patients were included. Velocity-diameter relationships are curvilinear, and well fit by a simple inverse power function. Middle cerebral arteries showed the highest power factors (>1.1, R²>0.9). Furthermore, velocity and diameter changed (P<0.033) consistent with the signature time course of cerebral vasospasm.

CONCLUSIONS

These results suggest that middle cerebral artery velocity-diameter relationships are most influenced by local fluid dynamics, which supports these vessels as preferred endpoints in Doppler detection of cerebral vasospasm. Other vessels showed less influence of local fluid dynamics, pointing to greater role of factors outside the local vessel segment in determining flow velocity.

AI-based analysis of 3D position and orientation of red blood cells using a digital in-line holographic microscopy

The morphological and mechanical characteristics of red blood cells (RBCs) largely vary depending on the occurrence of hematologic disorders. Variations in the rheological properties of RBCs affect the dynamic motions of RBCs, especially their rotational behavior. However, conventional techniques for measuring the orientation of biconcave-shaped RBCs still have some technical limitations, including complicated optical setups, complex post data processing, and low throughput. In this study, we propose a novel image-based technique for measuring 3D position and orientation of normal RBCs using digital in-line holographic microscopy (DIHM) and artificial intelligence (AI). Formaldehyde-fixed RBCs are immobilized in coagulated polydimethylsiloxane (PDMS). Holographic images of RBCs positioned at various out-of-plane angles are acquired by precisely manipulating the PDMS-trapped RBC sample attached to a 4-axis optical stage. With the aid of deep learning algorithms for data augmentation and regression analysis, the out-of-plane angle of RBCs is directly predicted from the captured holographic images. The 3D position and in-plane angle of RBCs are acquired by employing numerical reconstruction and ellipse detection methods. Combining these digital image processing techniques, the 3D positional and orientational information of each RBC recorded in a single holographic image is measured within 23.5 and 3.07 s, respectively. The proposed AI-based DIHM technique that can extract the 3D position, orientation, and morphology of individual RBCs would be utilized to analyze the dynamic translational and rotational motions of abnormal RBCs with hematologic disorders in shear flows through further research.

Rheological Behavior of DNP/HMX Melt-Cast Explosives with Bimodal and Trimodal Particle-Size Distributions

As a matrix for melt-cast explosives, 3,4-dinitropyrazole (DNP) is a promising alternative to 2,4,6-trinitrotoluene (TNT). However, the viscosity of molten DNP is considerably greater compared with that of TNT, thus, requiring the viscosity of DNP-based melt-cast explosive suspensions to be minimized. In this paper, the apparent viscosity of a DNP/HMX (cyclotetramethylenetetranitramine) melt-cast explosive suspension is measured using a Haake Mars III rheometer. Both bimodal and trimodal particle-size distributions are used to minimize the viscosity of this explosive suspension. First, the optimal diameter ratio and mass ratio (two crucial process parameters) between coarse and fine particles are obtained from the bimodal particle-size distribution. Second, based on the optimal diameter ratio and mass ratio, trimodal particle-size distributions are used to further minimize the apparent viscosity of the DNP/HMX melt-cast explosive suspension. Finally, for either the bimodal or trimodal particle-size distribution, if the original data between the apparent viscosity and solid content are normalized, the resultant plot of the relative viscosity versus reduced solid content collapses to a single curve, and the effect of the shear rate on this curve is further investigated.

Mechanism of Surface Wettability of Nanostructure Morphology Enhancing Boiling Heat Transfer: Molecular Dynamics Simulation

In this paper, the interaction mechanism between the solid–liquid–gas interface phenomenon caused by nanostructure and surface wettability and boiling heat transfer is described, and the heat transfer theory of single wettable nanostructure surface and mixed wettable nanostructure surface is proposed. Through molecular dynamics simulation, the thermodynamic model of the wettable surface of nanostructures is established. The nanostructures are set as four rectangular lattice structures with a height of 18 Å. The solid atoms are platinum atoms, and the liquid atoms are argon atoms. The simulation results show that with the increase of surface hydrophilicity of nanostructures, the fluid temperature increases significantly, and the heat transfer at the interface is enhanced. With the increase in surface hydrophobicity of nanostructures, the atoms staying on the surface of nanostructures are affected by the hydrophobicity, showing a phenomenon of exclusion, and the evaporation rate in the evaporation area of nanostructures is significantly increased. In addition, the mixed wettable surface is influenced by the atomic potential energy and kinetic energy of the solid surface, and when compared with the pure wettable surface under the nanostructure, it changes the diffusion behavior of argon atoms on the nanostructure surface, enhances the heat transfer phenomenon compared with the pure hydrophobic surface, and enhances the evaporation phenomenon compared with the pure hydrophilic surface. This study provides insights into the relationship between the vapor film and the heating surface with mixed wettability and nanostructures.

The power of super-resolution microscopy in modern biomedical science

Super-resolution microscopy (SRM) technology that breaks the diffraction limit has revolutionized the field of cell biology since its appearance, which enables researchers to visualize cellular structures with nanometric resolution, multiple colors and single-molecule sensitivity. With the flourishing development of hardware and the availability of novel fluorescent probes, the impact of SRM has already gone beyond cell biology and extended to nanomedicine, material science and nanotechnology, and remarkably boosted important breakthroughs in these fields. In this review, we will mainly highlight the power of SRM in modern biomedical science, discussing how these SRM techniques revolutionize the way we understand cell structures, biomaterials assembly and how assembled biomaterials interact with cellular organelles, and finally their promotion to the clinical pre-diagnosis. Moreover, we also provide an outlook on the current technical challenges and future improvement direction of SRM. We hope this review can provide useful information, inspire new ideas and propel the development both from the perspective of SRM techniques and from the perspective of SRM's applications.

Machine learning in additive manufacturing & Microfluidics for smarter and safer drug delivery systems

A new technological passage has emerged in the pharmaceutical field, concerning the management, application, and transfer of knowledge from humans to machines, as well as the implementation of advanced manufacturing and product optimisation processes. Machine Learning (ML) methods have been introduced to Additive Manufacturing (AM) and Microfluidics (MFs) to predict and generate learning patterns for precise fabrication of tailor-made pharmaceutical treatments. Moreover, regarding the diversity and complexity of personalised medicine, ML has been part of quality by design strategy, targeting towards the development of safe and effective drug delivery systems. The utilisation of different and novel ML techniques along with Internet of Things sensors in AM and MFs, have shown promising aspects regarding the development of well-defined automated procedures towards the production of sustainable and quality-based therapeutic systems. Thus, the effective data utilisation, prospects on a flexible and broader production of "on demand" treatments. In this study, a thorough overview has been achieved, concerning scientific achievements of the past decade, which aims to trigger the research interest on incorporating different types of ML in AM and MFs, as essential techniques for the enhancement of quality standards of customised medicinal applications, as well as the reduction of variability potency, throughout a pharmaceutical process.

Modular microfluidics for life sciences

The advancement of microfluidics has enabled numerous discoveries and technologies in life sciences. However, due to the lack of industry standards and configurability, the design and fabrication of microfluidic devices require highly skilled technicians. The diversity of microfluidic devices discourages biologists and chemists from applying this technique in their laboratories. Modular microfluidics, which integrates the standardized microfluidic modules into a whole, complex platform, brings the capability of configurability to conventional microfluidics. The exciting features, including portability, on-site deployability, and high customization motivate us to review the state-of-the-art modular microfluidics and discuss future perspectives. In this review, we first introduce the working mechanisms of the basic microfluidic modules and evaluate their feasibility as modular microfluidic components. Next, we explain the connection approaches among these microfluidic modules, and summarize the advantages of modular microfluidics over integrated microfluidics in biological applications. Finally, we discuss the challenge and future perspectives of modular microfluidics.

Effects of surface acoustic waves on droplet impact dynamics

HYPOTHESIS

Surface acoustic waves (SAW) propagating along a solid surface can significantly affect the dynamics of droplet impact. Although droplet impact in presence of SAW has been attempted recently, here, we investigate the effects of surface wettability, droplet size, impact velocity, and SAW power on the impact and spreading dynamics along with post-impact oscillation dynamics of a drop.

EXPERIMENTS

Here, we study droplet impact on a surface exposed to traveling SAW produced using an interdigitated electrode patterned on a piezoelectric substrate. The effects of Weber number (We), surface wettability, and SAW power on the impact and spreading dynamics and post-impact oscillation dynamics are studied.

FINDINGS

Our study unravels that the interplay between capillary and viscous forces, and inertia forces arising due to pre-impact kinetic energy and SAW-induced bulk acoustic streaming underpins the phenomena. Remarkably, we find that the effect of SAW on droplet impact dynamics is predominant in the case of a hydrophilic (HPL) substrate at a higher SAW power and smaller We and hydrophobic (HPB) substrate irrespective of SAW power. Our study reveals that the maximum droplet spreading diameter increases with SAW power at smaller We for an HPL surface whereas it is independent of SAW power at higher We. Post-impact oscillation of a droplet over an HPL surface is found to be overdamped with a smaller amplitude compared to an HPB substrate, and a faster decay in oscillation amplitude is observed in the case of an HPB surface and higher We. Our study provides an improved understanding of droplet impact on a surface exposed to SAW that may find relevance in various practical applications.

Diffusion of gold nanoparticles in porous silica monoliths determined by dynamic light scattering

HYPOTHESIS

The applicability of the dynamic light scattering method for the determination of particle diffusivity under confinement without applying refractive index matching was not adequately explored so far. The confinement effect on particle diffusion in a porous material which is relevant for particle chromatography has also not yet been fully characterized.

EXPERIMENTS

Dynamic light scattering experiments were performed for unimodal dispersions of 11-mercaptoundecanoic acid-capped gold nanoparticles. Diffusion coefficients of gold nanoparticles in porous silica monoliths were determined without limiting refractive index matching fluids. Comparative experiments were also performed with the same nanoparticles and porous silica monolith but applying refractive index matching.

FINDINGS

Two distinct diffusivities could be determined inside the porous silica monolith, both smaller than that in free media, showing a slowing-down of the diffusion processes of nanoparticles under confinement. While the larger diffusivity can be related to the slightly slowed-down diffusion of particles in the bulk of the pores and in the necks connecting individual pores, the smaller diffusivity might be related to the diffusion of particles near the pore walls. It shows that the dynamic light scattering method with a heterodyne detection scheme can be used as a reliable and competitive tool for determining particle diffusion under confinement.

A Novel Trans-Tracheostomal Retrograde Inhalation Technique Increases Subglottic Drug Deposition Compared to Traditional Trans-Oral Inhalation

Subglottic stenosis represents a challenging clinical condition in otolaryngology. Although patients often experience improvement following endoscopic surgery, recurrence rates remain high. Pursuing measures to maintain surgical results and prevent recurrence is thus necessary. Steroids therapy is considered effective in preventing restenosis. Currently, however, the ability of trans-oral steroid inhalation to reach and affect the stenotic subglottic area in a tracheotomized patient is largely negligible. In the present study, we describe a novel trans-tracheostomal retrograde inhalation technique to increase corticosteroid deposition in the subglottic area. We detail our preliminary clinical outcomes in four patients treated with trans-tracheostomal corticosteroid inhalation via a metered dose inhaler (MDI) following surgery. Concurrently, we leverage computational fluid-particle dynamics (CFPD) simulations in an extra-thoracic 3D airway model to gain insight on possible advantages of such a technique over traditional trans-oral inhalation in augmenting aerosol deposition in the stenotic subglottic region. Our numerical simulations show that for an arbitrary inhaled dose (aerosols spanning 1–12 µm), the deposition (mass) fraction in the subglottis is over 30 times higher in the retrograde trans-tracheostomal technique compared to the trans-oral inhalation technique (3.63% vs. 0.11%). Importantly, while a major portion of inhaled aerosols (66.43%) in the trans-oral inhalation maneuver are transported distally past the trachea, the vast majority of aerosols (85.10%) exit through the mouth during trans-tracheostomal inhalation, thereby avoiding undesired deposition in the broader lungs. Overall, the proposed trans-tracheostomal retrograde inhalation technique increases aerosol deposition rates in the subglottis with minor lower-airway deposition compared to the trans-oral inhalation technique. This novel technique could play an important role in preventing restenosis of the subglottis.

Efficient sonocatalytic degradation of orange II dye and real textile wastewater using peroxymonosulfate activated with a novel heterogeneous TiO2–FeZn bimetallic nanocatalyst

TiO₂–FeZn nanocatalyst combined with sonolysis were used to activate peroxymonosulfate (PMS) as a highly efficient advanced oxidation process (US/TiO₂–FeZn/PMS) for the decoloration of orange II dye (OII) and real textile wastewater. The characterization of the as-synthesized NPs was performed by SEM, FTIR, EDX and XRD analyses. Optimal experimental conditions of operational parameters were obtained: pH = 3, 15 mg/L initial OII concentration, 0.2 g/L PMS, 0.7 g/L nanocatalyst dosing, and 300 W ultrasonic power. The decolorization was observed to increase with increasing the dose of nanocatalyst and the ultrasonic power, and with decreasing pH (under acidic conditions). Under optimal experimental conditions, decolorization and COD removal of textile wastewater were 99.9% and 74.6%, respectively, at 40 min. The TiO₂–FeZn/PMS/US as a novel process exhibited a higher removal of OII (95%) than TiO₂ NPs/PMS/US process (54%). The OII removal efficiency by the different processes decreased in the following order: TiO₂–FeZn/US/PMS > TiO₂–FeZn/PMS > TiO₂–FeZn/US > TiO₂ /US/PMS > US/PMS > TiO₂–FeZn > PMS > US. The recyclability study revealed that the process could be reused up to three consecutive cycles. The current US/nanocatalyst/PMS system was concluded to be an efficient, reusable and stable nanocatalyst for the oxidation of textile dyes.

How can machine learning and multiscale modeling benefit ocular drug development?

The eyes possess sophisticated physiological structures, diverse disease targets, limited drug delivery space, distinctive barriers, and complicated biomechanical processes, requiring a more in-depth understanding of the interactions between drug delivery systems and biological systems for ocular formulation development. However, the tiny size of the eyes makes sampling difficult and invasive studies costly and ethically constrained. Developing ocular formulations following conventional trial-and-error formulation and manufacturing process screening procedures is inefficient. Along with the popularity of computational pharmaceutics, non-invasive in silico modeling & simulation offer new opportunities for the paradigm shift of ocular formulation development. The current work first systematically reviews the theoretical underpinnings, advanced applications, and unique advantages of data-driven machine learning and multiscale simulation approaches represented by molecular simulation, mathematical modeling, and pharmacokinetic (PK)/pharmacodynamic (PD) modeling for ocular drug development. Following this, a new computer-driven framework for rational pharmaceutical formulation design is proposed, inspired by the potential of in silico explorations in understanding drug delivery details and facilitating drug formulation design. Lastly, to promote the paradigm shift, integrated in silico methodologies were highlighted, and discussions on data challenges, model practicality, personalized modeling, regulatory science, interdisciplinary collaboration, and talent training were conducted in detail with a view to achieving more efficient objective-oriented pharmaceutical formulation design.

In silico modeling of patient-specific blood rheology in type 2 diabetes mellitus

Increased blood viscosity in type 2 diabetes mellitus (T2DM) is a risk factor for the development of insulin resistance and diabetes-related vascular complications; however, individuals with T2DM exhibit heterogeneous hemorheological properties, including cell deformation and aggregation. Using a multiscale red blood cell (RBC) model with key parameters derived from patient-specific data, we present a computational study of the rheological properties of blood from individual patients with T2DM. Specifically, one key model parameter, which determines the shear stiffness of the RBC membrane (μ) is informed by the high-shear-rate blood viscosity of patients with T2DM. At the same time, the other, which contributes to the strength of the RBC aggregation interaction (D₀), is derived from the low-shear-rate blood viscosity of patients with T2DM. The T2DM RBC suspensions are simulated at different shear rates, and the predicted blood viscosity is compared with clinical laboratory-measured data. The results show that the blood viscosity obtained from clinical laboratories and computational simulations are in agreement at both low and high shear rates. These quantitative simulation results demonstrate that the patient-specific model has truly learned the rheological behavior of T2DM blood by unifying the mechanical and aggregation factors of the RBCs, which provides an effective way to extract quantitative predictions of the rheological properties of the blood of individual patients with T2DM.

Effect of valve spacing on peristaltic pumping

Peristaltic fluid pumping due to a periodically propagating contraction wave in a vessel fitted with one-way elastic valves is investigated numerically. It is concluded that the valve spacing within the vessel relative to the contraction wavelength plays a critical role in providing efficient pumping. When the valve spacing does not match the wavelength, the valves open asynchronously and the volume of the vessel segments bounded by two consecutive valves changes periodically, thereby inducing volumetric fluid pumping. The volumetric pumping leads to higher pumping flowrate and efficiency against an adverse pressure gradient. The optimum pumping occurs when the ratio of valve spacing to contraction wavelength is about2/3. This pumping regime is characterized by a longer period during which the valves are open. The results are useful for further understanding the pumping features of lymphatic system and provide insight into the design of biomimetic pumping devices.

Multicellular structures in thin free liquid films induced by thermocapillary effect

HYPOTHESIS

Multicellular convective structures that are induced in a fluid exposed to temperature difference are commonly observed in nature and in daily life. Different types of basic flow patterns are induced in a free liquid film by thermocapillary effect, whereas the formation of such multicellular structures has not been hitherto unravelled.

EXPERIMENTS

A thin film of high-Prandtl-number liquid is prepared in a rectangular aperture of the order of 0.1 mm in thickness sustained by its surface tension. A designated temperature difference is imposed between the end surfaces of the aperture to generate a thermocapillary-driven convection in the free liquid film. We monitor the induced thermal flow patterns to evaluate the cell numbers and their wavelength by experimental and numerical approaches.

FINDINGS

The multicellular structure is established by the thermocapillary effect in the free liquid films. The cell number increases in a stepwise manner as the liquid-film width increases. When the cell number increases, another pair of the cells always newly emerges. We determine the wavelength in a non-dimensional manner, and present the variation of the wavelength against the aspect ratio corresponding to the liquid-film width. The results are compared to those of convectional Marangoni-Bénard convection.

Seagrass deformation affects fluid instability and tracer exchange in canopy flow

Monami is the synchronous waving of a submerged seagrass bed in response to unidirectional fluid flow. Here we develop a multiphase model for the dynamical instabilities and flow-driven collective motions of buoyant, deformable seagrass. We show that the impedance to flow due to the seagrass results in an unstable velocity shear layer at the canopy interface, leading to a periodic array of vortices that propagate downstream. Our simplified model, configured for unidirectional flow in a channel, provides a better understanding of the interaction between these vortices and the seagrass bed. Each passing vortex locally weakens the along-stream velocity at the canopy top, reducing the drag and allowing the deformed grass to straighten up just beneath it. This causes the grass to oscillate periodically even in the absence of water waves. Crucially, the maximal grass deflection is out of phase with the vortices. A phase diagram for the onset of instability shows its dependence on the fluid Reynolds number and an effective buoyancy parameter. Less buoyant grass is more easily deformed by the flow and forms a weaker shear layer, with smaller vortices and less material exchange across the canopy top. While higher Reynolds number leads to stronger vortices and larger waving amplitudes of the seagrass, waving amplitude is maximized at intermediate grass buoyancy. All together, our theory and computations develop an updated schematic of the instability mechanism consistent with experimental observations.

Tides regulate the flow and density of Antarctic Bottom Water from the western Ross Sea

Antarctic Bottom Water (AABW) stores heat and gases over decades to centuries after contact with the atmosphere during formation on the Antarctic shelf and subsequent flow into the global deep ocean. Dense water from the western Ross Sea, a primary source of AABW, shows changes in water properties and volume over the last few decades. Here we show, using multiple years of moored observations, that the density and speed of the outflow are consistent with a release from the Drygalski Trough controlled by the density in Terra Nova Bay (the "accelerator") and the tidal mixing (the "brake"). We suggest tides create two peaks in density and flow each year at the equinoxes and could cause changes of ~ 30% in the flow and density over the 18.6-year lunar nodal tide. Based on our dynamic model, we find tides can explain much of the decadal variability in the outflow with longer-term changes likely driven by the density in Terra Nova Bay.

Overflow Control for Sustainable Development by Superwetting Surface with Biomimetic Structure

Liquid flowing around a solid edge, i.e., overflow, is a commonly observed flow behavior. Recent research into surface wetting properties and microstructure-controlled overflow behavior has attracted much attention. Achieving controllable macroscale liquid dynamics by manipulating the micro-nanoscale liquid overflow has stimulated diverse scientific interest and fostered widespread use in practical applications. In this review, we outline the evolution of overflow and present a critical survey of the mechanism of surface wetting properties and microstructure-controlled liquid overflow in multilength scales ranging from centimeter to micro and even nanoscale. We summarize the latest progress in utilizing the mechanisms to manipulate liquid overflow and achieve macroscale liquid dynamics and in emerging applications to manipulate overflow for sustainable development in various fields, along with challenges and perspectives.

EHD stability of a cylindrical boundary separating double Reiner-Rivlin fluids

The major aim of this work is to achieve a mathematical technique to scrutinize the nonlinear instability of a vertical cylindrical boundary separation of two streaming Reiner-Rivlin liquids. The system is portrayed by an unchanged longitudinal electric strength. Furthermore, the action of mass and heat transfer (MHT) and permeable media are also considered. The problem is not only of methodological interest but also of scientific and practical interest. To shorten the mathematical analysis, Hsieh's modulation together with the viscous potential theory (VPT) is employed. The nonlinear diagram is contingent on tackling the governing linear mechanism along with the nonlinear applicable border restrictions. A non-dimensional process produces several non-dimensional physical numbers. A linear dispersion equation is attained and the stability standards are theoretically governed and numerically established. The nonlinear stability procedure reveals a Ginzburg-Landau formula. Consequently, nonlinear stability stipulations are accomplished. Furthermore, by way of the Homotopy perturbation approach, along with the expanded frequency concept, an accurate perturbed technique of surface deflection is attained theoretically and numerically. To validate the theoretical outcomes, the analytical expression is confirmed through the Rung-Kutta of the fourth order. The stable and unstable zones are signified graphically displaying the influences of several non-dimensional numbers.

Imaging Cerebral Arteries Tortuosity and Velocities by Transcranial Doppler Ultrasound Is a Reliable Assessment of Brain Aneurysm in Mouse Models

Background

During the past few decades, several pathophysiological processes contributing to intracranial aneurysm (IA) rupture have been identified, including irregular IA shape, altered hemodynamic stress within the IA, and vessel wall inflammation. The use of preclinical models of IA and imaging tools is paramount to better understand the underlying disease mechanisms.

Methods

We used 2 established mouse models of IA, and we analyzed the progression of the IA by magnetic resonance imaging, transcranial Doppler, and histology.

Results

In both models of IA, we observed, by transcranial Doppler, a significant decrease of the blood velocities and wall shear stress of the internal carotid arteries. We also observed the formation of tortuous arteries in both models that were correlated with the presence of an aneurysm as confirmed by magnetic resonance imaging and histology. A high grade of tortuosity is associated with a significant decrease of the mean blood flow velocities and a greater artery dilation.

Conclusions

Transcranial Doppler is a robust and convenient imaging method to evaluate the progression of IA. Detection of decreased blood flow velocities and increased tortuosity can be used as reliable indicators of IA.

Generalization of Young-Laplace, Kelvin, and Gibbs-Thomson equations for arbitrarily curved surfaces

The Young-Laplace, Kelvin, and Gibbs-Thomson equations form a cornerstone of colloidal and surface sciences and have found successful applications in many subfields of physics, chemistry, and biology. The Gibbs-Thomson effect, for example, predicts that small crystals are in equilibrium with their liquid melt at a lower temperature than large crystals and the positive interfacial energy increases the energy required to form small particles with a high curvature interface. In cases of liquids contained within porous media (confined geometry), the effect indicates decreasing the freezing/melting temperatures and the increment of the temperature is inversely proportional to the pore size. These phenomena can be reformulated for Gaussian maps of macromolecules and can be asked the following question: can one use the equations for predicting the melting temperature and shape of polymer chains in confined geometries? The answer is no, mainly because macromolecules form highly curved surfaces (Gaussian maps), and the equations hold only for simple geometries (sphere, plane, or cylinder). Here, we present general Young-Laplace, Kelvin, and Gibbs-Thomson equations for arbitrarily curved surfaces and apply them to predict temperature distribution on a few protein surfaces. Also, after increased interest toward liquid/liquid phase separation in biology, we derive generic Ostwald ripening and show that for shape-changing condensates, instead of a monotonic growing mechanism, a variety of processes are possible. Due to the generality of equations, we clarify that at appropriate internal/external pressure conditions systems, bounded by surfaces, may adopt any shape and thermal stability is strongly influenced by the geometries of confined spaces.

Innovative Methods of Centrifugal Separation

Considered are: (i) separation of gaseous molecules of different weight by the Ultra Centrifuge with application to large scale uranium enrichment to fuel nuclear power plants; (ii) separation of micron-sized particulate matter from fluids in mechanical devices in which use is made of inertial and centrifugal forces; (iii) separation of gaseous mixtures by fast expansion and cooling such that one of the gaseous components forms a mist of micron-sized droplets which are separated by centrifugation; and (iv) separation of components from gases by absorbing liquid films in small sized rotating channels. For each of these technologies we consider: physics of the separation process, assessment of the influence of fluid flow in the separation device, identification of leading parameters and their effect on design, experimental evidence, status of development, areas of application, position compared to other technologies, and economic value.

The Post-Shock Nonequilibrium Relaxation in a Hypersonic Plasma Flow Involving Reflection off a Thermal Discontinuity

The evolution in the post-shock nonequilibrium relaxation in a hypersonic plasma flow was investigated during a shock’s reflection off a thermal discontinuity. It was found that within a transitional period, the relaxation zone parameters past both the reflected and transmitted waves evolve differently compared to that in the incident wave. In a numerical example for the non-dissociating N2 gas heated to 5000 K/10,000 K across the interface and M = 3.5, the relaxation time determined for the transmitted wave is up to 50% shorter and the relaxation depth for both waves is significantly reduced, thus resulting in a weakened wave structure. The results of the extension into larger values of heating strength and the shock Mach numbers are discussed. The findings can be useful in the areas of research involving strong shocks interacting with optical discharges or other heated media on the scale where the shock structure becomes important.

Numerical Study of Particle Separation through Integrated Multi-Stage Surface Acoustic Waves and Modulated Driving Signals

The manipulation of biomedical particles, such as separating circulating tumor cells from blood, based on standing surface acoustic wave (SSAW) has been widely used due to its advantages of label-free approaches and good biocompatibility. However, most of the existing SSAW-based separation technologies are dedicated to isolate bioparticles in only two different sizes. It is still challenging to fractionate various particles in more than two different sizes with high efficiency and accuracy. In this work, to tackle the problems of low efficiency for multiple cell particle separation, integrated multi-stage SSAW devices with different wavelengths driven by modulated signals were designed and studied. A three-dimensional microfluidic device model was proposed and analyzed using the finite element method (FEM). In addition, the effect of the slanted angle, acoustic pressure, and the resonant frequency of the SAW device on the particle separation were systemically studied. From the theoretical results, the separation efficiency of three different size particles based on the multi-stage SSAW devices reached 99%, which was significantly improved compared with conventional single-stage SSAW devices.

Classification of Blood Rheological Models through an Idealized Symmetrical Bifurcation

The assumed rheological behavior of blood influences the hemodynamic characteristics of numerical blood flow simulations. Until now, alternative rheological specifications have been utilized, with uncertain implications for the results obtained. This work aims to group sixteen blood rheological models in homogeneous clusters, by exploiting data generated from numerical simulations on an idealized symmetrical arterial bifurcation. Blood flow is assumed to be pulsatile and is simulated using a commercial finite volume solver. An appropriate mesh convergence study is performed, and all results are collected at three different time instants throughout the cardiac cycle: at peak systole, early diastole, and late diastole. Six hemodynamic variables are computed: the time average wall shear stress, oscillatory shear index, relative residence time, global and local non-Newtonian importance factor, and non-Newtonian effect factor. The resulting data are analyzed using hierarchical agglomerative clustering algorithms, which constitute typical unsupervised classification methods. Interestingly, the rheological models can be partitioned into three homogeneous groups, whereas three specifications appear as outliers which do not belong in any partition. Our findings suggest that models which are defined in a similar manner from a mathematical perspective may behave substantially differently in terms of the data they produce. On the other hand, models characterized by different mathematical formulations may belong to the same statistical group (cluster) and can thus be considered interchangeably.

Wavelet image scattering based glaucoma detection

BACKGROUND

The ever-growing need for cheap, simple, fast, and accurate healthcare solutions spurred a lot of research activities which are aimed at the reliable deployment of artificial intelligence in the medical fields. However, this has proved to be a daunting task especially when looking to make automated diagnoses using biomedical image data. Biomedical image data have complex patterns which human experts find very hard to comprehend. Against this backdrop, we applied a representation or feature learning algorithm: Invariant Scattering Convolution Network or Wavelet scattering Network to retinal fundus images and studied the the efficacy of the automatically extracted features therefrom for glaucoma diagnosis/detection. The influence of wavelet scattering network parameter settings as well as 2-D channel image type on the detection correctness is also examined. Our work is a distinct departure from the usual method where wavelet transform is applied to pre-processed retinal fundus images and handcrafted features are extracted from the decomposition results. Here, the RIM-ONE DL image dataset was fed into a wavelet scattering network developed in the Matlab environment to achieve a stage-wise decomposition process called wavelet scattering of the retinal fundus images thereby, automatically learning features from the images. These features were then used to build simple and computationally cheap classification algorithms.

RESULTS

Maximum detection correctness of 98% was achieved on the held-out test set. Detection correctness is highly sensitive to scattering network parameter setting and 2-D channel image type.

CONCLUSION

A superficial comparison of the classification results obtained from our work and those obtained using a convolutional neural network underscores the potentiality of the proposed method for glaucoma detection.

A computational framework for investigating bacteria transport in microvasculature

Blood-borne bacteria disseminate in tissue through microvasculature or capillaries. Capillary size, presence of red blood cells (RBCs), and bacteria motility affect bacteria intracapillary transport, an important yet largely unexplored phenomenon. Computational description of the system comprising interactions between plasma, RBCs, and motile bacteria in 5-10 μm diameter capillaries pose several challenges. The Immersed Boundary Method (IBM) was used to resolve the capillary, deformed RBCs, and bacteria. The challenge of disparate coupled time scales of flow and bacteria motion are reconciled by a temporal multiscale simulation method. Bacterium-wall and bacterium-RBC collisions were detected using a hierarchical contact- detection algorithm. Motile bacteria showed a net outward radial velocity of 2.8 µm/s compared to -0.5 µm/s inward for non-motile bacteria; thus, exhibiting a greater propensity to escape the bolus flow region between RBCs and marginate for potential extravasation, suggesting motility enhances extravasation of bacteria from capillaries.