有限元单元计算子程序的OpenMP并行化

Intel和AMD双核乃至4核处理器的推出,使得并行计算已经普及到PC机。为了充分利用多核,需要对原有程序进行多线程改造,使其充分利用多核处理带来的性能提升。该文利用共享存储编程的工业标准OpenMP对有限元方法涉及的单元计算子程序进行了并行化实现。在机群的一个双CPU的SMP节点上的测试表明,共享并行化使得该单元子程序的性能提高了一倍。

有限元单元计算予程序的OpenMP并行化

Intel和AMD双核乃至4核处理器的推出，使得并行计算已经普及到PC机。为了充分利用多核，需要对原有程序进行多线程改造，使其充分利用多核处理带来的性能提升。该文利用共享存储编程的工业标准OpenMP对有限元方法涉及的单元计算子程序进行了并行化实现。在机群的一个双CPU的SMP节点上的测试表明，共享并行化使得该单元子程序的性能提高了一倍。

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

On the radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas

The radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas of finite depth is investigated. The analytical expressions for the radiated and diffracted potentials are derived as infinite series by use of the method of separation of variables. The unknown coefficients in the series are determined by the eigenfunction expansion matching method. The expressions for wave forces, hydrodynamic coefficients and reflection and transmission coefficients are given and verified by the boundary element method. Using the present analytical solution, the hydrodynamic influences of the angle of incidence, the submergence, the width and the thickness of the structure on the wave forces, hydrodynamic coefficients, and reflection and transmission coefficients are discussed in detail.

Reflection and transmission properties of double submerged rectangular blocks

An analytic method is used to study the reflection and transmission coefficients of the double submerged rectangular blocks (DSRBs) in oblique waves.. The scattering potentials are obtained by means of the eigenfunction expansion method, and expressions for the reflection and transmission coefficients are determined. The boundary element method is employed to verify the correctness of the present analytical method. The DSRBs have better performance than the single submerged rectangular block (SSRB) in certain cases. The reflection and transmission properties of the DSRBs are investigated for some specific cases, and the influences of the geometric parameters are also presented.

Scattering of linear water waves by a fixed and infinitely long rectangular structure parallel to a vertical wall in oblique seas

The scattering of linear water waves by an infinitely long rectangular structure parallel to a vertical wall in oblique seas is investigated. Analytical expressions for the diffracted potentials are derived using the method of separation of variables. The unknown coefficients in the expressions are determined through the application of the eigenfunction expansion matching method. The expressions for wave forces on the structure are given. The calculated results are compared with those obtained by the boundary element method. In addition, the influences of the wall, the angle of wave incidence, the width of the structure, and the distance between the structure and the wall on wave forces are discussed. The method presented here can be easily extended to the study of the diffraction of obliquely incident waves by multiple rectangular structures.

Transfer property of electron pulse in a femto-second electron diffraction system

An ultra-fast electron diffraction system has been designed. The static and dynamic characters of an electron pulse with 150 fs temporal dispersion are studied during its transmission in the whole ultra-fast electron diffraction system, including the size of the electron spot, temporal dispersion, distribution of azimuths and elevation angles. The initial status of the photoelectrons are put down by Monte Carlo method, both the two dimensional and three dimensional electric fields are calculated by finite difference method and the magnetic flux are. calculated by finite element method.

有限元子域法分析土石混合介质的饱和导水率

在黄土高原水蚀风蚀交错地区,碎石的存在影响着土壤的水力学特性,因此研究土石混合介质的渗流运动对该地区的植被恢复建设具有重要的指导意义。虽然用有限元法分析地下水渗流的工作很多,但对土石混合介质而言,由于有限元法的解题规模迅速增加,其研究甚少。本模型采用子域法矩形单元,假定土石体中石块含量增加不影响土的孔隙率且石块不透水分析了土石混合介质中的渗流过程。结果表明,影响土石混合介质饱和导水率的主要因素是土石介质的含石率;石块的大小基本不影响混合介质饱和导水率;石块分布方式增加水流路径时,会影响混合介质饱和导水率。这一数值分析结论与土石体渗流问题的部分理论结果和试验结果吻合良好。本模型子域法计算效率高,结果较准确,可用于较大规模的土石混合介质渗流分析,还可推广用于土石混合介质的三维分析和多级子域法有限元技术以获得更大规模的解题效果。

Theoretical and experimental investigations of flexural wave propagation in periodic grid structures designed with the idea of phononic crystals

The propagation characteristics of fiexural waves in periodic grid structures designed with the idea of phononic crystals are investigated by combining the Bloch theorem with the finite element method. This combined analysis yields phase constant surfaces, which predict the location and the extension of band gaps, as well as the directions and the regions of wave propagation at assigned frequencies. The predictions are validated by computation and experimental analysis of the harmonic responses of a finite structure with 11 × 11 unit cells. The fiexural wave is localized at the point of excitation in band gaps, while the directional behaviour occurs at particular frequencies in pass bands. These studies provide guidelines to designing periodic structures for vibration attenuation.

Design and numerical analysis for low birefringence silica on silicon waveguides

A new fabrication technology for three-dimensionally buried silica on silicon optical waveguide based on deep etching and thermal oxidation is presented. Using this method, a silicon layer is left at the side of waveguide. The stress distribution and effective refractive index are calculated by using finite element method and finite different beam propagation method, respectively. The results indicate that the stress of silica on silicon optical waveguide fabricated by this method can be matched in parallel and vertical directions and stress birefringence can be effectively reduced due to the side-silicon layer.

Design and optimization of highly nonlinear low-dispersion crystal fiber with high birefringence for four-wave mixing

We have proposed a novel type of photonic crystal fiber (PCF) with low dispersion and high nonlinearity for four-wave mixing. This type of fiber is composed of a solid silica core and a cladding with a squeezed hexagonal lattice elliptical airhole along the fiber length. Its dispersion and nonlinearity coefficient are investigated simultaneously by using the full vectorial finite element method. Numerical results show that the proposed highly nonlinear low-dispersion fiber has a total dispersion as low as +/- 2.5 ps nm(-1) km(-1) over an ultrabroad wavelength range from 1.43 to 1.8 mu m, and the corresponding nonlinearity coefficient and birefringence are about 150 W-1 km(-1) and 2.5 x 10(-3) at 1.55 mu m, respectively. The proposed PCF with low ultraflattened dispersion, high nonlinearity, and high birefringence can have important application in four-wave mixing. (C) 2010 Optical Society of America

Design of Topas photonic bandgap fiber with high birefringence and low confinement loss

A highly birefringent hollow-core photonic bandgap fiber based on Topas cyclic olefin copolymer is designed. The rhombic hollow-core with rounded corners is formed by omitting four central air holes of the cladding structure. The guided modes, birefringence and confinement loss of the fiber are investigated by using the full-vector finite element method. A high phase birefringence of the order of 10(-3), a group birefringence of the order of 10(-2) and confinement loss less than 0.1 dB/km are obtained at the central wavelength (1.55 mu m) range of the bandgap for fiber with seven rings of air holes in the cladding region. (C) 2010 Elsevier B.V. All rights reserved.

High birefringent rhombic-hole photonic crystal fibers

High birefringence induced by rhombic air-hole photonic crystal fibers (PCFs) is numerically analyzed by using the finite-element method. The birefringence of a few kinds of PCFs was investigated with different parameters related to rhombic holes, including the rhombic-hole shape, size, and spacing. It was found that the birefringence of the proposed rhombic-hole PCF in this study is relatively larger than that of an elliptical-hole PCF with the same air-filling fraction (f = 0.0375) when the ratio of the rhombic-hole diagonal length is equal to the elliptical-hole ellipticity. (C) 2010 Optical Society of America

The main factor influencing the tensile properties of surface nano-crystallized graded materials

The improved mechanical properties of surface nano-crystallized graded materials produced by surface severe plastic deformation ((SPD)-P-2) are generally owing to the effects of the refined structure, work-hardened region and compressive residual stress. However, during the (SPD)-P-2 process, residual stress is produced simultaneously with work-hardened region, the individual contribution of these two factors to the improved mechanical properties remains unclear. Numerical simulations are carried out in order to answer this question. It is found that work hardening predominates in improving the yield strength and the ultimate tensile strength of the surface nano-crystallized graded materials, while the influence of the residual stress mainly emerges at the initial stage of deformation and decreases the apparent elastic modulus of the surface nano-crystallized graded materials, which agrees well with the experimental results. (C) 2010 Elsevier B.V. All rights reserved.

Peeling behavior of a bio-inspired nano-film on a substrate

A peeling model is proposed to analyze the peeling properties of bio-mimetic nano-films using the finite element method (FEM) and theoretical approach. The influences of the nano-film's adhesion length, thickness, elastic modulus, roughness and peeling angle on the peeling force were considered as well as the effect of the viscoelastic behavior. It has been found that the effective adhesion length, at which the peeling force attained maximum, was much smaller than the real length of nano-films; and the shear force dominated in the case of smaller peeling angles, whereas, the normal force dominated at larger peeling angles. The total peeling force decreased with an increasing peeling angle. Two limiting values of the peeling-off force can be found in the viscoelastic model, which corresponds to the smaller and larger loading rate cases. The effects of nano-film thickness and Young's modulus on peeling behaviors were also discussed. The results obtained are helpful for understanding the micro-adhesion mechanisms of biological systems, such as geckos. (C) 2010 Elsevier Ltd. All rights reserved.