Algorithm of far-field centre estimation based on phase-only matched filter

Estimation of the far-field centre is carried out in beam auto-alignment. In this paper, the features of the far-field of a square beam are presented. Based on these features, a phase-only matched filter is designed, and the algorithm of centre estimation is developed. Using the simulated images with different kinds of noise and the 40 test images that are taken in sequence, the accuracy of this algorithm is estimated. Results show that the error is no more than one pixel for simulated noise images with a 99% probability, and the stability is restricted within one pixel for test images. Using the improved algorithm, the consumed time is reduced to 0.049 s.

Ground target detection, classification and sensor fusion in distributed fiber seismic sensor network - art. no. 683015

This paper describes the ground target detection, classification and sensor fusion problems in distributed fiber seismic sensor network. Compared with conventional piezoelectric seismic sensor used in UGS, fiber optic sensor has advantages of high sensitivity and resistance to electromagnetic disturbance. We have developed a fiber seismic sensor network for target detection and classification. However, ground target recognition based on seismic sensor is a very challenging problem because of the non-stationary characteristic of seismic signal and complicated real life application environment. To solve these difficulties, we study robust feature extraction and classification algorithms adapted to fiber sensor network. An united multi-feature (UMF) method is used. An adaptive threshold detection algorithm is proposed to minimize the false alarm rate. Three kinds of targets comprise personnel, wheeled vehicle and tracked vehicle are concerned in the system. The classification simulation result shows that the SVM classifier outperforms the GMM and BPNN. The sensor fusion method based on D-S evidence theory is discussed to fully utilize information of fiber sensor array and improve overall performance of the system. A field experiment is organized to test the performance of fiber sensor network and gather real signal of targets for classification testing.

Application of the Conjugate Gradient Optimization Algorithm of Neural Network on Laser Machining

在应用激光技术加工复杂曲面时,通常以采样点集为插值点来建立曲面函数,然后实现曲面上任意坐标点的精确定位。人工神经网络的BP算法能实现函数插值,但计算精度偏低,往往达不到插值精确要求,造成较大的加工误差。提出人工神经网络的共轭梯度最优化插值新算法,并通过实例仿真,证明了这种曲面精确定位方法的可行性,从而为激光加工的三维精确定位提供了一种良好解决方案。这种方法已经应用在实际中。

Simulating shock to detonation transition: Algorithm and results

An algorithm based on flux-corrected transport and the Lagrangian finite element method is presented for solving the problem of shock dynamics. It is verified through the model problem of one-dimensional strain elastoplastic shock wave propagation that the algorithm leads to stable, non-oscillatory results. Shock initiation and detonation wave propagation is simulated using the algorithm, and some interesting results are obtained. (C) 1999 Academic Press.

Ancient Chinese algorithm: The Ying Buzu Shu (method of surplus and deficiency) vs Newton iteration method

Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of double false position is actually a translation version of the ancient Chinese algorithm, a comparison with well-known Newton iteration method is also made. If derivative is introduced, the ancient Chinese algorithm reduces to the Newton method. A modification of the ancient Chinese algorithm is also proposed, and some of applications to nonlinear oscillators are illustrated.

Composite materials dynamic fracture studies by generalized Shmuely difference algorithm

The generalized Shmuely Difference Algorithm (GSDA) is presented here to analyze the dynamic fracture performance of orthogonal-anisotropic composite materials, such as glass fibre reinforced phenolplast. The difference recurrence Formulae and boundary condition difference extrapolation formulae are derived and programmed. The dynamic stress intensity factors (DSIF) of the isotropic and anisotropic centrally cracked plates are computed respectively using GSDA and compared with that published previously. GSDA is proved effective and reliable. Copyright (C) 1996 Elsevier Science Ltd.

Optimization of 2-stage gas gun for fusion refueling pellet injection with consideration of real effects

Analytic expression of pellet acceleration by constant base pressure with consideration of gas-wall friction, heat transfer and viscous dissipation that important for high speed injection is obtained. The process of compression stage is formulated by a set of governing equations and can be numerically integrated. Excellent confirmation with experiments is obtained and the ways to optimum match the compression stage with the launch stage are suggested.

Fusion energy-production from a deuterium-tritium plasma in the jet tokamak

Describes a series of experiments in the Joint European Torus (JET), culminating in the first tokamak discharges in deuterium-tritium fuelled mixture. The experiments were undertaken within limits imposed by restrictions on vessel activation and tritium usage. The objectives were: (i) to produce more than one megawatt of fusion power in a controlled way; (ii) to validate transport codes and provide a basis for accurately predicting the performance of deuterium-tritium plasmas from measurements made in deuterium plasmas; (iii) to determine tritium retention in the torus systems and to establish the effectiveness of discharge cleaning techniques for tritium removal; (iv) to demonstrate the technology related to tritium usage; and (v) to establish safe procedures for handling tritium in compliance with the regulatory requirements. A single-null X-point magnetic configuration, diverted onto the upper carbon target, with reversed toroidal magnetic field was chosen. Deuterium plasmas were heated by high power, long duration deuterium neutral beams from fourteen sources and fuelled also by up to two neutral beam sources injecting tritium. The results from three of these high performance hot ion H-mode discharges are described: a high performance pure deuterium discharge; a deuterium-tritium discharge with a 1% mixture of tritium fed to one neutral beam source; and a deuterium-tritium discharge with 100% tritium fed to two neutral beam sources. The TRANSP code was used to check the internal consistency of the measured data and to determine the origin of the measured neutron fluxes. In the best deuterium-tritium discharge, the tritium concentration was about 11% at the time of peak performance, when the total neutron emission rate was 6.0 × 10¹⁷ neutrons/s. The integrated total neutron yield over the high power phase, which lasted about 2 s, was 7.2 × 10¹⁷ neutrons, with an accuracy of ±7%. The actual fusion amplification factor, Q_DT was about 0.15

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

Algorithm for simulating fracture processes in concrete by lattice modeling

To simulate fracture behaviors in concrete more realistically, a theoretical analysis on the potential question in the quasi-static method is presented, then a novel algorithm is proposed which takes into account the inertia effect due to unstable crack propagation and meanwhile requests much lower computational efforts than purely dynamic method. The inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, results may become questionable if a fracture process including unstable cracking is simulated by the quasi-static procedure excluding completely inertia effects. However, it requires much higher computational effort to simulate experiments with not very high loading rates by the dynamic method. In this investigation which can be taken as a natural continuation, the potential question of quasi-static method is analyzed based on the dynamic equations of motion. One solution to this question is the new algorithm mentioned above. Numerical examples are provided by the generalized beam (GB) lattice model to show both fracture processes under different loading rates and capability of the new algorithm.

An Integrated Evolutionary Algorithm for Expensive Global Optimization

We propose an integrated algorithm named low dimensional simplex evolution extension (LDSEE) for expensive global optimization in which only a very limited number of function evaluations is allowed. The new algorithm accelerates an existing global optimization, low dimensional simplex evolution (LDSE), by using radial basis function (RBF) interpolation and tabu search. Different from other expensive global optimization methods, LDSEE integrates the RBF interpolation and tabu search with the LDSE algorithm rather than just calling existing global optimization algorithms as subroutines. As a result, it can keep a good balance between the model approximation and the global search. Meanwhile it is self-contained. It does not rely on other GO algorithms and is very easy to use. Numerical results show that it is a competitive alternative for expensive global optimization.

A block SPP algorithm for multidimensional tridiagonal equations with optimal message vector length

A parallel strategy for solving multidimensional tridiagonal equations is investigated in this paper. We present in detail an improved version of single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication cost. We show the resulting block SPP can achieve good speedup for a wide range of message vector length (MVL), especially when the number of grid points in the divided direction is large. Instead of only using the largest possible MVL, we adopt numerical tests and modeling analysis to determine an optimal MVL so that significant improvement in speedup can be obtained.

On optimal message vector length for block single parallel partition algorithm in a three-dimensional ADI solver

It has long been recognized that many direct parallel tridiagonal solvers are only efficient for solving a single tridiagonal equation of large sizes, and they become inefficient when naively used in a three-dimensional ADI solver. In order to improve the parallel efficiency of an ADI solver using a direct parallel solver, we implement the single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication costs. The measured performances show that the longest allowable message vector length (MVL) is not necessarily the best choice. To understand this observation and optimize the performance, we propose an improved model that takes the cache effect into consideration. The optimal MVL for achieving the best performance is shown to depend on number of processors and grid sizes. Similar dependence of the optimal MVL is also found for the popular block pipelined method.

Laser-confined fusion

An approach for producing a large quantity of neutrons is proposed. It involves compression of a fuel foil and confinement of the resulting plasma between two intense laser pulses. It is shown that two circularly polarized laser pulses of amplitude a=7 illuminating a deuterium-tritium foil of areal density 3.3 X 10(18) cm(-2) can produce about 4.2 X 10(6) neutrons per joule of the input laser energy.