Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Effects of geometric shape on the hydrodynamics of a self-propelled flapping foil

The hydrodynamics of a free flapping foil is studied numerically. The foil undergoes a forced vertical oscillation and is free to move horizontally. The effect of chord-thickness ratio is investigated by varying this parameter while fixing other ones such as the Reynolds number, the density ratio, and the flapping amplitude. Three different flow regimes have been identified when we increase the chord-thickness ratio, i.e., left-right symmetry, back-and-forth chaotic motion, and unidirectional motion with staggered vortex street. It is observed that the chord-thickness ratio can affect the symmetry-breaking bifurcation, the arrangement of vortices in the wake, and the terminal velocity of the foil. The similarity in the symmetry-breaking bifurcation of the present problem to that of a flapping body under constraint is discussed. A comparison between the dynamic behaviors of an elliptic foil and a rectangular foil at various chord-thickness ratios is also presented.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

Laser-field synchrotron radiation

An electron with an appropriate initial velocity injected into an oncoming, ultraintense circularly polarized laser pulse can execute a circular relativistic motion at the peak of the laser pulse. The circulating electron then radiates in the same manner as that in the storage ring of a conventional synchrotron source. Owing to the extremely small orbit radius, the laser-field synchrotron radiation thus generated can be a compact source of radiation pulses at short wavelength and short duration.

New doubly periodic waves of the (2+1)-dimensional double sine-Gordon equation

New exact solutions of the (2 + 1)-dimensional double sine-Gordon equation are studied by introducing the modified mapping relations between the cubic nonlinear Klein-Gordon system and double sine-Gordon equation. Two arbitrary functions are included into the Jacobi elliptic function solutions. New doubly periodic wave solutions are obtained and displayed graphically by proper selections of the arbitrary functions.

空间激光跟瞄中指向驱动电机运动特性研究

设计了一种深空非合作目标的激光扫描、捕获、跟踪地面实验装置,通过模拟深空同轨道运动的两颗卫星跟瞄过程,在理论上计算了跟瞄装置中光束指向驱动电机的最小加速度和其在跟踪过程中的运动特性。理论分析与仿真结果表明,当卫星偏离光斑中心一定距离时,指向驱动电机先加速后减速,补偿这个偏心,重新捕获跟踪卫星;重新捕获到跟踪所需时间受电机加速度和望远镜探测精度以及探测器响应处理时间影响,其中探测器精度对重新捕获到跟踪所需时间影响较大,探测器响应处理时间要减小到最小;为了使从捕获到跟踪过程中卫星始终在扫描光斑范围内,经纬仪驱动电机的最小角加速度为25.5°/s2。

外差干涉仪混频误差分析

共光路外差干涉仪具有很高的分辨率．但因为安装、调试误差会产生非线性误差．影响系统的测量精度。着重分析了在共光路外差干涉仪中由激光光源的椭圆偏振化和沃拉斯顿棱镜的安装方位角误差同时存在的情况下，引起的频率混叠综合误差的大小及变化规律。结果发现其造成的非线性误差可达2．2nm，同时还发现两者造成的误差在某些情况下存在一定程度的相互抵消作用。讨论了提高测量系统精度的有效措施，对正确设计和调试激光外差测试系统、提高测量系统精度具有重要意义。

Electrically controllable RKKY interaction in semiconductor quantum wires

We demonstrate in theory that it is possible to all-electrically manipulate the RKKY interaction in a quasi-one-dimensional electron gas embedded in a semiconductor heterostructure, in the presence of Rashba and Dresselhaus spin-orbit interaction. In an undoped semiconductor quantum wire where intermediate excitations are gapped, the interaction becomes the short-ranged Bloembergen-Rowland superexchange interaction. Owing to the interplay of different types of spin-orbit interaction, the interaction can be controlled to realize various spin models, e.g., isotropic and anisotropic Heisenberg-like models, Ising-like models with additional Dzyaloshinsky-Moriya terms, by tuning the external electric field and designing the crystallographic directions. Such controllable interaction forms a basis for quantum computing with localized spins and quantum matters in spin lattices.

Quantum tunneling through planar p-n junctions in HgTe quantum wells

We demonstrate that a p-n junction created electrically in HgTe quantum wells with inverted band structure exhibits interesting intraband and interband tunneling processes. We find a perfect intraband transmission for electrons injected perpendicularly to the interface of the p-n junction. The opacity and transparency of electrons through the p-n junction can be tuned by changing the incidence angle, the Fermi energy and the strength of the Rashba spin-orbit interaction (RSOI). The occurrence of a conductance plateau due to the formation of topological edge states in a quasi-one-dimensional (Q1D) p-n junction can be switched on and off by tuning the gate voltage. The spin orientation can be substantially rotated when the samples exhibit a moderately strong RSOI.

Electronic structures of wurtzite ZnO, BeO, MgO and p-type doping in Zn1-xYxO (Y = Mg, Be)

Using the density function theory within the generalized gradient approximation, the band structures of wurtzite ZnO, BeO and MgO have been calculated. The effective-mass parameters are fitted using the calculated eigenvalues. The Dresselhaus spin-orbit effect appears in the k[1 00] direction, and is zero in the high symmetry direction k[00 1]. The orderings of valence band split by the crystal-field and spin-orbit coupling in wurtzite ZnO, BeO and MgO are identified by analyzing the wave function characters calculated by projecting the wave functions onto p-state in the spherical harmonics. For wurtzite ZnO, the ordering of valence band is Still Gamma(7) > Gamma(9) > Gamma(7) due to the negative spin-orbit coupling splitting energy and the positive crystal-field splitting energy. Thus, the Thomas' conclusion is confirmed. For wurtzite BeO and MgO, although their orderings of valence bands are Gamma(7) > Gamma(9) > Gamma(7) too, the origins of their orderings are different from that of wurtzite ZnO. Zn1-x,YxO (Y = Mg, Be) doped with N and P atoms have been studied using first-principles method. The calculated results show that N atom doped in Zn1-x BexO has more shallow acceptor energy level with increasing the concentration of Be atom. (C) 2008 Elsevier B.V. All rights reserved.

Uniaxial Strain Effects on Optical Properties of c-plane Wurtzite GaN

We investigate the uniaxial strain effect in the c-plane on optical properties of wurtzite GaN based on k center dot p theory, the spin-orbit interactions are also taken into account. The energy dispersions show that the uniaxial strain in the c-plane gives an anisotropic energy splitting in the k(x) - k(y) plane, which can reduce the density of states. The uniaxial strain also results in giant in-plane optical polarization anisotropy, hence causes the threshold carrier density reduced. We clarify the relations between the uniaxial strain and the optical polarization properties. As a result, it is suggested that the compressive uniaxial strain perpendicular to the laser cavity direction in the c-plane is one of the preferable approaches for the effcient improvement of GaN-based laser performance.

Spin-dependent transport and spin polarization in coupled quantum wells

We theoretically investigate the electron transport and spin polarization of two coupled quantum wells with Dresselhaus spin-orbit interaction. In analogy with the optical dual-channel directional coupler, the resonant tunneling effect is treated by the coupled-mode equations. We demonstrate that spin-up and -down electrons can be completely separated from each other for the system with an appropriate system geometry and a controllable barrier. Our result provides a new approach to construct spin-switching devices without containing any magnetic materials or applying a magnetic field. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2981204]