Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

Visualization of turbulent wedges under favourable pressure gradients using both shear-sensitive and temperature-sensitive liquid crystals

Turbulent wedges induced by a 3D surface roughness placed in a laminar boundary layer over a flat plate were visualised for the first time using both shear-sensitive and temperature-sensitive liquid crystals. The experiments were carried out at three different levels of favourable pressure gradients. The purpose of this investigation was to examine the spreading angles of the turbulent wedges indicated by their associated surface shear stresses and heat transfer characteristics and hence obtain further insight about the difference in the behaviour of transitional momentum and thermal boundary layers when a streamwise pressure gradient exists. It was shown that under a zero pressure gradient the spreading angles indicated by the two types of liquid crystals are the same, but the difference increases as the level of favourable pressure gradient increases. The result from the present study could have an important implication to the transition modelling of thermal boundary layers over gas turbine blades.

Intermittency transport modeling of separated flow transition

An intermittency transport model is proposed for modeling separated-flow transition. The model is based on earlier work on prediction of attached flow bypass transition and is applied for the first time to model transition in a separation bubble at various degrees of free-stream turbulence. The model has been developed so that it takes into account the entrainment of the surrounding fluid. Experimental investigations suggest that it is this phenomena which ultimately determines the extent of the separation bubble. Transition onset is determined via a boundary layer correlation based on momentum thickness at the point of separation. The intermittent flow characteristic of the transition process is modeled via an intermittency transport equation. This accounts for both normal and streamwise variation of intermittency and hence models the entrainment of surrounding flow in a more accurate manner than alternative prescribed intermittency models. The model has been validated against the well established T3L semicircular leading edge flat plate test case for three different degrees of free-stream turbulence characteristic of turbomachinery blade applications.

Micro-ramp control for oblique shock wave / boundary layer interactions

Supersonic engine intakes operating supercritically feature shock wave / boundary layer interactions (SBLIs), which are conventionally controlled using boundary layer bleed. The momentum loss of bleed flow causes high drag, compromising intake performance. Micro-ramp sub-boundary layer vortex generators (SBVGs) have been proposed as an alternative form of flow control for oblique SBLIs in order to reduce the bleed requirement. Experiments have been conducted at Mach 2.5 to characterise the flow details on such devices and investigate their ability to control the interaction between an oblique shock wave and the naturally grown turbulent boundary layer on the tunnel floor. Micro-ramps of four sizes with heights ranging from 25% to 75% of the uncontrolled boundary layer thickness were tested. The flow over all sizes of microramp was found to be similar, featuring streamwise counter-rotating vortices which entrain high momentum fluid, locally reducing the boundary layer displacement thickness. When installed ahead of the shock interaction it was found that the positioning of the micro-ramps is of limited importance. Micro-ramps did not eliminate flow separation. However, the previously two-dimensional separation was broken up into periodic three-dimensional separation zones. The interaction length was reduced and the pressure gradient across the interaction was increased.

Compressor wake/leading-edge interactions at off design incidences

In this paper, the effects of wake/leading-edge interactions were studied at off-design conditions. Measurements were performed on the stator-blade suction surface at midspan. The leading-edge flow-field was investigated using hotwire micro-traverses, hotfilm surface shear-stress sensors and pressure micro-tappings. The trailing-edge flow-field was investigated using hotwire boundary-layer traverses. Unsteady CFD calculations were also performed to aid the interpretation of the results. At low flow coefficients, the time-averaged momentum thickness of the leading-edge boundary layer was found to rise as the flow coefficient was reduced. The time-resolved momentum-thickness rose due to the interaction of the incoming rotor wake. As the flow coefficient was reduced, the incoming wakes increased in pitch-wise extent, velocity deficit and turbulence intensity. This increased both the time-resolved rise in the momentum thickness and the turbulent spot production within the wake affected boundary-layer. Close to stall, a drop in the leading-edge momentum thickness was observed in-between wake events. This was associated with the formation of a leading-edge separation bubble in-between wake events. The wake interaction with the bubble gave rise to a shedding phenomenon, which produced large length scale disturbances in the surface shear stress. Copyright © 2008 by ASME.

Heat budget studies of the north Arabian Sea during summer and winter seasons, 1992

In this study heat budget components and momentum flux for August and January 1992 over the north Arabian Sea are computed. The marine meteorological data measured on board during the cruises of PAK-US joint project (NASEER) are used for the computation. Significant differences were found in the heat budget components as well as in the momentum flux during different monsoon periods over the north Arabian Sea. The latent heat flux was always positive and attributed to the large vapour pressure gradient. The computed moisture and latent heat fluxes in January were higher than August The highest value of latent heat flux 309 W/m2 at station 8 was evaluated. These higher latent heat fluxes were due to the large vapour pressure gradient, air-sea temperature difference, the wind speed, and the prevailing wind direction (from north and northeast). Negative values of sensible heat fluxes in both seasons indicate that the heat transfer was from the atmosphere to the ocean. The negative values of net heat gain indicate that the sea surface field became an energy sink: or the sea surface supplied more energy to the atmosphere than it received from it. Large variation in the momentum flux mainly attributed to the variation in the wind speed. Aerial averages of heat and momentum fluxes were also computed.

On computing the physical sources of jet noise

We derive a closed system of equations that relates the acoustically radiating flow variables to the sources of sound for homentropic flows. We use radiating density, momentum density and modified pressure as the dependent variables which leads to simple source terms for the momentum equations. The source terms involve the non-radiating parts of the density and momentum density fields. These non-radiating components are obtained by removing the radiating wavenumbers in the Fourier domain. We demonstrate the usefulness of this new technique on an axi-symmetric jet solution of the Navier-Stokes equations, obtained by direct numerical simulation (DNS). The dominant source term is proportional to the square of the non-radiating part of the axial momentum density. We compare the sound sources to that obtained by an acoustic analogy and find that they have more realistic physical properties. Their frequency content and amplitudes are consistent with. We validate the sources by computing the radiating sound field and comparing it to the DNS solution. © 2010 by S. Sinayoko, A. Agarwal.

Mandibular morphometric variation among Chinese cercopithecoids and the unique structure of the snub-nosed monkey (Rhinopithecus) mandible

In order to understand how mandibular structure differs among the Chinese cercopithecoids (Rhinopithecus, Trachypithecus and Macaca), particularly the uniqueness of the snub-nosed monkeys (Rhinopithecus), we analysed ten mandibular measurements by principal components analysis (PCA), and examined scaling patterns. The results provided by the PCA illustrated differences due to size among the cercopithecoids and the relationship between colobines (Trachypithecus and Rhinopithecus) and cercopithecines, in which macaques (Macaca) are included. Allometric analysis indicated that, biomechanically, there is not a marked difference between macaques and leaf-eating monkeys. This may be associated with the fact that both share some similar ecology and niches in south and southwest China. The snub-nosed monkeys exhibit a significantly more robust mandible, evident in the symphysis, corpus, condyle, and masticatory momentum arm. This supports the hypothesis, based on the study of dental structure, that Rhinopithecus is a unique group in Asian Old World monkeys (OWMs) and has developed some unique characteristics in order to adapt to the tough food available in the severe cold climate of the Plateaux of Qinghai-Tibet, Yun-Gui and Qingling in China.

COMPRESSOR BLADE

A blade (2) for a compressor comprising a pressure surface (6) and a suction surface (4), wherein the suction surface (4) comprises: a discontinuity in the chord-wise curvature of an intermediate portion of the suction surface between the blade leading edge and blade trailing edge; and a slot (20) arranged along at least a portion of the blade (2) in a substantially span-wise direction, the slot (2) being disposed at or near the said change in curvature, such that when in use lower momentum fluid near the suction surface (4) is aspirated into the slot (20).

Frictional flow of damp granular material in a conical centrifuge

An analysis is given of velocity and pressure-dependent sliding flow of a thin layer of damp granular material in a spinning cone. Integral momentum equations for steady state, axisymmetric flow are derived using a boundary layer approximation. These reduce to two coupled first-order differential equations for the radial and circumferential sliding velocities. The influence of viscosity and friction coefficients and inlet boundary conditions is explored by presentation of a range of numerical results. In the absence of any interfacial shear traction the flow would, with increasing radial and circumferential slip, follow a trajectory from inlet according to conservation of angular momentum and kinetic energy. Increasing viscosity or friction reduces circumferential slip and, in general, increases the residence time of a particle in the cone. The residence time is practically insensitive to the inlet velocity. However, if the cone angle is very close to the friction angle then the residence time is extremely sensitive to the relative magnitude of these angles. © 2011 Authors.

A model for diffusive and displacive phase transitions: Thermo-chemo-mechanical coupling effects

A diffuse interface phase field model is proposed for the unified analysis of diffusive and displacive phase transitions under nonisothermal conditions. Two order parameters are used for the description of the phenomena: one is related to the solute mass fraction and the other to the strain. The model governing equations come from the balance of linear momentum, the solute mass balance (which will lead to the Cahn-Hilliard equation) and the balance of internal energy. Thermodynamic restrictions allow to define constitutive relations for the thermodynamic forces and for the mechanical and chemical dissipations. Numerical tests carried out at different values of the initial temperature show that the model is able to describe the main features of both the displacive and the diffusive phase transitions, as well as their effect on the temperature. © 2010, Advanced Engineering Solutions.

Discrete element calculations of the impact of a sand column against rigid structures

Discrete particle simulations of column of an aggregate of identical particles impacting a rigid, fixed target and a rigid, movable target are presented with the aim to understand the interaction of an aggregate of particles upon a structure. In most cases the column of particles is constrained against lateral expansion. The pressure exerted by the particles upon the fixed target (and the momentum transferred) is independent of the co-efficient of restitution and friction co-efficient between the particles but are strongly dependent upon the relative density of the particles in the column. There is a mild dependence on the contact stiffness between the particles which controls the elastic deformation of the densified aggregate of particles. In contrast, the momentum transfer to a movable target is strongly sensitive to the mass ratio of column to target. The impact event can be viewed as an inelastic collision between the sand column and the target with an effective co-efficient of restitution between 0 and 0.35 depending upon the relative density of the column. We present a foam analogy where impact of the aggregate of particles can be modelled by the impact of an equivalent foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally and yet give predictions to within 5% of the full discrete particle calculations. They also suggest that "model" materials can be used to simulate the loading by an aggregate of particles within a laboratory setting. © 2012 Elsevier Ltd. All rights reserved.

Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

The deposition of small particles from a turbulent air flow in a curved duct

The paper describes an experimental and theoretical study of the deposition of small spherical particles from a turbulent air flow in a curved duct. The objective was to investigate the interaction between the streamline curvature of the primary flow and the turbulent deposition mechanisms of diffusion and turbophoresis. The experiments were conducted with particles of uranine (used as a fluorescent tracer) produced by an aerosol generator. The particles were entrained in an air flow which passed vertically downwards through a long straight channel of rectangular cross-section leading to a 90° bend. The inside surfaces of the channel and bend were covered with tape to collect the deposited particles. Following a test run the tape was removed in sections, the uranine was dissolved in sodium hydroxide solution and the deposition rates established by measuring the uranine concentration with a luminescence spectrometer. The experimental results were compared with calculations of particle deposition in a curved duct using a computer program that solved the ensemble-averaged particle mass and momentum conservation equations. A particle density-weighted averaging procedure was used and the equations were expressed in terms of the particle convective, rather than total, velocity. This approach provided a simpler formulation of the particle turbulence correlations generated by the averaging process. The computer program was used to investigate the distance required to achieve a fully-developed particle flow in the straight entry channel as well as the variation of the deposition rate around the bend. The simulations showed good agreement with the experimental results. © 2012 Elsevier Ltd.