Finite element modal solutions of planar photonic crystal fibers with rectangular air-holes

Modal solutions of planar photonic crystal waveguides with rectangular air-holes are presented by using a rigorous full-vectorial finite element-based approach. The effective indices, mode field profiles, spot-sizes, power confinements, modal hybridness, beat lengths and group velocity dispersions are shown for the fundamental and higher order modes of the quasi-TE and TM polarizations.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of x-ray fluorescence analysis to investigation of the composition of gunshot residues

We have used x-ray fluorescence analysis to study the elemental composition of gunshot residues from smooth-bore and rifled-bore guns. We have established that it is possible to differentiate between types of projectiles (jacketed/lead), types of primers (corrosive/noncorrosive), and types of propellant powders (black/smokeless) by analysis of the elemental composition of the gunshot residues. We have shown that the mercury content in gunshot residues on the surface of the object carrying the residues steadily decreases as storage time increases. Despite this fact, mercury can be preserved as a component of gunshot residues on the object (cotton cloth) under room conditions for more than 45 days. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 73, No. 3, pp. 320–325, May–June, 2006.     

利用有效液滴模型对超重元素α衰变半衰期的研究

利用有效液滴模型计算了偶偶超重核的α衰变半衰期,计算过程中采用了保持碎片体积守恒的不对称形状描述以及有效惯性系数计算Gamow势垒穿透因子.首先在质子数Z为88—98的区域检验了有效液滴模型,发现计算结果与实验符合得比较好.随后将此模型推广到Z≥100的情况,虽然只用了两个模型参数,计算结果与实验数据符合,说明有效液滴模型是计算偶偶超重核素α衰变半衰期的一个成功模型.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

A boundary element method for the calculation of noise barrier insertion loss in the presence of atmospheric turbulence

Atmospheric turbulence is an important factor that limits the amount of attenuation a barrier can provide in the outdoor environment. It is therefore important to develop a reliable method to predict its effect on barrier performance. The boundary element method (BEM) has been shown to be a very effective technique for predicting barrier insertion loss in the absence of turbulence. This paper develops a simple and efficient modification of the BEM formulation to predict the insertion loss of a barrier in the presence of atmospheric turbulence. The modification is based on two alternative methods: (1) random realisations of log-amplitude and phase fluctuations of boundary sources and (2) de-correlation of source coherence using the mutual coherence function (MCF). An investigation into the behaviours of these two methods is carried out and simplified forms of the methods developed. Some systematic differences between the predictions from the methods are found. When incorporated into the BEM formulation, the method of random realisations and the method of MCF de-correlation provide predictions that agree well with predictions by the parabolic equation method and by the scattering cross-section method on a variety of thin barrier configurations.     

溃坝问题的间断有限元方法

本文研究90年代初提出的Runge-Kutta间断Galerkin有限元方法，给出该方法的精度分析，通过经典算例验证该方法处理间断问题、捕捉锐利波形的能力，并将其推广到求解浅水问题．针对坝底无摩擦，无坡度的理想情形进行讨论，给出方溃坝和圆溃坝问题的数值模拟结果．     

THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

In this paper, we investigate the coupling of natural boundary element and finite element methods of exterior initial boundary value problems for hyperbolic equations. The governing equation is first discretized in time, leading to a time-step scheme, where an exterior elliptic problem has to be solved in each time step. Second, a circular artificial boundary FR consisting of a circle of radius R is introduced, the original problem in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formula are obtained in the infinite domainΩ2 outside circle of radius R. The coupled variational formulation is given. Only the function itself, not its normal derivative at artificial boundary ΓR, appears in the variational equation, so that the unknown numbers are reducedand the boundary element stiffness matrix has a few different elements. Such a coupled method is superior to the one based on direct boundary element method. This paper discusses finite element discretization for variational problem and its corresponding numerical technique, and the convergence for the numerical solutions. Finally, the numerical example is presented to illustrate feasibility and efficiency of this method.