Measurements of Flow Patterns and Velocities of Gas -Water Two -Phase Flows Using Electrical Resistance Imaging

This paper presents a measurement of flow patterns and flow velocities of gas-water two-phase flows based on the technique of electrical resistance tomography (ERT) in a 40m horizontal flow loop. A single-plane and dual-plane ERT sensor on conductive ring technique were used to gather sufficient information for the implementation of flow characteristics particularly flow pattern recognition and air cavity velocity measurement. A fast data collection strategy was applied to the dual-plane ERT sensor and an iterative algorithm was used for image reconstruction. Results, in respect to flow patterns and velocity maps, are reported.

Effect of Pipe Diameter on Flow Pattern Transition in Vertical Oil Water Flows

This paper investiges the effect of pipe diameter on flow pattern transition boundary in oil water vertical flows, and proposes a model to determine the maximum inner diameter (D_{infty s}) of a pipe in which the slug flow would not occur When pipe inner diameter D>D_{infty s}, only bubble flow exists, while D

Distributed Load Control in Multiphase Radial Networks

The current power grid is on the cusp of modernization due to the emergence of distributed generation and controllable loads, as well as renewable energy. On one hand, distributed and renewable generation is volatile and difficult to dispatch. On the other hand, controllable loads provide significant potential for compensating for the uncertainties. In a future grid where there are thousands or millions of controllable loads and a large portion of the generation comes from volatile sources like wind and solar, distributed control that shifts or reduces the power consumption of electric loads in a reliable and economic way would be highly valuable.

Load control needs to be conducted with network awareness. Otherwise, voltage violations and overloading of circuit devices are likely. To model these effects, network power flows and voltages have to be considered explicitly. However, the physical laws that determine power flows and voltages are nonlinear. Furthermore, while distributed generation and controllable loads are mostly located in distribution networks that are multiphase and radial, most of the power flow studies focus on single-phase networks.

This thesis focuses on distributed load control in multiphase radial distribution networks. In particular, we first study distributed load control without considering network constraints, and then consider network-aware distributed load control.

Distributed implementation of load control is the main challenge if network constraints can be ignored. In this case, we first ignore the uncertainties in renewable generation and load arrivals, and propose a distributed load control algorithm, Algorithm 1, that optimally schedules the deferrable loads to shape the net electricity demand. Deferrable loads refer to loads whose total energy consumption is fixed, but energy usage can be shifted over time in response to network conditions. Algorithm 1 is a distributed gradient decent algorithm, and empirically converges to optimal deferrable load schedules within 15 iterations.

We then extend Algorithm 1 to a real-time setup where deferrable loads arrive over time, and only imprecise predictions about future renewable generation and load are available at the time of decision making. The real-time algorithm Algorithm 2 is based on model-predictive control: Algorithm 2 uses updated predictions on renewable generation as the true values, and computes a pseudo load to simulate future deferrable load. The pseudo load consumes 0 power at the current time step, and its total energy consumption equals the expectation of future deferrable load total energy request.

Network constraints, e.g., transformer loading constraints and voltage regulation constraints, bring significant challenge to the load control problem since power flows and voltages are governed by nonlinear physical laws. Remarkably, distribution networks are usually multiphase and radial. Two approaches are explored to overcome this challenge: one based on convex relaxation and the other that seeks a locally optimal load schedule.

To explore the convex relaxation approach, a novel but equivalent power flow model, the branch flow model, is developed, and a semidefinite programming relaxation, called BFM-SDP, is obtained using the branch flow model. BFM-SDP is mathematically equivalent to a standard convex relaxation proposed in the literature, but numerically is much more stable. Empirical studies show that BFM-SDP is numerically exact for the IEEE 13-, 34-, 37-, 123-bus networks and a real-world 2065-bus network, while the standard convex relaxation is numerically exact for only two of these networks.

Theoretical guarantees on the exactness of convex relaxations are provided for two types of networks: single-phase radial alternative-current (AC) networks, and single-phase mesh direct-current (DC) networks. In particular, for single-phase radial AC networks, we prove that a second-order cone program (SOCP) relaxation is exact if voltage upper bounds are not binding; we also modify the optimal load control problem so that its SOCP relaxation is always exact. For single-phase mesh DC networks, we prove that an SOCP relaxation is exact if 1) voltage upper bounds are not binding, or 2) voltage upper bounds are uniform and power injection lower bounds are strictly negative; we also modify the optimal load control problem so that its SOCP relaxation is always exact.

To seek a locally optimal load schedule, a distributed gradient-decent algorithm, Algorithm 9, is proposed. The suboptimality gap of the algorithm is rigorously characterized and close to 0 for practical networks. Furthermore, unlike the convex relaxation approach, Algorithm 9 ensures a feasible solution. The gradients used in Algorithm 9 are estimated based on a linear approximation of the power flow, which is derived with the following assumptions: 1) line losses are negligible; and 2) voltages are reasonably balanced. Both assumptions are satisfied in practical distribution networks. Empirical results show that Algorithm 9 obtains 70+ times speed up over the convex relaxation approach, at the cost of a suboptimality within numerical precision.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

Directed and elliptic flow of charged particles in Cu plus Cu collisions at root s(NN)=22.4 GeV

This paper reports results for directed flow v(1) and elliptic flow v(2) of charged particles in Cu + Cu collisions at root s(NN) = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4-GeV Cu + Cu collisions the prior observation that v1 is independent of the system size at 62.4 and 200 GeV and also extend the scaling of v(1) with eta/y(beam) to this system. The measured v(2)(p(T)) in Cu + Cu collisions is similar for root s(NN) throughout the range 22.4 to 200 GeV. We also report a comparison with results from transport model (ultrarelativistic quantum molecular dynamics and multiphase transport model) calculations. The model results do not agree quantitatively with the measured v(1)(eta), v(2)(p(T)), and v(2)(eta).

A review of flow modeling for dense medium cyclones

A critical assessment is presented for the existing fluid flow models used for dense medium cyclones (DMCs) and hydrocyclones. As the present discussion indicates, the understanding of dense medium cyclone flow is still far from the complete. However, its similarity to the hydrocyclone provides a basis for improved understanding of fluid flow in DMCs. The complexity of fluid flow in DMCs is basically due to the existence of medium as well as the dominance of turbulent particle size and density effects on separation. Both the theoretical and experimental analysis is done with respect to two-phase motions and solid phase flow in hydrocyclones or DMCs. A detailed discussion is presented on the empirical, semiempirical, and the numerical models based upon both the vorticity-stream function approach and Navier-Stokes equations in their primitive variables and in cylindrical coordinates available in literature. The existing equations describing turbulence and multiphase flows in cyclone are also critically reviewed.

Multiphase modelling of the characteristics of close coupled gas atomization

Not withstanding the high demand of metal powder for automotive and High Tech applications, there are still many unclear aspects of the production process. Only recentlyhas supercomputer performance made possible numerical investigation of such phenomena. This thesis focuses on the modelling aspects of primary and secondary atomization. Initially two-dimensional analysis is carried out to investigate the influence of flow parameters (reservoir pressure and gas temperature principally) and nozzle geometry on final powder yielding. Among the different types, close coupled atomizers have the best performance in terms of cost and narrow size distribution. An isentropic contoured nozzle is introduced to minimize the gas flow losses through shock cells: the results demonstrate that it outperformed the standard converging-diverging slit nozzle. Furthermore the utilization of hot gas gave a promising outcome: the powder size distribution is narrowed and the gas consumption reduced. In the second part of the thesis, the interaction of liquid metal and high speed gas near the feeding tube exit was studied. Both axisymmetric andnon-axisymmetric geometries were simulated using a 3D approach. The filming mechanism was detected only for very small metal flow rates (typically obtained in laboratory scale atomizers). When the melt flow increased, the liquid core overtook the adverse gas flow and entered in the high speed wake directly: in this case the disruption isdriven by sinusoidal surface waves. The process is characterized by fluctuating values of liquid volumes entering the domain that are monitored only as a time average rate: it is far from industrial robustness and capability concept. The non-axisymmetric geometry promoted the splitting of the initial stream into four cores, smaller in diameter and easier to atomize. Finally a new atomization design based on the lesson learned from previous cases simulation is presented.

Development of multiphase and multiscale mathematical models for thermal spray process

High velocity oxyfuel (HVOF) thermal spraying is one of the most significant developments in the thermal spray industry since the development of the original plasma spray technique. The first investigation deals with the combustion and discrete particle models within the general purpose commercial CFD code FLUENT to solve the combustion of kerosene and couple the motion of fuel droplets with the gas flow dynamics in a Lagrangian fashion. The effects of liquid fuel droplets on the thermodynamics of the combusting gas flow are examined thoroughly showing that combustion process of kerosene is independent on the initial fuel droplet sizes. The second analysis copes with the full water cooling numerical model, which can assist on thermal performance optimisation or to determine the best method for heat removal without the cost of building physical prototypes. The numerical results indicate that the water flow rate and direction has noticeable influence on the cooling efficiency but no noticeable effect on the gas flow dynamics within the thermal spraying gun. The third investigation deals with the development and implementation of discrete phase particle models. The results indicate that most powder particles are not melted upon hitting the substrate to be coated. The oxidation model confirms that HVOF guns can produce metallic coating with low oxidation within the typical standing-off distance about 30cm. Physical properties such as porosity, microstructure, surface roughness and adhesion strength of coatings produced by droplet deposition in a thermal spray process are determined to a large extent by the dynamics of deformation and solidification of the particles impinging on the substrate. Therefore, is one of the objectives of this study to present a complete numerical model of droplet impact and solidification. The modelling results show that solidification of droplets is significantly affected by the thermal contact resistance/substrate surface roughness.

Volume averaging of multiphase flows with hydrate formation in subsea pipelines

In oil and gas pipeline operations, the gas, oil, and water phases simultaneously move through pipe systems. The mixture cools as it flows through subsea pipelines, and forms a hydrate formation region, where the hydrate crystals start to grow and may eventually block the pipeline. The potential of pipe blockage due to hydrate formation is one of the most significant flow-assurance problems in deep-water subsea operations. Due to the catastrophic safety and economic implications of hydrate blockage, it is important to accurately predict the simultaneous flow of gas, water, and hydrate particles in flowlines. Currently, there are few or no studies that account for the simultaneous effects of hydrate growth and heat transfer on flow characteristics within pipelines. This thesis presents new and more accurate predictive models of multiphase flows in undersea pipelines to describe the simultaneous flow of gas, water, and hydrate particles through a pipeline. A growth rate model for the hydrate phase is presented and then used in the development of a new three-phase model. The conservation equations of mass, momentum, and energy are formulated to describe the physical phenomena of momentum and heat transfer between the fluid and the wall. The governing equations are solved based on an analytical-numerical approach using a Newton-Raphson method for the nonlinear equations. An algorithm was developed in Matlab software to solve the equations from the inlet to the outlet of the pipeline. The developed models are validated against a single-phase model with mixture properties, and the results of comparative studies show close agreement. The new model predicts the volume fraction and velocity of each phase, as well as the mixture pressure and temperature profiles along the length of the pipeline. The results from the hydrate growth model reveal the growth rate and location where the initial hydrates start to form. Finally, to assess the impact of certain parameters on the flow characteristics, parametric studies have been conducted. The results show the effect of a variation in the pipe diameter, mass flow rate, inlet pressure, and inlet temperature on the flow characteristics and hydrate growth rates.

Flow Regime Identification Methodology with Mcnp-X Code And Artificial Neural Network

This paper presents flow regimes identification methodology in multiphase system in annular, stratified and homogeneous oil-water-gas regimes. The principle is based on recognition of the pulse height distributions (PHD) from gamma-ray with supervised artificial neural network (ANN) systems. The detection geometry simulation comprises of two NaI(Tl) detectors and a dual-energy gamma-ray source. The measurement of scattered radiation enables the dual modality densitometry (DMD) measurement principle to be explored. Its basic principle is to combine the measurement of scattered and transmitted radiation in order to acquire information about the different flow regimes. The PHDs obtained by the detectors were used as input to ANN. The data sets required for training and testing the ANN were generated by the MCNP-X code from static and ideal theoretical models of multiphase systems. The ANN correctly identified the three different flow regimes for all data set evaluated. The results presented show that PHDs examined by ANN may be applied in the successfully flow regime identification.

A two-fluid model for tissue growth within a dynamic flow environment

We study the growth of a tissue construct in a perfusion bioreactor, focussing on its response to the mechanical environment. The bioreactor system is modelled as a two-dimensional channel containing a tissue construct through which a flow of culture medium is driven. We employ a multiphase formulation of the type presented by G. Lemon, J. King, H. Byrne, O. Jensen and K. Shakesheff in their study (Multiphase modelling of tissue growth using the theory of mixtures. J. Math. Biol. 52(2), 2006, 571–594) restricted to two interacting fluid phases, representing a cell population (and attendant extracellular matrix) and a culture medium, and employ the simplifying limit of large interphase viscous drag after S. Franks in her study (Mathematical Modelling of Tumour Growth and Stability. Ph.D. Thesis, University of Nottingham, UK, 2002) and S. Franks and J. King in their study Interactions between a uniformly proliferating tumour and its surrounding: Uniform material properties. Math. Med. Biol. 20, 2003, 47–89). The novel aspects of this study are: (i) the investigation of the effect of an imposed flow on the growth of the tissue construct, and (ii) the inclusion of a chanotransduction mechanism regulating the response of the cells to the local mechanical environment. Specifically, we consider the response of the cells to their local density and the culture medium pressure. As such, this study forms the first step towards a general multiphase formulation that incorporates the effect of mechanotransduction on the growth and morphology of a tissue construct. The model is analysed using analytic and numerical techniques, the results of which illustrate the potential use of the model to predict the dominant regulatory stimuli in a cell population.

The multiphase starburst-driven galactic wind in NGC 5394

We present a detailed study of the neutral and ionized gas phases in the galactic wind for the nearby starburst galaxy NGC 5394 based on new integral field spectroscopy obtained with the INTEGRAL fibre system at the William Herschel Telescope. The neutral gas phase in the wind is detected via the interstellar Na I D doublet absorption. After a careful removal of the stellar contribution to these lines, a significant amount of neutral gas (∼10^7 M_⊙) is detected in a central region of ∼1.75 kpc size. This neutral gas is blueshifted by ∼165 km s^−1 with respect to the underlying galaxy. The mass outflow of neutral gas is comparable to the star formation rate of the host galaxy. Simultaneously, several emission lines (Hα, [N II], [S II]) are also analysed looking for the ionized warm phase counterpart of the wind. A careful kinematic decomposition of the line profiles reveals the presence of a secondary, broader, kinematic component. This component is found roughly in the same region where the Na I D absorption is detected. It presents higher [N II]/Hα and [S II]/Hα line ratios than the narrow component at the same locations, indicative of contamination by shock ionization. This secondary component also presents blueshifted velocities, although smaller than those measured for the neutral gas, averaging to ∼−30 km s^−1. The mass and mass outflow rate of the wind is dominated by the neutral gas, of which a small fraction might be able to escape the gravitational potential of the host galaxy. The observations in this system can be readily understood within a bipolar gas flow scenario.