Dynamics of Incoherent Photovoltaic Spatial Solitons

Propagation properties of bright and dark incoherent beams are numerically studied in photovoltaic-photorefractive crystal by using coherent density approach for the first time. Numerical simulations not only exhibit that bright incoherent photovoltaic quasi-soliton, grey-like incoherent photovoltaic soliton, incoherent soliton doublet and triplet can be established under proper conditions, but also display that the spatial coherence properties of these incoherent beams can be significantly affected during propagation by the photovoltaic field.

多相聚合物体系相分离动力学的理论和模拟

本文提出了将自洽场理论与多相流格子Boltzmann方法相结合的模型，并从该模型出发推导出了描述动量守恒的Navier-Stokes方程和扩散方程，验证了模型理论上的正确性。应用此模型，对聚合物的相分离过程进行了模拟。 首先证实了本模型最终能够得到正确的热力学平衡结果。对于高分子共混物和嵌段共聚物相分离的动力学过程。在分相各个阶段，对高分子共混物和嵌段共聚物都进行了验证。 其次，应用所提出的格子Boltzmann模型，分别对二元聚合物共混物和二嵌段共聚物的相分离后期相区增长过程进行了研究。 最后，通过模型的进一步扩展，实现了对具有复杂结构的嵌段共聚物和复杂共混物的模拟。

Potential light source for photonic integrated circuits based on the equilateral triangle resonator

The semiconductor microlasers based on the equilateral triangle resonator (ETR) can be fabricated from the edge-emitting laser wafer by dry-etching technique, and the directional emission can be obtained by connecting an output waveguide to one of the vertices of the ETR. We investigate the mode characteristics, especially the mode quality factor, for the ETR with imperfect vertices, which is inevitable in the real technique process. The numerical simulations show that the confined modes can still have a high quality factor in the ETR with imperfect vertices. We can expect that the microlasers is a suitable light source for photonic integrated circuits.

Thrust generation and wake structure of wiggling hydrofoil

Marine animals and micro-machines often use wiggling motion to generate thrust. The wiggling motion can be modeled by a progressive wave where its wavelength describes the flexibility of wiggling animals. In the present study, an immersed boundary method is used to simulate the flows around the wiggling hydrofoil NACA 65-010 at low Reynolds numbers. One can find from the numerical simulations that the thrust generation is largely determined by the wavelength. The thrust coefficients decrease with the increasing wavelength while the propulsive efficiency reaches a maximum at a certain wavelength due to the viscous effects. The thrust generation is associated with two different flow patterns in the wake: the well-known reversed Karman vortex streets and the vortex dipoles. Both are jet-type flows where the thrust coefficients associated with the reversed Karman vortex streets are larger than the ones associated with the vortex diploes.

Theoretical analysis of the relationships between hardness, elastic modulus, and the work of indentation for work-hardening materials

In our previous paper, the expanding cavity model (ECM) and Lame solution were used to obtain an analytical expression for the scale ratio between hardness (H) to reduced modulus (E-r) and unloading work (W-u) to total work (W-t) of indentation for elastic-perfectly plastic materials. In this paper, the more general work-hardening (linear and power-law) materials are studied. Our previous conclusions that this ratio depends mainly on the conical angle of indenter, holds not only for elastic perfectly-plastic materials, but also for work-hardening materials. These results were also verified by numerical simulations.

The main factor influencing the tensile properties of surface nano-crystallized graded materials

The improved mechanical properties of surface nano-crystallized graded materials produced by surface severe plastic deformation ((SPD)-P-2) are generally owing to the effects of the refined structure, work-hardened region and compressive residual stress. However, during the (SPD)-P-2 process, residual stress is produced simultaneously with work-hardened region, the individual contribution of these two factors to the improved mechanical properties remains unclear. Numerical simulations are carried out in order to answer this question. It is found that work hardening predominates in improving the yield strength and the ultimate tensile strength of the surface nano-crystallized graded materials, while the influence of the residual stress mainly emerges at the initial stage of deformation and decreases the apparent elastic modulus of the surface nano-crystallized graded materials, which agrees well with the experimental results. (C) 2010 Elsevier B.V. All rights reserved.

The Effects of Roughness on Adhesion Hysteresis

Adhesion hysteresis is defined as the difference between the work needed to separate two surfaces and that originally gained on bringing them together. Adhesion hysteresis is a common phenomenon in most surface/interface interactions. This paper studies the effects of surface roughness on adhesion hysteresis. We assumed that the surface asperity height distribution is Gaussian. Numerical simulations based on Fuller's model showed that adhesion hysteresis depended upon a single dimensionless parameter, the adhesion parameter, which represents the statistical average of a competition between the compressive forces exerted by the higher asperities, which are trying to separate the surfaces, and the adhesion forces of the lower asperities which are trying to hold the surfaces together. (C) Koninklijke Brill NV, Leiden, 2010

Mechanism Responsible for the Complete Suppression of Karman Vortex in Flows Past a Wavy Square-Section Cylinder

The Karman vortex shedding is totally suppressed in flows past a wavy square-section cylinder at a Reynolds number of 100 and the wave steepness of 0.025. Such a phenomenon is illuminated by the numerical simulations. In the present study, the mechanism responsible for it is mainly attributed to the vertical vorticity. The geometric disturbance on the rear surface leads to the appearance of spanwise flow near the base. The specific vertical vorticity is generated on the rear surface and convecting into the near wake. The wake flow is recirculated with the appearance of the pair of recirculating cells. The interaction between the upper and lower shear layers is weakened by such cells, so that the vortex rolls could not be formed and the near wake flow becomes stable.

Thermocapillary migrations of drops and bubbles

The transport phenomenon of drops or bubbles is a very important topic in fundamental hydrodynamics research and practical applications such as material processing and the chemical engineering. In microgravity environment, if drops or bubbles stay in a continuous phase with non-uniform temperature ¯eld, they will start to move as a result of the variance of the interface tension. This kind of movement is called the Marangoni migration. This review tries to sum up the main results in this ¯eld on theoretical analysis, numerical simulations and experiments. So far the theoretical analysis is still limited to the linear or weak nonlinear steady questions, while the current numerical simulations can already obtain the time- dependent process of the bubble/drop migration when the e®ect of heat convection is small. For strong heat convection problem, or when the Marangoni number is bigger than 100, no numerical result is in consistence with those of experiments so far. Some of the lastest numerical results are shown when heat convection is strong, and the main di®erence between strong and weak heat convection is analyzed. Finally, we also discuss the main unresolved problems in this ¯eld and some possible directions in the future.

MHD控制超声速边界层的理论研究和数值分析

对MHD(mechanisms of magnetohy drodynamics)控制超声速平板湍流边界层的机理进行了理论研究和数值模拟.理论上,采用等离子体低频近似碰撞频率模型,建立等离子体中电子和离子的力平衡方程,得到等离子体速度、极化电场以及边界层速度.数值上,通过空间HLLE格式、LU-SGS时间推进求解时均磁流体动力学湍流方程,其中湍流模型采用sst-kω双方程模型.研究结果表明:(1)边界层速度的理论结果和数值结果误差在7%范围内;(2)只有磁场而电场为零时,洛仑兹力起到减小摩阻的作用.施加电场后,洛仑兹力能够加速边界层低速区流体;(3)在边界层外层,越靠近壁面,作用参数越小;而在边界层近壁区黏性底层,虽然惯性力减小,但黏性力却迅速增加,因此越靠近壁面,作用参数反而越大,加速低速流的代价增加.

Heat Transfer Problems Induced By Multi-shocks Interaction

Numerical simulations of the multi-shock interactions observable around hypersonic vehicles were carried out by solving Navier-Stokes equations with the AUSMPW scheme and the new type of the IV interaction created by two incident shock waves was investigated in detail. Numerical results show that the intersection point of the second incident shock with the bow shock plays important role on the flow pattern, peak pressures and heat fluxes. In the case of two incident shocks interacting with the bow shock at the same position, the much higher peak pressure and more severe heat transfer rate are induced than the classical IV interaction. The phenomenon is referred to as the multi-shock interaction and higher requirements will be imposed on thermal protection systems.

Extinction of lean near-limit methane/air flames at elevated pressures under normal- and reduced-gravity

perimentally at evaluated pressures and under normal- and micro-gravity conditions utilizing the 3.5 s drop tower of the National Microgravity Laboratory of China. The results showed that under micro-gravity conditions the natural convection is minimized and the flames become more planar and symmetric compared to normal gravity. In both normal- and micro-gravity experiments and for a given strain rate and fuel concentration, the flame luminosity was found to enhance as the pressure increases. On the other hand, at a given pressure, the flame luminosity was determined to weaken as the strain rate decreases. At a given strain rate, the fuel concentration at extinction was found to vary non-monotonically with pressure, namely it first increases and subsequently decreases with pressure. The limit fuel concentration peaks around 3 and 4 atm under normal- and micro-gravity, respectively. The extinction limits measured at micro-gravity were in good agreement with predictions obtained through detailed numerical simulations but they are notably lower compared to the data obtained under normal gravity. The simulations confirmed the non-monotonic variation of flammability limits with pressure, in agreement with previous studies. Sensitivity analysis showed that for pressures between one and 5 atm, the near-limit flame response is dominated by the competition between the main branching, H + O2 ? OH + O, and the pressure sensitive termination, H+O2+M? HO2 + M, reaction. However, for pressures greater than 5 atm it was determined that the HO2 kinetics result in further chain branching in a way that is analogous to the third explosion limit of H2/O2 mixtures. 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Bubble Dynamics and Heat Transfer in Microgravity Pool Boiling

Boiling is an extremely complicated and illusive process. Microgravity experiments offer a unique opportunity to study the complex interactions without external forces, such as buoyancy, which can affect the bubble dynamics and the related heat transfer. Furthermore, they can also provide a means to study the actual influence of gravity on the boiling. Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments both in normal gravity and in short-term microgravity in the Drop Tower Beijing and numerical simulations have also been performed. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. It was found that the bubble dynamics in microgravity has a distinct difference from that in normal gravity, and that the heat transfer characteristic is depended upon the bubble dynamics. Lateral motions of bubbles on the heaters were observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drove it to detach from the heaters. Slight enhancement of heat transfer on wires is observed in microgravity, while diminution is evident for high heat flux in the plate case.

Thermal analysis of microstructural evolution of Al2O3 surface modified layers

Solidification behavior and microstructural evolution of surface modified layers in plasma cladding technique are studied via numerical simulations. Both the coupling effect of temperature and solid volume fraction are considered in the proposed thermal analytical model, by which the transient temperature distributions are calculated and the shape of melting pool is determined. Furthermore, we perform microscopic thermal analysis on the nucleation and growth behaviors of ceramic hardening phases and dendrites, as well as the kinetics of related two-phase flow systems. By comparing with experimental observations, the evolution mechanisms of the morphology of Al2O3 ceramic hardening layer are explained. Based on the above results, a relationship among the scanning velocity of plasma stream, dendritic growth rate and the advancing speed of solid/liquid interface is found, and an energy criterion is proposed for predicting the pushing/engulfing transition of ceramic particles by grain growth fronts. (C) 2009 Elsevier B.V. All rights reserved.

Controlling parameter for wave types of long flexible riser undergoing vortex-induced vibration

Submerged floating tunnel (SFT) is a popular concept of crossing waterways. The failure of the cable may occur due to vortex-induced-vibration (VIV), and the stability of the cable is crucial to the safety of the entire tunnel. Investigation results in recent years show that the vortex-induced vibration of the flexible cables with large aspect ratio reveals some new phenomena, for example, the vortex-induced wave, multi-mode competition, wide band random vibration, which have brought new challenges to the study of vortex-induced vibration of long flexible cables. In this paper, the dimensionless parameter controlling the wave types of dynamic response of slender cables undergoing vortex-induced vibration is investigated by means of dimensional analysis and finite element numerical simulations. Our results indicate that there are three types of response for a slender cable, i.e. standing wave vibration, traveling wave vibration and intermediate state. Based on dimensional analysis the controlling parameter is found to be related to the system damping including fluid damping and structural damping, order number of the locked-in modes and the aspect ratio of cable. Furthermore through numerical simulations and parameter regression, the expression and the critical value of controlling parameter is presented. At last the physical meaning of the parameter is analyzed and discussed.