Local Particle Research Articles

Measuring and Simulating the Transient Packing Density During Ultrasound Directed Self‐Assembly and Vat Polymerization Manufacturing of Engineered Materials

AbstractUltrasound‐directed self‐assembly (DSA) uses ultrasound waves to organize and orient particles dispersed in a fluid medium into specific patterns. Combining ultrasound DSA with vat photopolymerization (VP) enables manufacturing materials layer‐by‐layer, wherein each layer the organization and orientation of particles in the photopolymer is controlled, which enables tailoring the properties of the resulting composite materials. However, the particle packing density changes with time and location as particles organize into specific patterns. Hence, relating the ultrasound DSA process parameters to the transient local particle packing density is important to tailor the properties of the composite material, and to determine the maximum speed of the layer‐by‐layer VP process. This paper theoretically derives and experimentally validates a 3D ultrasound DSA model and evaluates the local particle packing density at locations where particles assemble as a function of time and ultrasound DSA process parameters. The particle packing density increases with increasing particle volume fraction, decreasing particle size, and decreasing fluid medium viscosity is determined. Increasing the particle size and decreasing the fluid medium viscosity decreases the time to reach steady‐state. This work contributes to using ultrasound DSA and VP as a materials manufacturing process.

Validation study on a coarse-grained DEM-CFD simulation in a bead mill

A bead mill is often employed in powder processes. In the bead mill, the bead-particles move rapidly under high shear rate due to the moving structures. The discrete element method (DEM) coupled with the computational fluid dynamics (CFD) has been applied to simulate the bead-particle behavior for the optimization of the bead mill design and operating conditions. However, the existing DEM-CFD method faces essential problems regarding the substantial limit of the number of computational particles, stable calculations, and modeling of rapidly moving structures. To solve these problems, we newly investigate a numerical model combining the coarse-grained DEM-CFD, an implicit algorithm for the drag force term, and a scalar field-based wall boundary in the calculation of a bead mill commonly employed in industrial powder processes. In this study, our proposed model is shown to successfully simulate the macroscopic behavior of the bead-particles through the validation test of the coarse-grained DEM-CFD simulation in the bead mill system. Specifically, the macroscopic characteristics such as the particle location, the particle velocity, and total kinetic energy, are shown to agree well between the original particle system and coarse-grained particle systems. Consequently, this study illustrates the combination of the coarse-grained DEM-CFD, an implicit algorithm for the drag force term, and the scalar field-based wall boundary is essential for simulating the typically industrial bead mill, where an enormous number of bead-particles should be calculated in liquid under high shear rate.

Phantom study of a fully automatic radioactive seed placement robot for the treatment of skull base tumours

BackgroundInterstitial brachytherapy is a form of intensive local irradiation that facilitates the effective protection of surrounding structures and the preservation of organ functions, resulting in a favourable therapeutic response. As surgical robots can perform needle placement with a high level of accuracy, our team developed a fully automatic radioactive seed placement robot, and this study aimed to evaluate the accuracy and feasibility of fully automatic radioactive seed placement for the treatment of tumours in the skull base.MethodsA fully automatic radioactive seed placement robot was established, and 4 phantoms of skull base tumours were built for experimental validation. All the phantoms were subjected to computed tomography (CT) scans. Then, the CT data were imported into the Remebot software to design the preoperative seed placement plan. After the phantoms were fixed in place, navigation registration of the Remebot was carried out, and the automatic seed placement device was controlled to complete the needle insertion and particle placement operations. After all of the seeds were implanted in the 4 phantoms, postoperative image scanning was performed, and the results were verified via image fusion.ResultsA total of 120 seeds were implanted in 4 phantoms. The average error of seed placement was (2.51 ± 1.44) mm.ConclusionThis study presents an innovative, fully automated radioactive particle implantation system utilizing the Remebot device, which can successfully complete automated localization, needle insertion, and radioactive particle implantation procedures for skull base tumours. The phantom experiments showed the robotic system to be reliable, stable, efficient and safe. However, further research on the needle-soft tissue interaction and deformation mechanism of needle puncture is still needed.

Effect of distribution shape on the melting transition, local ordering, and dynamics in a model size-polydisperse two-dimensional fluid

Abstract We study the effect of particle size polydispersity (δ) on the melting transition (T *), local ordering, solid–liquid coexistence phase and dynamics of two-dimensional Lennard–Jones fluids up to moderate polydispersity by means of computer simulations. The particle sizes are drawn at random from the Gaussian (G) and uniform (U) distribution functions. For these systems, we further consider two different kinds of particles, viz., particles having the same mass irrespective of size, and in the other case the mass of the particle scales with its size. It is observed that with increasing polydispersity, the value of T * initially increases due to improved packing efficiency (ϕ) followed by a decrease and terminates at δ ≈ 8% (U-system) and 14% (G-system) with no significant difference for both mass types. The interesting observation is that the particular value at which ϕ drops suddenly coincides with the peak of the heat capacity (CP ) curve, indicating a transition. The quantification of local particle ordering through the hexatic order parameter (Q 6), Voronoi construction and pair correlation function reveals that the ordering decreases with increasing δ and T. Furthermore, the solid–liquid coexistence region for the G-system is shown to be comparatively wider in the T–δ plane phase diagram than that for the U system. Finally, the study of dynamics reveals that polydisperse systems relax faster compared to monodisperse systems; however, no significant qualitative differences, depending on the distribution type and mass polydispersity, are observed.

Reducing direct exposure to exhaled aerosol through a portable desktop fan

Vulnerable individuals close to infected people emitting a respiratory cloud containing infectious load can inhale a pathogen dose, experiencing a more severe impact on their health compared to other individuals breathing the mixed air in the same room. In crowded spaces, this issue is crucial. Employing local airflow patterns can reduce the proximity risk of inhalation and subsequent transmission across short distances. This study proposes an experimental and numerical analysis of a novel personal and portable device creating a short-range air barrier to transmitting airborne pathogens in proximity. The portable device adopts V-shaped air blades affecting the trajectory of the particle-laden respiratory cloud emitted by the respiratory system of the infected individual. Experimental results, supported by CFD analysis, indicate that controlling local airflow through the V-shaped jet significantly reduces local particle concentrations by more than 60%, compared to typical scenarios without a local airflow control.

Structural and dynamical equilibrium properties of hard board-like particles in parallel confinement.

We performed Monte Carlo and dynamic Monte Carlo simulations to model the diffusion of monodispersed suspensions composed of impenetrable cuboidal particles, specifically hard board-like particles (HBPs), in the presence of parallel hard walls. The impact of the walls was investigated by adjusting the size of the simulation box while maintaining constant packing fractions, fixed at η = 0.150, for systems consisting of HBPs with prolate, dual-shaped, and oblate geometries. We observed that increasing the distance between the walls led to the recovery of an isotropic bulk phase, while local particle organization near the walls remained stable. Due to their shape, oblate HBPs exhibit more efficient anchoring at wall surfaces compared to prolate shapes. The formation of nematic-like particle assemblies near the walls, confirmed by theoretical calculations based on density functional theory, significantly influenced local particle dynamics. This effect was particularly pronounced to the extent that a modest portion of cuboids near the walls tended to diffuse exclusively in planes parallel to the confinement, even more efficiently than observed in the bulk regions.

Spatially resolved measurement of the electrostatic charge of turbulent powder flows

This article reports on measurements of the electrostatic charge of particles in a turbulent duct flow. In contrast to previous charge measurements, which do not apply to turbulent flows or give only the sum of all particles’ charges, the new method resolves the charge of a turbulent powder flow spatially. The experiment consists of a particle tracking velocimetry (PTV) system and electrode plates that generate an electric field. By comparing particle velocities and accelerations with and without the electric field, the time-averaged local particle charge profile is derived. Spatially resolving the charge profiles unveiled bipolar particle flow. The average of the charge profiles agreed well with a conventional Faraday pail measurement, demonstrating the accuracy of our measurements. However, the peak value of the charge profiles was 76 times higher than the average of the particles’ charge.

The fluidization behavior of particles in supercritical water fluidized bed: A region-dependent analysis

The distribution and fluctuation of particle concentration in a supercritical water fluidized bed (SCWFB) are studied using a specially-designed capacitance probe measurement system. The measurement error of capacitance probe system is 2.48%. Results show that the particle concentration at the axial line fluctuates largely with time and the fluctuation is weakened when increasing the distance from axial line or decreasing particle diameter. At Uf/Umb=1.5, the largest difference between particle concentrations is 0.026, 0.04, 0.082 when the particles with a diameter of 0.124 mm, 0.198 mm, 0.235 mm are used, indicating a larger non-uniformity of larger particles. The probability density of local particle concentration mainly follows a unimodal distribution with the peak at right side. The main frequency of SCWFB ranges from 1.39 Hz to 2.73 Hz in this experiment and the main frequency decreases when larger particles are used.

Topological Phase Diagram of an Interacting Kitaev Chain: Mean Field versus DMRG Study

In this work, we study the topological phase transitions of a Kitaev chain generalized by the addition of nearest-neighbor Coulomb interaction. We show the presence of a robust topological phase as a function of the interaction strength and of the on-site energy with associated non-zero energy Majorana states localized at the chain edges. We provide an effective mean-field model that allows for the self-consistent computation of the mean value of the local particle number operator, and we also perform Density Matrix Renormalization Group numerical simulations based on a tensor network approach. We find that the two methods show a good agreement in reporting the phase transition between trivial and topological superconductivity. Temperature robustness within a physically relevant threshold has also been demonstrated. These findings shed light on an entire class of topological interacting one-dimensional systems in which the effects of residual Coulomb interactions play a relevant role.

Improving depth uncertainty in plenoptic camera-based velocimetry

This work describes the development of a particle tracking velocimetry (PTV) algorithm designed to improve three-dimensional (3D), three-component velocity field measurements using a single plenoptic camera. Particular focus is on mitigating the longstanding depth uncertainty issues that have traditionally plagued plenoptic particle image velocimetry (PIV) experiments by leveraging the camera’s ability to generate multiple perspective views of a scene in order to assist both particle triangulation and tracking. 3D positions are first estimated via light field ray bundling (LFRB) whereby particle rays are projected into the measurement volume using image-to-object space mapping. Tracking is subsequently performed independently within each perspective view, providing a statistical amalgamation of each particle’s predicted motion through time in order to help guide 3D trajectory estimation while simultaneously protecting the tracking algorithm from physically unreasonable fluctuations in particle depth positions. A synthetic performance assessment revealed a reduction in the average depth errors obtained by LFRB as compared to the conventional multiplicative algebraic reconstruction technique when estimating particle locations. Further analysis using a synthetic vortex ring at a magnification of − 0.6 demonstrated plenoptic-PIV capable of maintaining the equivalent of 0.1–0.15 voxel accuracy in the depth domain at a spacing to displacement ratio of 5.3–10.5, an improvement of 84–89% compared to plenoptic-PIV. Experiments were conducted at a spacing to displacement ratio of approximately 5.8 to capture the 3D flow field around a rotor within the rotating reference frame. The resulting plenoptic-PIV/PTV vector fields were evaluated with reference to a fixed frame stereoscopic-PIV (stereo-PIV) validation experiment. A systematic depth-wise (radial) component of velocity directed toward the wingtip, consistent with observations from prior literature and stereo-PIV experiments, was captured by plenoptic-PTV at magnitudes similar to the validation data. In contrast, the plenoptic-PIV did not discern any coherent indication of radial motion. Our algorithm constitutes a significant advancement in enhancing the functionality and versatility of single-plenoptic camera flow diagnostics by directly addressing the primary limitation associated with plenoptic imaging.Graphical abstract

Bioinspired 1T-MoS2-based aerogel beads for efficient freshwater harvesting in harsh environments

Freshwater scarcity is one of the most critical issues worldwide, particularly in arid regions, stemming from population growth and climate change. Inspired by the hydrophilic bump structures of desert beetles, 1T-MoS2-based aerogel beads with porous structures and CaCl2-crystal loading (termed as MoAB-m@CaCl2-n) were prepared for freshwater harvesting. Metallic-phase MoS2 nanospheres exhibit excellent photothermal conversion abilities, facilitating solar-driven water desorption and evaporation. Owing to the synergistic effect of its localized surface features, hydrophilic groups, and dispersive CaCl2 particles, MoAB-2@CaCl2-2 efficiently harvests water from atmosphere with a superior moisture adsorption capacity (0.18–0.82 g g−1) at a wide range of relative humidity (10 %-70 %). Under one-sun illumination, MoAB-2@CaCl2-2 demonstrates an outstanding solar-driven water evaporation rate of 2.25 kg m-2h−1. The water evaporation rate from soil (water content = 20 %) is 1.19 kg m-2h−1, which is sufficient for sustainable freshwater generation from the soil in arid regions. More importantly, the multifunctional MoAB-2@CaCl2-2-based homemade freshwater generation prototype delivers a certain amount of water harvesting (0.99 g g−1 day−1) on a rainy day and provides an impressive daily freshwater yield (53.7 kg m−2) under natural sunlight. The integrated device exhibits excellent efficiency and practicality and offers a feasible method for freshwater harvesting in harsh environments.

Unlocking stray light mysteries in the CoRot baffle with the time-of-flight method

Stray light (SL) has emerged as a primary limiting factor for space telescopes. Pre-launch testing is essential for validating performance and identifying potential issues. However, traditional methods do not enable the decomposition and identification of individual SL contributors. Consequently, when problems arise, resolving them often involves a cumbersome and risky trial-and-error approach. The time-of-flight (ToF) method was recently introduced, employing a pulsed laser source and ultrafast sensor to characterize individual SL contributors. A proof of concept was achieved using a simple three-lens system. In this paper, we apply the ToF method to a real space optical system: the spare model of the CoRoT baffle. We successfully measured individual SL contributors over a dynamic range of 10−11, identifying direct scattering on vane edges and two-step scattering paths. Our results provide a performance breakdown, differentiating intrinsic baffle SL from contributions arising from experimental conditions. Notably, the ToF method allowed us to discriminate air scattering, eliminating the need for expensive vacuum testing. The ToF provides unparallel insights, including defects identification. For instance, we identified the presence of localized dust particles causing significant SL. These results confirm the utility of the ToF method even for the most challenging space systems.

Comparative evaluation of electrical resistance tomography, positron emission particle tracking and high-speed imaging for analysing horizontal particle-liquid flow in a pipe

We evaluate three experimental techniques - electrical resistance tomography (ERT), positron emission particle tracking (PEPT) and high-speed imaging (HSI) – for analysing the local particle velocity field and spatial distribution in a horizontal particle-liquid pipe flow under varying conditions of solid concentration. A new ERT methodology is devised for estimating particle velocity, circumventing the limitations of the conventional cross-correlation technique. Furthermore, an enhanced HSI approach is introduced and systematically compared with PEPT and ERT. Results show that, under all conditions, PEPT provides the most accurate particle velocity field followed by HSI, whilst ERT yields the most accurate concentration field, followed by HSI. The enhanced HSI emerges as a simple cost-effective option compared to PEPT and ERT. A combined measurement approach using PEPT for local particle velocity and ERT for local concentration, however, delivers the best comprehensive two-phase flow characterization, highlighting potential synergies between these methods for complex flow studies.

Evaluation of array capacitive sensor for local concentration measurement of gas–solid particles flow by coupled fields based on CFD-DEM

An array capacitive sensor adopting two excitation modes is designed for measuring local particle concentration in gas–solid two-phase flow. Considering the difficulty in obtaining instantaneous flow parameters under dynamic experiments, a 3D dynamic simulation model coupling gas–solid two-phase fluid and electrostatic fields is proposed to quantitatively evaluate the sensor performance. In the coupling model, Computational Fluid Dynamics and Discrete Element Method are applied to establish two-phase fluid field for particle micro-distribution. The method of defining particle shape and position is used to set permittivity distribution in electrostatic field for coupling gas–solid fluid field and electrostatic field of array capacitive sensor. Using the method of reconstructed image and data fitting, the performance of two excitation modes are evaluated under a typical dynamic flow with structure of stratified flow. Experimental result validates the simulation conclusions and proves the effectiveness of the coupling simulation model for evaluating capacitive sensors in gas–solid flow systems.

Topography prediction of high-aspect ratio features milled with overlapping abrasive slurry jet footprints considering fluid confinement effects

Overlapping adjacent passes of abrasive slurry jets can be used to fabricate micro-pockets or blind micro-slots. This study investigates the propagation of surfaces milled using overlapping passes of a high-pressure abrasive slurry jet on a ductile 6061-T6 target over a wide range of aspect ratios. Features such as shallow pockets with low aspect ratios to deep slots with aspect ratios ∼1.5 were fabricated to investigate the correlations between the aspect ratio, degree of pass overlap, and number of repeating passes with the machined surface evolution characteristics. A comprehensive model was presented that combined computational fluid dynamics (CFD) with surface evolution to predict the evolving machined feature topography after each of the overlapping and repeating passes, and to provide insights into experimental observations. The model explicitly considered the CFD-predicted local particle impact angles and velocities as well as effects due to secondary impacts, high sidewall slopes, and stagnation effects. For larger degrees of overlap, measurements revealed a more than linear rate of increase in depth after each repeat of two adjacent overlapped passes, as opposed to a linear depth increase for smaller overlaps and blind pockets. Large overlaps were found to result in asymmetric features. Up to a critical aspect ratio of ∼0.65, the model predicted the surface evolution of the features to within <8.6% of those measured. Beyond that, randomness in the machined feature shape and size made the process challenging to control and the surface evolution difficult to predict. Nevertheless, the model was able to elucidate the reasons for the randomness and other observed phenomena such as the more than linear growth in etch rate, and the occurrence of feature asymmetry. The findings emphasize that flow confinement, increases in the jet's turbulent kinetic energy, and the formation of vortices at the bottom of the machined features are critical factors influencing the machining process, and that in these cases explicit modeling approaches like the one presented in this study must be used.

Exploring the nonclassical dynamics of the “classical” Schrödinger equation

The introduction of nonlinearities into the Schrödinger equation has been considered in the literature as an effective manner to describe the action of external environments or mean fields. Here, we explore the nonlinear effects induced by subtracting a term proportional to Bohm’s quantum potential to the usual (linear) Schrödinger equation, which generates the so-called “classical” Schrödinger equation. Although a simple nonlinear transformation allows us to recover the well-known classical Hamilton–Jacobi equation, by combining a series of analytical results (in the limiting cases) and simulations (whenever the analytical treatment is unaffordable), we find an analytical explanation to why the dynamics in the nonlinear “classical” regime is still strongly nonclassical. This is even more evident by establishing a one-to-one comparison between the Bohmian trajectories associated with the corresponding wave function and the classical trajectories that one should obtain. Based on these observations, it is clear that the transition to a fully classical regime requires extra conditions in order to remove any trace of coherence, which is the truly distinctive trait of quantum mechanics. This behavior is investigated in three paradigmatic cases, namely, the dispersion of a free propagating localized particle, the harmonic oscillator, and a simplified version of Young’s two-slit experiment.

The study on modulation mechanism for anti-corrosion and low-friction of R0-MoS2 by Ti-/Pb-doped model

Based on the first principles, the mechanism of the corrosion and friction failure of 2S-MoS2 (MoS2 stands for R0-MoS2 structure), of the anti-corrosion and friction modulation of Ti-/Pb-doped MoS2 are revealed in the atmosphere environment. For 2S-MoS2, O-doped and H2O-adsorbed 2S-MoS2 interface are formed with interlayer electron reconfiguration, which form the local MoO3 particles and increases their sliding barrier. Although the anti-corrosion role can be shown by trapping interlayer O2/H2O in 2Ti(23)-MoS2 +O2/H2O or by forming a clean Pb-doped interlayer in 2Pb(14)-MoS2 +O2/H2O, the tribological modulation effect occurs only 2Pb(14)-MoS2 +O2/H2O. Compared with 2S-MoS2 +O2/H2O, the tribological performance of 2Pb(14)-MoS2 +O2/H2O is improved by 31.4%/35.2%. These are attributed to 2Pb(14)-MoS2 +O2/H2O with O-passivated structure (S-O/O-S)/H2O-intercalated structure (S/H2O/S) that presents lower shear strength. However, 2Ti(23)-MoS2 +O2/H2O with the cross-linked structure (Ti-O-O-Ti/Ti-OH-Ti) shows higher shear strength. In addition, the ultra-low friction (improved by 26.3%/32.8%, compared with 2S-MoS2 in vacuum) and load-carrying capacity of 2Pb(14)-MoS2 +O2/H2O are also presented and 2Pb(14)-MoS2 +H2O is more excellent (improved by 8.8%, compared with 2Pb(14)-MoS2 +O2). In summary, the tribological modulation of Ti/Pb-doped MoS2 in the atmosphere environment is relative to the structure of their anti-corrosion interface.

Numerical extraction of three-dimensional acoustic polarizabilities for acoustically small scatterers

Polarizability is a convenient descriptor for scattering from low-ka inhomogeneities in both electromagnetism, where it originated, and in acoustics, where it is increasingly used in the design of metamaterial elements. In two or three dimensions, the polarizability is represented by a block matrix that couples the local pressure and particle velocity fields to the lowest-order components of the multipole expansion, usually truncated to dipole order. A fully dense acoustic polarizability matrix has components that couple spatially uniform pressure and velocity oscillations to dipole and monopole scattering, respectively. This unique coupling leads to acoustic bianisotropy, also known as Willis coupling, for a collection of scatterers. The design of Willis materials necessitates computationally efficient methods to extract all components of the polarizability matrix and, at present, only two-dimensional analytical extraction methods are found in the literature. In this work, we present an algorithm to extract all components of the three-dimensional polarizability matrix from a scatterer with arbitrary geometry and composition. The method presented here is numerically efficient and yields results that are in agreement with an analytically obtained benchmark polarizability, providing a useful tool for acoustic metamaterial design.

Mode entanglement in fermionic and bosonic Harmonium

Mode entanglement in many-body quantum systems is an active area of research. It provides crucial insight into the suitability of many-body systems for quantum information processing tasks. Local super-selection rules must be taken into account when assessing the amount of physically accessible entanglement. This requires amending well-established entanglement measures by incorporating local parity and local particle number constraints. In this paper, we report on mode entanglement present in the analytically solvable system of N-Harmonium. To the knowledge of the authors, this is the first analytic study of the physically accessible mode and mode-mode entanglement of an interacting many-body system in a continuous state space. We find that super-selection rules dramatically reduce the amount of physically accessible entanglement, which vanishes entirely in some cases. Our results strongly suggest the need to re-evaluate intra and inter-mode entanglement in other fermionic and bosonic systems.

Active acoustic metamaterials with independently programmable bulk modulus and full mass density tensor

Active acoustic metamaterials consist of unit cells with sensor-driver pairs that produce a coherent response to incident waves. The effective acoustic properties of the metamaterials depend on the gains programmed between the sensing and driving components. The strength of the monopole response to the local pressure determines the effective bulk modulus and the dipole response to the local particle velocity determines the effective mass density. The unit cells are controlled individually, so in theory the metamaterials can be scaled to an arbitrary size and geometry, but prior realizations were limited to only a few cells with 1D operation. Here, we present an active acoustic metamaterial of nine cells in 2D with programmable bulk modulus and mass density tensor. We demonstrate the ability to independently program the property components (bulk modulus and four mass density elements) for any desired set within stable limits, including previously unattainable acoustic properties necessary for the fabrication of transformation acoustics devices. We also demonstrate complex effective properties, specifically a lossy layer with impedance matched to the background, such that incident waves are absorbed with no reflection. The effective properties of the active metamaterial are validated by comparing the experimental total and scattered fields to ideal simulation results.