New solutions to the fragmentation equation

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. The impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

Development and analysis of moments preserving finite volume schemes for multi-variate nonlinear breakage model

Modelling and simulation of collisional particle breakage mechanisms are crucial in several physical phenomena (asteroid belts, molecular clouds, raindrop distribution etc.) and process industries (chemical, pharmaceutical, material etc.). This paper deals with the development and analysis of schemes to numerically solve the multi-dimensional nonlinear collisional fragmentation model. Two numerical techniques are presented based on the finite volume discretization method. It is shown that the proposed schemes are consistent with the hypervolume conservation property. Moreover, the number preservation property law also holds for one of them. Detailed mathematical discussions are presented to establish the convergence analysis and consistency of the multi-dimensional schemes under predefined restrictions on the kernel and initial data. The proposed schemes are shown to be second-order convergent. Finally, several numerical computations (one-, two- and three-dimensional fragmentation) are performed to validate the numerical schemes.

Ultrasound frequency sonication facilitates high-throughput and uniform dissociation of cellular aggregates and tissues

A sample preparation step involving dissociation of tissues into their component cells is often required to conduct analysis of nucleic acids and other constituents from tissue samples. Frequently, the extracellular matrix and cell-cell adhesions are disrupted via treatment with a chemical dissociating reagent or various mechanical forces. In this work, a new, high-throughput, multiplexed method of dissociating tissues and cellular aggregates into single cells using ultrasound frequency bath sonication is explored and characterized. Different operating parameters are evaluated, and a treatment protocol with potential for uniform, high-throughput tissue dissociation is compared to the existing best chemical and orbital plate shaking protocol. Metrics such as percent dissociation, cellular recovery, average aggregate size, proportion of various aggregate sizes, membrane circularity, and cellular viability are subsequently assessed and found to be favorable. In optimized conditions, 53 ± 8% of 1 mm biopsy cores are dissociated within 30 min using sonication alone, surpassing leading high-throughput orbital plate shaking techniques five-fold. Chemical digestion is also 2 times more effective when complexed with sonication rather than orbital plate shaking. RNA content, quality, and expression are found to be superior to the standard protocol in terms of transcriptional preservation.

The influence of non-stationarity and interphase curvature on the growth dynamics of spherical crystals in a metastable liquid

This manuscript is concerned with the theory of nucleation and evolution of a polydisperse ensemble of crystals in metastable liquids during the intermediate stage of a phase transformation process. A generalized growth rate of individual crystals is obtained with allowance for the effects of their non-stationary evolution in unsteady temperature (solute concentration) field and the phase transition temperature shift appearing due to the particle curvature (the Gibbs-Thomson effect) and atomic kinetics. A complete system of balance and kinetic equations determining the transient behaviour of the metastability degree and the particle-radius distribution function is analytically solved in a parametric form. The coefficient of mutual Brownian diffusion in the Fokker-Planck equation is considered in a generalized form defined by an Einstein relation. It is shown that the effects under consideration substantially change the desupercooling/desupersaturation dynamics and the transient behaviour of the particle-size distribution function. The asymptotic state of the distribution function (its 'tail'), which determines the relaxation dynamics of the concluding (Ostwald ripening) stage of a phase transformation process, is derived. This article is part of the theme issue 'Transport phenomena in complex systems (part 1)'.

Nucleation and growth dynamics of ellipsoidal crystals in metastable liquids

When describing the growth of crystal ensembles from metastable solutions or melts, a significant deviation from a spherical shape is often observed. Experimental data show that the shape of growing crystals can often be considered ellipsoidal. The new theoretical models describing the transient nucleation of ellipsoidal particles and their growth with and without fluctuating rates at the intermediate stage of bulk phase transitions in metastable systems are considered. The nonlinear transport (diffusivity) of ellipsoidal crystals in the space of their volumes is taken into account in the Fokker-Planck equation allowing for fluctuating growth rates. The complete analytical solutions of integro-differential models of kinetic and balance equations are found and analysed. Our solutions show that the desupercooling dynamics is several times faster for ellipsoidal crystals as compared to spherical particles. In addition, the crystal-volume distribution function is lower and shifted to larger particle volumes when considering the growth of ellipsoidal crystals. What is more, this function is monotonically increasing to the maximum crystal size in the absence of fluctuations and is a bell-shaped curve when such fluctuations are taken into account. This article is part of the theme issue 'Transport phenomena in complex systems (part 1)'.

Modeling Food Particle Systems: A Review of Current Progress and Challenges

For many years, food engineers have attempted to describe physical phenomena such as heat and mass transfer in food via mathematical models. Still, the impact and benefits of computer-aided engineering are less established in food than in most other industries today. Complexity in the structure and composition of food matrices are largely responsible for this gap. During processing of food, its temperature, moisture, and structure can change continuously, along with its physical properties. We summarize the knowledge foundation, recent progress, and remaining limitations in modeling food particle systems in four relevant areas: flowability, size reduction, drying, and granulation and agglomeration. Our goal is to enable researchers in academia and industry dealing with food powders to identify approaches to address their challenges with adequate model systems or through structural and compositional simplifications. With advances in computer simulation capacity, detailed particle-scale models are now available for many applications. Here, we discuss aspects that require further attention, especially related to physics-based contact models for discrete-element models of food particle systems.

Conservative Finite Volume Schemes for Multidimensional Fragmentation Problems

In this article, a new numerical scheme for the solution of the multidimensional fragmentation problem is presented. It is the first that uses the conservative form of the multidimensional problem. The idea to apply the finite volume scheme for solving one-dimensional linear fragmentation problems is extended over a generalized multidimensional setup. The derivation is given in detail for two-dimensional and three-dimensional problems; an outline for the extension to higher dimensions is also presented. Additionally, the existing one-dimensional finite volume scheme for solving conservative one-dimensional multi-fragmentation equation is extended to solve multidimensional problems. The accuracy and efficiency of both proposed schemes is analyzed for several test problems.

Statistical Mechanics of Discrete Multicomponent Fragmentation

We formulate the statistics of the discrete multicomponent fragmentation event using a methodology borrowed from statistical mechanics. We generate the ensemble of all feasible distributions that can be formed when a single integer multicomponent mass is broken into fixed number of fragments and calculate the combinatorial multiplicity of all distributions in the set. We define random fragmentation by the condition that the probability of distribution be proportional to its multiplicity, and obtain the partition function and the mean distribution in closed form. We then introduce a functional that biases the probability of distribution to produce in a systematic manner fragment distributions that deviate to any arbitrary degree from the random case. We corroborate the results of the theory by Monte Carlo simulation, and demonstrate examples in which components in sieve cuts of the fragment distribution undergo preferential mixing or segregation relative to the parent particle.

On the approximate solutions of fragmentation equations

A numerical model based on the finite volume scheme is proposed to approximate the binary breakage problems. Initially, it is considered that the particle fragments are characterized by a single property, i.e. particle’s volume. We then investigate the extension of the proposed model for solving breakage problems considering two properties of particles. The efficiency to estimate the different moments with good accuracy and simple extension for multi-variable problems are the key features of the proposed method. Moreover, the mathematical convergence analysis is performed for one-dimensional problems. All mathematical findings and numerical results are validated over several test problems. For numerical validation, we propose the extension of Bourgade & Filbet (2008 Math. Comput. 77 , 851–882. ( doi:10.1090/S0025-5718-07-02054-6 )) model for solving two-dimensional pure breakage problems. In this aspect, numerical treatment of the two-dimensional binary breakage models using finite volume methods can be treated to be the first instance in the literature.

A note on the self-similar solutions to the spontaneous fragmentation equation

We provide a method to compute self-similar solutions for various fragmentation equations and use it to compute their asymptotic behaviours. Our procedure is applied to specific cases: (i) the case of mitosis, where fragmentation results into two identical fragments, (ii) fragmentation limited to the formation of sufficiently large fragments, and (iii) processes with fragmentation kernel presenting a power-like behaviour.

Renormalization of the fragmentation equation: exact self-similar solutions and turbulent cascades

Using an approach developed earlier for renormalization of the Boltzmann collision integral [Saveliev and Nanbu, Phys. Rev. E 65, 051205 (2002)], we derive an exact divergence form for the fragmentation operator. Then we reduce the fragmentation equation to the continuity equation in size space, with the flux given explicitly. This allows us to obtain self-similar solutions and to find the integral of motion for these solutions (we call it the bare flux). We show how these solutions can be applied as a description of cascade processes in three- and two-dimensional turbulence. We also suggested an empirical cascade model of impact fragmentation of brittle materials.

A study of the collisional fragmentation problem using the Gamma distribution approximation

The nonlinear fragmentation population balance formulation has been elevated in recent years from a prototype for studying nonlinear integro-differential equations to a vehicle for analyzing and understanding several physicochemical processes of technological interest. The so-called pure collisional fragmentation, which is the particular mode of nonlinear fragmentation induced by collisions between particles, is studied here. It is shown that the corresponding population balance equation admits large time asymptotic (self-similarity) solutions for homogeneous fragmentation and collision functions (kernels). The self-similar solutions are given in closed form for some simple kernels. Based on the shape of the self-similar solutions the method of moments with Gamma distribution approximation is employed for transient solution (from initial state to establishment of the asymptotic shape) of the collisional fragmentation equation. These solutions are presented for several sets of parameters and their behavior is discussed rather extensively. The present study is similar to the one has already been performed for the case of the much simpler linear fragmentation equation [G. Madras, B.J. McCoy, AIChE J. 44 (1998) 647].

The equal-size binary breakage problem: evolution toward a steady shape or periodic behavior?

During the last few years, the self-similar particle size distribution for a particle population undergoing breakage in equal size fragments has been derived using approximating, numerical, and analytical means. But very recently it was shown [N.V. Mantzaris, J. Phys. A Math. Gen. 38 (2005) 5111] through transient simulation of the breakage process that the particle size distribution in case of breakage in two equal fragments, never attains a steady shape, i.e., a self-similar form. The new results give rise to questions about the real meaning and utility of the previously derived self-similar distributions for these systems. The scope of the present work is to answer these questions and it is attempted using only analytical (exact) means for the solution of the transient breakage problem. In doing so, the very interesting and rich underlying structure and properties of the solutions of the equal size breakage problem (seemingly, very simple) are revealed. It appears that the utility of the known self-similar distributions for this particular problem has to be redefined but yet not entirely abandoned.

On the self-similar solution of fragmentation equation: Numerical evaluation with implications for the inverse problem

It is well known that the fragmentation equation admits self-similar solutions for evolving particle-size distributions (PSD); i.e., if the shape of PSD is independent of time after an initial transient period. Although an analytical derivations of the self-similar PSD cases have been studied extensively, results for cases requiring numerical solutions are rare. The aim of the present work is to fill this gap for the case of homogeneous breakage functions. The known analytical and approximate solutions for the self-similar PSD are reviewed and a general algorithm for the numerical solution is proposed. Results for a broad range of breakage functions (kernel and rate) are presented. Further, the work is focused on the sensitivity of the relation between self-similar PSD and breakage kernel and its influence on the inverse breakage problem, i.e., that of estimating the breakage kernel from experimental self-similar PSDs. Useful suggestions are made for tackling the inverse problem.