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.

Langmuir turbulence and deeply penetrating jets in an unstratified mixed layer

The influence of surface waves and an applied wind stress is studied in an ensemble of large eddy simulations to investigate the nature of deeply penetrating jets into an unstratified mixed layer. The influence of a steady monochromatic surface wave propagating parallel to the wind direction is parameterized using the wave-filtered Craik-Leibovich equations. Tracer trajectories and instantaneous downwelling velocities reveal classic counterrotating Langmuir rolls. The associated downwelling jets penetrate to depths in excess of the wave's Stokes depth scale, δs. Qualitative evidence suggests the depth of the jets is controlled by the Ekman depth scale. Analysis of turbulent kinetic energy (tke) budgets reveals a dynamical distinction between Langmuir turbulence and shear-driven turbulence. In the former, tke production is dominated by Stokes shear and a vertical flux term transports tke to a depth where it is dissipated. In the latter, tke production is from the mean shear and is locally balanced by dissipation. We define the turbulent Langmuir number Lat = (v*/Us)0.5 (v* is the ocean's friction velocity and Us is the surface Stokes drift velocity) and a turbulent anisotropy coefficient Rt = /( + ). The transition between shear-driven and Langmuir turbulence is investigated by varying external wave parameters δs and Lat and by diagnosing Rt and the Eulerian mean and Stokes shears. When either Lat or δs are sufficiently small the Stokes shear dominates the mean shear and the flow is preconditioned to Langmuir turbulence and the associated deeply penetrating jets.

How to save a bad element with weak boundary conditions

The P-1-P-1 finite element pair is known to allow the existence of spurious pressure (surface elevation) modes for the shallow water equations and to be unstable for mixed formulations. We show that this behavior is strongly influenced by the strong or the weak enforcement of the impermeability boundary conditions. A numerical analysis of the Stommel model is performed for both P-1-P-1 and P-1(NC)-P-1 mixed formulations. Steady and transient test cases are considered. We observe that the P-1-P-1 element exhibits stable discrete solutions with weak boundary conditions or with fully unstructured meshes. (c) 2005 Elsevier Ltd. All rights reserved.

The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model (HadGEM1). Part II: Aspects of variability and regional climate

The performance of the atmospheric component of the new Hadley Centre Global Environmental Model (HadGEM1) is assessed in terms of its ability to represent a selection of key aspects of variability in the Tropics and extratropics. These include midlatitude storm tracks and blocking activity, synoptic variability over Europe, and the North Atlantic Oscillation together with tropical convection, the Madden-Julian oscillation, and the Asian summer monsoon. Comparisons with the previous model, the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3), demonstrate that there has been a considerable increase in the transient eddy kinetic energy (EKE), bringing HadGEM1 into closer agreement with current reanalyses. This increase in EKE results from the increased horizontal resolution and, in combination with the improved physical parameterizations, leads to improvements in the representation of Northern Hemisphere storm tracks and blocking. The simulation of synoptic weather regimes over Europe is also greatly improved compared to HadCM3, again due to both increased resolution and other model developments. The variability of convection in the equatorial region is generally stronger and closer to observations than in HadCM3. There is, however, still limited convective variance coincident with several of the observed equatorial wave modes. Simulation of the Madden-Julian oscillation is improved in HadGEM1: both the activity and interannual variability are increased and the eastward propagation, although slower than observed, is much better simulated. While some aspects of the climatology of the Asian summer monsoon are improved in HadGEM1, the upper-level winds are too weak and the simulation of precipitation deteriorates. The dominant modes of monsoon interannual variability are similar in the two models, although in HadCM3 this is linked to SST forcing, while in HadGEM1 internal variability dominates. Overall, analysis of the phenomena considered here indicates that HadGEM1 performs well and, in many important respects, improves upon HadCM3. Together with the improved representation of the mean climate, this improvement in the simulation of atmospheric variability suggests that HadGEM1 provides a sound basis for future studies of climate and climate change.

An optical study of multiple NEIAL events driven by low energy electron precipitation

Optical data are compared with EISCAT radar observations of multiple Naturally Enhanced Ion-Acoustic Line (NEIAL) events in the dayside cusp. This study uses narrow field of view cameras to observe small-scale, short-lived auroral features. Using multiple-wavelength optical observations, a direct link between NEIAL occurrences and low energy (about 100 eV) optical emissions is shown. This is consistent with the Langmuir wave decay interpretation of NEIALs being driven by streams of low-energy electrons. Modelling work connected with this study shows that, for the measured ionospheric conditions and precipitation characteristics, growth of unstable Langmuir (electron plasma) waves can occur, which decay into ion-acoustic wave modes. The link with low energy optical emissions shown here, will enable future studies of the shape, extent, lifetime, grouping and motions of NEIALs.

Scattering of electromagnetic waves by rough interfaces and inhomogeneous layers

We consider a two-dimensional problem of scattering of a time-harmonic electromagnetic plane wave by an infinite inhomogeneous conducting or dielectric layer at the interface between semi-infinite homogeneous dielectric half-spaces. The magnetic permeability is assumed to be a fixed positive constant. The material properties of the media are characterized completely by an index of refraction, which is a bounded measurable function in the layer and takes positive constant values above and below the layer, corresponding to the homogeneous dielectric media. In this paper, we examine only the transverse magnetic (TM) polarization case. A radiation condition appropriate for scattering by infinite rough surfaces is introduced, a generalization of the Rayleigh expansion condition for diffraction gratings. With the help of the radiation condition the problem is reformulated as an equivalent mixed system of boundary and domain integral equations, consisting of second-kind integral equations over the layer and interfaces within the layer. Assumptions on the variation of the index of refraction in the layer are then imposed which prove to be sufficient, together with the radiation condition, to prove uniqueness of solution and nonexistence of guided wave modes. Recent, general results on the solvability of systems of second kind integral equations on unbounded domains establish existence of solution and continuous dependence in a weighted norm of the solution on the given data. The results obtained apply to the case of scattering by a rough interface between two dielectric media and to many other practical configurations.

Extracting planetary waves from geomagnetic time series using Empirical Mode Decomposition

Empirical Mode Decomposition is presented as an alternative to traditional analysis methods to decompose geomagnetic time series into spectral components. Important comments on the algorithm and its variations will be given. Using this technique, planetary wave modes of 5-, 10-, and 16-day mean periods can be extracted from magnetic field components of three different stations in Germany. In a second step, the amplitude modulation functions of these wave modes can be shown to contain significant contribution from solar cycle variation through correlation with smoothed sunspot numbers. Additionally, the data indicate connections with geomagnetic jerk occurrences, supported by a second set of data providing reconstructed near-Earth magnetic field for 150 years. Usually attributed to internal dynamo processes within the Earth's outer core, the question of who is impacting whom will be briefly discussed here.

Interaction of equatorial waves through resonance with the diurnal cycle of tropical heating

In this work, a new theoretical mechanism is presented in which equatorial Rossby and inertio-gravity wave modes may interact with each other through resonance with the diurnal cycle of tropical deep convection. We have adopted the two-layer incompressible equatorial primitive equations forced by a parametric heating that roughly represents deep convection activity in the tropical atmosphere. The heat source was parametrized in the simplest way according to the hypothesis that it is proportional to the lower-troposphere moisture convergence, with the background moisture state function mimicking the structure of the ITCZ. In this context, we have investigated the possibility of resonant interaction between equatorially trapped Rossby and inertio-gravity modes through the diurnal cycle of the background moisture state function. The reduced dynamics of a single resonant duo shows that when this diurnal variation is considered, a Rossby wave mode can undergo significant amplitude modulations when interacting with an inertio-gravity wave mode, which is not possible in the context of the resonant triad non-linear interaction. Therefore, the results suggest that the diurnal variation of the ITCZ can be a possible dynamical mechanism that leads the Rossby waves to be significantly affected by high frequency modes.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Sediment deformation in Miocene and post-Miocene strata, Northeastern Brazil: Evidence for paleoseismicity in a passive margin

Stratigraphic intervals characterized by varied and complex styles of soft-sediment deformation structures are well preserved in Miocene and Late Pleistocene to Holocene deposits of a sedimentary basin located in Northeastern Brazil. The Miocene strata, represented by the Barreiras Formation, record only brittle structures, including numerous faults and fractures with straight and high angle-dipping planes that are often filled with sands derived from overlying beds. Folds consisting of broad anticlines and synclines are also present in this unit. The late Pleistocene to Holocene deposits, named Post-Barreiras Sediments, contain an indurated sandy package with a large variety of ductile and brittle deformation structures (i.e., massive sandstones with isolated sand fragments and breccias, undulatory strata, sand dykes and diapirs, sinks and bowls, pebbly pockets, plunged sediment mixtures, fitted sand masses, cone-shaped cracks, fault grading and sedimentary enclaves). These features, confined to sharp-based stratigraphic horizons that progressively grade downward into undisturbed deposits, are related to seismic shocks of high surface-wave magnitude (i.e., Ms>5 or 6). Amalgamated seismites suggest that previously formed seismites were affected by subsequent seismic-wave propagation. Seismic waves caused by activity along one, or most likely, several tectonic structures would have propagated throughout the depositional environment, producing laterally extensive seismites. The close proximity to earthquake epicenters would have promoted pervasive re-sedimentation due to pore overpressure, resulting high volumes of massive sandstones and breccia. The similarity between deposits with correlatable strata from many other areas along the Brazilian coast allows raise the hypothesis that the seismic episodes might have affected sedimentation patterns in a large (i.e., extension of several hundreds of kilometers) geographic area. Thus, the modern seismicity recorded along Northeastern Brazil was recurrent during the Quaternary and, perhaps, also in the Pliocene. The estimated high magnitude of the seismic events and the great regional extent of the affected area demonstrate that the Brazilian coast experienced tectonic stress through the last geological episodes of its evolution, which would have favored sediment accumulation and penecontemporaneous re-sedimentation. This geological context is unexpected in a passive margin, inducing to revisit the debate on how active is a passive margin. (C) 2010 Elsevier B.V. All rights reserved.

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN(+) and CNN(-): a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN(+) and CNN(-) were investigated at a high level of theory. Very accurate structural parameters for the states X(3)Sigma(-) and A(3)Pi of CNN, and X(2)Pi of both CNN(+) and CNN(-) were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Pi states. At the UCCSD(T)/CBS(T-5) level of theory, Delta(f)H(0 K) = 138.89 kcal/mol and Delta(f)H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS(T-5) approach for the states X(3)Sigma(-) A(3)Pi, and a(1)Delta of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b(1)Sigma(+), c(1)Pi, d(1)Sigma(-), and B(3)Sigma(-), and excitation energies (T(e)) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X (3)Sigma(-) <- E (3)Pi.

Raman and photoluminescence of Er3+-doped SnO2 obtained via the sol-gel technique from solutions with distinct pH

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Arcabouço geofísico, isostasia e causas do magmatismo cenozóico da província Borborema e de sua margem continental (Nordeste do Brasil)

The Borborema Province (BP) is a geologic domain located in Northeastern Brazil. The BP is limited at the south by the São Francisco craton, at the west by the Parnaíba basin, and both at the north and east by coastal sedimentary basins. Nonetheless the BP surface geology is well known, several key aspects of its evolution are still open, notably: i)its tectonic compartmentalization established after the Brasiliano orogenesis, ii) the architecture of its cretaceous continental margin, iii) the elastic properties of its lithosphere, and iv) the causes of magmatism and uplifting which occurred in the Cenozoic. In this thesis, a regional coverage of geophysical data (elevation, gravity, magnetic, geoid height, and surface wave global tomography) were integrated with surface geologic information aiming to attain a better understanding of the above questions. In the Riacho do Pontal belt and in the western sector of the Sergipano belt, the neoproterozoic suture of the collision of the Sul domain of the BP with the Sanfranciscana plate (SFP) is correlated with an expressive dipolar gravity anomaly. The positive lobule of this anomaly is due to the BP lower continental crust uplifting whilst the negative lobule is due to the supracrustal nappes overthrusting the SFP. In the eastern sector of the Sergipano belt, this dipolar gravity anomaly does not exist. However the suture still can be identified at the southern sector of the Marancó complex arc, alongside of the Porto da Folha shear zone, where the SFP N-S geophysical alignments are truncated. The boundary associated to the collision of the Ceará domain of the BP with the West African craton is also correlated with a dipolar gravity anomaly. The positive lobule of this anomaly coincides with the Sobral-Pedro II shear zone whilst the negative lobule is associated with the Santa Quitéria magmatic arc. Judging by their geophysical signatures, the major BP internal boundaries are: i)the western sector of the Pernambuco shear zone and the eastern continuation of this shear zone as the Congo shear zone, ii) the Patos shear zone, and iii) the Jaguaribe shear zone and its southwestern continuation as the Tatajuba shear zone. These boundaries divide the BP in five tectonic domains in the geophysical criteria: Sul, Transversal, Rio Grande do Norte, Ceará, and Médio Coreaú. The Sul domain is characterized by geophysical signatures associated with the BP and SFP collision. The fact that Congo shear zone is now proposed as part of the Transversal domain boundary implies an important change in the original definition of this domain. The Rio Grande do Norte domain presents a highly magnetized crust resulted from the superposition of precambrian and phanerozoic events. The Ceará domain is divided by the Senador Pompeu shear zone in two subdomains: the eastern one corresponds to the Orós-Jaguaribe belt and the western one to the Ceará-Central subdomain. The latter subdomain exhibits a positive ENE-W SW gravity anomaly which was associated to a crustal discontinuity. This discontinuity would have acted as a rampart against to the N-S Brasiliano orogenic nappes. The Médio Coreaú domain also presents a dipolar gravity anomaly. Its positive lobule is due to granulitic rocks whereas the negative one is caused by supracrustal rocks. The boundary between Médio Coreaú and Ceará domains can be traced below the Parnaíba basin sediments by its geophysical signature. The joint analysis of free air anomalies, free air admittances, and effective elastic thickness estimates (Te) revealed that the Brazilian East and Equatorial continental margins have quite different elastic properties. In the first one 10 km < Te < 20 km whereas in the second one Te ≤ 10 km. The weakness of the Equatorial margin lithosphere was caused by the cenozoic magmatism. The BP continental margin presents segmentations; some of them have inheritance from precambrian structures and domains. The segmentations conform markedly with some sedimentary basin features which are below described from south to north. The limit between Sergipe and Alagoas subbasins coincides with the suture between BP and SFP. Te estimates indicates concordantly that in Sergipe subbasin Te is around 20 km while Alagoas subbasin has Te around 10 km, thus revealing that the lithosphere in the Sergipe subbasin has a greater rigidity than the lithosphere in the Alagoas subbasin. Additionally inside the crust beneath Sergipe subbasin occurs a very dense body (underplating or crustal heritage?) which is not present in the crust beneath Alagoas subbasin. The continental margin of the Pernambuco basin (15 < Te < 25 km) presents a very distinct free air edge effect displaying two anomalies. This fact indicates the existence in the Pernambuco plateau of a relatively thick crust. In the Paraíba basin the free air edge effect is quite uniform, Te ≈ 15 km, and the lower crust is abnormally dense probably due to its alteration by a magmatic underplating in the Cenozoic. The Potiguar basin segmentation in three parts was corroborated by the Te estimates: in the Potiguar rift Te ≅ 5 km, in the Aracati platform Te ≅ 25 km, and in the Touros platform Te ≅ 10 km. The observed weakness of the lithosphere in the Potiguar rift segment is due to the high heat flux while the relatively high strength of the lithosphere in the Touros platform may be due to the existence of an archaean crust. The Ceará basin, in the region of Mundaú and Icaraí subbasins, presents a quite uniform free air edge effect and Te ranges from 10 to 15 km. The analysis of the Bouguer admittance revealed that isostasy in BP can be explained with an isostatic model where combined surface and buried loadings are present. The estimated ratio of the buried loading relative to the surface loading is equal to 15. In addition, the lower crust in BP is abnormally dense. These affirmations are particularly adequate to the northern portion of BP where adherence of the observed data to the isostatic model is quite good. Using the same above described isostatic model to calculate the coherence function, it was obtained that a single Te estimate for the entire BP must be lower than 60 km; in addition, the BP north portion has Te around 20 km. Using the conventional elastic flexural model to isostasy, an inversion of crust thickness was performed. It was identified two regions in BP where the crust is thickened: one below the Borborema plateau (associated to an uplifting in the Cenozoic) and the other one in the Ceará domain beneath the Santa Quitéria magmatic arc (a residue associated to the Brasiliano orogenesis). On the other hand, along the Cariri-Potiguar trend, the crust is thinned due to an aborted rifting in the Cretaceous. Based on the interpretation of free air anomalies, it was inferred the existence of a large magmatism in the oceanic crust surrounding the BP, in contrast with the incipient magmatism in the continent as shown by surface geology. In BP a quite important positive geoid anomaly exists. This anomaly is spatially correlated with the Borborema plateau and the Macaú-Queimadas volcanic lineament. The integrated interpretation of geoid height anomaly data, global shear velocity model, and geologic data allow to propose that and Edge Driven Convection (EDC) may have caused the Cenozoic magmatism. The EDC is an instability that presumably occurs at the boundary between thick stable lithosphere and oceanic thin lithosphere. In the BP lithosphere, the EDC mechanism would have dragged the cold lithospheric mantle into the hot asthenospheric mantle thus causing a positive density contrast that would have generated the main component of the geoid height anomaly. In addition, the compatibility of the gravity data with the isostatic model, where combined surface and buried loadings are present, together with the temporal correlation between the Cenozoic magmatism and the Borborema plateau uplifting allow to propose that this uplifting would have been caused by the buoyancy effect of a crustal root generated by a magmatic underplating in the Cenozoic

Theory of small aspect ratio waves in deep water

In the limit of small values of the aspect ratio parameter (or wave steepness) which measures the amplitude of a surface wave in units of its wave-length, a model equation is derived from the Euler system in infinite depth (deep water) without potential flow assumption. The resulting equation is shown to sustain periodic waves which on the one side tend to the proper linear limit at small amplitudes, on the other side possess a threshold amplitude where wave crest peaking is achieved. An explicit expression of the crest angle at wave breaking is found in terms of the wave velocity. By numerical simulations, stable soliton-like solutions (experiencing elastic interactions) propagate in a given velocities range on the edge of which they tend to the peakon solution. (c) 2005 Elsevier B.V. All rights reserved.

Sea waves generated by extratropical cyclones in the South Atlantic Ocean: Hindcast and validation against altimeter data

This work presents analyses of the atmospheric conditions and the hindcast of the surface wave field when six extratropical cyclones formed and displaced over the South Atlantic Ocean (10degreesN, 60degreesS; 75degreesW, 15degreesE) between April and September 1999. These events caused high sea waves associated with hazardous conditions along the south and southeast coast of Brazil. The meteorological composite fields for these cyclones show a strong near-surface wind velocity (up to 14 m s(-1)) during its mature phase. The sea-state wave hindcast was obtained using a third-generation wave model forced by the 10-m above ground level wind field from the National Centers for Environmental Prediction-National Center for Atmospheric Research reanalysis dataset. Closer to the south and southeast Brazilian coast, the hindcast results showed significant wave heights of up to 5 m in some of the events. The wave hindcast results for the significant wave height were compared against satellite altimeter data at 6 h intervals. The statistical index showed a systematic underestimation of the significant wave height by 0.5 m. The correlation between wave hindcast and altimeter measurements was greater than 90%, showing a good phase reproduction by the wave model.