绕流直接数值模拟误差分析

对于无边界绕流问题的计算流体力学模拟通常是将物体置于“足够大”的槽道中，而通过不断改变槽道尺寸以及离散网格密度，后验对比方式来检查模拟误差。本文结合多种经典流场理论，提出一种简单的先验误差估计方法来确定槽道尺寸以及相应的网格分布。在此方法中，对于槽道尺寸的确定基于线性叠加原理（即在极小雷诺数下采用Stokes理论解叠加，而在其他雷诺数条件下采用势流理论解叠加），来估计槽道尺寸对绕流结果的影响。而对网格尺寸与分布，则是使用多项式逼近中的基本误差分析工具，应用到速度边界层，远场势流，以及Rankine涡等简单流动，从而确定整个绕流问题中的离散误差。为了验证前面的理论分析结果，本文模拟了相当大雷诺数范围内的二维翼型以及三维圆球绕流，所得数值结果非常好地验证了理论分析。结果表明，对于Stokes流动问题，槽道尺寸需要大约100倍于物体特征尺寸来保证其结果与无边界绕流相差不超过1%；而在雷诺数超过大约100时，槽道尺寸只需10倍（二维绕流）或者5倍（三维绕流）于物体特征尺寸来达到同等精度。在此先验误差估计方法可应用于一般化的绕流问题。     

具有弹性支承的圆管在内外部流激发下动力特性及实验研究

本文建立了具有弹性支承的圆管在内外部流激发下的力学模型.推导了内部流与静止外部流作用下圆管的耦联振动方程.提出了确定弹性系数的方法.采用振型叠加法分析圆管动力特性问题.对内部流与静止外部流情况下圆管固有频率进行了计算和测量,计算值与实验值吻合较好.另外,对内外流同时激发下圆管的固有频率进行了测量,得到若干对工程实际有用的结论.     

机械载荷作用下单边裂纹载流薄板的应力场

采用坐标变换的方式,将单边裂纹载流薄板通电瞬间由温度产生的应力场表达式中的各应力分量分离,并用极坐标进行表示．给出了Ⅰ型穿透裂纹尖端附近的应力场的表达式．最后将温度产生的应力场与单向拉伸载荷作用产生的应力场相叠加,推导出用极坐标表示的机械载荷作用下单边裂纹载流薄板的应力场的表达式,并给出算例．     

欧拉-伯努利梁运动场的正交性及其能量传导特性分析

从无阻尼欧拉—伯努利梁振动方程解析解出发, 推导了有限长梁的关于谱系数的时间—空间平均能量和功率流表达式. 在此基础上, 从泛函分析观点, 探讨了弯曲运动场: 衰减振动、行波模式分解关于能量、功率泛函的正交性. 结果表明: 弯曲衰减振动模式和行波模式关于功率流、机械能时间—空间平均是相互独立的, 即关于场能和场功率互不干涉, 满足叠加原理; 衰减振动场导能与行波场导能的重要区别在于功率流关于右、左衰振动模式分解不满足叠加原理, 即弯曲衰减振动场间的相互"干涉"是使其具有能量传导能力的内在原因. 通过右端集中阻尼器支撑的梁的稳态功率流仿真分析计算, 表明低频区振动导能不可忽略, 同时, 衰减振动场和行波场间存在一定的能量交换现象, 但随着频率升高, 振动场传导能量不断下降, 同时能量传导效率也不断下降.      

蜿蜒边界下平面流的线性流动稳定性

为便于数值分析,蜿蜒河流水动力和演变模型中一般隐性假设二次时均流-二次涡的关系与明渠流时均流-明渠湍流的关系相同,但由于高雷诺数下的DNS算力限制和实验尺度限制,这种隐含假设是否成立目前尚无相关湍流研究来支撑.文章试图通过分析明渠湍流和二次湍流发展初期的研究,侧面揭示其湍流结构的异同.通过对曲线正交坐标系下的平面二维NS方程使用双参数摄动的方法,建立了一种求解蜿蜒边界弱非线性层流的摄动解法,并推导得出一个适用于蜿蜒边界的EOS方程以及其特征值问题的解法.蜿蜒边界下弱非线性层流解为一系列蜿蜒谐波分量的叠加,其中线性部分使得两壁产生流速差,非线性部分随着雷诺数增大呈指数增长.水流的扰动增长率特征谱的第一模态与直道流相似,由3条曲线、4个波段合成,但其长波段和短波段的扰动流场与直道流不同,所有短波段的扰动流速近似于KH涡.蜿蜒边界对内部水流扰动有一定的选择性.偏角幅值越大扰动增长越快;蜿蜒波数的影响则为先增后减,有一个使扰动增长最快的蜿蜒波数.扰动流场由一个典型的TS波和一对波包形式的二次涡叠加而成,波包只有纵向流速分量,包络线由蜿蜒波数控制,波包内是与直道扰动波参数相同的TS波.     

可变形儒可夫斯基翼型非定常气动力的研究

对于翼面变形法向运动速度远小于来流速度的儒可夫斯基机翼,将解析解和离散涡方法相结合计算变形机翼的流场及非定常气动力,较详细地分析了变形机翼升力系数的准定常计算方法的误差来源,并给出修正方法.计算结果显示脱落涡尾迹对升力系数和机翼绕流环量的影响很小,变形机翼升力系数准定常计算方法的误差丰要来源于流体非定常运动引起的虚拟质量力,该非定常附加升力仅与当前时刻飞行姿态及翼犁形状和变形速率有关,与具体的变形历史过程无关,变形机翼的升力近似等于准定常计算结果叠加上相应的虚拟质量力.     

变截面悬臂梁的刚度叠加法

利用梁的弯曲刚度与位移成反比的关系, 可较为方便地计算变截面悬臂梁端截面的位移．该 方法是区别于载荷叠加法和逐段变形效应叠加法的叠加方法．     

缓坡方程的推广

为了描述水波和强烈的环境流在非平整海底上的相互作用,运用无旋运动的Lagrangian变分原理,对经典的Berkhoff缓坡方程进行了改进。假定水流沿水深方向基本上保持均匀性,这正如潮流运动的特征。海底地形由慢变、快变两个分量叠加构成;慢变分量满足缓坡逼近假定,快变分量的波长与表面波波长为同一量级,但其振幅小于表面波的振幅。在以上假定条件下,得到了适用于非平整海底的推广型浅水方程和应用性更加广泛的波－流－非平整海底相互作用的一般缓坡方程,并且从理论上证明一般缓坡方程包含了以下3种著名的缓坡型方程：经典的Berkhoff缓坡方程;波－流相互作用的Kirby缓坡方程、Dingemans关于沙纹海底的缓坡方程。最后,通过与Bragg反射实验数据的比较,表明该模型可以准确地反映快变海底的典型地貌特征。     

多边形板高频振动的能量虚源法研究

本研究的目的是将室内声学中的虚源法引入高频横向点激励作用下多边形板的能量辐射传递模型中,以预测结构弯曲振动的能量响应,其核心是如何构建多边形板的虚源云场。在多边形板的能量辐射传递模型框架中,以实源的镜像虚源替代结构的边界,以此来表示边界处能量的连续反射。基于直接场的能量密度和功率流强度核函数,结构内部任意一点的能量可由实源产生的能量和无限多的镜像源产生的能量叠加得到。对直角三角形板和方形板求解,获得了它们的能量密度与功率流强度分布特征。数值算例结果与能量辐射传递法以及解析解对比,验证了模型的准确性。最后讨论了阻尼、频率和虚源的阶数对计算结果收敛性的影响。     

水流冲击管道内滞留气团现象的VOF模型仿真分析

针对复杂管道系统内水流冲击滞留气团现象,采用VOF模型(Volume of Fluid Model)进行了数值模拟计算,并与一维模型进行了比较计算分析,结果表明:系统的最大压力并不总是气团的最大压力,有可能还会叠加水体对管壁的撞击而形成的突然升高压力.与实验实测结果的比较分析表明:采用VOF模型,能够较精细地仿真水流冲击滞留气团现象的气团形态、流场结构以及压力分布等的变化过程,其压力数值计算结果与实验实测基本吻合,其计算误差明显小于现有一维模型的计算误差,是深入研究该复杂瞬变流现象的有效方法.     

Stagnation-point flow of a viscoelastic fluid towards a stretching surface

An analysis is made of the steady two-dimensional stagnation-point flow of an incompressible viscoelastic fluid over a flat deformable surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. It is shown that for a viscoelastic fluid of short memory (obeying Walters’ B′ model), a boundary layer is formed when the stretching velocity of the surface is less than the inviscid free-stream velocity and velocity at a point increases with increase in the elasticity of the fluid. On the other hand, an inverted boundary layer is formed when the surface stretching velocity exceeds the velocity of the free stream and the velocity decreases with increase in the elasticity of the fluid. A novel result of the analysis is that the flow near the stretching surface is that corresponding to an inviscid stagnation-point flow when the surface stretching velocity is equal to the velocity of the free stream. Temperature distribution in the boundary layer is found when the surface is held at constant temperature and surface heat flux is determined. It is found that temperature at a point decreases with increase in the elasticity of the fluid.     

Deformation of an ellipsoidal ferrogel sample in a uniform magnetic field

The elongation of a ferroelastic material sample (whose initial shape is a sphere or an ellipsoid of revolution) under the action of an external magnetic field is studied in an in approximation of small strains. For a sphere, there is a classical estimate obtained under the assumption that elongating in the direction of the field, it becomes a spheroid and the stress and strain fields remain uniform. In the present calculation, it is assumed that the body is an ellipsoid (a sphere in a particular case) only in the absence of an external field; the shape of the sample in the presence of a field is not specified in advance but is found from the condition of balance of surface forces (elastic and magnetic). For the spherical case, the problem is solved exactly: it is shown, that the contour of the deformed body is described by a third-order algebraic equation. The case where the initial configuration is an ellipsoid of revolution is studied numerically. It is shown that in all versions, the refined solution leads to an appreciable increase in the elongation of the sample compared to the classical estimate. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 153–164, May–June, 2005.     

Nonlinear and chaotic analysis of a financial complex system

In this paper,determination of the characteristics of futures market in China is presented by the method of the phase-randomized surrogate data.There is a significant difference in the obtained critical values when this method is used for random timeseries and for nonlinear chaotic timeseries.The singular value decomposition is used to reduce noise in the chaotic timeseries.The phase space of chaotic timeseries is decomposed into range space and null noise space.The original chaotic timeseries in range space is restructured.The method of strong disturbance based on the improved general constrained randomized method is further adopted to re-deternination.With the calculated results,an analysis on the trend of futures market of commodity is made in this paper.The results indicate that China＇s futures market of commodity is a complicated nonlinear system with obvious nonlinear chaotic characteristic.     

Momentum and heat transfer in the magnetohydrodynamic stagnation-point flow of a viscoelastic fluid toward a stretching surface

An analysis is made of the steady two-dimensional stagnation-point flow of an incompressible viscoelastic fluid over a flat deformable surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. It is shown that for a viscoelastic conducting fluid of short memory (obeying Walters’ Bʹ model), a boundary layer is formed when the stretching velocity of the surface is less than the inviscid free-stream velocity and velocity at a point increases with increase in the Hartmann number. On the other hand an inverted boundary layer is formed when the surface stretching velocity exceeds the velocity of the free stream and the velocity decreases with increase in the Hartmann number. A novel result of the analysis is that the flow near the stretching surface is that corresponding to an inviscid stagnation-point flow when the surface stretching velocity is equal to the velocity of the free stream. Temperature distribution in the boundary layer is found when the surface is held at constant temperature and surface heat flux is determined. It is found that in the absence of viscous and Ohmic dissipation and strain energy in the flow, temperature at a point decreases with increase in the Hartmann number.     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

汶川地震触发陕西省境内地质灾害灾险情特征

在实地调查、汶川地震震后次生地质灾害隐患排查成果的基础上,结合县市地质灾害调查数据,从灾情和险情两个方面,对汶川地震触发地质灾害灾险情的致灾灾险情、灾种、类型等特征进行了较系统的分析对比。进一步比较了陕西在汶川地震触发地质灾害受灾省份中的情况,并对地质灾害趋势做了初步预测。研究表明,灾险情集中分布在龙门山构造带向陕西境内的延伸部分,险情比灾情分布范围广。灾险情的灾种都以崩塌为主,其次为滑坡,地震对地质体的震动触发作用明显。机关学校是汶川地震触发陕西境内地质灾害的主要受灾对象。陕西受汶川地震触发地质灾害的灾险情轻于四川、甘肃,三省地质灾害震前均较为严重,地震触发地质灾害的影响具有长期性,未来陕西境内地质灾害发生的频率可能增大。     

Hydraulic analysis of free-surface flows with solidification

Summary The paper is concerned with a one-dimensional analysis of plane open-channel flow with continuous solidification. The process is of relevance for recent developments in the casting of steel and other metals. The bottom of the channel consists of a rotating casting roll and a horizontal cooling table, where the solidified material is withdrawn with given velocity. The study is restricted to the region downstream of the top of the casting roll. Surface tension is neglected. In the main part of the analysis inviscid fluid flow is considered since the Reynolds number is very large in the applications. It is found that the steady-state solutions are nonunique in a certain parameter range. In addition to a continuous solution, there are two solutions including hydraulic jumps, with one hydraulic jump being located on the casting roll, the other one on the cooling table. Regarding the stability of the non unique solutions, the evolution of disturbances is investigated numerically as an initial-value problem. It is concluded that the hydraulic jump on the cooling table is unstable, while the other discontinuous solution as well as the continuous solution are stable for sufficiently small disturbances. Which stable solution is attained in the steady state, depends on the history of the process. Friction at the liquid/solid interface is taken into account in the last part of the analysis. A constant friction coefficient is assumed. It is found that the history of the process determines the steady-state solution if, and only if, the friction coefficient is sufficiently small. For larger values of the friction coefficient, the steady-state solution is unique and independent of the history of the transient process. Furthermore, for sufficiently large friction coefficients, stable hydraulic jumps are found, in contrast to the inviscid case, also on the cooling table. Received 19 March 1999; accepted for publication 3 May 1999     

Performance of a flapping foil flow energy harvester in shear flows

Through numerical simulations, we investigate the energy harvesting performance of a heaving/pitching foil in shear flow. With two-dimensional Navier–Stokes simulations, we examined the energy harvesting efficiencies of such a system in linear shear flows and compared the results with those in uniform flows. It is found that in low shear rates, the performance of the system in linear shear flow is slightly higher than that in uniform flow, whereas the energy harvesting efficiency is greatly diminished if the shear rate is sufficiently high (this effect is more pronounced in higher frequencies). This is attributed to the effects of linear shear on the vorticity generation and the synchronization between fluid forcing and foil motion – when a strong shear flow is introduced the lift force induced by the leading edge vortex that is in phase with the heaving motion of the foil is diminished. Furthermore, by studying the instability of the wake behind the foil, we confirm that the optimal performance of the foil in linear shear flows is associated with the same physical mechanism that controls the performance of the foil in uniform flows, i.e. the excitation of the most unstable modes in the wake when the oscillation frequency of the foil is close to the frequencies of these modes.     

新疆下坂地水库坝基防渗墙试验研究

新疆下坂地水利枢纽工程地处高寒、高海拔的山区,坝基覆盖层厚度达147.95m,透水性强,易坍塌和漏浆,块石坚硬,坝基防渗处理难度罕见。针对这一工程问题进行了坝基防渗墙试验研究,结果表明,以大深度混凝土防渗墙为主墙,下接帷幕灌浆的墙幕结合方案。本文针对新疆下坂地深厚覆盖层坝基防渗墙进行了试验研究,找出了防渗墙有效的施工深度及合理的施工方法,从而为坝基防渗墙设计施工提供了可靠的依据。     

岩溶塌陷的类型、成因机制及防治途径──兼论武汉市岩溶塌陷勘查工作要点

按岩溶塌陷的类型, 对其成因机制作了分析。在此基础上, 提出了武汉市岩溶塌陷勘查与防治工作要点。