流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

浮力效应对射流扩散火焰转捩特性影响的实验研究

对不同空气伴流速度下丙烷层流射流火焰向湍流火焰的转捩过程进行实验观测，分析伴流对火焰转捩行为及稳定性的影响.相对于静止环境中的射流火焰，较大速度的伴流可以减小浮力效应对射流扩散火焰转捩过程的影响，使火焰发生转捩的临界Reynolds数（Re_cr）增大，即火焰推迟转捩.但当伴流速度较小时，Re_cr保持不变，转捩过程中的射流火焰发生周期性振荡，振荡幅度随着伴流速度的增大而减小，继续增大伴流速度，火焰振荡的周期性最终消失，转而呈现随机性.实验还发现，喷嘴直径较大的扩散火焰的Re_cr更大.考虑到火焰对燃料射流局部流动状态的影响，对此现象进行了解释.      

Cahn-Hilliard-Schlgl反应-扩散方程分叉问题

本文利用Lyapunov-Schmidt约化方法、奇异性理论及摄动方法,对服从Cahn-Hilliard广义扩散定律的带有Schlgl反应项的Cahn-Hilliard-Schlgl反应-扩散方程的分叉情况进行了研究。研究结果表明,C-H-S方程可以产生单模态的空间周期一次和二次分叉。特别地,经小扰动的退化的C-H-S方程,在不同条件下,还可能产生不同形式的更为复杂的双模态混合的二次分叉。本文还给出了分叉类型和一次分叉解的近似表达式。     

煤粉燃烧过程的数值模拟

采用计算流体动力学软件对煤粉实际燃烧过程进行数值计算并结合其热态试验数据进行对比分析,已成为验证数学模型和指导工程实践的一种重要研究手段.应用商业软件FLUENT对煤粉燃烧及流场进行了数值模拟分析,采用混合分数/概率密度函数法模拟湍流燃烧,用P-1辐射模型开展辐射传热模拟,利用拉格朗日/欧拉法处理气固两相间的动量、质量和能量交换,对挥发份的析出采用双速率竞争模型,采用动力/扩散反应速率模型模拟煤粉颗粒的表面燃烧.并对模拟结果进行分析与对比,计算结果与实际燃烧过程有较好的一致性.数值模拟计算结果揭示了挥发分释放与焦炭燃烧的过程,为改善和优化煤粉的燃烧提供了重要的参考依据.     

湍流燃烧反应机理的简化

采用速率敏感度分析法分别对燃料为CO/H2和CH3OH的钝体绕流湍流扩散火焰的燃烧机理进行了简化.判断出了冗余组分,并通过对敏感度矩阵F的主成分分析,找出了次要反应,使原始的复杂湍流燃烧机理得到了大幅度的简化.简化机理与详细机理的计算结果吻合得很好.结果表明,速率敏感度分析的方法能大大简化复杂的湍流反应机理,且计算效率高,精度满足要求.     

基于非Favre过滤的大涡模拟

将直接过滤的可压缩Navier-Stokes(NS)方程组作为可压缩湍流大涡模拟方程组.方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,得到可压缩张量扩散模型.与传统的采用质量加权过滤(Favre过滤)的可压缩大涡模拟方程组相比,由于这种方法没有密度信息的损失,所以能反映真实的密度脉动信息;方程所求量为流场物理量的过滤值,物理意义明确,使大涡模拟结果与直接数值模拟或实验结果更具有可比性.用这种大涡模拟方程组对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验、Baldwin-Lomax(BL)模型及Smagorinsky亚格子模型的结果进行了对比.      

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

翼型湍流尾缘噪声半经验预测公式改进

改进了传统的翼型湍流边界层尾缘噪声BPM半经验预测公式。传统的BPM半经验湍流边界层尾缘噪声预测公式对高攻角和厚翼型在高频范围的预测结果大于实验结果,通过分析比较传统BPM半经验预测公式和Howe翼型尾缘噪声理论模型发现:这主要是由于传统BPM半经验预测公式对压力面声源噪声辐射高估引起的。因此将压力面声源噪声辐射与吸力面声源噪声辐射的幅值比由原来的边界层位移厚度一次方比值改进为二次方比值,进而得到了改进后的BPM半经验预测公式;使用改进后的BPM半经验预测公式对NACA0012翼型在不同来流不同攻角下的噪声辐射进行了预测比较,发现对于NACA0012翼型,改进后的BPM半经验预测公式具有较高精度;另外也预测了较厚的风力机翼型DU-96-W-180,预测结果明显改善。     

Conditional sea surface statistics and their impact on geophysical sea surface parameters retrieved from SAR altimetry signals

In this study we present an analytical formulation of synthetic-aperture radar (SAR) altimetry signals including narrow banded nonlinear wave fields and conditional statistics between wave elevation displacements, horizontal wave slopes and vertical wave particle velocities. Considering the wave elevation displacements coskewness with respect to horizontal slopes leads to an analytical formulation of the electromagnetic bias within a SAR-mode altimeter stack. This formulation can be either parametrized by the significant wave height (SWH) and mean wave steepness, or in terms of the variance of vertical wave velocities. The effect of conditional vertical wave particle velocity variances with respect to the observed horizontal wave slopes close to nadir incidence angles leads to an effective reduction of the azimuth blurring of SAR-mode stacks. We present here a formulation of this effect by examining JONSWAP ocean wave spectra. In most cases this effect reduces the azimuth blurring by 10% to 30%. Additionally we investigate the effect of a nonlinear wave elevation displacement probability density function (PDF) on estimated geophysical parameters. We were able to show that including an elevation displacement skewness of 0.13 improves significantly the SWH consistency between altimetry and ECMWF Reanalysis v5 ERA5 results.All of these effects are validated with respect to ERA5 model data in the North East Atlantic region and in situ data located in the German Bight and Baltic Sea.The developed model can be used in both SAR and conventional altimetry retrackers.     

微重力效应的物理解释及其应用

讨论了微重力效应及其分子物理学解释，万有引力的单极性决定了它的累积性和宏观性，重力对自由分子的影响是微科其微的，然而在流钵，特别是在液体内部的具体条件下，重力对流体内部分子集团产生明显影响，以宏观的静压强形式表现出来。而当重力变小，变微，流体内部的静压强也就趋于消失，由静压强产生的二级现象，如异相沉降或浮泛。同相的浮泛对流，液体自约束成球形，也就消失，这些就是微重力的一，二级物理效应，本文认为，Bossinesq近似对于流体而言实质上是线性热力学假设或近平衡态假设，可以用初始平衡态的物性参量表征临界点时的无量纳数，提出一种重力消失诱发非平衡态向平衡态蜕变的猜想，从文中导出的瑞利数Rα的表达式可以看出，重力加速度的消失，可以使瑞利数的值从远离平衡态的湍流蜕变到稳定流动的瑞利－贝纳对流胞，甚至是没有宏观流动的平衡态，最后，将本文给出的理论解释模型应用到空间材料加工所需微重力水平的估算、窨材料加工工艺的微重力利用准则和空间装置中气体自然对流传热状态模拟的低真空和微重力的互换性方面。     

On the heating mechanism of magnetic flux loops in the solar atmosphere

Hα filtergrams and magnetograms indicate that bright features (such as plages and granulation boundaries) correspond to areas of strong vertical magnetic fields and dark features (such as fibrils and filaments) are associated with strong horizontal magnetic field. It was suggested by /1/ that there is an excess dissipation of waves, available for heating, in regions of vertical magnetic fields. With this suggestion in mind, we have investigated the physical heating mechanism due to ponderomotive forces exerted by turbulent waves along curved magnetic flux loops. Results show that the temperature difference (ΔT) between the inside and outside of the flux loop can be classified into three parts; ΔT = ΔT₁ + ΔT₂ + ΔT₃; in which ΔT₁ and ΔT₃ represent the heating or cooling effect from the ponderomotive force, and ΔT₂ is the heating effect due to conversion of turbulent energy from the localized plasma. The specific physical mechanism (i.e., the ponderomotive forces exerted by turbulent waves), is used to illustrate solar atmospheric heating via an example leading to the formulation of plages.     

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

环形桁架式天线展开不同步现象分析

从国外相关文献所反映的地面试验结果看,环形桁架式天线周边的桁架结构由于阻尼等因素在展开过程中存在明显的不同步现象。文章利用绝对节点坐标法建立了天线展开过程的多柔体展开动力学模型,获得了连续展开过程中桁架各部分承受的载荷和变形情况,并重点对展开过程中的不同步现象的产生原因及其影响进行了分析。结果表明,环形桁架式天线展开的不同步现象是由展开过程中桁架的变形与摩擦引起的,这种现象的出现导致桁架结构承受更高的载荷。该研究为环形桁架式天线的设计和验证提供了参考。     

后掠翼身干扰区流动特性及改善措施研究

利用流动显示及表面压力测量方法研究了后掠翼身干扰区的流动特性,并研究了用小边条等措施改善干扰区的流动特性的效果.结果表明,随着不同机翼后掠角、不同迎角及不同Re数对干扰区流动特性的影响,流态可以从一涡系变成多涡系,由定常变成非定常,而且在一定的Re数以后涡系会紊流化;翼身干扰区上游的的逆压梯度是导致边界层分离的物理原因,利用面积很小的边条可以降低干扰区局部的逆压梯度,可以导致干扰区的旋涡很弱,甚至不出现,这是很有实际意义的.     

利用热力学与统计物理方法研究电离层E层的形成

将太阳光子作为反应物引入电离层化学反应中, 运用化学热力学方法, 获得低电离层光化学反应的反应度及其化学平衡常数, 进而得到低电离层的电子密度, 并通过统计物理方法推导出相应的化学平衡常数. 结果显示, 运用热力学方法得到的低电离层电子密度与Chapman理论几乎完全一致, 表明利用热力学方法研究电离层的形成是可行的. 分别用热力学方法和统计物理方法得到的平衡常数之间存在一定的差异, 分析表明, 这是因为光化学过程中并非所有光子能量都转化成了电离能, 其中有一部分能量转化为热能或使粒子处于激发态, 这种差异的存在使利用统计物理方法推导电离层等离子体温度成为可能.      

氢氧发动机推力室化学反应流场计算

采用弱耦合点隐式方法的MacCormack格式对氢氧火箭发动机推力室化学反应粘性流场进行数值模拟.数值模拟时采用6种组分、8个反应有限速率的化学反应模型,湍流模型采用Baldwin-Lomax代数湍流模型.数值模拟得到了流场参数在燃烧室和喷管中的分布.结果分析表明,得到的数值模拟结果与理论分析一致,说明结果是可靠的.本文的工作为氢氧火箭发动机喷管设计提供了依据,并为进一步开展火箭发动机推力室有化学反应的两相流动的数值模拟打下了基础.