In November 2000, the 4th Workshop on Measurement Techniques for Stationary and Transient Multiphase Flows took place in Rossendorf. Three previous workshops of this series were national meetings; this time participants from different countries took part. The programme comprised 14 oral presentations, 9 of which are included in these proceedings in full length. A special highlight of the meeting was the main lecture "Ultrasonic doppler method for bubbly flow measurement" of Professor Masanori Aritomi, Dr. Hiroshige Kikura and Dr. Yumiko Suzuki, which was read by Dr. Hiroshige Kikura. The workshop again dealt with high-resolution phase distribution and phase velocity measurement techniques based on electrical conductivity, ultrasound, laser light and high-speed cinematography. A number of presentations were dedicated to the application of wire-mesh sensors developed by FZR for different applications used by the Technical Universities of Delft and Munich and the Tokyo Institute of Technology. The presentations were in particular: M. Aritomi, H. Kikura, Y. Suzuki (Tokyo Institute of Technology): Ultrasonic doppler method for bubbly flow measurement V. V. Kontelev, V. I. Melnikov (TU Nishny Novgorod): An ultrasonic mesh sensor for two-phase flow visualisation A. V. Duncev (TU Nishny Novgorod): Waveguide ultrasonic liquid level transducers for power generating equipment H.-M. Prasser, E. Krepper, D. Lucas, J. Zschau (FZR), D. Peters, G. Pietzsch, W. Taubert, M. Trepte (Teletronic Ingenieurbüro GmbH), Fast wire-mesh sensors for gas-liquid flows and decomposition of gas fraction profiles according to bubble size classes D. Scholz, C. Zippe (FZR): Validation of bubble size measurements with wire-mesh sensors by high-speed video observation A. Manera, H. Hartmann, W.J.M. de Kruijf, T.H.J.J. van der Hagen, R.F. Mudde, (TU Delft, IRI): Low-pressure dynamics of a natural-circulation two-phase flow loop H. Schmidt, O. Herbst, W. Kastner, W. Köhler (Siemens AG KWU): Measuring methods for the investigation of the flow phenomena during external pressure vessel cooling of the boiling water reactor SWR1000 A. Traichel, W. Kästner, S. Schefter, V. Schneider, S. Fleischer, T. Gocht, R. Hampel (HTWS Zittau/Görlitz - IPM): Verification of simulation results of mixture level transients and evaporation processes in level measurement systems using needle-shaped probes S. Richter, M. Aritomi (Tokyo Institute of Technology): Methods for studies on bubbly flow characteristics applying a new electrode-mesh tomograph

In November 2000, the 4th Workshop on Measurement Techniques for Stationary and Transient Multiphase Flows took place in Rossendorf. Three previous workshops of this series were national meetings; this time participants from different countries took part. The programme comprised 14 oral presentations, 9 of which are included in these proceedings in full length. A special highlight of the meeting was the main lecture "Ultrasonic doppler method for bubbly flow measurement" of Professor Masanori Aritomi, Dr. Hiroshige Kikura and Dr. Yumiko Suzuki, which was read by Dr. Hiroshige Kikura. The workshop again dealt with high-resolution phase distribution and phase velocity measurement techniques based on electrical conductivity, ultrasound, laser light and high-speed cinematography. A number of presentations were dedicated to the application of wire-mesh sensors developed by FZR for different applications used by the Technical Universities of Delft and Munich and the Tokyo Institute of Technology. The presentations were in particular: M. Aritomi, H. Kikura, Y. Suzuki (Tokyo Institute of Technology): Ultrasonic doppler method for bubbly flow measurement V. V. Kontelev, V. I. Melnikov (TU Nishny Novgorod): An ultrasonic mesh sensor for two-phase flow visualisation A. V. Duncev (TU Nishny Novgorod): Waveguide ultrasonic liquid level transducers for power generating equipment H.-M. Prasser, E. Krepper, D. Lucas, J. Zschau (FZR), D. Peters, G. Pietzsch, W. Taubert, M. Trepte (Teletronic Ingenieurbüro GmbH), Fast wire-mesh sensors for gas-liquid flows and decomposition of gas fraction profiles according to bubble size classes D. Scholz, C. Zippe (FZR): Validation of bubble size measurements with wire-mesh sensors by high-speed video observation A. Manera, H. Hartmann, W.J.M. de Kruijf, T.H.J.J. van der Hagen, R.F. Mudde, (TU Delft, IRI): Low-pressure dynamics of a natural-circulation two-phase flow loop H. Schmidt, O. Herbst, W. Kastner, W. Köhler (Siemens AG KWU): Measuring methods for the investigation of the flow phenomena during external pressure vessel cooling of the boiling water reactor SWR1000 A. Traichel, W. Kästner, S. Schefter, V. Schneider, S. Fleischer, T. Gocht, R. Hampel (HTWS Zittau/Görlitz - IPM): Verification of simulation results of mixture level transients and evaporation processes in level measurement systems using needle-shaped probes S. Richter, M. Aritomi (Tokyo Institute of Technology): Methods for studies on bubbly flow characteristics applying a new electrode-mesh tomograph

Fluid flow network simulation codes are commonly used as a design and analysis tool for many engineering problems such as gas distribution networks, power plants and heat pumps. Two formulations of conservation of momentum have been widely applied in fluid flow network simulation models namely those based on static pressure and those based on total pressure. The total pressure formulations are convenient in that they eliminate the difficulties associated with the calculation of the convective terms and components such as pipe junctions are treated in a straightforward manner based on total pressure losses. However, the different formulations of total pressure for compressible and incompressible flow require different formulations of the momentum conservation equation, which is inconvenient for implementation in a generic network simulation code. In this thesis a united total pressure formulation is first derived which is valid for all fluids and therefore eliminates the inconvenience of switching between the compressible and incompressible formulations. A non-iterative method for the solution of the non-isothermal discretised equations based on the total pressure formulation is then introduced and consistency is illustrated. The method appears to be very stable for subsonic flows, while rapid steady state convergence is observed. A systematic comparison is also done with traditional static pressure based methods and the similarities and differences between the two formulations are illuminated. The different time scales involved in the simulation of transient flow in fluid networks are problematic when conventional fixed time step methods are used for time-wise integration. The time scales associated with acoustic and kinematic wave phenomena as well as storage effects can differ by orders in magnitude. This thesis also presents a simple adaptive time step algorithm which can be readily used in conjunction with all the commonly used first order methods for fluid flow networks. Two test problems are selected to demonstrate the efficiency and savings obtained with this procedure. The adaptive time step algorithm correctly selects appropriate time steps for all phenomena and significant computational savings are observed for accurate integration. In addition, a procedure is implemented which automatically selects the appropriate integration method. The resulting algorithm is a fully adaptive algorithm which switches between a fully implicit method and a semi-implicit method. Two test problems are once again used to demonstrate the efficiency and savings. The fully adaptive algorithm correctly selects appropriate methods for all phenomena and significant additional computational savings are observed. / Thesis (PhD (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2013

Fluid flow network simulation codes are commonly used as a design and analysis tool for many engineering problems such as gas distribution networks, power plants and heat pumps. Two formulations of conservation of momentum have been widely applied in fluid flow network simulation models namely those based on static pressure and those based on total pressure. The total pressure formulations are convenient in that they eliminate the difficulties associated with the calculation of the convective terms and components such as pipe junctions are treated in a straightforward manner based on total pressure losses. However, the different formulations of total pressure for compressible and incompressible flow require different formulations of the momentum conservation equation, which is inconvenient for implementation in a generic network simulation code. In this thesis a united total pressure formulation is first derived which is valid for all fluids and therefore eliminates the inconvenience of switching between the compressible and incompressible formulations. A non-iterative method for the solution of the non-isothermal discretised equations based on the total pressure formulation is then introduced and consistency is illustrated. The method appears to be very stable for subsonic flows, while rapid steady state convergence is observed. A systematic comparison is also done with traditional static pressure based methods and the similarities and differences between the two formulations are illuminated. The different time scales involved in the simulation of transient flow in fluid networks are problematic when conventional fixed time step methods are used for time-wise integration. The time scales associated with acoustic and kinematic wave phenomena as well as storage effects can differ by orders in magnitude. This thesis also presents a simple adaptive time step algorithm which can be readily used in conjunction with all the commonly used first order methods for fluid flow networks. Two test problems are selected to demonstrate the efficiency and savings obtained with this procedure. The adaptive time step algorithm correctly selects appropriate time steps for all phenomena and significant computational savings are observed for accurate integration. In addition, a procedure is implemented which automatically selects the appropriate integration method. The resulting algorithm is a fully adaptive algorithm which switches between a fully implicit method and a semi-implicit method. Two test problems are once again used to demonstrate the efficiency and savings. The fully adaptive algorithm correctly selects appropriate methods for all phenomena and significant additional computational savings are observed. / Thesis (PhD (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2013

With the high costs of flight testing, especially at hypersonic speeds, ground based facility testing of scramjets becomes an attractive option. The expansion tube is the only facility currently that can offer full flight property duplication at the total pressures and total enthalpies required, while maintaining correct chemical composition. Due to difficulties with short test times and unsteady flow phenomena, scramjet testing in these facilities has not been thoroughly investigated. This study examines these issues, in order to explore the practicality of testing a full ’tip to tail’ scramjet engine at a true flight replication condition in an expansion tube facility. An investigation was initially undertaken on the large X3 expansion tunnel facility to maximise test time and core flow, aimed at producing a 30 km altitude, Mach 10 flow condition. This was identified as the limitation point of the T4 reflected shock tunnel, which has generally been accepted to produce reliable scramjet data for propulsion tests. Using a condition that is also able to be produced in the T4 facility, will permit direct comparison of data between the two facilities in the future, providing confidence in results from expansion tube facilities. Both experimental measurements and numerical calculations showed that the limitation of the test time was due to large boundary layer growth after transition, which engulfed the entire core flow 200 μs into the test time. This phenomenon is likely to affect all scramjet duplication conditions in expansion tubes, as the flow properties are conducive to boundary layer transition occurring. Two solutions where proposed and investigated in order to overcome the flow disruption caused by boundary layer transition; the use of a steady expansion nozzle at the acceleration tube exit; the use of hydrogen as an accelerator gas. Since the smaller X2 facility had a Mach 10 steady expansion nozzle and X3 was decommissioned for the free piston driver to be upgraded, the investigation was shifted to X2. Due to a restricted test time of 550 μs, the static pressure of the flow condition was increased to allow a reduction in the length of the scramjet (pressure-length scaling). A combination of experimental and numerical calculations of the facility was used to define the flow properties. With the confidence of overcoming the phenomenon associated with boundary transition in the X2 facility, numerical modelling of the X3 facility with a steady expansion nozzle was then undertaken to show a 1 ms condition could be produced. Although initially promising, the hydrogen accelerator gas solution requires further investigation. A two dimensional scramjet was designed with upstream injection for testing in X2. This was a three shock inlet with a constant area combustor and a planar thrust surface. Since the flow condition involved changes in flow properties during the test time, aninvestigation of the appropriateness of a quasi steady analysis was undertaken. Using a fuel off simulation of the scramjet duct with the transient inflow properties from the X2 facility nozzle exit, the convective terms for pressure were shown to be two orders of magnitude larger than local terms indicating the dominance of the convective terms change in flow properties at any location allowing quasi-steady flow to be assumed. A normalisation procedure was developed to deal with the transient nature of the data and to accurately represent the axial progression of the gas through the duct. The numerical simulations were also used to show that both flow establishment was achieved and that impulsive starting of the intake would occur. Experimentation with the scramjet using static pressure measurements throughout the body side of the engine provided verification of supersonic combustion. This was verified by the doubling of the static pressure from the start to the end of the combustor for an air test gas, whereas with a nitrogen test gas no significant change in pressure occurred. Effects of fuel equivalence ratio, injector size and cowl position were also investigated. A net inviscid thrust was predicted, using the quasi-steady flow analysis, indicating a specific impulse of 183 s. This work provides evidence to validate the use of expansion tube facilities for experimental testing of scramjets at flight duplication conditions. Limitations due to boundary layer transition flow effects has been shown to be avoidable. Numerical simulations of the facilities showed good agreement with experimental measurements, allowing definition of freestream properties and can now be applied to further scramjet conditions with confidence. Stable, supersonic combustion was shown to be produced for these expansion tube conditions. Coupling the transient simulation of the flow condition with a numerical calculation of the fuel off experimental scramjet has been useful in both verification of the design and performance predictions. Appropriate techniques have been presented to analyse scramjet pressure and thrust measurements where transient effects are present in the freestream.

O ensaio de palheta de campo em material siltoso levanta problemáticas relacionadas à sua execução e à interpretação dos resultados. O presente trabalho apresenta uma modelagem em poroelasticidade do ensaio de palheta. O modelo de ensaio é definido pelo problema de rotação de um cilindro infinito em um solo poroso. A solução do problema é buscada adotando um comportamento poroelástico não linear fictício tal que a resposta é localmente equivalente àquela do comportamento plástico perfeito. O modelo considera que a rotação do cilindro gera deformações volumétricas não desprezíveis e que a solução fechada de poropressão é garantida por um módulo de cisalhamento equivalente. As soluções do campo de tensões e deformações descritas por equações analíticas são obtidas numericamente por meio do método de diferenças finitas. O modelo é avaliado através de uma comparação com os resultados de simulação em elementos finitos e de solução do modelo com o uso de software de álgebra computacional. Os critérios de Tresca e de Drucker-Prager são considerados nas avaliações. As condições de drenagem foram estudadas através da curva característica de drenagem no espação grau de drenagem U versus velocidade normalizada V. A influência da rigidez e da resistência do material sobre o fenômeno de dissipação foram interpretados no mesmo espaço U ×V . Também foi demonstrado que a dissipação é sensível à definição da zona de influência. O parâmetro numérico de discretização do domínio para o método de diferenças finitas para a obtenção de resultados de boa precisão foi identificado. O modelo foi aplicado para a modelagem do ensaio de palheta em resíduo de zinco e para a interpretação dos resultados desse ensaio. Demonstrou-se que o modelo proposto permite a identificações dos padrões de ensaio que garantem os comportamentos desejados, além de possibilitar o estudo da sensibilidade do processo de dissipação em relação à rigidez e resistência do material. / The Field Vane Test in silty materials raises problematics related to its execution and interpretation of results. The work presents a model for the vane test based on poroelasticity. The modeling of the test is characterized by the problem of the rotation of an infinite cylinder in a porous soil. The solution of the problem is sought by adopting a fictitious non-linear poroelastic behavior such that the answer is locally equivalent to the one corresponding to a perfect plastic behavior. The revised model assumes that rotation of the cylinder does generate volumetric deformation which one is not negligible and the close form of solution for pore pressure is guaranteed by an equivalent shear modulus. The solutions of stresses and displacement field are obtained numerically using the finite difference method. The model is evaluated for materials characterized by two criteria Tresca and Drucker-Prager through a comparison of results obtained by simulation in finite element model and by simulation using computer algebra software. The drainage conditions have been studied through the drainage characteristic curve U ×V . The influence of the stiffness and strength of the material on the dissipation phenomena were interpreted in the same space U ×V . It was also demonstrated that the dissipation process is sensitive to the definition of the influence zone. The numerical parameters to obtain good precision results were identified. The model was applied to the modeling Vane Test in zinc residue and the interpretation of experimental results. It has been shown that the proposed model allows the identification of test patterns that ensures the desired drainage behavior and allows the study of the sensitivity of the dissipation process for stiffness and strength of the material.

Cryogenic carbon capture process can be used to efficiently eliminate CO2 emissions from fossil-fueled power plants. The energy-storing embodiment of cryogenic carbon capture (ES-CCC) integrates energy storage with cryogenic carbon capture and uses natural gas as a refrigerant. ES-CCC captures CO2 from slowly varying or steady-state sources even as it absorbs and replaces large amounts of energy on the grid for energy storage. These large transients occur in the LNG generation as the process moves through energy storing to energy recovery operations. Additionally, raw natural gas often includes CO2 that forms an azeotrope with ethane. Breaking this azeotrope and separating CO2 from other hydrocarbons to meet natural gas pipeline and liquefied natural gas (LNG) standards is very energy intensive. The purpose of this work is to (a) describe a dynamic heat exchanger that reduces the heat exchanger performance and efficiency losses experienced under transient conditions and (b) introduce an alternative extractive distillation system for CO2 separation from ethane that requires less capital and has a lower operating cost than the conventional system for the same purification. This investigation demonstrates theoretically and experimentally that the dynamic heat exchangers can absorb sudden and large changes in flow rates and other properties without compromising either the heat exchanger efficiency or creating thermal or other stresses. These heat exchangers play an essential role in the ES-CCC process. Designs for retrofitting existing heat exchangers and for replacing existing heat exchangers with new designs are both theoretically and experimentally tested. The ES-CCC process requires natural gas processing to meet pipeline and LNG standards in many applications, depending primarily on the CO2 content of locally available NG. The energy, cost, and dynamic response of such processing hinges primarily on the most difficult step, breaking the CO2-ethane azeotrope. This project proposes and analyzes an alternative process for breaking this azeotrope and a control scheme that dramatically improves the dynamic response of natural gas processing plants, including steady and transient control scheme and processing simulations. These contributions to the ES-CCC capture process all have much broader applications in many chemical and energy processes. These contributions to ES-CCC and other industrial processes improve energy efficiency and dynamic performance of many processes and are ready for larger scale demonstration.

Field-scale transport through unsaturated soil is influenced by surface and subsurface boundary conditions, and the spatial variability of state soil variables. The objective of this thesis is to examine the relative importance of the spatial redistribution of surface water versus spatial variability of soil properties on long-term transient water flow and transport under semi-arid conditions. The field-scale transport (34 yr) of a surface applied tracer (chloride), spatial variability of other pedogenic tracers, and surface water redistribution over a 19 mo fallow period were measured in a catchment basin. In 1966 and 1971, a chloride tracer (KCl) was surface applied to plots (6.1 m x 90 m, Chernozemic soil) near Saskatoon, Saskatchewan. In 2000 and 2001, 262 soil cores were taken along and perpendicular to one KCl strip. Soil layering at each core was recorded and samples were analysed for chloride concentration, electrical conductivity, bulk density and water content. Sulphate and nitrate concentrations were measured on selected cores. The site is level by common definitions, with a very slight concave depression (1.8% grade) midway along the KCl strip and a slight grade (¡Ü2.1%) perpendicular to the KCl strip. Measured water recharge indicated slight differences in surface slope had a marked effect on redistribution of water and spatial distribution of the chloride tracer. An estimated 90% of redistributed water was subsequently used by plants and 10% resulted in an increase in deep drainage. A varved layer had a strong influence on the subsurface redistribution of water and chloride below the root zone. There were sharp horizontal transitions between areas of slow and faster transport, which corresponded to sharp increases in catchment area and water recharge. Small surface depressions, which controlled pedogenic transport and soil formation, have been filled in by tillage translocation. Spatial variability of soil horizon thickness (and associated hydraulic properties) had little effect on transport of chloride after 34 yr. Computer simulations also suggest substantial surface redistribution of precipitation and snowmelt. In contrast to the measured chloride data, the model was sensitive to changes in hydraulic properties and horizon thickness in the root zone. Surface water redistribution was the primary factor controlling long-term transport.

Field-scale transport through unsaturated soil is influenced by surface and subsurface boundary conditions, and the spatial variability of state soil variables. The objective of this thesis is to examine the relative importance of the spatial redistribution of surface water versus spatial variability of soil properties on long-term transient water flow and transport under semi-arid conditions. The field-scale transport (34 yr) of a surface applied tracer (chloride), spatial variability of other pedogenic tracers, and surface water redistribution over a 19 mo fallow period were measured in a catchment basin. In 1966 and 1971, a chloride tracer (KCl) was surface applied to plots (6.1 m x 90 m, Chernozemic soil) near Saskatoon, Saskatchewan. In 2000 and 2001, 262 soil cores were taken along and perpendicular to one KCl strip. Soil layering at each core was recorded and samples were analysed for chloride concentration, electrical conductivity, bulk density and water content. Sulphate and nitrate concentrations were measured on selected cores. The site is level by common definitions, with a very slight concave depression (1.8% grade) midway along the KCl strip and a slight grade (¡Ü2.1%) perpendicular to the KCl strip. Measured water recharge indicated slight differences in surface slope had a marked effect on redistribution of water and spatial distribution of the chloride tracer. An estimated 90% of redistributed water was subsequently used by plants and 10% resulted in an increase in deep drainage. A varved layer had a strong influence on the subsurface redistribution of water and chloride below the root zone. There were sharp horizontal transitions between areas of slow and faster transport, which corresponded to sharp increases in catchment area and water recharge. Small surface depressions, which controlled pedogenic transport and soil formation, have been filled in by tillage translocation. Spatial variability of soil horizon thickness (and associated hydraulic properties) had little effect on transport of chloride after 34 yr. Computer simulations also suggest substantial surface redistribution of precipitation and snowmelt. In contrast to the measured chloride data, the model was sensitive to changes in hydraulic properties and horizon thickness in the root zone. Surface water redistribution was the primary factor controlling long-term transport.

Horizontal wells have been used to increase reservoir recovery, especially in unconventional reservoirs, and hydraulic fracturing has been applied to further extend the contact with the reservoir to increase the efficiency of development. In the past, many models, analytical or numerical, were developed to describe the flow behavior in horizontal wells with fractures. Source solution is one of the analytical/semi-analytical approaches. To solve fractured well problems, source methods were advanced from point sources to volumetric source, and pressure change inside fractures was considered in the volumetric source method. This study aims at developing a method that can predict horizontal well performance and the model can also be applied to horizontal wells with multiple fractures in complex natural fracture networks. The method solves the problem by superposing a series of slab sources under transient or pseudosteady-state flow conditions. The principle of the method comprises the calculation of semi-analytical response of a rectilinear reservoir with closed outer boundaries. A statistically assigned fracture network is used in the study to represent natural fractures based on the spacing between fractures and fracture geometry. The multiple dominating hydraulic fractures are then added to the natural fracture system to build the physical model of the problem. Each of the hydraulic fractures is connected to the horizontal wellbore, and the natural fractures are connected to the hydraulic fractures through the network description. Each fracture, natural or hydraulically induced, is treated as a series of slab sources. The analytical solution of superposed slab sources provides the base of the approach, and the overall flow from each fracture and the effect between the fractures are modeled by applying superposition principle to all of the fractures. It is assumed that hydraulic fractures are the main fractures that connect with the wellbore and the natural fractures are branching fractures which only connect with the main fractures. The fluid inside of the branch fractures flows into the main fractures, and the fluid of the main fracture from both the reservoir and the branch fractures flows to the wellbore. Predicting well performance in a complex fracture network system is extremely challenged. The statistical nature of natural fracture networks changes the flow characteristic from that of a single linear fracture. Simply using the single fracture model for individual fracture, and then adding the flow from each fracture for the network could introduce significant error. This study provides a semi-analytical approach to estimate well performance in a complex fracture network system.