Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Solar flare particle propagation--comparison of a new analytic solution with spacecraft measurements

A new analytic solution has been obtained to the complete Fokker-Planck equation for solar flare particle propagation including the effects of convection, energy-change, corotation, and diffusion with ĸ_r = constant and ĸ_Ɵ ∝ r². It is assumed that the particles are injected impulsively at a single point in space, and that a boundary exists beyond which the particles are free to escape. Several solar flare particle events have been observed with the Caltech Solar and Galactic Cosmic Ray Experiment aboard OGO-6. Detailed comparisons of the predictions of the new solution with these observations of 1-70 MeV protons show that the model adequately describes both the rise and decay times, indicating that ĸ_r = constant is a better description of conditions inside 1 AU than is ĸ_r ∝ r. With an outer boundary at 2.7 AU, a solar wind velocity of 400 km/sec, and a radial diffusion coefficient ĸ_r ≈ 2-8 x 10²⁰ cm²/sec, the model gives reasonable fits to the time-profile of 1-10 MeV protons from "classical" flare-associated events. It is not necessary to invoke a scatter-free region near the sun in order to reproduce the fast rise times observed for directly-connected events. The new solution also yields a time-evolution for the vector anisotropy which agrees well with previously reported observations.

In addition, the new solution predicts that, during the decay phase, a typical convex spectral feature initially at energy T_o will move to lower energies at an exponential rate given by T_KINK = T_oexp(-t/Ƭ_KINK). Assuming adiabatic deceleration and a boundary at 2.7 AU, the solution yields Ƭ_KINK ≈ 100h, which is faster than the measured ~200h time constant and slower than the adiabatic rate of ~78h at 1 AU. Two possible explanations are that the boundary is at ~5 AU or that some other energy-change process is operative.

The Gaussian mixture MCMC particle algorithm for dynamic cluster tracking

We present a novel filtering algorithm for tracking multiple clusters of coordinated objects. Based on a Markov chain Monte Carlo (MCMC) mechanism, the new algorithm propagates a discrete approximation of the underlying filtering density. A dynamic Gaussian mixture model is utilized for representing the time-varying clustering structure. This involves point process formulations of typical behavioral moves such as birth and death of clusters as well as merging and splitting. For handling complex, possibly large scale scenarios, the sampling efficiency of the basic MCMC scheme is enhanced via the use of a Metropolis within Gibbs particle refinement step. As the proposed methodology essentially involves random set representations, a new type of estimator, termed the probability hypothesis density surface (PHDS), is derived for computing point estimates. It is further proved that this estimator is optimal in the sense of the mean relative entropy. Finally, the algorithm's performance is assessed and demonstrated in both synthetic and realistic tracking scenarios. © 2012 Elsevier Ltd. All rights reserved.

Occurrence of particle debris field during focused Ga ion beam milling of glassy carbon

To explore the machining characteristics of glassy carbon by focused ion beam (FIB), particles induced by FIB milling on glassy carbon have been studied in the current work. Nano-sized particles in the range of tens of nanometers up to 400 nm can often be found around the area subject to FIB milling. Two ion beam scanning modes - slow single scan and fast repetitive scan - have been tested. Fewer particles are found in single patterns milled in fast repetitive scan mode. For a group of test patterns milled in a sequence, it was found that a greater number of particles were deposited around sites machined early in the sequence. In situ EDX analysis of the particles showed that they were composed of C and Ga. The formation of particles is related to the debris generated at the surrounding areas, the low melting point of gallium used as FIB ion source and the high contact angle of gallium on glassy carbon induces de-wetting of Ga and the subsequent formation of Ga particles. Ultrasonic cleaning can remove over 98% of visible particles. The surface roughness (Ra) of FIB milled areas after cleaning is less than 2 nm. © 2010.

A study of the suppression of the high-temperature helium embrittlement in an oxide-particle dispersion strengthened alloy

In this paper, an investigation on the micro-structure of an Fe-base oxide-dispersion-strengthened (ODS) alloy irradiated with high-energy Ne-20 ions to different doses at a temperature around 0.5T(m) (T-m is the melting point of the alloy) is presented. Investigation with the transmission electron microscopy found that the accelerated growth of voids at grain-boundaries, which is usually a concern in conventional Fe-base alloys under conditions of inert-gas implantation, was not observed in the ODS alloy irradiated even to the highest dose (12000 at.ppm Ne). The reason is ascribed to the enhanced recombination of point defects and strong trapping of Ne atoms at the interfaces of the nano-scale oxide particles in grains. The study showed that ODS alloys have good resistance to the high-temperature inter-granular embrittlement due to inert-gas accumulation, exhibiting prominence of application in harsh situations of considerable helium production at elevated temperatures like in a fusion reactor.

Identified particle production, azimuthal anisotropy, and interferometry measurements in Au plus Au collisions at root s(NN)=9.2 GeV

We present the first measurements of identified hadron production, azimuthal anisotropy, and pion interferometry from Au + Au collisions below the nominal injection energy at the BNL Relativistic Heavy-Ion Collider (RHIC) facility. The data were collected using the large acceptance solenoidal tracker at RHIC (STAR) detector at root s(NN) = 9.2 GeV from a test run of the collider in the year 2008. Midrapidity results on multiplicity density dN/dy in rapidity y, average transverse momentum < p(T)>, particle ratios, elliptic flow, and Hanbury-Brown-Twiss (HBT) radii are consistent with the corresponding results at similar root s(NN) from fixed-target experiments. Directed flow measurements are presented for both midrapidity and forward-rapidity regions. Furthermore the collision centrality dependence of identified particle dN/dy, < p(T)>, and particle ratios are discussed. These results also demonstrate that the capabilities of the STAR detector, although optimized for root s(NN) = 200 GeV, are suitable for the proposed QCD critical-point search and exploration of the QCD phase diagram at RHIC.

Volumetric consequences of particle loss by grading entropy.

Chemical and biological processes, such as dissolution in gypsiferous sands and biodegradation in waste refuse, result in mass or particle loss, which in turn lead to changes in solid and void phase volumes and grading. Data on phase volume and grading changes have been obtained from oedometric dissolution tests on sand–salt mixtures. Phase volume changes are defined by a (dissolution-induced) void volume change parameter (Λ). Grading changes are interpreted using grading entropy coordinates, which allow a grading curve to be depicted as a single data point and changes in grading as a vector quantity rather than a family of distribution curves. By combining Λ contours with pre- to post-dissolution grading entropy coordinate paths, an innovative interpretation of the volumetric consequences of particle loss is obtained. Paths associated with small soluble particles, the loss of which triggers relatively little settlement but large increase in void ratio, track parallel to the Λ contours. Paths associated with the loss of larger particles, which can destabilise the sand skeleton, tend to track across the Λ contours.

Particle Image Velocimetry Investigation of a finite amplitude pressure wave

Particle image velocimetry is used to study the motion of gas within a duct subject to the passage of a finite amplitude pressure wave. The wave is representative of the pressure waves found in the exhaust systems of internal combustion engines. Gas particles are accelerated from stationary to 150 m/s and then back to stationary in 8 ms. It is demonstrated that gas particles at the head of the wave travel at the same velocity across the duct cross section at a given point in time. Towards the tail of the wave viscous effects are plainly evident causing the flow profile to tend towards parabolic. However, the instantaneous mean particle velocity across the section is shown to match well with the velocity calculated from a corresponding measured pressure history using 1D gas dynamic theory. The measured pressure history at a point in the duct was acquired using a high speed pressure transducer of the type typically used for engine research in intake and exhaust systems. It is demonstrated that these are unable to follow the rapid changes in pressure accurately and that they are prone to resonate under certain circumstances.

An exactly solvable four transition point model

Heavy particle collisions, in particular low-energy ion-atom collisions, are amenable to semiclassical JWKB phase integral analysis in the complex plane of the internuclear separation. Analytic continuation in this plane requires due attention to the Stokes phenomenon which parametrizes the physical mechanisms of curve crossing, non-crossing, the hybrid Nikitin model, rotational coupling and predissociation. Complex transition points represent adiabatic degeneracies. In the case of two or more such points, the Stokes constants may only be completely determined by resort to the so-called comparison- equation method involving, in particular, parabolic cylinder functions or Whittaker functions and their strong-coupling asymptotics. In particular, the Nikitin model is a two transition-point one-double-pole problem in each half-plane corresponding to either ingoing or outgoing waves. When the four transition points are closely clustered, new techniques are required to determine Stokes constants. However, such investigations remain incomplete, A model problem is therefore solved exactly for scattering along a one-dimensional z-axis. The energy eigenvalue is b(2)-a(2) and the potential comprises -z(2)/2 (parabolic) and -a(2) + b(2)/2z(2) (centrifugal/centripetal) components. The square of the wavenumber has in the complex z-plane, four zeros each a transition point at z = +/-a +/- ib and has a double pole at z = 0. In cases (a) and (b), a and b are real and unitarity obtains. In case (a) the reflection and transition coefficients are parametrized by exponentials when a(2) + b(2) > 1/2. In case (b) they are parametrized by trigonometrics when a(2) + b(2) <1/2 and total reflection is achievable. In case (c) a and b are complex and in general unitarity is not achieved due to loss of flux to a continuum (O'Rourke and Crothers, 1992 Proc. R. Sec. 438 1). Nevertheless, case (c) coefficients reduce to (a) or (b) under appropriate limiting conditions. Setting z = ht, with h a real constant, an attempt is made to model a two-state collision problem modelled by a pair of coupled first-order impact parameter equations and an appropriate (T) over tilde-tau relation, where (T) over tilde is the Stueckelberg variable and tau is the reduced or scaled time. The attempt fails because (T) over tilde is an odd function of tau, which is unphysical in a real collision problem. However, it is pointed out that by applying the Kummer exponential model to each half-plane (O'Rourke and Crothers 1994 J. Phys. B: At. Mel. Opt. Phys. 27 2497) the current model is in effect extended to a collision problem with four transition points and a double pole in each half-plane. Moreover, the attempt in itself is not a complete failure since it is shown that the result is a perfect diabatic inelastic collision for a traceless Hamiltonian matrix, or at least when both diagonal elements are odd and the off-diagonal elements equal and even.

Modulation of surface charge, particle size and morphological properties of chitosan-TPP nanoparticles intended for gene delivery

This work investigates the polyanion initiated gelation process in fabricating chitosan-TPP (tripolyphosphate) nanoparticles in the size range of 100-250 nm intended to be used as carriers for the delivery of gene or protein macromolecules. It demonstrates that ionic gelation of cationic chitosan molecules offers a flexible and easily controllable process for systematically and predictably manipulating particle size and surface charge which are important properties in determining gene transfection efficacy if the nanoparticles are used as non-viral vectors for gene delivery, or as delivery carriers for protein molecules. Variations in chitosan molecular weight, chitosan concentration, chitosan to TPP weight ratio and solution pH value were examined systematically for their effects on nanoparticle size, intensity of surface charge, and tendency of particle aggregation so as to enable speedy fabrication of chitosan nanoparticles with predetermined properties. The chitosan-TPP nanoparticles exhibited a high positive surface charge across a wide pH range, and the isoelectric point (IEP) of the nanoparticles was found to be at pH 9.0. Detailed imaging analysis of the particle morphology revealed that the nanoparticles possess typical shapes of polyhedrons (e.g., pentagon and hexagon), indicating a similar crystallisation mechanism during the particle formation and growth process. This study demonstrates that systematic design and modulation of the surface charge and particle size of chitosan-TPP nanoparticles can be readily achieved with the right control of critical processing parameters, especially the chitosan to TPP weight ratio. (c) 2005 Elsevier B.V. All rights reserved.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

Comparison of a novel spray congealing procedure with emulsion-based methods for the micro-encapsulation of water-soluble drugs in low melting point triglycerides

The particle size characteristics and encapsulation efficiency of microparticles prepared using triglyceride materials and loaded with two model water-soluble drugs were evaluated. Two emulsification procedures based on o/w and w/o/w methodologies were compared to a novel spray congealing procedure. After extensive modification of both emulsification methods, encapsulation efficiencies of 13.04% tetracycline HCl and 11.27% lidocaine HCl were achievable in a Witepsol (R)-based microparticle. This compares to much improved encapsulation efficiencies close to 100% for the spray congealing method, which was shown to produce spherical particles of similar to 58 mu m. Drug release studies from a Witepsol (R) formulation loaded with lidocaine HCl showed a temperature-dependent release mechanism, which displayed diffusion-controlled kinetics at temperatures similar to 25 degrees C, but exhibited almost immediate release when triggered using temperatures close to that of skin. Therefore, such a system may find application in topical semi-solid formulations, where a temperature-induced burst release is preferred.

Nonlinear theory of solitary waves associated with longitudinal particle motion in lattices - Application to longitudinal grain oscillations in a dust crystal

The nonlinear aspects of longitudinal motion of interacting point masses in a lattice are revisited, with emphasis on the paradigm of charged dust grains in a dusty plasma (DP) crystal. Different types of localized excitations, predicted by nonlinear wave theories, are reviewed and conditions for their occurrence (and characteristics) in DP crystals are discussed. Making use of a general formulation, allowing for an arbitrary (e.g. the Debye electrostatic or else) analytic potential form phi(r) and arbitrarily long site-to-site range of interactions, it is shown that dust-crystals support nonlinear kink-shaped localized excitations propagating at velocities above the characteristic DP lattice sound speed v(0). Both compressive and rarefactive kink-type excitations are predicted, depending on the physical parameter values, which represent pulse- (shock-)like coherent structures for the dust grain relative displacement. Furthermore, the existence of breather-type localized oscillations, envelope-modulated wavepackets and shocks is established. The relation to previous results on atomic chains as well as to experimental results on strongly-coupled dust layers in gas discharge plasmas is discussed.

Rao-blackwellized particle filter for human appearance and position tracking

In human motion analysis, the joint estimation of appearance, body pose and location parameters is not always tractable due to its huge computational cost. In this paper, we propose a Rao-Blackwellized Particle Filter for addressing the problem of human pose estimation and tracking. The advantage of the proposed approach is that Rao-Blackwellization allows the state variables to be splitted into two sets, being one of them analytically calculated from the posterior probability of the remaining ones. This procedure reduces the dimensionality of the Particle Filter, thus requiring fewer particles to achieve a similar tracking performance. In this manner, location and size over the image are obtained stochastically using colour and motion clues, whereas body pose is solved analytically applying learned human Point Distribution Models.

Particle swarm optimization applied to integrated demand response resources scheduling

The concept of demand response has a growing importance in the context of the future power systems. Demand response can be seen as a resource like distributed generation, storage, electric vehicles, etc. All these resources require the existence of an infrastructure able to give players the means to operate and use them in an efficient way. This infrastructure implements in practice the smart grid concept, and should accommodate a large number of diverse types of players in the context of a competitive business environment. In this paper, demand response is optimally scheduled jointly with other resources such as distributed generation units and the energy provided by the electricity market, minimizing the operation costs from the point of view of a virtual power player, who manages these resources and supplies the aggregated consumers. The optimal schedule is obtained using two approaches based on particle swarm optimization (with and without mutation) which are compared with a deterministic approach that is used as a reference methodology. A case study with two scenarios implemented in DemSi, a demand Response simulator developed by the authors, evidences the advantages of the use of the proposed particle swarm approaches.