货运动车组装载门凹陷对隧道内会车压力波的影响

随着我国技术经济水平的提高，电子商务产业快速发展，货运动车组应运而生。采用有限体积方法和重叠网格方法，以及SST k-ω高雷诺数湍流模型，数值模拟研究时速350 km的货运动车组在隧道内交会的压力波，探究货运动车组车体结构即装载门凹陷对隧道内交会压力波的影响。实车试验结果验证了本文数值模拟方法的准确性。研究结果表明，凹陷的装载门使平直车身处的最大正压值增大约20%~30%，且有无凹陷的装载门中心最大正压值的差异百分比维持在35%~80%，最大负压值差异百分比在20%~25%。     

A new endocytic model based on the co-rotational grid method

Endocytosis plays important roles in many cellular physiological processes, such as metabolism and immune. Many theoretical models have been proposed to study the endocytic process, but little has considered the tensile deformation of the membrane and the actin forces. In this paper, a new endocytic model is proposed based on the co-rotational grid method. Our model gives a direct estimation of the in-plane strain of the plasma membrane and provides a basis for the calculation of further scission process of the vesicle. The results fit well with experimental data in the literature. Moreover, it is suggested that the active forces of actin at the endocytic site is the main mechanism driving the invagination of the plasma membrane.     

Characteristic behaviour of plasma fluctuations inside a plasma bubble in presence of a magnetic field due to the formation of a potential well

Experiments were carried out to observe the effect of a magnetic field and grid biasing voltage in presence of a plasma bubble in a magnetized, filamentary discharge plasma system. A spherical mesh grid of 80% optical transparency was negatively biased and introduced into the plasma for creating a plasma bubble. Diagnostics via an electrical Langmuir probe and a hot emissive probe were extensively used for scanning the plasma bubble. Plasma floating potential fluctuations were measured at three different positions of the plasma bubble. The instability in the pattern showed the dynamic transition from periodic to chaotic for increasing magnetic fields. Time scale analysis using continuous wavelet transform was carried out to identify the presence of non‐linearity from the contour plots. The mechanisms of the low‐frequency instabilities along with the transition to chaos could be qualitatively explained. Non‐linear techniques such as fast Fourier transform, phase space plot, and recurrence plot were used to explore the dynamics of the system appearing during plasma fluctuations. In order to demonstrate the observed chaotic phenomena in this study, characteristics of chaos such as the Lyapunov exponent were obtained from experimental time series data. The experimentally observed potential structure is confirmed with numerical analysis based on fluid hydrodynamics.     

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理，给出了八节点六面体等参实体单元的基本公式，通过投影方式将空间曲面六面体转换为投影六面体，采用参考面的逆向运动求解节点纯变形，通过单元形函的虚功方程计算节点内力；针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序，并进行工程结构算例分析。结果表明，所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析，验证了本文理论和程序的正确性和实用性。     

New algorithms for discrete vector optimization based on the Graef-Younes method and cone-monotone sorting functions

Abstract

The well-known Jahn-Graef-Younes algorithm, proposed by Jahn in 2006, generates all minimal elements of a finite set with respect to an ordering cone. It consists of two Graef-Younes procedures, namely the forward iteration, which eliminates a part of the non-minimal elements, followed by the backward iteration, which is applied to the reduced set generated by the previous iteration. Without using the backward iteration, we develop new algorithms that also compute all minimal elements of the initial set, by combining the forward iteration with certain sorting procedures based on cone-monotone functions. In particular, when the ordering cone is polyhedral, computational results obtained in MATLAB allow us to compare our algorithms with the Jahn-Graef-Younes algorithm, within a bi-objective optimization problem.     

A 3D mesh deformation technique for irregular in‐flight ice accretion

Aircraft holding around busy airports may be requested to sustain as much as 45 min of icing before landing or being diverted to another airport. In this paper, a three‐dimensional mesh deformation scheme, based on a structural frame analogy, is proposed for the numerical simulation of ice accretion during extended exposure to adverse weather conditions. The goal is to provide an approach that is robust and efficient enough to delay or altogether avoid re‐meshing while preserving (enforcing) nearly orthogonal elements at the highly distorted ice surface. Robustness is achieved by suitably modifying the axial and torsional stiffness components of the frame elements in order to handle large and irregular grid displacements typical of in‐flight icing. Computational efficiency is obtained by applying the mesh displacement to an automatically selected small subset of the entire computational domain. The methodology is validated first in the case of deformations typical of fluid‐structure interaction problems, including wing bending, a helicopter rotor in forward flight, and the twisting of a high‐lift wing configuration. The approach is then assessed for aero‐icing on two swept wings and compared against experimental measurements where available. Copyright © 2015 John Wiley & Sons, Ltd.     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法，采用涡粘性模型求解湍流熵产方程，系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时，0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性，验证了算法及程序的正确性。结果表明，选择不同湍流模型时，翼型流场熵产的计算结果存在差异，湍流耗散是引起流场熵产的主要原因；翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域，流场熵产率与翼型阻力系数线性相关；当产生分离涡时，粘性耗散引起的熵产下降。     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象，发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式，避免在单元内进行速度平均，从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比，网格的边界通量由节点处的多维黎曼求解器提供，克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增，给出稀疏波出现的合理判据，从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象，同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.     

基于改进格林元的压裂水平井动态分析方法

以建立高效的动态分析方法为出发点，以边单元作为求解点，改进传统的格林元方法，减少未知数和求解矩阵维度；并提出基于改进格林元的加密网格加密方法，保证考虑复杂裂缝网络的压裂水平井动态模拟的早期精度.退化模型与半解析解、数值模拟结果进行对比，验证本文基于加密网格的改进格林元方法的准确性和动态分析的高效性.最后进行动态响应的敏感性分析，结果表明：①格林元方法是一种高精度的动态模拟方法，将求解节点设置在网格的边上可以提高压裂水平井动态模拟的速度；②改进格林元方法的加密基于叠加原理，不需要通过插值近似，其求解精度高.在相同加密网格条件下，基于本文改进格林元方法的加密效果比有限差分加密效果更佳；③复杂裂缝导流能力、改造区渗透率提高倍数、改造区大小等参数对压裂水平井动态特征影响较大，在动态分析和参数反演时，应着重考虑这些因素的影响.