Mathematical Modelling of the Dynamics of Granular Materials in a Rotating Cylinder

The current study is aimed at the development of a theoretical simulation tool based on Discrete Element Method (DEM) to 'interpret granular dynamics of solid bed in the cross section of the horizontal rotating cylinder at the microscopic level and subsequently apply this model to establish the transition behaviour, mixing and segregation.The simulation of the granular motion developed in this work is based on solving Newton's equation of motion for each particle in the granular bed subjected to the collisional forces, external forces and boundary forces. At every instant of time, the forces are tracked and the positions velocities and accelarations of each partcle is The software code for this simulation is written in VISUAL FORTRAN 90 After checking the validity of the code with special tests, it is used to investigate the transition behaviour of granular solids motion in the cross section of a rotating cylinder for various rotational speeds and fill fraction.This work is hence directed towards a theoretical investigation based on Discrete Element Method (DEM) of the motion of granular solids in the radial direction of the horizontal cylinder to elucidate the relationship between the operating parameters of the rotating cylinder geometry and physical properties ofthe granular solid.The operating parameters of the rotating cylinder include the various rotational velocities of the cylinder and volumetric fill. The physical properties of the granular solids include particle sizes, densities, stiffness coefficients, and coefficient of friction Further the work highlights the fundamental basis for the important phenomena of the system namely; (i) the different modes of solids motion observed in a transverse crosssection of the rotating cylinder for various rotational speeds, (ii) the radial mixing of the granular solid in terms of active layer depth (iii) rate coefficient of mixing as well as the transition behaviour in terms of the bed turnover time and rotational speed and (iv) the segregation mechanisms resulting from differences in the size and density of particles.The transition behaviour involving its six different modes of motion of the granular solid bed is quantified in terms of Froude number and the results obtained are validated with experimental and theoretical results reported in the literature The transition from slumping to rolling mode is quantified using the bed turnover time and a linear relationship is established between the bed turn over time and the inverse of the rotational speed of the cylinder as predicted by Davidson et al. [2000]. The effect of the rotational speed, fill fraction and coefficient of friction on the dynamic angle of repose are presented and discussed. The variation of active layer depth with respect to fill fraction and rotational speed have been investigated. The results obtained through simulation are compared with the experimental results reported by Van Puyvelde et. at. [2000] and Ding et at. [2002].The theoretical model has been further extended, to study the rmxmg and segregation in the transverse direction for different particle sizes and their size ratios. The effect of fill fraction and rotational speed on the transverse mixing behaviour is presented in the form of a mixing index and mixing kinetics curve. The segregation pattern obtained by the simulation of the granular solid bed with respect to the rotational speed of the cylinder is presented both in graphical and numerical forms. The segregation behaviour of the granular solid bed with respect to particle size, density and volume fraction of particle size has been investigated. Several important macro parameters characterising segregation such as mixing index, percolation index and segregation index have been derived from the simulation tool based on first principles developed in this work.

Study on Granular Dynamics in Vertically Vibrated Beds using Tracking Technique

The study of granular material is of great interest to many researchers in both engineering and science communities. The importance of such a study derives from its complex rheological character and also its significant role in a wide range of industrial applications, such as coal, food, plastics, pharmaceutical, powder metallurgy and mineral processing. A number of recent reports have been focused on the physics of non-cohesive granular material submitted to vertical vibration in either experimental or theoretical approaches. Such a kind of system can be used to separate, mix and dry granular materials in industries. It exhibits different instability behaviour on its surface when under vertical vibration, for example, avalanching, surface fluidization and surface wave, and these phenomena have attracted particular interest of many researchers. However, its fundamental understanding of the instability mechanism is not yet well-understood. This paper is therefore to study the dynamics of granular motion in such a kind of system using Positron Emission Particle Tracking (PEPT), which allows the motion of a single tracer particle to be followed in a non-invasive way. Features of the solids motion such as cycle frequency and dispersion index were investigated via means of authors’ specially-written programmes. Regardless of the surface behaviour, particles are found to travel in rotational movement in horizontal plane. Particle cycle frequency is found to increase strongly with increasing vibration amplitude. Particle dispersion also increased strongly with vibration amplitude. Horizontal dispersion is observed to always exceed vertical dispersion.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

Molecular motion in nanocomposites of poly(ethylene oxide) and montmorillonite

Motion of chains of poly(ethylene oxide) within the interlayer spacing of 2:1 phyllosilicate/montmorillonite was studied with H-1 and C-13 NMR spectroscopy. Measurements of the H-1 NMR line widths and relaxation times across a large temperature range were used to determine the effect of bulk thermal transitions on polymer chain motion within the nanocomposites. The results were consistent with previous reports of low apparent activation energies of motion. Details of the frequency and geometry of motion were obtained from a comparison of the C-13 cross-polarity/magic-angle spinning spectra and relaxation times of the nanocomposite with those of the pure polymer. (C) 2001 John Wiley & Sons, Inc.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

Evaluation of the sonically induced narrowing of nuclear magnetic resonance spectra of solids

The Nuclear Magnetic Resonance (NMR) spectra of liquids contain a wealth of quantitative information that may be derived, for instance, from chemical shifts and spin-spin couplings. The available information depends on the incoherent rapid molecular motion that causes complicating effects present in the solid state to average to zero. Whereas liquid state NMR spectra show narrow lines, the corresponding NMR spectra from the solid state are normally composed of exceedingly broad resonance lines due to highly restricted molecular motion. It is, therefore, difficult to obtain directly as detailed information from the spectra of solids as from those derived from the liquid state. Studies on a new technique (SINNMR, the sonically induced narrowing of the NMR spectra of solids) to remove line broadening effects in the NMR spectra of the solid state are reported within this thesis. SINNMR involves narrowing the NMR absorptions from solid particles by irradiating them with ultrasound when they are suspended in a support liquid. It is proposed that ultrasound induces incoherent motion of the suspended particles, producing motional characteristics of the particles similar to those of rather large molecules. The first report of apparently successful experiments involving SINNMR[1] emphasised both the irreproducibility of the technique and the uncertainty regarding its true origin. If SINNMR can be made reproducible and the effect definitively attributed to the sonically induced incoherent motional averaging of particles, the technique could offer a simple alternative to the now classical magic-angle spinning (MAS) NMR[2] and the recently reported dynamic angle spinning (DAS)[3] and double rotation (DOR)[4] techniques. Evidence is presented in this thesis to support the proposal that ultrasound may be used to narrow the NMR spectral resonances from solids by inducing incoherent motion of particles suspended in support liquids and, additionally, for some solids, by inducing rotational motion of molecular constituents in the lattices of solids. Successful SINNMR line narrowing using 20 kHz ultrasound is reported for a variety of samples: including trisodium orthophosphate, polytetrafluoroethylene and aluminium alloys. Investigations of SINNMR line narrowing in trisodium phosphate have revealed the relationship between ultrasonic power, particle size and support liquid density for the production of optimum SINNMR conditions. It is also proposed that the incoherent motion of particles induced by 20 kHz ultrasound can originate from interactions between acoustically induced cavitation microjets and particles.

Sonically induced narrowing of the NMR spectra of solids and other novel NMR techniques

A typical liquid state NMR spectrum is composed of a number of discrete absorptions which can be readily interpreted to yield detailed information about the chemical environment of the nuclei found within the sample. The same cannot be said about the spectra of solid samples. For these the absorptions are typically broad, featureless and yield little information directly. This situation may be further exacerbated by the characteristically long T1 values of nuclei bound within a solid lattice which, consequently, require long inter-sequence delays that necessitate lengthy experiments. This work attempts to address both of these inherent problems. Classically, the resolution of the broad-line spectra of solids into discrete resonances has been achieved by imparting to the sample coherent rotation about specific axes in relation to the polarising magnetic field, as implemented in the magic-angle spinning (MAS) [1], dynamic angle spinning (DAS) [2] and double rotation (DOR) [3] NMR experiments. Recently, an alternative method, sonically induced narrowing of the NMR spectra of solids (SINNMR) [4], has been reported which yields the same well resolved solid-state spectra as the classic solid-state NMR experiments, but which achieves the resolution of the broad-line spectra through the promotion of incoherent motion in a suspension of solid particles. The first part of this work examines SINNMR and, in particular, concentrates on ultrasonically induced evaluation, a phenomenon which is thought to be essential to the incoherent averaging mechanism. The second part of this work extends the principle of incoherent motion, implicit in SINNMR, to a new genre of particulate systems, air fluidized beds, and examines the feasibility of such systems to provide well resolved solid state NMR spectra. Samples of trisodium phosphate dodecahydrate and of aluminium granules are examined using the new method with partially resolved spectra being reported in the case of the latter.

Pain Intensity And Cervical Range Of Motion In Women With Myofascial Pain Treated With Acupuncture And Electroacupuncture: A Double-blinded, Randomized Clinical Trial.

Acupuncture stimulates points on the body, influencing the perception of myofascial pain or altering physiologic functions. The aim was to evaluate the effect of electroacupuncture (EAC) and acupuncture (AC) for myofascial pain of the upper trapezius and cervical range of motion, using SHAM acupuncture as control. Sixty women presenting at least one trigger point at the upper trapezius and local or referred pain for more than six months were randomized into EAC, AC, and SHAM groups. Eight sessions were scheduled and a follow-up was conducted after 28 days. The Visual Analog Scale assessed the intensity of local and general pain. A fleximeter assessed cervical movements. Data were analyzed using paired t or Wilcoxon's tests, ANOVA or Friedman or Kruskal-Wallis tests and Pearson's correlation (α=0.05). There was reduction in general pain in the EAC and AC groups after eight sessions (P<0.001). A significant decrease in pain intensity occurred for the right trapezius in all groups and for the left trapezius in the EAC and AC groups. Intergroup comparisons showed improvement in general pain in the EAC and AC groups and in local pain intensity in the EAC group (P<0.05), which showed an increase in left rotation (P=0.049). The AC group showed increases in inclination (P=0.005) sustained until follow-up and rotation to the right (P=0.032). EAC and AC were effective in reducing the pain intensity compared with SHAM. EAC was better than AC for local pain relief. These treatments can assist in increasing cervical range of motion, albeit subtly.

Análise do desempenho muscular do quadríceps e dos isquiotibiais em função da série temporal e da amplitude de movimento de atletas amadoras de futsal feminino = : Analysis of muscle performance of the quadriceps and hamstrings as a function of the times series and range of motion of amateurfemale futsal athletes

Universidade Estadual de Campinas . Faculdade de Educação Física

Calibração de câmeras para análises subaquáticas de movimentos = : Cameras calibration for underwater motion analysis

Universidade Estadual de Campinas . Faculdade de Educação Física

Proposição de um método de análise de movimentos de jogadores de tênis de campo a partir de vídeos televisivos = : Proposition of a method for tennis players motion analysis on broadcast videos

Universidade Estadual de Campinas . Faculdade de Educação Física

Determination of Three-Dimensional Left Ventricle Motion to Analyze Ventricular Dyssyncrony in SPECT Images

A method to compute three-dimension (3D) left ventricle (LV) motion and its color coded visualization scheme for the qualitative analysis in SPECT images is proposed. It is used to investigate some aspects of Cardiac Resynchronization Therapy (CRT). The method was applied to 3D gated-SPECT images sets from normal subjects and patients with severe Idiopathic Heart Failure, before and after CRT. Color coded visualization maps representing the LV regional motion showed significant difference between patients and normal subjects. Moreover, they indicated a difference between the two groups. Numerical results of regional mean values representing the intensity and direction of movement in radial direction are presented. A difference of one order of magnitude in the intensity of the movement on patients in relation to the normal subjects was observed. Quantitative and qualitative parameters gave good indications of potential application of the technique to diagnosis and follow up of patients submitted to CRT.

A frequency approach to identifying asteroid families II. Families interacting with nonlinear secular resonances and low-order mean-motion resonances

Aims. In an earlier paper we introduced a new method for determining asteroid families where families were identified in the proper frequency domain (n, g, g + s) ( where n is the mean-motion, and g and s are the secular frequencies of the longitude of pericenter and nodes, respectively), rather than in the proper element domain (a, e, sin(i)) (semi-major axis, eccentricity, and inclination). Here we improve our techniques for reliably identifying members of families that interact with nonlinear secular resonances of argument other than g or g + s and for asteroids near or in mean-motion resonant configurations. Methods. We introduce several new distance metrics in the frequency space optimal for determining the diffusion in secular resonances of argument 2g - s, 3g - s, g - s, s, and 2s. We also regularize the dependence of the g frequency as a function of the n frequency (Vesta family) or of the eccentricity e (Hansa family). Results. Our new approaches allow us to recognize as family members objects that were lost with previous methods, while keeping the advantages of the Carruba & Michtchenko (2007, A& A, 475, 1145) approach. More important, an analysis in the frequency domain permits a deeper understanding of the dynamical evolution of asteroid families not always obtainable with an analysis in the proper element domain.