An adaptive hierarchical domain decomposition method for parallel contact dynamics simulations of granular materials

A fully parallel version of the contact dynamics (CD) method is presented in this paper. For large enough systems, 100% efficiency has been demonstrated for up to 256 processors using a hierarchical domain decomposition with dynamic load balancing. The iterative scheme to calculate the contact forces is left domain-wise sequential, with data exchange after each iteration step, which ensures its stability. The number of additional iterations required for convergence by the partially parallel updates at the domain boundaries becomes negligible with increasing number of particles, which allows for an effective parallelization. Compared to the sequential implementation, we found no influence of the parallelization on simulation results.Comment: 19 pages, 15 figures, published in Journal of Computational Physics (2011

Parallel computational strategies for modelling transient Stokes fluid flow.

The present work is centred on two main research areas; the development of finite element techniques for the modelling of transient Stokes flow and implementation of an effective parallel system on distributed memory platforms for solving realistic large-scale Lagrangian flow problems. The first part of the dissertation presents the space-time Galerkin / least-square finite element implicit formulation for solving incompressible or slightly compressible transient Stokes flow with moving boundaries. The formulation involves a time discontinuous Galerkin method and includes least-square terms in the variational formulation. Since the additional terms involve the residual of the Euler- Lagrangian equations evaluated over element interiors, it prevents numerical oscillation on the pressure field when equal lower order interpolation functions for velocity and pressure fields are used, without violating the Babuska-Brezzi stability condition. The space-time Galerkin / least-square formulation has been successfully extended into the finite element explicit analysis, in which the penalty based discrete element contact algorithm is adopted to simulate fiuid-structure or fluid-fluid particle contact. The second part of the dissertation focuses on the development of an effective parallel processing technique, using the natural algorithm concurrency of finite element formulations. A hybrid iterative direct parallel solver is implemented into the ELFEN/implicit commercial code. The solver is based on a non-overlapping domain decomposition and sub-structure approach. The modified Cholesky factorisation is used to eliminate the unknown variables of the internal nodes at each subdomain and the resulting interfacial equations are solved by a Krylov subspace iterative method. The parallelization of explicit fluid dynamics is based on overlapping domain decomposition and a Schwarz alternating procedure. Due to the dual nature of the overlapping domain decomposition a buffer zone between any two adjacent subdomains is introduced for handling the inter-processor communication. Both solvers are tested on a PC based interconnected network system and its performances are judged by the parallel speed-up and efficiency

The effect of chemical treatment, fibre length, fibre content and injection moulding parameters to uv irradiation resistance of oil palm fibre reinforced composites

There are many types of polymer used in engineering materials expose to UV irradiation, such as automotive parts (car body, bumper, dashboard etc.) which can cause material degradation. Some polymers are used in pure polymer and some of them in composite material. This study has investigated the composite material degradation. In this study, polypropylene was used as a matrix of the composite material samples, while oil palm fibre as reinforcement. The effect of the fibre length, fibre content, fibre treatment, coupling agent and injection moulding parameter to ultraviolet (UV) light resistance of this composite and also the optimum setting of them were investigated. The UV resistance was examined via the change of mechanical properties after UV exposured in the UV accelerated weathering tester. The Linear regression models were generated for tensile strength, strain at maximum stress, break stress, break strain and Charpy Impact strength based on six different UV exposure time i.e.: 0, 96, 336, 504, 1008, and 1512 hours. The significance of the regression models were tested by Analysis Of Variance (ANOVA) and verified by two expose time i.e: 168 and 672 hours. This study found that all of the mechanical properties decrease after UV irradiated. The largest property decrease was break strain of the samples, which was decrease in the range 44.54% to 79.21% after 1512 hours UV irradiated. The lowest decrease was break stress in the range 27.38 % to 63.82%. It was also found in this study that fibre content, and UV irradiation time, significantly affect all properties. Coupling agent and alkali treatment significantly affect all properties except strain at maximum stress. Whereas fibre length and injection moulding parameter only significantly affect the Impact strength of the specimens. All of the regression models are significant which are signed by the p value of each of regression models were lower than 0.05. The equation for predicting the lifetime of UV exposed of oil palm fibre reinforced composite have been generated. The lifetime UV irradiated specimen can be predicted using this equation

Parallelization of an object-oriented FEM dynamics code: influence of the strategies on the Speedup

This paper presents an implementation in C++ of an explicit parallel finite element code dedicated to the simulation of impacts. We first present a brief overview of the kinematics and the explicit integration scheme with details concerning some particular points. Then we present the OpenMP parallelization toolkit used in order to parallelize our FEM code, and we focus on how the parallelization of the DynELA FEM code has been conducted for a shared memory system using OpenMP. Some examples are then presented to demonstrate the efficiency and accuracy of the proposed implementations concerning the Speedup of the code. Finally, an impact simulation application is presented and results are compared with the ones obtained by the commercial Abaqus explicit FEM code

FullSWOF_Paral: Comparison of two parallelization strategies (MPI and SKELGIS) on a software designed for hydrology applications

In this paper, we perform a comparison of two approaches for the parallelization of an existing, free software, FullSWOF 2D (http://www. univ-orleans.fr/mapmo/soft/FullSWOF/ that solves shallow water equations for applications in hydrology) based on a domain decomposition strategy. The first approach is based on the classical MPI library while the second approach uses Parallel Algorithmic Skeletons and more precisely a library named SkelGIS (Skeletons for Geographical Information Systems). The first results presented in this article show that the two approaches are similar in terms of performance and scalability. The two implementation strategies are however very different and we discuss the advantages of each one.Comment: 27 page

Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model

We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for the space discretization, our implementation provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner (LMW) benchmark [1]

A multidomain spectral method for solving elliptic equations

We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three distinct features. First, the combined problem of solving the PDE, satisfying the boundary conditions, and matching between different subdomains is cast into one set of equations readily accessible to standard linear and nonlinear solvers. Second, touching as well as overlapping subdomains are supported; both rectangular blocks with Chebyshev basis functions as well as spherical shells with an expansion in spherical harmonics are implemented. Third, the code is very flexible: The domain decomposition as well as the distribution of collocation points in each domain can be chosen at run time, and the solver is easily adaptable to new PDEs. The code has been used to solve the equations of the initial value problem of general relativity and should be useful in many other problems. We compare the new method to finite difference codes and find it superior in both runtime and accuracy, at least for the smooth problems considered here.Comment: 31 pages, 8 figure