Dynamic stress concentration in a ribbed plate using the acoustical wave propagator technique

This paper presents an explicit acoustical wave propagator technique to investigate the flexural wave scattering and dynamic stress concentration in a ribbed plate based on Mindlin plate theory. A scheme combining Chebyshev polynomial expansion and fast Fourier transformation is introduced to implement the operation of the acoustical wave propagator. The exact analytical solutions are also presented to demonstrate the validity of the present technique. The wave propagation patterns are used to investigate the effect of rib on the flexural wave scattering and distribution of dynamic stress concentration. Dynamic stress concentration factors around the discontinuities are examined in detail.     

基于声传播算子的声源定位实验模拟方法研究

声传播算子是一种高效的时域声场计算方法,它能够很方便地计算出给定系统参数下任意时刻任意位置的声场变化情况,本文采用这种方法计算所得的二维房间声场信息进行传声器阵列的声源定位仿真实验。计算结果表明,用该方法获取的阵列数据能有效地应用于阵列信号处理算法中,准确地估计出初始高斯脉冲声源的方向。声传播算子声场计算方法能为传声器阵列声源定位的实验提供方便,使得传声器阵列声源定位算法在不同混响时间的鲁棒性实验研究变得更加简捷。     

An improved shooting algorithm and its application to high-power fiber lasers

An improved shooting algorithm for high-power fiber laser is proposed based on the relation of the two-end boundary laser powers. By the improved shooting algorithm, the evolutions of the pump and laser powers along the fiber position for Yb-doped and Tm-doped high-power double-clad fiber lasers with the laser scattering loss are analyzed. The results prove that the improved shooting algorithm can be efficient even if the pump power is up to kilowatt level, and the initial value can all be fast convergent as long as the initial guessed laser power is smaller than the truth value. So, this improved shooting algorithm can be used in high-power double-clad fiber lasers efficiently.     

The concept of cyclic sound intensity and its application to acoustical imaging

This paper demonstrates how to take advantage of the cyclostationarity property of engine signals to define a new acoustical quantity, the cyclic sound intensity, which displays the instantaneous flux of acoustical energy in the angle-frequency domain during an average engine cycle. This quantity is attractive in that it possesses the ability of being instantaneous and averaged at the same time, thus reconciling two conflicting properties into a rigourous and unambiguous framework. Cyclic sound intensity is a rich concept with several original ramifications. Among other things, it returns a unique decomposition into instantaneous active and reactive parts. Associated to acoustical imaging techniques, it allows the construction of sound radiation movies that evolve within the engine cycle and whose each frame is a sound intensity map calculated at a specific time - or crankshaft angle - in the engine cycle. This enables the accurate localisation of sources in space, in frequency and in time (crankshaft angle). Furthermore, associated to cyclic Wiener filtering, this methodology makes it possible to decompose the overall radiated sound into several noise source contributions whose cyclic sound intensities can then be analysed independently.     

An extended Jacobi elliptic function rational expansion method and its application to ()-dimensional dispersive long wave equation

With the aid of computerized symbolic computation, a new elliptic function rational expansion method is presented by means of a new general ansatz, in which periodic solutions of nonlinear partial differential equations that can be expressed as a finite Laurent series of some of 12 Jacobi elliptic functions, is more powerful than exiting Jacobi elliptic function methods and is very powerful to uniformly construct more new exact periodic solutions in terms of rational formal Jacobi elliptic function solution of nonlinear partial differential equations. As an application of the method, we choose a (2+1)-dimensional dispersive long wave equation to illustrate the method. As a result, we can successfully obtain the solutions found by most existing Jacobi elliptic function methods and find other new and more general solutions at the same time. Of course, more shock wave solutions or solitary wave solutions can be gotten at their limit condition.     

一种改进的高性能Lorenz系统构造及其应用

Lorenz系统是一种最具有代表性、典型性的混沌模型之一, 一直被众多学者深入研究和广泛应用.为了获取结构和动力学行为更为复杂的混沌吸引子, 不断改善Lorenz系统已成为混沌动力系统研究中的重要课题之一. 为此, 本文提出了一个具有复杂系统动力学行为的改进的Lorenz系统, 并将其用于图像信息安全保护. 在现有各种改进的Lorenz系统的基础上, 首先通过增加Lorenz系统的控制参数和改变非线性项相结合的方法构造出一种新的Lorenz 混沌系统; 其次采用微分动力系统方法深入研究该系统并获得与Lorenz系统、Bao系统、Tee系统和Y系统等具有相似的耗散性、对称性、稳定性, 以及更加复杂的混沌特性和动力学行为, 同时分析该系统所产生随机序列具有良好的相关性和复杂性; 最后将其所产生的离散伪随机序列用于图像置乱和扩散加密, 通过对图像加密结果的相邻像素相关性分析、灰度空间相关特性不确定性分析、抗差分攻击以及密钥敏感性测试, 表明本文所构造的改进的Lorenz系统应用于图像加密能获得相对较高的安全性.     

A propagator theory applied to wave mechanics in phase space

The fact that the classical Liouville equation can be analyzed as a dynamical equation in Hilbert-Koopman (HK.) space is used in order to develop a perturbative method for the wave mechanics in phase space: an explicit solution of the Liouville equation inqp representation is exhibited. The connection between the solution obtained and the dynamics of correlations is established by computing theqp-kp transformation function in HK space. To elucidate the method, an application is presented and the result compared to that available in the literature.     

声波在固体板中的多径传播及其时间反转声场

采用无穷长条形阵元在固体板自由表面激发,研究了固体板中声波(纵波和横波)传播的多径效应,在射线近似下深入分析了柱面纵波和横波在固体板中的多次反射及波型转换特性,给出了简明解析表达式;将时间反转法应用于固体板中声波的传播,在理论和实验上分析了时间反转声场规律,证明了从不同阵元发出的沿不同路径传播的声波时间反转后同时同相到达原接收点,表明时间反转法能自动补偿固体板中由于多径效应造成的波形畸变;还通过聚焦增益和主副瓣比定量地分析了时间反转声场的自适应聚焦过程,考查了焦点位置和换能器阵列孔径对聚焦效果的影响,得到了理论和实验相符的结果。     

等离子体的SO-FDTD算法和对电磁波反射系数的计算分析

将一类色散介质的介电常数写成有理分式函数形式，进而导出FDTD中电位移矢量D和电场强度E之间的关系，形成SO-FDTD方法. 应用该方法计算了冷等离子体平板对电磁波的反射系数，通过与解析结果的比较，验证了该算法的高效性和高精度，同时，应用SO-FDTD方法，计算了等离子体层对垂直入射电磁波的反射系数，结果表明：等离子体厚度、电子密度、电子密度的分布形式和入射波频率是影响反射系数的重要因素. 关键词： 等离子体 电磁波 FDTD方法     

Fractional-harmonic components in an acoustical wave in a liquid

The amplitudes and frequencies of the components of the acoustical spectrum of an acoustical wave, generated in a non-viscous liquid filling a cavity, are calculated by solving a boundary value problem. Apart from forced oscillations, the frequencies of the acoustical spectrum are equal to the frequency of the fundamsntal oscillation of the liquid in the cavity and its higher harmonics. If the frequency of the driving acoustical wave coincides with one of these (proper) frequencies, thefractional harmonic components appear. The amplitudes of the component oscillations decrease monotonically as the absolute value of the difference between the frequency of the driving acoustical wave and the frequency of the respective oscillation is increased. The derived relations are compared with the results of some published measurements.     

驻波实验仪的一种改进

设计了一种新型的驻波实验装置,该装置采用振动频率可随外加交流电压频率变化的振动子作为振源,再利用单片机及功率驱动电路产生各种所需的不同频率的电压,来控制振源的振动频率.用它可以验证传统的驻波实验装置不能验证的驻波频率f与驻波波长λ之间的关系.     

An order-N electronic structure theory with generalized eigenvalue equations and its application to a ten-million-atom system

A linear algebraic theory called the 'multiple Arnoldi method' is presented and realizes large-scale (order-N) electronic structure calculations with generalized eigenvalue equations. A set of linear equations, in the form of (zS - H)x = b, are solved simultaneously with multiple Krylov subspaces. The method is implemented in a simulation package ELSES (www.elses.jp) with tight-binding-form Hamiltonians. A finite-temperature molecular dynamics simulation is carried out for metallic and insulating materials. A calculation with 10(7) atoms was realized by a workstation. The parallel efficiency is shown up to 1024 CPU cores.     

An improved model for soil ionization around grounding system and its application to stratified soil

An improved model for taking into account the effect of the soil ionization around grounding system under lightning strike is proposed in this paper. In this model, the soil ionization region is assumed to retain 7% of its pre-ionization resistivity, which is consistent with the experimental results on soil ionization found in literature (Trans. SA Inst. Electr. Eng. (1988) 63; AIEE Trans. 61 (1942) 349; Proc. IEE 121(2) (1974) 123) and our own laboratory experiments (Time domain modelling of the response of grounding systems subjected to lightning currents, Licenciate Thesis, Uppsala University, 2003). Compared with modelling the soil ionization as an increase in the size of the ground conductor, the model presented here will not overestimate the beneficial influence of the soil ionization in reducing the ground potential rise, especially in high resistivity soil. The model is also applied to study the transient behaviour of grounding conductors in stratified soil under lightning strike including soil ionization. It shows that making the grounding conductor to penetrate the lower resistivity soil layer could help to decrease the ground potential rise at the injection point several times.     

An efficient physical device-circuit simulator and its application to accurate design of millimeter wave oscillators

A rigorous one-dimensional physical device-circuit simulator is developed to accurately determine the transport properties and the electrical performance of semiconductor devices embedded into a passive circuit. The simulator is well suited to study the effects of device geometry, doping level, bias voltage, and mounting structure on the performance and to compare the accuracy and the computational efficiency of different physical models. The model is applied to determine the application limits of the drift-diffusion and hydrodynamic models when they are used to characterize the operation and to optimize the structure of cm- and mm-wave oscillators with two-terminal devices, and to produce a set of design curves which indicate the performance limits of different devices under various operating conditions.     

An improved variational method

In order to improve the unitarity of the S-matrix, an improved variational formulism is derived by proposing new generating functionals and adopting proper asymptotic boundary conditions for trial relative wave functions. The formulas with the weighted line-column balance for the single-channel and multi-channel scatterings, where the non-central interaction is implicitly considered, are presented. A numerical check is performed with a soluble model in a four coupled channel scattering problem. The result shows that the high accuracy and the unitarity of the S-matrix are reached.     

An improved dynamic subgrid-scale model and its application to large eddy simulation of rotating channel flows

Rotating turbulence occurs extensively in nature and engineering circumstances. Meanwhile, understanding physical mechanisms of the rotating turbulence is important to the fundamental research of turbulence. The turbulent flow in rotating frames undergoes two kinds of Coriolis force effects. First, a secondary flow is induced in the case that there is a mean vorticity component perpendicular to the rotating axis. Second, there are augmenting or suppressing effects on the turbulence if there i…     

An improved variational method

In order to improve the unitarity of the S-matrix, an improved variational formulism is derived by proposing new generating functionals and adopting proper asymptotic boundary conditions for trial relative wave functions. The formulas with the weighted line-column balance for the single-channel and multi-channel scatterings, where the non-central interaction is implicitly considered, are presented. A numerical check is performed with a soluble model in a four coupled channel scattering problem. The result shows that the high accuracy and the unitarity of the S-matrix are reached.     

具有不规则自由表面的粘弹性介质中声波模拟的有限差分算法

声波在具有不规则自由表面的粘弹性介质传播的数值模拟工作在地震勘探,地震预测中非常重要,特别是当模型有起伏的自由界面和较强的衰减特性时更是如此。基于镜像方法和直接方法,本文发展了一种二维有限差分算法,可以模拟不规则自由表面引起的声波散射问题。该方法将自由表面条件与速度应力方程结合求解粘弹性波动方程,在垂直和水平自由段及拐点处设置相应的边界条件。该方法假设自由表面穿过剪切应力和相应参数的网格点。为了提高计算精度,分别计算了水平和垂直方向上应力镜像值,而对质点速度,采用先水平方向后垂直方向分析进行镜像计算和更新。将粘弹性水平自由表面镜像方法和不同倾斜度平滑自由表面的改进方法的数值结果进行对比,并验证了算法的精度。     

离散Arnold变换改进及其在图像置乱加密中的应用

为了改善传统二维Arnold变换用于图像置乱加密的效果,提出了离散Arnold变换的改进方法,并将其用于图像置乱加密测试研究.该方法利用现有离散标准映射的构造思想,将传统离散二维Arnold变换表达式中第一个变换表达式所对应变换结果非线性融入第二个变换表达式,实现经典离散二维Arnold变换的非线性去拟仿射化修改,以便快速改善图像置乱加密效果.数学证明改进方法不再保持现有离散二维Arnold变换所具有的拟仿射不变性,但是改进变换仍是一种具有周期性的可逆映射,将其用于图像置乱加密时,利用其周期性或逆变换能恢复置乱前原图像.大量实验结果表明,本文所建议的改进方法是有效的,相比现有的离散Arnold变换更具有实用价值意义.