Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.

A GPU-Based Parallel Algorithm for 2D Large Deformation Contact Problems Using the Finite Particle Method

Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation from total motion in large deformation problems.In addition,the decoupled procedures of the FPM make it suitable for parallel computing,which may provide an approach to solve time-consuming issues.In this study,a graphics processing unit(GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems.The fundamentals of the FPM for planar solids are first briefly introduced,including the equations of motion of particles and the internal forces of quadrilateral elements.Subsequently,a linked-list data structure suitable for parallel processing is built,and parallel global and local search algorithms are presented for contact detection.The contact forces are then derived and directly exerted on particles.The proposed method is implemented with main solution procedures executed in parallel on a GPU.Two verification problems comprising large deformation frictional contacts are presented,and the accuracy of the proposed algorithm is validated.Furthermore,the algorithm’s performance is investigated via a large-scale contact problem,and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4,respectively,relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit(CPU).The contact calculation time percentage of the total calculation time is only 18%with the FPM,much smaller than that(50%)with Abaqus/Explicit,demonstrating the efficiency of the proposed method.

FINITE PART INTEGRALS OF LOCAL BIVARIATE C^1 QUASI-INTERPOLATING SPLINES

A local bivariate C1 quasi-interpolating spline operator with a four directional mesh is considered and studied. Based on the above operator we present cubature formulas for 2-D singular integrals, defined in the Hadamard finite part sense. Convergence results are obtained for a wide class of functions. Moreover numerical tests are given.

Form Finding and Collapse Analysis of Cable Nets Under Dynamic Loads Based on Finite Particle Method

This paper presents form finding and collapse analysis of cable net structure under strong wind using the finite particle method(FPM).As a kind of particle method,the theoretical fundamentals of the FPM are given.Methods to handle geometric and material nonlinearities of cable element are proposed.The fracture criterion and model for cable element are built to simulate the failure of cable nets.The form-finding and load analysis of two cable nets are then performed in order to initialize the successive of nonlinear analysis.The failure progress of cable nets under dynamic loads is simulated,and the dynamic responses of the typical fracture element are given in details.Analyses of the energy variations during the collapse process also show the failure mechanisms of cable nets,which is useful for the structure collapse resistance design.The numerical applications highlight the capability of the proposed procedure to solve complicate collapse problems with the FPM.

AWeighted Average Finite Difference Scheme for the Numerical Solution of Stochastic Parabolic Partial Differential Equations

In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.

A dynamic large-deformation particle finite element method for geotechnical applications based on Abaqus

In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.

Effects of layer interactions on instantaneous stability of finite Stokes flows

The stability analysis of a finite Stokes layer is of practical importance in flow control. In the present work, the instantaneous stability of a finite Stokes layer with layer interactions is studied via a linear stability analysis of the frozen phases of the base flow. The oscillations of two plates can have different velocity amplitudes, initial phases, and frequencies. The effects of the Stokes-layer interactions on the stability when two plates oscillate synchronously are analyzed. The growth rates of two most unstable modes when δ < 0.12 are almost equal, and δ = δ*/h*, where δ*and h*are the Stokes-layer thickness and the half height of the channel, respectively. However, their vorticities are different. The vorticity of the most unstable mode is symmetric, while the other is asymmetric. The Stokes-layer interactions have a destabilizing effect on the most unstable mode when δ < 0.68, and have a stabilizing effect when δ > 0.68. However, the interactions always have a stabilizing effect on the other unstable mode. It is explained that one of the two unstable modes has much higher dissipation than the other one when the Stokes-layer interactions are strong. We also find that the stability of the Stokes layer is closely related to the inflectional points of the base-flow velocity profile. The effects of inconsistent velocity-amplitude, initial phase, and frequency of the oscillations on the stability are analyzed. The energy of the most unstable eigenvector is mainly distributed near the plate of higher velocity amplitude or higher oscillation frequency. The effects of the initial phase difference are complicated because the base-flow velocity is extremely sensitive to the initial phase.

A Deep Learning Approach to Shape Optimization Problems for Flexoelectric Materials Using the Isogeometric Finite Element Method

A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.

Indentation behavior of a semi-infinite piezoelectric semiconductor under a rigid flat-ended cylindrical indenter

This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.

Profiling Electronic and Phononic Band Structures of Semiconductors at Finite Temperatures: Methods and Applications

Semiconductor devices are often operated at elevated temperatures that are well above zero Kelvin,which is the temperature in most first-principles density functional calculations.Computational approaches to com-puting and understanding the properties of semiconductors at finite temperatures are thus in critical demand.In this review,we discuss the recent progress in computationally assessing the electronic and phononic band structures of semiconductors at finite temperatures.As an emerging semiconductor with particularly strong temperature-induced renormalization of the electronic and phononic band structures,halide perovskites are used as a representative example to demonstrate how computational advances may help to understand the band struc-tures at elevated temperatures.Finally,we briefly illustrate the remaining computational challenges and outlook promising research directions that may help to guide future research in this field.

Application of the finite analytic numerical method to a flowdependent variational data assimilation

An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.

Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

Finite Element Method Simulation of Wellbore Stability under Different Operating and Geomechanical Conditions

The variation of the principal stress of formations with the working and geo-mechanical conditions can trigger wellbore instabilities and adversely affect the well completion.A finite element model,based on the theory of poro-elasticity and the Mohr-Coulomb rock damage criterion,is used here to analyze such a risk.The changes in wellbore stability before and after reservoir acidification are simulated for different pressure differences.The results indicate that the risk of wellbore instability grows with an increase in the production-pressure difference regardless of whether acidification is completed or not;the same is true for the instability area.After acidizing,the changes in the main geomechanical parameters(i.e.,elastic modulus,Poisson’s ratio,and rock strength)cause the maximum wellbore instability coefficient to increase.

Finite Deformation, Finite Strain Nonlinear Dynamics and Dynamic Bifurcation in TVE Solids with Rheology

This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.

Collapse Behavior of Pipe-Framed Greenhouses with and without Reinforcement under Snow Loading:A 3-D Finite Element Analysis

The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are considered.Three snow load distribution patterns related to the wind-driven snow particle movement are used in the analysis.It is found that snow load distribution affects the deformation and collapse behavior of the pipe-framed greenhouse significantly.The results obtained in this study are consistent with the actual damage observed.Next,discussion is made of the effects of reinforcements by adding members to the basic frame on the strength of the whole structure,in which seven kinds of reinforcement methods are examined.A buckling analysis is also carried out.The results indicate that the most effective reinforcement method depends on the snow load distribution pattern.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Study of Axisymmetric Infinite Guide Lined with Locally Reacting Material without Flow Using DtN Operators

The present work proposed a new method for the modeling by the finite element method of the acoustic propagation problems in infinite axisymmetric cylindrical guides lined with locally reacting absorbent materials without flow. The method deals with the development of an efficient transparent boundary condition based on DtN operators. The method developed in this study is successfully applied to a straight axisymmetric lined guide by imposing a mode on one of the artificial boundaries of the truncated guide. The results are in good agreement with analytical solutions. Applying the method for a non-uniform axisymmetric lined guide which is a complex case, proved its effectiveness and the results compared to those of PML layers are in very good agreement.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

On Some Results in the Theory of Finite Partially Soluble Groups

This article provides an overview of some recent results and ideas relatedto the study of finite groups depending on the restrictions on some systems of theirsections.In particular,we discuss some properties of the lattice of all subgroups ofa finite group related with conditions of permutability and generalized subnormality for subgroups.The paper contains more than 30 open problems which were posed,atdifferent times,by some mathematicians working in the discussed direction.

Texture evolution and inhomogeneous deformation of polycrystalline Cu based on crystal plasticity finite element method and particle swarm optimization algorithm

Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm optimization(PSO) algorithm. The texture-based representative volume element(TBRVE) is used in the crystal plasticity finite element model, where a given number of crystallographic orientations are obtained by means of discretizing the orientation distribution function(ODF) based on electron backscattered diffraction(EBSD) experiment data. Three-dimensional grains with different morphologies are generated on the basis of Voronoi tessellation. The PSO algorithm plays a significant role in identifying the material parameters and saving computational time. The macroscopic stress–strain curve is predicted based on CPFEM, where the simulation results are in good agreement with the experimental ones. Therefore, CPFEM is a powerful candidate for capturing the texture evolution and clarifying the inhomogeneous plastic deformation of polycrystalline Cu. The simulation results indicate that the <110> fiber texture is generated finally with the progression of plastic deformation. The inhomogeneous distribution of rotation angles lays the foundation for the inhomogeneous deformation of polycrystalline Cu in terms of grain scale.