Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Effects of adiabatic exponent on Richtmyer-Meshkov instability

We present a systematical numerical study of the effects of adiabatic exponent gamma on Richtmyer-Meshkov instability (RMI) driven by cylindrical shock waves, based on the gamma model for the multi-component problems and numerical simulation with high-order and high-resolution method for compressible Euler equations. The results show that the RMI of different gamma across the interface exhibits different evolution features with the case of single gamma. Moreover, the large gamma can hold back the development of nonlinear structures, such as spikes and bubbles.

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算法的应用范围.

7 fs超快强激光驱动Ar原子产生支持单个阿秒脉冲的高次谐波连续谱

采用脉冲宽度为7 fs，脉冲能量为0．4mJ的超快强激光脉冲与气体盒子中Ar原子作用获得了高次谐波截止区连续谱，并发现当驱动激光稳定在不同的载波包络相位时，高次谐波的谱结构、谱调制深度和连续谱的带宽都有很大区别。在某些载波包络相位时获得了平滑的连续谱，调制深度小于17％，连续带宽达10eV，从而支持时域上获得变换极限500as的单个阿秒脉冲。

Single attosecond pulse generation using two-color polarized time-gating technique

It is proposed that single attosecond pulses be generated via high-order harmonic generation by using a two-color pump pulse with time dependent ellipticity. The two-color pump pulse is created by the fundamental field and its second harmonic: the fundamental field is left-circularly polarized and the second harmonic is right-circularly polarized. Numerical simulations show that single attosecond pulses can be produced in the cut-off region by using the synthesis of 20 fs left-hand and right-hand circularly polarized pulses with a pulse delay of 20 fs. The attosecond pulses produced this way are much stronger than that produced by a few-cycle linear polarized pulse of comparable intensity. (c) 2005 Optical Society of America

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

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

极紫外阿秒脉冲在高次谐波产生过程中引起的非偶极效应

当红外强激光和极紫外（XUV）阿秒脉冲共同作用于原子分子时，电离出去的电子通常会吸收和辐射激光光子而发生能量扩展．讨论了由于XUV阿秒脉冲的短波长与扩展后的电子波包尺度可相比拟时在高次谐波产生过程中引起的非偶极效应．采用彤作为模型分子，并把分子轴置于激光场的传播方向，通过解二维含时薛定谔方程并比较考虑非偶极效应和采用偶极近似两种方法计算得到的结果，两者相比，前者的谐波强度降低，谐波频率向低级次稍有移动，电子能谱的能带内出现了更多的光电子峰．在相同的光电子能量处，两种方法计算得到的信号强度相差2—5倍．并且

Control of single attosecond pulse generation from the reflection of a synthesized relativistic laser pulse on a solid surface

The influence of the carrier-envelope phase (CEP) of the driving laser pulse on the generation of single attosecond (as) pulses from surface harmonics by using the polarization gating technique is investigated in detail. It is found that the modulation depth of the high-order harmonic spectrum depends on the CEP, and a strong single 68 as pulse can be generated when the CEP is stable and has the proper value. The physical origin of the influence of the CEP is explained in terms of the oscillating mirror model. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2997342]