展向喷嘴湍流横向射流大涡模拟及其统计分析

采用大涡模拟方法数值模拟了展向椭圆喷嘴的湍流横向射流,对其大尺度结构的时空演化和湍流脉动速度场的时间序列分析、频谱分析、PDF分析以及时、空截面上的统计平均特性进行分析.结果表明,在射流出口附近的下游核心区中速度脉动剧烈,显现出明显的湍流特征.除了三维涡环脱落、扭曲、变形、摆动所对应频率之外,还存在很宽的湍流基频,它与在喷嘴出口附近产生的三维涡环的时空演化过程密切相关.由于展向椭圆喷嘴的湍流横向射流中的三维涡环快速脱落和强相互作用导致射流尾迹中的强湍流脉动,展向椭圆喷嘴湍流横向射流的PDF空间演化特征结构复杂.在射流核心区的湍流偏应力变化平缓,其统计平均值分布接近左右对称.展向椭圆喷嘴的湍流横向射流脉动速度场具有极为复杂的统计行为,与流向椭圆喷嘴相比具有更好的掺混能力.     

应用保角变换法的亚音速矩形湍流射流理论研究

为求解亚音速矩形湍流射流的速度场结构,采用复变函数中的保角变换法将圆形喷口平面转化为矩形喷口平面;进而将已有的关于轴对称喷口射流流场的相关理论拓展到矩形喷口射流研究中,避免了在数值模拟中引入辅助方程、离散化、设置边界条件来求解三维N-S方程。通过该方法得到了矩形喷口宽高比AR分别为2、4、8时对射流结构的定量影响;求出了喷口尺寸为1×5mm2、喷口马赫数Ma=0.3时,喷口对称中心的宽、窄对称面内的速度云图。结果表明:宽高比AR越大,射流核心段长度越长,轴向衰减越慢;喷口面积为10cm2、20cm2时,加大宽高比可以缩短核心段长度且可加快轴向速度衰减,在射流宽、窄对称面都存在射流核心区与衰减区。首次提出速度分布的"马鞍面"主要集中在射流核心段内。通过与现有实验研究结果对比发现:核心段长度误差最小达到30%;轴向速度分布与实验结果吻合度也很高,误差在5%以内。因此,本文理论研究的正确性得到验证。本文所得结论可为进一步研究亚音速矩形射流结构与矩形喷口参数优化设计提供参考。     

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理, 完善等离子体激励器数学模型, 提升等离子体激励器扰动能力, 采用粒子图像测速技术, 在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究. 实验时, 将非对称布局激励器布置在平板模型上, 随后将带有激励器的模型放置在有机玻璃箱内, 从而避免环境气流对测试结果的影响. 基于激励器诱导流场, 分析了激励电压对诱导射流特性的影响, 揭示了较高电压下诱导射流近壁区的拟序结构, 获得了卷起涡、二次涡等拟序结构的演化发展过程, 计算了卷起涡脱落频率, 阐述了卷起涡与启动涡的区别, 初步探索了卷起涡的耗散机制. 结果表明: (1)层流射流不能完全概括等离子体诱导射流特性, 激励电压是影响射流特性的重要参数. 当电压较低时, 诱导射流为层流射流; 当电压较高时, 诱导射流的雷诺数提高, 射流剪切层不稳定, 层流射流逐渐发展为湍流射流. (2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构; 在固定电压下, 这些涡结构存在恒定的卷起频率. (3)当激励电压较高时, 流动不稳定使得卷起涡发生了拉伸、变形, 引起了流场湍动能增大, 从而加速了卷起涡的耗散. 研究结果为全面认识激励器射流特性, 进一步挖掘激励器卷吸掺混能力, 提升激励器控制能力积累基础.      

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理,完善等离子体激励器数学模型,提升等离子体激励器扰动能力,采用粒子图像测速技术,在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究.实验时,将非对称布局激励器布置在平板模型上,随后将带有激励器的模型放置在有机玻璃箱内,从而避免环境气流对测试结果的影响.基于激励器诱导流场,分析了激励电压对诱导射流特性的影响,揭示了较高电压下诱导射流近壁区的拟序结构,获得了卷起涡、二次涡等拟序结构的演化发展过程,计算了卷起涡脱落频率,阐述了卷起涡与启动涡的区别,初步探索了卷起涡的耗散机制.结果表明:(1)层流射流不能完全概括等离子体诱导射流特性,激励电压是影响射流特性的重要参数.当电压较低时,诱导射流为层流射流;当电压较高时,诱导射流的雷诺数提高,射流剪切层不稳定,层流射流逐渐发展为湍流射流.(2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构;在固定电压下,这些涡结构存在恒定的卷起频率.(3)当激励电压较高时,流动不稳定使得卷起涡发生了拉伸、变形,引起了流场湍动能增大,从而加速了卷起涡的耗散.研究结果为全面认识激励器射流特性,进一步挖掘激励器卷吸掺混能力,提升激励器控制能力积累基础.     

深海多立柱浮式平台涡激运动特性研究

基于不可压缩纳维-斯托克斯(Navier-Stokes)方程和改进的延迟脱体涡模拟方法(improved delayed detached edd ysimulation method,IDDES),数值研究了深海多立柱浮式平台在海流作用下的涡激运动响应特性.以张力腿平台为对象,计算了该平台在0°,22.5°和45°流向角下的横向和艏摇涡激运动响应,分析了涡激运动响应幅度、频率比随约化速度的变化规律,研究了涡激运动响应能量的分布趋势.数值预报结果与模型实验数据吻合良好,证实了数值模型的有效性;研究发现,当约化速度介于7.0和14.0之间时,横向运动发生锁频,运动幅值稳定在0.2D~0.4D(D为立柱宽)之间,而艏摇涡激运动和约化速度呈线性递增关系;在横向运动锁频区间内,由于艏摇激振力矩主要受升力主导,艏摇频率与横向运动频率相同;相对于0°来流,22.5°和45°流向下的涡激运动频率更高,但艏摇运动能量仅为0°流向角下的10%.基于计算结果,进一步分析了多立柱平台涡激运动中的三维流场结构.      

具有壁面射流激励的圆管内气相流动大涡模拟

采用大涡模拟方法和单方程亚格子模式对小尺度量进行模拟。研究了不同强度壁面射流激励对圆管内气相流动的影响,模拟结果给出了射流对瞬态拟序结构发展、时平均流向速度分布的影响。随着射流强度的增加,射流入口附近流体的回流现象增强。射流强度足够大时可以减小管壁处的切应力值,同时会减小壁面附近流动速度,这种速度分布会导致气体夹带颗粒的能力下降,从而在实际两相流动中容易造成壁面附近的气粒返混现象。     

热激励器对超声速圆管射流的控制机理

采用大涡模拟方法研究了热激励器对马赫1.3 超声速圆管射流涡结构的影响, 采用加入净热源的方法模拟了热激励器的热效应, 分析了不同激励模态(m =±1 和m =±4) 的降噪和增加射流掺混的效果. 研究发现: (1)射流在受到激励作用时, 产生了更大的径向和周向速度扰动, 这有利于流向涡的形成, m =±1 激励模态下流向涡的发展更快, 特别是在摆动面上会有更大的增长速度; (2) 气动噪声分布的频带很宽, 激励作用对抑制高频气动噪声有较强的作用; (3) 对于在喷管壁面处产生的固壁噪声, m =±4 模态下的激励作用一定程度上加强了固壁噪声.      

逆壁射流中射流剪切层内湍流结构特性的实验研究

采用粒子图像测速技术对逆壁射流全流场进行了实验测量,射流与主流的速度比为8.89,基于射流圆管内径的雷诺数为对射流中心线上不同流向位置的脉动速度场统计分析发现:在占据主导地位.在驻点附近(尺度进行分析,在x/D=46～51总尺度向射流下游发展呈减小趋势.在x/D=35～41,参考点下游尺度大于上游尺度.在本征正交分解方法对湍流结构进行了定量分析,发现模态能量集中在低频,流场中能量最大的模态频率为f D/Uj=0.000 5,出现在再循环区.频率为产生了湍流结构,并且沿再循环区外围输运.高频结构的构型是类似的,均位于射流剪切层内,且频率越高,越接近射流出口,尺度越小.     

高雷诺数湍流风场大涡模拟的并行直接求解方法

在环境流体力学中,风场是风沙流、风雪流等自然环境特性问题研究的动力源和基础.通常采用壁湍流模型进行风场大涡模拟(large eddy simulation, LES)计算,但受到计算规模的限制使得高雷诺数风场的模拟计算难以实现.并行计算技术是解决大规模高雷诺数风场大涡模拟的关键技术之一.在不可压湍流风场的LES模拟中,压力泊松方程的并行计算技术是进行规模并行计算的困难点.根据风场流动模拟计算的特点,采用水平网格等距而垂直于地面网格非等距,在解决规模并行计算中求解压力泊松方程的难点问题时,利用FFT解耦三维泊松方程使其变为垂向的一维三对角方程,并利用可并行的三对角方程PDD求解技术,可建立三维泊松方程的直接并行求解技术.结合其它容易并行的动量方程计算,本文建立风场LES模拟的并行直接求解方法 (parallel direct method-LES, PDM-LES).在超级计算机上对新方法进行并行效率测试,并行计算效率达到90%.新的方法可用于进行湍流风场大涡模拟的大规模并行计算.计算结果表明,湍流风场瞬时速度分布近壁面存在条带状的拟序结构,平均场的速度分布符合速度对数律特性,风场湍流特性基本合理.     

微型射流涡流器对湍流边界层控制的数值研究

运用数值方法,模拟出展向分布的同向倾斜微型射流列与平板湍流边界层相互作用形成流向涡列的流场结构,验证了利用其来对湍流边界层进行控制的可能性.随射流间距减小,流向涡列控制作用流向渗透能力增强,但作用区域减小;随射流速度提高,流向涡列控制作用增强,但过大的射流速度反而会导致流向涡列在局部区域内控制作用的下降;随射流俯仰角减小、倾斜角增大,流向涡列初始控制作用增强,但过小的俯仰角、过大的倾斜角会导致流向涡列流向控制区域明显缩小.要保证流向涡列具有较强的湍流边界层控制作用,必须通过合理配置射流列各主要参数,在保证各流向涡具有一定强度的同时,还要确保各流向涡在形成时部分嵌入边界层内部.     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

两自由度振动系统的斜碰撞分析

研究斜碰撞振动系统动力学的一个关键问题是对系统在碰撞前后的状态进行合理描述和正确计算．针对两弹性体斜碰撞问题，基于瞬间碰撞假设，提出了采用步进冲量来分析和求解斜碰撞前后的状态关系；并以弹簧摆和振子组成的两自由度斜碰撞振动系统为例，具体介绍了该算法如何实现．用解析方法讨论了该系统在斜碰撞过程中可能出现的各种力学现象，将冲量步进算法得到的数值解与解析结果进行对比，取得了完全一致的结果．该数值方法能适应多种斜碰撞问题的计算．     

Application of the theory of branching of solutions of nonlinear equations in the estimation of pole perturbance of linear closed-loop systems

When the state and input matrices of a multivariable linear time-invariant system are perturbed, the problem of the estimation of pole perturbance of closed-loop system is considered by making use of the theory of branching of solutions for nonlinear equations.     

The distribution of zeroes of solutions of neutral equations

1.IntroductionThestudyofoscillator}'behaviorofneutraldifferentialequationsisarelativealnewfieldandisveryinterestinginapplications.Inrecentyears.theoscillatorytheoryfortheseequationshasbeenextensivelydeveloped.WerefertotherecentbookbyGyoriandLadaSI'I,Baino…     

CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

SWT-120风洞稳定段的性能测量

在普通超音速风洞中,由于受到噪声干扰,很难进行有效的边界层转捩特性试验研究,针对超声速流动特点发展较低噪声风洞十分必要,而稳定段设计的好坏直接影响到下游试验段噪声水平。本文介绍一座低噪声风洞稳定段的结构和性能测量结果。先对稳定段的结构设计做了简单介绍,然后对实验结果进行分析,实验结果表明在大角度扩散段内装置孔锥,稳定段安装消音夹层和阻尼网组等部件后,气流的速度脉动和压力脉动明显降低,其中压力脉动降低一个量级,速度脉动为1%。进一步优化设计和改进工艺,速度脉动还可进一步降低。测量结果表明SWT-120稳定段的设计是成功的,对我们以后发展更高性能的静风洞有借鉴和参考价值。     

综合管廊早龄期混凝土温度应力仿真模拟

以汉中市综合管廊为工程背景,通过有限元软件,首先模拟分析综合管廊早龄期温度场,随后将温度场结果附加到管廊结构中求解温度应力,从而分析综合管廊早龄期混凝土温度应力与浇筑长度的关系。模拟分析发现:综合管廊在混凝土浇筑早期,温度场沿长度方向变化不显著,且随着浇筑长度的增加,温度场不变;在墙体温度上升阶段,墙体内部沿长度方向的应力主要为压应力,当墙体混凝土内部温度达到最高,混凝土开始降温,结构内变形由膨胀逐渐变为收缩,内力由压应力逐渐变为拉应力,随着时间的变化,墙体内外温差趋于稳定,墙体温度应力也逐渐平稳,此时温度应力主要受外界环境变化影响;当综合管廊浇筑长度在60m范围内,浇筑长度越大,温度应力越大,但当浇筑长度超过60m,温度应力不再增加。     

One type of solutions of conduction of nervous pulses and the stability of action potential

This paper reports a type of laws which governs action poten-tial of nervous impulses,and it is discussed by general form-nonlinear dispersive process.We find that the nervous wave is aslowly varying amplitude solitary wave in the small dispersivecase.we prove that the solitary wave is not generated in the or-dinary dispersion,but a travelling wave with varying amplitudesmay be obtained.The stability of various possible action poten-tials and bifurcation in overdamped case are also discussed inthis paper.     

Evolution of viscoelastic properties of polyurethane in the course of curing

Simultaneous measurements of the components of dynamic modulus at the set of frequencies were made by Fourier-Transform Mechanical Spectroscopy (FTMS) to investigate the curing process of a segmented polyurethane.The whole process of curing is treated as consisting of three stages. In the first stage, the material can flow and this stage is complete at the gel-point. In the second stage, the curing proceeds in the network and ends with phase separation of aggregating segments (blocks) of polymeric chains. In the third stage, the curing continues in the two-phase system. The object of investigation was the transition from the second to the third stage. It was shown that the position of the maximum of loss tangent depends on frequency and does not coincide with different characteristic points marking the time dependencies of G and G. This is due to the formation of micro-particles of the newly formed phase. The FTMS method is used for separation of the modulus changes into parts related to formation of chemical and physical networks.Special means of treating calorimetric measurements identified a weak maximum heat output at the moment corresponding to the transition point.