Search

Predicting particle segregation has remained challenging due to the lack of a general model for the segregation velocity that is applicable across a range of granular flow geometries. Here, a segregation velocity model for dense granular flows is developed by exploiting momentum balance and recent advances in particle-scale modelling of the segregation driving and drag forces over a wide range of particle concentrations, size and density ratios, and flow conditions. This model is shown to correctly predict particle segregation velocity in a diverse set of idealized and natural granular flow geometries simulated using the discrete element method. When incorporated in the well-established advection-diffusion-segregation formulation, the model has the potential to accurately capture segregation phenomena in many relevant industrial application and geophysical settings.

Currently, Si vacancy ($\mathrm{V}_\mathrm{Si}$) color centers in SiC are of significant interest due to their potential applications in quantum sensing and quantum communication. Meanwhile, the qualities of laser-induced color centers are well guaranteed. Femtosecond laser processing suffices for increasing the yield of $\mathrm{V}_\mathrm{Si}$ color centers in bulk materials and forms crater-shaped enriched regions on the surface. However, there is a notable absence of existing simulation methods to explain the mechanisms behind laser-assisted $\mathrm{V}_\mathrm{Si}$ color center generation. In this work, we design a three-dimensional molecular dynamics (3D-MD) model using an integral hemi-ellipsoidal shell mathematical model to simulate the interaction of Gaussian laser beams with bulk materials. Furthermore, we calculate the transmittance, absorption coefficient, refractive index, and reflectivity of $4H$-SiC. Then, the absorptance of a 1030 nm laser in 350 {\mu}m-thick $4H$-SiC material is abtained to simulate the energy loss during the actual processing. Finally, the study analyzes the movement trajectories of $\mathrm{V}_\mathrm{Si}$ color centers and explains the source of $\mathrm{V}_\mathrm{Si}$ on the surface. This analysis explains the reasons for the enrichment of color centers in the crater-shaped regions formed after laser deposition. Our work provides an effective 3D-MD modeling approach to study the processing mechanisms of laser interaction with semiconductor materials, offering insights into efficient $\mathrm{V}_\mathrm{Si}$ color center creation processes., Comment: 7 pages, 5 figures

The universal approximation property (UAP) of neural networks is a fundamental characteristic of deep learning. It is widely recognized that a composition of linear functions and non-linear functions, such as the rectified linear unit (ReLU) activation function, can approximate continuous functions on compact domains. In this paper, we extend this efficacy to the scenario of dynamical systems with controls. We prove that the control family $\mathcal{F}_1 = \mathcal{F}_0 \cup \{ \text{ReLU}(\cdot)\} $ is enough to generate flow maps that can uniformly approximate diffeomorphisms of $\mathbb{R}^d$ on any compact domain, where $\mathcal{F}_0 = \{x \mapsto Ax+b: A\in \mathbb{R}^{d\times d}, b \in \mathbb{R}^d\}$ is the set of linear maps and the dimension $d\ge2$. Since $\mathcal{F}_1$ contains only one nonlinear function and $\mathcal{F}_0$ does not hold the UAP, we call $\mathcal{F}_1$ a minimal control family for UAP. Based on this, some sufficient conditions, such as the affine invariance, on the control family are established and discussed. Our result reveals an underlying connection between the approximation power of neural networks and control systems., Comment: 19 pages

Particle segregation in dense flowing size-disperse granular mixtures is driven by gravity and shear, but predicting the associated segregation force due to both effects has remained an unresolved challenge. Here, a model of the combined gravity- and kinematics-induced segregation force on a single intruder particle is integrated with a model of the concentration dependence of the gravity-induced segregation force. The result is a general model of the net particle segregation force in flowing size-bidisperse granular mixtures. Using discrete element method simulations for comparison, the model correctly predicts the segregation force for a variety of mixture concentrations and flow conditions in both idealized and natural shear flows.

The study of universal approximation properties (UAP) for neural networks (NN) has a long history. When the network width is unlimited, only a single hidden layer is sufficient for UAP. In contrast, when the depth is unlimited, the width for UAP needs to be not less than the critical width $w^*_{\min}=\max(d_x,d_y)$, where $d_x$ and $d_y$ are the dimensions of the input and output, respectively. Recently, \cite{cai2022achieve} shows that a leaky-ReLU NN with this critical width can achieve UAP for $L^p$ functions on a compact domain ${K}$, \emph{i.e.,} the UAP for $L^p({K},\mathbb{R}^{d_y})$. This paper examines a uniform UAP for the function class $C({K},\mathbb{R}^{d_y})$ and gives the exact minimum width of the leaky-ReLU NN as $w_{\min}=\max(d_x,d_y)+\Delta (d_x, d_y)$, where $\Delta (d_x, d_y)$ is the additional dimensions for approximating continuous functions with diffeomorphisms via embedding. To obtain this result, we propose a novel lift-flow-discretization approach that shows that the uniform UAP has a deep connection with topological theory., Comment: Include errata of the previous versions

Deep learning has made significant applications in the field of data science and natural science. Some studies have linked deep neural networks to dynamic systems, but the network structure is restricted to the residual network. It is known that residual networks can be regarded as a numerical discretization of dynamic systems. In this paper, we back to the classical network structure and prove that the vanilla feedforward networks could also be a numerical discretization of dynamic systems, where the width of the network is equal to the dimension of the input and output. Our proof is based on the properties of the leaky-ReLU function and the numerical technique of splitting method to solve differential equations. Our results could provide a new perspective for understanding the approximation properties of feedforward neural networks.

In dense flowing bidisperse particle mixtures varying in size or density alone, smaller particles sink (driven by percolation) and lighter particles rise (driven by buoyancy). But when the particle species differ from each other in both size and density, percolation and buoyancy can either enhance (large/light and small/heavy) or oppose (large/heavy and small/light) each other. In the latter case, a local equilibrium condition can exist in which the two segregation mechanisms balance and particles remain mixed: this allows the design of minimally-segregating mixtures by specifying particle size ratio, density ratio, and mixture concentration. Using experimentally validated DEM simulations, we show that mixtures specified by the methodology remain relatively well-mixed in the thin rapid surface flows characteristic of heaps and tumblers commonly used in industry. Furthermore, minimally-segregating particle mixtures prepared in a fully segregated state in a tumbler mix over time and eventually reach a state of nearly uniform species concentration., Comment: 19 pages, 17 figures, submitted to AIChE Journal

Using simulations and a virtual-spring-based approach, we measure the segregation force, Fseg, over a range of size-bidisperse mixture concentrations, particle size ratios, and shear rates to develop a model for Fseg that extends its applicability from the well-studied non-interacting intruders regime to finite-concentration mixtures where cooperative phenomena occur. The model predicts the concentration below which the single intruder assumption applies and provides an accurate description of the pressure partitioning between species.

In bidisperse particle mixtures varying in size or density alone, large particles rise (driven by percolation) and heavy particles sink (driven by buoyancy). When the two particle species differ from each other in both size and density, the two segregation mechanisms either enhance (large/light and small/heavy) or oppose (large/heavy and small/light) each other. In the latter case, an equilibrium condition exists in which the two segregation mechanisms balance and the particles no longer segregate. This leads to a methodology to design non-segregating particle mixtures by specifying particle size ratio, density ratio, and mixture concentration to achieve the equilibrium condition. Using DEM simulations of quasi-2D bounded heap flow, we show that segregation is significantly reduced for particle mixtures near the equilibrium condition. In addition, the rise-sink transition for a range of particle size and density ratios matches the combined size and density segregation model predictions.

Flowing granular materials segregate due to differences in particle size (driven by percolation) and density (driven by buoyancy). Modelling the segregation of mixtures of large/heavy particles and small/light particles is challenging due to the opposing effects of the two segregation mechanisms. Using discrete element method (DEM) simulations of combined size and density segregation we show that the segregation velocity is well described by a model that depends linearly on the local shear rate and quadratically on the species concentration. Concentration profiles predicted by incorporating this segregation velocity model into a continuum advection-diffusion-segregation transport model match DEM simulation results well for a wide range of particle size and density ratios. Most surprisingly, the DEM simulations and the segregation velocity model both show that the segregation direction for a range of size and density ratios depends on the local species concentration. This leads to a methodology to determine the combination of particle size ratio, density ratio, and particle concentration for which a bidisperse mixture will not segregate.

Individual constituent balance equations are often used to derive expressions for species-specific segregation velocities in flows of dense granular mixtures. We propose a semiempirical expression for the interspecies momentum exchange in density-bidisperse granular flows as an extension of ideas from kinetic theory and compare it to a previous viscous drag approach that is analogous to particles settling in a fluid. The proposed model expands the range of the granular kinetic theory from short-duration binary collisions to the multiple enduring contacts characteristic of dense shear flows and incorporates the effects of particle friction, concentration ratio, and local flow conditions. The segregation velocities derived from the momentum balance equation using both interspecies drag models match the downward and upward segregation velocities of heavy and light particles obtained from DEM simulations through the flowing layer depth for different density ratios and constituent concentrations in confined shear flows. Predictions of the kinetic theory inspired approach are additionally compared to results from free surface heap flow simulations, and, again, a close match is observed.

Kinetic theory (KT) has been successfully used to model rapid granular flows in which particle interactions are frictionless and near elastic. However, it fails when particle interactions become frictional and inelastic. For example, the KT is not able to accurately predict the free cooling process of a vibrated granular medium that consists of inelastic frictional particles under microgravity. The main reason that the classical KT fails to model these flows is due to its inability to account for the particle surface friction and its inelastic behavior, which are the two most important factors that need be considered in modeling collisional granular flows. In this study, we have modified the KT model that is able to incorporate these two factors. The inelasticity of a particle is considered by establishing a velocity-dependent expression for the restitution coefficient based on many experimental studies found in the literature, and the particle friction effect is included by using a tangential restitution coefficient that is related to the particle friction coefficient. Theoretical predictions of the free cooling process by the classical KT and the improved KT are compared with the experimental results from a study conducted on an airplane undergoing parabolic flights without the influence of gravity [Y. Grasselli, G. Bossis, and G. Goutallier, EPL (Europhysics Letters) 86, 60007 (2009)]. Our results show that both the velocity- dependent restitution coefficient and the particle surface friction are important in predicting the free cooling process of granular flows; the modified KT model that integrates these two factors is able to improve the simulation results and led to a better agreement with the experimental results., Comment: Accepted for Physical Review E

Granular materials are characterized by large collections of discrete particles of sizes larger than one micron, where the particle-particle interactions are significantly more important than the particle-fluid interactions. These flows can be successfully modeled by the existing Kinetic Theory (KT) models when they are in the dilute regime with low particle-particle collision frequencies, yielding results that agree well with the simulation results of the event-driven hard sphere model or the more sophisticated soft-sphere Discrete Element Method (DEM). However, these KT models become less accurate for granular flows with soft particles (low particle stiffness) at high particle-particle collision frequencies when the predicted collision interval (the time of free flight for a particle prior to the next collision) is comparable to the collision duration; there is a large discrepancy between the results of these KT models and those from the DEM models. In this work we develop a new KT model that could be used to model granular systems of high collision frequencies with a finite particle stiffness. This is done by modifying the fluctuation energy dissipation rate to incorporate the ratio of collision duration to collision interval, a parameter that is determined by both the collision frequency and particle stiffness. We use a linear-spring-dashpot collision scheme to model the elastic potential energy in the system and to uncover the relationship between the constitutive relations of KT and the ratio of collision duration to collision interval.