随机噪声对经验模态分解非线性信号的影响

采用Monte-Carlo随机模拟方法来研究外部噪声对经验模态分解非线性信号的影响.结果表明:噪声对低阶特征模态函数(IMF)影响较为明显,对高阶IMF影响较小;白噪声强度系数越大,分解出的IMF纯噪声分量阶数越多;用含噪声信号减去经验模态分解后的主要IMF噪声分量,可较为明显地削弱噪声的影响;含噪声响应的最大Lyapunov指数比不含噪声响应的最大Lyapunov指数小     

The influence of cell geometry on the Godunov scheme applied to the linear wave equation

By studying the structure of the discrete kernel of the linear acoustic operator discretized with a Godunov scheme, we clearly explain why the behaviour of the Godunov scheme applied to the linear wave equation deeply depends on the space dimension and, especially, on the type of mesh. This approach allows us to explain why, in the periodic case, the Godunov scheme applied to the resolution of the compressible Euler or Navier–Stokes system is accurate at low Mach number when the mesh is triangular or tetrahedral and is not accurate when the mesh is a 2D (or 3D) cartesian mesh. This approach confirms also the fact that a Godunov scheme remains accurate when it is modified by simply centering the discretization of the pressure gradient.     

Effect of shot peening on the microstructure of laser hardened 17-4PH

In order to investigate the influence of shot peening on microstructure of laser hardened steel and clarify how much influence of initial microstructure induced by laser hardening treatment on final microstructure of laser hardened steel after shot peening treatment, measurements of retained austenite, measurements of microhardness and microstructural analysis were carried out on three typical areas including laser hardened area, transitional area and matrix area of laser hardened 17-4PH steel. The results showed that shot peening was an efficient cold working method to eliminate the retained austenite on the surface of laser hardened samples. The surface hardness increased dramatically when shot peening treatments were carried out. The analyses of microstructure of laser hardened 17-4PH after shot peening treatment were carried out in matrix area and laser hardened area via Voigt method. With the increasing peening intensity, the influence depth of shot peening on hardness and microstructure increased but the surface hardness and microstructure did not change when certain peening intensity was reached. Influence depth of shot peening on hardness was larger than influence depth of shot peening on microstructure due to the kinetic energy loss along the depth during shot peening treatment. From the microstructural result, it can be shown that the shot peening treatment can influence the domain size and microstrain of treated samples but laser hardening treatment can only influence the microstrain of treated samples.     

Two-electron entanglement in elliptically deformed quantum dots

Entropic entanglement measures of a two-dimensional system of two Coulombically interacting particles confined in an anisotropic harmonic potential are discussed in dependence on the anisotropy and the interaction strength. The harmonic approximation appears exact in the strong interaction limit, allowing determination of the asymptotic expression for the linear entropy. Entanglement properties are dramatically influenced by the anisotropy of the confining potential in the strong-correlation regime.     

Unconditionally Stable Complex Envelope Wave Equation PML Algorithm for Band Limited FDTD Simulations

Unconditionally stable complex envelope (CE) perfectly matched layer (PML) absorbing boundary conditions (ABCs) are presented for truncating the scalar wave-equation finite difference time domain (WE-FDTD) grids. The formulations are based on incorporating the alternating direction implicit (ADI) scheme into the CE FDTD implementations of the scalar wave-equation derived in the PML region at the domain boundaries. Numerical example carried out in two dimensional domain shows that the proposed formulations are more accurate than the classical ADI scalar wave equation PML formulations when it is used for modelling band limited electromagnetic applications.     

Exact Solution of a Test Particle in Presence of Thick Domain Walls

We have obtained the exact solution of the equations of motion of a test particle near a thick domain walls. From the solution it has been shown that the domain walls have repulsive gravitational fields.     

On the spurious eigenvalues for a concentric sphere in BIEM

In this paper, the null-field integral equation method is employed to study the occurring mechanism of spurious eigenvalues for a concentric sphere. By expanding the fundamental solution into degenerate kernels and expressing the boundary density in terms of spherical harmonics, all boundary integrals can be analytically determined. It is noted that our null-field integral formulation can locate the collocation point on the real boundary thanks to the degenerate kernel. In addition, the spurious eigenvalues are parasitized in the formulations while true eigensolutions are dependent on the boundary condition such as the Dirichlet or Neumann problem. By using the updating term and updating document of singular value decomposition (SVD) technique, true and spurious eigenvalues can be extracted out, respectively. Besides, true and spurious boundary eigenvectors are obtained in the right and left unitary vectors in the SVD structure of the influence matrices. This finding agrees with that of 2D cases.     

Fast evaluation of Helmholtz potential on graphics processing units (GPUs)

This paper presents a parallel algorithm implemented on graphics processing units (GPUs) for rapidly evaluating spatial convolutions between the Helmholtz potential and a large-scale source distribution. The algorithm implements a non-uniform grid interpolation method (NGIM), which uses amplitude and phase compensation and spatial interpolation from a sparse grid to compute the field outside a source domain. NGIM reduces the computational time cost of the direct field evaluation at N observers due to N co-located sources from O(N²) to O(N) in the static and low-frequency regimes, to O(N log N) in the high-frequency regime, and between these costs in the mixed-frequency regime. Memory requirements scale as O(N) in all frequency regimes. Several important differences between CPU and GPU implementations of the NGIM are required to result in optimal performance on respective platforms. In particular, in the CPU implementations all operations, where possible, are pre-computed and stored in memory in a preprocessing stage. This reduces the computational time but significantly increases the memory consumption. In the GPU implementations, where handling memory often is a critical bottle neck, several special memory handling techniques are used to accelerate the computations. A significant latency of the GPU global memory access is hidden by implementing coalesced reading, which requires arranging many array elements in contiguous parts of memory. Contrary to the CPU version, most of the steps in the GPU implementations are executed on-fly and only necessary arrays are kept in memory. This results in significantly reduced memory consumption, increased problem size N that can be handled, and reduced computational time on GPUs. The obtained GPU–CPU speed-up ratios are from 150 to 400 depending on the required accuracy and problem size. The presented method and its CPU and GPU implementations can find important applications in various fields of physics and engineering.     

基于改进格林元的压裂水平井动态分析方法

以建立高效的动态分析方法为出发点,以边单元作为求解点,改进传统的格林元方法,减少未知数和求解矩阵维度;并提出基于改进格林元的加密网格加密方法,保证考虑复杂裂缝网络的压裂水平井动态模拟的早期精度.退化模型与半解析解、数值模拟结果进行对比,验证本文基于加密网格的改进格林元方法的准确性和动态分析的高效性.最后进行动态响应的敏感性分析,结果表明：①格林元方法是一种高精度的动态模拟方法,将求解节点设置在网格的边上可以提高压裂水平井动态模拟的速度;②改进格林元方法的加密基于叠加原理,不需要通过插值近似,其求解精度高.在相同加密网格条件下,基于本文改进格林元方法的加密效果比有限差分加密效果更佳;③复杂裂缝导流能力、改造区渗透率提高倍数、改造区大小等参数对压裂水平井动态特征影响较大,在动态分析和参数反演时,应着重考虑这些因素的影响.     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象,发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式,避免在单元内进行速度平均,从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比,网格的边界通量由节点处的多维黎曼求解器提供,克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增,给出稀疏波出现的合理判据,从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象,同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.