Particle Migration Research Articles

Ultra-Stretchable Microfluidic Devices for Optimizing Particle Manipulation in Viscoelastic Fluids.

Viscoelastic microfluidics leverages the unique properties of non-Newtonian fluids to manipulate and separate micro- or submicron particles. Channel geometry and dimension are crucial for device performance. Traditional rigid microfluidic devices require numerous iterations of fabrication and testing to optimize these parameters, which is time-consuming and costly. In this work, we developed a flexible microfluidic device using ultra-stretchable and biocompatible Flexdym material to overcome this issue. Our device allows for simultaneous modification of channel dimensions by external stretching. We fabricated a stretchable device with an initial square microchannel (30 μm × 30 μm), and the channel aspect ratio can be adjusted from 1 to 5 by external stretching. Next, the effects of aspect ratio, particle size, flow rate, and poly(ethylene oxide) (PEO) concentration that make the fluid viscoelastic on particle migration were investigated. Finally, we demonstrated the feasibility of our approach by testing channels with an aspect ratio of 3 for the separation of both particles and cells.

Surface Self-Diffusion Induced Sintering of Nanoparticles.

Despite the critical role of sintering phenomena in constraining the long-term durability of nanosized particles, a clear understanding of nanoparticle sintering has remained elusive due to the challenges in atomically tracking the neck initiation and discerning different mechanisms. Through the integration of in situ transmission electron microscopy and atomistic modeling, this study uncovers the atomic dynamics governing the neck initiation of Pt-Fe nanoparticles via a surface self-diffusion process, allowing for coalescence without significant particle movement. Real-time imaging reveals that thermally activated surface morphology changes in individual nanoparticles induce significant surface self-diffusion. The kinetic entrapment of self-diffusing atoms in the gaps between closely spaced nanoparticles leads to the nucleation and growth of atomic layers for neck formation. This surface self-diffusion-driven sintering process is activated at a relatively lower temperature compared to the classic Ostwald ripening and particle migration and coalescence processes. The fundamental insights have practical implications for manipulating the morphology, size distribution, and stability of nanostructures by leveraging surface self-diffusion processes.

Soil fluidisation induced by fine particles migration: Insights from the Shenzhen 2015 landfill landslide

Naturally completely decomposed granite (CDG) soil typically exhibits strain-hardening behavior under undrained shear conditions. Nevertheless, flow-type landslides are not uncommon in CDG landfills. This paper endeavors to address the observed contradiction by conducting a case study of the 2015 Shenzhen landslides. Based on field investigations, we propose a hypothesis for the initiation and evolution of flow-type landslides in CDG landfill slopes, termed ‘clay particle argillization, mud-water migration, and static liquefaction’. This hypothesis was verified by element-scale internal erosion tests and triaxial tests, and further elucidated by microscale particle analysis. It was observed that the internal erosion-induced removal of plastic fine particles and retention of low-plasticity fine particles from CDG soil promotes the sliding and reorganization of coarse granules under shear stress, thereby increasing the soil's susceptibility to fluidization under undrained conditions. The proposed hypothesis and experimental findings provide new insights into the instability and subsequent extensive runout of CDG landfills and analogous broadly graded landslides.

Study on the Migration Mechanism of Temporary Plugging Agent in Artificial Fractures of Deep Geothermal Reservoirs Based on the Euler-Euler Multiphase Flow Model

The successful use of temporary plugging diverting fracturing technology requires an understanding of the migration and plugging processes of temporary plugging agents into artificial fractures under high temperature settings. In this study, a multiphase flow model for the migration of temporary plugging agents in artificial fractures was developed using the Euler-Euler framework, and numerical simulations were conducted at elevated temperatures. Various factors, including plugging agent injection velocity, concentration, carrying fluid viscosity, wall temperature, and fracture width, were systematically analyzed to assess their impact on the agent’s migration behavior. Detailed analyses, using cloud diagrams of particle volume fraction, velocity, and turbulence intensity, clarified the underlying mechanisms influencing the migration process. The results indicate that as the injection velocity increases, the height of blockages near the wellbore decreases, while the blockage length initially increases before declining. Increasing the concentration of the plugging agent leads to a rise in blockage height and a shift in the front edge toward the injection point. Enhancing the viscosity of the carrying fluid enables the plugging agent to migrate deeper into the fracture, improving deep plugging effectiveness. While changes in wall temperature have limited impact on blockage morphology, temperatures exceeding the critical threshold of 573K significantly intensify particle migration. Moreover, increasing fracture width enhances both the height and length of blockages, with the optimal plugging effect observed when the plugging agent diameter is approximately one-third of the fracture width.

Mass and momentum balance during particle migration in the pressure-driven flow of frictional non-Brownian suspensions

The transient shear-induced particle migration of frictional non-Brownian suspensions is studied using particle-resolved simulations. The numerical method – the fictitious domain method – is well suited to heterogeneous flows thanks to a frame-invariant formulation of the subgrid (lubrication) corrections that does not involve the ambient flow (Orsi et al., J. Comput. Phys., vol. 474, 2023, 111823). The paper aims to give an accurate quantitative picture of the mass and momentum balance during the flow. The various assumptions and local constitutive laws that together form the suspension balance model (SBM) are thoroughly examined. To this purpose, the various quantities of interest are locally averaged in space and time, and their profile across the channel is extensively studied, with specific attention to the time evolution of the different contributions, either hydrodynamic in nature or from contact interactions, to the shear and normal stresses. The latter, together with the velocity gradient in the wall-normal direction and the volume fraction profile, yield the local constitutive laws, which are compared with their counterpart obtained in homogeneous shear flow. A fair agreement is observed except in a layering area at the boundaries and at the very centre of the channel. In addition, the main assumption of the SBM, i.e. the local relation between the hydrodynamic force on the particles and the particle flux, is meticulously investigated. The hydrodynamic force is found to be mainly a drag, except in the lower range of the probed volume fractions, where a non-drag contribution is observed.

Soil behavior near the PVD filter and interpretation of clogging during vacuum preloading

Clogging near the prefabricated vertical drain filter during vacuum preloading reduces the efficiency of ground treatment, and existing research has yet to agree on the cause. This study conducted a small model test to comprehensively examine the soil's water content and particle size distribution at different distances from the filter and tailwater for different preloading durations. The results showed that the particle size distribution of the soil at different selected points exhibited a slightly irregular difference. This indicates that no measurable migration of particles occurs, which was further confirmed by the particle size distribution of the soil after the one-dimensional consolidation test. The application of vacuum pressure did not significantly change the soil particle composition and, thus, its compressibility and permeability properties. The particle size distribution of tailwater revealed that a small number of fine particles with a mean size of 3–4 μm rapidly discharged with tailwater during the initial stage of vacuum preloading, which was much less than the mean size of the slurry and pore size of the filter. The decrease in water content and permeability indicates non-uniform consolidation, and the increase in the permeability coefficient ratio suggests the formation and development of the clogging zone.

Approximate Solution of Fuzzy Sedimentary Ocean Basin Boundary Value Problem Using Variational Iteration Method

The problem that is considered in this paper is developed from the migration of the coastline over time in sedimentary ocean basins that may be considered as a moving boundary problem with a variable latent heat transfer equation as a governing equation if one of the boundaries is unknown, with boundary and initial conditions are given. This model discusses the shoreline movement within sedimentary basins. The main objective of the paper is to use fuzzy logic to formulate this problem when modelling sand particles migration in saltwater by utilizing fuzzy numbers whenever replacing some of the problem's parameters. In reality, the proposed model is more general than the nonfuzzy or crisp models. By employing the variational iteration method, the approximate numerical solution of the proposed problem has been found, and computer software written in Mathematica 11 have been used to obtain the results.

Research on energy utilization of electron beam melting for silicon purification

As a key material for solar photovoltaic power generation, the purity of silicon plays a crucial role in the photovoltaic conversion efficiency and performance. The electron beam melting (EBM), which is irreplaceable metallurgical methods, has been widely used for silicon purification. However, high energy consumption result in relatively high cost, which limits its widespread application in melting industry. This study was centered on the energy distribution that proposes a new approach to reducing energy consumption during EBM. The scanning form of electron beam is most important factor concerning energy consumption. Different scanning forms can affect the characteristics of molten pool which is closely correlated with the purification process. In this study, efficient and energy-saving preparation processes of EBM are being explored through the numerical model and experimental. The relationship of energy distribution form and purification efficiency is studied. The research results indicate that electron beam scanning under uniform distribution can achieve more uniform temperature and flow field distribution in the molten pool, which has a higher energy utilization rate during silicon purification by EBM. The mathematical model established in this study can effectively predict the energy consumption level under different melting parameters, providing important reference for the optimization and energy conservation of EBM process for other materials. In addition, evaporation characteristics of volatile impurity, migration of the second phase particles in the melt, controlling of melt pool characteristics and the crystal growth control all can be conducted based on this study.

Dramatic droplet deformation through interfacial particles jamming

Droplets of one fluid in a second, immiscible fluid are typically spherical in shape due to the interfacial tension between the two fluids. Shear forces can lead to droplet deformation when they are subjected to flow, and these effects can be further modified when the droplet is stabilized by a surfactant due to a flow-induced gradients in the surfactant concentration. An alternative method of stabilizing droplets is through the use of colloidal particles, whose stabilization behavior is intrinsically different from molecular surfactants. Under the same flow condition, a gradient of particle concentration can result in the jamming of particles in regions with a high packing density, making the interface solid-like, albeit only under compression and not tension. However, how this asymmetry in the surfactant properties alters the droplet shape under shear is unknown. Here, we show that shear of particle-stabilized droplets can lead to a remarkable array of shape deformations as the droplets flow through a constrained microchannel. The shear-induced migration of particles on the surface results in the formation of an elastic shell at the back of the droplet, which can wrinkle and invaginate, ultimately leading to a unique core-shell structure. The shapes depend on the Peclet number of the flow, reflecting the balance of shear forces that drive the particles and diffusion that randomizes them. These findings highlight the consequences of the asymmetry in the forces between the particles and provide a unique method to controllably create droplets with a vast array of different shapes.

Aggregate migration and self-organisation properties of asphalt mixtures during compaction process

The objective of this investigation is to analyze the migration behaviours and self-organisation properties of aggregate-asphalt particle system. The three-dimensional discrete element dynamic compaction model was established, and the accuracy and reliability were validated by self-developed particle migration tracking test. The evaluation indices of aggregate particles was measured to investigate the migration behaviours of aggregates and its evolution law. The self-organisation behaviours of aggregate-asphalt particle system were explored combined with the contact behaviours. The results indicate that the migration process of aggregate particles can be delimited into three stages of compaction degree lower than 50%, from 50% to 70%, and above 70%. The larger the aggregate size, the more significant the rotation migration. Under the combined effect of migration and contact, the coarse aggregates migrate to find stable contact points to self-organize the skeleton structure. Asphalt mortar fills the voids eliminated to self-organize the dense structure.

Lateral migration of a deformable fluid particle in a square channel flow of viscoelastic fluid

Cross-stream migration of a deformable fluid particle is investigated computationally in a pressure-driven channel flow of a viscoelastic fluid via interface-resolved simulations. Flow equations are solved fully coupled with the Giesekus model equations using an Eulerian–Lagrangian method and extensive simulations are performed for a wide range of flow parameters to reveal the effects of particle deformability, fluid elasticity, shear thinning and fluid inertia on the particle migration dynamics. Migration rate of a deformable particle is found to be much higher than that of a solid particle under similar flow conditions mainly due to the free-slip condition on its surface. It is observed that the direction of particle migration can be altered by varying shear thinning of the ambient fluid. With a strong shear thinning, the particle migrates towards the wall while it migrates towards the channel centre in a purely elastic fluid without shear thinning. An onset of elastic flow instability is observed beyond a critical Weissenberg number, which in turn causes a path instability even for a nearly spherical particle. An inertial path instability is also observed once particle deformation exceeds a critical value. Shear thinning is found to be suppressing the path instability in a viscoelastic fluid with a high polymer concentration whereas it reverses its role and promotes path instability in a dilute polymer solution. It is found that migration of a deformable particle towards the wall induces a secondary flow with a velocity that is approximately an order of magnitude higher than the one induced by a solid particle under similar flow conditions.

Particle migration due to non-uniform laminar flow

Abstract Using methods of granular mechanics in the quasi-static limit, with inter-particle interactions derived from the lubrication limit, the intensity of velocity fluctuations in the slurry is associated with fluctuations in the local distribution of inter-particle distances. These are shown to consist of a vector intensity and a scalar intensity; the former couples to the first velocity gradient, the latter (which is associated with solidosity fluctuations) couples to the second velocity gradient. Rheologies for both are presented, as is the rheology that links the particle pressure to the intensity of the velocity fluctuations (also known as the ‘granular temperature’) to the dispersive pressure. The rheologies are informed by experimental results. The granular temperature profile, modified from previous work, is responsible for axial particle migration (Bagnold effect). Two broad categories are assessed: symmetrical vertical and non-symmetrical lateral flow. For the latter the roughness of the boundary walls and a non-zero density contrast are important; this case is studied for a system in which flow effects are confined to the immediate vicinity of the boundary. Sensitivity analysis reveals several key variables including the parameters that control a slipping boundary condition and the mean solidosity in the conduit. For lateral flow, a sedimentary deposit with a solidosity profile may develop near the upper or lower boundary. The theory predicts an approximate relation between the fluid-particle density contrast and sediment thickness as a function of the mean flow rate, conduit width, the mean particle diameter and fluid viscosity that has utility in a range of engineering and geological situations where particulate matter is transported in the laminar flow regime.

The impact of Microbially Induced Calcite Precipitation (MICP) on sand internal erosion resistance: A microfluidic study

Fine particle internal erosion involves the detachment, transport, and deposition of fine particles within soil, significantly impacting agriculture, engineering, and environmental protection. Microbially induced calcite precipitation (MICP) has proven to be an effective method for controlling internal erosion. To optimize MICP protocols for effective erosion control, understanding the microscopic mechanisms by which MICP reduces fine particle erosion is essential but remains unclear. In this study, microfluidic chip experiments were conducted to observe the characteristics of calcium carbonate crystals and fine particles before and after MICP reinforcement and erosion. Through quantitative analysis of the calcium carbonate produced by MICP and eroded fine particles, the effects of bacterial density, concentration of cementation solution, and flow rate on the efficiency of MICP in resisting soil internal erosion were investigated. In addition, three primary mechanisms by which MICP reduces fine particle erosion were identified. Firstly, in situ sand stabilization occurs when calcium carbonate generated by MICP bonds and encapsulates fine particles, forming larger particles that remain stable under erosive flow conditions. Secondly, regional sand stabilization is achieved as MICP-produced calcium carbonate crystals narrow or block the flow channels, preventing extensive migration of fine particles. Lastly, adjacent particle stabilization is facilitated by calcium carbonate crystals, which remain stable during water flow erosion and alter the erosion flow lines, creating zones adjacent to the crystals where fine particle movement is minimized. These findings provide a deeper understanding of the role of MICP in erosion control and can inform the development of more effective erosion mitigation strategies.

Examining the Mid to Long‐Term Variability in Saturated Hydraulic Conductivity of Sandy Soils and Its Influencing Factors Under Constant Head Test in the Laboratory

AbstractSaturated hydraulic conductivity (Ks) is a crucial parameter that influences water flow in saturated soils, with applications in various fields such as surface water runoff, soil erosion, drainage, and solute transport. However, accurate determination of Ks is challenging due to temporal and spatial uncertainties. This study addresses the knowledge gap regarding the long‐term behavior of Ks in sandy soils with less than 10% fine particles. The research investigates the changes in Ks over a long period of constant head tests and examines the factors influencing its variation. Two sandy samples were tested using a hydraulic conductivity cell, and the hydraulic head and discharge were recorded for over 50 days. The results show a general decline in Ks throughout the test, except for brief periods of increase. At the end of both tests, there are noticeable reductions in the saturated hydraulic conductivities of the samples, with one sample being 96% and the other sample 91% less than the maximum recorded saturated hydraulic conductivity during the tests. Furthermore, the relationship between flow rate and hydraulic head gradient does not follow the expected linear correlation from Darcy's law, highlighting the complex nature of sandy soil saturated hydraulic conductivity. The investigation of soil properties in three different sections of the samples before and after the tests revealed a decrease in the percentage of fine particles and a shift in specific gravity from the bottom to the top of the sample, suggesting particle migration along the flow direction. Factors such as clogging by fine particles and pore pressure variation contribute to the changes in Ks. The findings of this research show the importance of considering changes of saturated hydraulic conductivity during constant‐head laboratory tests. Therefore, this study provides evidence for the requirement to further assess the laboratory methods for measurement of the saturated hydraulic conductivity in sandy soil mixtures.

Experimental study on fluid-Induced vibration caused by backward erosion piping based on acoustic emission

Abstract Double-layer dike foundation is one of the most common stratum types with the highest probability of catastrophic failure in dike engineering, and the most important danger is backward erosion piping. The occurrence and development of piping need to be monitored by reliable technology. Therefore, the influence of different piping outlet forms on acoustic emission monitoring is explored by using a self-designed test device. The results show that the overall evolution process of ringing count and accumulated energy can be divided into two stages according to the time of piping. In the first and middle stages, the particle collision is caused by sand boiling at the outlet, and the measured value is small. In the later stage, piping channels appear, particle migration occurs, and the measured value is large. Moreover, the maximum impact value decreases with the increase of the piping outlet, and the cumulative energy growth rate of the three groups of tests in the early stage is the same, while the head difference required for the rapid growth stage is inversely proportional to the piping outlet.

Influence of concentration on thermophoresis of spherical aerosol particles within a Brinkman medium

Abstract We are examining the thermophoretic movement of a uniform mixture of spherical aerosol particles, all with the same properties, as they are situated within a porous material. These particles can have various thermal conductivity and surface characteristics. This analysis focuses on situations where the Péclet and Reynolds numbers are small. The influence of particle interactions is carefully considered by using a unit cell model, a well-established method known for its accurate predictions in the context of sedimentation for monodisperse suspensions of spherical particles. The porous medium is represented as a Brinkman fluid characterized by a Darcy permeability, which can be determined directly from experimental observations. This medium is considered to be uniform and isotropic, and the solid matrix is in thermal equilibrium with the fluid flowing through the voids of the medium. The Knudsen number is assumed to be low, enabling the description of fluid flow through the porous medium using a continuum model that includes temperature jump, thermal creep, frictional slip, and thermal stress slip at the aerosol particle’s surface. The conservation equations for energy and momentum are individually tackled within each cell. In this model, each cell represents a spherical particle enclosed by a concentric shell of surrounding fluid. The thermophoretic particle migration velocity is determined across different cases. We derive analytical expressions for this average particle velocity, expressing it in terms of the particle volume fraction. It is observed that different cell models yield somewhat varied results for particle velocity. Generally, with a fixed permeability parameter characterizing the porous medium, an increase in the thermal stress slip coefficient tends to decrease the normalized thermophoretic velocity across the different cell models. The results are in good agreement with the available data as documented in the existing literature. Additionally, a parallel examination of aerosol sphere sedimentation is provided.

A new evaluation method for sand control performance of slotted screen

Evaluation method of slotted screen performance is a critical concern for sand control. Despite its significance, comprehensive investigations into quantitative study of sand retaining during the complex flow process need further in-depth study, particularly in numerical simulation research model of the slotted screen. This paper adopts the two-way coupling calculation method of computational fluid dynamics (CFD) and discrete element method (DEM) for the particle movement and fluid flow at the slotted screen and establishes a numerical model of the sand particle migration process through the slotted screen. The sand particle deposition and plugging dynamics on the surface of the sand retaining medium and the micro sand plugging mechanism and blockage pattern are studied. The fractal dimension of formation sand is introduced to characterize different distribution of formation sand. The effects of fluid viscosity, width of the slotted screen, sand production concentration, fluid velocity and particle size distribution on sand production, sand passing rate and particle size distribution of formation sand passing through the slot and plugging behavior are studied. The sand plugging process of slotted screen can be divided into three stages: the beginning of blockage, the intensification of blockage, and the balance of blockage. Under simulated conditions, the initial sand plugging rate is 87%, and as the blockage intensifies, the sand plugging rate gradually increases to 96.5%. Based on the analysis of numerical simulation results, the index of total seepage resistance damage amplitude and velocity, general sand retention ratio and passed-sand size were defined to construct the integrated evaluation method of slotted screen media, which involves the performance of sand retention, anti-plugging and flow ability. A brand new set of ideas on quantifying the integrated evaluation method were introduced to calculate the individual and integrated performance indicators of the slotted screen based on simulated data. The quantitative evaluation of the sand control performance of the slotted screen provides guidance and reference for the parameter design of the slotted screen in the target oil region and the accuracy optimization of the sand control medium.

Damage Evaluation of Unconsolidated Sandstone Particle Migration Reservoir Based on Well–Seismic Combination

The primary factor constraining the performance of unconsolidated sandstone reservoirs is blockage from particle migration, which reduces the capacity of liquid production. By utilizing logging, seismic, core–testing, and oil–well production data, the reservoir damage induced by particle migration in the Bohai A oilfield was characterized and predicted through combined well–seismic methods. This research highlights the porosity, permeability, median grain diameter, and pore structure as the primary parameters influencing reservoir characteristics. Based on their permeability differences, reservoirs can be categorized into Type I (permeability ≥ 800 mD), Type II (400 mD < permeability < 800 mD), and Type III (permeability ≤ 400 mD). The results of the core displacement experiments revealed that, compared to their initial states, the permeability change rates for Type I and Type II reservoirs exceeded 50%, whereas the permeability change rate for Type III reservoirs surpassed 200%. Furthermore, by combining this quantitative relationship model with machine learning techniques and well–seismic methods, the distribution of permeability change rates caused by particle migration across the entire region was successfully predicted and validated against production data from three oil wells. In addition, to build a reliable deep learning model, a sensitivity analysis of the hyperparameters was conducted to determine the activation function, optimizer, learning rate, and neurons. This method enhances the prediction efficiency of reservoir permeability changes in offshore oilfields with limited coring data, providing important decision support for reservoir protection and field development.

Integrated neural network based simulation of thermo solutal radiative double-diffusive convection of ternary hybrid nanofluid flow in an inclined annulus with thermophoretic particle deposition

Numerical simulation of magnetohydrodynamic radiative double-diffusive convective heat and mass transfer of a ternary hybrid nanofluid with thermophoretic particle deposition in an inclined annulus is studied. Both sides of cylinders are preserved at uniform temperatures, while the remaining sides are thermally insulated. The coupled nonlinear partial differential equations are solved using a finite difference approach. The detailed computational results of heat and mass transfer rate, temperature, concentration and flow fields are presented and discussed for a distinct range of significant physical parameters. The results demonstrated that increasing the angle of the inclination factor increases fluid flow. In light of buoyancy, when the annular cavity is set downward, the enhanced shear velocity is transported upwards, where it reaches its optimal temperature. The higher temperature difference leads to augmented dynamic fluid motion, particularly in the radial direction, as the system seeks to balance the thermal energy through convective transfer. The Brownian motion parameter has abrupt mobility of nanoparticles in liquid, which promotes particle collisions with liquid molecules, yielding kinetic energy. Increased values of thermophoretic parameters augment the thermophoresis force. Thermophoretic particle deposition increases the concentration of nanoparticles in an annulus due to the migration of particles from hot to cold regions under the influence of temperature gradient. The heat and mass transfer characteristics of ternary hybrid nanofluid are forecasted through an artificial neural network-based Levenberg–Marquardt backpropagated algorithm. The built model shows the heat and mass transfer rate root mean squared values through the Levenberg–Marquardt algorithm as one and mean squared error values as 1e-07 and 1e-08 respectively.

The effect of MnO2 additive on the microstructure and mechanical properties of magnesium aluminate spinel

AbstractIn this study, varying amounts of MnO2 up to 5 wt.% were added to magnesium aluminate spinel (MA) bodies using a solid‐state sintering method at 1200–1600°C. The effect of MnO2 addition on the phase composition, microstructure, distribution of elements, and ionic valence of MA was investigated via X‐ray diffraction, scanning electron microscopy, energy‐dispersive spectroscopy, and X‐ray photoelectron spectroscopy, respectively. The results showed that Mg2+ ions in MA crystals were replaced by Mn2+ ions, resulting in the formation of the (Mg1‐xMnx)Al2O4 solid solution. The distorted crystal structures promoted the sintering reactions, and the mechanical characteristics of MA were greatly improved by the solid solution strengthening process. When the additive amount of MnO2 was 5 wt.% and the sintered temperature reached at 1600°C, excess manganese ions hardly dissolved into the lattice of MA. And these ions were only distributed at the grain boundaries of MgAl2O4, forming a “barrier” that hindered the migration and diffusion of particles, thereby suppressing the sintering process and weakening the mechanical strength of MA.