用扰动涡方法模拟叶轮机内动静叶的相互作用

发展了一种研究叶轮机内动、静叶间的相互作用的新方法———扰动涡方法，它利用全三维的定常解为基础解，并由此给出非定常扰动场的初始解．为计算叶片对扰动场的响应过程，采用拉格朗日方法追踪扰动涡团的对流流动过程，用确定性涡方法来描述流体的粘性扩散过程．发展了代数湍流模型（Ｂａｌｄｗｉｎ Ｌｏｍａｘ湍流模型）在尾迹中的应用方法，克服了其它数值方法中无法准确捕捉尾迹中心线的运动轨迹，以及计算出的边界层外的湍流涡粘性系数偏大的缺陷．利用该方法计算轴流叶轮机内由于动、静叶间的相互作用而引起的非定常流动过程，与实验的对照表明，模拟结果与实验数据吻合得相当好，从而说明本文发展的方法是可信的，为更直观地描述尾迹等非定常因素的流动及叶轮机内的掺混问题提供依据．     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

水力机械非定常流动的三维涡方法计算

利用拉格朗日型涡方法计算了混流式水轮机转轮内部的非定常流动. 采用高精度的边界元法计算域内势流速度以满足固体边界不可穿透条件, 通过在固体边界上导入新生涡满足边界的无滑移条件. 利用拉格朗日法计算涡元的运动, 而离散涡元之间的诱导速度则通过毕奥-萨伐尔公式结合快速多级子算法和自适应树结构算法获得. 通过对某混流式转轮水轮机在设计工况及非设计工况下的数值模拟, 以及模拟结果与实验测量的对比, 验证了该拉格朗日涡方法在计算转轮内部高雷诺数非定常流动中的有效性.     

三个共轴涡环运动混沌特性的研究

本言语利用离散涡环方法对三个共轴涡环运动进行了数值模拟，在此基础上，采用非线性动力学中混沌特征量的计算方法，得到了涡环呈规则运动和混沌运动的结论。结果表明，三个共轴涡环在运动时，总是会产主规划与混沌的现象。而这些现象的产生取决于三个涡环的初一文析得出对于利用共轴涡环模拟轴对称射流场涡结构的演变具有重要的参考价值。     

气固两相混合层二维涡配对的数值研究

采用双向耦合模型中的速度耦合模型,数值模拟了气固两相混合层中涡的配对、合并过程,文中采用拟谱方法对流场进行了直接数值模拟,用Lagrange模型跟踪固粒,结果发现,在双向耦合过程中Stokes数仍然是重要的参数,但不是唯一影响流场的参数,流场的发展还与固粒的体积浓度、固粒的相对密度以及固粒大小等因素有关。固粒抑制流场的发展,阻碍涡的配对、合并,加快了涡量的扩散;小St数的固粒仍然跟随流体运动,大St数的固粒趋向于低涡量区的趋势减弱。     

气固两相边界层中固粒与拟序结构相互作用的研究

本文对气固两相边界层中固体颗粒与拟序结构的相互作用进行了研究，建立了基于速度修正的双向耦合模型，提出了计算固粒对流场反作用的新方法以及大大减少计算量的快速涡方法，并据此得到了边界气固两相之间的相互作用结果。     

不可压缩湍流非对称脱体涡流场数值分析

采用拟压缩性方法，Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式数值求解三维定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。对Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，采用Ｄｅｇａｎｉ－ｓｃｈｉｆｆ修正。计算绕尖头正切拱型旋成柱体的大迎角大雷诺数脱体涡流场，计算结果中的非对称脱体涡与实验相符．     

全附体潜艇粘性流场的RANS模拟及场量和涡量的校验分析

采用相同拓扑结构和相近网格质量的4套网格和5种湍流模型,对全附体Suboff潜艇粘性流场进行RANS模拟,分析了网格密度、节点空间分布规律和湍流模型对计算精度的影响,详细校验了其力积分量、速度场量和涡量特征。结果表明:网格密度最大的G4网格(140万)计算精度最高,总阻力较实验值误差为0.723%,其采用SST湍流模型时最优。计算得到的压力系数和剪切应力系数分布均与实验值吻合很好;桨盘面速度等值线分布计算精度与文献相当,轴向相对速度0.9以上的计算半径稍大于实验值,其余半径与实验吻合较好;桨盘面上0.25倍半径处速度分量沿周向分布计算精度较文献高,轴向分量与实验值吻合较好,径向分量峰值稍小于实验值,但峰值所处周向位置与实验值一致。成功捕捉到了附体端面绕流诱导对旋涡、附体叶根截面下游处项链形涡对、尾翼端面尾缘上方附着涡蹄、附体马蹄涡系、尾翼截面通道流体挤压作用诱导涡以及桨盘面涡量汇集的潜艇涡量场特征,且围壳端面绕流诱导对旋涡沿流动方向持续稳定,不影响桨盘面涡量场,均与文献中由大涡模拟模拟得到的定性结论一致。研究表明,在网格密度较大、节点分布合理、网格质量较高、湍流模型选取适当和壁面函数使用有效的条件下,RANS模拟潜艇粘性流场的场量和涡量特征同样具有很高的计算精度,能够在工程应用中有力支撑新型艇型设计与性能分析。     

随动曲线坐标能量法模拟平面塑性压缩

本文给出平面压缩时随动曲线坐标系下运动可能速度场和相应的力学张量表达式.利用电算法优化可变参数使总能量泛函达最小值,从而获得流动速度场、位移场、应变速率和应力场的数值解.描述出不同宽高比时坯料产生单鼓形和双鼓形等现象,计算结果与实验吻合较好.     

流固耦合分析的一种改进CBS有限元算法

针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.     

隔水管弯曲对钻柱振动影响的计算与分析

海洋油气资源钻探中隔水管的弯曲对钻柱振动以及钻进特性有特别的影响。为得到隔水管弯曲对钻柱振动的影响规律,对南海已钻深水井使用非线性有限元软件建立全井钻井数值计算模型,研究获得了不同垂深时隔水管弯曲对钻柱振动特性的影响规律。研究表明:隔水管弯曲会加剧钻柱的振动,钻柱振动加剧会导致钻井能耗上升、钻头切削能力下降并且会加快钻柱疲劳;当隔水管的弯曲达到某临界值,钻柱与隔水管间的接触力会陡增;井口的钩载越大,隔水管弯曲带来的井口钩载波动量越大;井越深,隔水管弯曲对全井钻柱最大弯矩和钻头切削能力的影响越小。     

进气道入口形状对冲压发动机性能影响数值研究

进气道是冲压发动机设计中非常重要的环节,其入口设计会直接影响冲压发动机的总压恢复系数σ和流量系数φ等重要性能。本文采用二阶迎风隐式TVD格式,内外流分区耦合求解可压缩N-S方程,数值模拟了掺混段出口压力与来流压力比为P/P∞=4.2超临界状态下,二元方形截面“X”型布局进气道、弹体和掺混段一体化通气模型复杂流场。计算比较了两种不同二元进气道入口形状流场和冲压发动机总压恢复系数σ、流量系数φ的结果,对其产生的原因进行了分析。     

On the existence of the derivative of the volume average

A general theorem on the derivative of the volume average is formulated and proved. Conditions for the existence of the derivative are presented and discussed. This is done in order to give a better base to the theory of spatial averaging.Latin Letters E ₃ three-dimensional vector space over the field of real numbers - K, K(x) averaging domain - G, G _w, G_s open sets in E ₃; components of the two-phase system - C ¹(G) the set of functions 1-times continuously differentiable in G - W1/2(G) Sobolev space - V volume of the domain K - f function defined in G, G _w - K _infi ^sup* (x), K _infi ^sup– (x) special parts of K(x) Greek Letters boundary of G, G _w, G_s; w-s interface - _ij Kronecker delta - v unit outward normal of G, G _w - _j j-dimensional Lebesgue measure Other M closure of a set M in the metric space E ₃ - f phase average of f for the w-phase - (u, v) scalar product of u, v in E ₃ - one-sided derivatives     

表面织构对刀具切削性能及前刀面摩擦特性的影响

为了研究刀具表面织构对切削性能的影响,通过控制前刀面上微凹坑的直径、深度和面积占有率等参数制备了两组表面高度算术平均值Sa相同的织构刀具.使用Talysurf CCI Lite非接触式三维光学轮廓仪对织构表面进行测量,采用ISO25178参数对测量表面进行表征;通过SL200KS接触角测量仪测量织构刀具前刀面湿润性的变化;随后在CA6140车床上进行45钢切削试验,研究不同织构对切削力、前刀面摩擦特性、工件表面粗糙度和切屑形变的影响.结果表明:合理的表面织构能够改善刀具的切削性能;切削液在前刀面铺展较快的织构刀具获得了较好的切削性能;凹坑直径、深度和面积占有率的协同作用对改善刀具切削性能至关重要,为了协调这三个织构参数,进一步研究三维表征参数Sku、Ssk、Sv、Vvv、Vvc、Sdv、Sda、Spd与刀具切削性能的关联性,为刀具表面设计提供了基础.     

On the derivation for the gradients of the principal invariants

In this notes, it is shown that the method carried out by Truesdell and Noll [1] for the gradients of the principal invariants is available to be used for arbitrary second-order A. Finally, some trace identities for the derivatives of invariants are given.     

On the study of the initiation of the micro crack on the smooth surface of polycrystalline

Crystallographic plasticity was applied to study the initiation of micro cracks on the smooth surface of polycrystalline under uniform applied stress. Even under the uniform external stress, due to the different crystallographic orientations of the grains in the polycrystalline, there is un-uniform stress distribution and the deformation is also not uniform. Under the fatigue loading, the roughness increases with the number of fatigue, and deformation will localize in some places, where micro cracks form.     

地球与月球关系的研究

基于鹦鹉螺外壳生长线条数的变化是月亮绕地球运动周期随年代变化反映的假设, 采用外推 法得出两个结论: (1) 月球最初被地球俘获绕地球运动时距地球的距离为目前月地距离的0.2 倍左右; (2) 月球最初被地球俘获绕地球运动时的时间距今为7亿年左右. 这 为进一步研究地球与月球的关系提供了有益的启示.     

On the criterion for the absolute stability of the control system

In this paper, firstly we give the criterion for the absolute stability of the second canonical form for the control system, including the equation of the longitudinal motions of a plane as a particular example. The corresponding result in [8], [9] is a particular example given in this paper. Secondly, we give the criteria for the absolute stability of the first canonical form in the usual case and in the critical case. Finally, we give some criteria for the absolute stability of the general form for the direct control system.All the results in this paper merely depend upon the relations between the parameters of the system itself to give an explicit algebraic discriminant.     

捷联式惯性导航系统惯性元件的设置与可靠性

本文对惯性元件的最佳冗余设置及其在飞行器上的布局与可靠性之间的关系进行了研究，建立了相应的可靠性数学模型，探讨了３／６（Ｇ）冗余惯性元件出现第三个故障时故障的检测与识别的方法。所得结果可供惯性导航系统可靠性分析及工程上应用。