快速数值差分递推算法在偏振模色散作用下光脉冲传输中的应用研究

利用快速数值差分递推算法对非线性薛定谔方程进行推导分析,得到一个同时考虑群速色散效应、非线性效应以及偏振模色散效应的时域快速数值差分递推关系,并将该方法应用于研究偏振模色散作用下的光脉冲传输,得到的计算结果分别与解析结果和传统的分步傅里叶算法所得结果进行对比,表明这种快速数值差分递推算法拥有很高的计算准确度,并且物理模型更加合理,完全适合分析计算偏振模色散作用下的光脉冲传输问题.利用该方法研究了光脉冲的畸变以及脉冲展宽问题,获得的计算结果为光纤传输系统的设计提供了参考.

海岸河口二维潮流可视化数学模型

经过多年的不间断的研究，逐步改进和完善了在海岸河口非常适用的、基于不规则三角形网格显式差分法的二维潮流数学模型，研制了潮流场可视化成品实用软件，将二者结合起来即构成了二维潮流可视化数学模型。可视化数学模型不仅使数学模型可视化、直观化，而且加快了数学模型的调试速度．方便了各种工程方案的比选，从而使用数学模型解决实际工程问题的效率大为提高。

不同物理和几何参数条件下滑坡要素的地质雷达探测响应研究

利用改进的时域有限差分算法,对不同物理和几何参数条件下各种滑坡要素组成的综合模型进行数值模拟,研究滑带厚度、滑带充填介质、滑体岩土类型、滑坡裂缝等的地质雷达探测响应效果。研究结果表明,具有各自物理和几何参数的滑坡要素与地质雷达响应特征之间存在一一对应的数学关联,通过这种关联可以从野外实测数据中有效提取和判读滑坡要素的类型及物理几何特性,为滑坡稳定性评价提供地球物理依据。为了展示该研究成果的适用性,以三峡库区将军滩滑坡的地质雷达野外实测数据为例进行解释推断,成功识别出滑坡体的滑带埋深及分布、滑体的裂缝发育程度等情况,为评价滑坡稳定性,合理优化滑坡治理方案提供了科学参考。

齐大山铁矿北采边坡治理稳定性分析

在充分认识齐大山铁矿北采下盘滑坡特征、工程和水文地质特征的基础上,用有限差分数值法(FLAC)和极限平衡法(SLOPE2000)对滑坡治理的多个方案进行了稳定性综合分析评价,从而优化出治理方案。

NUMERICAL SIMULATIONS OF EVOLUTIONS FROM WALL PULSE TO TURBULENT COHE RENT STRUCTURES

The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.

Study on fast calibration of a ULA's gain and phase error

Array calibration is important in engineering practice. In this paper, fast calibration methods for a ULA's gain and phase errors both in far and near fields are proposed. In the far field, using a single sound source without exact orientation, this method horizontally rotates the array exactly once, performs eigen value decomposition for the covariance matrix of received data, then computes the gain and phase error according to the formulas. In the near field, using the same single sound source, it is necessary to rotate the array horizontally at most three times, build equations according to geometric relations, then solve them. Using the formula proposed in this paper, spherical waves are modified into plane waves. Then eigen values decomposition is performed. These two calibration methods were shown to be valid by simulation and are fast, accurate and easy to use. Finally, an analysis of factors influencing estimation precision is given.

Improved Cotes Formula and Error Analysis

The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.

Effects of liquefaction-induced large lateral ground deformation on pile foundations

The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.

Two Unconditional Stable Schemes for Simulation of Heat Equation on Manifold Using DEC

To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.

Determining the Effective Refractive Index of AIGaAs-GaAs Slab Waveguide Based on Analytical and Finite Difference Method

Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics needs such numerical methods. These methods must be precise and useful in terms of memory capacity and time duration. In this paper, we study basic analytical and finite difference methods to determine the effective refractive index of AIGaAs-GaAs slab waveguide. Also, appropriate effective refractive index value is obtained with respect to number of grid points and number of matrix sizes. Finally, the validity of the obtained values by both methods is compared to using waveguide type.

Mechanical analysis of fixed geosynthetic technique of GRPS embankment

To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data.

MATLAB fitting for overcurrent relay and a new digital algorithm study

The overcurrent (OC) protection limit is set usually accorging to a OC protection setting table on digital integrated protection equipment in mine explode isolation high voltage (HV) vacuum switch. For digital integrated protection equipment, OC protection setting table must be converted to be a microcomputer algorithm. This paper first intro-duced a method of the fitting OC protection setting table to be OC relay inverse time characteristics equations using MATLAB least square fitting. On the basis of analyzing these fitting equations, a notion, “integral limit rate” was put forward initially and a OC in-verse time digital algorithm was developed. MATLAB simulation results and a digital signal processor (DSP) based digital integrated protection equipment running test indicate that this algorithm has less calculation amount, less taking up memory, high control accuracy, implements the no-grade setting of OC delay values, suits for all kinds of low-middle mi-crocomputer system implementation.

Modeling of low-frequency seismic waves in a shallow sea using the staggered grid difference method

Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships. To obtain the propagation characteristics of ship seismic waves, an algorithm for calculating Seismic waves at the seafloor is presented based on the staggered-grid finite difference method. The accuracy of the algorithm was tested by comparison with analytical solutions. Numerical simulation of seismic waves generated by a low-frequency point sotmd source in a typical shallow sea environment was carried out. Using various source frequencies and locations in the numerical simulation, we show that the seismic waves in the near field are composed mostly of transmitted S-waves and interface waves while transmitted P-waves are weak near the seafloor. However, in the far field, the wave components of the seismic wave are mainly normal modes and interface waves, with the latter being relatively strong in the waveforms, As the source frequency decreases, the normal modes become smaller and the interface waves dominate the time series of the seismic waves.

A numerical error analysis method and its application on a 4RRR parallel kinematic machine

To guarantee the accuracy of error analysis and evaluate the manufacturing tolerance s influence,anumerical error analysis method for parallel kinematic machines (PKMs) is presented in this paper.Quasi-Newton method and genetic algorithm are introduced for the forward kinematic solution.Based onthe inverse and forward kinematic solutions,the end-effector s error calculation procedure is developed.To solve the accuracy problem caused by the length and angular parameters' different units,a normalizationmethod is proposed based on the manufacturing tolerance.Comparison between the error analysis resultscalculated by the traditional method and the numerical method for a 4RRR PKM shows that,this numericalerror analysis method is more accurate,simpler,and can evaluate the machine s real error basedon the manufacturing tolerance.

Improved modal truncation error in the directly analytical method for damage identification of frame structures

The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equations were solved by the minimal least square method （viz. general inverse method）. It demonstrates that the feasibility and the accuracy of the present approach were impoved significantly, compared with the result of unimproved damage identification.

Numerical Homogenization of the Elliptic Problem with Rapidly Varying Tensor or Boundary Condition

Homogenization is concerned with obtaining the average properties of a material. The problem on its own has no easy solution, except in cases like the periodic case, when it can be obtained in closed form. In this paper we consider a numerical solution of the elliptic homogenization problem for the case of rapidly varying tensor or boundary conditions. The method makes use of an adaptive finite element method to correctly capture the rapid change in the tensor or boundary condition. In the numerical experiments we vary the mesh size and do a posteriori error analysis on test problems.

Numerical study of detonation shock dynamics using generalized finite difference method

The generalized finite difference method (GFDM) used for irregular grids is first introduced into the numerical study of thelevel set equation, which is coupled with the theory of detonation shock dynamics (DSD) describing the propagation of thedetonation shock front. The numerical results of a rate-stick problem, a converging channel problem and an arc channel prob-lem for specified boundaries show that GFDM is effective on solving the level set equation in the irregular geometrical domain.The arrival time and the normal velocity distribution of the detonation shock front of these problems can then be obtainedconveniently with this method. The numerical results also confirm that when there is a curvature effect, the theory of DSDmust be considered for the propagation of detonation shock surface, while classic Huygens construction is not suitable anymore.

Joint semiparametric mean-covariance model in longitudinal study

Semiparametric regression models and estimating covariance functions are very useful for longitudinal study. To heed the positive-definiteness constraint, we adopt the modified Cholesky decomposition approach to decompose the covariance structure. Then the covariance structure is fitted by a semiparametric model by imposing parametric within-subject correlation while allowing the nonparametric variation function. We estimate regression functions by using the local linear technique and propose generalized estimating equations for the mean and correlation parameter. Kernel estimators are developed for the estimation of the nonparametric variation function. Asymptotic normality of the the resulting estimators is established. Finally, the simulation study and the real data analysis are used to illustrate the proposed approach.

An unconditionally energy stable finite difference scheme for a stochastic Cahn-Hilliard equation

In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation, is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional.For the non-stochastic case, we develop an unconditionally energy stable difference scheme which is proved to be uniquely solvable. For the stochastic case, by adopting the same splitting of the energy functional, we construct a similar and uniquely solvable difference scheme with the discretized stochastic term. The resulted schemes are nonlinear and solved by Newton iteration. For the long time simulation, an adaptive time stepping strategy is developed based on both first- and second-order derivatives of the energy. Numerical experiments are carried out to verify the energy stability, the efficiency of the adaptive time stepping and the effect of the stochastic term.