Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

A generalized 1D particle transport method for convection diffusion reaction model

This work is devoted to the development of numerical method to deal with convection diffusion dominated problem with reaction term, non - stiff chemical reaction and stiff chemical reaction. The technique is based on the unifying Eulerian - Lagrangian schemes (particle transport method) under the framework of operator splitting method. In the computational domain, the particle set is assigned to solve the convection reaction subproblem along the characteristic curves created by convective velocity. At each time step, convection, diffusion and reaction terms are solved separately by assuming that, each phenomenon occurs separately in a sequential fashion. Moreover, adaptivities and projection techniques are used to add particles in the regions of high gradients (steep fronts) and discontinuities and transfer a solution from particle set onto grid point respectively. The numerical results show that, the particle transport method has improved the solutions of CDR problems. Nevertheless, the method is time consumer when compared with other classical technique e.g., method of lines. Apart from this advantage, the particle transport method can be used to simulate problems that involve movingsteep/smooth fronts such as separation of two or more elements in the system.

Influence of particle size on the properties of Pt-Ru/C catalysts prepared by a microemulsion method

Pt-Ru/C nanocatalysts were prepared by a microemulsion method using different values of water/surfactant molar ratio in order to get different particle sizes. Crystallite sizes and structural properties were determined by X-ray diffraction. Particle size and distribution were characterized by transmission electron microscopy and average composition was determined by energy-dispersive X-ray analysis. Thermogravimetric analysis was used to estimate the amount of supported metals. Differential scanning calorimetry measurements indicated the presence of hydrous ruthenium oxides in the as-prepared catalysts. Results for the oxidation of adsorbed CO as well as for methanol oxidation revealed significant differences in the behavior of the prepared catalysts. All together, the results demonstrate that the variation of particle size produces changes in other properties of the Pt-Ru/C catalysts and that to establish direct correlations between electrocatalytic activity and particle size is not possible because the effects of the different parameters cannot be separated. (c) 2007 the Electrochemical Society.

A new method to control particle size and particle size distribution of SnO2 nanoparticles for gas sensor applications

The controlled growth of SnO2 nanoparticles for gas sensor applications is reported by these authors. Nb2O5 additive is used to control nucleation and growth of the SnO2 (see Figure), which is synthesized by the polymeric precursor method. Preliminary gas sensing measurements are performed and it is demonstrated that the response time of the Nb2O5-doped SnO2 is faster than that of the undoped material.

Effect of reaction time on the particle size of ZnO and ZnO : Ce obtained by a sol-gel method

In this work particles of ZnO of size range 33-56 Angstrom were prepared by a sol-gel method. The effect of reaction time on the particle size of ZnO or ZnO:Ce was investigated by transmission electron microscopy measurements, UV-vis absorption and luminescence spectroscopy. A linear increase of the mean particle size is observed as a function of reaction time. The cerium-doped particles are bigger than the pure ZnO ones obtained at the same reaction time. A shift to lower energy at the maximum of the bands is observed in all absorption, emission and excitation spectra as a function of particle growth. From the absorption spectra the optical energy gap values (Eg) for these particles were determined. In the quantum size regime, Eg was found to decrease with particle growth.

Particle growth during calcination of polycation oxides synthesized by the polymeric precursors method

The particle-growth kinetics of sodium niobate and zirconium titanate powders that were processed by the polymeric precursors method were studied. The growth kinetics that were studied for the particle, in the final stage of crystallization, showed that the growth process occurs in two different stages. For temperatures <800°C, the particle-growth mechanism is associated with surface diffusion, with an activation energy in the range of 40-80 KJ/mol. For temprratures >800°C, particle growth is controlled by densification of the nanometric particle cluster and by a neck-size-controlled particle-growth mechanism. The results suggest that this behavior was typical of the synthesis method, because two different polycation oxides presented the same behavior.

Moderating effect of ammonia on particle growth and stability of quasi-monodisperse silver nanoparticles synthesized by the Turkevich method

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

Airborne-particle Abrasion Is a Minimally Invasive Method to Pre-Treat Sclerotic, Non-Carious Cervical Dentin

Objectives: To investigate substance loss and bond strength capacity of sclerotic, non-carious cervical dentin after airborne-particle abrasion or diamond bur preparation. Methods: Fifteen non-sclerotic dentin specimens were made from crowns of extracted human incisors of which the labial surfaces had been ground with silicon carbide papers (non-sclerotic control; Group 1). Forty-five sclerotic dentin specimens (n=15/group) were made from the labial, non-carious cervical root part of extracted human incisors and underwent either no pre-treatment (sclerotic control; Group 2), pre-treatment with airborne-particle abrasion (CoJet Prep [3M ESPE] and 50 µm aluminium oxide; Group 3), or with diamond bur preparation (40 µm grit size; Group 4). Substance loss after pre-treatment was measured in Groups 3 and 4. Subsequently, Scotchbond Universal (3M ESPE) and resin composite (CeramX [DENTSPLY DeTrey]) were applied on the treated dentin surfaces. The specimens were stored at 37°C and 100% humidity for 24 h. After storage, shear bond strength (SBS) was measured and data analyzed with nonparametric ANOVA followed by Wilcoxon rank sum tests. Results: Substance loss (medians) was 19 µm in Group 3 and 113 µm in Group 4. SBS-values (MPa; medians) in Group 2 (9.24) were significantly lower than in Group 1 (13.15; p=0.0069), Group 3 (13.05; p=0.01), and Group 4 (13.02; p=0.0142). There were no significant differences in SBS between Groups 1, 3, and 4 (p≥0.8063). Conclusion: Airborne-particle abrasion and diamond bur preparation restored bond strength of Scotchbond Universal to sclerotic dentin to the level of non-sclerotic dentin, with airborne-particle abrasion being less invasive than diamond bur preparation.

A Semi-Lagrangian Particle Level Set Finite Element Method for Interface Problems

We present a quasi-monotone semi-Lagrangian particle level set (QMSL-PLS) method for moving interfaces. The QMSL method is a blend of first order monotone and second order semi-Lagrangian methods. The QMSL-PLS method is easy to implement, efficient, and well adapted for unstructured, either simplicial or hexahedral, meshes. We prove that it is unconditionally stable in the maximum discrete norm, � · �h,∞, and the error analysis shows that when the level set solution u(t) is in the Sobolev space Wr+1,∞(D), r ≥ 0, the convergence in the maximum norm is of the form (KT/Δt)min(1,Δt � v �h,∞ /h)((1 − α)hp + hq), p = min(2, r + 1), and q = min(3, r + 1),where v is a velocity. This means that at high CFL numbers, that is, when Δt > h, the error is O( (1−α)hp+hq) Δt ), whereas at CFL numbers less than 1, the error is O((1 − α)hp−1 + hq−1)). We have tested our method with satisfactory results in benchmark problems such as the Zalesak’s slotted disk, the single vortex flow, and the rising bubble.

A packing computational method relating fractal particle size distribution and void fraction in granular media

The study of granular systems is of great interest to many fields of science and technology. The packing of particles affects to the physical properties of the granular system. In particular, the crucial influence of particle size distribution (PSD) on the random packing structure increase the interest in relating both, either theoretically or by computational methods. A packing computational method is developed in order to estimate the void fraction corresponding to a fractal-like particle size distribution.

A Packing Computational Method Relating Fractal Particle Size Distribution and Void Fraction in Granular Media

The study of granular systems is of great interest to many fields of science and technology. The packing of particles affects to the physical properties of the granular system. In particular, the crucial influence of particle size distribution (PSD) on the random packing structure increase the interest in relating both, either theoretically or by computational methods. A packing computational method is developed in order to estimate the void fraction corresponding to a fractal-like particle size distribution.

Coherent method of acceleration of particle bunches : paper presented at CERN symposium, 1956 /

"European Organization for Nuclear Research."