Laminar free convection heat transfer from isothermal convex bodies of arbitrary shape: a new dynamic model

Calculation of free convection from bodies of arbitrary shape has been investigated previously. The Body Gravity Function (BGF) which accounts for the geometry of each body shape was considered to be a constant value. In the present study, it is shown that BGF is not a constant value in a wide range of Rayleigh number. Instead, its value changes as Rayleigh number increases. Therefore, by analytical modeling of Dynamic BGF and derivation of a new parameter called Body Fluid Function, a novel method is proposed to calculate laminar free convection heat transfer from isothermal convex bodies of arbitrary shape. Results for 24 different body shapes are compared with the available experimental and numerical data. Excellent agreement shows that the present simple method accurately predicts laminar free convection heat transfer from isothermal convex bodies of arbitrary shape in the whole range of laminar flow and for fluids of any Prandtl number.     

Nichtisotherme ebene Freistrahlen unter Schwerkrafteinfluß

Higher-order boundary layer theory is used to study the behaviour of nonisothermal laminar and turbulent free jet flows. In addition to the Prandtl boundary layer equations, an equation is used to describe the equilibrium of forces normal to the flow direction. This equilibrium exists between the buoyancy forces caused by gravity and the centrifugal forces resulting from the curvature in the flow. The proper selection of reference values permits the characteristics of the jet flow to be expressed as universal functions in which only the initial jet orientation and the Prandtl number in the case of laminar flow are input parameters. When the volume flow is given in addition to the momentum and thermal energy, the characteristic parameter are the Archimedes number for turbulent flow and the modified Archimedes number for laminar flow. The jet flow is calculated using an integral method in which the eddy viscosity and the turbulent Prandtl number are given as functions of the local Archimedes number. Comparison of experimental data from the literature and from our laboratory on nonisothermal free jets with the theoretical results, show satisfactory agreement. The universal diagrams given in the paper are valid forall plane laminar (Pr=0.7) and turbulent nonisothermal jets.     

Numerical investigation of three-dimensional viscous incompressible flows in divergent curved channels and turbulent model study

In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density- weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flows in divergent curved channels was carried out by means of parabolized computation method. The calculating results were used to analyze and investigate the aerodynamic performance of stator cascades in compressors preliminarily.     

TVD格式在超音速喷管三维粘性流动求解中的应用

详细给出了任意三维曲线坐标系中Ｎｏｖｉｅｒ－Ｓｔｏｋｅｓ方程的对流项ＴＶＤ格式的构造过程，建立了数值求解三维粘性流动的计算方法，应用该方法对三维超音速喷管中有激波及无激波情况下的两种工况的层流流场进行了数值求解，并与实验做了对比。结果表明本文建立的计算方法具有较高的精度，同时也证明ＴＶＤ格式具有分辩率高，稳定收敛等优点，为进一步开展叶栅流场及紊流的研究打下了基础。     

Water hammer prediction and control: the Green’s function method

By Green’s function method we show that the water hammer (WH) can be analytically predicted for both laminar and turbulent flows (for the latter,with an eddy viscosity depending solely on the space coordinates),and thus its hazardous effect can be rationally controlled and minimized.To this end,we generalize a laminar water hammer equation of Wang et al.(J.Hydrodynamics,B2,51,1995) to include arbitrary initial condition and variable viscosity,and obtain its solution by Green’s function method.The predicted characteristic WH behaviors by the solutions are in excellent agreement with both direct numerical simulation of the original governing equations and,by adjusting the eddy viscosity coefficient,experimentally measured turbulent flow data.Optimal WH control principle is thereby constructed and demonstrated.     

Turbulent free convection heat transfer to power-law fluids from arbitrary geometric configurations

The problem of turbulent free convection heat transfer from curved surfaces to non-Newtonian power-law fluids has been investigated using the Nakayama-Koyama solution methodology. The scheme is designed to deal with bodies of arbitrary geometric configurations and hence can be viewed as a generalized version of the Shenoy-Mashelkar approach for turbulent free convection heat transfer from a flat vertical plate to a power-law fluid. The surface wall temperature is allowed to vary in the streamwise direction in an arbitrary fashion, and calculations are carried out for the turbulent free convection about the horizontal circular cylinder and sphere for illustrative purposes. Available theoretical and experimental data have been compared with the predictions of the present analysis and the comparison of results has been found to be reasonably good.     

Prediction of periodic boundary layers

A relatively simple, yet efficient and accurate finite difference method is developed for the solution of the unsteady boundary layer equations for both laminar and turbulent flows. The numerical procedure is subjected to rigorous validation tests in the laminar case, comparing its predictions with exact analytical solutions, asymptotic solutions, and/or experimental results. Calculations of periodic laminar boundary layers are performed from low to very high oscillation frequencies, for small and large amplitudes, for zero as well as adverse time-mean pressure gradients, and even in the presence of significant flow reversal. The numerical method is then applied to predict a relatively simple experimental periodic turbulent boundary layer, using two well-known quasi-steady closure models. The predictions are shown to be in good agreement with the measurements, thereby demonstrating the suitability of the present numerical scheme for handling periodic turbulent boundary layers. The method is thus a useful tool for the further development of turbulence models for more complex unsteady flows.     

A PRESSURE CORRECTION METHOD FOR UNSTRUCTURED MESHES WITH ARBITRARY CONTROL VOLUMES

A pressure correction procedure for general unstructured meshes is presented. It is a cell-centred, collocated finite volume method and the pressure–velocity coupling is treated using SIMPLEC. The cells can have an arbitrary number of grid points (cell vertices). In the present study the number of faces on the cells varies between three and six. The discretized equations are solved using either a symmetric Gauss–Seidel solver or a conjugate gradient solver with a preconditioner. The method is applied to three two-dimensional test cases in which the flow is incompressible and laminar. The extension to three dimensions as well as to turbulent flow using transport models is straightforward. It can also be extended to handle compressible flow.     

Study on correction coefficients of liminar and turbulent entrance region effect in round pipe

In this paper, a universal method for calculation of the inlet length and correction coefficient for the pressure losses and the flow in the case of laminar and turbulent flow in a smooth round pipe is put forward.With the aid of concrete examples, theoretical formulas are presented for calculation of the correction coefficient for pressure losses and flow in the case of laminar and turbulent flow. And a method for calculation of the flow is also presented here in consideration of the inlet effect of turbulent flow.By comparing the theoretical results with the experimental data, formulas presented in this paper are very simple and reliable.     

Steady flow separation along a straight streamwise corner

This note considers the steady incompressible flow of a laminar or turbulent boundary layer along a straight streamwise corner, formed by two flat plates intersecting each other at an arbitrary included angle. It is assumed that the undisturbed freestream is parallel or roughly parallel to the cornerline and that the intersecting corner-walls are impermeable. Attention is focussed on the necessary and sufficient conditions for steady flow separation from the sharp cornerline.     

Turbulent boundary layer on a catalytic surface in a nonequilibrium dissociating gas

The majority of the studies which consider the flow of a dissociating gas in a turbulent boundary layer are devoted to the investigation of either frozen or equilibrium flows on a flat plate.The frozen turbulent boundary layer has been studied by Dorrance [1], Kutateladze and Leont'ev [2], and Lapin and Sergeev [3]. A study of the effect of catalytic recombination processes at the plate surface on the heat transfer in a frozen turbulent boundary layer was made by Lapin [4].Kosterin and Koshmarov [5], Ginzburg [6], Dorrance [7], and Lapin [8] have studied the turbulent boundary layer on a plate in equilibrium dissociating gas.The calculation of the heat transfer in a turbulent boundary layer on a catalytic plate surface with nonequilibrium dissociation was made by Kulgein [9]. In this study the nonequilibrium nature of the dissociation process was taken into account only in the laminar sublayer, while the flow in the turbulent core was considered frozen. The solution was found numerically using a computer by means of a laborious iteration process.The present paper reports a method for calculating the turbulent boundary layer on a flat catalytic plate with arbitrary dissociation rate. The method, constructed using the assumptions customary for turbulent boundary layer theory, is a successive approximation method. Good convergence of the method is assured by the fact that the effect of the nonequilibrium nature of the dissociation process on the parameter distribution in the boundary layer and, consequently, on the friction and heat transfer may be allowed for merely by finding corrections, usually relatively small, to the distribution of these parameters in the equilibrium or frozen flows. The basis of the study is the two-layer scheme of the turbulent boundary layer. The Prandtl and Schmidt numbers and also their turbulent analogs are taken equal to unity. As the model of the dissociating gas we use the Lighthill model of the ideal dissociating gas [10], extended by Freeman [11] to nonequilibrium flows.     

An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A matrix‐free preconditioned Newton/GMRES method for unsteady Navier–Stokes solutions

The unsteady compressible Reynolds‐averaged Navier–Stokes equations are discretized using the Osher approximate Riemann solver with fully implicit time stepping. The resulting non‐linear system at each time step is solved iteratively using a Newton/GMRES method. In the solution process, the Jacobian matrix–vector products are replaced by directional derivatives so that the evaluation and storage of the Jacobian matrix is removed from the procedure. An effective matrix‐free preconditioner is proposed to fully avoid matrix storage. Convergence rates, computational costs and computer memory requirements of the present method are compared with those of a matrix Newton/GMRES method, a four stage Runge–Kutta explicit method, and an approximate factorization sub‐iteration method. Effects of convergence tolerances for the GMRES linear solver on the convergence and the efficiency of the Newton iteration for the non‐linear system at each time step are analysed for both matrix‐free and matrix methods. Differences in the performance of the matrix‐free method for laminar and turbulent flows are highlighted and analysed. Unsteady turbulent Navier–Stokes solutions of pitching and combined translation–pitching aerofoil oscillations are presented for unsteady shock‐induced separation problems associated with the rotor blade flows of forward flying helicopters. Copyright © 2000 John Wiley & Sons, Ltd.     

Visualizing supersonic inlet duct unstart using planar laser Rayleigh scattering

Planar laser Rayleigh scattering (PLRS) from condensed CO₂ particles is used to visualize flow structure in a Mach 5 wind tunnel undergoing unstart. Detailed flow features such as laminar/turbulent boundary layers and shockwaves are readily illustrated by the technique. A downstream transverse air jet, inducing flow unchoking downstream of the jet, is injected into the free stream flow of the tunnel, resulting in tunnel unstart. Time sequential PLRS images reveal that the boundary layer growth/separation on a surface with a thick turbulent boundary layer, initiated by the jet injection, propagates upstream and produces an oblique unstart shock. The tunnel unstarts upon the arrival of the shock at the inlet. In contrast, earlier flow separation on the opposite surface, initially supporting a thin laminar boundary layer, is observed when a jet induced bow shock strikes that surface. The resulting disturbance to this boundary layer also propagates upstream and precedes the formation of an unstart shock.     

Forced convective film condensation inside vertical tubes

A theoretical study of forced convective film condensation inside vertical tubes is presented. We propose a unified procedure for predicting the pressure gradient and condensation heat transfer coefficient of a vapor flowing turbulently in the core and associated with laminar or turbulent film on the tube wall. The analysis for the vapor flows is performed under the condition that the velocity profiles are locally self-similar. The laminar and turbulent film models equate the gravity, pressure and viscous forces, and consider the effect of interfacial shear. The transition from laminar to turbulent film depends not only on the liquid Reynolds number but also on the interfacial shear stress. In this work we also proposed a new eddy viscosity model which is divided into three regions: the inner region in liquid condensate near the wall; the interface region including both liquid and vapor; and the outer region for the vapor core. Comparisons of the theory with some published experimental data showed good agreement.     

Time-dependent liquid metal flows with free convection and a deformable free surface

The finite element method is employed to investigate time-dependent liquid metal flows with free convection, free surfaces and Marangoni effects. The liquid circulates in a two-dimensional shallow trough with differentially heated vertical walls. The spatial formulation incorporates mixed Lagrangian approximations to the velocity, pressure, temperature and free surface position. The time integration is performed with the backward Euler and trapezoid rule methods with step size control. The Galerkin method is used to reduce the problem to a set of non-linear equations which are solved with the Newton–Raphson method. Calculations are performed for conditions relevant to the electron beam vaporization of refractory metals. The Prandtl number is 0·015 and Grashof number are in the transition range between laminar and turbulent flow. The results reveal the effects of flow intensity, surface tension gradients, mesh refinement and time integration strategy.     

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.     

考虑法向压力梯度的三维旋转边界层计算

本文以量级分析为基础,建立了一般曲线坐标系上的三维旋转边界层方程。对旋转在边界层中的影响进行分析之后,提出了一个能够处理壁面法向压力梯度不为零问题的压力梯度迭代方法。在传统的Box法的基础上发展了一套完整的求解三维旋转边界层的方法和程序,并对螺旋面、压气机转子叶面以及圆柱面上的旋转边界层进行了计算,与他人的计算和实验的对比分析表明,该方法和程序是正确的,可用于求解任意几何物面上的三维旋转边界层。     

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.     

Conjugated laminar forced convection in ducts with periodic variation of inlet temperature

Conjugated laminar forced convection with parabolic velocity is studied for flow inside both parallel-plate and circular ducts subjected to a periodically varying inlet temperature. The analysis of this classs of problems lead to a Sturm-Liouville type complex eigenvalue problem for which no known solution is available. In this work, a new methodology is presented for a direct solution of such complex eigenvalue problems in the complex domain using the shooting method along with the Runge-Kutta method. The methodology is applicable for solving both laminar and turbulent flow problems; here we consider only the laminar flow with a parabolic velocity profile. The results clearly show that the slug-flow assumption overestimates the phase lag and the Nusselt number. Some benchmark results are also presented for the eigenquantities in tabular form.