During the production and transportation of oil and gas, flow assurance issues may occur due to the solid deposits that are formed and carried by the flowing fluid. Solid deposition may cause serious damage and possible failure to production equipment in the flow lines. The major flow assurance problems that are faced in the fields are concerned with asphaltene, wax and scale deposition, as well as hydrate formations. Hydrates, wax and asphaltene deposition are mostly addressed in deep-water environments, where fluid flows through a long path with a wide range of pressure and temperature variations (Hydrates are generated at high pressure and low temperature conditions). In fact, a large change in the thermodynamic condition of the fluid yields phase instability and triggers solid deposit formations. In contrast, scales are formed in aqueous phase when some incompatible ions are mixed. Among the different flow assurance issues in hydrocarbon reservoirs, asphaltenes are the most complicated one. In fact, the difference in the nature of these molecules with respect to other hydrocarbon components makes this distinction. Asphaltene molecules are the heaviest and the most polar compounds in the crude oils, being insoluble in light n-alkenes and readily soluble in aromatic solvents. Asphaltene is attached to similarly structured molecules, resins, to become stable in the crude oils. Changing the crude oil composition and increasing the light component fractions destabilize asphaltene molecules. For instance, in some field situations, CO₂ flooding for the purpose of enhanced oil recovery destabilizes asphaltene. Other potential parameters that promote asphaltene precipitation in the crude oil streams are significant pressure and temperature variation. In fact, in such situations the entrainment of solid particulates in the flowing fluid and deposition on different zones of the flow line yields serious operational challenges and an overall decrease in production efficiency. The loss of productivity leads to a large number of costly remediation work during a well life cycle. In some cases up to $5 Million per year is the estimated cost of removing the blockage plus the production losses during downtimes. Furthermore, some of the oil and gas fields may be left abandoned prematurely, because of the significance of the damage which may cause loss about $100 Million. In this dissertation, we developed a robust wellbore model which is coupled to our in-house developed compositional reservoir model (UTCOMP). The coupled wellbore/reservoir simulator can address flow restrictions in the wellbore as well as the near-wellbore area. This simulator can be a tool not only to diagnose the potential flow assurance problems in the developments of new fields, but also as a tool to study and design an optimum solution for the reservoir development with different types of flow assurance problems. In addition, the predictive capability of this simulator can prescribe a production schedule for the wells that can never survive from flow assurance problems. In our wellbore simulator, different numerical methods such as, semi-implicit, nearly implicit, and fully implicit schemes along with blackoil and Equation-of-State compositional models are considered. The Equation-of-State is used as state relations for updating the properties and the equilibrium calculation among all the phases (oil, gas, wax, asphaltene). To handle the aqueous phase reaction for possible scales formation in the wellbore a geochemical software package (PHREEQC) is coupled to our simulator as well. The governing equations for the wellbore/reservoir model comprise mass conservation of each phase and each component, momentum conservation of liquid, and gas phase, energy conservation of mixture of fluids and fugacity equations between three phases and wax or asphaltene. The governing equations are solved using finite difference discretization methods. Our simulation results show that scale deposition is mostly initiated from the bottom of the wellbore and near-wellbore where it can extend to the upper part of the well, asphaltene deposition can start in the middle of the well and the wax deposition begins in the colder part of the well near the wellhead. In addition, our simulation studies show that asphaltene deposition is significantly affected by CO₂ and the location of deposition is changed to the lower part of the well in the presence of CO₂. Finally, we applied the developed model for the mechanical remediation and prevention procedures and our simulation results reveal that there is a possibility to reduce the asphaltene deposition in the wellbore by adjusting the well operation condition. / text

<div>Many two-phase flow systems might undergo flow instabilities even if the system is adiabatic but operates near the saturation conditions, especially in vertical flow conditions. Such instabilities are caused by <i>flashing</i> of the fluid in flow. Flashing is a sudden phase change in the fluid caused when local saturation enthalpy falls below the fluid enthalpy and the excess energy is used as latent heat for gas generation.</div><div> In the current analysis, a mathematical model is presented for analysis of such instability analytically. The conservation equations have been obtained by statistical averaging in time and space. Then, the concerned system is divided into various regions based on flow conditions, and these averaged equations are used to describe the flow. For flashing-based instability, two parameters are derived from constitutive relationships for the fluid. These two parameters are <i>Flashing Boundary</i> and <i>Gas Generation due to Flashing</i>. These parameters provide for the closure of the mathematical model. Some simple models for flashing have been developed and discussed.</div><div> The mathematical model is then solved analytically for <i>Uniform Heat</i> and <i>Flat Model</i> for the heater and flashing region respectively. The solution is in terms of the characteristic equation which is used to predict the onset of instability caused by flashing. The results are then plotted on the Subcooling-Phase Change number plane. It is observed that inlet and outlet restrictions in the flow does <b>not</b> affect the onset of flashing induced instability as the flow rate is coupled with the pressure drop of the system. This is important as these restrictions play a major role in other two-phase flow instabilities such as <i>Density Wave Oscillations</i></div><div> Finally, the stability boundary in the stability plane is compared to experimental data present for flashing. The comparison was made with data of S. Shi, A. Dixit, and F. Inada. The stability boundary satisfactorily agrees with the experimental data thus corroborating the present mathematical model and analysis.</div>

[pt] Na produção de gás e petróleo em campos de águas profundas são comumente encontrados trechos verticais de dutos (risers) na aproximação final à plataforma. Nesta configuração, podem ocorrer escoamentos bifásicos no regime de golfadas severas (severe slug) que gera alternância na produção da fase gasosa e líquida. Esta alternância é caracterizada por períodos de produção de gás sem líquido seguido de altas taxas de produção de ambas as fases. O regime severo de golfadas é geralmente descrito em quatro fases: formação da golfada, produção da golfada, rompimento da golfada pela fase gás e fluxo reverso do que restou da fase líquida. Este regime induz o escoamento a condições mais extremas do que um outro regime, visto que resultam em um aumento de pressão no duto durante a formação da golfada e em um aumento na taxa de produção durante a expulsão da mesma. O presente trabalho consiste da simulação numérica do regime de golfadas severas para um trecho de tubulação horizontal seguido de outro vertical, assim como apresentar uma análise de um regime estatisticamente permanente. A previsão do escoamento é obtida utilizando-se uma formulação unidimensional baseada no modelo de Drift. A freqüência das golfadas é comparada com outros estudos numéricos da literatura, obtendo-se uma concordância bastante satisfatória. / [en] In the gas and oil offshore deep water production is usual to find risers in Production Unit final approach. Regarding to this configuration, two-phase flows can evolve to a severe slug regime that create gas and liquid alternate production. This cyclic behavior is characterized by periods of gas production followed by very high liquid and gas flow rates. The severe slug flow regime is normally described as occurring in four phases: slug formation, slug production, blowout, and liquid fallback. This flow regime introduces new conditions that can be found in other regimes because of the pressure increase during the slug formation and the large flow rates during the slug production. The present work consists of the numerical simulation of the severe slug flow regime into a horizontal pipeline section followed by a vertical section, as well as to present an statistically steady state analysis. The flow prediction is obtained through a one-dimensional formulation based on the Drift Flux Model. The slug frequency is compared with other numerical studies available in the literature, and a very satisfactory agreement is obtained.

O escoamento bifásico é muito importante para vários ramos industriais. As misturas bifásicas podem escoar em diversas configurações, as quais são chamadas de padrões de escoamento, e que, ao longo dos anos, receberam diversas classificações. Neste trabalho foi estudado o padrão bolhas dispersas em duto anular. Os trabalhos sobre escoamento bifásico em geometria anular são mais escassos, especialmente quando se trata de dutos anulares de grande dimensão. Dentre os modelos existentes para a modelagem do escoamento bifásico, trabalhamos com o modelo de mistura (Drift Flux), pois, apesar das altas velocidades superficiais, o escoamento em bolhas não possui um comportamento homogêneo, já que a fase gasosa ainda escoa na região central e com uma velocidade maior do que a velocidade da mistura. Foi utilizado o modelo de mistura unidimensional para a modelagem do escoamento bifásico água-ar, pois, além de bastante preciso, é de baixo custo computacional e fácil implementação. Neste trabalho empregamos, a princípio, equações elaboradas para dutos circulares, as quais foram adaptadas para a geometria anular fazendo-se uso do conceito de diâmetro hidráulico. As equações foram implementadas no software Mathematica® e as previsões para o holdup e queda de pressão foram comparadas com dados experimentais próprios. O trabalho experimental foi realizado no laboratório de escoamentos multifásicos do Núcleo de Engenharia Térmica e de Fluidos (NETeF) da EESC-USP, o qual conta com uma instalação experimental em estado operacional para a simulação de escoamentos bifásicos vertical e inclinados em duto anular. Foram colhidas medidas de holdup e queda de pressão bifásica. A concordância entre os dados experimentais e as previsões foi muito satisfatória. / Two-phase flow is present in a wide range of industrial processes. Such flows occur in many geometrical configurations known as flow patterns. In this work is presented a study about bubbly flow in annular duct. Works on two-phase flow in large annular geometry are very scarce. Among the existing models for modeling two-phase flow we chose the drift flux model because the bubbly flow does not behave as a homogeneous mixture, despite of the high superficial velocities, since the gas phase flows faster in the core region in comparison with the mixture velocity. The onedimensional drift flux model besides to be accurate is of low computational cost and easy implementation. The goal of this work is to analyze the prediction of holdup and pressure drop in annular ducts of large diameter. The set of equations used was adapted from those of circular ducts by using the hydraulic diameter concept. The equations were implemented in Mathematica software and holdup and pressure drop prediction were compared with experimental data taken in the experimental facility of the Thermal-fluids Engineering Laboratory of EESC-USP. It was collected measure of holdup and pressure drop. The agreement between experimental data and model predictions is encouraging.

O escoamento bifásico é muito importante para vários ramos industriais. As misturas bifásicas podem escoar em diversas configurações, as quais são chamadas de padrões de escoamento, e que, ao longo dos anos, receberam diversas classificações. Neste trabalho foi estudado o padrão bolhas dispersas em duto anular. Os trabalhos sobre escoamento bifásico em geometria anular são mais escassos, especialmente quando se trata de dutos anulares de grande dimensão. Dentre os modelos existentes para a modelagem do escoamento bifásico, trabalhamos com o modelo de mistura (Drift Flux), pois, apesar das altas velocidades superficiais, o escoamento em bolhas não possui um comportamento homogêneo, já que a fase gasosa ainda escoa na região central e com uma velocidade maior do que a velocidade da mistura. Foi utilizado o modelo de mistura unidimensional para a modelagem do escoamento bifásico água-ar, pois, além de bastante preciso, é de baixo custo computacional e fácil implementação. Neste trabalho empregamos, a princípio, equações elaboradas para dutos circulares, as quais foram adaptadas para a geometria anular fazendo-se uso do conceito de diâmetro hidráulico. As equações foram implementadas no software Mathematica® e as previsões para o holdup e queda de pressão foram comparadas com dados experimentais próprios. O trabalho experimental foi realizado no laboratório de escoamentos multifásicos do Núcleo de Engenharia Térmica e de Fluidos (NETeF) da EESC-USP, o qual conta com uma instalação experimental em estado operacional para a simulação de escoamentos bifásicos vertical e inclinados em duto anular. Foram colhidas medidas de holdup e queda de pressão bifásica. A concordância entre os dados experimentais e as previsões foi muito satisfatória. / Two-phase flow is present in a wide range of industrial processes. Such flows occur in many geometrical configurations known as flow patterns. In this work is presented a study about bubbly flow in annular duct. Works on two-phase flow in large annular geometry are very scarce. Among the existing models for modeling two-phase flow we chose the drift flux model because the bubbly flow does not behave as a homogeneous mixture, despite of the high superficial velocities, since the gas phase flows faster in the core region in comparison with the mixture velocity. The onedimensional drift flux model besides to be accurate is of low computational cost and easy implementation. The goal of this work is to analyze the prediction of holdup and pressure drop in annular ducts of large diameter. The set of equations used was adapted from those of circular ducts by using the hydraulic diameter concept. The equations were implemented in Mathematica software and holdup and pressure drop prediction were compared with experimental data taken in the experimental facility of the Thermal-fluids Engineering Laboratory of EESC-USP. It was collected measure of holdup and pressure drop. The agreement between experimental data and model predictions is encouraging.

<p>Flerfasestrømning er av betydning i en lang rekke anvendelser, blant annet i olje- og gassindustrien, i den kjemiske og i prosessindustrien, inkludert i varmepumpende systemer, samt i sikkerhetsanalyse av kjernekraftverk. Denne avhandlingen analyserer modeller for tofasestrømning, og metoder for numerisk løsning av disse modellene. Den er derfor ett bidrag til å utvikle pålitelige ingeniørverktøy for flerfase-anvendelser. Slike verktøy trengs og forventes av ingeniører i industrien.</p><p>Den tilnærmede Riemann-løseren framsatt av Roe har blitt studert. Roe-skjema for tre ulike modeller for tofasestrømning har blitt implementert i rammen av en standard numerisk algoritme for løsning av hyperbolske bevaringslover. Disse skjemaene har blitt analysert ved hjelp av referanse-regnestykker fra litteraturen, og ved sammenlikning med hverandre.</p><p>Et Roe-skjema for den fire-liknings ettrykks tofluid-modellen har blitt implementert, og det har blitt vist at en andreordens utvidelse basert på bølge-dekomponering og fluksdifferanse-splitting virker godt, og gir forbedrede resultater sammenliknet med det førsteordens skjemaet.</p><p>Det har blitt foreslått et Roe-skjema for en fem-liknings totrykks tofluid-modell med trykkrelaksering. Bruken av analoge numeriske metoder for fire-liknings- og fem-liknings-modellene gjorde det mulig med en direkte sammenlikning av en metode med og uten trykkrelaksering. Numeriske eksperiment demonstrerte at de to framgangsmåtene konvergerte til samme resultat, men at den fem-liknings trykkrelakserings-metoden var betydelig mer dissipativ, særlig for kontakt-diskontinuiteter. Videre viste beregninger at selv om fem-liknings-modellen har reelle egenverdier, så produserte den oscillasjoner for tilfeller der fire-liknings-modellen hadde komplekse egenverdier.</p><p>Et Roe-skjema har blitt konstruert for driftfluks-modellen med generelle lukningslover. Roe-matrisen er helt analytisk for det tilfellet at man kan anvende Zuber-Findlay-slipp-loven som beskriver boblestrømning. Dermed er dette Roe-skjemaet mer effektivt enn tidligere fullt numeriske Roe-skjema for driftfluks-modellen.</p><p>En isentropisk diskret-nivå-flerfasemodell har blitt presentert. En diskusjon av hvordan man kan ta hensyn til ulike interfase-trykkmodeller har blitt gitt. Med de passende modellene for interfase-trykk og -fart, var samsvaret svært godt mellom diskret-nivå-modellen og det fem-liknings Roe-skjemaet.</p><p>Multi-steg- (MUSTA) metoden har som siktemål å komme nær oppstrøms-metodene i nøyaktighet, samtidig som den bevarer enkelheten til sentrerte skjema. Her har metoden blitt brukt på driftfluks-modellen. Når antallet steg økes, nærmer resultatene fra MUSTA-metoden seg det man får med Roe-metoden. De gode resultatene til MUSTA-metoden er avhengige av at man bruker et stort nok lokalt grid. Derfor er hovedfordelen med MUSTA-metoden at den er enkel, snarere enn at man sparer regnetid.</p><p>En karakteristikk-basert metode for å spesifisere grensebetinglser for flerfase-modeller har blitt testet, og funnet å virke godt for transiente problem.</p> / <p>This thesis analyses models for two-phase flows and methods for the numerical resolution of these models. It is therefore one contribution to the development of reliable design tools for multiphase applications. Such tools are needed and expected by engineers in a range of fields, including in the oil and gas industry.</p><p>The approximate Riemann solver of Roe has been studied. Roe schemes for three different two-phase flow models have been implemented in the framework of a standard numerical algorithm for the solution of hyperbolic conservation laws. The schemes have been analysed by calculation of benchmark tests from the literature, and by comparison with each other.</p><p>A Roe scheme for the four-equation one-pressure two-fluid model has been implemented, and a second-order extension based on wave decomposition and flux-difference splitting was shown to work well and to give improved results compared to the first-order scheme. The convergence properties of the scheme were tested on smooth and discontinuous solutions.</p><p>A Roe scheme has been proposed for a five-equation two-pressure two-fluid model with pressure relaxation. The use of analogous numerical methods for the five-equation and four-equation models allowed for a direct comparison of a method with and without pressure relaxation. Numerical experiments demonstrated that the two approaches converged to the same results, but that the five-equation pressure-relaxation method was significantly more dissipative, particularly for contact discontinuities. Furthermore, even though the five-equation model with instantaneous pressure relaxation has real eigenvalues, the calculations showed that it produced oscillations for cases where the four-equation model had complex eigenvalues.</p><p>A Roe scheme has been constructed for the drift-flux model with general closure laws. For the case of the Zuber-Findlay slip law describing bubbly flows, the Roe matrix is completely analytical. Hence the present Roe scheme is more efficient than previous fully numerical Roe schemes for the drift-flux model.</p><p>An isentropic discrete-equation multiphase model has been presented. The incorporation of different interfacial-pressure models was discussed, and examples were given. With the adequate models for the interfacial pressure and velocity, the agreement was very good between the discrete-equation model and the five-equation Roe scheme.</p><p>The flux-limiter centred (FLIC) scheme was tested for the four-equation two-fluid model. Only the first-order version (FORCE) of the scheme was found to work well, but it was rather diffusive. The purpose of the multi-stage (MUSTA) method is to come close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. Here it has been applied to the drift-flux model. As the number of stages was increased, the results of the MUSTA scheme approached those of the Roe method. The good results of the MUSTA scheme were dependent on the use of a large-enough local grid. Hence, the main advantage of the MUSTA scheme is its simplicity.</p><p>A multiphase characteristic-based boundary-condition method has been tested, and it was shown to be workable for transient problems.</p>

No description available.

Flerfasestrømning er av betydning i en lang rekke anvendelser, blant annet i olje- og gassindustrien, i den kjemiske og i prosessindustrien, inkludert i varmepumpende systemer, samt i sikkerhetsanalyse av kjernekraftverk. Denne avhandlingen analyserer modeller for tofasestrømning, og metoder for numerisk løsning av disse modellene. Den er derfor ett bidrag til å utvikle pålitelige ingeniørverktøy for flerfase-anvendelser. Slike verktøy trengs og forventes av ingeniører i industrien. Den tilnærmede Riemann-løseren framsatt av Roe har blitt studert. Roe-skjema for tre ulike modeller for tofasestrømning har blitt implementert i rammen av en standard numerisk algoritme for løsning av hyperbolske bevaringslover. Disse skjemaene har blitt analysert ved hjelp av referanse-regnestykker fra litteraturen, og ved sammenlikning med hverandre. Et Roe-skjema for den fire-liknings ettrykks tofluid-modellen har blitt implementert, og det har blitt vist at en andreordens utvidelse basert på bølge-dekomponering og fluksdifferanse-splitting virker godt, og gir forbedrede resultater sammenliknet med det førsteordens skjemaet. Det har blitt foreslått et Roe-skjema for en fem-liknings totrykks tofluid-modell med trykkrelaksering. Bruken av analoge numeriske metoder for fire-liknings- og fem-liknings-modellene gjorde det mulig med en direkte sammenlikning av en metode med og uten trykkrelaksering. Numeriske eksperiment demonstrerte at de to framgangsmåtene konvergerte til samme resultat, men at den fem-liknings trykkrelakserings-metoden var betydelig mer dissipativ, særlig for kontakt-diskontinuiteter. Videre viste beregninger at selv om fem-liknings-modellen har reelle egenverdier, så produserte den oscillasjoner for tilfeller der fire-liknings-modellen hadde komplekse egenverdier. Et Roe-skjema har blitt konstruert for driftfluks-modellen med generelle lukningslover. Roe-matrisen er helt analytisk for det tilfellet at man kan anvende Zuber-Findlay-slipp-loven som beskriver boblestrømning. Dermed er dette Roe-skjemaet mer effektivt enn tidligere fullt numeriske Roe-skjema for driftfluks-modellen. En isentropisk diskret-nivå-flerfasemodell har blitt presentert. En diskusjon av hvordan man kan ta hensyn til ulike interfase-trykkmodeller har blitt gitt. Med de passende modellene for interfase-trykk og -fart, var samsvaret svært godt mellom diskret-nivå-modellen og det fem-liknings Roe-skjemaet. Multi-steg- (MUSTA) metoden har som siktemål å komme nær oppstrøms-metodene i nøyaktighet, samtidig som den bevarer enkelheten til sentrerte skjema. Her har metoden blitt brukt på driftfluks-modellen. Når antallet steg økes, nærmer resultatene fra MUSTA-metoden seg det man får med Roe-metoden. De gode resultatene til MUSTA-metoden er avhengige av at man bruker et stort nok lokalt grid. Derfor er hovedfordelen med MUSTA-metoden at den er enkel, snarere enn at man sparer regnetid. En karakteristikk-basert metode for å spesifisere grensebetinglser for flerfase-modeller har blitt testet, og funnet å virke godt for transiente problem. / This thesis analyses models for two-phase flows and methods for the numerical resolution of these models. It is therefore one contribution to the development of reliable design tools for multiphase applications. Such tools are needed and expected by engineers in a range of fields, including in the oil and gas industry. The approximate Riemann solver of Roe has been studied. Roe schemes for three different two-phase flow models have been implemented in the framework of a standard numerical algorithm for the solution of hyperbolic conservation laws. The schemes have been analysed by calculation of benchmark tests from the literature, and by comparison with each other. A Roe scheme for the four-equation one-pressure two-fluid model has been implemented, and a second-order extension based on wave decomposition and flux-difference splitting was shown to work well and to give improved results compared to the first-order scheme. The convergence properties of the scheme were tested on smooth and discontinuous solutions. A Roe scheme has been proposed for a five-equation two-pressure two-fluid model with pressure relaxation. The use of analogous numerical methods for the five-equation and four-equation models allowed for a direct comparison of a method with and without pressure relaxation. Numerical experiments demonstrated that the two approaches converged to the same results, but that the five-equation pressure-relaxation method was significantly more dissipative, particularly for contact discontinuities. Furthermore, even though the five-equation model with instantaneous pressure relaxation has real eigenvalues, the calculations showed that it produced oscillations for cases where the four-equation model had complex eigenvalues. A Roe scheme has been constructed for the drift-flux model with general closure laws. For the case of the Zuber-Findlay slip law describing bubbly flows, the Roe matrix is completely analytical. Hence the present Roe scheme is more efficient than previous fully numerical Roe schemes for the drift-flux model. An isentropic discrete-equation multiphase model has been presented. The incorporation of different interfacial-pressure models was discussed, and examples were given. With the adequate models for the interfacial pressure and velocity, the agreement was very good between the discrete-equation model and the five-equation Roe scheme. The flux-limiter centred (FLIC) scheme was tested for the four-equation two-fluid model. Only the first-order version (FORCE) of the scheme was found to work well, but it was rather diffusive. The purpose of the multi-stage (MUSTA) method is to come close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. Here it has been applied to the drift-flux model. As the number of stages was increased, the results of the MUSTA scheme approached those of the Roe method. The good results of the MUSTA scheme were dependent on the use of a large-enough local grid. Hence, the main advantage of the MUSTA scheme is its simplicity. A multiphase characteristic-based boundary-condition method has been tested, and it was shown to be workable for transient problems.

No description available.

<p>  </p> <p>The study of two-phase flow in different orientations can allow for greater understanding of the fundamentals of two-phase flow dynamics. While a large amount of work has been performed for vertical flows and recent work has been done for horizontal flows, limited research has been done studying inclined upward two-phase flows between those two orientations. Studying two-phase flows at various inclinations is important for developing physical models and simulations of two-phase flow systems and understanding the changes between what is observed for symmetric vertical flows and asymmetric horizontal flows. The present work seeks to systematically characterize the effects of inclination on adiabatic concurrent air-water two-phase flows in straight pipes. An experimental database is established for local and global two-phase flow parameters in a novel inclinable 25.4 mm inner diameter test facility using four-sensor conductivity probes, high speed video capabilities, a ring-type impedance meter, a pressure transducer, and a gamma densitometer. Rotatable measurement ports are employed to allow for local conductivity probe measurements across the flow profile to capture asymmetric parameter distributions during experiments without stopping the flow. Some of the major effects of inclination are investigated, including the effects on flow regime transition, bubble distribution, frictional pressure loss, and relative motion between the two phases. Flow visualization and machine-learning methods are employed to identify the transitions between flow regimes for inclined orientations, and these transitions are compared against existing theoretical flow regime transition criteria proposed in literature. The theoretical transitions in literature agree well with both methods for vertical flow, but additional work is necessary for angles between 0 degrees and 60 degrees. The effect of inclination on two-phase frictional pressure drop is explored, and a novel adaption of the Lockhart-Martinelli pressure drop correlation is proposed, which is able to predict the pressure drop for the conditions investigated with an absolute percent difference of 2.6%. To explore the relationships between orientation, void fraction, and relative motion, one-dimensional drift flux analyses are performed for the data at each angle investigated. It is observed that the relative velocity between phases decreases as the angle is reduced, with a relative velocity near zero at some intermediate angles and a negative relative velocity for near-horizontal orientations.  Existing modeling capabilities that have been developed for vertical and horizontal flows are evaluated based on the local two-phase parameters collected at multiple orientations. The performance of the one-dimensional interfacial area transport equation for vertical and horizontal flows is tested against experimental data and a novel model for horizontal and inclined-upward bubbly flows is proposed. Finally, an evaluation of existing momentum transfer relations is performed for the two-fluid model using three-dimensional computational fluid dynamics tools for horizontal and inclined. The prediction of the void fraction distribution and gas velocity profiles are compared against experimental data, and improvements to the lift force model are identified based on changes in the relative velocity between phases. </p>