Assessment of an improved turbulence model in simulating the unsteady flows around a D-shaped cylinder and an open cavity

Most engineering flows are still predicted by the conventional Reynolds-averaged Navier-Stokes method because of the low requirements of the computational quantities. However, the resolution capability of Reynolds-averaged Navier-Stokes models is still open to deliberation, especially in the recirculation and wake regions, where the vortical flows dominate. In the present work, an improved turbulence model derived from the original shear stress transport k-ω model is proposed and its superiority is assessed by our modeling the unsteady flows around a D-shaped cylinder and an open cavity, corresponding to two different Reynolds numbers. The results are compared with results from experiments and other turbulence models in terms of the flow morphology and mean velocity profiles. This shows that the predictive accuracy of the modified turbulence model is increased significantly in the bluff body wake flows and in the shear layer and separation flows of the cavity. Some special vortex structures can be captured in the open cavity, in which the secondary vortex emerging from the shear layer and the separation vortex near the trailing edge can induce large flow instability, and this phenomenon should be eliminated in engineering applications. It is believed that this improved turbulence model can be used for the more complex turbomachinery flows with better prediction of the hydrodynamic/aerodynamic performance and the unsteady vortical flows, which can provide some guidelines to design or optimize rotating machines.     

Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

The effect of gravity on flow past a semi-circular cylinder with a constant pressure wake

The flow past a semi-cylinder with a trailing wake region is considered. In the absence of gravity the only known high Reynolds number solutions have tangential separation from the body and a cusped shape at the back of the wake. This flow can be a simple model for several situations, including the classical approximation of a constant pressure wake and the flow past an object with a region of trapped fluid of different density (or an air cavity) attached on the downstream side. Here we relax the assumption of high flow speeds to examine the effects of gravity. It is shown that there are situations in which a stagnation point can form either on the body or at the tail of the wake and that there is a minimum velocity beneath which a cavity will not form. Non-uniqueness in the parameter space is found in certain cases.     

Dual solutions in mixed convection boundary layer flow of micropolar fluids

The steady mixed convection boundary layer flow over a vertical surface immersed in an incompressible micropolar fluid is considered in this paper. Employing suitable similarity transformations, the governing partial differential equations are transformed into ordinary differential equations, and the transformed equations are solved numerically by the Keller-box method. Numerical results are obtained for the skin friction coefficient and the local Nusselt number as well as the velocity, angular velocity and temperature profiles. Both cases of assisting and opposing buoyant flows are considered. It is found that dual solutions exist for the assisting flow, besides that usually reported in the literature for the opposing flow. Moreover, in contrast to the classical boundary layer theory, the separation point of the boundary layer is found to be distinct from the point of vanishing skin friction.     

Species separation in narrow gaps—a simple model

The problem of the separation of Uranium isotopes from the fluid-dynamic point of view is considered. The technique based on the separator nozzle effects is discussed. As this system works at low pressure, the analysis is performed in the slip flow regime. The boundary layer, transition and asymptotic regions for the velocity and concentration profiles in the nozzle are studied.     

Mixed convection on a vertical circular cylinder

The problem of the mixed convection boundary-layer flow past an isothermal vertical circular cylinder is considered in both the cases when the buoyancy forces aid and oppose the development of the boundary layer. A series solution is obtained, valid near the leading edge, and this is extended by a numerical solution of the full equations, which in the aiding case, becomes inaccurate downstream. An approximate solution is also derived which gives a good estimate for the heat transfer near the leading edge and has the correct asymptotic form well downstream. In the opposing case, the boundary layer is seen to separate at a finite distance downstream, with, for moderate values of the buoyancy parameter, the numerical solution indicating a regular behaviour near separation.     

Nonequilibrium jet boundary layer in a polyatomic gas

The two-dimensional nonequilibrium hypersonic free jet boundary layer gas flow in the near wake of a body is studied using a closed system of macroscopic equations obtained (as a thin-layer version) from moment equations of kinetic origin for a polyatomic single-component gas with internal degrees of freedom. (This model is can be used to study flows with strong violations of equilibrium with respect to translational and internal degrees of freedom.) The solution of the problem under study (i.e., the kinetic model of a nonequilibrium homogeneous polyatomic gas flow in a free jet boundary layer) is shown to be related to the known solution of the well-studied simpler problem of a Navier-Stokes free jet boundary layer, and a method based on this relation is proposed for solving the former problem. It is established that the gas flow velocity distribution along the separating streamline in the kinetic problem of a free jet boundary layer coincides with the distribution obtained by solving the Navier-Stokes version of the problem. It is found that allowance for the nonequilibrium nature of the flow with respect to the internal and translational degrees of freedom of a single-component polyatomic gas in a hypersonic free jet boundary layer has no effect on the base pressure and the wake angle.     

三维横向流体诱发直管振动的数值模拟

基于有限体积法和有限元法,结合动网格控制技术,建立了横向流体作用下三维弹性直管流致振动计算的数值模型,实现了计算结构动力学与计算流体力学之间的联合仿真．首先,通过对刚性管的静止绕流计算,研究了网格离散方式和不同湍流模型对圆柱类结构静止绕流流场特征的影响和预测能力,得到了适用于双向耦合分析的CFD模型;其次,利用基于双向流固耦合方法的流致振动模型,计算并分析了流体力与结构位移间的相位关系,指出流体力与位移间的相位差是由流体力引起的,同时对双向耦合和单向耦合进行了比较分析;最后通过对直管流致振动的数值计算,联合管表面压力、尾流区时均速度、分离角等时均量,分析了尾流区的流场特征.     

A second order boundary layer solution for separated flow past a bluff body

In this investigation second order boundary layer solutions are obtained for subsonic, plane flow over a circular cylinder with a wake. The outer flow due to displacement thickness is obtained by extending a first order solution, which includes wake effects. This enables realistic solutions for the displacement effect to be calculated as well as for the curvature effect. It is shown that as the first order separation point is approached the magnitudes of the second order velocity components become large and the surface shear becomes large and negative. There is also a significant upstream shift in the position of the separation point.

---

Résumé Dans ce travail, on obtient des solutions du second ordre des équations de la couche limite pour un écoulement subsonique plan autour d'un cylindre circulaire avec un sillage. On obtient l'écoulement extérieur dû à l'épaisseur de déplacement en étendant une solution du premier ordre qui comprend les effets du sillage. Ceci permet de calculer des solutions réalistes pour l'effet de déplacement ainsi que pour l'effet de courbure. On montre que les composantes de la vitesse du second ordre de même que la tension tangentielle deviennent très grandes et cette dernière négative au voisinage du point de décollement du premier ordre. Il y a aussi un déplacement significatif vers l'amont de la position du point de décollement.

     

Nonsimilar incompressible laminar boundary layers with magnetic field

The solution of the steady laminar incompressible nonsimilar magneto-hydrodynamic boundary layer flow and heat transfer problem with viscous dissipation for electrically conducting fluids over two-dimensional and axisymmetric bodies with pressure gradient and magnetic field has been presented. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for flow over a cylinder and a sphere. The results indicate that the magnetic field tends to delay or prevent separation. The heat transfer strongly depends on the viscous dissipation parameter. When the dissipation parameter is positive (i.e. when the temperature of the wall is greater than the freestream temperature) and exceeds a certain value, the hot wall ceases to be cooled by the stream of cooler air because the ‘heat cushion’ provided by the frictional heat prevents cooling whereas the effect of the magnetic field is to remove the ‘heat cushion’ so that the wall continues to be cooled. The results are found to be in good agreement with those of the local similarity and local nonsimilarity methods except near the point of separation, but they are in excellent agreement with those of the difference-differential technique even near the point of separation.     

Receptivity of the boundary layer with separation bubble to external sound waves

The receptivity problem of laminar separated two-dimensional boundary-layer under influence of external acoustic waves is considered. Basic features of this theory were formulated by Ruban (1984) and then Goldstein (1985), but these investigations were limited by small perturbations of steady parallel flow only. The current paper takes into account essentially nonparallel regimes of base flow, including attached flows, flow with separated bubbles and the base flow at marginal value of parameter. The steady flow near a corner point of profile is of considerable interest and a question naturally arises about the sensitivity of separation bubble inside the boundary layer to external disturbances such as sound waves. Another question arises about the sensitivity of near critical base flow, namely if corner is concave there is a marginal angle which limits existence of such base flow. This fact usually is linked with sudden reconstruction of full flow pattern past a profile. Development of perturbations inside the interaction region may give some answers about the appearance of turbulence. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Long Eddies in Sheared Flows

The characteristic feature of the wide variety of hydraulic shear flows analyzed in this study is that they all contain a critical level where some of the fluid is turned relative to the ambient flow. One example is the flow produced in a thin layer of fluid, contained between lateral boundaries, during the passage of a long eddy. The boundaries of the layer may be rigid, or flexible, or free; the fluid may be either compressible or incompressible. A further example is the flow produced when a shear layer separates from a rigid boundary producing a region of recirculating flow. The equations used in this study are those governing inviscid hydraulic shear flows. They are similar in form to the classical boundary layer equations with the viscous term omitted. The main result of the study is to show that when the hydraulic flow is steady and contained between lateral boundaries, the variation of vorticity ω(ψ) cannot be prescribed at any streamline which crosses the critical level. This variation is, in fact, determined by (1) the vorticity distribution at all streamlines which do not cross the critical level, by (2) the auxiliary conditions which must be satisfied at the boundaries of the fluid layer, and by (3) the dimensions of the region containing the turned flow. If at some instant the vorticity distribution is specified arbitrarily at all streamlines, generally the subsequent flow will be unsteady. In order to emphasize this point, a class of exact solutions describing unsteady hydraulic flows are derived. These are used to describe the flow produced by the passage of a long eddy which distorts as it is convected with the ambient flow. They are also used to describe the unsteady flow that is produced when a shear layer separates from a boundary. Examples are given both of flows in which the shear layer reattaches after separation and of flows in which the shear layer does not reattach. When the shear layer vorticity distribution has the form ωαyⁿ, where y is a distance measure across the layer, the steady flows are of Falkner-Skan type inside, and adjacent to, the separation region. The unsteady flows described in this paper are natural generalizations of these Falkner-Skan flows. One important result of the analysis is to show that if the unsteady flow inside the separation region is strongly sheared, then the boundary of the separation region moves upstream towards the point of separation, forming large transverse currents. Generally, the assumption of hydraulic flow becomes invalid in a finite time. On the other hand, if the flow inside the separation region is weakly sheared, this region is swept downstream and the flow becomes self-similar.     

Analytical solutions of laminar swirl decay in a straight pipe

In this work, the laminar swirl flow in a straight pipe is revisited and solved analytically by using prescribed axial flow velocity profiles. Based on two axial velocity profiles, namely a slug flow and a developed parabolic velocity profiles, the swirl velocity equation is solved by the separation of variable technique for a rather general inlet swirl velocity distribution, which includes a forced vortex in the core and a free vortex near the wall. The solutions are expressed by the Bessel function for the slug flow and by the generalized Laguerre function for the developed parabolic velocity. Numerical examples are calculated and plotted for different combinations of influential parameters. The effects of the Reynolds number, the pipe axial distance, and the inlet swirl profiles on the swirl velocity distribution and the swirl decay are analyzed. The current results offer analytical equations to estimate the decay rate and the outlet swirl intensity and velocity distribution for the design of swirl flow devices.     

Solution of the Falkner–Skan equation for wedge by Adomian Decomposition Method

The Adomian Decomposition Method is employed in the solution of the two dimensional laminar boundary layer of Falkner–Skan equation for wedge. This work aims at the solution of momentum equation in the case of accelerated flow and decelerated flow with separation. The Adomian Decomposition Method is provided an analytical solution in the form of an infinite power series. The effect of Adomian polynomials terms is considered on accuracy of the results. The velocity profiles in boundary layer are obtained. Results show a good accuracy compared to the exact solution.     

Turbulent Poiseuille flow with near-critical wall transpiration

An incompressible, pressure-driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress with the first and the second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary-value problem for a second-order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. The case of near-critical transpiration, when the shear stress at the injection wall vanishes, is considered. It is shown that the maximum point on the mean velocity profile lies in a thin sublayer near the suction wall in this case. A formula for the position of the maximum point as a function of the transpiration factor is obtained. The mean velocity profiles near the suction wall are calculated. A friction law for Poiseuille flow with near-critical transpiration is found, which makes it possible to describe the relation between the shear stress at the wall, the Reynolds number, and the transpiration velocity by a single function of one variable. Direct numerical simulation of the flow for some transpiration factors is performed. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of laminar flow past a circular cylinder

The present paper focuses on the analysis of two- and three-dimensional flow past a circular cylinder in different laminar flow regimes. In this simulation, an implicit pressure-based finite volume method is used for time-accurate computation of incompressible flow using second order accurate convective flux discretisation schemes. The computation results are validated against measurement data for mean surface pressure, skin friction coefficients, the size and strength of the recirculating wake for the steady flow regime and also for the Strouhal frequency of vortex shedding and the mean and RMS amplitude of the fluctuating aerodynamic coefficients for the unsteady periodic flow regime. The complex three dimensional flow structure of the cylinder wake is also reasonably captured by the present prediction procedure.     

Shock Interactions in Nonequilibrium Hypersonic Flow

A shock interaction problem is solved with finite difference methods for a hypersonic flow of air with chemical reactions. If a body has two concave corners, a secondary shock is formed in the shock layer and it meets the main shock later. As the two shocks meet, the flow becomes singular at the interaction point, and a new main shock, a contact discontinuity and an expansion wave appear as a result of interaction between the two shocks. Therefore, the problem is very complicated. Using proper combinations of implicit and explicit finite difference schemes according to the property of the equations and the boundary conditions, we compute the flow behind the interaction point successfully.     

粘弹性流体流经竖直表面时的瞬时自然对流

静止流体中,在一个竖直的、不可渗透的等温表面附近,研究粘弹性边界层的流动及其热传导.得到其控制方程,并利用MackCormak技术对其进行数值求解.与先前发表的关于该问题特例的结果相比较,有着很好的一致性.对于不同的粘弹性参数值,图示了速度和温度分布、边界层厚度、Nusselt数、局部摩擦因数等典型结果.一般而言,粘弹性流体与Newton流体相比较,由于拉应力的促进作用,流体动力边界层里的速度是增加的,热边界层里的温度是下降的.粘弹性参数值越高,摩擦因数和传热系数越高.     

Simulation of wake effects of windfarms using an Actuator Disk implementation

Wind turbines extract energy from the approaching flow field resulting in reduced wind speeds, increased turbulence and a wake downstream of the wind turbine. Wakes extending up to 20 km downstream of a windfarm are reported [1]. In [2] a multi-scale numerical method is proposed that addresses length scales of the boundary layer, windfarm and meso-scale. As a step towards full realization of the multi scale approach the implementation of an Actuator Disk (AD) model into OpenFOAM including a proper atmospheric boundary layer (ABL) is presented. For a meso-scale area of open sea, the standard ABL conditions of OpenFOAM can not be used. The pre implemented ABL boundary condition for OpenFOAM is restricted to roughness lengths greater than 0.001 m. The boundary condition are modified in order to achieve adequate homogeneous flow properties throughout the considered domain enabling almost constant velocity profiles. A new AD class is implemented in OpenFOAM to establish an arbitrary number and location of windturbines. The class uses the standard k-ϵ turbulence model. The wake effect of windfarms composed of windturbines with different heights and different lateral spacing is investigated. Specifically the windfarm Trianel Borkum in the North Sea of Germany is simulated. Within this work wakes extending up to 20 km in length could be observed, suggesting to study the ABL-windfarm interaction in more detail. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)