垂直管道中密相两相流动的解析解

本文依据文献[1]的密相两相流动的数学模型,对垂直圆管中密相两相流动进行了解析求解,分别得到了连续相和分散相的速度解析表达式.在相间阻力与相间速度差成比例时,除了在离管壁面很近的薄区之外,管道流动规律与达西渗流定律完全一致.本文验证了文献[1]的密相两相流动数学模型的假定在本文情形下是合理的.     

Modelling of an Homogeneous Equilibrium Mixture Model (HEM)

We present here a model for two phase flows which is simpler than the 6-equations models (with two densities, two velocities, two temperatures) but more accurate than the standard mixture models with 4 equations (with two densities, one velocity and one temperature). We are interested in the case when the two-phases have been interacting long enough for the drag force to be small but still not negligible. The so-called Homogeneous Equilibrium Mixture Model (HEM) that we present is dealing with both mixture and relative quantities, allowing in particular to follow both a mixture velocity and a relative velocity. This relative velocity is not tracked by a conservation law but by a closure law (drift relation), whose expression is related to the drag force terms of the two-phase flow. After the derivation of the model, a stability analysis and numerical experiments are presented.     

Numerical simulation of scour around a submarine pipeline using computational fluid dynamics and discrete element method

Scour under a submarine pipeline can lead to structural failure; hence, a good understanding of the scour mechanism is paramount. Various numerical methods have been proposed to simulate scour, such as potential flow theory and single-phase and two-phase turbulent models. However, these numerical methods have limitations such as their reliance on calibrated empirical parameters and inability to provide detailed information. This paper investigates the use of a coupled computational fluid dynamics-discrete element method (CFD-DEM) model to simulate scour around a pipeline. The novelty of this work is to use CFD-DEM to extract detailed information, leading to new findings that enhance the current understanding of the underlying mechanisms of the scour process. The simulated scour evolution and bed profile are found to be in good agreement with published experimental results. Detailed results include the contours of the fluid velocity and fluid pressure, particle motion and velocity, fluid forces on the particles, and inter-particle forces. The sediment transport rate is calculated using the velocity of each single particle. The quantitative analysis of the bed load layer is also presented. The numerical results reveal three scour stages: onset of scour, tunnel erosion, and lee-wake erosion. Particle velocity and force distributions show that during the tunnel erosion stage, the particle motion and particle–particle interactive forces are particularly intense, suggesting that single-phase models, which are unable to account for inter-particle interactions, may be inadequate. The fluid pressure contours show a distinct pressure gradient. The pressure gradient force is calculated and found to be comparable with the drag force for the onset of scour and the tunnel erosion. However, for the lee-wake erosion, the drag force is shown to be the dominant mechanism for particle movements.     

Active Control of Separating Boundary Layer

Flow control refers to the ability to alter flows with the aim of achieving a desired effect: examples include drag reduction, noise attenuation, improved mixing or increased combustion efficiency among many other industrial applications. The reduction and control of the viscous drag force exerted on bodies moving in a fluid is of great technical interest. Several active and passive methods to achieve a delay of separation in the boundary layer have been developed and are being developed. In this paper we present a new concept for boundary layer separation control that is based on the synthetic jet concept which converts acoustic oscillations into mean fluid motions. We use synthetic jets as methods for Coanda effect's amplification. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Maximization of the lift/drag ratio of airfoils with a turbulent boundary layer: Sharp estimates,approximation, and numerical solutions

The lift/drag ratio of an airfoil placed in an incompressible attached flow is maximized taking into account the viscosity in the boundary-layer approximation. An exact solution is constructed. The situation when the resulting solutions are not in the admissible class of univalent flows is discussed. A procedure is proposed for determining physically feasible airfoils (with a univalent flow region) with a high lift/drag ratio. For this purpose, a class of airfoils is constructed that are determined by a twoparameter function approximating the found exact solution to the variational problem. For this class, the ranges of free parameters leading to physically feasible flows are found. The results are verified by computing a turbulent boundary layer using Eppler’s method, and airfoils with a high lift/drag ratio in an attached flow are detected.     

Empirical validation of a queueing approach to uninterrupted traffic flows

In this paper, the use of queueing theory for modeling uninterrupted traffic flows is evaluated. Empirical data on speeds and flows are used to evaluate speeds generated by the different queueing models. Using the Theil inequality coefficient as evaluation criterion, the speeds generated by the queueing models are compared to the empirical speeds. Queueing models that best fit the observed speeds are obtained. It appears that traffic flow on a highway during non-congested hours is best described using a M/G/1 queueing model. During the congested hours however, the state dependent queueing GI/G/z models are more realistic. Because the queueing models describe the empirical data well, they can also be used to evaluate potential improvements in existing traffic conditions. Received: April 2005 / Revised version: June 2005 AMS classification: 60K30, 68M20     

Object-oriented modelling and simulation of heat exchangers with finite element methods

The paper describes the derivation of finite-element models of one-dimensional fluid flows with heat transfer in pipes, using the Galerkin/least-squares approach. The models are first derived for one-phase flows, and then extended to homogeneous two-phase flows. The resulting equations have then been embedded in the context of object-oriented system modelling; this allows one to combine the fluid flow model with a model for other phenomena such as heat transfer, as well as with models of other discrete components such as pumps or valves, to obtain complex models of heat exchangers. The models are then validated by simulating a typical heat exchanger plant.     

Active control of cylinder wake

The objective of this paper is to develop an efficient active control algorithm for manipulating wake flows past a solid cylinder in an electrically low-conducting fluid (e.g. seawater). The intent is to avoid both vortex shedding and flow separation from the body. It is expected to reduce the mean drag significantly. This is achieved through the introduction of a Lorentz force in the azimuthal direction generated by an array of permanent magnets and electrodes located on the solid structure. With the use of a symmetric and static Lorentz force over the entire surface of the cylinder, the vortex shedding behind the cylinder weakens and eventually disappears completely when the Lorentz force is sufficiently large. The localized Lorentz force along the rear surface of the cylinder was also used to control the vortex shedding behind the cylinder. In this case, numerical results show that the efficiency of the localized Lorentz force in controlling the flow is to that of the Lorentz force distributed over the whole surface.     

Advanced continuum modelling of gas-particle flows beyond the hydrodynamic limit

The accurate prediction of dilute gas-particle flows using Euler–Euler models is challenging because particle–particle collisions are usually not dominant in such flows. In other words, in dilute flows the particle Knudsen number is not small enough to justify a Chapman–Enskog expansion about the collision-dominated near-equilibrium limit. Moreover, due to the fluid drag and inelastic collisions, the granular temperature in gas-particle flows is often small compared to the mean particle kinetic energy, implying that the particle-phase Mach number can be very large. In analogy to rarefied gas flows, it is thus not surprising that two-fluid models fail for gas-particle flows with moderate Knudsen and Mach numbers. In this work, a third-order quadrature-based moment method, valid for arbitrary Knudsen number, coupled with a fluid solver has been applied to simulate dilute gas-particle flow in a vertical channel with particle-phase volume fractions between 0.0001 and 0.01. In order to isolate the instabilities that arise due to fluid-particle coupling, a fluid mass flow rate that ensures that turbulence would not develop in a single phase flow (Re = 1380) is employed. Results are compared with the predictions of a two-fluid model with standard kinetic theory based closures for the particle phase. The effect of the particle-phase volume fraction on flow instabilities leading to particle segregation is investigated, and differences with respect to the two-fluid model predictions are examined. The influence of the discretization on the solution of both models is investigated using three different grid resolutions. Radial profiles of phase velocities and particle concentration are shown for the case with an average particle volume fraction of 0.01, showing the flow is in the core-annular regime.     

Toward design and fabrication of wind-driven vehicles: Procedure to optimize the threshold of driving forces

Wind energy has been continuously considered as a green, available, and economical alternative source of energy. For centuries, the transformed wind energy to drag-force has been used for transportation in watercrafts. With improvement of aerodynamics, the airfoil was invented to create and use a higher magnitude aerodynamic force, lift-force, in order to elevate airplanes. Later, the lift-force was horizontally applied as the thrust force in land/water wind-crafts. Whereas in airplanes horizontal airfoils (wing) create a vertical lift-force, installed vertical airfoils (wing-sail) produce a horizontal lift-force in wind-crafts. Therefore, this force can be used as thrust (driving) force in lift-based ice, water, and land vehicles. If the prevailing wind is constantly available, the vehicle speed can even exceed the wind velocity. Due to the complex kinematics of such vehicles, however, it should be noted that there would be always an optimum for the thrust force in order to control and navigate the vehicle to the destination point, and to avoid the severe undesired side-forces. This optimum is calculated in wind-craft trajectory software (WTS) which requires many inputs, including variable and constant parameters. Variable parameters consist of wind direction and magnitude in addition to vehicle’s position, velocities, and accelerations. On the other hand, design characteristics of the wind-driven vehicle are known as constant parameters. The land-yacht body’s drag is an unknown constant parameter which alters according to the relative wind. This implies that several wind tunnel experiment in different wind directions and speeds are required in order to obtain the drag coefficients.Therefore in order to bypass the wind tunnel measurements, this study aims to propose a fast and economical procedure to find the aforementioned drag coefficient by integration of a measurement and by a simulation approach. The obtained data can be later used in the optimization and control module of the WTS. The performance of this procedure has been investigated using a case study. For this purpose, a 1:4 prototype three-wheel land-yacht is first designed and fabricated. The land-yacht comprises of three major parts; horizontal airfoil (axle), vertical airfoil, and body. The dimensions of these elements are obtained after development of a code based on kinematics of the land-yacht. The axle is designed to increase the stability of the land-yacht, whereas the shape of the body is intended to produce a low drag coefficient in various directions. Furthermore, a set of experiments has been conducted to measure the body drag of the land-yacht in a direction parallel to the relative wind. This experiment is later used to develop and validate a computational fluid dynamics (CFD) model in order to estimate the drag of the land-yacht body in its various directions against the relative wind. The results show the adequate efficacy of this procedure to provide the required data for the optimization and control module of the WTS.     

Numerical investigation of body and hole effects on the cavitating flow behind a disk cavitator at extremely low cavitation numbers

Supercavitation has been recently proposed as an effective method for drag reduction of underwater vehicles. For this purpose, a cavitator is used to generate a large continuous bubble to cover the vehicle. Extremely low cavitation numbers are required in these applications to provide a sufficiently voluminous cavity. In this study, the effects of body presence inside the cavity of a disk-shaped cavitator are studied numerically. An element based finite volume approach is used for solving Reynolds-averaged Navier-Stokes equations. In order to validate the numerical scheme, turbulence and cavitation models, the cavitating flow over a hemispherical nosed cylinder body and also behind a disk-shaped cavitator has been simulated and the numerical results are validated against existing analytical and experimental data. Next, the effects of body presence inside the cavity of the disk-cavitator on the cavity characteristics are investigated. Results show that the cavity length is slightly smaller when it closes on the body, in comparison with the freely-closing cavity. Furthermore, the effects of adding a concentric hole with various diameters to the disk-cavitator on the cavity features are examined. The degree of drag reduction and changes of the cavity dimensions are determined at various cavitation numbers.     

Skew normal measurement error models

In this paper we define a class of skew normal measurement error models, extending usual symmetric normal models in order to avoid data transformation. The likelihood function of the observed data is obtained, which can be maximized by using existing statistical software. Inference on the parameters of interest can be approached by using the observed information matrix, which can also be computed by using existing statistical software, such as the Ox program. Bayesian inference is also discussed for the family of asymmetric models in terms of invariance with respect to the symmetric normal distribution showing that early results obtained for the normal distribution also holds for the asymmetric family. Results of a simulation study and an analysis of a real data set analysis are provided.     

Time of complete displacement of an immiscible compressible fluid by water in porous media: Application to gas migration in a deep nuclear waste repository

A system of evolutionary partial differential equations (PDEs) describing the two-phase flow of immiscible fluids, such as water–gas, through porous media is studied. In this formulation, the wetting and nonwetting phases are treated to be incompressible and compressible, respectively. This treatment is indeed necessary when a compressible nonwetting phase is subjected to compression during confinement. The system of PDEs consists of an evolution equation for the wetting-phase saturation and an evolution equation for the pressure in the nonwetting phase. This system is applied to the problem of unsaturated flows to assess gas migration and two-phase flow through engineered and geological barriers for a deep repository for radioactive waste. This paper is primarily concerned with the large time behavior of solutions of this system. Under some realistic assumptions on the data, we derive estimates of the speed of propagation of the gas by water in porous media. Namely, we establish estimates of time stabilization for the water saturation to a constant limit profile. The analysis is based on the energy methods whose main idea involves deriving and studying suitable ordinary differential inequalities. We show that the time of complete displacement of a gas by water may be at most infinite or finite depending essentially on the power parameters defining the capillary pressure and the relative permeabilities. This result is then illustrated with two examples in the context of gas migration in a deep nuclear waste repository. We consider Van Genuchten’s and Brooks–Corey’s models for a two-phase water–gas system.     

Numerical modeling of unsteady flow in steam turbine stage

This work deals with numerical solution of unsteady flow in turbine stage. We use models of compressible single-phase flow of air and two-phase flow of wet steam. Presented numerical methods are based on different stator-rotor matching algorithms, as well as different numerical schemes. Numerical results achieved by both methods and flow models are discussed.     

Slow viscous flow through a membrane built up from porous cylindrical particles with an impermeable core

This paper concerns the slow viscous flow through an aggregate of concentric clusters of porous cylindrical particles with Happel boundary condition. An aggregate of clusters of porous cylindrical particles is considered as a hydro-dynamically equivalent to solid cylindrical core with concentric porous cylindrical shell. The Brinkman equation inside the porous cylindrical shell and the Stokes equation outside the porous cylindrical shell in their stream function formulations are used. The drag force acting on each porous cylindrical particle in a cell is evaluated. In certain limiting cases, drag force converges to pre-existing analytical results, such as, the drag on a porous circular cylinder and the drag on a solid cylinder in a Happel unit cell. Representative results are then discussed and presented in graphical forms. The hydrodynamic permeability of the membrane built up from porous particles is evaluated. The variation of hydrodynamic permeability with different parameters is graphically presented. Some new results are reported for flow pattern in the porous region. Being in resemblance with the model of colloid particles with a coating of porous layers due to adsorption phenomenon, results obtained through this model can be useful to study the membrane filtration process.     

Numerical simulation of cavitation around a two-dimensional hydrofoil using VOF method and LES turbulence model

In this paper simulation of cavitating flow over the Clark-Y hydrofoil is reported using the large eddy simulation (LES) turbulence model and volume of fluid (VOF) technique. We applied an incompressible LES modelling approach based on an implicit method for the subgrid terms. To apply the cavitation model, the flow has been considered as a single fluid, two-phase mixture. A transport equation model for the local volume fraction of vapour is solved and a finite rate mass transfer model is used for the vapourization and condensation processes. A compressive volume of fluid (VOF) method is applied to track the interface of liquid and vapour phases. This simulation is performed using a finite volume, two phase solver available in the framework of the OpenFOAM (Open Field Operation and Manipulation) software package. Simulation is performed for the cloud and super-cavitation regimes, i.e., σ = 0.8, 0.4, 0.28. We compared the results of two different mass transfer models, namely Kunz and Sauer models. The results of our simulation are compared for cavitation dynamics, starting point of cavitation, cavity’s diameter and force coefficients with the experimental data, where available. For both of steady state and transient conditions, suitable accuracy has been observed for cavitation dynamics and force coefficients.     

Finite element simulation of compressible particle-laden gas flows

A macroscopic two-fluid model of compressible particle-laden gas flows is considered. The governing equations are discretized by a high-resolution finite element method based on algebraic flux correction. A multidimensional limiter of TVD type is employed to constrain the local characteristic variables for the continuous gas phase and conservative fluxes for a suspension of solid particles. Special emphasis is laid on the efficient computation of steady state solutions at arbitrary Mach numbers. To avoid stability restrictions and convergence problems, the characteristic boundary conditions are imposed weakly and treated in a fully implicit manner. A two-way coupling via the interphase drag force is implemented using operator splitting. The Douglas-Rachford scheme is found to provide a robust treatment of the interphase exchange terms within the framework of a fractional-step solution strategy. Two-dimensional simulation results are presented for a moving shock wave and for a steady nozzle flow.     

A network approach to the modelling of active magnetic bearings

The modelling of active magnetic bearings based on a network approach is considered. Unlike in the standard modelling approach, where a linearization of the current-force relation for the centred shaft position is used, network models permit to include the position dependence of the bearing force in the force model. This becomes necessary when model based controllers are used to stabilize a magnetically supported shaft in tracking applications.

The approach is based on the well known application of network models to magnetic circuits. Further simplifying assumptions are discussed which allow one to obtain a network with a limited number of lumped parameters describing the magnetic behaviour of a magnetic bearing. The modelling of a combined radial and axial bearing serves as an example for the application of the proposed approach. Furthermore, the fitting of the network based model to measured characteristic force curves is discussed. In this context, a method for including saturation effects in the model is sketched.     

Numerical study of wave overtopping of a seawall supported by porous structures

This paper presents a VOF-based two-phase flow model for the simulation of wave interactions with seawall supported by a porous terrace. Firstly, the model was verified against laboratory data in a simple case for wave overtopping of a vertical wall. Comparison of computed and measured wave properties showed reasonably good agreement. The model was then applied to study the interactions of waves and a seawall protected by porous structures with a permeable terrace. The application results showed that the overtopping rate was strongly related to the energy dissipation through the drag force; the porous reef and terrace were very effective to produce a low crest type seawall. It is concluded that there exist two optimum values of porosity of the submerged reef, about of 0.25 and 0.7, that give minimum overtopping rates. Whereas, there is an effective range of porosity of the permeable terrace varying from 0.4 to 0.65 for significantly reducing the overtopping rate. The verification results confirm that the VOF-based two-phase flow model is sufficient robust to simulate the wave overtopping of coastal structures with reasonable accuracy.     

Computation of forces acting on bodies in plane and axisymmetric cavitation flow problems

Plane and axisymmetric cavitation flow problems are considered using Riabouchinsky’s scheme. The incoming flow is assumed to be irrotational and steady, and the fluid is assumed to be inviscid and incompressible. The flow problems are solved by applying the boundary element method with quadrature formulas without saturation. The free boundary is determined using a gradient descent technique based on Riabouchinsky’s principle. The drag force acting on the cavitator is expressed in terms of the Riabouchinsky functional. As a result, for small cavitation numbers, the force is calculated with a fairly high accuracy. Dependences of the drag coefficient are investigated for variously shaped cavitators: a wedge, a cone, a circular arc, and a spherical segment.