Centrifugal Convection in Rapid Rotation of Bodies Made of Cellular-Porous Materials

Gas flows inside and around rapidly rotating bodies made of cellular-porous materials are studied numerically and experimentally. Within the framework of the previously proposed physicomathematical model, an appropriate numerical algorithm is developed and tested. Internal flows and a conjugate problem with the external flow are considered. The calculated moment and dynamic pressure are in good agreement with experimentally measured characteristics of a rotating porous disk on a solid substrate. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 46–57, January–February, 2006.     

Some unique features of supersonic unsteady gas flow around bodies

We consider unsteady supersonic gas flow about bodies for small Strouhal numbers. The amplitude of the angles of attack is assumed arbitrary under the condition that the bow shock remains attached, and the flow behind the shock is supersonic. A criterion is formulated which permits the comparison of the damping characteristic for small and large amplitudes of the disturbed motion. A comparison is made, using the wedge as an example, of the damping characteristics calculated by various theories, and the hypersonic similarity law is verified. Simple asymptotic equations are presented for the rotary derivatives of the thin wedge in a hypersonic gas stream.     

Aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials

In this study, a thin-walled beam made of functionally graded material (FGM) which is used as rotating blades in turbomachinery under aerothermoelastic loading is investigated. The governing equations, which are based on first-order shear deformation theory, include the effects of the presetting angle, the secondary warping, temperature gradient through the wall thickness of the beam and also the rotational speed. Moreover, quasi-steady aerodynamic pressure loadings are determined using first-order piston theory, and steady beam surface temperature is obtained from gas dynamics theory. Then, the blade partial differential equations are transformed into a set of ordinary differential equations using the extended Galerkin method. Finally, having solved the resulting structural–fluid–thermal eigenvalue system of equations, the effects of Mach number and geometric parameters on natural frequencies are presented. The results demonstrate that the natural frequencies decrease under aerothermoelastic loading at high Mach numbers.     

On the flow of a fluid–particle mixture between two rotating cylinders, using the theory of interacting continua

Mixture theory is used to develop a model for a flowing mixture of solid particulates and a fluid. Equations describing the flow of a two-component mixture consisting of a Newtonian fluid and a granular solid are derived. These relatively general equations are then reduced to a system of coupled ordinary differential equations describing Couette flow between concentric rotating cylinders. The resulting boundary value problem is solved numerically and representative results are presented.     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

On solving the problem of an ideal incompressible fluid flow around large-aspect-ratio axisymmetric bodies

Based on Babenko’s fundamental mathematical ideas, principally new (unsaturated) algorithms are developed for the numerical solution of problems of a potential axisymmetric ideal fluid flow around bodies of revolution, in particular, an ellipsoid of revolution with an aspect ratio equal to 1000. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 56–67, September–October, 2006.     

Internal flow numerical simulation of a cross‐flow fan in refrigerant operating condition using user‐defined functions

In the paper, a cross‐flow fan in refrigerant operating condition is systematically simulated using user‐defined functions. Three‐dimensional simulations are acquired with Navier–Stokes equations coupled with k–ε turbulence model, and internal flow characteristics of an indoor split‐type air conditioner are obtained, which is mainly composed of cross‐flow fan and heat exchanger. It has systematically been simulated in the isothermal flow condition that the performance of cross‐flow fan may be reduced easily with dry or humid air, and in the refrigerant operating condition in which user‐defined functions are applied to the humid air, considered as a mixture of dry air and vapor. A density‐modulated function is adopted to deal with the condensation of the vapor at the heat‐transfer region approximately. The results show flow mechanism of the two gas‐phase flow, including phase‐vary process. The distribution of the parameters is not uniform at the inlet of the machine, the intensity and position of pressure and velocity vary along the axial direction of the fan, the distribution of vapor volume fraction and turbulent intensity in heat‐transfer region is obtained, and the external characteristic data of the indoor machine are obtained and analyzed. Compared with the experimental data, the calculated characteristic curves and designed parameters are on target. © British Crown Copyright 2010/MOD. Reproduced with permission. Published by John Wiley & Sons, Ltd.     

Simulation of rapidly varying flow using an efficient TVD–MacCormack scheme

An efficient numerical scheme is outlined for solving the SWEs (shallow water equations) in environmental flow; this scheme includes the addition of a five‐point symmetric total variation diminishing (TVD) term to the corrector step of the standard MacCormack scheme. The paper shows that the discretization of the conservative and non‐conservative forms of the SWEs leads to the same finite difference scheme when the source term is discretized in a certain way. The non‐conservative form is used in the solution outlined herein, since this formulation is simpler and more efficient. The time step is determined adaptively, based on the maximum instantaneous Courant number across the domain. The bed friction is included either explicitly or implicitly in the computational algorithm according to the local water depth. The wetting and drying process is simulated in a manner which complements the use of operator‐splitting and two‐stage numerical schemes. The numerical model was then applied to a hypothetical dam‐break scenario, an experimental dam‐break case and an extreme flooding event over the Toce River valley physical model. The predicted results are free of spurious oscillations for both sub‐ and super‐critical flows, and the predictions compare favourably with the experimental measurements. Copyright © 2006 John Wiley & Sons, Ltd.     

高超声速非定常流动的数值模拟与气动热计算

高超声速飞行器研究中的一个重点问题是飞行器表面的气动加热,它对飞行器的气动、热特性及安全性有重要的影响.受到当前实验技术的限制,地面实验无法准确模拟真实飞行条件,所以采用数值模拟研究气动加热问题成为目前重要的研究手段.本文采用数值方法求解三维N-S方程,得到钝头体再入模型绕流的瞬态流场,驻点温度及表面热流沿轨道变化规律.计算中采用变边界条件模拟沿轨道飞行的非定常性.     

Large eddy simulation in a fully developed turbulent flow in a channel and comparison of subgrid eddy viscosity models

The accuracy and computational efficiency are compared for a number of models of subgrid eddy viscosity (Smagorinsky model, renormalization group model, and dynamic and one-parameter models). Space-filtered Navier-Stokes equations are solved numerically by the control-volume approach on a nonuniform grid with the use of high-resolution schemes in time and space. The numerical data are compared with the results of a physical experiment and direct numerical simulation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 31–42, May–June, 2006.     

Numerical simulation of a flow around an impulsively started radially deforming circular cylinder

The development of a two‐dimensional viscous incompressible flow generated by a deformable circular cylinder impulsively started into rectilinear motion is studied numerically for the Reynolds numbers equal to 550 and 3000. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the streamfunction is solved by a Fourier–Galerkin method in the direction of the flow that is assumed to remain symmetrical and a second‐order finite difference for the radial direction. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The computed results are compared qualitatively with experimental and numerical results done before in the particular non‐deformable case. The comparison is found to be satisfactory. The influence of the deformation of the cylinder on the flow structure and the drag coefficient is then analyzed. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

基于分区拼接网格技术高升力装置流场数值模拟

针对高升力装置构型模型结构复杂、流场变化剧烈等特点,本文采用分区拼接网格的思想分别按照流场和结构拓扑特点对高升力装置进行了网格分区。在分区的基础上逐块生成网格,减小了增升装置网格的生成难度,提高了网格质量,减少了网格数目。首先,研究了高升力装置的分区策略及流场特点;接着以MD30P-30N多段翼型为研究算例研究了网格比例和插值方法对计算结果的影响。经过分析对比可知:外部区域与近壁区域之间的比例不宜小于1:5;内部域网格比例不宜超过1:1.8,最好保持在1:1左右;计算中应该采用高阶精度插值以保证计算精度。采用某四段翼型进行了验证;最后采用NASA标准高升力装置进行了三维高升力装置流场数值模拟并与相应风洞实验数据及对接网格计算结果进行了比较与分析,验证了拼接网格技术的高效性与可靠性。同时分析研究了绕三维增升装置的流动及其周围复杂的粘性流动现象。     

Some features of the buckling of stringer shells

A technique for stability analysis of stringer shells is proposed. It is used to analyze the minimum critical stresses. The dependence of the dimensionless parameters σ_cr/σ_cl on the number of stringers is plotted. The linear and nonlinear theories of ribbed shells are used to examine the features of how stringer shells lose stability. It is shown that the minimum critical stresses determined using the theory of ribbed shells and a structurally orthotropic model are close within the range of stiffness parameters considered __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 2, pp. 59–64, February 2006.     

Features of the turbulent flow around symmetric elongated bluff bodies

The flow fields around three elongated bluff bodies with the same chord-to-thickness ratios but distinct leading and trailing edge details were measured at a Reynolds number of 3×10⁴. These models each represent a case where: leading edge shedding dominates, trailing edge shedding dominates and a case where there is a balance between the two. The results show that the vortex street parameters vary between the models, and in particular, the shedding frequencies are significantly altered by the geometry. However, contrary to the current understanding for shorter bluff bodies, the scale of the recirculation region is found to be similar for each model, even though the shedding frequency changes within the range from 0.15 to 0.24. Also, the base pressure does not follow trends with shedding frequency expected from shorter bluff bodies. A force balance of the recirculation region shows that the near wake of each body is significantly affected by the Reynolds shear stress distribution and the resultant force due to the pressure field in the mean recirculation region. These differences infer that the distinct vortex formation characteristics depend on the state of the trailing edge shear layers. The boundary layers at the trailing edge have been quantified, as have the leading edge separation bubbles, and the marked differences in the wake details are shown to depend on the leading edge separation.     

风口区空气流动的有限元模拟与实验研究

通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.     

Effects of body shape and viscoelasticity on the slow flow around an obstacle

A study was made of the motion of axisymmetrical objects in viscous and viscoelastic fluids within a cylindrical tank with the assumption of negligible inertial effects. A numerical treatment based on the Stokes equations of motion and an optimization technique enabled the details of the velocity and rate-of-deformation fields for a Newtonian fluid to be predicted. The influence of the shape of various bodies, some even with concave surfaces, was examined. The corresponding experiments were carried out with viscous and viscoelastic fluids using a visualization technique. A correlation between the main flow characteristics and the rheological behaviour of the fluids was established.     

Factors affecting the stability of viscous axial flow in annuli with a rotating inner cylinder

Hot-wire measurements are presented of the onset of instability in developed axial flow and in both developing and developed tangential flow caused by inner cylinder rotation in concentric annuli of radius ratio N of 0.909, 0.809 and 0.565 for axial-flow Reynolds numbers (Re) between 86 and 2000. Within assessed uncertainty intervals, the consistency of marginal stability measurements, at four azimuthal locations 90° apart, indicates insensitivity to small variations in gap width; the measurements also confirm the destabilisation of nearly-developed and developed tangential flow identified by Takeuchi and Jankowski¹ with the occurrence at increasing Re of three-dimensional initial disturbances of spiral-vortex form. Comparison with earlier measurements suggests that in particular annuli, destabilisation may be delayed to higher Re by high values of certain geometrical factors, including radius ratio and the resultant end-effects parameter. Stability may also be restored or improved at high Re by reversion to developing tangential flow in which the initial instability is not of spiral-vortex form and where, for given N, the critical Taylor number appears uniquely related to the dimensionless axial co-ordinate. Stability is then generally greatest at low N.     

AN EXTENDED k-ε MODEL FOR NUMERICAL SIMULATION OF WIND FLOW AROUND BUILDINGS

It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold’s normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.