Análise da evolução termomecânica da Bacia do Amazonas

Os dados geológicos e geofísicos escolhidos para o tema de estudo pertencem a Bacia do Amazonas, na região centro-norte do Brasil. A Bacia do Amazonas é uma bacia intracratônica com cerca de 500.000 km. A mesma está limitada ao norte pelo Escudo das Guianas e ao sul pelo Escudo Brasileiro. O limite oeste com a Bacia do Solimões é marcado pelo Arco de Purus, ao passo que o Arco de Gurupá constitui seu limite leste. Possui características inerentes a uma bacia intracratônica paleozóica, com uma longa história evolutiva, marcada por discordâncias expressivas e com uma cunha sedimentar relativamente rasa se comparada às bacias cretáceas brasileiras, o que levanta controvérsia a respeito da suficiência do soterramento para a eficiência de geração de hidrocarboneto. Podem ser reconhecidas nos 5000 m do preenchimento sedimentar da Bacia do Amazonas, duas seqüências de primeira ordem: uma paleozóica, intrudida por diques e soleiras de diabásio, na passagem do Triássico para o Jurássico, e uma mesozóica-cenozóica que representam um aspecto importante na evolução térmica da matéria orgânica que ocorre na primeira seqüência. Com relação à exploração de petróleo, apesar do fomento exploratório ocorrido nos últimos anos, a bacia ainda é considerada pouco explorada sendo sua maior reserva a da província de Urucu. Um dos fatores que dificultam bastante a exploração desta bacia assim como a bacia do Solimões a oeste é o acesso restrito, pois estão situadas em áreas remotas e florestadas, de difícil acesso, com muitas reservas indígenas e florestais, o que causa restrições logísticas, operacionais e legais. O efeito térmico das intrusões ígneas é considerado como o responsável pelo acréscimo de calor necessário à maturação da matéria orgânica e conseqüente geração de hidrocarbonetos. Este trabalho contribui com a reconstrução da história térmica desta bacia a partir da modelagem das variáveis termais e da história de soterramento. Para isso, foram utilizados modelos consagrados na literatura, que permitem, de forma simples, a estimativa do fluxo térmico através do embasamento e da seqüência sedimentar. Na análise da influência de intrusões ígneas na estrutura térmica da bacia, o modelo bidimensional desenvolvido pelo método de diferenças finitas se mostrou apropriado. Utilizou-se o fluxo térmico basal calculado nas condições de contorno da modelagem da influência térmica das ígneas. Como resultado obteve-se a estruturação térmica da bacia e a historia maturação de suas rochas geradoras

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.

HEAVY VEHICLE VIBRATION - A CASE STUDY.

A method is presented for the digital simulation of multiple degrees-of-freedom lumped parameter vibrating systems with arbitrary constitutive elements in an inertial frame of reference. The geometry of the system is treated independently of the constitutive elements and as a result nonlinear (time domain) or linearised (frequency domain) calculations may be performed using a single input description. The method is used to simulate a 3-axle rigid heavy commercial vehicle for harsh vibrating conditions. Some of the assumptions to which the calculations are sensitive are examined. Agreement between the response of a 3-dimensional whole vehicle model and measurements on the test vehicle is satisfactory.

A note on application of hills procedure for growth-curve analysis of Indian major carp, Labeo rohita (Ham.), during the gonodal development period at Shillong, Meghalaya

The present study deals with the length increment data of 15 adult Labeo rohita (Ham.) over a period of five months by the applicatin of finite difference method at an altitude of 1496 m above mean sea level at Shilllong, Meghalaya.

Adjoint linearised Euler solver in the frequency domain for jet noise modelling

The paper is devoted to extending the new efficient frequency-domain method of adjoint Green's function calculation to curvilinear multi-block RANS domains for middle and farfield sound computations. Numerical details of the method such as grids, boundary conditions and convergence acceleration are discussed. Two acoustic source models are considered in conjunction with the method and acoustic modelling results are presented for a benchmark low-Reynolds-number jet case.

A comparison of methods for the coupled analysis of floating structures

The dynamic analysis of a deepwater floating platform and the associated mooring/riser system should ideally be fully coupled to ensure a reliable response prediction. It is generally held that a time domain analysis is the only means of capturing the various coupling and nonlinear effects accurately. However, in recent work it has been found that for an ultra-deepwater floating system (2000m water depth), the highly efficient frequency domain approach can provide highly accurate response predictions. One reason for this is the accuracy of the drag linearization procedure over both first and second order motions, another reason is the minimal geometric nonlinearity displayed by the mooring lines in deepwater. In this paper, the aim is to develop an efficient analysis method for intermediate water depths, where both mooring/vessel coupling and geometric nonlinearity are of importance. It is found that the standard frequency domain approach is not so accurate for this case and two alternative methods are investigated. In the first, an enhanced frequency domain approach is adopted, in which line nonlinearities are linearized in a systematic way. In the second, a hybrid approach is adopted in which the low frequency motion is solved in the time domain while the high frequency motion is solved in the frequency domain; the two analyses are coupled by the fact that (i) the low frequency motion affects the mooring line geometry for the high frequency motion, and (ii) the high frequency motion affects the drag forces which damp the low frequency motion. The accuracy and efficiency of each of the methods are systematically compared. Copyright © 2007 by ASME.

Monitoring tunnel deformation induced by close-proximity bored tunneling using distributed optical fiber strain measurements

The paper describes the use of optical fiber Brillouin Optical Time Domain Reflectometry (BOTDR) to monitor the strain distribution in an existing tunnel while a twin tunnel was bored at close-proximity. The twin circular bored tunnels between Serangoon and Bartley stations on the new Circle Line Stage 3 subway in Singapore were constructed at close-proximity to avoid underpinning the foundations of adjacent buildings. The minimum clear separation of the two tunnels is 2.3m (0.4 times the tunnel diameter). The Outer Tunnel was constructed first, followed by the Inner Tunnel, with the earth-pressure balance tunnel boring machines maintained at a minimum of 100m apart. In this trial application of BOTDR, the strain distribution along the Outer Tunnel was measured, in order to monitor its deformation due to the boring of the Inner Tunnel at close-proximity. The aim of the trial application was to determine the practicality of this monitoring method for future use in 'live' tunnels. This paper compares the measurements obtained from optical fiber BOTDR with conventional methods of tunnel monitoring and describes preliminary installation and workmanship guidelines derived from lessons learnt during this trial. © 2007 ASCE.

Fibre optics monitoring of clay cuttings and embankments along London's ring motorway

Although a wide range of techniques exist for slope monitoring, the task of monitoring slopes is sometimes complicated by the extensive nature and unpredictability of slope movements. The Brillouin optical time-domain reflectometer (BOTDR) is a distributed optical fiber strain measurement technology utilising Brillouin scattering. This method measures continuous strain along a standard optical fibre over a distance up to 10 km and hence has potential to detect deformations and diagnose problems along large sections of slopes and embankments. This paper reports the demonstration of BOTDR method for monitoring surface ground movements of clay cuttings and embankments along London's ring M25 motorway. A field trial investigating varying methods of onsite fibre optic installations was conducted. The surrounding ground was artificially moved by excavating a 3 m deep trench perpendicular to the instrumented sections. Results obtained from onsite installations after slope movement demonstrate a half-pipe covered fibre optic installed on wide (200mm) Tensar ™SS20 geogrid gives the most consistent recorded strain change profile. Initial conclusions suggest this method best represents induced ground motion at the surface and hence is recommended for implementation in future sitework. Copyright ASCE 2008.

The use of spatial impulse responses to characterise flexible forming processes with mobile tools

A novel test method for the characterisation of flexible forming processes is proposed and applied to four flexible forming processes: Incremental Sheet Forming (ISF), conventional spinning, the English wheel and power hammer. The proposed method is developed in analogy with time-domain control engineering, where a system is characterised by its impulse response. The spatial impulse response is used to characterise the change in workpiece deformation created by a process, but has also been applied with a strain spectrogram, as a novel way to characterise a process and the physical effect it has on the workpiece. Physical and numerical trials to study the effects of process and material parameters on spatial impulse response lead to three main conclusions. Incremental sheet forming is particularly sensitive to process parameters. The English wheel and power hammer are strongly similar and largely insensitive to both process and material parameters. Spinning develops in two stages and is sensitive to most process parameters, but insensitive to prior deformation. Finally, the proposed method could be applied to modelling, classification of existing and novel processes, product-process matching and closed-loop control of flexible forming processes. © 2012 Elsevier B.V.

Prediction of tread block forces for a free-rolling tyre in contact with a rough road

This paper describes a new approach to model the forces on a tread block for a free-rolling tyre in contact with a rough road. A theoretical analysis based on realistic tread mechanical properties and road roughness is presented, indicating partial contact between a tread block and a rough road. Hence an asperity-scale indentation model is developed using a semi-empirical formulation, taking into account both the rubber viscoelasticity and the tread block geometry. The model aims to capture the essential details of the contact at the simplest level, to make it suitable as part of a time-domain dynamic analysis of the coupled tyre-road system. The indentation model is found to have a good correlation with the finite element (FE) predictions and is validated against experimental results using a rolling contact rig. When coupled to a deformed tyre belt profile, the indentation model predicts normal and tangential force histories inside the tyre contact patch that show good agreement with FE predictions. © 2012 Elsevier B.V..

Robust F0 estimation based on a multi-microphone periodicity function for distant-talking speech

This work addresses the problem of deriving F0 from distanttalking speech signals acquired by a microphone network. The method here proposed exploits the redundancy across the channels by jointly processing the different signals. To this purpose, a multi-microphone periodicity function is derived from the magnitude spectrum of all the channels. This function allows to estimate F0 reliably, even under reverberant conditions, without the need of any post-processing or smoothing technique. Experiments, conducted on real data, showed that the proposed frequency-domain algorithm is more suitable than other time-domain based ones.

Stiffness and strength of tridimensional periodic lattices

This paper presents a method for the linear analysis of the stiffness and strength of open and closed cell lattices with arbitrary topology. The method hinges on a multiscale approach that separates the analysis of the lattice in two scales. At the macroscopic level, the lattice is considered as a uniform material; at the microscopic scale, on the other hand, the cell microstructure is modelled in detail by means of an in-house finite element solver. The method allows determine the macroscopic stiffness, the internal forces in the edges and walls of the lattice, as well as the global periodic buckling loads, along with their buckling modes. Four cube-based lattices and nine cell topologies derived by Archimedean polyhedra are studied. Several of them are characterized here for the first time with a particular attention on the role that the cell wall plays on the stiffness and strength properties. The method, automated in a computational routine, has been used to develop material property charts that help to gain insight into the performance of the lattices under investigation. © 2012 Elsevier B.V.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Shock dynamics of random structures.

Predicting the response of a structure following an impact is of interest in situations where parts of a complex assembly may come into contact. Standard approaches are based on the knowledge of the impulse response function, requiring the knowledge of the modes and the natural frequencies of the structure. In real engineering structures the statistics of higher natural frequencies follows those of the Gaussian Orthogonal Ensemble, this allows the application of random point process theory to get a mean impulse response function by the knowledge of the modal density of the structure. An ensemble averaged time history for both the response and the impact force can be predicted. Once the impact characteristics are known in the time domain, a simple Fourier Transform allows the frequency range of the impact excitation to be calculated. Experimental and numerical results for beams, plates, and cylinders are presented to confirm the validity of the method.

On the time delay in continuous in-cylinder sampling from IC engines

When gas sample is continuously drawn from the cylinder of an internal combustion engine, the sample that appears at the end of the sampling system corresponds to the in-cylinder content sometime ago because of the finite transit time which is a function of the cylinder pressure history. This variable delay causes a dispersion of the sample signal and makes the interpretation of the signal difficult An unsteady flow analysis of a typical sampling system was carried out for selected engine loads and speeds. For typical engine operation, a window in which the delay is approximately constant may be found. This window gets smaller with increase in engine speed, with decrease in load, and with the increase in exit pressure of the sampling system.