On dispersion-controlled principles for non-oscillatory shock-capturing schemes

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time. It is only during the last two decades that extensive studies on the dispersion-controlled dissipative (DCD) schemes were reported. The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations. The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do. Research progresses on the dispersion controlled principles are reviewed in this paper, including the exploration of the role of dispersions in numerical simulations, the development of the dispersion-controlled principles, efforts devoted to high-order dispersion-controlled dissipative schemes, the extension to both the finite volume and the finite element methods, scheme verification and solution validation, and comments on several aspects of the schemes from author's viewpoint.

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.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Fracture analysis for 3-dimensional bodies with a surface crack

This paper deals with fracture analyses in 3-dimensional bodies containing a surface crack. A general solution of stress-strain fields at crack tip is proposed. Based on the stress-strain fields obtained, a high-order 3-dimensional special element is established to calculate the stress intensity factors in a plate with a surface crack. The variation of stress intensity factors with geometric parameters is investigated.

高精度非定常激波装配法

为正确模拟高超声速绕流中,来流小扰动与弓形激波之间的干扰对流动特征的影响,将弓形激波作为动边界,利用非定常特征关系处理激波处的边界条件.应用五阶精度迎风紧致格式和六阶精度的对称格式与三阶精度的R-K方法相结合,建立高精度非定常激波装配方法.采用该方法数值模拟钝锥高超声速定常流场和二维抛物外形高超声速边界层流动的感受性问题,数值模拟来流小扰动与弓形激波干扰激波后非定常扰动流场,研究扰动波进入边界层产生边界层不稳定波的特征.

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

高阶精度CE/SE算法及其应用

对时-空守恒元解元算法(CE/SE)的网格设置做较大改进,提出一种新的六面体解元和元定义;同时在解元中对物理量进行高阶Taylor展开,给出一种在时间和空间上均具有高阶精度CE/SE算法.在此基础上,把新型的高阶精度CE/SE算法推广应用于高速流动捕捉激波间断、气相化学反应流动、计及固体动态效应的流体-弹塑性流动和非稳态多相不可压缩粘性流动中.数值实践表明,提出的新型网格结构上的高阶精度CE/SE算法具有算法简单、计算精度高、计算效率和计算效果好的优点,并大大改进和拓展了CE/SE算法的应用范围.

超短激光脉冲在不同色散参量光子晶体光纤中传输的数值模拟

为了使得数值模拟更为精确, 采用广义非线性薛定谔方程(GNSE)描述超短激光脉冲在光子晶体光纤中的传输演化过程, 并利用二阶分步傅里叶方法通过求解方程, 数值计算了相同脉宽和能量的超短脉冲在不同色散参量的光子晶体光纤中非线性传输和超连续谱的产生。比较了超短脉冲在光纤不同色散区传输时, 高阶色散和非线性效应对超连续谱的产生以及对脉冲波形演化的影响。结果表明, 相对于超短脉冲中心波长位于光子晶体光纤的正常和反常色散区, 可以相应获得短波波段和长波波段的超连续谱输出, 当超短脉冲中心波长位于零色散波长点时, 通

Phase control of higher spectral components in the presence of a static electric field

We investigate the higher spectral component generations driven by a few-cycle laser pulse in a dense medium when a static electric field is present. Our results show that, when assisted by a static electric field, the dependence of the transmitted laser spectrum on the carrier-envelope phase (CEP) is significantly increased. Continuum and distinct peaks can be achieved by controlling the CEP of the few-cycle ultrashort laser pulse. Such a strong variation is due to the fact that the presence of the static electric field modifies the waveform of the combined electric field, which further affects the spectral distribution of the generated higher spectral components.

Quantum path control in few-optical-cycle regime

We theoretically show that selection of a single quantum path in high-order harmonics generation can be realized in a few-optical-cycle regime with two-color schemes. We also demonstrate, in theory as well, the generation of spectrally smooth and ultrabroad extreme ultraviolet supercontinuum in argon gas which can produce single similar to 79 as pulses with currently available ultrafast laser sources. Our finding can be beneficial for generating isolated sub-100 as extreme ultraviolet pulses.

较高阶非线性抽运-探测过程中的空间效应

在飞秒抽运一探测光谱技术中，空间分辨的探测光信号反映了在不同空间位置的材料的非线性效应。当抽运光强度增大时，探测光信号中会出现明显的高阶特别是五阶非线性效应。利用劈裂算子方法直接解决了一维非线性传播方程的问题。在数值模拟中，研究了在不同抽运强度和位置下的抽运一探测过程中的五阶非线性效应。在足够高的抽运场下，探测信号出现清晰的振荡，显示了三阶和五阶非线性效应之间的干涉。当空间位置离抽运场中心足够远时，五阶比三阶非线性效应的衰减快得多，对其物理机制和趋势进行了定性的讨论。

高阶群速度色散引起的高斯超短脉冲宽度的展宽和形变

基于动量守恒和光参变过程中的三波耦合波方程, 和负单轴非线性光学晶体CsLiB6O10的色散方程, 研究了在光参变效应中超短激光脉冲由于群速度色散引起的展宽和形变。数值模拟显示, 在超短脉冲波形为双曲正割形和无啁啾调制时, 高阶群速度色散引起的超短脉冲为50 fs时, 晶体长度为10 mm, 紫外光213 nm作为基波入射时的脉冲展宽是波长为532 nm绿光在同等条件下的1.6倍。脉冲展宽程度与入射波长和晶体长度有关, 波长越短和晶体长度越长则脉冲展宽和波形变化越严重,高阶色散引起的超短高斯脉冲展宽, 将

XPM-induced pulse compression in birefringent dispersion decreasing fiber

The pulse compression induced by cross-phase modulation in birefringent dispersion decreasing fiber is discussed theoretically by solving the coupled Schrodinger equations which include the contribution of the high-order non-linear effects, and third-order dispersion. In particular, it is found that a high quality compressed signal pulse can be obtained by a pump pulse of low intense through the technique. The dependence of optimum compression on the non-linear factor N, time delay tau(d) and the dispersive ratio f is also discussed in detail. (c) 2005 Elsevier B.V. All rights reserved.

In situ surface topography measurement method of granite base in scanning wafer stage with laser interferometer

Topography of a granite surface has an effect on the vertical positioning of a wafer stage in a lithographic tool, when the wafer stage moves on the granite. The inaccurate measurement of the topography results in a bad leveling and focusing performance. In this paper, an in situ method to measure the topography of a granite surface with high accuracy is present. In this method, a high-order polynomial is set up to express the topography of the granite surface. Two double-frequency laser interferometers are used to measure the tilts of the wafer stage in the X- and Y-directions. From the sampling tilts information, the coefficients of the high-order polynomial can be obtained by a special algorithm. Experiment results shows that the measurement reproducibility of the method is better than 10 nm. (c) 2006 Elsevier GmbH. All rights reserved.