A conceptual model for sheet-flow drawn from rapid granular flow theories

33rd IAHR Congress: Water Engineering for a Sustainable Environment

Granular flow and design studies in flighted rotating drums

Magdeburg, Univ., Fak. für Verfahrens- und Systemtechnik, Diss., 2013

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.

Cluster identification and characterization in the riser of a circulating fluidized bed from numerical simulation results

A methodology of identification and characterization of coherent structures mostly known as clusters is applied to hydrodynamic results of numerical simulation generated for the riser of a circulating fluidized bed. The numerical simulation is performed using the MICEFLOW code, which includes the two-fluids IIT`s hydrodynamic model B. The methodology for cluster characterization that is used is based in the determination of four characteristics, related to average life time, average volumetric fraction of solid, existing time fraction and frequency of occurrence. The identification of clusters is performed by applying a criterion related to the time average value of the volumetric solid fraction. A qualitative rather than quantitative analysis is performed mainly owing to the unavailability of operational data used in the considered experiments. Concerning qualitative analysis, the simulation results are in good agreement with literature. Some quantitative comparisons between predictions and experiment were also presented to emphasize the capability of the modeling procedure regarding the analysis of macroscopic scale coherent structures. (c) 2007 Elsevier Inc. All rights reserved.

The thickness of the transport layer in stratified geomorphic flows

River Flow 2010

Geomorphic dam-break flows. Part I: conceptual model

Proceedings of the Institution of Civil Engineers - Water Management 163 Issue WM6

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Litoestratigrafia e deformação cenzóica na região de Icapuí, Ceará, e implicações para a estruturação de campos de petróleo na borda ocidental da bacia Potiguar(NE do Brasil)

This thesis deals with the sedimentological/stratigraphic and structural evolution of the sedimentary rocks that occur in the NW continental border of the Potiguar Basin. These rocks are well exposed along coastal cliffs between the localities of Lagoa do Mato and Icapuí, Ceará State (NE Brazil). The sedimentological/stratigraphic study involved, at the outcrop scale, detailed facies descriptions, profile mapping of the vertical succession of different beds, and columnar sections displaying inferred lateral relationships. The approach was complemented by granulometric and petrographic analyses, including the characterization of heavy mineral assemblages. The data set allowed to recognize two kinds of lithological units, a carbonate one of very restricted occurrence at the base of the cliffs, and three younger, distinct siliciclastic units, that predominate along the cliffs, in vertical and lateral extent. The carbonate rocks were correlated to the late Cretaceous Jandaíra Formation, which is covered by the siliciclastic Barreiras Formation. The Barreiras Formation occurs in two distinct structural settings, the usual one with nondeformed, subhorizontal strata, or as tilted beds, affected by strong deformation. Two lithofacies were recognized, vertically arranged or in fault contacts. The lower facies is characterized by silty-argillaceous sandstones with low-angle cross bedding; the upper facies comprises medium to coarse grained sandstones, with conglomeratic layers. The Tibau Formation (medium to coarse-grained sandstones with argillite intercalations) occurs at the NW side of the studied area, laterally interlayered with the Barreiras Formation. Eolic sediments correlated to the Potengi Formation overly the former units, either displaying an angular unconformity, or simply an erosional contact (stratigraphic unconformity). Outstanding structural features, identified in the Barreiras Formation, led to characterize a neocenozoic stress field, which generated faults and folds and/or reactivated older structures in the subjacent late cretaceous (to paleogene, in the offshore basin) section. The structures recognized in the Barreiras Formation comprise two distinct assemblages, namely a main extensional deformation between the localities of Ponta Grossa and Redonda, and a contractional style (succeeded by oblique extensional structures) at Vila Nova. In the first case, the structural assemblage is dominated by N-S (N±20°Az) steep to gently-dipping extensional faults, displaying a domino-style or listric geometry with associated roll-over structures. This deformation pattern is explained by an E-W/WNW extension, contemporaneous with deposition of the upper facies of the Barreiras Formation, during the time interval Miocene to Pleistocene. Strong rotation of blocks and faults generated low-angle distensional faults and, locally, subvertical bedding, allowing to estimate very high strain states, with extension estimates varying between 40% up to 200%. Numerous detachment zones, parallel to bedding, help to acommodate this intense deformation. The detachment surfaces and a large number of faults display mesoscopic features analoguous to the ones of ductile shear zones, with development of S-C fabrics, shear bands, sigmoidal clasts and others, pointing to a hydroplastic deformation regime in these cases. Local occurrences of the Jandaíra limestone are controled by extensional faults that exhume the pre-Barreiras section, including an earlier event with N-S extension. Finally, WNWtrending extensional shear zones and faults are compatible with the Holocene stress field along the present continental margin. In the Vila Nova region, close to Icapuí, gentle normal folds with fold hinges shallowly pluging to SSW affect the lower facies of the Barreiras Formation, displaying an incipient dissolution cleavage associated with an extension lineation at high rake (a S>L fabric). Deposition of the upper facies siliciclastics is controlled by pull-apart graben structures, bordered by N-NE-trending sinistral-normal shear zones and faults, characterizing an structural inversion. Microstructures are compatible with tectonic deformation of the sedimentary pile, burried at shallow depths. The observed features point to high pore fluid pressures during deformation of the sediments, producing hydroplastic structures through mechanisms of granular flow. Such structures are overprinted by microfractures and microfaults (an essentially brittle regime), tracking the change to microfracturing and frictional shear mechanisms accompanying progressive dewatering and sediment lithification. Correlation of the structures observed at the surface with those present at depth was tested through geophysical data (Ground Penetrating Radar, seismics and a magnetic map). EW and NE-trending lineaments are observed in the magnetic map. The seismic sections display several examples of positive flower structures which affect the base of the cretaceous sediments; at higher stratigraphic levels, normal components/slips are compatible with the negative structural inversion characterized at the surface. Such correlations assisted in proposing a structural model compatible with the regional tectonic framework. The strong neogenepleistocene deformation is necessarily propagated in the subsurface, affecting the late cretaceous section (Açu and Jandaíra formations), wich host the hydrocarbon reservoirs in this portion of the Potiguar Basin. The proposed structural model is related to the dextral transcurrent/transform deformation along the Equatorial Margin, associated with transpressive terminations of E-W fault zones, or at their intersections with NE-trending lineaments, such as the Ponta Grossa-Fazenda Belém one (the LPGFB, itself controlled by a Brasiliano-age strike-slip shear zone). In a first step (and possibly during the late Cretaceous to Paleogene), this lineament was activated under a sinistral transpressional regime (antithetic to the main dextral deformation in the E-W zones), giving way to the folds in the lower facies of the Barreiras Formation, as well as the positive flower structures mapped through the seismic sections, at depth. This stage was succeeded (or was penecontemporaneous) by the extensional structures related to a (also sinistral) transtensional movement stage, associated to volcanism (Macau, Messejana) and thermal doming processes during the Neogene-Pleistocene time interval. This structural model has direct implications to hydrocarbon exploration and exploitation activities at this sector of the Potiguar Basin and its offshore continuation. The structure of the reservoirs at depth (Açu Formation sandstones of the post-rift section) may be controlled (or at least, strongly influenced) by the deformation geometry and kinematics characterized at the surface. In addition, the deformation event recognized in the Barreiras Formation has an age close to the one postulated for the oil maturation and migration in the basin, between the Oligocene to the Miocene. In this way, the described structural cenario represents a valid model to understand the conditions of hydrocarbon transport and acummulation through space openings, trap formation and destruction. This model is potentially applicable to the NW region of the Potiguar Basin and other sectors with a similar structural setting, along the brazilian Equatorial Atlantic Margin

The grain-scale distribution and behaviour of melt and fluid in crystalline analogue systems

The production, segregation and migration of melt and aqueous fluids (henceforth called liquid) plays an important role for the transport of mass and energy within the mantle and the crust of the Earth. Many properties of large-scale liquid migration processes such as the permeability of a rock matrix or the initial segregation of newly formed liquid from the host-rock depends on the grain-scale distribution and behaviour of liquid. Although the general mechanisms of liquid distribution at the grain-scale are well understood, the influence of possibly important modifying processes such as static recrystallization, deformation, and chemical disequilibrium on the liquid distribution is not well constrained. For this thesis analogue experiments were used that allowed to investigate the interplay of these different mechanisms in-situ. In high-temperature environments where melts are produced, the grain-scale distribution in “equilibrium” is fully determined by the liquid fraction and the ratio between the solid-solid and the solid-liquid surface energy. The latter is commonly expressed as the dihedral or wetting angle between two grains and the liquid phase (Chapter 2). The interplay of this “equilibrium” liquid distribution with ongoing surface energy driven recrystallization is investigated in Chapter 4 and 5 with experiments using norcamphor plus ethanol liquid. Ethanol in contact with norcamphor forms a wetting angle of about 25°, which is similar to reported angles of rock-forming minerals in contact with silicate melt. The experiments in Chapter 4 show that previously reported disequilibrium features such as trapped liquid lenses, fully-wetted grain boundaries, and large liquid pockets can be explained by the interplay of the liquid with ongoing recrystallization. Closer inspection of dihedral angles in Chapter 5 reveals that the wetting angles are themselves modified by grain coarsening. Ongoing recrystallization constantly moves liquid-filled triple junctions, thereby altering the wetting angles dynamically as a function of the triple junction velocity. A polycrystalline aggregate will therefore always display a range of equilibrium and dynamic wetting angles at raised temperature, rather than a single wetting angle as previously thought. For the deformation experiments partially molten KNO3–LiNO3 experiments were used in addition to norcamphor–ethanol experiments (Chapter 6). Three deformation regimes were observed. At a high bulk liquid fraction >10 vol.% the aggregate deformed by compaction and granular flow. At a “moderate” liquid fraction, the aggregate deformed mainly by grain boundary sliding (GBS) that was localized into conjugate shear zones. At a low liquid fraction, the grains of the aggregate formed a supporting framework that deformed internally by crystal plastic deformation or diffusion creep. Liquid segregation was most efficient during framework deformation, while GBS lead to slow liquid segregation or even liquid dispersion in the deforming areas.