Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非特殊振荡的直接原因是在于非均一的波群的群速度，因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式，以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例，证明了文中所构造的格式具有方法简便，物理含义清楚，精度高，网格基架点小，和捕捉激波能力较强的优点。

The Heterogeneous Growth of P(3MeTh)--An Ellipsometric Study

Successive thicker P(3MeTh) layers are analysed by ex situ conventional and imaging ellipsometry. Thin films display a smooth surface, are compact and homogeneous while for a growth charge above 20 mC cm(-2) the polymer structure modifies to a still uniform but less dense layer. A two-layer model is used and a mathematical procedure is developed to obtain, simultaneously, from the experimental ellipsometric parameters, Delta and Psi, the thickness and the complex refractive index of P(3MeTh) films grown up to 80 mC cm(-2). Thicker polymer layers are disordered and present a high degree of surface roughness. 

迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jaeobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer-Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.

Crack Tip Field and J-Integral with Strain Gradient Effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventional J_2 deformation theory No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory Two typical crack Problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-I K-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.

Explosive consolidation of amorphous cobalt-based alloys

This paper addresses the explosive consolidation of amorphous cobalt-based alloys. Using the experimental setup introduced in the present paper, specimens with high compact density, excellent magnetic properties and great wearability have been made. In comparison with permalloy and ferrite, the present specimens exhibit superior magnetic properties. Therefore, the compact is deemed as being a promising material for magnetic recording heads.

Growth Mechanism And Corrosion Behavior Of Ceramic Coatings On Aluminum Produced By Autocontrol Ac Pulse Peo

Ceramic coatings are produced on aluminum alloy by autocontrol AC pulse Plasma Electrolytic Oxidation (PEO) with stabilized average current. Transient signal gathering system is used to study the current, voltage, and the transient wave during the PEO process. SEM, OM, XRD and EDS are used to study the coatings evolution of morphologies, composition and structure. TEM is used to study the micro profile of the outer looser layer and inner compact layer. Polarization test is used to study the corrosion property of PEO coatings in NaCl solution. According to the test results, AC pulse PEO process can be divided into four stages with different aspects of discharge phenomena, voltage and current. The growth mechanism of AC PEO coating is characterized as anodic reaction and discharge sintering effect. PEO coating can increase the corrosion resistance of aluminum alloy by one order or two; however, too long process time is not necessarily needed to increase the corrosion resistance. In condition of this paper, PEO coating at 60 min is the most protective coating for aluminum alloy substrate. (C) 2008 Elsevier B.V. All fights reserved.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

Numerical analysis of spatial evolution of the small signal gain in a chemical oxygen-iodine laser operating without primary buffer gas

A chemical oxygen iodine laser (COIL) that operates without primary buffer gas has become a new way of facilitating the compact integration of laser systems. To clarify the properties of spatial gain distribution, three-dimensional (3-D) computational fluid dynamics (CFD) technology was used to study the mixing and reactive flow in a COIL nozzle with an interleaving jet configuration in the supersonic section. The results show that the molecular iodine fraction in the secondary flow has a notable effect on the spatial distribution of the small signal gain. The rich iodine condition produces some negative gain regions along the jet trajectory, while the lean iodine condition slows down the development of the gain in the streamwise direction. It is also found that the new configuration of an interleaving jet helps form a reasonable gain field under appropriate operation conditions. (c) 2007 Elsevier Ltd. All rights reserved.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

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.