纳米压痕过程的三维有限元数值试验研究

采用有限元方法模拟了纳米压痕仪的加、卸载过程，三维有限元模型考虑了纳米压痕仪的标准Berkovich压头．介绍了有限元模型的几何参数、边界条件、材料特性与加载方式，讨论了摩擦、滑动机制、试件模型的大小对计算结果的影响，进行了计算结果与标准试样实验结果的比较，证实了模拟的可靠性．在此基础上，重点研究了压头尖端曲率半径对纳米压痕实验数据的影响．对比分析了尖端曲率半径r=O与r=100nm两种压头的材料压痕载荷-位移曲线．结果表明，当压头尖端曲率半径r≠O时，基于经典的均匀连续介质力学本构理论、传统的实验手段与数据处理方法，压痕硬度值会随着压痕深度的减小而升高．

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Diffusive waves in a channel with concentrated lateral inflow

In the current paper an analytical solution for diffusive wave equation with the concentrate-distributed lateral inflow is yielded. Finite-difference numerical method is also employed to validate this model. The backwater effects drawn from lateral inflow on the mainstream are examined finally.

微尺度压痕实验的数值模拟及相关问题研究

目前国际上占主导地位的纳米压痕技术是由Oliver与Pharr提出并发展,目前的纳米压痕可以给出整个加、卸载过程的载荷—位移曲线以及硬度与弹性模量随压痕深度变化的曲线,从而提供了丰富的、比较精确的信息,为利用它探索材料比较完整的力学特性提供了可能.为达到该目的,就必须对压痕实验的加、卸过程进行较为深入的研究.作为主要的研究工具,有限元方法模拟微压痕过程在探讨通过实验数据得到更多、更准确的材料表层力学性能参数以及解释实验现象等方面发挥着重要作用.基于计算机速度与容量的原因,较早进行微压痕过程有限元模拟的BhattacharyaandNix、LaursenandSino都使用圆锥压头模拟维氏显微硬度标准正四棱锥Vicker压头与纳米压痕仪标准正三棱锥Berkovich压头,因为圆锥压头具有旋转对称性,可用二维旋转对称单元(二维实体单元)进行计算从而降低计算规模.即便如此,以当时大型计算机的水平,对规模为400~2000个四节点矩形单元的有限元模型进行一次完整的加、卸载过程也需要1~2天.到目前为止,微尺度压痕实验的数值模拟沿用二维模型.事实上,由于加工工艺的限制,微尺度压痕仪的压头如Berkovich与Vicker压头均不个旋转对称性;就微观尺度而言,实际的表层材料都是非均匀的.这些特征均不能由二维模拟体现,所以该文首先建立三维有限元模型,模拟带滑动接触的微尺度压痕加、卸载过程.在此基础上重点讨论了压头几何效应的问题,如二维模拟与三维模拟的关系、显微硬度与纳米的压痕硬度的关系、不同压头下材料的应力应变场、压痕间距与压痕边界的效应等,最后针对微尺度压痕实验中出现的压痕硬度随压痕深度减小而升高的现象,讨论了影响不同压痕深度硬度值的因素.

若干流体力学问题的数学性质研究

该文利用高智的扩散抛物化方程组理论及流体力学基本方程组的特征次特征理论,流体大小尺度（LSS）方程组理论以及摄动有限差分（PFD）方法,研究若干流体力学问题的数学性质.该文得到的主要结论有：1.利用湍流大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组,并证明两个关于湍流大小尺度涡量的命题,从而得到湍流封闭大小尺度涡量（CLSSV）方程组,并对已有的近程相互作用命题进行推广.2.根据扩散抛物化方程组理论和流体力学层次结构方程组的特征和次特征方法,研究了抛物化稳定性方程组（PSE）的特征和次特征以及消除PSE的剩余椭圆特性的问题.3.利用摄动有限差分（PFD）方法得到对流扩散反应方程的变步长摄动有限差分格式,是等步长摄动有限差分格式的推广.

Propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems

We investigate the propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems using an iterative predictor-corrector finite-difference time-domain method. It is shown that when the initial effective area is equal to 2 pi, the effective area will remain invariant during the course of propagation, and a complete Rabi oscillation can be achieved. However, for an elliptically polarized few-cycle ultrashort laser pulse, polarization conversion can occur. Eventually, the laser pulse will evolve into two separate circularly polarized laser pulses with opposite helicities.

Propagation of a few-cycle laser pulse in a V-type three-level system

Propagation of a few-cycle laser pulse in a V-type three-level system (fine structure levels of rubidium) is investigated numerically. The full three-level Maxwell-Bloch equations without the rotating wave approximation and the standing slowly varying envelope approximation are solved by using a finite-difference time-domain method. It is shown that, when the usual unequal oscillator strengths are considered, self-induced transparency cannot be recovered and higher spectral components can be produced even for small-area pulses. (c) 2005 Pleiades Publishing, Inc.

Polarization evolution of a few-cycle ultrashort laser pulse propagating in a degenerate three-level medium

The propagation of an arbitrary polarized few-cycle ultrashort laser pulse in a degenerate three-level medium is investigated by using an iterative predictor-corrector finite-difference time-domain method. It is found that the polarization evolution of the ultrashort laser pulse is dependent not only on the initial atomic coherence of the medium but also on the polarization condition of the incident laser pulse. When the initial effective area is equal to 2 pi, complete linear-to-circular and circular-to-linear polarization conversion of few-cycle ultrashort laser pulses can be achieved due to the quantum interference effects between the two different transition paths.

Polarization-dependent Talbot effect

The term "polarization-dependent Talbot effect" means that the Talbot self-imaging intensity of a high-density grating is different for TE and TM polarization modes. Numerical simulations with the finite-difference time-domain method show that the polarization dependence of the Talbot images is obvious for gratings with period d between 2 lambda and 3 lambda. Such a polarization-dependent difference for TE and TM polarization of, a high-density grating of 630 lines/mm (corresponding to d/lambda = 2.5) is verified through experiments with the scanning near-field optical microscopy technique, in which a He-Ne laser is used as its polarization is changed from the TE mode to the TM mode. The polarization-dependent Talbot effect should help us to understand more clearly the diffraction behavior of a high-density grating in nano-optics and contribute to wide application of the Talbot effect. (c) 2006 Optical Society of America.

Talbot effect of a grating with different kinds of flaws

The Talbot effect of a grating with different kinds of flaws is analyzed with the finite-difference time-domain (FDTD) method. The FDTD method can show the exact near-field distribution of different flaws in a high-density grating, which is impossible to obtain with the conventional Fourier transform method. The numerical results indicate that if a grating is perfect, its Talbot imaging should also be perfect; if the grating is distorted, its Talbot imaging will also be distorted. Furthermore, we evaluate high-density gratings by detecting the near-field distribution with the scanning near-field optical microscopy technique. Experimental results are also given. (c) 2005 Optical Society of America.

FDTD analysis of Talbot effect of a high density grating

Talbot effect of a grating with different flaws is analyzed with the finite-difference time-domain (FDTD) method. The FDTD method can show the exact near-field distribution of different flaws in a high-density grating, which is impossible to obtain with the conventional Fourier transform method. The numerical results indicate that if a grating is perfect, its Talbot imaging should also be perfect; if the grating is distorted, its Talbot imaging would also be distorted. Furthermore, we can evaluate high density gratings by detecting the near-field distribution.

Phase interpretation for polarization-dependent near-field images of high-density gratings

It has been described that the near-field images of a high-density grating at the half self-imaging distance could be different for TE and TM polarization states. We propose that the phases of the diffraction orders play an important role in such polarization dependence. The view is verified through the coincidence of the numerical result of finite-difference time-domain method and the reconstructed results from the rigorous coupled-wave analysis. Field distributions of TE and TM polarizations are given numerically for a grating with period d = 2.3 lambda, which are verified through experiments with the scanning near-field optical microscopy technique. The concept of phase interpretation not only explains the polarization dependence at the half self-imaging distance of gratings with a physical view, but also, it could be widely used to describe the near-field diffraction of a variety of periodic diffractive optical elements whose feature size comparable to the wavelength. (C) 2008 Elsevier B.V. All rights reserved.

Absolute left-handed behaviors in a triangular elliptical-rod photonic crystal

In this paper we theoretically study the left-handed behaviors in a two-dimensional triangular photonic crystal made of elliptical rods in air. An absolute left-handed region is found in the second photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Typical left-handed behaviors such as negative refraction, flat superlensing and plano-concave lensing are demonstrated by the finite-difference time-domain simulations. These behaviors are also compared with the quasi-negative refraction and the resulted focusing effects in a square-lattice two-dimensional photonic crystal. (c) 2005 Optical Society of America

Subwavelength imaging with two symmetrical interfaces by dielectric-tube photonic crystals

As distinct from coated photonic crystals, in this paper we propose a novel one that is made of dielectric tubes arranged in a close-packet square lattice. Without metallic cores, this structure is low-loss and convenient to fabricate. A left-handed frequency region is found in the second band by dispersion characteristic analysis. Without inactive modes for the transverse electric mode, negative refraction and subwavelength imaging are demonstrated by the finite-difference time-domain simulations with two symmetrical interfaces, i.e. Gamma X and Gamma M.

Rigorous transmittance analysis of a sub-wavelength Sb thin film lens

We quantitatively analysed the factors contributing to the optical transmission enhancement of a sub-wavelength Sb thin film lens, using the finite-difference time-domain (FDTD) method. The results show that the transmission enhancement of the dielectric with a Gaussian distributed refractive index loaded in a sub-wavelength circular hole is not only due to the high refractive index dielectric, but also due to the specific distributions of refractive index. It is the first study about the effects of the refractive index distribution on the transmission of a sub-wavelength aperture. This kind of lens has practical applications in the very small aperture lasers and for near-field optical storage and lithography.