Vortex-induced vibrations of a rigid cylinder on elastic supports with endstops. Part 1: Experimental Results

This paper describes an experimental investigation into the effect of restricting the vortex-induced vibrations of a spring-mounted rigid cylinder by means of stiff mechanical endstops. Cases of both asymmetric and symmetric restraint are investigated. Results show that limiting the amplitude of the vibrations strongly affects the dynamics of the cylinder, particularly when the offset is small. Fluid-structure interaction is profoundly affected, and the well-known modes of vortex shedding observed with a linear elastic system are modified or absent. There is no evidence of lock-in, and the dominant impact frequency corresponds to a constant Strouhal number of 0.18. The presence of an endstop on one side of the motion can lead to large increases in displacements in the opposite direction. Attention is also given to the nature of the developing chaotic motion, and to impact velocities, which in single-sided impacts approach the maximum velocity of a cylinder with linear compliance undergoing VIV at lock-in. With symmetrical endstops, impact velocities were about one-half of this. Lift coefficients are computed from an analysis of the cylinder’s motion between impacts.

Nearshore oscillating wave surge converters and the development of Oyster

Oscillating wave surge converters (OWSCs) are a class of wave power technology that exploits the enhanced horizontal fluid particle movement of waves in the nearshore coastal zone with water depths of 10–20 m. OWSCs predominantly oscillate horizontally in surge as opposed to the majority of wave devices, which oscillate vertically in heave and usually are deployed in deeper water. The characteristics of the nearshore wave resource are described along with the hydrodynamics of OWSCs. The variables in the OWSC design space are discussed together with a presentation of some of their effects on capture width, frequency bandwidth response and power take-off characteristics. There are notable differences between the different OWSCs under development worldwide, and these are highlighted. The final section of the paper describes Aquamarine Power’s 315kW Oyster 1 prototype, which was deployed at the European Marine Energy Centre in August 2009. Its place in the OWSC design space is described along with the practical experience gained. This has led to the design of Oyster 2, which was deployed in August 2011. It is concluded that nearshore OWSCs are serious contenders in the mix of wave power technologies. The nearshore wave climate has a narrower directional spread than the offshore, the largest waves are filtered out and the exploitable resource is typically only 10–20% less in 10m depth compared with 50m depth. Regarding the devices, a key conclusion is that OWSCs such as Oyster primarily respond in the working frequency range to the horizontal fluid acceleration; Oyster is not a drag device responding to horizontal fluid velocity. The hydrodynamics of Oyster is dominated by inertia with added inertia being a very significant contributor. It is unlikely that individual flap modules will exceed 1MW in installed capacity owing to wave resource, hydrodynamic and economic constraints. Generating stations will be made up of line arrays of flaps with communal secondary power conversion every 5–10 units.

How does Oyster work? The simple interpretation of Oyster mathematics

Oyster® is a surface-piercing flap-type device designed to harvest wave energy in the nearshore environment. Established mathematical theories of wave energy conversion, such as 3D point-absorber and 2D terminator theory, are inadequate to accurately describe the behaviour of Oyster, historically resulting in distorted conclusions regarding the potential of such a concept to harness the power of ocean waves. Accurately reproducing the dynamics of Oyster requires the introduction of a new reference mathematical model, the “flap-type absorber”. A flap-type absorber is a large thin device which extracts energy by pitching about a horizontal axis parallel to the ocean bottom. This paper unravels the mathematics of Oyster as a flap-type absorber. The main goals of this work are to provide a simple–yet accurate–physical interpretation of the laws governing the mechanism of wave power absorption by Oyster and to emphasise why some other, more established, mathematical theories cannot be expected to accurately describe its behaviour.

An 40Ar/39Ar, Rb/Sr, and stable isotope study of micas in low-grade fold-and-thrust belt: An example from the Swiss Helvetic Alps

White micas in carbonate-rich tectonites and a few other rock types of large thrusts in the Swiss Helvetic fold-and-thrust belt have been analyzed by Ar-40/Ar-39 and Rb/Sr techniques to better constrain the timing of Alpine deformation for this region. Incremental Ar-40/Ar-39 heating experiments of 25 weakly metamorphosed (anchizone to low greenschist) samples yield plateau and staircase spectra. We interpret most of the staircase release spectra result from variable mixtures of syntectonic (neoformed) and detrital micas. The range in dates obtained within individual spectra depends primarily on the duration of mica nucleation and growth, and relative proportions of neoformed and detrital mica. Rb/Sr analyses of 12 samples yield dates of ca. 10-39 Ma (excluding one anomalously young sample). These dates are slightly younger than the Ar-40/Ar-39 total gas dates obtained for the same samples. The Rb/ Sr dates were calculated using initial Sr-87/Sr-86 ratios obtained from the carbonate-dominated host rocks, which are higher than normal Mesozoic carbonate values due to exchange with fluids of higher Sr-87/Sr-86 ratios (and lower O-18/O-16 ratios). Model dates calculated using Sr-87/Sr-86 values typical of Mesozoic marine carbonates more closely approximate the Ar-40/Ar-39 total gas dates for most of the samples. The similarities of Rb/Sr and Ar-40/Ar-39 total gas dates are consistent with limited amounts of detrital mica in the samples. The delta(18)O values range from 24-15%. (VSMOW) for 2-6 mum micas and 27-16parts per thousand for the carbonate host rocks. The carbonate values are significantly lower than their protolith values due to localized fluid-rock interaction and fluid flow along most thrust surfaces. Although most calcite-mica pairs are not in oxygen isotope equilibrium at temperatures of ca. 200-400 degreesC, their isotopic fractionations are indicative of either 1) partial exchange between the minerals and a common external fluid, or 2) growth or isotopic exchange of the mica with the carbonate after the carbonate had isotopically exchanged with an external fluid. The geological significance of these results is not easily or uniquely determined, and exemplifies the difficulties inherent in dating very fine-grained micas of highly deformed tectonites in low-grade metamorphic terranes. Two generalizations can be made regarding the dates obtained from the Helvetic thrusts: 1) samples from the two highest thrusts (Mt. Gond and Sublage) have all of their Ar-40/Ar-39 steps above 20 Ma, and 2) most samples from the deepest Helvetic thrusts have steps (often accounting for more than 80% of Ar-39 release) between 15 and 25 Ma. These dates are consistent with the order of thrusting in the foreland-imbricating system and increase proportions of neoformed to detrital mica in the more metamorphosed hinterland and deeply buried portions of the nappe pile. Individual thrusts accommodated the majority of their displacement during their initial incorporation into the foreland-imbricating system, and some thrusts remained active or were reactivated down to 15 Ma.

Une nouvelle mise en oeuvre de la méthode IIM pour les équations de Navier-Stokes en présence d'une force singulière

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

Simulation de la nage anguilliforme

Ce document traite premièrement des diverses tentatives de modélisation et de simulation de la nage anguilliforme puis élabore une nouvelle technique, basée sur la méthode de la frontière immergée généralisée et la théorie des poutres de Reissner-Simo. Cette dernière, comme les équations des fluides polaires, est dérivée de la mécanique des milieux continus puis les équations obtenues sont discrétisées afin de les amener à une résolution numérique. Pour la première fois, la théorie des schémas de Runge-Kutta additifs est combinée à celle des schémas de Runge-Kutta-Munthe-Kaas pour engendrer une méthode d’ordre de convergence formel arbitraire. De plus, les opérations d’interpolation et d’étalement sont traitées d’un nouveau point de vue qui suggère l’usage des splines interpolatoires nodales en lieu et place des fonctions d’étalement traditionnelles. Enfin, de nombreuses vérifications numériques sont faites avant de considérer les simulations de la nage.

Investigations on Structural response of Water Backed Perforated plates

The study envisaged herein contains the numerical investigations on Perforated Plate (PP) as well as numerical and experimental investigations on Perforated Plate with Lining (PPL) which has a variety of applications in underwater engineering especially related to defence applications. Finite element method has been adopted as the tool for analysis of PP and PPL. The commercial software ANSYS has been used for static and free vibration response evaluation, whereas ANSYS LS-DYNA has been used for shock analysis. SHELL63, SHELL93, SOLID45, SOLSH190, BEAM188 and FLUID30 finite elements available in the ANSYS library as well as SHELL193 and SOLID194 available in the ANSYS LS-DYNA library have been made use of. Unit cell of the PP and PPL which is a miniature of the original plate with 16 perforations have been used. Based upon the convergence characteristics, the utility of SHELL63 element for the analysis of PP and PPL, and the required mesh density are brought out. The effect of perforation, geometry and orientation of perforation, boundary conditions and lining plate are investigated for various configurations. Stress concentration and deflection factor are also studied. Based on these investigations, stadium geometry perforation with horizontal orientation is recommended for further analysis.Linear and nonlinear static analysis of PP and PPL subjected to unit normal pressure has been carried out besides the free vibration analysis. Shock analysis has also been carried out on these structural components. The analytical model measures 0.9m x 0.9m with stiffener of 0.3m interval. The influence of finite element, boundary conditions, and lining plate on linear static response has been estimated and presented. Comparison of behavior of PP and PPL in the nonlinear strain regime has been made using geometric nonlinear analysis. Free vibration analysis of the PP and PPL has been carried out ‘in vacuum’ condition and in water backed condition, and the influence of water backed condition and effect of perforation on natural frequency have been investigated.Based upon the studies on the vibration characteristics of NPP, PP and PPL in water backed condition and ‘in vacuum’ condition, the reduction in the natural frequency of the plate in immersed condition has been rightly brought out. The necessity to introduce the effect of water medium in the analysis of water backed underwater structure has been highlighted.Shock analysis of PP and PPL for three explosives viz., PEK, TNT and C4 has been carried out and deflection and stresses on plate as well as free field pressure have been estimated using ANSYS LS-DYNA. The effect of perforations and the effect of lining plate have been predicted. Experimental investigations of the measurement of free field pressure using PPL have been conducted in a shock tank. Free field pressure has been measured and has been validated with finite element analysis results. Besides, an experiment has been carried out on PPL, for the comparison of the static deflection predicted by finite element analysis.The distribution of the free field pressure and the estimation of differential pressure from experimentation and the provision for treating the differential pressure as the resistance, as a part of the design load for PPL, has been brought out.

Numerical Methods for the Unsteady Compressible Navier-Stokes Equations

We consider numerical methods for the compressible time dependent Navier-Stokes equations, discussing the spatial discretization by Finite Volume and Discontinuous Galerkin methods, the time integration by time adaptive implicit Runge-Kutta and Rosenbrock methods and the solution of the appearing nonlinear and linear equations systems by preconditioned Jacobian-Free Newton-Krylov, as well as Multigrid methods. As applications, thermal Fluid structure interaction and other unsteady flow problems are considered. The text is aimed at both mathematicians and engineers.

Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional- and shear-wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse-echo method. The measurements were made both in vacuum-dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin-impregnated polished thin sections, X-ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional- and shear-wave velocities (V-p and V-s, respectively) decrease with increasing porosity and that V-p decreases approximately twice as fast as V-s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore-structure-dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot-Gassmann velocity values are greater than the measured velocity values due to the rock-fluid interaction. This is not accounted for in the Biot-Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time-average relationships overestimated the measured velocities even more than the Biot model.

HYDROTHERMAL REE FLUOROCARBONATE MINERALIZATION AT BARRA DO ITAPIRAPUA, A MULTIPLE STOCKWORK CARBONATITE, SOUTHERN BRAZIL

The Barra do Itapirapua ( BIT) carbonatites in southern Brazil belong to the final stages of the Early Cretaceous alkaline rock - carbonatite magmatism of the Ponta Grossa Arch Province. The BIT complex is a dyke and vein stockwork in which four main carbonatitic phases are recognized, mainly magnesiocarbonatites and ferrocarbonatites. These carbonatites are generally overprinted by pervasive hydrothermal events. The C-O stable isotopic data indicate re-equilibration under hydrothermal conditions at temperatures between 375 and 80 degrees C. Significant amounts of REE fluorocarbonate minerals, relatively Sr- and Th-rich, were deposited. Syntaxy between synchysite-(Ce) and parisite-(Ce) is very common owing to the similarity in structures, with alternating (001) layers of (CeF), (CO3) and (Ca). However, bastnasite-(Ce) occurs as individual crystals, overgrown by the synchysite and parisite polycrystals. Textural and chemical reactions between the REE fluorocarbonates provide insights into the mobility of rare-earth elements during fluid-rock interaction. The BIT complex is considered to be of potential economic interest for production of the rare-earth concentrates.

Unraveling an antique subduction process from metamorphic basement around Medellin city, Central Cordillera of Colombian Andes

In the surroundings of Caldas and El Retiro cities (Colombia) metamorphic rocks derived from basic and pelitic protoliths comprise the Caldas amphibole schist and the Ancon schist respectively. Subordinated metamorphosed granite bodies (La Miel gneiss) are associated to these units, and The El Retiro amphibolites, migmatites and granulites crops out eastwards of these units, separated by shear zones. The Caldas amphibole schist and the Ancon schist protoliths could have been formed in a distal marine reduced environment and amalgamated to the South American continent in an apparent Triassic subduction event. The El Retiro rocks are akin to a continental basement and possible include impure metasediments of continental margin, whose metamorphism originated granulite facies rocks and migmatites as a result of the anatexis of quartz-feldspathic rocks. The metamorphism was accompanied by intense deformation, which has juxtaposed both migmatites and granulite blocks. Afterward, heat and fluid circulation associated with the emplacement of minor igneous intrusions resulted in intense fluid-rock interaction, variations in the grain size of the minerals and, especially, intense retrograde metamorphic re-equilibrium. Thermobarometric estimations for the Caldas amphibole schist indicate metamorphism in the Barrovian amphibolite fades. The metamorphic path is counter-clockwise, but retrograde evolution could not be precisely defined. The pressures of the metamorphism in these rocks range from 6.3 to 13.5 kbar, with narrow temperature ranging from 550 to 630 degrees C. For the Ancon schist metapelites the P-T path is also counter-clockwise, with a temperature increase evidenced by the occurrence of sillimanite and the cooling by later kyanite. The progressive metamorphism event occurred at pressures of 7.6-7.2 kbar and temperatures of 645-635 degrees C for one sample and temperature between 500 and 600 degrees C under constant pressure of 6 kbar. The temperature estimated for these rocks varies between 400 and 555 degrees C at pressures of 5-6 kbar in the retrograde metamorphic path. The El Retiro rocks evidence strong decompression with narrow variation in temperature, showing pressure values between 8.7 and 2.7 kbar at temperatures of 740-633 degrees C. These metamorphic fragments of the basement in the Central Cordillera of the Colombian Andes could represent a close relationship with an antique subduction zone. (C) 2011 Elsevier Ltd. All rights reserved.

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

Cubic interactions and quantum criticality in dimerized antiferromagnets

In certain Mott-insulating dimerized antiferromagnets, triplet excitations of the paramagnetic phase display both three-particle and four-particle interactions. When such a magnet undergoes a quantum phase transition into a magnetically ordered state, the three-particle interaction becomes part of the critical theory provided that the lattice ordering wave vector is zero. One microscopic example is the staggered-dimer antiferromagnet on the square lattice, for which deviations from O(3) universality have been reported in numerical studies. Using both symmetry arguments and microscopic calculations, we show that a nontrivial cubic term arises in the relevant order-parameter quantum field theory, and we assess its consequences using a combination of analytical and numerical methods. We also present finite-temperature quantum Monte Carlo data for the staggered-dimer antiferromagnet which complement recently published results. The data can be consistently interpreted in terms of critical exponents identical to that of the standard O(3) universality class, but with anomalously large corrections to scaling. We argue that the cubic interaction of critical triplons, although irrelevant in two spatial dimensions, is responsible for the leading corrections to scaling due to its small scaling dimension.

Micro-raman spectroscopic characterization of a CR-39 detector

Characterization by micro-Raman spectroscopy of polymeric materials used as nuclear track detectors reveals physico-chemical and morphological information on the material's molecular structure. In this work, the nuclear track detector poly(allyl diglycol carbonate), or Columbia Resin 39 (CR-39), was characterized according to the fluence of alpha particles produced by a 226Ra source and chemical etching time. Therefore, damage of the CR-39 chemical structure due to the alpha-particle interaction with the detector was analyzed at the molecular level. It was observed that the ionization and molecular excitation of the CR-39 after the irradiation process entail cleavage of chemical bonds and formation of latent track. In addition, after the chemical etching, there is also loss of polymer structure, leading to the decrease of the group density C-O-C (∼888 cm-1), CH=CH (∼960 cm -1), C-O (∼1110 cm-1), C-O-C (∼1240 cm -1), C-O (∼1290 cm-1), C-O (∼1741 cm -1), -CH2- (∼2910 cm-1), and the main band -CH2- (∼2950 cm-1). The analyses performed after irradiation and chemical etching led to a better understanding of the CR-39 molecular structure and better comprehension of the process of the formation of the track, which is related to chemical etching kinetics. Copyright © 2013 Society for Applied Spectroscopy.