State Key Laboratory of Nonlinear Mechanics

We report significant back stress strengthening and strain hardening in gradient structured (GS) interstitial-free (IF) steel. Back stress is long-range stress caused by the pileup of geometrically necessary dislocations (GNDs). A simple equation and a procedure are developed to calculate back stress basing on its formation physics from the tensile unloading-reloading hysteresis loop. The gradient structure has mechanical incompatibility due to its grain size gradient. This induces strain gradient, which needs to be accommodated by GNDs. Back stress not only raises the yield strength but also significantly enhances strain hardening to increase the ductility. [GRAPHICS] .

An analysis based on the energy conservation is presented for the self-propelled droplet during coalescence of two droplets of the same size over a superhydrophobic rough surface. The self-propelled behavior occurs only for the coalescence of droplets with a certain range of radius. An analytical relation is established among the coalescence-induced velocity, surface energy, viscous dissipation, and droplet size if gravity is negligible. The coalescence-induced velocity increases with increasing droplet size to a maximum and then decreases with the size, which is in good accord with the experimental observation reported in the literature.

Inertial microfluidics has emerged as an important tool for manipulating particles and cells. For a better design of inertial microfluidic devices, we conduct 3D direct numerical simulations (DNS) and experiments to determine the complicated dependence of focusing behaviour on the particle size, channel aspect ratio, and channel Reynolds number. We find that the well-known focusing of the particles at the two centers of the long channel walls occurs at a relatively low Reynolds number, whereas additional stable equilibrium positions emerge close to the short walls with increasing Reynolds number. Based on the numerically calculated trajectories of particles, we propose a two-stage particle migration which is consistent with experimental observations. We further present a general criterion to secure good focusing of particles for high flow rates. This work thus provides physical insight into the multiplex focusing of particles in rectangular microchannels with different geometries and Reynolds numbers, and paves the way for efficiently designing inertial microfluidic devices.

Viscoelasticity-induced particle migration has recently received increasing attention due to its ability to obtain high-quality focusing over a wide range of flow rates. However, its application is limited to low throughput regime since the particles can defocus as flow rate increases. Using an engineered carrier medium with constant and low viscosity and strong elasticity, the sample flow rates are improved to be 1 order of magnitude higher than those in existing studies. Utilizing differential focusing of particles of different sizes, here, we present sheathless particle/cell separation in simple straight microchannels that possess excellent parallelizability for further throughput enhancement. The present method can be implemented over a wide range of particle/cell sizes and flow rates. We successfully separate small particles from larger particles, MCF-7 cells from red blood cells (RBCs), and Escherichia coli (E. coli) bacteria from RBCs in different straight microchannels. The proposed method could broaden the applications of viscoelastic microfluidic devices to particle/cell separation due to the enhanced sample throughput and simple channel design.

Cerebral aneurysms and brain tumors are leading life-threatening diseases worldwide. By deliberately occluding the target lesion to reduce the blood supply, embolization has been widely used clinically to treat cerebral aneurysms and brain tumors. Conventional embolization is usually performed by threading a catheter through blood vessels to the target lesion, which is often limited by the poor steerability of the catheter in complex neurovascular networks, especially in submillimeter regions. Here, we propose magnetic soft microfiberbots with high steerability, reliable maneuverability, and multimodal shape reconfigurability to perform robotic embolization in submillimeter regions via a remote, untethered, and magnetically controllable manner. Magnetic soft microfiberbots were fabricated by thermal drawing magnetic soft composite into microfibers, followed by magnetizing and molding procedures to endow a helical magnetic polarity. By controlling magnetic fields, magnetic soft microfiberbots exhibit reversible elongated/aggregated shape morphing and helical propulsion in flow conditions, allowing for controllable navigation through complex vasculature and robotic embolization in submillimeter regions. We performed in vitro embolization of aneurysm and tumor in neurovascular phantoms and in vivo embolization of a rabbit femoral artery model under real-time fluoroscopy. These studies demonstrate the potential clinical value of our work, paving the way for a robotic embolization scheme in robotic settings.

This letter addresses the issue of deformation mechanisms and mechanical tensile behavior of the twinned metal nanowires using atomistic simulations. Free surfaces are always the preferential dislocation nucleation sites in the initial inelastic deformation stage, while with further plastic deformation, twin boundary interfaces will act as sources of dislocations with the assistance of the newly formed defects. The smaller the twin boundary spacing, the higher the yielding stresses of the twinned nanowires. Twin boundaries, which serve both as obstacles to dislocation motion and dislocation sources, can lead to hardening effects and contribute to the tensile ductility. This work illustrates that the mechanical properties of metal nanowires could be controlled by tailoring internal growth twin structures.

For the first time, the enhanced recovery of confined methane (CH4) with carbon dioxide (CO2) is investigated through molecular dynamics simulations. The adsorption energy and configuration of CH4 and CO2 on the carbon surface were compared, which shows that CO2 is a good candidate in displacing confined CH4. The energy barrier required for displacing CH4 by CO2 injection was found to depend on the displacement angle. When CO2 approached vertically to the carbon surface, the displacement of CH4 occurred most easily. The curvature and size effects of the carbon nanopores on CH4 recovery were revealed and indicated that there exists an optimum pore size making the displacement occur most efficiently. The underlying mechanisms of these phenomena were uncovered. Our findings and related analyses may help to understand CO2 enhanced gas recovery from the atomic level and assist the future design in engineering.

Nanostructured (NS) metals typically fail soon after yielding, starting with the formation of narrow shear bands. Here we report the observation of shear band delocalization in gradient metals. Shear bands were nucleated and delocalized in the NS layers by propagating along the gage length soon after yielding, converting the shear band into a localized strain zone (LSZ). Synergistic work hardening was developed in the LSZ by regaining dislocation hardening capability, and by back-stress hardening from the strain gradients in the axial and depth directions, which helped with enhancing global ductility.

Spherical Janus particles are one of the most prominent examples for active Brownian objects. Here, we study the diffusiophoretic motion of such microswimmers in experiment and in theory. Three stages are found: simple Brownian motion at short times, superdiffusion at intermediate times, and finally diffusive behavior again at long times. These three regimes observed in the experiments are compared with a theoretical model for the Langevin dynamics of self-propelled particles with coupled translational and rotational motion. Besides the mean square displacement also higher displacement moments are addressed. In particular, theoretical predictions regarding the non-Gaussian behavior of self-propelled particles are verified in the experiments. Furthermore, the full displacement probability distribution is analyzed, where in agreement with Brownian dynamics simulations either an extremely broadened peak or a pronounced double-peak structure is found, depending on the experimental conditions.

Inertial microfluidics has been widely used in high-throughput manipulation of particles and cells by hydrodynamic forces, without the aid of externally applied fields. The performance of inertial microfluidic devices largely relies on precise prediction of particle trajectories that are determined by inertial lift acting on particles. The only way to accurately obtain lift forces is by direct numerical simulation (DNS); however, it is burdensome when applied to practical microchannels with complex geometries. Here, we propose a fitting formula for inertial lift on a sphere drawn from DNS data obtained in straight channels. The formula consists of four terms that represent the shear-gradient-induced lift, the wall-induced lift, the slip-shear lift, and the correction of the shear-gradient-induced lift, respectively. Notably, as a function of the parameters of a local flow field, it possesses good adaptability to complex channel geometries. This generalized formula is further implemented in the Lagrangian particle tracking method to realize fast prediction of particle trajectories in two types of widely used microchannels: a long serpentine and a double spiral microchannel, demonstrating its ability to efficiently design and optimize inertial microfluidic devices.

Cracking and fracture of electrodes under diffusion during lithiation and delithiation is one of the main factors responsible for short life span of lithium based batteries employing high capacity electrodes. Coupling effects among lithium diffusion, stress evolution and crack propagation have a significant effect on dynamic processes of electrodes during cycling. In this paper, a phase field model coupling lithium diffusion and stress evolution with crack propagation is established. Then the model is applied to a silicon thin film electrode to explore the coupling effects on diffusion and crack propagation paths. During lithiation, simulation results show that lithium accumulates at crack tips and the lithium accumulation further reduces the local hydrostatic stress. Single and multiple crack geometries are considered to elucidate some of the crack patterns in thin film electrodes as a consequence of coupling effects and crack interactions.

The snap-through and pull-in instabilities of the micromachined arch-shaped beams under an electrostatic loading are studied both theoretically and experimentally. The pull-in instability that results in a system collision with an electrode substrate may lead to a system failure and, thus, limits the system maximum displacement. The beam/plate structure with a flat initial configuration under an electrostatic loading can only experience the pull-in instability. With the different arch configurations, the structure may experience either only the pull-in instability or the snap-through and pull-in instabilities together. As shown in our computation and experiment, those arch-shaped beams with the snap-through instability have the larger maximum displacement compared with the arch-shaped beams with only the pull-in stability and those with the flat initial configuration. The snap-through occurs by exerting a fixed load, and the structure experiences a discontinuous displacement jump without consuming power. Furthermore, after the snap-through jump, the structures are demonstrated to have the capacity to withstand further electrostatic loading without pull-in. Those properties of consuming no power and increasing the structure deflection range without pull-in is very useful in microelectromechanical systems design, which can offer better sensitivity and tuning range.

Most deformation twins in nanocrystalline face-centered cubic (fcc) metals have been observed to form from grain boundaries. The growth of such twins requires the emission of Shockley partials from the grain boundary on successive slip planes. However, it is statistically improbable for a partial to exist on every slip plane. Here we propose a dislocation reaction and cross-slip mechanism on the grain boundary that would supply a partial on every successive slip plane for twin growth. This mechanism can also produce a twin with macrostrain smaller than that caused by a conventional twin.

Confined diffusion is ubiquitous in nature. Ever since the "anomalous yet Brownian" motion was observed, the non-Gaussianity in confined diffusion has been unveiled as an important issue. In this Letter, we experimentally investigate the characteristics and source of non-Gaussian behavior in confined diffusion of nanoparticles suspended in polymer solutions. A time-varied and size-dependent non-Gaussianity is reported based on the non-Gaussian parameter and displacement probability distribution, especially when the nanoparticle's size is smaller than the typical polymer mesh size. This non-Gaussianity does not vanish even at the long-time Brownian stage. By inspecting the displacement autocorrelation, we observe that the nanoparticle-structure interaction, indicated by the anticorrelation, is limited in the short-time stage and makes little contribution to the non-Gaussianity in the long-time stage. The main source of the non-Gaussianity can therefore be attributed to hopping diffusion that results in an exponential probability distribution with the large displacements, which may also explain certain processes dominated by rare events in the biological environment.

Covalent adaptable networks are polymers that can alter the arrangement of network connections by bond exchange reactions where an active unit attaches to an existing bond then kicks off its pre-existing peer to form a new bond. When the polymer is stretched, bond exchange reactions lead to stress relaxation and plastic deformation, or the so-called reforming. In addition, two pieces of polymers can be rejoined together without introducing additional monomers or chemicals on the interface, enabling welding and reprocessing. Although covalent adaptable networks have been researched extensively in the past, knowledge about the macromolecular level network alternations is limited. In this study, molecular dynamics simulations are used to investigate the macromolecular details of bond exchange reactions in a recently reported epoxy system. An algorithm for bond exchange reactions is first developed and applied to study a crosslinking network formed by epoxy resin DGEBA with the crosslinking agent tricarballylic acid. The trace of the active units is tracked to show the migration of these units within the network. Network properties, such as the distance between two neighboring crosslink sites, the chain angle, and the initial modulus, are examined after each iteration of the bond exchange reactions to provide detailed information about how material behaviors and macromolecular structure evolve. Stress relaxation simulations are also conducted. It is found that even though bond exchange reactions change the macroscopic shape of the network, microscopic network characteristic features, such as the distance between two neighboring crosslink sites and the chain angle, relax back to the unstretched isotropic state. Comparison with a recent scaling theory also shows good agreement.

The microfluidic passive control of microparticles largely relies on the hydrodynamic effects of the carrier media such as Newtonian fluids and viscoelastic fluids. Yet the viscoelastic/Newtonian interfacial effect has been scarcely investigated, especially for high-resolution particle separation. Here we report a microfluidic co-flow of Newtonian (water or PBS) and viscoelastic fluids (PEO) for the size-dependent separation of microparticles. The co-flow condition generates a stable viscoelastic/Newtonian interface, giving rise to the wall-directed elastic lift forces that compete with the center-directed lift forces, and efficiently hinders the migration of microparticles from the Newtonian to the viscoelastic fluid in a size-dependent manner. An almost complete separation of a binary mixture of 1 μm and 2 μm polystyrene particles is achieved by the co-flow of water and a very dilute PEO solution (100 ppm), whereas the sole use of water or PEO could not lead to an efficient separation. This co-flow microfluidic system is also applied for the separation of Staphylococcus aureus (1 μm) from platelets (2-3 μm) with >90% efficiencies and purities.

Liquid crystal elastomers (LCE) and magnetic soft materials are promising active materials in many emerging fields, such as soft robotics. Despite the high demand for developing active materials that combine the advantages of LCE and magnetic actuation, the lack of independent programming of the LCE nematic order and magnetization in a single material still hinders the desired multi-responsiveness. In this study, a ferromagnetic LCE (magLCE) ink with nematic order and magnetization is developed that can be independently programmed to be anisotropic, referred to as "dual anisotropy", via a customized 3D-printing platform. The magLCE ink is fabricated by dispersing ferromagnetic microparticles in the LCE matrix, and a 3D-printing platform is created by integrating a magnet with 3-DoF motion into an extrusion-based 3D printer. In addition to magnetic fields, magLCEs can also be actuated by heating sources (either environmental heating or photo-heating of the embedded ferromagnetic microparticles) with a high energy density and tunable actuation temperature. A programmed magLCE strip robot is demonstrated with enhanced adaptability to complex environments (different terrains, magnetic fields, and temperatures) using a multi-actuation strategy. The magLCE also has potential applications in mechanical memory, as demonstrated by the multistable mechanical metastructure array with remote writability and stable memory.

Viscoelastic microfluidics becomes an efficient and label-free hydrodynamic technology to enrich and separate micrometer-scale particles, including blood cells, circulating tumor cells, and bacteria. However, the manipulation of nanoscale particles by viscoelastic microfluidics remains a major challenge, because the viscoelastic force acting on the smaller particle decreases dramatically. In contrast to the commonly used polymer solutions of high molecular weight, herein we utilize the aqueous solutions of poly(ethylene oxide) (PEO) of low molecular weight with minimized shear thinning but sufficient elastic force for high-quality focusing and separation of various nanoparticles. The focusing efficiencies of 100 nm polystyrene (PS) nanoparticles and λ-DNA molecules are 84% and 85%, respectively, in a double spiral microchannel, without the aid of sheath flows. Furthermore, we demonstrate the size-based viscoelastic separation of two sets of binary mixtures-100/2000 nm PS particles and λ-DNA molecules/blood platelets-all achieving separation efficiencies of >95% in the same device. Our proposal technique would be a promising approach for enrichment/separation of the nanoparticles encountered in applications of analytical chemistry and nanotechnology.

Direct numerical simulations of turbulent channel flow at Reτ = 205 and 932 have been carried out to examine Taylor’s “frozen turbulence” hypothesis. The terms in Taylor’s hypothesis appear in the transport equation for instantaneous momentum (Navier-Stokes) in this flow. The additional terms, i.e., the additional convective acceleration term and the pressure gradient and viscous force terms, act to diminish the validity of Taylor’s hypothesis when they are relatively large compared to the Taylor’s hypothesis terms and are not in balance. A similar analysis has been applied to the transport equation for instantaneous vorticity. The additional terms in this equation, namely, the additional convective rates of change of vorticity terms, the stretching/compression/rotation of vorticity terms, and the viscous diffusion of vorticity terms, similarly act to diminish the validity of Taylor’s hypothesis when they are relatively large compared to the terms in the hypothesis and are not in balance. Where in the channel flow this diminishment occurs, and to what degree, and which of the non-Taylor’s hypothesis terms in the momentum and vorticity equations contribute most to this diminishment are unraveled here.

The dynamic spreading of a liquid droplet on micropillar-arrayed surfaces is experimentally investigated. A theoretical model is proposed to include energy dissipations raised from both the viscous resistance at mesoscale and the molecular friction at microscale in the triple-phase region. The scaling laws and spreading shape of the droplet change with the variation of the liquid viscosity because of the competition between these two mechanisms of energy dissipations at the moving contact line. The Laplace pressures at the interior corner and at the wavy contact line are the answers to the excess driving energy and the superwetting on pillar-arrayed surfaces. The formation and evolution of the bulk and the fringe are also analyzed in detail. Our results may help to understand the wetting dynamics on microtextured surfaces and assist the future design of engineered surfaces in practical applications.