非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧,提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步,在粗网格上采用Oseen迭代法求解非线性问题,在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散,时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.     

沙丘背风侧不同铁路结构形式对风沙环境的适应性分析

当铁路穿越大风沙漠地区时,风沙灾害会对铁路工程及其正常运营产生严重威胁,而设计一种合理的铁路结构形式能够减小风沙沉积对铁路工程的危害. 在本文中,以敦煌至格尔木铁路沙山沟段落为研究对象,采用多相流的方法对越过沙丘的风沙运动过程进行数值模拟,分别讨论了风沙运动对位于沙丘背风坡的铁路路基工程和桥梁工程的影响. 主要的模拟结果显示:路基工程明显降低了风速并且将沙丘后的回流区分成了两部分,而桥梁工程的导流效应则压缩了沙丘背风坡的回流区;轨道间的道碴增大了铁路表面的粗糙度,在轨道间有少量沙粒沉积,而路基工程两侧则有大量积沙;铁路表面的积沙量与摩阻风速呈现出非线性关系,随着摩阻风速的增大,路基工程沙粒沉积的增加速度大于风蚀能力的增加速度,而桥梁工程则正好相反. 在防止风沙危害铁路方面,设置桥梁工程明显优于路基工程. 本研究为风沙运动对铁路工程的影响提供了理论支持,也为今后的铁路工程设计提供了新的思路与研究工具.      

攒尖四坡屋面风致雪漂移的数值模拟

为获得攒尖四坡屋面的风致雪漂移规律,基于欧拉-欧拉方法和风雪单向耦合假定,运用计算流体动力学软件,选取Mixture模型分别对立方体周边和高低屋面上的风致雪漂移运动进行数值模拟,将模拟结果与两者的实地观测数据对比,探讨分析数值风洞的关键技术和参数设置,验证数值模拟方法的合理性与可靠性。依据攒尖四坡房屋的使用功能要求,设计分析模型与分析工况,在试算的基础上对屋面进行分区。以风速5 m/s,7 m/s,9 m/s,11 m/s,13 m/s和15 m/s,风向角0°,15°,30°和45°以及屋面坡度25°,30°,34°,40°和45°为分析参数,对攒尖四坡房屋的120种工况进行数值模拟,得到屋面各分区侵蚀沉积的基本规律,提出可用于抗雪设计的屋面积雪分布系数。研究表明,风向角的改变会使屋面积雪分布状态发生较大程度的变化,风速和屋面坡度的变化对屋面整体积雪量有较大影响。     

燃气透平叶片表面颗粒沉积特性数值研究进展

污染物颗粒在燃气透平叶片表面及冷却通道内的沉积将影响透平叶片的冷却特性和安全服役寿命。国外许多研究机构对燃气透平叶片表面污染物颗粒沉积特性及其与叶片气膜冷却的相互作用机制进行了持续研究。本文首次系统总结了近年来燃气透平叶片表面污染物颗粒沉积特性数值研究进展,重点介绍了燃气透平叶栅流道内污染物颗粒的沉积机制、颗粒沉积和脱离物理模型以及颗粒沉积与透平叶片冷却特性相互作用的数值模拟方法和研究成果。基于当前研究热点和发展趋势,结合作者的研究经历,指出了国内在研发先进燃气轮机透平冷却结构时应注意的问题和方向,为高效、燃料适应性更广的燃气轮机技术自主开发奠定基础。     

An Insight into the Diffusion of Unsaturated Polyester (UP) Resin into Expanded Graphite (EG) to Improve the Mechanical Properties of the UP/EG Composites

Abstract

Polymer/expanded graphite (EG) nanocomposites have great importance in many industrial applications mainly due to their high electrical/thermal conductivity or flame retardancy. However, to fully employ the benefits of polymer/EG nanocomposites one must consider the high degree of porosity of EG. The high degree of porosity of EG can deteriorate the composites’ mechanical properties if the polymer chains cannot diffuse completely into the EG pores. In this article, an insight is given into the diffusion of unsaturated isophthalic polyester (UP) resin, consisting of a combination of maleic anhydride and isophthalic anhydride in the resin backbone, with two viscosities, into the pores of the EG particles of various degrees of porosity. The diffusion experiments were carried out on compressed EG tablets with the same density but different porosity due to the different porosity of the EG particles. The results showed that the diffusion rate of the UP resin with higher viscosity slightly decreased when the EG porosity decreased but, in the opposite way, it strongly increased for the low viscosity UP resin. The EG nanocomposites samples were molded at varying pressures. The micrographs of the fractured surfaces of the EG nanocomposites showed that the EG pores were not filled with resin, thus the EG nanocomposites had residual pores. It was found that composites containing EGs with higher expansion ratio and larger particles and pores showed larger residual pores. Furthermore, the composites prepared with the more viscous UP resin showed more residual pores. By applying a pressure of 10?bar instead of 1?bar, a reduction of 7–20% in the residual pores of the nanocomposites was observed which led to improved mechanical properties by up to 20% in flexural strength for the EG with the highest expansion ratio.     

吸收峰混叠的太赫兹光谱区间组合特征提取算法

太赫兹光谱是物质识别的前沿方法之一。由于不同物质的分子组成或结构各异,许多物质的太赫兹吸收谱会在特定频率上出现吸收峰,可以作为混合物成分检测的重要特征。有效准确地提取这些吸收峰的参数,是提高识别率的关键。多峰拟合算法将光谱曲线拟合成若干个标准峰函数之和,能够同时提取到吸收峰的频率、峰高、峰宽等信息。但是该算法以寻峰算法结果为基础确定吸收峰的大致位置和数量,寻峰结果不一定是最优的拟合结果,而且很难准确识别定位混叠状态的吸收峰。为了提高混叠光谱中吸收峰的识别定位精度,提出以大幅度平滑后的曲线波谷为分界点,将预处理后的光谱分成若干个子区间。然后将子区间组合起来进行多峰拟合,通过遗传算法得到最优的拟合子区间组合和吸收峰频率近似值,拟合时每个子区间中通过峰数递增最优化方法确定拟合的吸收峰数,最后微调优化得到最优的吸收峰频率、峰高值。为了实现物质的识别,通过密度聚类算法得到同一类纯净物在多次测量中的共同吸收峰,以此作为标准数据,通过提出的基于吸收峰特征的光谱匹配算法实现了纯净物和不同含量混合物的快速识别。对10类纯净物的实际光谱数据进行拟合聚类,得到其吸收峰参数,结果与太赫兹光谱数据库一致。通过识别算法对纯净物测试集进行识别的识别率为100%,证明了特征提取和物质识别算法的有效性。对于含有混叠峰的混合物光谱,二阶导数法对葡萄糖-乳糖混合物光谱中被掩盖吸收峰(1.280 THz)的识别率仅为70%,提取到的频率平均值为1.316 THz;而该算法提高识别率至95%,频率平均值为1.281 THz,该算法提高了对混叠峰的分辨能力,能够精确定位混叠峰。对10类纯净物构成的6类不同程度混叠的二元混合物前二、三识别率分别达到90.8%和98.3%,提取到的特征能够有效应用于混合物的成分检测。该算法能够以纯净物数据为标准数据实现成分各异的混合物成分检测,对于太赫兹光谱混合物成分检测有重要意义。     

秋千参数自激振动的数学模型与数值模拟

将秋千视为变长度单摆,基于归一化高斯函数构造了摆球相对摆杆的速度表达式.在此基础上建立了秋千的非线性动力学微分方程,并据此数学模型对秋千运动进行了数值模拟,得到了秋千的运动规律和功能转化的机制.     

Modified Gerchberg–Saxton (G-S) Algorithm and Its Application

The Gerchberg–Saxton (G-S) algorithm is a phase retrieval algorithm that is widely used in beam shaping and optical information processing. However, the G-S algorithm has difficulty obtaining the exact solution after iterating, and an approximate solution is often obtained. In this paper, we propose a series of modified G-S algorithms based on the Fresnel transform domain, including the single-phase retrieval (SPR) algorithm, the double-phase retrieval (DPR) algorithm, and the multiple-phase retrieval (MPR) algorithm. The analysis results show that the convergence of the SPR algorithm is better than that of the G-S algorithm, but the exact solution is not obtained. The DPR and MPR algorithms have good convergence and can obtain exact solutions; that is, the information is recovered losslessly. We discuss the security advantages and verification reliability of the proposed algorithms in image encryption. A multiple-image encryption scheme is proposed, in which n plaintexts can be recovered from n ciphertexts, which greatly improves the efficiency of the system. Finally, the proposed algorithms are compared with the current phase retrieval algorithms, and future applications are discussed. We hope that our research can provide new ideas for the application of the G-S algorithm.     

Dynamics of a Predator-Prey Model with Allee Effect and Herd Behavior

This paper deals with dynamics of a predator-prey model with Allee effect and herd behavior. We first study the stability of non-negative constant solutions for such system. We also establish the existence of Hopf bifurcation solutions for such predator-prey model. The stability and bifurcation direction of Hopf bifurcation solution in the case of spatial homogeneity are further discussed. At the same time, several examples are given by MATLAB. Finally, the numerical simulations of the system are carried out through MATLAB, which intuitively verifies and supplements the theoretical analysis results.     

Structural Entropy of the Stochastic Block Models

With the rapid expansion of graphs and networks and the growing magnitude of data from all areas of science, effective treatment and compression schemes of context-dependent data is extremely desirable. A particularly interesting direction is to compress the data while keeping the “structural information” only and ignoring the concrete labelings. Under this direction, Choi and Szpankowski introduced the structures (unlabeled graphs) which allowed them to compute the structural entropy of the Erdős–Rényi random graph model. Moreover, they also provided an asymptotically optimal compression algorithm that (asymptotically) achieves this entropy limit and runs in expectation in linear time. In this paper, we consider the stochastic block models with an arbitrary number of parts. Indeed, we define a partitioned structural entropy for stochastic block models, which generalizes the structural entropy for unlabeled graphs and encodes the partition information as well. We then compute the partitioned structural entropy of the stochastic block models, and provide a compression scheme that asymptotically achieves this entropy limit.