Numerical study of swirling one- and two-phase turbulent flows in a cylindrical channel

Translated from Zhurnal Prikladnoi Mekhanikl I Tekhnicheskoi Fiziki, No. 2, pp. 51–58, March–April, 1988.     

Mathematical modeling of jet interaction with a high-enthalpy flow in an expanding channel

Results of modeling the interaction of a plane supersonic jet with a supersonic turbulent high-enthalpy flow in a channel are reported. The problem is solved in a two-dimensional formulation at external flow Mach numbers M_∞ = 2.6 and 2.8 and at high values of the total temperature of the flow T ₀ = 1800–2000 K. The mathematical model includes full averaged Navier-Stokes equations supplemented with a two-equation turbulence model and an equation that describes the transportation of the injected substance. The computations are performed by using the ANSYS Fluent 12.1 software package. Verification of the computational technique is performed against available experimental results on transverse injection of nitrogen and helium jets. The computed and experimental results are demonstrated to agree well. For the examined problems, in addition to surface distributions of characteristics, fields of flow parameters are obtained, which allow one to reproduce specific features that can be hardly captured in experiments. Parametric studies show that an increase in the angle of inclination and the mass flow rate of the jet leads to an increase in the depth of jet penetration into the flow, but more intense separated flows and shock waves are observed in this case.     

Slow two-phase flow through a sinusoidal channel

The two-phase flow through a symmetric sinusoidal channel is studied by means of a regular perturbation analysis, where the small parameter is defined as the ratio between the amplitude of variation of the channel wall and the average thickness of the non-wetting phase. Results are valid for Reynolds numbers of the same order of magnitude as that of the expansion parameter. It is thus found that the fluid-fluid interface presents a wavy shape characterized by an amplitude and a phase-shift with respect to the fixed solid-fluid interface. Instabilities of the two-phase flow can arise for large values of the viscosity, flow rate and phase thickness ratios. Results are expected to be a first step towards the understanding of the hydrodynamics of trickle bed reactors, where several flow regimes are possible.     

Boundary conditions on a rigid surface in a two-phase flow

A study is made of the boundary conditions on a rigid surface in a two-component disperse flow. Appropriate boundary conditions are obtained for the kinetic equation and macroscopic equations of a pseudogas of solid particles proposed in [1–3]. The reasons for the occurrence of bubbles in two-phase systems are discussed. On the basis of the similitude parameters of the kinetic equation of the pseudogas, disperse systems are classified generally on the basis of the concentration of solid particles and their diameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 46–51, May–June, 1980.I thank V. P. Myasnikov for suggesting the problem and for a helpful discussion.     

Investigation on two-phase flow-induced vibrations of a piping structure with an elbow

The dynamic behaviors of a horizontal piping structure with an elbow due to the two-phase flow excitation are experimentally investigated. The effects of flow patterns and superficial velocities on the pressure pulsations and vibration responses are evaluated in detail. A strong partition coupling algorithm is used to calculate the flow-induced vibration (FIV) responses of the pipe, and the theoretical values agree well with the experimental results. It is found that the lateral and axial vibration responses of the bend pipe are related to the momentum flux of the two-phase flow, and the vibration amplitudes of the pipe increase with an increase in the liquid mass flux. The vertical vibration responses are strongly affected by the flow pattern, and the maximum response occurs in the transition region from the slug flow to the bubbly flow. Moreover, the standard deviation (STD) amplitudes of the pipe vibration in three directions increase with an increase in the gas flux for both the slug and bubbly flows. The blockage of liquid slugs at the elbow section is found to strengthen the vibration amplitude of the bend pipe, and the water-blocking phenomenon disappears as the superficial gas velocity increases.     

Stratification of boiling two-phase flow in a single horizontal channel

The objective of this study is to investigate experimentally the stratification phenomena of boiling two-phase flow in a uniformly heated horizontal channel. Two-phase flow stratification due to gravity effects, and consequently its thermal and hydrodynamic behavior, under steady state conditions, have been determined by measuring 16 top and 16 bottom wall temperatures. Six distinct wall temperature profiles are found, and the corresponding flow patterns are discussed. A dimensionless number has been formulated for the prediction of the occurrence of different flow patterns.     

Investigation of turbulence modification in a non-reactive two-phase flow

In a two-phase flow the influence of a dispersed phase on the turbulence properties of a continuous phase, known as turbulence modification, is investigated. An experimental approach is discussed that is suitable for studying the decay of grid-generated turbulence in a vertically orientated wind tunnel with a cross-section large enough to avoid influences from walls. Phase Doppler anemometry is used to characterize both single and two-phase flow by measuring mean axial and radial velocity components, velocity fluctuations, turbulent kinetic energy, and integral time scales. By direct comparison of results from single- and two-phase flows, the feedback of the dispersed phase on the continuous phase can be isolated. The data is used to deduce a source term for particle-induced turbulence production appropriate for a numerical simulation of the flow, based on the Reynolds-averaged Navier-Stokes equations. Although of special importance for a detailed understanding of turbulent two-phase combustion, additional complexity introduced by evaporation and chemical reactions is avoided by using glass beads as dispersed phase.     

Dynamic instabilities of boiling two-phase flow in a single horizontal channel

Dynamic instabilities of two-phase flow associated with refrigerant R-11 in a uniformly heated horizontal in-tube boiling system were experimentally investigated. An experimental setup was designed and built to work in a wide range of mass fluxes G [75–1050 kg/(m² s)], heat fluxes q (0–100 kW/m²), and fluid inlet temperatures T_inlet (2–24°C). Dynamic instability data were obtained under various working conditions. The dependence of oscillation amplitude and period on system parameters is discussed, and the boundaries of various oscillations are located on the steady-state characteristic curves.     

Effect of inlet subcooling on two-phase flow oscillations in a vertical boiling channel

In this work, results of experimental research to investigate the effects of heat transfer augmentation and inlet subcooling on two-phase flow instabilities are presented. For this purpose, a simple set-up was designed and built. The effect of inlet subcooling was investigated using different heat transfer surfaces and inlet temperatures at constant heat input of 415 W. Freon-11 has been used as the test fluid, and the experiments were carried out for six heater tubes having different heat transfer surfaces. Inlet temperatures were in the range of –9.8°C to 38°C. The results indicate that, in the range of present experiments, the system becomes more stable, that it's instability boundary moves into lower mass flow rates, with increase in the inlet subcooling. However, the amplitudes and the periods of the oscillations increase with increase in the inlet subcooling. For some of the tested surfaces there was a particular inlet subcooling, above and below which the system's stability decreased.

---

Der Einfluß der Eintrittsunterkühlung auf die Oszillationen einer Zwei-Phasen-Strömung in einem senkrechten Siede-Kanal

Zusammenfassung In dieser Arbeit werden die Ergebnisse einer experimentellen Untersuchung des Einflusses von Verbesserungen der Wärmeübertragung und der Eintrittsunterkühlung auf Instabilitäten der Zweiphasenströmung dargestellt. Zu diesem Zweck wurde ein einfacher Aufbau konstruiert und aufgebaut. Der Einfluß der Eintrittsunterkühlung wurde unter Benutzung verschiedener Wärmeübertragungsoberflächen und Eintrittstemperaturen bei einer konstanten Wärmezufuhr von 415 W untersucht. Als Testfluid wurde Freon-11 benutzt. Die Experimente wurden für sechs Heizrohre mit verschiedenen Wärmeübertragungsoberflächen durchgeführt. Die Eintrittstemperaturen lagen in diesem Bereich von –9.8°C bis 38°C. Die Ergebnisse zeigen, daß in dem Bereich der vorliegenden Untersuchungen das System mit einem Ansteigen der Eintrittsunterkühlung stabiler wird und daß dessen Instabilitätsgrenze sich zu niedrigeren Massenstromdichten bewegt. Die Amplituden und Perioden der Schwingungen steigen mit Zunahme der Eintrittsunterkühlung an. Für einige der getesteten Oberflächen existierte eine bestimmte Eintrittsunterkühlung, über bzw. unter welcher die Stabilität des Systems abnahm.

     

Analytic solution to the micropolar-fluid flow through a semi-porous channel with an expanding or contracting wall

The flow of a micropolar fluid in a semi-porous channel with an expanding or contracting wall is investigated. The governing equations are reduced to ordinary ones by using similar transformations. To get the analytic solution to the problem, the homotopy analysis method （HAM） is employed to obtain the expressions for velocity fields. Graphs are sketched and discussed for various parameters, especially the effect of the expansion ratio on velocity and micro-rotation fields.     

The effects of inlet conditions and expansion ratio on the onset of flow reversal in swirling flow in a sudden expansion

Experiments are reported in which the minimum swirl intensity required to produce a central recirculation zone in a swirling sudden expansion flow is determined as a function of expansion ratio and inlet conditions. Using a swirl generator which allows for independent variation of velocity profile shape and swirl number, it is shown that an inlet tangential velocity distribution with a large solid body vortex core or an axial velocity profile with a maximum on the axis will lead to a higher critical swirl.     

Investigation of nanofluid flow in a vertical channel considering polynomial boundary conditions by Akbari-Ganji's method

In this research, a vertical channel containing a laminar and fully developed nanofluid flow is investigated. The channel surface's boundary conditions for temperature and volume fraction functions are considered qth-order polynomials. The equations related to this problem have been extracted and then solved by the AGM and validated through the Runge-Kutta numerical method and another similar study. In the study, the effect of parameters, including Grashof number, Brownian motion parameter, etc., on the motion, velocity, temperature, and volume fraction of nanofluids have been analyzed. The results demonstrate that increasing the Gr number by 100% will increase the velocity profile function by 78% and decrease the temperature and fraction profiles by 20.87% and 120.75%. Moreover, rising the Brownian motion parameter in five different sizes (0.1, 0.2, 0.3, 0.4, and 0.5) causes lesser velocity, about 24.3% at first and 4.35% at the last level, and a maximum 52.86% increase for temperature and a 24.32% rise for Ψ occurs when N_b rises from 0.1 to 0.2. For all N_t values, at least 55.44%, 18.69%, for F(η), and Ω(η), and 20.23% rise for Ψ(η) function is observed. Furthermore, enlarging the N_r parameter from 0.25 to 0.1 leads F(η) to rise by 199.7%, fluid dimensionless temperature, and dimensional volume fraction to decrease by 18% and 92.3%. In the end, a greater value of q means a more powerful energy source, amplifying all velocity, temperature, and volume fraction functions. The main novelty of this research is the combined convection qth-order polynomials boundary condition applied to the channel walls. Moreover, The AMG semi-analytical method is used as a novel method to solve the governing equations.     

Calculation of the swirling flow of an ideal gas in a laval nozzle

The numerical solution of the problem of the motion of a swirling flow of an ideal gas in a Laval nozzle in axisymmetric formulation is obtained by the method of stabilization. As a result, a number of effects appear that are essentially not one-dimensional, in particular, the drawing-in of the sonic line into the nozzle, an effect that leads to a decrease in the nozzle's expansion coefficient. The dependence of this coefficient on the intensity of the swirling is obtained. A number of problems connected with the control of the expansion of a gas through a Laval nozzle and with variation of the thrust of a nozzle can be solved successfully in cases where a rotary motion is imparted to the flow of gas exhausted from the nozzle. Investigation of such a swirling flow in [1, 2] and a number of other papers are based on a one-dimensional model of gas flow, which makes it possible in principle to obtain integrated characteristics of the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 72–76, September–October, 1971.     

Numerical investigation of a perturbed swirling annular two-phase jet

A swirling annular gas–liquid two-phase jet flow system has been investigated by solving the compressible, time-dependent, non-dimensional Navier–Stokes equations using highly accurate numerical methods. The mathematical formulation for the flow system is based on an Eulerian approach with mixed-fluid treatment while an adjusted volume of fluid method is utilised to account for the gas compressibility. Surface tension effects are captured by a continuum surface force model. Swirling motion is applied at the inlet while a small helical perturbation is also applied to initiate the instability. Three-dimensional spatial direct numerical simulation has been performed with parallelisation of the code based on domain decomposition. The results show that the flow is characterised by a geometrical recirculation zone adjacent to the nozzle exit and by a central recirculation zone further downstream. Swirl enhances the flow instability and vorticity and promotes liquid dispersion in the cross-streamwise directions. A dynamic precessing vortex core is developed demonstrating that the growth of such a vortex in annular configurations can be initiated even at low swirl numbers, in agreement with experimental findings. Analysis of the averaged results revealed the existence of a geometrical recirculation zone and a swirl induced central recirculation zone in the flow field.     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

Effect of heat transfer augmentation on two-phase flow instabilities in a vertical boiling channel

The effect of different heater surface configurations on two-phase flow instabilities has been investigated in a single channel, forced convection, open loop, up-flow system. Freon-11 is used as the test fluid, and six different heater tubes with various inside surface configurations have been tested at five different heat inputs. In addition to temperature and pressure recordings, high speed motion pictures of the two-phase flow were taken for some of the experiments to study the two-phase flow behavior at different operating points. Experimental results are shown on system pressure drop versus mass flow rate curves, and stability boundaries are also indicated on these curves. Comparison of different heater tubes is made by the use of the stability boundary maps and the plots of inlet throttling necessary to stabilize the system versus mass flow rate. Tubes with internal springs were found to be more stable than the other tubes.     

Spatial instability of a swirling flow

The instability of a swirling flow of an inviscidand incompressible fluid is studied on the assumption that the wavenumber k=k_r + ik_i of the disturbance is complex while its frequency ω is real. This implies that the disturbance grows with distance along the axis of the swirling flow, but it does not grow with time. The occurrence of such disturbance is called spatial instability, in contrast to the temporal instability, in which k is a real number and ω=ω_r + iω_i is a complex one. The results show that spatial instability analysis is a useful tool for the comprehensive understanding of the instability behaviours of a swirling flow.