Identifying Multiquark Hadrons from Heavy Ion Collisions

Identifying hadronic molecular states and/or hadrons with multiquark components either with or without exotic quantum numbers is a long-standing challenge in hadronic physics. We suggest that studying the production of these hadrons in relativistic heavy ion collisions offers a promising resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their structures. Using the coalescence model for hadron production, we find that, compared to the case of a nonexotic hadron with normal quark numbers, the yield of an exotic hadron is typically an order of magnitude smaller when it is a compact multiquark state and a factor of 2 or more larger when it is a loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic states could be produced and realistically measured in these experiments.

K*(0) production in Cu plus Cu and Au plus Au collisions at root s(NN)=62.4 GeV and 200 GeV

We report on K*(0) production at midrapidity in Au + Au and Cu + Cu collisions at root s(NN) = 62.4 and 200 GeV collected by the Solenoid Tracker at the Relativistic Heavy Ion Collider detector. The K*(0) is reconstructed via the hadronic decays K*(0) -> K(+)pi(-) and (K*(0)) over bar -> K(+)pi(-). Transverse momentum, p(T), spectra are measured over a range of p(T) extending from 0.2 GeV/c up to 5 GeV/c. The center-of-mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, < p(T)>, are presented. The measured N(K*(0))/N(K) and N(phi)/N(K*(0)) ratios favor the dominance of rescattering of decay daughters of K*(0) over the hadronic regeneration for the K*(0) production. In the intermediate p(T) region (2.0 < p(T) < 4.0 GeV/c), the elliptic flow parameter, v(2), and the nuclear modification factor, R(CP), agree with the expectations from the quark coalescence model of particle production.

D (sJ) (2317) meson production at RHIC

We study the production of D (sJ) (2317) mesons in relativistic heavy ion collisions using the quark coalescence model. The predicted D (sJ) (2317) abundance depends sensitively on the quark structure of the D (sJ) (2317) meson. We have also evaluated the absorption cross sections of the D (sJ) (2317) meson by pi, rho, kaon and K* in a phenomenological hadronic model. We find that the final yield of D (sJ) (2317) mesons remains sensitive to its initial number produced from the quark-gluon plasma, providing thus the possibility of studying the quark structure of the D (sJ) (2317) meson and its production mechanism in relativistic heavy ion collisions.

A BEM model applied to failure analysis of multi-fractured structures

Due to manufacturing or damage process, brittle materials present a large number of micro-cracks which are randomly distributed. The lifetime of these materials is governed by crack propagation under the applied mechanical and thermal loadings. In order to deal with these kinds of materials, the present work develops a boundary element method (BEM) model allowing for the analysis of multiple random crack propagation in plane structures. The adopted formulation is based on the dual BEM, for which singular and hyper-singular integral equations are used. An iterative scheme to predict the crack growth path and crack length increment is proposed. This scheme enables us to simulate the localization and coalescence phenomena, which are the main contribution of this paper. Considering the fracture mechanics approach, the displacement correlation technique is applied to evaluate the stress intensity factors. The propagation angle and the equivalent stress intensity factor are calculated using the theory of maximum circumferential stress. Examples of multi-fractured domains, loaded up to rupture, are considered to illustrate the applicability of the proposed method. (C) 2011 Elsevier Ltd. All rights reserved.

Modeling the kinetics of the coalescence of water droplets in crude oil emulsions subject to an electric field, with the cellular automata technique

In this study, the concept of cellular automata is applied in an innovative way to simulate the separation of phases in a water/oil emulsion. The velocity of the water droplets is calculated by the balance of forces acting on a pair of droplets in a group, and cellular automata is used to simulate the whole group of droplets. Thus, it is possible to solve the problem stochastically and to show the sequence of collisions of droplets and coalescence phenomena. This methodology enables the calculation of the amount of water that can be separated from the emulsion under different operating conditions, thus enabling the process to be optimized. Comparisons between the results obtained from the developed model and the operational performance of an actual desalting unit are carried out. The accuracy observed shows that the developed model is a good representation of the actual process. (C) 2010 Published by Elsevier Ltd.

Coalescence of Water Droplets in Crude Oil Emulsions: Analytical Solution

In petroleum refineries, water is used in desalting units to remove the salt contained in crude oil. Typically, 7% of the volume of hot crude oil is water, forming a water-and-oil emulsion. The emulsion flows between two electrodes and is subjected to an electric field. The electrical forces promote the coalescence of small droplets of water dispersed in crude oil, and these form bigger droplets. This paper calculates the forces acting on the droplets, highlighting particularly the mechanisms proposed for droplet-droplet coalescence under the influence of an applied electric field. Moreover, a model is developed in order to calculate the displacement speed of the droplets and the time between droplet collisions. Thus, it is possible to simulate and optimize the process by changing the operational variables (temperature, electrical field, and water quantity). The main advantage of this study is to show that it is feasible to increase the volume of water recycled in desalting processes, thus reducing the use of freshwater and the generation of liquid effluents in refineries.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

Population balance modelling of granulation with a physically based coalescence kernel

It was previously published by the authors that granules can either coalesce through Type I (when granules coalesce by viscous dissipation in the surface liquid layer before their surfaces touch) or Type II (when granules are slowed to a halt during rebound, after their surfaces have made contact) (AIChE J. 46 (3) (2000) 529). Based on this coalescence mechanism, a new coalescence kernel for population balance modelling of granule growth is presented. The kernel is constant such that only collisions satisfying the conditions for one of the two coalescence types are successful. One constant rate is assigned to each type of coalescence and zero is for the case of rebound. As the conditions for Types I and II coalescence are dependent on granule and binder properties, the coalescence kernel is thus physically based. Simulation results of a variety of binder and granule materials show good agreement with experimental data. (C) 2002 Elsevier Science Ltd. All rights reserved.

Experimental study of drop-interface coalescence in the presence of polymer stabilisers

An experimental study has been carried out to characterise the performance of polymer stabilisers, partially hydrolysed polyvinyl acetate (PVAc), used in suspension polymerisation processes. The stabilisers are ranked by their ability to stabilise the dispersion characterised by the median coalescence time of a single drop with its homophase at a planar liquid/liquid interface. Results show that the stability of the dispersion relates closely to the molecular properties of the PVAcs. Other conditions being equal, PVAcs with higher molecular weights or lower degrees of hydrolysis can better stabilise a liquid-liquid dispersion. The stability of the dispersion also depends strongly on where the PVAc resides. The presence of a PVAc in the dispersed phase significantly reduces stability. Consistent with results reported in the literature, considerable scatter has been observed on the coalescence times of identical drops under the same conditions. An explanation for the scatter is also proposed in the paper, based on the classical Reynolds model for film thinning. (C) 2002 Elsevier Science B.V. All rights reserved.

A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Étude de la coalescence de nanogouttelettes par dynamique moléculaire

Ce travail est consacré à l’étude de la coalescence de gouttelettes liquides à l’échelle du nanomètre. Nous nous sommes intéressés principalement à l’évolution du changement topologique des gouttelettes à partir de la rupture des surfaces au moment du contact initial jusqu’à la coalescence complète. Nous utilisons la dynamique moléculaire afin de simuler plusieurs types de gouttelettes soit en utilisant le potentiel empirique de type Stillinger-Weber pour des gouttelettes de silicium (l-Si) en 3 dimensions et le modèle Embedded Atom Method pour des gouttelettes de cuivre liquide (l-Cu) en 2d, quasi-2d (disques) et 3 dimensions. Qualitativement, toutes les simulations démontrent une coalescence similaire indépendamment de la dimension de calcul (2d à 3d), de la taille et de la température initiale des gouttelettes. La coalescence évolue par une déformation rapide des surfaces sans mixage important entre les atomes des deux gouttelettes initiales. De plus, nous étudions l’évolution du col de coalescence formé lors du contact initial entre les gouttelettes et, pour les systèmes en 3d, nous observons une transition claire d’un régime visqueux vers un régime inertiel du rayon de ce col, tel que suggéré par des modèles théoriques. Pour les gouttelettes de cuivre nous observons exactement le comportement des prédictions analytiques et confirmons que le premier régime suit un comportement visqueux sans aplatissement local des gouttelettes. La situation est différente pour les gouttelettes de l-Si où nous observons un effet plus grand, par rapport aux prédictions analytiques, du rayon et de la température initiale des gouttelettes sur l’évolution du col de coalescence. Nous suggérons que les paramètres décrivant l’évolution de la coalescence dépendent des propriétés des matériaux utilisés contrairement à la théorie universelle couramment utilisée.

Processus de coalescence dans une population subdivisée avec possibilité de coalescences multiples

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Phonon confinement model to measure the average sizes of anatase nanoparticles synthesized by a solvothermal method using H2O2

In this work, crystalline titanium dioxide (TiO2) nanoparticles with variable average crystallite sizes (e.g., 8 nm) and surface areas (e.g., 192 m² g-1) were synthesized in pure anatase phase using H2O2 to reduce the hydrolysis rate of the titanium ions. An isopropanol (IP) solution was employed as the reaction medium. The TiO2 nanoparticles were characterized by powder X-ray diffraction analysis (XRD), Raman spectroscopy and transmission electron microscopy (TEM). By changing the synthesis parameters it was possible to control nanoparticle size and avoid the coalescence process. A dependence of the Raman wavenumber on the nanocrystal sizes was determined, which is quite useful for a quick check of the size of TiO2 nanocrystals.

Subcritical Propagation and Coalescence of Oil-Filled Cracks: Getting the Oil Out of Low-Permeability Source Rocks

We use a fracture mechanics model to study subcritical propagation and coalescence of single and collinear oil-filled cracks during conversion of kerogen to oil. The subcritical propagation distance, propagation duration, crack coalescence and excess oil pressure in the crack are determined using the fracture mechanics model together with the kinetics of kerogen-oil transformation. The propagation duration for the single crack is governed by the transformation kinetics whereas the propagation duration for the multiple collinear cracks may vary by two orders of magnitude depending on initial crack spacing. A large amount of kerogen (>90%) remains unconverted when the collinear cracks coalesce and the new, larger cracks resulting from coalescence will continue to propagate with continued kerogen-oil conversion. The excess oil pressure on the crack surfaces drops precipitously when the collinear cracks are about to coalesce, and crack propagation duration and oil pressure on the crack surfaces are strongly dependent on temperature. Citation: Jin, Z.-H., S. E. Johnson, and Z. Q. Fan (2010), Subcritical propagation and coalescence of oil-filled cracks: Getting the oil out of low-permeability source rocks, Geophys. Res. Lett., 37, L01305, doi:10.1029/2009GL041576.

A fracture model for pearlitic steel bars using a cohesive model

The fracture of ductile materials, such as metals, is usually explained with the theory of nucleation, growth and coalescence of microvoids. Based on this theory, many numerical models have been developed, with a special mention to Gurson-type models. These models simulate mathematically the physical growth of microvoids, leading to a progressive development of the internal damage that takes place during a tensile test. In these models, the damage starts to develop in very early stages of the test. Tests carried out by the authors suggest that, in the case of some eutectoid steels such as those used for manufacturing prestressing steel wires, the internal damage that takes place as a result of the growth of microvoids is only noticeable in very late stages of the tensile test. In the authors’ opinion, using a cohesive model as a failure criterion may be interesting in this case; a cohesive model only requires two parameters to be defined, with the fracture energy being one of them, which can be obtained experimentally. In addition to this, given that it is known that the stress triaxiality has a strong influence on the fracture of ductile materials, a cohesive model whose parameters are affected by the value of the stress triaxiality can be considered. This work presents a fracture model for steel specimens in a tensile test, based on a cohesive behaviour and taking into account the effect of stress triaxiality, which is different at each point of the fracture plane.