不连续问题的无网格分域算法

提出一种无网格方法中采用分域的思想处理材料和位移不连续问题的方法。该方法将求解域沿不连续面进行分域,通过使用两种转换矩阵使子域交界面上的位移连续性得到满足;采用分块矩阵法计算转换后的刚度矩阵,所得刚度矩阵仍具有稀疏、带状性。可采用与有限元耦合的方式施加本质边界条件。编制了该算法的计算程序,通过对材料不连续悬臂梁弯曲问题的分析和单边裂纹板裂纹张开位移的计算,验证了该算法的正确性和有效性。

非规则复杂域等值填充图的快速绘制方法

针对非规则复杂区域填充等值线图的绘制问题,提出了一种非规则的、复杂区域填充等值线图绘制算法。算法基本思想是:首先应用环形矩形域分割数据点;然后分区逐步插入新点快速生成二维约束Delaunay三角网格化;最后应用三叉树递归原理,快速等值剖分Delaunay三角形,颜色填充绘制等值域。通过研究实例表明,该方法具有很好的实时显示与应用效果。

Optimal Boundary Control Method for Domain Decomposition Algorithm

To study the domain decomposition algorithms for the equations of elliptic type, the method of optimal boundary control was used to advance a new procedure for domain decomposition algorithms and regularization method to deal with the ill posedness of the control problem. The determination of the value of the solution of the partial differential equation on the interface——the key of the domain decomposition algorithms——was transformed into a boundary control problem and the ill posedness of the control problem was overcome by regularization. The convergence of the regularizing control solution was proven and the equations which characterize the optimal control were given therefore the value of the unknown solution on the interface of the domain would be obtained by solving a series of coupling equations. Using the boundary control method the domain decomposion algorithm can be carried out.

ANALYSIS OF WAVEGUIDE PROBLEMS USING A RELAXED ITERATIVE DOMAIN DECOMPOSITION METHOD COMBINED WITH MULTIFRONTAL ALGORITHM

In this paper, an absorbing Fictitious Boundary Condition (FBC) is presented to generate an iterative Domain Decomposition Method (DDM) for analyzing waveguide problems.The relaxed algorithm is introduced to improve the iterative convergence. And the matrix equations are solved using the multifrontal algorithm. The resulting CPU time is greatly reduced.Finally, a number of numerical examples are given to illustrate its accuracy and efficiency.

A Bregman adaptive sparse-spike deconvolution method in the frequency domain

To improve the anti-noise performance of the time-domain Bregman iterative algorithm,an adaptive frequency-domain Bregman sparse-spike deconvolution algorithm is proposed.By solving the Bregman algorithm in the frequency domain,the influence of Gaussian as well as outlier noise on the convergence of the algorithm is effectively avoided.In other words,the proposed algorithm avoids data noise effects by implementing the calculations in the frequency domain.Moreover,the computational efficiency is greatly improved compared with the conventional method.Generalized cross validation is introduced in the solving process to optimize the regularization parameter and thus the algorithm is equipped with strong self-adaptation.Different theoretical models are built and solved using the algorithms in both time and frequency domains.Finally,the proposed and the conventional methods are both used to process actual seismic data.The comparison of the results confirms the superiority of the proposed algorithm due to its noise resistance and self-adaptation capability.

Flotation bubble image segmentation based on seed region boundary growing

Segmenting blurred and conglutinated bubbles in a flotation image is done using a new segmentation method based on Seed Region and Boundary Growing（SRBG）.Bright pixels located on bubble tops were extracted as the seed regions.Seed boundaries are divided into four curves：left-top,right-top,right-bottom, and left-bottom.Bubbles are segmented from the seed boundary by moving these curves to the bubble boundaries along the corresponding directions.The SRBG method can remove noisy areas and it avoids over- and under-segmentation problems.Each bubble is segmented separately rather than segmenting the entire flotation image.The segmentation results from the SRBG method are more accurate than those from the Watershed algorithm.

Parallel computation of unified finite-difference time-domain for underwater sound scattering

In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.

Calculation and application of full-wave airborne transient electromagnetic data in electromagnetic detection

Airborne electromagnetic transient method enjoys the advantages of high-efficiency and the high resolution of electromagnetic anomalies,especially suitable for mining detection around goaf areas and deep exploration of minerals.In this paper,we calculated the full-wave airborne transient electromagnetic data,according to the result of numerical research,the advantage of switch-off time response in electromagnetic detection was proofed via experiments.Firstly,based on the full-wave airborne transient electromagnetic system developed by Jilin University(JLU-ATEMI),we proposed a method to compute the full-waveform electromagnetic(EM)data of 3D model using the FDTD approach and convolution algorithm,and verify the calculation by the response of homogenous half-space.Then,through comparison of switch-off-time response and off-time response,we studied the effect of ramp time on anomaly detection.Finally,we arranged two experimental electromagnetic detection,the results indicated that the switch-off-time response can reveal the shallow target more effectively,and the full-waveform airborne electromagnetic system is an effective technique for shallow target detection.

Parallel programming characteristics of a DSP-based parallel system

This paper firstly introduces the structure and working principle of DSP-based parallel system, parallel accelerating board and SHARC DSP chip. Then it pays attention to investigating the system’s programming characteristics, especially the mode of communication, discussing how to design parallel algorithms and presenting a domain-decomposition-based complete multi-grid parallel algorithm with virtual boundary forecast (VBF) to solve a lot of large-scale and complicated heat problems. In the end, Mandelbrot Set and a non-linear heat transfer equation of ceramic/metal composite material are taken as examples to illustrate the implementation of the proposed algorithm. The results showed that the solutions are highly efficient and have linear speedup.

Multi-path routing algorithm in WSN using an improvedparticle swarm optimization

To slove the problems of constrained energy and unbalanced load of wireless sensor network(WSN)nodes,a multipath load balancing routing algorithm based on neighborhood subspace cooperation is proposed.The algorithm adopts the improved particle swarm optimization(PSO)algorithm,takes the shortest distance and minimum energy consumption as optimization target and divides the nodes in one-hop neighborhood near the base station area into different regions.Furthermore,the algorithm designs a fitness function to find the best node in each region as a relay node and forward the data in parallel through the different paths of the relay nodes.The simulation results show that the proposed algorithm can reduce energy consumption and average end-to-end delay,balance network load and prolong network lifetime effectively.

Image Segmentation Using an Improved Watershed Algorithm

As watershed algorithm suffers from over-segmentation problem, this paper presented an efficient method to resolve this problem. First, pre-process of the image using median filter is made to reduce the effect of noise. Second, watershed algorithm is employed to provide initial regions. Third, regions are merged according to the information between the region and boundary. In the merger processing based on the region information, an adaptive threshold of the difference between the neighboring regions is used as the region merge criteria, which is based on the human visual character. In the merger processing on the boundary information, the gradient is used to judge the true boundary of the image to avoid merging the foreground with the background regions. Finally, post-process to the regions using mathematical morphology open and close filter is done to smooth object boundaries. The experimental results show that this method is very efficient.

Expert Classification Method Based on Patch-Based Neighborhood Searching Algorithm

The expert knowledge has been widely used to improve the remotely sensed classification accuracy. Generally, the ex-pert classification system mainly depends on DEM and some thematic maps. The spatial relationship information in pixel level was commonly introduced into the expert classification. Because the geographic objects were found spatially dependent relationship to a certain degree, the commonly used basic unit of spatial relationship information in pixel greatly limited the efficiency of spatial in-formation. A patch-based neighborhood searching algorithm was proposed to implement the expert classification. The homogene-ous spectral unit, patch, was used as the basic unit in the spatial object granularity, and different types of patches' relationship in-formation were obtained through a spatial neighborhood searching algorithm. And then the neighborhood information and DEM data were added into the expert classification system and used to modify the primitive classification errors. In this case, the classi-fication accuracies of wetland, grassland and cropland were obviously improved. In this work, water was used as base object, and different types of water extraction methods were tested to get a result in a high accuracy.

The Novel Spectrum Allocation Algorithm in Cognitive Radio System

In order to reduce the disturbance on an authorizing user and lower the competition between cognitive users, assure the normal communication of a cognitive radio system, reliability theory is applied to describe if a channel can be used by a cognitive user or not and the probability that the channel is continually used for a period. Three aspects including space, time domain and frequency domain are united for the research on the distribution of frequency spectrum. The simulation result shows that, in the space domain, time domain, frequency domain algorithm, the transmitted data volume and the total throughput of the system are superior to those in greedy algorithm and time domain—frequency domain algorithm, the novel algorithm is helpful to reduce the disturbance caused by a cognitive user to an authorizing user and lower the competition between cognitive users, this simulation result shows that the proposed algorithm is effective.

OPTIMAL DESIGN OF SMART ANTENNA ARRAY

This letter investigates an efficient design procedure integrating the Genetic Algorithm (GA) with the Finite Difference Time Domain (FDTD) for the fast optimal design of Smart Antenna Arrays (SAA). The FDTD is used to analyze SAA with mutual coupling. Then,on the basis of the Maximal Signal to Noise Ratio (MSNR) criteria, the GA is applied to the optimization of weighting elements and structure of SAA. Finally, the effectiveness of the analysis is evaluated by experimental antenna arrays.

High performance parallel computing of large eddy simulation of the flow in a curved duct with square cross section

Large eddy simulation(LES) cooperated with a high performance parallel computing method is applied to simulate the flow in a curved duct with square cross section in the paper. The method consists of parallel domain decomposition of grids, creation of virtual diagonal bordered matrix, assembling of boundary matrix, parallel LDL^T decomposition, parallel solving of Poisson Equation, parallel estimation of convergence and so on. The parallel computing method can solve the problems that are difficult to solve using traditional serial computing. Furthermore, existing microcomputers can be fully used to resolve some large-scale problems of complex turbulent flow.

Prediction of Instability Separation Modes and Its Application in Practical Dynamic Security Region

The transient critical boundary of dynamic security region (DSR) can be approximated by a few hyper planes correlated with instability separation modes. A method to fast predict instability separation modes is proposed for DSR calculation in power injection space. The method identifies coherent generation groups by the developed K-medoids algorithm, taking a similarity matrix derived from the reachability Grammian as the index. As an experimental result, reachability Grammian matrices under local injections are approximately invariant. It indicates that the generator coherency identifications are nearly consistent for different injections. Then instability separation modes can be predicted at the normal operating point, while average initial acceleration is considered as the measure of the critical generator group to amend the error. Moreover, based on these predicted instability separation modes, a critical point search strategy for DSR calculation is illustrated in the reduced injection space of the critical generators. The proposed method was evaluated using New England Test System, and the computation accuracy and speed in determining the practical DSR were improved.

Dispersion Characteristics of the Microstrip Bend at Millimeter Wave Frequencies

The effect of the microstrip bend on the propagation characteristics of a microstrip line was investigated based on the Finite-Difference Time-Domain(FDTD)method.The dispersive characteristics of a microstrip line with a bend are quite different from that of a uniform straight microstrip line.The effect of bend discontinuity on propagation constants deceases exponentially with the distance form the corner of the bend.This can be explained by the fact that the higher order modes excited by the bend discontinuities have intrinsic properties of exponentially decay with the distance from the bend discontinuities.This effect can be negligible(less than five percent)when the distance is beyond ten times of conductor strip width.The S parameters and propagation constants comparison between unmitered and mitered bends are also presented.The proposed method which takes advantage of the FDTD algorithm and the Root Mean Square Deviation(RMSD)analysis is also an effective way for analyzing the dispersion characteristics of other planar transmission lines.

NEW DESIGN AND ANALYSIS OF FINITE DIFFERENCE DOMAIN DECOMPOSITION ALGORITHM FOR THE TWO-DIMENSIONAL HEAT EQUATION

It is presented in this paper that the new design and its analysis of finite difference domain decomposition algorithms for the two-dimensional heat equation, and the numerical results have shown the stability and accuracy of the algorithms, where SauFyev asymmetric schemes have been used at the interface points. The Algorithm II in this paper has further extended those developed by Dawson and the others, Zhang and Shen.

Joint spectrogram segmentation and ridge-extraction method for separating multimodal guided waves in long bones

Ultrasonic guided waves (GWs) can be used to evaluate long bones effectively because of the ability to provide the information of the whole bone. In this study, a joint spectrogram segmentation and ridge-extraction (JSSRE) method was proposed to separate multiple modes in long bones. First, the Gabor time-frequency transform was applied to obtain the spectrogram of multimodal signals. Then, a multi-class image segmentation algorithm was used to find the corresponding region of each mode in the spectrogram, including an improved watershed transform and a region growing procedure. Finally, the ridges were extracted and the time domain signals representing individual modes were reconstructed from these ridges in each region. The validations of this method were discussed by simulated multimodal signals with different signal-to-noise ratios (SNR). The correlation coefficients between the original signals without noise and the reconstructed signals were calculated to analyze the results quantitatively. The results showed that the extracted ridges were in good agreement with generated theoretical dispersion curves, and the reconstructed signals were highly related to the original signals, even under the SNR=3 dB situation.

Light intensification effect of trailing indent crack in fused silica subsurface

A finite-difference time-domain algorithm was applied to solve Maxwell＇s equations to obtain the redistribution of an electro- magnetic plane wave in the vicinity of a trailing indent crack （TIC）. The roles of five geometrical parameters playing in light intensification were calculated numerically under the irradiation of a 355-nm normal incidence laser. The results show that the light intensity enhancements between the nearest neighbor pits were remarkable, which may lead to damage, The calculated results reveal that the light intensity enhancement factor （LIEF） can be up to 11.2 when TIC is on the rear-surface. With the increase of the length as well as the depth of pits, LIEF increased. Conversely, with the increase of the axis of pits, LIEF grad- ually declined to a stable status. It was observed that there exists an optima width or gap, which enables LIEF to be increased dramatically and then decreased gently. By comparison, results suggest that the worst cases occur when the depth and the length are both very large, especially if the width equals to 22 and the gap equals the width. This work provides a recom- mended theoretical criterion for defect inspection and classification.