Particle Trajectory Simulations Research Articles

Simulation and reconstruction of particle trajectories in the CEPC drift chamber

The circular electron-positron collider (CEPC) is designed to precisely measure the properties of the Higgs boson, study electroweak interactions at the Z-boson peak, and search for new physics beyond the Standard Model. As a component of the 4th\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$4^{\ ext {th}}$$\\end{document} conceptual CEPC detector, the drift chamber facilitates the measurement of charged particles. This study implemented a Geant4-based simulation and track reconstruction for the drift chamber. For the simulation, detector construction and response were implemented and added to the CEPC simulation chain. The development of track reconstruction involves track finding using the combinatorial Kalman filter method and track fitting using the tool of GenFit. Using the simulated data, the tracking performance was studied. The results showed that both the reconstruction resolution and tracking efficiency satisfied the requirements of the CEPC experiment.

Design of a multi-carrier X-ray source for communication with energy modulation information.

X-ray communication is a kind of space communication technology which uses X-ray as information carrier. In order to improve the information transmission capacity, communication rate and anti-interference ability of X-ray communication, we proposes to design a novel multi-target X-ray source. The source is composed of a fast switching module of light channels based on FPGA technology and four photoelectric X-ray tubes with different target materials: Cr, Fe, Ni, and Cu. Using Geant4 software, we determined the optimal target thickness for each material, which enabled us to fully leverage the characteristic X-rays for multi-channel signal modulation transmission. Moreover, using CST software for particle trajectory simulation and optimization of the electron beam revealed that at a tube voltage of 20 kV, the focus area measures approximately 1.2 mm×1.2 mm. The simulations show that four kinds of spectra with high distinctiveness can be generated from the Cr, Fe, Ni, and Cu targets. Within a single modulation period, these spectra can be combined in various ways to create 16 different X-ray spectra signals, thereby increasing the number of communication elements and enhancing the information transmission rate.

Diffusion of noiseless active particles in a planar convection array.

We investigated, both analytically and numerically, the dynamics of a noiseless overdamped active particle in a square lattice of planar counter-rotating convection rolls. Below a first threshold of the self-propulsion speed, a fraction of the simulated particle's trajectories spatially diffuse around the convection rolls, whereas the remaining trajectories remain trapped inside the injection roll. We detected two chaotic diffusion regimes: (i) below a second, higher threshold of the self-propulsion speed, the particle performs a random motion characterized by asymptotic normal diffusion. Long superdiffusive transients were observed for vanishing small self-propulsion speeds. (ii) above that threshold, the particle follows chaotic running trajectories with speed and orientation close to those of the self-propulsion vector at injection and its dynamics is superdiffusive. Chaotic diffusion disappears in the ballistic limit of extremely large self-propulsion speeds.

Effect of eddy currents on metal fine particles in high gradient magnetic separation using a superconducting magnet

Abstract Dry high gradient magnetic separation is one of the key technologies expected in removing magnetic metal fine particles as contaminations in the production process of various powders. However, the effect of eddy currents on such metal fine particles in dry high gradient magnetic separation has not been adequately explored in previous research. In order to understand this phenomenon more comprehensively, we conducted particle trajectory simulations using a finite element analysis with numerical methods and actual magnetic separation experiments to investigate the influence of eddy currents on magnetic separation efficiency. The results show that for both paramagnetic and diamagnetic metal particles, eddy currents can favourably influence the capture efficiency. Furthermore, these effects depend on the particle size, electrical conductivity, and the rate of change in the magnetic field. Therefore, by taking eddy currents into account, magnetic separation systems can be better optimized and their performance further improved.

The motion of Brownian particles suspended in a non-uniform fluid flow

The effect of the type of fluid flow velocity profile on the motion of a Brownian particle in a circular tube has been studied using a Monte Carlo technique to simulate particle trajectories. Three types of fluid flow were studied: uniform (UF), parabolic (PF), and an initial uniform flow which evolves and may eventually attain the parabolic velocity profile, i.e. transient flow (TF). For UF and PF, radial and axial particle diffusion were simultaneously considered in some cases; for TF, where fluid dynamics equations had to be solved in addition to the Monte Carlo simulations, axial diffusion was neglected for fluid and particles, so that calculations had to be restricted to tube aspect ratios (radius/length) smaller than 0.1, a common scenario in most practical aerosol applications. A first series of simulations were performed to determine the overall particle penetration through the tube. For fixed tube geometry and particle diffusion coefficient, penetration in TF lied between those in PF (highest) and UF (lowest), in such a manner that, for a given value of the diffusion coefficient, the larger the degree of flow development (i.e., the smaller the values of the tube aspect ratio and the Reynolds number), the larger the penetration. These findings are coherent with those found in a former investigation dealing with the particle residence time in the tube. In a second series of simulation runs, calculations were done for selected initial locations of the particle within the tube entrance cross section, to examine the effect that the proximity of the particle starting position to the tube wall or to the tube axis has on variables such as penetration, mean first hitting time and length, mean first exit time, and mean axial velocity of lost and survived particles. When the fluid flow is not uniform but varies from place to place, the mean particle axial velocity is larger than that of the fluid. This counterintuitive fundamental result allows interpretation of the trends observed in penetration and residence time.

Fundamental Study on Recovery of Spilled Oil by Magnetic Pickering Emulsification

Heavy oil spills from stranded oil tankers are a global problem. The recovery of spilled oil must have a low environmental impact and must be capable of processing large quantities at high speeds. In this study, we investigated the application of magnetic separation for the recovery of spilled oil. Pickering emulsions were formed by heavy oil and magnetite, and recovery was carried out using a drum‐type permanent magnet. Pickering emulsions were successfully formed by stirring heavy oil A and magnetite, and magnetic separation experiments showed that up to 90% of the heavy oil A could be recovered. To recycle magnetite, a small amount of heavy oil was added to the recovered Pickering emulsions, and then a strong magnetic field of a superconducting bulk magnet was applied, which resulted in easy solid–liquid separation. The possibility of recycling magnetite after separation was also demonstrated. Finally, a small‐scale system using drum‐type permanent magnets was developed, and magnetic separation experiments were performed. The results showed that approximately 77% of the heavy oil A could be recovered. Particle trajectory simulations also demonstrated the recoverability of heavy oil at the sea surface. © 2024 Institute of Electrical Engineers of Japan and Wiley Periodicals LLC.

Research on particle separation microfluidic chip based on standing surface acoustic wave

Microparticle and cell separation is a key process in the study of cell properties, clinical diagnosis and disease treatment. In this paper, a particle separation method based on the combination of surface acoustic wave and microfluidic technology is proposed, which uses standing surface acoustic wave to achieve the separating of polystyrene microparticles of different sizes twice in a continuous flow. Through numerical simulation of the device model and particle trajectory simulation, the optimum particle separation parameters were determined. The particle separation microfluidic chip has successfully realized the particle separation. The separation yield of 5 μm particles is about 92.1%. The separation rate of 1 μm particles is about 83.3%.

On the parallelization of multipacting simulation codes for the design of particle accelerator components

Particle trajectory and collision simulation is a critical step of the design and construction of novel particle accelerator components. However it requires a huge computational effort which can slow down the design process. We started from a sequential simulation program which is used to study an event called Multipacting. Our work explains the physical problem that is simulated and the implications it can have on the behavior of the components. Then we analyze the original program’s operation to find the best options for parallelization. We first developed a parallel version of the Multipacting simulation and were able to accelerate the execution up to ∼35×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 35\ imes $$\\end{document} with 48 or 56 cores. In the best cases, parallelization efficiency was maintained up to 16 cores (∼95\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 95$$\\end{document}%) and the speed-up plateaus at around 40–48 cores. When this first parallelization effort was tried for multi-power simulations, we found that parallelism was severely limited with a maximum of 20×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$20\ imes $$\\end{document} speed-up. For this reason, we introduced a new method to improve the parallelization efficiency for this second use case. This method uses a shared processor pool for all simulations of electrons (OnePool). OnePool improved scalability by pushing the speed-up to over 32×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$32\ imes $$\\end{document}.

Ion energy distribution measurement device using a capillary plate with high-aspect ratio

An ion energy measurement method using a capillary plate with a high-aspect ratio is proposed to obtain ion energy distribution (IED). When a capillary plate is used, an electron-shading (ES) effect is introduced in which a potential difference occurs between the top and bottom of the capillary plate because of a difference in the mobility of ions and electrons. Electrons are repelled by the potential difference, and only ions reach the bottom of the capillary plate. Additionally, the measured IED at the bottom of the capillary plate is similar to the value at the bottom of the contact hole in the etch process. This means that the actual ion energy reaching the inside of the contact hole can be predicted in the etch process. In our experiment, the IED was obtained using the ES effect at various discharge conditions, such as RF bias power, pressure, and aspect ratio (L/D), of the capillary plate. The ES effect was verified by a particle trajectory simulation, and the ion current measured at the collector was compared with the increasing aspect ratio. The obtained IEDs were in good agreement with previous works.

Mechanisms controlling giant sea salt aerosol size distributions along a tropical orographic coastline

Abstract. Sea salt aerosol (SSA) is a naturally occurring phenomenon that arises from the breaking of waves and consequent bubble bursting on the ocean's surface. The resulting particles exhibit a bimodal distribution spanning orders of magnitude in size that introduces significant uncertainties when estimating the total annual mass of SSA on a global scale. Although estimates of mass and volume are significantly influenced by the presence of giant particles (dry radius >1 µm), effectively observing and quantifying these particles proves to be challenging. Additionally, uncertainties persist regarding the contribution of SSA production along coastlines, but preliminary studies suggest that coastal interactions may increase SSA particle concentrations by orders of magnitude. Moreover, our knowledge regarding the vertical distribution of SSA particles in the marine boundary layer remains limited, resulting in significant gaps in understanding the vertical mixing of giant aerosol particles and specific environmental conditions facilitating their dispersion. By addressing these uncertainties, particularly in regions where SSA particles constitute a substantial percentage of total aerosol loading, we can enhance our comprehension of the complex relationships between the air, sea, aerosols, and clouds. A case study conducted on the Hawaiian island of O`ahu offers insight into the influence of coastlines and orography on the production and vertical distribution of giant SSA size distributions. Along the coastline, the frequency of breaking waves is accelerated, serving as an additional source of SSA production. Furthermore, the steep island orography generates strong and consistent uplift under onshore trade wind conditions, facilitating vertical mixing of SSA particles along windward coastlines. To investigate this phenomenon, in situ measurements of SSA size distributions for particles with dry radii (rd) ≥ 2.8 µm were conducted for various altitudes, ranging from approximately 80 to 650 m altitude along the windward coastline and from 80 to 250 m altitude aboard a ship offshore. Comparing size distributions onshore and offshore confirmed significantly higher concentrations along the coastline, with 2.7–5.4 times greater concentrations than background open-ocean concentrations for supermicron particles. These size distributions were then analyzed in relation to environmental variables influencing SSA production and atmospheric dynamics. It was found that significant wave height exhibited the strongest correlation with changes in SSA size distributions. Additionally, simulated sea salt particle trajectories provided valuable insight into how production distance from the coastline impacts the horizontal and vertical advection of SSA particles of different sizes under varying trade wind speeds. Notably, smaller particles demonstrated reduced dependence on local wind speeds and production distance from the coastline, experiencing minimal dry deposition and high average maximum altitudes relative to larger particles. This research not only highlights the role of coastlines in enhancing the presence and vertical mixing potential of giant SSA particles, but also emphasizes how important it is to consider the influence of local factors on aerosol observations at different altitudes.

Automated design of a 3D passive microfluidic particle sorter.

Microfluidic chips that can sort mixtures of cells and other particles have important applications in research and healthcare. However, designing a sorter chip for a given application is a slow and difficult process, especially when we extend the design space from 2D into a 3D scenario. Compared to the 2D scenario, we need to explore more geometries to derive the appropriate design due to the extra dimension. To evaluate sorting performance, the simulation of the particle trajectory is needed. The 3D scenario brings particle trajectory simulation more challenges of runtime and collision handling with irregular obstacle shapes. In this paper, we propose a framework to design a 3D microfluidic particle sorter for a given application with an efficient 3D particle trajectory simulator. The efficient simulator enables us to simulate more samples to ensure the robustness of the sorting performance. Our experimental result shows that the sorter designed by our framework successfully separates the particles with the targeted size.

Feasibility Analysis of Calcium Carbonate Particle Trajectory Simulation in a Dual Horizontal Shaft Mixer.

This article aims to investigate the feasibility of using discrete element software EDEM 2022.0 to simulate the trajectory of artificial marble patterns in a dual horizontal shaft mixer. Research was conducted on the mixing uniformity of particles in the mixing chamber, and the optimal speed range for particle mixing was established. By simulating the trajectory of pigment particles, the trajectories of the particles at different positions of the stirring paddle were obtained, and the trajectories were compared with the measured results. In the study of uniform particle mixing, the Lacey index at different speeds was compared, and the optimal speed range was established between 40 RPM and 60 RPM. Based on this, the particle trajectory simulation found that the motion trajectories of particles at different positions of the stirring paddle varied significantly. The particles in the stirring paddle rod exhibit a gradual trend, in which they gradually decrease as they approach the head of the stirring paddle. Finally, the feasibility of this method was established by comparing the simulated and actual patterns through proportional replication of the mixing process, and it was discovered that the two were similar.

Particle focusing in a wavy channel

It is known that inertial lift forces can lead to particle focusing in channel flows; yet oscillatory straining effects have also been suggested as a mechanism for particle focusing in wavy channels. To explore the synergy between these two mechanisms, we analytically and experimentally investigate the focusing behaviour of rigid neutrally buoyant particles in a wavy channel. We decompose the particle-free channel flow into a primary Poiseuille flow and secondary eddies induced by the waviness. We calculate the perturbation of the particle on the particle-free flow and the resulting lateral lift force exerted on the particle using the method of matched asymptotic expansions. This yields a zeroth-order lift force arising from the Poiseuille flow and a first-order lift force due to the waviness of the channel. We further incorporate the inertial lift force into the Maxey–Riley equation and simulate particle trajectories in wavy channels. The interactions between the zeroth-order lift force and the particle-free flow largely determine the focusing locations. Experiments in wavy channels with varying amplitudes at channel Reynolds numbers ranging from 5 to 250 are consistent with the predictions of the focusing locations, which are mainly governed by the channel Reynolds number as well as the competition between the inertial lift and the oscillatory straining effects.

The effect of relative humidity on deposition pattern in inertial impactors: The role of particle elasticity and surface attraction

The role of relative humidity on particle bounce in an inertial impactor was investigated through a program of experiments and particle trajectory simulations. Inertial impactors are devices used to obtain particle size distributions by passing particle-laden air through a nozzle and collecting particles having sufficient inertia on a flat surface placed directly below the nozzle. Herein monodisperse hygroscopic particles impacted the flat hydrophobic surface of a single stage inertial impactor while varying the relative humidity of the flow. The results show that a circular deposition pattern occurs when the relative humidity is high. When the relative humidity is low secondary deposits beyond the circular deposit also occur. Particle trajectory simulations were performed where particle/surface interactions were quantified via the Hamaker constant (A) and the loss of kinetic energy via the coefficient of restitution (e). In a first-of-its-kind approach, (A, e) were iteratively adjusted until prominent features of the particle deposits in the experiments and simulations agreed, thereby providing values of (A,e) for each relative humidity. The results show that the observed deposition behavior is due to (i) the increase with relative humidity of the repulsion between the hygroscopic particles and the hydrophobic surface, (ii) multiple bounces along the particle trajectory, and (iii) kinetic energy loss at each particle bounce.

Fundamental Study of the Removal of Microplastic Fibers Using Swirling Flow and Magnetic Field

Microplastic fibers (MPFs) derived from laundry wastewater and fishing nets, etc. have been found in the environment around the world, raising concerns about their ecological impact. Particularly for laundry wastewater, some of the MPFs are not removed by sewage plants, and are discharged into the environment. We have studied the removal method of MPFs by magnetic separation in sewage plants, but there were some problems for practical applications such as the desorption of magnetite from MPFs and blockage of magnetic filters. To solve the problems, we proposed a novel magnetic separation system combining swirling flow and opposing magnets. In this study, we investigated whether MPFs can be concentrated by particle trajectory simulations and lab-scale experiments. The simulations and experiments showed similar trends that MPFs can be concentrated by the proposed magnetic separation system. The magnetic separation performance will be improved by using a higher magnetic field of opposed superconducting bulk magnets.

Removal of SUS304 fine particles transformed into martensite by the high gradient magnetic separation under dry condition

The interfusion of impurities such as metallic wear debris has been a major issue in the manufacturing process of foods, medicines and industrial products. Such debris originates with wearing of stainless-steel pipe joints or mechanical moving parts. Since the debris shows ferromagnetic properties by undergoing the strain-induced martensitic transformation of non-magnetic SUS304 (X5CrNi18-10), the magnetic separation system is much efficient to remove such debris from raw materials. In order to study the magnetic separation properties for the martensitic transformation fine particles, the several kinds of magnetic powder with different magnetization were prepared by controlling the amount of martensitic transformation by heat-treating SUS304 powders. The magnetic separation performance was evaluated by the high gradient magnetic separation (HGMS) experiments using the multilayer unidirectional magnetic wire filter under dry condition. The experimental results and FEM particle trajectory simulations revealed that the SUS304 powders transformed into martensite were removed at high speed of 0.1 m/s by the superconducting HGMS in a relatively low magnetic field.

Particle trajectory simulation using python and spreadsheet as an online learning alternative

Education today is required to utilize technological knowledge and skills in preparation for global competition. Along with the rapid development of technology, educators are required to develop learning alternatives. The purpose of this research is to create a particle trajectory simulation that is used as an alternative to online learning. The simulation uses Python programming language and Origin Pro assisted spreadsheet. Simulation in Python programming uses the Euler Cromer method to describe particle trajectories affected by electric and magnetic fields. This paper has successfully simulated particle trajectories affected by electric and magnetic fields with the Python programming language and Spreadsheet. The case where the motion of a charged particle is affected by a combination of electric and magnetic fields is when a positively charged particle moves perpendicular to the magnetic field, it will form a helical trajectory. However, when the electric field is in the direction of the magnetic field, the motion in the direction of the magnetic field will be accelerated by the electric force in the direction of the magnetic field which causes the helix to increase in width.

Inertial active ratchet: Simulation versus theory.

We present the inertial active dynamics of an Ornstein-Uhlenbeck particle in a piecewise sawtooth ratchet potential. Using the Langevin simulation and matrix continued fraction method (MCFM), the particle transport, steady-state diffusion, and coherence in transport are investigated in different parameter regimes of the model. Spatial asymmetry is found to be a key criterion for the possibility of directed transport in the ratchet. The MCFM results for net particle current of overdamped dynamics of the particle agree well with the simulation results. The simulated particle trajectories for the inertial dynamics and the corresponding position and velocity distribution functions reveal that the system passes through an activity-induced transition in the transport from the running phase to the locked phase of the dynamics. This is further corroborated by the mean square displacement (MSD) calculations, where the MSD gets suppressed with increase in the persistent duration of activity or self-propulsion in the medium and finally approaches zero for a very large value of self propulsion time. The nonmonotonic behavior of the particle current and Péclet number with self-propulsion time confirms that the particle transport and its coherence can be enhanced or reduced by fine tuning the persistent duration of activity. Moreover, for intermediate ranges of self-propulsion time as well as mass of the particle, even though the particle current shows a pronounced unusual maximum with mass, there is no enhancement in the Péclet number, instead the Péclet number decreases with mass, confirming the degradation of coherence in transport.

Simulation of Particle Trajectories in Gas Turbine Components and Assessment of Unsteady Effects Using an Efficient Eulerian-Lagrangian Technique

In recent years, CFD has proven to be a very useful asset to help with predicting complex flows in a wide range of situations, including multiphase and gas-particle flows. On this track, numerical modelling of particle-laden flows in multistage turbomachinery has become an important step in helping to analyse the behaviour of a discrete phase in gas turbines. Furthermore, unsteady effects due, for example, to rotor–stator interaction may have an effect on trajectories and capture efficiencies of the discrete phase. Unfortunately, computational times for transient simulations can be exceedingly high, especially if a discrete-phase needs also to be simulated. For this reason, this work reports a new method for the efficient and accurate simulation of particle-laden flows in gas turbine engines components. The Harmonic Balance Method is exploited to gain orders of magnitude speedup exploiting the idea that once the flow field has been embedded in the spectral basis, it can be reconstructed at any desired time. In this way, not only can the computational time needed to reach convergence of the flow field be dramatically reduced, but there is also no need to keep simulating the flow field during particle tracking. On the contrary, the continuous phase field can be retrieved at any desired time through flow reconstruction. This technique is conceptually simple, but, to the authors’ knowledge, has never been applied so far in particle-laden flow simulations and represents a novelty in the field. First, the implementation of the method is described, and details are given on how phase-lagged boundary conditions can be applied to flow and particles to further speed up the calculation. Then, some relevant case studies are presented to highlight the performance of the method.

Effects of Eulerian current, Stokes drift and wind while simulating surface drifter trajectories in the Baltic Sea

The simulation of Lagrangian drift is an important task in applications such as dispersion of pollutants, larvae and search and rescue activities. In this study, the Eulerian current, Stokes drift and wind effect on the simulation of observed drifters were analysed. The Lagrangian OceanParcels model was used, and the surface trajectories were assessed by comparison with 9 GPS drifter trajectories in the Gulf of Finland, Gulf of Riga and Lithuanian coast. The Normalised Cumulative Lagrangian Separation (NCLS) distance between the simulated and the satellite-tracked drifter trajectories, and the mean absolute error (MAE) were used as comparison metrics. The present study suggests the need to consider the Stokes drift and the wind factor in addition to the modelled Eulerian currents to obtain a better description of the trajectories of particles. By making these considerations, the OceanParcels model could adequately simulate particle trajectories in the sub-basins within the Baltic Sea. The realized model tests showed that motion of surface drifters are strongly controlled by the Stokes drift when the significant wave height is >1 m, whereas the wind component and the Eulerian currents are crucial when the significant wave height is <0.6 m or the wave (Stokes drift) directions do not match the wind direction.