高频感应等离子体风洞的流场研究

高频感应等离子体风洞适合于基础研究特别是用于防热材料的催化特性的研究，本文用实验测量和数值模拟了该风洞的温度场和速度场。温度场是用OSA （Optical Spectrum Analyzer）进行光谱分析得到的，同时用水冷总压和静压探针对风洞的流场进行测量得到了速度场的空间分布。数值模拟采用完全二维的高频感应等离子模型，我们通过求解相关的能量方程、流体动力学方程和电磁场方程得到了相应的温度场和流场。比较了数值模拟结果与实现测量结果，并给出了不同运行参数下的典型数值结果。

基于统计光学的无透镜鬼成像数值模拟与实验验证

作为量子信息领域分支的鬼成像,由于物体的像将出现在不包含物体的光路上的特点,使得这一领域的研究引人入胜。一度认为,只有基于纠缠态双光子的纠缠光源,才能实现鬼成像;但近年来的研究表明,经典热光场也能实现这一过程。从经典统计光学入手,建立了热光场的数值模型,模拟符合热光特性的光场变化、光场传播、以及物体透射函数对热光场的调制,进而从光强度起伏的关联函数中,分别重现振幅型物体和纯相位型物体的傅里叶变换图像;通过与真实实验结果的对比,表明基于统计光学原理的该数值模型所预测的实验结果,与真实的实验结果完全一致。

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

Influence of Coriolis force on bulk flows during horizontal Bridgman growth on a centrifuge: a numerical study

Three-dimensional and time-dependent numerical simulations are performed For melt convection in horizontal Bridgman crystal growth tinder high gravity conditions by means of a centrifuge. The numerical results show that Coriolis Force can cause a stabilizing effect on the fluctuations of the melt flow under a specific relation direction and relation rates of the centrifuge as reported in previous experiments (Ma et al., Materials Processing in High Gravity, Plenum Press, New York, 1994, p. 61). The present simulation provides details of the now features associated with the effect of the Coriolis force. There are also some differences between the present three-dimensional and former two-dimensional numerical solutions particularly in the prediction of the critical conditions and flow patterns.

On the influence zone and the prediction of tensile strength of particulate polymeric composites

An intended numerical investigation is carried out. The results indicate that, even if a perfect adhesive bond is preserved between the particles and matrix materials, the two-phase element cell model is unable to predict the strength increment of the particulate polymeric composites (PPC). To explore the main reinforcing mechanism, additional microscopic experiment is performed. An ''influence zone'' was observed around each particle which is measured about 2 to 10 micrometers in thickness for a glass-polyethylene mixture. Then, an improved computational model is presented to include the ''influence zone'' effect and several mechanical behaviors of PPC are well simulated through this new model.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Prediction and suppression of strong dispersive coupling in microracetrack channel drop filters

For a four-port microracetrack channel drop filter, unexpected transmission characteristics due to strong dispersive coupling are demonstrated by the light tunneling between the input-output waveguides and the resonator, where a large dropping transmission at off-resonance wavelengths is observed by finite-difference time-domain simulation. It causes a severe decline of the extinction ratio and finesse. An appropriate decrease of the coupling strength is found to suppress the dispersive coupling and greately increase the extinction ratio and finesse, a decreased coupling strength can be realized by the application of an asymmetrical coupling waveguide structure. In addition, the profile of the coupling dispersion in the transmission spectra can be predicted based on a coupled mode theory analysis of an equivalent system consisting of two coupling straight waveguides. The effects of structure parameters on the transmission spectra obtained by this method agree well with the numerical results. It is useful to avoid the strong dispersive coupling region in the filter design. (c) 2007 Optical Society of America.

Prediction of simultaneously large and opposite generalized Goos-Hanchen shifts for TE and TM light beams in an asymmetric double-prism configuration

It is predicted that large and opposite generalized Goos-Hanchen (GGH) shifts may occur simultaneously for TE and TM light beams upon reflection from an asymmetric double-prism configuration when the angle of incidence is below but near the critical angle for total reflection, which may lead to interesting applications in optical devices and integrated optics. Numerical simulations show that the magnitude of the GGH shift can be of the order of beam's width.

An eigenfunction expansion-variational method in prediction of the transverse thermal conductivity of fiber reinforced composites considering interfacial characteristics

An eigenfunction expansion-variational method based on a unit cell is developed to deal with the steady-state heat conduction problem of doubly-periodic fiber reinforced composites with interfacial thermal contact resistance or coating. The numerical results show a rapid convergence of the present method. The present solution provides a unified first-order approximation formula of the effective thermal conductivity for different interfacial characteristics and fiber distributions. A comparison with the present high-order results, available experimental data and micromechanical estimations demonstrates that the first-order approximation formula is a good engineering closed-form formula. An engineering equivalent parameter reflecting the overall influence of the thermal conductivities of the matrix and fibers and the interfacial characteristic on the effective thermal conductivity, is found. The equivalent parameter can greatly simplify the complicated relation of the effective thermal conductivity to the internal structure of a composite. (c) 2010 Elsevier Ltd. All rights reserved.

PREDICTION OF VORTEX-INDUCED VIBRATION OF FLEXIBLE RISER USING AN IMPROVED WAKE-OSCILLATOR MODEL

Based on improving the wake-oscillator model, an analytical model for vortex-induced vibration (VIV) of flexible riser under non-uniform current is presented, in which the variation of added mass at lock-in and the nonlinear relationship between amplitude of response and reduced velocity are considered. By means of empirical formula combining iteration computation, the improved analytical model can be conveniently programmed into computer code with simpler and faster computation process than CFD so as to be suitable to application of practical engineering. This model is validated by comparing with experimental result and numerical simulation. Our results show that the improved model can predict VIV response and lock-in region more accurately. At last, illustrative examples are given in which the amplitude of response of flexible riser experiencing VIV under action of non-uniform current is calculated and effects of riser tension and flow distribution along span of riser are explored. It is demonstrated that with the variation of tension and flow distribution, lock-in region of mode behaves in different way, and thus the final response is a synthesis of response of locked modes.

Numerical simulation of mesoscopic mechanical behaviors of gradual multi-fiber-reinforced polymer matrix composites

To simulate the deformation and the fracture of gradual multi-fiber-reinforced matrix composites, a numerical simulation method for the mesoscopic mechanical behaviors was developed on the basis of the finite element and the Monte Carlo methods. The results indicate that the strength of a composite increases if the variability of statistical fiber strengths is decreased.

Prediction of GC retention values under various column temperature conditions from temperature programmed data

A numerical approach has been developed for the correlation of retention limes (total retention lime) with temperature in gas chromatography, which allows the calculation of retention parameters including retention index from data acquired under two or more different temperature program conditions. By using this procedure the optimization of temperature condition can be further achieved, especially when a temperature-programmed run is the most suitable mode in the preliminary development of an analytical method for the analysis of an unknown sample.

Numerical Analysis of a Dual Polarization Mode-Locked Laser with a Quarter Wave Plate

Dynamical behaviors and frequency characteristics of an active mode-locked laser with a quarter wave plate (QWP) are numerically studied by using a set pf vectorial laser equation. Like a polarization self-modulated laser, a frequency shift of half the cavity mode spacing exists between the eigen-modes in the two neutral axes of QWP. Within the active medium, the symmetric gain and cavity structure maintain the pulse's circular polarization with left-hand and right-hand in turn for each round trip. Once the left-hand or right-hand circularly polarized pulse passes through QWP, its polarization is linear and the polarized direction is in one of the directions of i45o with respect to the neutral axes of QWP. The output components in the directions of i45" from the mirror close to QWP are all linearly polarized with a period of twice the round-trip time.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.