Coupled deformation-flow analysis for methane hydrate production by depressurized wells

A formulation for coupled flow-deformation analysis of methane-hydrate extraction problems is presented. By assuming that the hydrate does not flow, a two phase flow formulation is considered, based on Darcy's law and capillary pressure relation. The formulation is implemented in the finite difference code FLAC. The code was used to investigate the stability of a methane extraction well by depressurizing the well. © 2005 Taylor & Francis Group, London.

Bubble-to-Slug Transition for Low-Re Two-Phase Flows

A critical review on the mechanism and models on the bubble-to-slug transition of two-phase gas-liquid flows are presented in the present paper. It is shown that the most possible mechanism controlled the bubble-to-slug transition is the bubble coalescence. Focusing on the bubble-to-slug transition for the low-Re two-phase flow, a simple Monte Carlo method is used to simulate the influence of the initial bubble size on the bubble-to-slug transition. Some secondary factors, such as the liquid viscosity, the surface tension, and the relative slip between the two phases, are ignored in the present study. It is found that the locus of the dimensionless rate of collision is a universal curve. Based on this curve, it is determined that the bubble initial size can affect the phase distribution and flow pattern when its dimensionless value is in the range from 0.03 to 0.4. A simple relationship between the critical void fraction and the initial bubble size is proposed, which agrees very well with the experimental data.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Multiphase Flow Measurement in the Offshore Oil Industry

In recent years, considerable research has been conducted into the development of a three-phase flowmeter suitable for use in an offshore environment, and oil/gas/water three-phase metering becomes an important aspect in multiphase flow measurement. This paper discusses the importance of three-phase flow measurement in offshore oil industry, describes the current development in this area, and points out the principal strategies which may be used to meter three-phase flow.

Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.

I. The effects of variations in certain parameters upon the efficiency of in vivo hemodialysis with a Kiil dialyzer. II. Interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube.

Part I

Regression analyses are performed on in vivo hemodialysis data for the transfer of creatinine, urea, uric acid and inorganic phosphate to determine the effects of variations in certain parameters on the efficiency of dialysis with a Kiil dialyzer. In calculating the mass transfer rates across the membrane, the effects of cell-plasma mass transfer kinetics are considered. The concept of the effective permeability coefficient for the red cell membrane is introduced to account for these effects. A discussion of the consequences of neglecting cell-plasma kinetics, as has been done to date in the literature, is presented.

A physical model for the Kiil dialyzer is presented in order to calculate the available membrane area for mass transfer, the linear blood and dialysate velocities, and other variables. The equations used to determine the independent variables of the regression analyses are presented. The potential dependent variables in the analyses are discussed.

Regression analyses were carried out considering overall mass-transfer coefficients, dialysances, relative dialysances, and relative permeabilities for each substance as the dependent variables. The independent variables were linear blood velocity, linear dialysate velocity, the pressure difference across the membrane, the elapsed time of dialysis, the blood hematocrit, and the arterial plasma concentrations of each substance transferred. The resulting correlations are tabulated, presented graphically, and discussed. The implications of these correlations are discussed from the viewpoint of a research investigator and from the viewpoint of patient treatment.

Recommendations for further experimental work are presented.

Part II

The interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube in the wavy flow regime has been measured using resistance wave gages. The median water depth, r.m.s. wave height, wave frequency, extrema frequency, and wave velocity have been measured as functions of air and water flow rates. Reynolds numbers, Froude numbers, Weber numbers, and bulk velocities for each phase may be calculated from these measurements. No theory for wave formation and propagation available in the literature was sufficient to describe these results.

The water surface level distribution generally is not adequately represented as a stationary Gaussian process. Five types of deviation from the Gaussian process function were noted in this work. The presence of the tube walls and the relatively large interfacial shear stresses precludes the use of simple statistical analyses to describe the interfacial structure. A detailed study of the behavior of individual fluid elements near the interface may be necessary to describe adequately wavy two-phase flow in systems similar to the one used in this work.

On the transport properties of fluid-particle flow

The hydrodynamic forces acting on a solid particle in a viscous, incompressible fluid medium at low Reynolds number flow is investigated mathematically as a prerequisite to the understanding of transport processes in two-phase flow involving solid particles and fluid. Viscous interaction between a small number of spherical particles and continuous solid boundaries as well as fluid interface are analyzed under a “point-force” approximation. Non-spherical and elastic spherical particles in a simple shear flow area are then considered. Non-steady motion of a spherical particle is briefly touched upon to illustrate the transient effect of particle motion.

A macroscopic continuum description of particle-fluid flow is formulated in terms of spatial averages yielding a set of particle continuum and bulk fluid equations. Phenomenological formulas describing the transport processes in a fluid medium are extended to cases where the volume concentration of solid particles is sufficiently high to exert an important influence. Hydrodynamic forces acting on a spherical solid particle in such a system, e.g. drag, torque, rotational coupling force, and viscous collision force between streams of different sized particles moving relative to each other are obtained. Phenomenological constants, such as the shear viscosity coefficient, and the diffusion coefficient of the bulk fluid, are found as a function of the material properties of the constituents of the two-phase system and the volume concentration of solid. For transient heat conduction phenomena, it is found that the introduction of a complex conductivity for the bulk fluid permits a simple mathematical description of this otherwise complicated process. The rate of heat transfer between particle continuum and bulk fluid is also investigated by means of an Oseen-type approximation to the energy equation.

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.

Poroelastic indentation analysis for hydrated biological tissues

Indentation techniques are employed for the measurement of mechanical properties of a wide range of materials. In particular, techniques focused at small length-scales, such as nanoindentation and AFM indentation, allow for local characterization of material properties in heterogeneous materials including natural tissues and biomimetic materials. Typical elastic analysis for spherical indentation is applicable in the absence of time-dependent deformation, but is inappropriate for materials with time-dependent responses. Recent analyses for the viscoelastic indentation problem, based on elastic-viscoelastic correspondence, have begun to address the issue of time-dependent deformation during an indentation test. The viscoelastic analysis has been shown to fit experimental indentation data well, and has been demonstrated as useful for characterization of viscoelasticity in polymeric materials and in hydrated mineralized tissues. However, a viscoelastic analysis is not necessarily sufficient for multi-phase materials with fluid flow. In the current work, a poroelastic analysis-based on fluid motion through a porous elastic network-is used to examine spherical indentation creep responses of hydrated biological materials. Both analytical and finite element approaches are considered for the poroelastic Hertzian indentation problem. Modeling results are compared with experimental data from nanoindentation of hydrated bone immersed in water and polar solvents (ethanol, methanol, acetone). Baseline (water-immersed) bone responses are characterized using the poroelastic model and numerical results are compared with altered hydration states due to polar solvents. © 2007 Materials Research Society.

Characteristics of blood flow in microchannels and relevant impact on modelling blood behaviour in biochip separators

This paper reports a perspective investigation of computational modelling of blood fluid in microchannel devices as a preparation for future research on fluid-structure interaction (FSI) in biofluid mechanics. The investigation is carried out through two aspects, respectively on physical behaviours of blood flow in microchannels and appropriate methodology for modelling. The physics of blood flow is targeted to the challenges for describing blood flow in microchannels, including rheology of blood fluid, suspension features of red blood cells (RBCs), laminar hydrodynamic influence and effect of surface roughness. The analysis shows that due to the hyperelastic property of RBC and its comparable dimension with microchannels, blood fluid shows complex behaviours of two phase flow. The trajectory and migration of RBCs require accurate description of RBC deformation and interaction with plasma. Following on a discussion of modelling approaches, i.e. Eulerian method and Lagrangian method, the main stream modelling methods for multiphase flow are reviewed and their suitability to blood flow is analysed. It is concluded that the key issue for blood flow modelling is how to describe the suspended blood cells, modelled by Lagrangian method, and couple them with the based flow, modelled by Eulerian method. The multiphase flow methods are thereby classified based on the number of points required for describing a particle, as follows: (i) single-point particle methods, (ii) mutli-point particle methods, (iii) functional particle methods, and (iv) fluid particle methods. While single-point particle methods concentrate on particle dynamic movement, multipoint and functional particle methods can take into account particle mechanics and thus offer more detailed information for individual particles. Fluid particle methods provide good compromise between two phases, but require additional information for particle mechanics. For furthermore detailed description, we suggest to investigate the possibility using two domain coupling method, in which particles and base flow are modelled by two separated solvers. It is expected that this paper could clarify relevant issues in numerical modelling of blood flow in microchannels and induce some considerations for modelling blood flow using multiphase flow methods. © 2012 IEEE.

Unsteady numerical study of wet steam flow in a low pressure steam turbine

In steam power plants condensation already starts in the flow path of the low pressure part of the steam turbine, which leads to a complex three-dimensional two-phase flow. Wetness losses are caused due to thermodynamic and mechanical relaxation processes during condensation and droplet transport. The present investigation focuses on the unsteady effects due to rotor-stator interaction on the droplet formation process. Results of unsteady three dimensional flow simulations of a two-stage steam turbine are presented, whereby this is the first time that non-equilibrium condensation is considered in such simulations. The numerical approach is based on RANS equations, which are extended by a wet steam specific nucleation and droplet growth model. Despite the use of a high performance cluster the unsteady simulation has a considerably high simulation time of approximately 60 days by use of 48 CPUs. © Springer-Verlag Berlin Heidelberg 2012.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Transient flow patterns and bubble slug lengths in parallel microchannels with oxygen gas bubbles produced by catalytic chemical reactions

Transient flow patterns and bubble slug lengths were investigated with oxygen gas (O-2) bubbles produced by catalytic chemical reactions using a high speed camera bonded with a microscope. The microreactor consists of an inlet liquid plenum, nine parallel rectangular microchannels followed by a micronozzle, using the MEMS fabrication technique. The etched surface was deposited by the thin platinum film, which is acted as the catalyst. Experiments were performed with the inlet mass concentration of the hydrogen peroxide from 50% to 90% and the pressure drop across the silicon chip from 2.5 to 20.0 kPa. The silicon chip is directly exposed in the environment thus the heat released via the catalytic chemical reactions is dissipated into the environment and the experiment was performed at the room temperature level. It is found that the two-phase flow with the catalytic chemical reactions display the cyclic behavior. A full cycle consists of a short fresh liquid refilling stage, a liquid decomposition stage followed by the bubble slug flow stage. At the beginning of the bubble slug flow stage, the liquid slug number reaches maximum, while at the end of the bubble slug flow stage the liquid slugs are quickly flushed out of the microchannels. Two or three large bubbles are observed in the inlet liquid plenum, affecting the two-phase distributions in microchannels. The bubble slug lengths, cycle periods as well as the mass flow rates are analyzed with different mass concentrations of hydrogen peroxide and pressure drops. The bubble slug length is helpful for the selection of the future microreactor length ensuring the complete hydrogen peroxide decomposition. Future studies on the temperature effect on the transient two-phase flow with chemical reactions are recommended.

Pressure drops of single and two-phase flows through T-type microchannel mixers

In this paper, preliminary experimental results are presented on pressure drop characteristics of single and two-phase flows through two T-type rectangular microchannel mixers with hydraulic diameters of 528 and 333 mum, respectively. It is shown that both N-2 and water single-phase laminar flows in microchannels, with consideration of experimental uncertainties, are consistent with classic theory, if additional effects, such as entrance effects that will interfere with the interpretation of experimental results, are eliminated by carefully designing the experiments. The obtained pressure drop data of N-2-water two-phase flow in micromixers are analyzed and compared with existing flow pattern-independent models. It is found that the Lockhart-Martinelli method generally underpredicts the frictional pressure drop. Thereafter, a modified correlation of C value in the Chisholm's equation based on linear regression of experimental data is proposed to provide a better prediction of the two-phase frictional pressure drop. Also among the homogeneous flow models investigated, the viscosity correlation of McAdams indicates the best performance in correlating the frictional pressure drop data (mean deviations within +/-20% for two micromixers both). Finally it is suggested that systematic studies are still required to accurately predict two-phase frictional performance in microchannels. (C) 2004 Elsevier B.V. All rights reserved.

ELECTROCATALYTIC OXIDATION AND FLOW DETECTION OF HYDRAZINE COMPOUNDS IN LIQUID-CHROMATOGRAPHY AT A VITAMIN-B-12 ADSORBED GLASSY-CARBON ELECTRODE

A vitamin B-12 chemically modified electrode (CME) was constructed by adsorption of vitamin B-12 onto a glassy carbon surface. The electrode catalyzes the electrooxidation of hydrazine compounds over a wide pH range. The electrocatalytic behavior of hydrazines is elucidated with respect to the CME preparation conditions, solution pH, operating potential, mobile phase flow rate, and other variables. When applied to liquid chromatographic detection of the analytes, the vitamin B-12 CME yielded a linear response range over 2 orders of magnitude, and detection limits at the picomole level. The vitamin B-12 CME offers acceptable catalytic stability in both batch and flow systems.