The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

Three-dimensional analysis of the flow in the cyclones

Sincethecycloneisextensivelyusedinindustry,themechanismofitsstreamhasbeeninvestigatedandanalysed.Buttheprocessofthepreviousdiscussionstaysattheleveloftwo-dimensionedmodelanalysis,withtheaddedconsiderationabouttheaxialIysymmetry.lntheseyears,somespecialistshavediscusseditinthelightofthetotalthree-dimensionview.Thispapertriestogivesystematicanddetai]edexp1anationsandtopresentacorrectresultofRef[1].I.GeneralEquationEstablishedFirstweestablishmathematicalexpressionofflowinthecyclone.Wecanimaget…     

Verification and validation of impinging round jet simulations using an adaptive FEM

This paper illustrates the use of an adaptive finite element method as a means of achieving verification of codes and simulations of impinging round jets, that is obtaining numerical predictions with controlled accuracy. Validation of these grid‐independent solution is then performed by comparing predictions to measurements. We adopt the standard and accepted definitions of verification and validation (Technical Report AIAA‐G‐077‐1998, American Institute of Aeronautics and Astronautics, 1998; Verification and Validation in Computational Science and Engineering. Hermosa Publishers: Albuquerque, NM, 1998). Mesh adaptation is used to perform the systematic and rigorous grid refinement studies required for both verification and validation in CFD. This ensures that discrepancies observed between predictions and measurements are due to deficiencies in the mathematical model of the flow. Issues in verification and validation are discussed. The paper presents an example of code verification by the method of manufactured solution. Examples of successful and unsuccessful validation for laminar and turbulent impinging jets show that agreement with experiments is achieved only with a good mathematical model of the flow physics combined with accurate numerical solution of the differential equations. The paper emphasizes good CFD practice to systematically achieve verification so that validation studies are always performed on solid grounds. Copyright © 2004 John Wiley & Sons, Ltd.     

纳米颗粒多相流体动力学研究及应用

纳米颗粒多相流研究是目前多相流研究中新的研究方向及重点发展领域. 为探索纳米尺度多相流相间作用机理及内部存在机制,采用理论分析及数值计算手段,对一般动力学方程的封闭处理、颗粒碰撞率宏观模型的有效构建、颗粒凝并系统动力学演变特性的机理分析、非稀相问题碰撞率的求取、双变量问题求解方法的建立以及一些实际应用进行了系统研究,提出了新的针对纳米尺度颗粒动力学演变的一般动力学方程求解方法,并将其应用于实际工业过程问题的研究. 该文对上述研究工作进行了综述.     

Attenuation analysis of hydraulic transients with laminar-turbulent flow alternations

An improved compound mathematical model is established to simulate the attenuation of hydraulic transients with laminar-turbulent alternations,which usually occur when the pipeline flow velocity fluctuates near the critical velocity.The laminar friction resistance and the turbulent friction resistance are considered respectively in this model by applying different resistance schemes to the characteristics method of fluid transient analysis.The hydraulic transients of the valve closing process are simulated using the model.A more reasonable attenuation of hydraulic transients is obtained.The accurate attenuation is more distinct than that obtained from the traditional mathematical model.The research shows that the hydraulic transient is a type of energy waves,and its attenuation is caused by the friction resistance.The laminar friction resistance is more important than the turbulent friction resistance if the flow velocity is smaller than the critical velocity.Otherwise the turbulent friction resistance is more important.The laminar friction resistance is important in the attenuation of hydraulic transients for the closing process.Thus,it is significant to consider the different resistances separately to obtain more accurate attenuation of hydraulic transients.     

Extension of some numerical schemes to the analysis of gas and particle mixtures

In this paper, several numerical schemes are extended to obtain approximate solutions to the system of equations encountered in the analysis of multiphase mixtures of gas and particles. Both dense and dilute mixtures are studied, the gas is modelled as a perfect gas and the solid is considered incompressible. Although the tests employed throughout this work for studying the behaviour of the schemes are essentially one dimensional, the finite volume method developed permits its application to multidimensional problems in unstructured grids. Copyright © 2007 John Wiley & Sons, Ltd.     

MHD squeezing flow of second‐grade fluid between two parallel disks

This paper examines the unsteady two‐dimensional flow of a second‐grade fluid between parallel disks in the presence of an applied magnetic field. The continuity and momentum equations governing the unsteady two‐dimensional flow of a second‐grade fluid are reduced to a single differential equation through similarity transformations. The resulting differential system is computed by a homotopy analysis method. Graphical results are discussed for both suction and blowing cases. In addition, the derived results are compared with the homotopy perturbation solution in a viscous fluid (Math. Probl. Eng., DOI: 10.1155/2009/603916 ). Copyright © 2011 John Wiley & Sons, Ltd.     

Non-isothermal polymer melt flow in sudden expansions

This paper presents numerical results for laminar, incompressible and non-isothermal polymer melt flow in sudden expansions. The mathematical model includes the mass, momentum and energy conservation laws within the framework of a generalized Newtonian formulation. Two constitutive relations are adopted to describe the non-Newtonian behavior of the flow, namely Cross and Modified Arrhenius Power-Law models. The governing equations are discretized using the finite difference method based on central, second-order accurate formulas for both convective and diffusive terms. The pressure–velocity coupling is treated by solving a Poisson equation for pressure. The results are presented for two commercial polymers and demonstrate that important flow parameters, such as pressure drop and viscosity distribution, are strongly affected by heat transfer features.     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

Visualisation and modelling studies of churn flow in a vertical pipe

Churn flow is an important flow regime intermediate between slug flow and annular flow. A feature of this regime is the occurrence of very large waves travelling upwards over a liquid film substrate which may intermittently travel downwards. These waves are often formed close to the liquid inlet where their behaviour is usually difficult to observe. This paper describes a series of experiments using a test section with a specially constructed transparent liquid inlet. High-speed video recordings show clearly the process of wave formation and analysis of the recordings gave data on wave frequencies and typical velocities. Also, predictions of velocity and distance travelled by the waves were obtained via the application of a simple mathematical model that takes into account the forces acting on a circumferentially coherent wave.     

Evaluation of Smagorinsky‐based subgrid‐scale models in a finite‐volume computation

Smagorinsky‐based models are assessed in a turbulent channel flow simulation at Re_b=2800 and Re_b=12500. The Navier–Stokes equations are solved with three different grid resolutions by using a co‐located finite‐volume method. Computations are repeated with Smagorinsky‐based subgrid‐scale models. A traditional Smagorinsky model is implemented with a van Driest damping function. A dynamic model assumes a similarity of the subgrid and the subtest Reynolds stresses and an explicit filtering operation is required. A top‐hat test filter is implemented with a trapezoidal and a Simpson rule. At the low Reynolds number computation none of the tested models improves the results at any grid level compared to the calculations with no model. The effect of the subgrid‐scale model is reduced as the grid is refined. The numerical implementation of the test filter influences on the result. At the higher Reynolds number the subgrid‐scale models stabilize the computation. An analysis of an accurately resolved flow field reveals that the discretization error overwhelms the subgrid term at Re_b=2800 in the most part of the computational domain. Copyright © 2002 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

On the global asymptotic stability of periodic solutions of a class of ecological system

In this paper, we have made reaearches on the mathematical models which have three populations of mutual action: We have obtained the sufficient conditions respectively for the systems(*)and(**)for existence and uniqueness of single positive periodic solutions which are globally asymptotically stable.     

Similar solution for three‐dimensional flow in an Oldroyd‐B fluid over a stretching surface

The present investigation deals with the three‐dimensional flow of an Oldroyd‐B fluid over a stretching surface. The governing equations for the three‐dimensional flow are developed. Similarity transformations are invoked for the conversion of nonlinear partial differential equations into the coupled system of ordinary differential equations. Computations for the series solution are presented through implementation of homotopy analysis method. The salient features of the involved parameters have been pointed out. Copyright © 2011 John Wiley & Sons, Ltd.     

An analytic solution of dense two-phase flow in a vertical pipeline

According to a mathematical model for dense two-phase flows presented in theprevious paper,a dense two-phase flow in a vertical pipeline is analytically solved,and theanalytic expressions of velocity of each continuous phase and dispersed phase arerespectively derived The results show that when the drag force between two phases dependslinearly on their relative velocity,the relative velocity profile in the pipeline coincides withDarcy’s law except for the thin layer region near the pipeline wall,and that the theoreticalassumptions in the dense two-phase flow theory mentioned are reasonable.     

Modelling a recirculating density-driven turbulent flow

A mathematical model of turbulent density-driven flows is presented and is solved numerically. A form of the k–? turbulence model is used to characterize the turbulent transport, and both this non-linear model and a sediment transport equation are coupled with the mean-flow fluid motion equations. A partitioned, Newton–Raphson-based solution scheme is used to effect a solution. The model is applied to the study of flow through a circular secondary sedimentation basin.     

Three‐dimensional flow of upper‐convected Maxwell (UCM) fluid

An analysis has been carried out for three‐dimensional fluid over a stretching surface. Similarity transformations are invoked for the conversion of nonlinear partial differential equations to the ordinary differential equations. Computations for the series solution are made by using homotopy analysis method. Graphical results are obtained. Attention has been particularly given to the influence of Deborah number. Copyright © 2010 John Wiley & Sons, Ltd.