Development of a PC Tool for Estimating Pressure Drops in Forced-convection Boiler Tubes

Kaksifaasivirtauksen kuvaamiseen käytettävät mallit, ja menetelmät kaksifaasivirtauksen painehäviön määrittämiseksi kehittyvät yhä monimutkaisimmiksi. Höyrystinputkissa tapahtuvien painehäviöiden arvioinnin vaatiman laskennan suorittamiseksi tietokoneohjelman kehittäminen on välttämätöntä. Tässä työssä on kehitetty itsenäinen PC-ohjelma painehäviöiden arvioimiseksi pakotetulle konvektiovirtaukselle pystysuorissa höyrykattilan höyrystinputkissa. Veden ja vesihöyryn aineominaisuuksien laskentaan käytetään IAPWS-IF97 –yhtälökokoelmaa sekä muita tarvittavia IAPWS:n suosittelemia yhtälöitä. Höyrystinputkessa kulloinkin vallitsevan virtausmuodon määrittämiseen käytetään sovelluskelpoisia virtausmuotojen välisiä rajoja kuvaavia yhtälöitä. Ohjelmassa käytetään painehäviön määritykseen kirjallisuudessa julkaistuja yhtälöitä, virtausmuodosta riippuen, alijäähtyneelle virtaukselle, kupla-, tulppa- ja rengasvirtaukselle sekä tulistetun höyryn virtaukselle. Ohjelman laskemia painehäviöarvioita verrattiin kirjallisuudesta valittuihin mittaustuloksiin. Laskettujen painehäviöiden virhe vaihteli välillä –19.5 ja +23.9 %. Virheiden itseisarvojen keskiarvo oli 12.8 %.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

CFD modeling of dispersion water feed in wastewater cleaning application

Fluid particle breakup and coalescence are important phenomena in a number of industrial flow systems. This study deals with a gas-liquid bubbly flow in one wastewater cleaning application. Three-dimensional geometric model of a dispersion water system was created in ANSYS CFD meshing software. Then, numerical study of the system was carried out by means of unsteady simulations performed in ANSYS FLUENT CFD software. Single-phase water flow case was setup to calculate the entire flow field using the RNG k-epsilon turbulence model based on the Reynolds-averaged Navier-Stokes (RANS) equations. Bubbly flow case was based on a computational fluid dynamics - population balance model (CFD-PBM) coupled approach. Bubble breakup and coalescence were considered to determine the evolution of the bubble size distribution. Obtained results are considered as steps toward optimization of the cleaning process and will be analyzed in order to make the process more efficient.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Johtajan jäljillä – Johtaminen varhaiskasvatuksen hajautetuissa organisaatioissa laadun ja pedagogisen tuen näkökulmasta

Tutkimus sijoittuu varhaiskasvatuksen hajautetun organisaation kontekstiin, mutta tulokset ovat siirrettävissä muihinkin suomalaisiin kasvatus- ja opetustoimen organisaatioihin. Hajautettujen organisaatioiden tutkimus on ollut varhaiskasvatuksen kentällä vielä vähäistä, vaikka organisaatiomallin vaikutukset johtajuuden toteuttamiselle ovat merkittävät. Hajautetulla organisaatiolla varhaiskasvatuksessa tarkoitetaan sitä, että yhden johtajan alaisuudessa on monta eri päiväkotia tai erilaisia päivähoitomuotoja. Tämä organisaatiomalli on yhä enenevässä määrin kasvava suomalaisessa varhaiskasvatuksessa. Varhaiskasvatuksen hajautettujen organisaatioiden tutkimuksessa on aiemmin tarkasteltu johtajan ja työntekijöiden ja työntekijöiden keskinäisiä ammatillisia suhteita. Tässä tutkimuksessa näkökulma painottuu johtamiseen ja työskentelyyn hajautetuissa organisaatiossa sinänsä sekä myös laadunarviointiin sekä pedagogiikkaan. Viitekehyksenä tutkimuksessa on LMX-teoria (leader-member-exchange, johtajuuden vaihtoteoria), jossa tarkastellaan esimies-alaissuhdetta ja siihen kiinteästi liittyvää luottamuksen käsitettä. Luottamuksen merkitys hajautetuissa organisaatioissa korostuu, koska esimies ei ole fyysisesti päivittäin läsnä työntekijöiden arjessa. Tutkimuksessa tarkastellaan hajautetuissa varhaiskasvatuksen organisaatioissa työskentelyä seuraavien tutkimuskysymysten avulla: 1) Mitkä ovat varhaiskasvatuksen hajautettujen organisaatioiden johtamisen erityispiirteet? 2) Miten eri työntekijäryhmät kokevat hajautetussa organisaatiossa työskentelyn? 3) Millaisia kokemuksia esimiehillä ja työntekijöillä on heidän yksiköissään toteutetusta laadunarvioinnista? 4) Millaiseksi työntekijät ja esimiehet kokevat esimieheltään saadun tuen? Tutkimuksessa oli kolme eri aineistoa. Ensimmäinen aineisto koostui 11 hajautetun organisaation johtajan haastattelusta. Toinen aineisto (n = 223) sisälsi haastateltujen esimiesten lomakevastausten lisäksi heidän alaisuudessaan toimivien työntekijöiden, 10 esimieskoulutukseen osallistuneen johtajan sekä kolmen erillisyksikön työntekijöiden vastaukset. Kolmas aineisto oli kerätty pääkaupunkiseudulta varhaiskasvatuksen johtajilta lomakekyselynä (n = 112). Aineistoa on analysoitu teorialähtöisen ja aineistolähtöisen sisällönanalyysin ja tilastollisten analyysien avulla Tulokset osoittavat, että johtajat kokivat hallinnollisten töiden vievän paljon aikaa. Esimiehen kanssa eri työpaikassa työskentelevät työntekijät hahmottivat koko organisaation selkeämmin kuin esimiehen kanssa fyysisesti samassa paikassa työskentelevät. Esimiesten käsitysten mukaan laadunarviointia suoritettiin enemmän kuin mitä työntekijöiden mukaan. Työntekijät kaipasivat esimiehiltään tukea yhteistyöhön ja vuorovaikutukseen, pedagogiseen ohjaukseen, kehittämiseen ja toiminnan resursseihin liittyen. Erillisyksikössä työskentelevät kokivat saavansa enemmän tukea kuin esimiehen kanssa fyysisesti samassa yksikössä työskentelevät työntekijät. Sekä esimieheltä saadun pedagogisen tuen että luottamuksen kokemukset kiinnittävät tämän tutkimuksen tulosten mukaan huomion rakenteiden merkitykseen hajautetuissa organisaatioissa. Arviointiin, pedagogiseen tukeen ja tiedonkulkuun liittyvien rakenteiden huomioiminen helpottaa hajautetussa organisaatiossa johtamista. Edellisten lisäksi johtajan selkeä visio omasta johtamistyöstään ja jaetun johtajuuden hyödyntäminen edesauttavat työn hallinnan kokemuksia.

Reconsidérer les interactions entre l’écoulement, le transport de sédiments en charge de fond et la morphologie en rivière à lit de graviers : approches, échelles et analyses

L'objectif ultime en géomorphologie fluviale est d'expliquer les formes des cours d'eau et leur évolution temporelle et spatiale. La multiplication des études nous a mené à la réalisation que les systèmes géomorphologiques sont complexes. Les formes observées sont plus que la somme des processus individuels qui les régissent en raison d’interactions et de rétroactions non-linéaires à de multiples échelles spatiales et temporelles. Dans ce contexte, le but général de la thèse est de proposer et de tester de nouvelles avenues de recherche afin de mieux appréhender la complexité des dynamiques fluviales en utilisant des approches méthodologiques et analytiques mettant l’accent sur les interactions entre l’écoulement, le transport de sédiments en charge fond et la morphologie du lit en rivière graveleuse. Cette orientation découle du constat que les paradigmes actuels en géomorphologie fluviale n’arrivent pas à expliquer adéquatement la variabilité naturelle du transport en charge de fond ainsi que des formes du lit qui en résultent. Cinq pistes de réflexion sont développées sous forme d’articles basés sur des études de cas : 1. L'intégration des échelles de variation de l'écoulement permet d’insérer la notion de structures turbulentes dans des pulsations de plus grande échelle et d'améliorer la compréhension de la variabilité du transport de sédiments. 2. La quantification des taux de changement de l’écoulement (accélération /décélération) au cours d’une crue permet d’expliquer la variabilité des flux de transport en charge fond autant que la magnitude de l’écoulement. 3. L’utilisation de techniques de mesures complémentaires révèle une nouvelle dynamique du lit des rivières graveleuses, la dilatation et la contraction du lit suite à une crue. 4. La remise en cause du fait généralement accepté que le transport en charge de fond est corrélé positivement à l'intensité des modifications morphologiques en raison d’un problème associé aux échelles différentes des processus en cause. 5. L’approche systémique des dynamiques fluviales par l’utilisation d’analyses multivariées permet d’appréhender la complexité des dynamiques de rétroactions linéaires et non-linéaires dans l’évolution d’un chenal et d’illustrer l’importance de l’historique récent des changements géomorphologiques en réponse aux crues. Cette thèse se veut une avancée conceptuelle issue d'une profonde réflexion sur les approches classiques que l'on utilise en géomorphologie fluviale depuis plusieurs décennies. Elle est basée sur un jeu de données unique récolté lors du suivi intensif de 21 évènements de crue dans un petit cours d’eau à lit de graviers, le ruisseau Béard (Québec). Le protocole expérimental axé sur la simultanéité des mesures de l’écoulement, de la morphologie du lit et du transport de sédiments en charge de fond a permis de centrer la recherche directement sur les interactions entre les processus plutôt que sur les processus individuels, une approche rarement utilisée en géomorphologie fluviale. Chacun des chapitres illustre un nouveau concept ou une nouvelle approche permettant de résoudre certaines des impasses rencontrées actuellement en géomorphologie fluviale. Ces travaux ont des implications importantes pour la compréhension de la dynamique des lits de rivières et des habitats fluviaux et servent de point de départ pour de nouveaux développements.

Bubble Behavior in a Taylor Vortex

We present an experimental study on the behavior of bubbles captured in a Taylor vortex. The gap between a rotating inner cylinder and a stationary outer cylinder is filled with a Newtonian mineral oil. Beyond a critical rotation speed (ω[subscript c]), Taylor vortices appear in this system. Small air bubbles are introduced into the gap through a needle connected to a syringe pump. These are then captured in the cores of the vortices (core bubble) and in the outflow regions along the inner cylinder (wall bubble). The flow field is measured with a two-dimensional particle imaging velocimetry (PIV) system. The motion of the bubbles is monitored by using a high speed video camera. It has been found that, if the core bubbles are all of the same size, a bubble ring forms at the center of the vortex such that bubbles are azimuthally uniformly distributed. There is a saturation number (N[subscript s]) of bubbles in the ring, such that the addition of one more bubble leads eventually to a coalescence and a subsequent complicated evolution. Ns increases with increasing rotation speed and decreasing bubble size. For bubbles of non-uniform size, small bubbles and large bubbles in nearly the same orbit can be observed to cross due to their different circulating speeds. The wall bubbles, however, do not become uniformly distributed, but instead form short bubble-chains which might eventually evolve into large bubbles. The motion of droplets and particles in a Taylor vortex was also investigated. As with bubbles, droplets and particles align into a ring structure at low rotation speeds, but the saturation number is much smaller. Moreover, at high rotation speeds, droplets and particles exhibit a characteristic periodic oscillation in the axial, radial and tangential directions due to their inertia. In addition, experiments with non-spherical particles show that they behave rather similarly. This study provides a better understanding of particulate behavior in vortex flow structures.

Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

On the structure of Langmuir turbulence

The Stokes drift induced by surface waves distorts turbulence in the wind-driven mixed layer of the ocean, leading to the development of streamwise vortices, or Langmuir circulations, on a wide range of scales. We investigate the structure of the resulting Langmuir turbulence, and contrast it with the structure of shear turbulence, using rapid distortion theory (RDT) and kinematic simulation of turbulence. Firstly, these linear models show clearly why elongated streamwise vortices are produced in Langmuir turbulence, when Stokes drift tilts and stretches vertical vorticity into horizontal vorticity, whereas elongated streaky structures in streamwise velocity fluctuations (u) are produced in shear turbulence, because there is a cancellation in the streamwise vorticity equation and instead it is vertical vorticity that is amplified. Secondly, we develop scaling arguments, illustrated by analysing data from LES, that indicate that Langmuir turbulence is generated when the deformation of the turbulence by mean shear is much weaker than the deformation by the Stokes drift. These scalings motivate a quantitative RDT model of Langmuir turbulence that accounts for deformation of turbulence by Stokes drift and blocking by the air–sea interface that is shown to yield profiles of the velocity variances in good agreement with LES. The physical picture that emerges, at least in the LES, is as follows. Early in the life cycle of a Langmuir eddy initial turbulent disturbances of vertical vorticity are amplified algebraically by the Stokes drift into elongated streamwise vortices, the Langmuir eddies. The turbulence is thus in a near two-component state, with suppressed and . Near the surface, over a depth of order the integral length scale of the turbulence, the vertical velocity (w) is brought to zero by blocking of the air–sea interface. Since the turbulence is nearly two-component, this vertical energy is transferred into the spanwise fluctuations, considerably enhancing at the interface. After a time of order half the eddy decorrelation time the nonlinear processes, such as distortion by the strain field of the surrounding eddies, arrest the deformation and the Langmuir eddy decays. Presumably, Langmuir turbulence then consists of a statistically steady state of such Langmuir eddies. The analysis then provides a dynamical connection between the flow structures in LES of Langmuir turbulence and the dominant balance between Stokes production and dissipation in the turbulent kinetic energy budget, found by previous authors.

Detection of coherent airstreams using cluster analysis: application to an extratropical cyclone

Flow in geophysical fluids is commonly summarized by coherent streams, for example conveyor belt flows in extratropical cyclones or jet streaks in the upper troposphere. Typically, parcel trajectories are calculated from the flow field and subjective thresholds are used to distinguish coherent streams of interest. This methodology contribution develops a more objective approach to distinguish coherent airstreams within extratropical cyclones. Agglomerative clustering is applied to trajectories along with a method to identify the optimal number of cluster classes. The methodology is applied to trajectories associated with the low-level jets of a well-studied extratropical cyclone. For computational efficiency, a constraint that trajectories must pass through these jet regions is applied prior to clustering; the partitioning into different airstreams is then performed by the agglomerative clustering. It is demonstrated that the methodology can identify the salient flow structures of cyclones: the warm and cold conveyor belts. A test focusing on the airstreams terminating at the tip of the bent-back front further demonstrates the success of the method in that it can distinguish fine-scale flow structure such as descending sting jet airstreams.

Simulação do escoamento multifásico no interior de bombas de cavidades progressivas metálicas

The progressing cavity pumping (PCP) is one of the most applied oil lift methods nowadays in oil extraction due to its ability to pump heavy and high gas fraction flows. The computational modeling of PCPs appears as a tool to help experiments with the pump and therefore, obtain precisely the pump operational variables, contributing to pump s project and field operation otimization in the respectively situation. A computational model for multiphase flow inside a metallic stator PCP which consider the relative motion between rotor and stator was developed in the present work. In such model, the gas-liquid bubbly flow pattern was considered, which is a very common situation in practice. The Eulerian-Eulerian approach, considering the homogeneous and inhomogeneous models, was employed and gas was treated taking into account an ideal gas state. The effects of the different gas volume fractions in pump volumetric eficiency, pressure distribution, power, slippage flow rate and volumetric flow rate were analyzed. The results shown that the developed model is capable of reproducing pump dynamic behaviour under the multiphase flow conditions early performed in experimental works

Numerical simulation of a radial diffuser turbulent airflow

In the present work are presented results from numerical simulations performed with the ANSYS-CFX (R) code. We have studied a radial diffuser flow case, which is the main academic problem used to study the flow behavior on flat plate valves. The radial flow inside the diffuser has important behavior such as the turbulence decay downstream and recirculation regions inside the valve flow channel due to boundary layer detachment. These flow structures are present in compressor reed valve configurations, influencing to a greater extent the compressor efficiency. The main target of the present paper was finding the simulation set-up (computational domain, boundary conditions and turbulence model) that better fits with experimental data published by Tabatabai and Pollard. The local flow turbulence and velocity profiles were investigated using four different turbulence models, two different boundary conditions set-up, two different computational domains and three different flow conditions (Re-in - Reynolds number at the diffuser inlet). We used the Reynolds stress (BSL); the k-epsilon; the RNG k-epsilon; and the shear stress transport (SST) k-omega turbulence models. The performed analysis and comparison of the computational results with experimental data show that the choice of the turbulence model, as well as the choice of the other computational conditions, plays an important role in the results physical quality and accuracy. (c) 2007 Elsevier B.V. All rights reserved.

Informação visual no controle postural de crianças: efeito das características do estímulo

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

Opalas gemológicas do Piauí: gênese revelada por microtermometria e minerais associados

As opalas de Pedro II e Buriti dos Montes, no estado do Piauí, constituem as mais importantes ocorrências brasileiras dessa gema, tanto em termos de volume quanto pela qualidade gemológica, que é comparável à das famosas opalas australianas. No entanto, a informalidade na extração e comercialização destas opalas, assim como a falta de informações quanto à gênese destes depósitos não permitem a prospecção por novas jazidas e o estabelecimento de um certificado de procedência para as opalas do Piauí que permitisse sua inserção formal no mercado gemológico internacional. Alguns autores têm se dedicado ao estudo dessas opalas, revelando fortes evidências de sua origem hidrotermal, mas até então, nenhum trabalho abordou as características físico-químicas dos fluidos que teriam originado esses depósitos de opalas. Diante disso, o principal objetivo deste trabalho foi entender o sistema hidrotermal responsável pela gênese das opalas do Piauí, ou seja, caracterizar os fluidos que originaram a mineralização e mostrar sua relação com o contexto geológico da região. Os municípios de Pedro II e Buriti dos Montes se localizam na porção nordeste do estado do Piauí, a aproximadamente 230 km a leste da capital Teresina, e as ocorrências de opala se encontram na porção basal da Bacia do Parnaíba, constituindo veios e vênulas nos arenitos dos grupos Serra Grande (Buriti dos Montes) e Canindé (Pedro II), os quais são seccionados por soleiras e diques de diabásio da Formação Sardinha. Elas também ocorrem cimentando brechas e como depósitos coluvionares e de paleocanal. Associados às opalas, localmente encontram-se veios de quartzo, calcedônia, barita e hematita (ou goethita). De maneira geral, as opalas de Pedro II apresentam jogo de cores, são predominantemente brancas ou azuladas com aspecto leitoso, semitranslúcidas a opacas e com inclusões sólidas pouco aparentes. Em contrapartida, as opalas de Buriti dos Montes não apresentam jogo de cores, a cor varia entre amarelo claro e vermelho amarronzado, são semitransparentes a translúcidas e contêm grande variedade de inclusões sólidas. Os dados obtidos revelam que as opalas de Pedro II são tipicamente do tipo amorfo (opala-A), enquanto as opalas de Buriti dos Montes variam entre amorfas e cristobalita-tridimita (opala-CT). Na opala preciosa, o típico jogo de cores é causado pelo arranjo regular das esferas de sílica que as constituem. A ausência de cimento opalino entre as esferas reforça a beleza desse efeito. Em contrapartida, as opalas laranja não apresentam jogo de cores, mas têm maior transparência devido ao diminuto tamanho das esferas. As inclusões sólidas também produzem belos efeitos nas opalas estudadas, principalmente na variedade laranja, que é mais transparente. Além disso, o conjunto de inclusões sólidas revela características intrínsecas aos processos hidrotermais que originaram as opalas estudadas. Agregados botrioidais, dendríticos e nodulares são exemplos de inclusões formadas por fragmentos dos arenitos hospedeiros carreados pelos fluidos hidrotermais que geraram as opalas. As inclusões sólidas também têm relação direta com a cor das opalas. Nas opalas de Buriti dos Montes, os tons de vermelho, laranja e amarelo são produzidos pela dissolução parcial das inclusões constituídas por oxihidróxidos de Fe. De maneira semelhante, a cor verde nas opalas preciosas está relacionada aos microcristais de Co-pentlandita inclusos nas mesmas. O conjunto de minerais associados às opalas conduz a uma assinatura mineralógicogeoquímica marcada pelos elevados teores de Fe e Al nas opalas com inclusões de hematita/goethita e caulinita, e assim também com aumento considerável dos teores de elementos terras raras nas opalas em que se concentram as inclusões de caulinita e apatita. Entre os elementos-traço, Ba é o mais abundante, e provavelmente foi incorporado pelo fluido hidrotermal, tendo em vista que veios de barita são encontrados com frequência nessa região da Bacia do Parnaíba. Várias feições como estruturas de fluxo nas opalas, corrosão e dissolução parcial dos cristais de quartzo hialino e de inclusões mineralógicas, vênulas de quartzo hidrotermal sobrecrescidas aos grãos detríticos, e zoneamento dos cristais de quartzo confirmam que essas opalas têm origem hidrotermal. A ruptura do Gondwana teria provocado um vasto magmatismo básico fissural, que por sua vez foi responsável pelo aporte de calor que gerou as primeiras células convectivas de fluidos quentes. A água contida nos arenitos certamente alimentou o sistema e se enriqueceu em sílica através da dissolução parcial ou total dos próprios grãos de quartzo dos arenitos. Este fluido hidrotermal foi posteriormente aprisionado em sistemas de fraturas e nelas se resfriou, precipitando a opala e minerais associados.

Geochemistry of Pleistocene volcanic rocks in the Mariana Forearc

Geophysical surveys of the Mariana forearc, in an area equidistant from the Mariana Trench and the active Mariana Island Arc, revealed a 40-m-deep graben about 13 km northwest of Conical Seamount, a serpentine mud volcano. The graben and its bounding horst blocks are part of a fault zone that strikes northwest-southeast beneath Conical Seamount. One horst block was drilled during Leg 125 of the Ocean Drilling Program (Site 781). Three lithologic units were recovered at Site 781: an upper sedimentary unit, a middle basalt unit, and a lower sedimentary unit. The upper unit, between 0 and 72 mbsf, consists of upper Pliocene to Holocene diatomaceous and radiolarian-bearing silty clay that grades down into vitric silty clay and vitric clayey silt. The middle unit is a Pleistocene vesicular, porphyritic basalt, the top of which corresponds to a high-amplitude reflection on the reflection profiles. The lower unit is a middle to upper (and possibly some lower) Pliocene vitric silty clay and vitric clayey silt similar to the lower part of the upper unit. The thickness of the basalt unit can only be estimated to be between 13 and 25 m because of poor core recovery (28% to 55%). The absence of internal flow structures and the presence of an upper glassy chilled zone and a lower, fine-grained margin suggest that the basalt unit is either a single lava flow or a near-surface sill. The basalt consists of plagioclase phenocrysts with subordinate augite and olivine phenocrysts and of plagioclase-augite-olivine glomerocrysts in a groundmass of plagioclase, augite, olivine, and glass. The basalt is an island arc tholeiite enriched in large-ion-lithophile elements relative to high-field-strength elements, similar to the submarine lavas of the southern arc seamounts. In contrast, volcanic rocks from the active volcanoes on Pagan and Agrigan islands, 100 km to the west of the drill site, are calc-alkaline. The basalt layer, the youngest in-situ igneous layer reported from the Izu-Bonin and Mariana forearcs, is enigmatic because of its location more than 100 km from the active volcanic arc. The sediment layers above and below the basalt unit are late Pliocene in age (about 2.5 Ma) and normally magnetized. The basalt has schlierenlike structures, reverse magnetization, and a K-Ar age of 1.68±0.37 Ma. Thus, the basalt layer is probably a sill fed by magma intruded along a fault zone bounding the horst and graben in the forearc. The geochemistry of the basalt is consistent with a magma source similar to that of the active island arc and from a mantle source above the subducting Pacific plate.