Numerical simulation of coherent structure in two-dimensional compressible mixing layers

The coherent structure in two-dimensional mixing layers is simulated numerically with the compressible Navier-Stokes equations. The Navier-Stokes equations are discretized with high-order accurate upwind compact schemes. The process of development of flow structure is presented: loss of stability, development of Kelvin-Helmholtz instability, rolling up and pairing. The time and space development of the plane mixing layer and influence of the compressibility are investigated.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting. and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However. the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C(t) parameters. so that crack growth rates correlate rather well with C(t) and C(t). A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.

A study of the mechanism of consolidating metal-powder under explosive-implosive shock-waves

A study of the two-dimensional flow pattern of particles in consolidation process under explosive-implosive shock waves has been performed to further understand the mechanism of shock-wave consolidation of metal powder, in which bunched low-carbon steel wires were used instead of powder. Pressure in the compact ranges from 6 to 30 GPa. Some wires were electroplated with brass, some pickled. By this means, the flow pattern at particle surfaces was observed. The interparticle bonding and microstructure have been investigated systematically for the consolidated specimens by means of optical and electron microscopy, as well as by microhardness. The experimental results presented here are qualitatively consistent with Williamson's numerical simulation result when particle arrangement is close packed, but yield more extensive information. The effect of surface condition of particle on consolidation quality was also studied in order to explore ways of increasing the strength of the compacts. Based on these experiments, a physical model for metal powder shock consolidation has been established.

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法.首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract,1993,62:203～218)相比较,表明该方法具有很高的精度和很好的计算效率.随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则.进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.

迎风紧致格式的混淆误差分析及其同谱方法的比较

对运用迎风紧致格式求解非线性方程时混淆误差产生的机理进行了研究，通过算例对五阶迎风紧致格式与谱方法进行了比较，发现在混淆误差的处理上迎风紧致格式优于谱方法。

粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型，将高精度对称紧致格式引入边值法数值线性稳定分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法，实现了时间模式和空间模式的统一求解，并将扰动特征及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性，涉及二维/三维扰动波、粘性/无粘扰动波、第一/第二模态、特征函数、伪特征值谱等。研究表明，压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用；在Mc = 1附近，从高波数段开始，粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡。

Formation and Structure of Composite Coating of HAD and Micro-Plasma Oxidation on A3 Steel

Composite coatings were obtained on A3 steel by hot dipping aluminum(HAD) at 720 degreesC for 6 min and micro-plasma oxidation (MPO) in alkali electrolyte. The surface morphology, element distribution and interface structure of composite coatings were studied by means of XRD, SEM and EDS. The results show that the composite coatings obtained through HAD/MPO on A3 steel consist of four layers. From the surface to the substrate, the layer is loose Al2O3 ceramic, compact Al2O3 ceramic, At and FeAl intermetallic compound layer in turn. The adhesions among all the layers are strengthened because the ceramic layer formed at the At surface originally, FeAl intermetallic compound layer and substrate are combined in metallurgical form through mutual diffusion during HAD process. Initial experiment results disclose that the anti-corrosion performance and wear resistance of composite coating are obviously improved through HAD/MPO treatment.

Structural Basis for Feed-Forward Transcriptional Regulation of Membrane Lipid Homeostasis in Staphylococcus aureus

11 p.

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

Generating cluster submodels from a multistage stochastic mixed integer optimization model using break stage

We present a scheme to generate clusters submodels with stage ordering from a (symmetric or a nonsymmetric one) multistage stochastic mixed integer optimization model using break stage. We consider a stochastic model in compact representation and MPS format with a known scenario tree. The cluster submodels are built by storing first the 0-1 the variables, stage by stage, and then the continuous ones, also stage by stage. A C++ experimental code has been implemented for reordering the stochastic model as well as the cluster decomposition after the relaxation of the non-anticipativiy constraints until the so-called breakstage. The computational experience shows better performance of the stage ordering in terms of elapsed time in a randomly generated testbed of multistage stochastic mixed integer problems.

可压缩流气动声场的数值模拟

气动声学是一门流动力学和声学之间的交叉学科，主要研究流动及其与物体相互作用产生噪声的机理。动用计算技术研究气动声学问题的手段称为计算气动声学。本文的目的是，基于高精度数值算法的研究，分别运用Lighthill比拟理论、Kirchhoff积分和直接数值模拟等方法，针对翼型绕流、激波-涡干扰和轴对称射流，研究了物面非定常脉动压力、涡脱落、激波-涡干扰以及涡对并等产生噪声的机理。首先针对声场与主流场在能级和特征尺度等方面的差异，从空间离散角度分析了几种差分格式，表明迎风紧致格式/对称紧致格式有较小的数值色散、耗散和各向异性误差，因而适用于气动噪声的计算。以Runge-Kutta格式为例，对时间离散带来的误差进行了分析。指出对声波计算来说，仅考虑格式稳定性是不够的，时间步长还受到允许色散误差和耗散误差的限制。基于保色戎关系的思想，构造了优化Runge-Kutta格式。处例显示优化Runge-Kutta格式相对于经典格式有更高的计算效率。采用3阶迎风紧致格式和3阶Runge-Kutta格式数值模拟了NACA0012翼型的可压缩非定常绕流流场，并将此流场作为近场声源，运用声学比拟理论对偶极子声和四极子声进行研究。结果指出，主流速度对远场声压有决定性影响，在来流马赫数较大时，四极子噪声和偶极子噪声具有相同量级，不能被忽略，表明了可压缩效应对声场的影响。采用5阶迎风紧致格式和4阶Runge-Kutta格式求解非定常可压缩Navier-Stokes方程，对激波-单涡/双涡干扰导致的声场进行了直接数值模拟。详细研究了激波-涡干扰产生噪声的机理，指出噪声的产生及其性质和激波变形密切相关。研究了近场噪声衰减和传播距离r的关系，发现噪声衰减大致和r~(4/5)而不是r~(1/2)成反比关系，提出这种差异是由流场的非线性效应引起的。构造了Kirchhoff积分和非定常流动计算相结合的算法。采用5阶迎风紧致格式和3阶Runge-Kutta格式对亚声速轴对称射流进行直接数值模拟。将射流流场作为近场声源，结合Kirchhoff方法求解远场 气动噪声。数值结果表明远场噪声具有方向性，噪声声压在离开对称轴20°处达到最大值。随着传播距离增大，噪声方向性逐渐减弱。

Sampling and learning the Mallows and Generalized Mallows models under the Cayley distance

[EN]The Mallows and Generalized Mallows models are compact yet powerful and natural ways of representing a probability distribution over the space of permutations. In this paper we deal with the problems of sampling and learning (estimating) such distributions when the metric on permutations is the Cayley distance. We propose new methods for both operations, whose performance is shown through several experiments. We also introduce novel procedures to count and randomly generate permutations at a given Cayley distance both with and without certain structural restrictions. An application in the field of biology is given to motivate the interest of this model.

槽道湍流的直接数值模拟

通过直接数值模拟(DNS)研究槽道湍流的性质和机理。包含五个部分：1）湍流直接数值模拟的差分方法研究。2）求解不可压N－S方程的高效算法和不可压槽道湍流的直接数值模拟。3）可压缩槽道湍流的直接数值模拟和压缩性机理分析。4）“二维湍流”的机理分析。5）槽道湍流的标度律分析。1．针对壁湍流计算网格变化剧烈的特点，构造了基于非等距网格的的迎风紧致格式。该方法直接针对计算网格构造格式中的系数，克服了传统方法采用 Jacobian 变换因网格变化剧烈而带来的误差。针对湍流场的多尺度特性分析了差分格式的精度、网格尺度与数值模拟能分辨的最小尺度的关系，给出不同差分格式对计算网格步长的限制。同时分析了计算中混淆误差的来源和控制方法，指出了迎风型紧致格式能很好地控制混淆误差。2．将上述格式与三阶精度的Adams半隐格式相结合，构造了不可压槽道湍流直接数值模拟的高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项，避免了压力边界条件处理的困难。利用FFT对方程中的隐式部分进行解耦，解耦后的方程采用追赶法（LU分解法）求解，大大减少了计算量。为了检验该方法，进行了三维不可压槽道湍流的直接数值模拟，得到了Re＝2800的充分发展不可压槽道湍流，并对该湍流场进行了统计分析。包括脉动速度偏斜因子在内的各阶统计量与实验结果及Kim等人的计算结果吻合十分理想，说明本方法是行之有效的。3．进行了三维充分发展的可压缩槽道湍流的直接数值模拟。得到了 Re=3300，Ma=0.8的充分发展可压槽道湍流的数据库。流场的统计特征（如等效平均速度分布，“半局部”尺度无量纲化的脉动速度均方根）和他人的数值计算结果吻合。得到了可压槽道湍流的各阶统计量，其中脉动速度的偏斜因子和平坦因子等高阶统计量尚未见其他文献报道。同时还分析了压缩性效应对壁湍流影响的机理，指出近壁处的压力－膨胀项将部分湍流脉动的动能转换成内能，使得可压湍流近壁速度条带结构更加平整。4．模拟了二维不可压槽道流动的饱和态（所谓“二维湍流”），分析了“二维槽道湍流”的非线性行为特征。分析了流场中的上抛－下扫和间歇现象，研究了“二维湍流”与三维湍流的区别。指出“二维湍流”反映了三维湍流的部分特征，同时指出了展向扰动对于湍流核心区发展的重要性。5．首次对可压缩槽道湍流及“二维槽道湍流”标度律进行了分析，得出了以下结论：a）槽道湍流中，在槽道中心线附近较宽的区域，存在标度律。b）该区域流场存在扩展自相似性（ESS）。c）在Mach数不是很高时，压缩性对标度指数影响不大。本文结果同SL标度律的理论值吻合较好，有效支持了该理论。对“二维槽道湍流”也有相似的结论，但与三维湍流不同的是，“二维槽道湍流”存在标度律的区域更宽，近壁处的标度指数比中心处有所升高。