并列双方柱绕流的大涡模拟及频谱分析

应用二阶全展开ETG有限元离散格式与大涡模拟相结合的方法,对间距比为1.5情况下的并列双方柱绕流进行了数值模拟．由王小华、何钟怡提出的二阶全展开ETG有限元方法通过对N-S方程中的时变项进行Taylor展开,从而把时间导数用空间导数来代替,其作用相当于引入了人工粘性．计算得到了不同时刻的流线图,给出了两方柱的阻力系数、升力系数以及两对称点上流向速度随时间的变化历程,并采用谱分析的方法研究了对称边界条件下并列双方柱绕流的频谱对称性问题．为了消除初始条件的影响,在所取样本中去除了计算中初始段的数据,分别分析了阻力系数、升力系数以及两对称点上流向速度的频谱．结果表明:对称边界条件下,双方柱绕流运动参量的时域过程虽然是不对称的,但频域过程却是对称的．     

功能梯度压电圆柱体的旋转分析及其应用——内压和热荷载作用下的力学传感器

对旋转的功能梯度压电(FGP)圆柱体,进行了精确的热弹性分析．圆柱体同时承受电、热和力学荷载的作用,这是一个旋转传感器或者调节器的简化模型．应用能量法得到系统的控制微分方程．为了正确地评估能量函数,引进了一个被称为“附加能量”的新项．在两种边界条件下:自由旋转的圆柱体和受内压作用的圆柱体,求解所得到的控制方程．研究角速度对各个物理量沿半径分布的影响．所研究的结构也可以认为是,在压力及热荷载作用下,测量圆柱体角速度的一个传感器．结果表明,电势与角速度成正比例关系．     

带有球凹/球凸的旋转通道内流动结构和强化换热的数值研究

对带有凹坑和凸包的内流通道在不同旋转数下的对流换热特性进行了数值分析,探讨了Coriolis力对通道中流场和换热特征的影响．研究发现,随着旋转数增加,通道前缘呈现出较弱的流动冲击,但存在较大的尾迹和延迟的流动再附着,后缘凹坑内部有一较小旋涡和较强射流使得后缘传热得到强化,最高可达60%．总体Nusselt数随着旋转数的增加先减小而后增大.     

无限圆柱体旋转运动时的热应力

研究旋转对确定边界条件下无限圆柱体的影响.当热荷载沿径向作用时,给出了旋转圆柱体中热应力、位移和温度的分析过程.当无限弹性圆柱体部分弯曲界面有常温作用,而其余界面维持零温度时,讨论其热动应力的分布.圆柱体表面绝缘材料熔化时出现这种情况.得到了应力分量、位移分量和温度的解和数值结果.提出的半解析法所得到的结果,与早期采用方法所得到的结果比较,发现两者显示出很好的一致性.     

圆管纤维悬浮湍流场中粒子运动的研究

利用从细长体理论出发得到的三维分段积分法和湍流简化方法模拟了大量纤维粒子在圆管湍流内的运动．统计了不同Re数下计算区域内的纤维的取向分布,计算结果与实验结果基本吻合,结果表明湍流的脉动速度导致纤维取向趋于无序,且随着Re数的增加,纤维取向的分布越来越趋于均匀．其后又考虑了纤维速度和角速度的脉动,二者都充分体现了流体速度脉动的影响,且纤维速度的脉动在流向上的强度大于横向,而其角速度的脉动在流向上的强度小于横向．最后统计了纤维在管道截面上的位置分布,说明Re数的增加加速了纤维在管道截面上的位置扩散．     

非Darcy多孔介质中的非线性对流

利用温度-浓度-密度关系,研究非Darcy多孔介质中的自由对流问题.对于不同的惯性参数、传递参数、Rayleigh数、Lewis数、Soret数和Dufour数,分析了非线性温度参数和浓度参数对非线性对流的影响.浮力对对流起着辅助的附加作用,当惯性作用不计时,切向速度随着非线性温度和浓度的增加而急剧地增加.然而,当惯性效应不为0时,非线性温度和浓度对切向速度的影响是有限的.对两个传递参数、惯性影响参数以及控制非线性温度和浓度的其他参数,取不同的数值时,浓度分布有点儿变化,并在不同的范围内传播.随着非线性温度和浓度的增加,传热/传质在很大的范围内变化,这取决于是Dacry多孔介质,还是非Darcy多孔介质.当所有的影响(惯性的影响、两个传递系数的影响、Soret和Dufour的影响)同时为0/不为0,在非线性温度/浓度参数以及浮力的共同作用下,分析了传热/传质的变化.发现在Darcy多孔介质中,温度和浓度以及它们的交叉扩散,对传热/传质的影响,要比非Darcy多孔介质要大.发现了浮力的负面作用,随着非线性温度系数的增加,传热/传质率是提高的,而随着非线性浓度系数的增加,传热/传质率是下降的.     

粘性导电流体在磁场作用下流过多孔通道时传热传质振荡流的数值解——用于病理状态动脉中血液的流动

当血管内壁出现多孔性结构时,流过多孔性血管的血液将作不稳定的MHD流动.研究血液在其中的传热传质问题,考虑了与时间相关的渗透率和振荡引起的吸入速度,并数值地求解该问题.对分析中出现的参数取不同数值时,图形给出了速度、温度、浓度场,以及表面摩擦因数、Nusselt数和Sherwood数的计算结果.研究表明,血液流动受磁场和Grashof数的影响明显.     

固壁近旁Stokes流中粘性液滴的运动和变形

发展了一种模拟固壁近旁轴对称Stokes流中粘性液滴的运动和变形及直接计算固壁上应力的边界积分方法．用此方法对不同的液滴-固壁初始相对间距、粘度比、表面张力和浮力联合参数以及环境流动参数情况进行了数值实验．数值结果显示,由于环境流动和浮力的作用,随着时间的推进,液滴在轴向压缩,在径向拉伸．当环境流动的作用弱于浮力作用时,随着时间的推移,液滴上升并向上弯,固壁上由液滴运动所引起的应力不断减小．当环境流动的作用强于浮力作用时,随着时间的推移,液滴变得越来越扁．在这种情形,当大初始间距时,壁面上的应力随液滴的演变而增大;当小初始间距时,由环境流动、浮力及壁面对流动的较强作用的联合影响,此应力随液滴的演变而减小．由于液滴运动所引起的壁面应力的有效作用仅限于对称轴附近的一个小范围内,且此范围随液滴与固壁的初始间距增大而增大．应力的大小随初始间距增大而大为减小．表面张力对液滴变形有阻止作用．液滴粘性会减小液滴的变形和位置迁移．     

在柱坐标系下研究由拉伸缸引起的Casson纳米流体流动与传热传质现象

首次利用柱坐标研究速度滑移和对流表面边界条件下,由拉伸缸引起的稳态层流Casson纳米流体流动、传热及传质现象.采用恰当的相似变换将偏微分控制方程转化为高阶非线性耦合常微分方程,并通过打靶法进行数值求解,图示并详细分析了不同物理参数对速度、温度及浓度分布的影响.结果显示,速度受滑移参数的影响较大,温度和浓度分别受Biot数和Lewis数的影响较大;随着Casson参数的增大,速度下降而温度和浓度都增加;温度随着Brown(布朗)运动参数或热泳参数的增加而上升;浓度随着Brown运动参数的增大而减小,随着热泳参数的增大而增大,当热泳参数较大时,浓度出现了"回流"现象.     

Stokes Flow in Cylindrical Containers With Rotating Ends北大核心CSCD

该文以端部旋转的圆柱形容器内的Stokes流为研究对象,根据流动的特点,将轴向坐标模拟为时间,则问题归结为Hamilton对偶方程的本征值和本征解问题.利用本征解空间的完备性和本征解之间的共轭辛正交关系,给出了问题解的展开形式,并建立了展开系数的数值求解方法.采用该方法研究了单端旋转、两端以相同或相反角速度旋转时不同外形比(容器的高度与半径之比)时圆柱形容器内流动速度和应力的分布情况,展示了不同边界条件下流场的一些特点.     

两并列方柱绕流相互干扰的数值研究

用数值方法研究了绕流两并列方柱的气动力及其相互干扰.使用了陈素琴等人提出的对流项用三阶逆风格式的带多重网格的改进的MAC方法,对两并列方柱的流场进行了模拟.结果表明:两并列方柱在不同的间距比时其干扰特征有很大的不同,小间距比时出现双稳态偏流,流动并不对称于双方柱间隙的中心轴线,而是偏向其中的一个方柱.     

Analysis of viscoelastic fluid flow with temperature dependent properties in plane Couette flow and thin annuli

The steady state flow in very thin annuli has been studied analytically for the case where the annular gap is much smaller than the radius of the inner cylinder and for the outer cylinder rotating at constant angular speed and the inner cylinder at rest. The cylinders were subjected to two different thermal boundary conditions. The exponential effect of temperature on the relaxation time and the viscosity coefficient was accounted into the governing differential equations using Nahme’s law. Effects of viscous dissipation as well as εDe² (viscoelastic index for SPTT constitutive equation) on the dimensionless velocity and temperature profiles have been investigated. Results show that while the properties of the fluid depend on temperature, the velocity and temperature profiles are different compared to those obtained with constant physical properties. The Nahme–Griffith number increases whereas εDe² as a viscoelastic index decreases when temperature dependent physical properties are considered. In addition, the results indicate that the viscous dissipation has a sensible effect on heat transfer and the Nusselt number decreases with an increase in the Nahme–Griffith number.     

MILU-CG方法和绕旋转圆柱流动的数值研究

本文采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),结合高阶隐式差分格式,改进了作者(1992)提出的基于区域分解、有限差分法与涡法杂交的数值方法(HDV).系统地研究了雷诺数Re=1000,200,旋转速度比α∈(0.5,3.25)范围内,绕旋转圆柱从突然起动到充分发展,长时间内尾流旋涡结构和阻力、升力系数的变化规律.计算所得流线与实验流场显示相比,完全吻合.首次揭示了临界状态时的旋涡结构特性,并指出最佳升阻比就在该状态附近得到.     

On the Influence of Thermal Convection in Classical Taylor‐Couette Flows

The presented work deals with the instabilities that occur in the flow of a viscous fluid between axisymmetric cylinders with a rotating inner and stationary outer cylinder. The results of a numerical study of convective flows are presented. The inner cylinder is rotating and heated from within, while the outer cylinder is stationary and cooled outside. Stationary horizontal endplates are used to seal the annulus, forming an enclosure. The working fluid is silikon oil M3. The flow of oil was rendered visible by injecting aluminium powder. By increasing the Reynolds number with angular velocity of the driving inner cylinder, the flow bifurcates into different types of instabilities. Investigation was aimed to find the values of critical Reynolds and Rayleigh numbers corresponding to the critical speeds and temperature differences at which these instabilities set in. The three‐dimensional problem was modelled numerically using software package FLUENT in which discretization is performed by means of finite volume techniques. Computational grid was created in preprocessor Gambit. Numerical experiments are conducted to determine the interdependence between the heat transfer mechanism and the structure of secondary flows     

Forced convective flow and heat transfer past an unconfined blunt headed cylinder at different angles of incidence

In the present paper, a numerical investigation has been carried out to study the forced convective flow and heat transfer characteristics past a blunt-headed cylinder in crossflow. Employing air as an operating fluid, calculations are carried out for a range of Reynolds number (Re) from 40 to 160. The angle of incidence is varied in the range of 0^∘ ≤ α ≤ 180 ^∘. The thermofluid features of flow and heat transport are analysed in detail for different angles of incidence. To analyse the aerodynamic characteristics, several parameters such as drag and lift coefficients, moment coefficient, Strouhal number, recirculation length, and local time-averaged vorticity flux have been calculated. Furthermore, a stability analysis has been undertaken by using the Stuart Landau equation to enumerate the critical Reynolds number at each angle of incidence. Heat transfer characteristics are studied by computing local and time-averaged values of Nusselt numbers. When compared to a rectangular cylinder, a blunt-headed cylinder exhibits an enhanced heat transfer rate. In the end, an entropy generation analysis has been carried out to study the effects of Re and angle of incidence on the efficiency of thermofluid transport characteristics.     

On the accuracy of the viscous form in simulations of incompressible flow problems

We present a numerical study of drag/lift and flux estimates using two forms of Navier‐Stokes equations (NSE) that are equivalent in the continuum formulation but not in the discrete finite element formulation. The two investigated forms of the NSE differ in the viscous term that is represented in one form by νΔ u with ν being the viscosity and 2ν?·?^S u in the other form where ?^S represents the deformation tensor. The study consists of numerical analysis of the two forms and computations of drag/lift, pressure drop on the cylinder problem and computations of flux for the Poiseuille flow. The main objective is to provide a clear comparison of the reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time for the two forms of NSEs. Our computational results of the reference values do not differ significantly between the two forms, but the differences are there. For the Poiseuille flow, the differences in the flux computations were much smaller, and this agreed with the computationally obtained results of the divergence of the velocity field. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 523–541, 2012     

Buoyancy-assisted flow of yield stress fluids past a cylinder: Effect of shape and channel confinement

The aiding-buoyancy mixed convection heat transfer in Bingham plastic fluids from an isothermal cylinder of elliptical and circular shape in a vertical adiabatic channel is numerically investigated. For a fixed shape of the elliptical cylinder E = 2 (ratio of major to minor axes), the effect of confinement is studied for three values of blockage ratio, B, defined as the ratio of the channel width to the circumference of the cylinder/π, as 6.5, 2.17 and 1.3. In order to delineate the role of cross-section of the cylinder, results are also presented here for a circular cylinder of the same heat transfer area as the elliptical cylinder. The results presented herein span the range of conditions as: Bingham number, 0 ≤ Bn ≤ 100, Reynolds number, 1 ≤ Re ≤ 40, and Prandtl number, 1 ≤ Pr ≤ 100 over the range of Richardson number Ri = 0 (pure forced convection) to Ri = 10. Extensive results on drag coefficient, local and surface averaged values of the Nusselt number and yield surfaces are presented herein to elucidate the combined effects of buoyancy, blockage ratio and fluid yield stress. The morphology of the yield surfaces shows that the unyielded plug regions formed upstream and downstream of the cylinder grow faster at low Reynolds numbers with the increasing yield stress effects under the weak buoyancy forces, i.e., small values of Grashof or Richardson number. The heat transfer enhancement is observed with the increasing channel-confinement due to the sharpening of the temperature gradients near the surface of the cylinder. The average Nusselt number shows a positive dependence on the Reynolds number, Prandtl number and Richardson number irrespective of the shape of the cylinder or the type of fluid. By employing the modified definitions of the dimensionless parameters (based on the two choices of the overall effective fluid velocity), predictive correlations have been established for estimating the value of the average Nusselt number in a new application.     

Analytic solutions for a rotating radial fin of rectangular and various convex parabolic profiles

The homotopy analysis method (HAM) is used to develop an analytical solution for the thermal performance of a radial fin of rectangular and various convex parabolic profiles mounted on a rotating shaft and losing heat by convection to its surroundings. The convection heat transfer coefficient is assumed to be a function of both the radial coordinate and the angular speed of the shaft. Results are presented for the temperature distribution, heat transfer rate, and the fin efficiency illustrating the effect of thickness profile, the ratio of outer to inner radius, and the angular speed of the shaft. Comparison of HAM results with the direct numerical solutions shows that the analytic results produced by HAM are highly accurate over a wide range of parameters that are likely to be encountered in practice.     

On low Reynolds number optimization of non-linear partial differential equations for convergent–divergent engulfment: A numerical result

The flow field in a convergent–divergent engulfment along with the installation of infinite cylinders as an obstacle results in non-linear partial differential equations and the scientific computation in this regard remains a challenging task. The present attempt is the numerical motivation in this direction to evaluate the flowing liquid stream in the convergent–divergent channel at a low Reynolds number. From the left wall, the liquid stream move with the parabolic profile and have interaction with the case-wise installation of infinite cylinders in the left vicinity of the convergent–divergent throat. The differential system is constructed for the flow field in the channel and hybrid meshed finite element method is utilized to report the numerical solution. A comparative study is enclosed for the hydrodynamic forces faced by obstructions in the left region of the convergent–divergent throat. The drag coefficient for a triangular cylinder acting as an obstruction is higher than that of a circular hitch. In comparison to both triangular and circular hitches, the square-shaped obstacle suffered the most drag force. Considering drag coefficient one can extend this work to obtain information for the real behavior of the vehicle toward air flow and may conclude findings toward reduction of fuel consumption.