Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Coupling of melting-solidification phase-change convection with floating-zone convection under microgravity

A finite element algorithm is used to analyze the process of floating zone crystal growth under microgravity. The effect of phase change convection coupled with surface tension convection is considered. The results show that the rate of crystal growth is very important. The single-crystal-melt interface is steeper than the feed-melt interface during the process of crystal growth. When the rate exceeds a critical value, the Marangoni vortex near the feed-melt interface will become so large that a secondary vortex will exist.

Preliminary-results of gaas single-crystal growth under high gravity conditions

GaAs single crystals have been grown under high gravity conditions, up to 9g0, by a recrystallization method with decreasing temperature. The impurity striations in GaAs grown under high gravity become weak and indistinct with smaller striation spacings. The dislocation density of surcharge-grown GaAs increases with increase of centrifugal force. The cathodoluminescence results also show worse perfection in the GaAs grown at high gravity than at normal earth gravity.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

SiC_p/6151Al复合材料的微结构效应

利用分离式Hopkinson压杆和MTS通用材料试验机研究了SiC_p/6151Al颗粒增强复合材料在不同应变率下的变形行为和增强颗粒的尺寸对复合材料微结构及变形行为的影响。结果表明，对于在不同应变率下的SiC_p/6151Al复合材料，增强颗粒尺寸大小的流动应力高于增强颗粒尺寸的流动应力。根据位错强化理论中的Hall-Petch关系对这个结果进行了解释。首次在实验上观测到增强颗粒对复合材料微损伤-微带形成的影响，并根据微带（microband）形成的双位错墙理论（double dislocation walls）, 分析了增强颗粒对复合材料微带损伤及力学性能影响的微结构效应。

高压扭转铜试样的微观组织与压缩性能

通过高压扭转对铜试样施加不同程度的变形，研究了样品扭转面（ND面）和纵截面（TD面）上微观组织特征．对ND面，在较小的剪应变下，原始晶粒形貌模糊，晶粒内部形成等轴状的位错胞及亚晶结构；随变形量的增大，亚晶间取向差及亚晶内部的位错密度增大，最后形成亚微米尺度的等轴晶粒．对TD面，变形初期原始晶粒被拉长，晶粒内部为位错墙分割成的层状结构，层内为拉长的位错胞；随变形程度的增大，拉长晶粒的宽度减小，与剪切方向的夹角减小，晶内层状组织间距减小，并逐渐演化成拉长的亚晶组织；进一步增大变形，晶粒拉长痕迹消失，变形组织与ND面相似，为等轴状亚微米晶粒．压缩实验表明，经16圈扭转后，整个试样上的压缩性能基本均匀，σ0．2达到385MPa，应变率敏感性指数增大至0．021．

水中悬浮隧道锚索在剪切流中的涡激响应

以Iwan改进的尾流振子模型为基础，给出了剪切流中水中悬浮隧道锚索横向涡激振动的工程分析方法，分析了悬浮隧道重浮比以及剪切流特性对锚索涡激响应的影响．计算结果表明，悬浮隧道重浮比的改变，使锚索各阶模态的频锁区域分布发生变化，从而改变了结构的涡激振动响应；剪切流陡度参数的增大会使锚索涡激响应的幅值减小，若使用均匀流场代替剪切流场计算悬浮隧道锚索的涡激振动，会过高估计其响应幅值．

中等雷诺数圆柱尾流旋涡脱落的控制

介绍通过实验对圆柱尾流旋涡脱落进行抑制的方法及其结果。实验模型的展径比为38，实验的雷诺数范围为3 * 10~2～1.6 * 10~3。抑制方法是在圆柱（直径为D）表面沿展向每隔一定间距伸出一直径0.18D、长度为1.5D的小棒，实验结果表明，当棒间距小于3D，棒与来流夹角在30 °～ 90 °范围内，可有效抑制旋涡脱落。

湍流大小尺度涡量方程组及其封闭

利用流体大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组，给出两个关于湍流大小尺度涡量的命题，从而得到湍流封闭大小尺度涡量（CLSSV）方程组。同时，对近程相互作用命题进行了推广。

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Nanocrystalline Formation during Mechanical Attrition of Cobalt

The nanocrystalline (nc) formation was studied in cobalt (a mixture of c (hexagonal close packed) and gamma (face-centered cubic) phases) subjected to surface mechanical attrition treatment. Electron microscopy revealed the operation of {10(1) over bar 0}< 11(2) over bar 0 > prismatic and {0001}< 11(2) over bar 0 > basal slip in the E phase, leading to the successive subdivision of grains to nanoscale. In particular, the dislocation splitting into the stacking faults was observed to occur in ultrafine and nc grains. By contrast, the planar dislocation arrays, twins and martensites were evidenced in the gamma phase. The strain-induced gamma ->epsilon martensitic transformation was found to progress continuously in ultrafine and nc grains as the strain increased. The nc formation in the gamma phase was interpreted in terms of the martensitic transformation and twinning.

Physical Modeling of Current-Induced Seabed Scour around a Vibrating Submarine Pipeline

Most of the existing researches either focus on vortex-induced-vibrations (VIV) of a pipeline near a rigid boundary, or on seabed scour around a fixed pipeline. In the fields, pipeline vibration and seabed scour are actually always coupled. Based on the similarity analysis, a series of tests were conducted with a hydro-elastic facility to investigate the influence of pipe vibration on the local scour and the effects of scour process on the pipeline dynamic responses. Experimental results indicate that, there exist two phases in the process of sand scouring around the pipeline with small embedment, i.e. Phase I: scour beneath pipe without VIV, and Phase II: scour with VIV of pipe. It is also found that the gap-to-diameter ratio (e/D) has much effect upon the scour depth for the fixed pipes. For a given value of e/D, the vibrating pipes with close proximity to seabed may induce a deeper scour hole than the fixed ones. Within the examined gap-to-diameter ratio range (425 < e/D < 0.75), the influences of gap-to-diameter ratio on the maximum values of scour-depth for the case of vibrating pipes are not as much as those for the case of fixed pipes.

Characteristics of Wake Flow in Long Bluff Bodies

Vortex shedding is the main characteristics of bluff bodies, which will cause the bluff bodies to vibrate and sometimes result in the structures failure. In this paper the wake flow characteristics of 21 bluff bodies with rectangular, rounded and angular profiles and the length-to-width ratio in the range of 4~12 were deeply studied by Micro ADV. Two parameters, namely the relative intensity of the load due to Karman vortices and the large scale vortex intensity, were introduced to measure the wake flow intensity. Generally, the values of these parameters for different bluff bodies are consistent with each other. The experiment results showed that the key factor affecting the wake flow characteristics is the bluff edge, especially the leading edge geometry. The wake flow in bluff bodies with rounded edge profiles has more regular vortices and becomes more periodic than that in bluff bodies with rectangular ones. A bluff body with angular edged profile was witnessed to have not only small wake loading but small hydraulic resistance also.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.