Dynamic Parallel Method for Eulerian Method of Elastoplastic Hydrodynamics in Three Dimension

Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such as detonation and penetration, the dynamic parallel method （DPM） is designed to adjust the computational domain dynamically to get better load balance. Dynamic parallel method can be separated into two parts： one is division of initial computational domain and location of the data, the other is expansion of the computational domain and adjustment of the data location. DPM program can greatly shorten computational time and be preferable in simulating actual problems. The speedup of the DPM program is linear in parallel test. DPM can be popularized to parallel program of other multi-component high dimension Eulerian methods naturally.

A parallel method for numerical solution of delay differential equations

A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.

An Approximate Time-Parallel Method for the Fast and Accurate Computation of Particle Trajectories in a Magnetic Field

A temporal multiscale hybridization method is presented that carefully couples coarse scale gyrokinetic models with exact charged particle solution trajectories (that is, with full phase information) in a magnetic field. The approach is based on the careful approximation of a sum, generally employed for time-parallel (TP) computing applications. While the hybridization method presented is highly parallelizable, a computational efficiency gain is seen from considering serial computations only. A complete numerical method is only presented for the aforementioned charged particle application, however, the general approach depicted likely has relevance to a wide swath of challenging multiscale/multiphysics problems. Additionally, the approach has obvious relevance to TP computing applications (such as variable selection on which to perform TP calculations and fine scale sampling strategies).

Fast Parallel Method for Polynomial Evaluation at Points in Arithmetic Progression

We present a fast method for polynomial evaluation at points in arithmetic progression. By dividing the progression into m new ones and evaluating the polynomial at each point of these new progressions recursively,this method saves most of the multiplications in the price of little increase of additions comparing to Horner's method, while their accuracy are almost the same. We also introduce vector structure to the recursive process making it suitable for parallel applications.

PARALLEL INTERVAL MATRIX MULTISPLITTING AOR METHODS AND THEIR CONVERGENCE

This paper proposes a class of parallel interval matrix multisplitting AOR methods far solving systems of interval linear equations and discusses their convergence properties under the conditions that the coefficient matrices are interval H-matrices.

ON THE CONVERGENCE OF PARALLEL BFGS METHOD

According to the sequential BFGS method, in this paper we present an asynchronous parallel BFGS method in the case when the gradient information about the function is inexact. We assume that we have p + q processors, which are divided-into two groups, the first group has p processors, the second group has q processors, the two groups are asynchronous. parallel, If we assume the objective function is twice continuously differentiable and uniformly convex, we prove the iteration converge globally to the solution, and under some additional conditions we show the method is superlinearly convergent. Finally, we show the numerical results of this algorithm.

PARALLEL REGION PRESERVING MULTISECTION METHOD FOR SOLVING GENERALIZED EIGENPROBLEM

The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of symmetric tridiagonal matrix. The multisection method for solving the generalized eigenproblem applied significantly in many science and engineering domains has not been studied. The parallel region preserving multisection method (PRM for short) for solving generalized eigenproblems of large sparse and real symmetric matrix is presented in this paper. This method not only retains the advantages of the conventional determinant search method (DS for short), but also overcomes its disadvantages such as leaking roots and disconvergence. We have tested the method on the YH 1 vector computer, and compared it with the parallel region preserving determinant search method the parallel region preserving bisection method (PRB for short). The numerical results show that PRM has a higher speed up, for instance, it attains the speed up of 7.7 when the scale of the problem is 2 114 and the eigenpair found is 3, and PRM is superior to PRB when the scale of the problem is large.

LOCAL AND PARALLEL FINITE ELEMENT METHOD FOR THE MIXED NAVIER-STOKES/DARCY MODEL WITH BEAVERS-JOSEPH INTERFACE CONDITIONS

In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.

Coupling analysis of transmission lines excited by space electromagnetic fields based on time domain hybrid method using parallel technique

We present a time domain hybrid method to realize the fast coupling analysis of transmission lines excited by space electromagnetic fields, in which parallel finite-difference time-domain (FDTD) method, interpolation scheme, and Agrawal model-based transmission line (TL) equations are organically integrated together. Specifically, the Agrawal model is employed to establish the TL equations to describe the coupling effects of space electromagnetic fields on transmission lines. Then, the excitation fields functioning as distribution sources in TL equations are calculated by the parallel FDTD method through using the message passing interface (MPI) library scheme and interpolation scheme. Finally, the TL equations are discretized by the central difference scheme of FDTD and assigned to multiple processors to obtain the transient responses on the terminal loads of these lines. The significant feature of the presented method is embodied in its parallel and synchronous calculations of the space electromagnetic fields and transient responses on the lines. Numerical simulations of ambient wave acting on multi-conductor transmission lines (MTLs), which are located on the PEC ground and in the shielded cavity respectively, are implemented to verify the accuracy and efficiency of the presented method.

Global SH-wave propagation in a 2D whole Moon model using the parallel hybrid PSM/FDM method

We present numerical modeling of SH-wave propagation for the recently proposed whole Moon model and try to improve our understanding of lunar seismic wave propagation. We use a hybrid PSM/FDM method on staggered grids to solve the wave equations and implement the calculation on a parallel PC cluster to improve the computing efficiency. Features of global SH-wave propagation are firstly discussed for a 100-km shallow and900-km deep moonquakes, respectively. Effects of frequency range and lateral variation of crust thickness are then investigated with various models. Our synthetic waveforms are finally compared with observed Apollo data to show the features of wave propagation that were produced by our model and those not reproduced by our models. Our numerical modeling show that the low-velocity upper crust plays significant role in the development of reverberating wave trains. Increasing frequency enhances the strength and duration of the reverberations.Surface multiples dominate wavefields for shallow event.Core–mantle reflections can be clearly identified for deep event at low frequency. The layered whole Moon model and the low-velocity upper crust produce the reverberating wave trains following each phases consistent with observation. However, more realistic Moon model should be considered in order to explain the strong and slow decay scattering between various phases shown on observation data.

Estimates of Convergence Rate of Parallel Multisplitting Itertive Methods

This paper givers an estimated formula of convergence rate for parallel multisplitting iterative method.Using the formula,we can simplify and unify the proof of convergence of PMI_method.

Space decomposition based parallelization solutions for the combined finiteediscrete element method in 2D

The combined finiteediscrete element method （FDEM） belongs to a family of methods of computationalmechanics of discontinua. The method is suitable for problems of discontinua, where particles aredeformable and can fracture or fragment. The applications of FDEM have spread over a number of disciplinesincluding rock mechanics, where problems like mining, mineral processing or rock blasting canbe solved by employing FDEM. In this work, a novel approach for the parallelization of two-dimensional（2D） FDEM aiming at clusters and desktop computers is developed. Dynamic domain decompositionbased parallelization solvers covering all aspects of FDEM have been developed. These have beenimplemented into the open source Y2D software package and have been tested on a PC cluster. Theoverall performance and scalability of the parallel code have been studied using numerical examples. Theresults obtained confirm the suitability of the parallel implementation for solving large scale problems. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

CONTINUATION METHOD APPLIED IN KINEMATICS OF PARALLEL ROBOT

Continuation method solving forward kinematics problem of parallel robot was discussed. And through a coefficient-parameter continuation method the efficiency and feasibility of continuation method were improved. Using this method all forward solutions of a new parallel robot model which was put forward lately by Robot Open Laboratory of Science Institute of China were obtained. Therefore it provided the basis of mechanism analysis and real-time control for new model.

A Class of Real-Time Parallel Combined Methods ofDigital Simulation for Large Systems

In this paper, a class of real-time parallel combined methods (RTPCM) of the digital simulation for a partitioned large system is presented. By means of combination of the parallelism across the system with the parallelism across the method, stiff and non-stiff subsystems are solved in parallel on parallel computer by a parallel Rosenbrock method and a parallel RK method, respectively. Their construction, convergence and numerical stability are discussed, and the digitalsimulation experiments are conducted.

Class of modified parallel combined methods of real-time numerical simulation for a stiff system

A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.

The Model of Asynchronous Parallel Nonlinear Multisplitting Method on Shared Memory System

Nonlinear multisplitting method is known as parallel iterative methods for solving a large-scale system of nonlinear equations F(x) = 0. We extend the idea of nonlinear multisplitting and consider a new model ill which the iteration is executed asynchronously: Each processor calculate the solution of an individual nonlinear system belong to its nonlinear multisplitting and can update the global approximation residing in the shared memory at any time. A local convergence analysis of this model is presented. Finally, we give a uumerical example which shows a 'strange' property that speedup Sp > p and efficiency Ep > 1.

THE PARALLEL MULTISPLITTING METHOD FOR CONSISTENT SYMMETRIC POSITIVE(SEMI-)DEFINITE SYSTEMS

Main resultsTheorem 1 Let A be an n×n symmetric positive semidefinite matrix and let

Review of the algebraic linear methods and parallel implementation in numerical simulation of groundwater flow

The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.

Stability analysis method considering non-parallelism:EPSE method and its application

The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory （LST） based on the local parallel hypothesis. Considering the non-parallelism effect, the parabolized stability equation （PSE） method lacks local characteristic of stability analysis. In this paper, a local stability analysis method considering non-parallelism is proposed, termed as EPSE since it may be considered as an expansion of the PSE method. The EPSE considers variation of the shape function in the streamwise direction. Its local characteristic is convenient for stability analysis. This paper uses the EPSE in a strong non-parallel flow and mode exchange problem. The results agree well with the PSE and the direct numerical simulation （DNS）. In addition, it is found that the growth rate is related to the normalized method in the non-parallel flow. Different results can be obtained using different normalized methods. Therefore, the normalized method must be consistent.

Type Synthesis of 4-DOF Parallel Kinematic Mechanisms Based on Grassmann Line Geometry and Atlas Method

Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms（PKMs）, but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom（DOFs） and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R（T： translational DOF; R： rotational DOF）, are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs＇ type synthesis. Moreover, this method can provide theoretical guidance for other PKMs＇ type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods＇ problems in terms of complicated, not intuitive and unsuitable for practical application.