Análise e modelagem do escoamento transitório de fluidos de perfuração

Drilling fluids present a thixotropic behavior and they usually gel when at rest. The sol-gel transition is fundamental to prevent the deposit of rock fragments, generated by drilling the well, over the drill bit during eventual stops. Under those conditions, high pressures are then required in order to break-up the gel when circulation is resumed. Moreover, very high pressures can damage the rock formation at the bottom of the well. Thus, a better understanding of thixotropy and the behavior of thixotropic materials becomes increasingly important for process control. The mechanisms that control thixotropy are not yet well defined and modeling is still a challenge. The objective of this work is to develop a mathematical model to study the pressure transmission in drilling fluids. This work presents a review of thixotropy and of different mathematical models found in the literature that are used to predict such characteristic. It also shows a review of transient flows of compressible fluids. The problem is modeled as the flow between the drillpipe and the annular region (space between the wall and the external part of the drillpipe). The equations that describe the problem (mass conservation, momentum balance, constitutive and state) are then discretized and numerically solved by using a computational algorithm in Fortran. The model is validated with experimental and numerical data obtained from the literature. Comparisons between experimental data obtained from Petrobras and calculated by three viscoplastic and one pseudoplastic models are conducted. The viscoplastic fluids, due to the yield stress, do not fully transmit the pressure to the outlet of the annular space. Sensibility analyses are then conducted in order to evaluate the thixotropic effect in pressure transmission.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Analytical solution for transient one-dimensional couette flow considering constant and time-dependent pressure gradients

This paperaims to determine the velocity profile, in transient state, for a parallel incompressible flow known as Couette flow. The Navier-Stokes equations were applied upon this flow. Analytical solutions, based in Fourier series and integral transforms, were obtained for the one-dimensional transient Couette flow, taking into account constant and time-dependent pressure gradients acting on the fluid since the same instant when the plate starts it´s movement. Taking advantage of the orthogonality and superposition properties solutions were foundfor both considered cases. Considering a time-dependent pressure gradient, it was found a general solution for the Couette flow for a particular time function. It was found that the solution for a time-dependent pressure gradient includes the solutions for a zero pressure gradient and for a constant pressure gradient.

Transient raise of endothelin-1 plasma level and reduction of ocular blood flow in a patient with optic neuritis

PURPOSE: To analyze how far an ischemic component might have been involved in optic neuritis. METHODS: Case report: a 32-year-old man with symptoms characteristic for optic neuritis underwent extensive clinical, laboratory/serological and vascular examination for systemic associations and vascular involvement. RESULTS: The patient was found to have a temporary ocular blood flow dysregulation and increased plasma endothelin-1 levels which decreased after the acute phase of the optic nerve. CONCLUSIONS: We conclude that there might be an ischemic component in this patient with optic neuritis and hypothesize that this ischemic component is at least in part due to a temporarily increased endothelin-1 level.

Transient viscous flow around an elliptic cylinder /

Originally presented as the author's thesis (M.S.), University of Illinois at Urbana-Champaign.

Interim report on transient heat-transfer measurements of catalytic recombination in a step-function flow of atomic oxygen.

"AFOSR 3203. Contract no. AF 49(638)-782."

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Linear and sliding-mode control design for matrix converter-based unified power flow controllers

This paper presents the design and compares the performance of linear, decoupled and direct power controllers (DPC) for three-phase matrix converters operating as unified power flow controllers (UPFC). A simplified steady-state model of the matrix converter-based UPFC fitted with a modified Venturini high-frequency pulse width modulator is first used to design the linear controllers for the transmission line active (P) and reactive (Q) powers. In order to minimize the resulting cross coupling between P and Q power controllers, decoupled linear controllers (DLC) are synthesized using inverse dynamics linearization. DPC are then developed using sliding-mode control techniques, in order to guarantee both robustness and decoupled control. The designed P and Q power controllers are compared using simulations and experimental results. Linear controllers show acceptable steady-state behaviour but still exhibit coupling between P and Q powers in transient operation. DLC are free from cross coupling but are parameter sensitive. Results obtained by DPC show decoupled power control with zero error tracking and faster responses with no overshoot and no steady-state error. All the designed controllers were implemented using the same digital signal processing hardware.