Load balancing for aspect ratio in adaptive finite element simulations

In parallel adaptive finite element simulations the work load on the individual processors may change frequently. To (re)distribute the load evenly over the processors a load balancing heuristic is needed. Common strategies try to minimise subdomain dependencies by optimising the cutsize of the partitioning. However for certain solvers cutsize only plays a minor role, and their convergence is highly dependent on the subdomain shapes. Degenerated subdomain shapes cause them to need significantly more iterations to converge. In this work a new parallel load balancing strategy is introduced which directly addresses the problem of generating and conserving reasonably good subdomain shapes in a dynamically changing Finite Element Simulation. Geometric data is used to formulate several cost functions to rate elements in terms of their suitability to be migrated. The well known diffusive method which calculates the necessary load flow is enhanced by weighting the subdomain edges with the help of these cost functions. The proposed methods have been tested and results are presented.

Dynamic multi-partitioning for parallel finite element applications

The central product of the DRAMA (Dynamic Re-Allocation of Meshes for parallel Finite Element Applications) project is a library comprising a variety of tools for dynamic re-partitioning of unstructured Finite Element (FE) applications. The input to the DRAMA library is the computational mesh, and corresponding costs, partitioned into sub-domains. The core library functions then perform a parallel computation of a mesh re-allocation that will re-balance the costs based on the DRAMA cost model. We discuss the basic features of this cost model, which allows a general approach to load identification, modelling and imbalance minimisation. Results from crash simulations are presented which show the necessity for multi-phase/multi-constraint partitioning components.

A Posteriori Error Analysis of hp-Version Discontinuous Galerkin Finite Element Methods for Second-Order Quasilinear Elliptic Problems

We develop the a-posteriori error analysis of hp-version interior-penalty discontinuous Galerkin finite element methods for a class of second-order quasilinear elliptic partial differential equations. Computable upper and lower bounds on the error are derived in terms of a natural (mesh-dependent) energy norm. The bounds are explicit in the local mesh size and the local degree of the approximating polynomial. The performance of the proposed estimators within an automatic hp-adaptive refinement procedure is studied through numerical experiments.

Appraisal of open software for finite element simulation of 2D metal sheet laser cut.

FEA simulation of thermal metal cutting is central to interactive design and manufacturing. It is therefore relevant to assess the applicability of FEA open software to simulate 2D heat transfer in metal sheet laser cuts. Application of open source code (e.g. FreeFem++, FEniCS, MOOSE) makes possible additional scenarios (e.g. parallel, CUDA, etc.), with lower costs. However, a precise assessment is required on the scenarios in which open software can be a sound alternative to a commercial one. This article contributes in this regard, by presenting a comparison of the aforementioned freeware FEM software for the simulation of heat transfer in thin (i.e. 2D) sheets, subject to a gliding laser point source. We use the commercial ABAQUS software as the reference to compare such open software. A convective linear thin sheet heat transfer model, with and without material removal is used. This article does not intend a full design of computer experiments. Our partial assessment shows that the thin sheet approximation turns to be adequate in terms of the relative error for linear alumina sheets. Under mesh resolutions better than 10e−5 m , the open and reference software temperature differ in at most 1 % of the temperature prediction. Ongoing work includes adaptive re-meshing, nonlinearities, sheet stress analysis and Mach (also called ‘relativistic’) effects.

A class of domain decomposition preconditioners for hp-discontinuous Galerkin finite element methods

In this article we address the question of efficiently solving the algebraic linear system of equations arising from the discretization of a symmetric, elliptic boundary value problem using hp-version discontinuous Galerkin finite element methods. In particular, we introduce a class of domain decomposition preconditioners based on the Schwarz framework, and prove bounds on the condition number of the resulting iteration operators. Numerical results confirming the theoretical estimates are also presented.

High-order hp-adaptive discontinuous Galerkin finite element methods for compressible fluid flows

This article is concerned with the construction of general isotropic and anisotropic adaptive strategies, as well as hp-mesh refinement techniques, in combination with dual-weighted-residual a posteriori error indicators for the discontinuous Galerkin finite element discretization of compressible fluid flow problems.

Slope stability analysis for fill slopes using finite element limit analysis

This paper investigates the stability of fill slopes often found in embankment cases where frictional fill materials are placed on purely cohesive undrained clay with increasing strength. By using finite element upper and lower bound limit analysis for this investigation, the limit load can be truly bounded. It is known that two-dimensional analysis yields a more conservative result due to plain strain condition when compared to three-dimensional analysis. Therefore, this paper will focus on three-dimensional (3D) slope stability analysis and for comparison purposes two-dimensional analysis results will be employed. In fact, the final results are presented in the form of comprehensive chart solutions for the convenience of practicing engineers during preliminary slope design. The failure mechanism will also be discussed in order to further illustrate the situation during failure. It should be highlighted that the failure mechanisms are obtained through the numerical method itself and no prior assumptions are required, therefore, are more realistic and able to provide a better understanding for the slope failure surfaces.

A study of guided wave propagation in timber pole using spectral finite element method

Timber is one of the most widely used structural material all over the world. Round timbers can be seen as a structural component in historical buildings, jetties, short span bridges and also as piles for foundation and poles for electrical and power distribution. To evaluate the current condition of these cylindrical type timber structures, guided wave has a great potential. However, the difficulties associated with the guided wave propagation in timber materials includes orthotropic behaviour of wood, moisture contents, temperature, grain direction, etc. In addition, the effect of fully or partially filled surrounding media, such as soil, water, etc. causes attenuation on the generated stress wave. In order to investigate the effects of these parameters on guided wave propagation, extensive numerical simulation is required to conduct parametric studies. Moreover, due to the presence of multi modes in guided wave propagation, dispersion curves are of great importance. Even though conventional finite element method (FEM) can determine dispersion curves along with wave propagation in time domain, it is highly computationally expensive. Furthermore, incorporating orthotropic behaviour and surrounding media to model a thick cylindrical wave (large diameter cylindrical structures) make conventional FEM inefficient for this purpose. In contrast, spectral finite element method (SFEM) is a semi analytical method to model the guided wave propagation which does not need fine meshes compared to the other methods, such as FEM or finite difference method (FDM). Also, even distribution of mass and stiffness of structures can be obtained with very few elements using SFEM. In this paper, the suitability of SFEM is investigated to model guided wave propagation through an orthotropic cylindrical waveguide with the presence of surrounding soil. Both the frequency domain analysis (dispersion curves) and time domain reconstruction for a multi-mode generated input signal are presented under different loading location. The dispersion curves obtained from SFEM are compared against analytical solution to verify its accuracy. Lastly, different numerical issues to solve for the dispersion curves and time domain results using SFEM are also discussed.

SIMULATING LARGE DEFORMATION OF INTRANEURAL GANGLION CYST USING FINITE ELEMENT METHOD

Intraneural Ganglion Cyst is disorder observed in the nerve injury, it is still unknown and very difficult to predict its propagation in the human body so many times it is referred as an unsolved history. The treatments for this disorder are to remove the cystic substance from the nerve by a surgery. However these treatments may result in neuropathic pain and recurrence of the cyst. The articular theory proposed by Spinner et al., (Spinner et al. 2003) considers the neurological deficit in Common Peroneal Nerve (CPN) branch of the sciatic nerve and adds that in addition to the treatment, ligation of articular branch results into foolproof eradication of the deficit. Mechanical modeling of the affected nerve cross section will reinforce the articular theory (Spinner et al. 2003). As the cyst propagates, it compresses the neighboring fascicles and the nerve cross section appears like a signet ring. Hence, in order to mechanically model the affected nerve cross section; computational methods capable of modeling excessively large deformations are required. Traditional FEM produces distorted elements while modeling such deformations, resulting into inaccuracies and premature termination of the analysis. The methods described in research report have the capability to simulate large deformation. The results obtained from this research shows significant deformation as compared to the deformation observed in the conventional finite element models. The report elaborates the neurological deficit followed by detail explanation of the Smoothed Particle Hydrodynamic approach. Finally, the results show the large deformation in stages and also the successful implementation of the SPH method for the large deformation of the biological organ like the Intra-neural ganglion cyst.

Development of a Coupling Model for Fluid-Structure Interaction using the Mesh-free Finite Element Method and the Lattice Boltzmann Method

In the presented thesis work, the meshfree method with distance fields was coupled with the lattice Boltzmann method to obtain solutions of fluid-structure interaction problems. The thesis work involved development and implementation of numerical algorithms, data structure, and software. Numerical and computational properties of the coupling algorithm combining the meshfree method with distance fields and the lattice Boltzmann method were investigated. Convergence and accuracy of the methodology was validated by analytical solutions. The research was focused on fluid-structure interaction solutions in complex, mesh-resistant domains as both the lattice Boltzmann method and the meshfree method with distance fields are particularly adept in these situations. Furthermore, the fluid solution provided by the lattice Boltzmann method is massively scalable, allowing extensive use of cutting edge parallel computing resources to accelerate this phase of the solution process. The meshfree method with distance fields allows for exact satisfaction of boundary conditions making it possible to exactly capture the effects of the fluid field on the solid structure.

Characterisation and investigation of local variations in mechanical behaviour in cast aluminium using gradient solidification, Digital Image Correlation and finite element simulation

Due to design and process-related factors, there are local variations in the microstructure and mechanical behaviour of cast components. This work establishes a Digital Image Correlation (DIC) based method for characterisation and investigation of the effects of such local variations on the behaviour of a high pressure, die cast (HPDC) aluminium alloy. Plastic behaviour is studied using gradient solidified samples and characterisation models for the parameters of the Hollomon equation are developed, based on microstructural refinement. Samples with controlled microstructural variations are produced and the observed DIC strain field is compared with Finite Element Method (FEM) simulation results. The results show that the DIC based method can be applied to characterise local mechanical behaviour with high accuracy. The microstructural variations are observed to cause a redistribution of strain during tensile loading. This redistribution of strain can be predicted in the FEM simulation by incorporating local mechanical behaviour using the developed characterization model. A homogeneous FEM simulation is unable to predict the observed behaviour. The results motivate the application of a previously proposed simulation strategy, which is able to predict and incorporate local variations in mechanical behaviour into FEM simulations already in the design process for cast components.

Finite element modeling of tube deformation during cold pilgering

A three-dimensional finite element model of cold pilgering of stainless steel tubes is developed in this paper. The objective is to use the model to increase the understanding of forces and deformations in the process. The focus is on the influence of vertical displacements of the roll stand and axial displacements of the mandrel and tube. Therefore, the rigid tools and the tube are supported with elastic springs. Additionally, the influences of friction coefficients in the tube/mandrel and tube/roll interfaces are examined. A sensitivity study is performed to investigate the influences of these parameters on the strain path and the roll separation force. The results show the importance of accounting for the displacements of the tube and rigid tools on the roll separation force and the accumulative plastic strain.

Experimental and finite element analysis of machinability of AL-6XN super austenitic stainless steel

This paper presents a finite element cutting modelbased on physical microstructure to investigate the thermomechanicalbehaviour of AL-6XN Super AusteniticStainless Steel in the primary shear zone. Frozen chip rootsamples were created under dry turning operation to observethe plasticity behaviour occurring in the shear zones to comparewith the model for analysis. Chip samples were generatedunder cutting velocities at 65 and 94 m/min, feed rate at0.2 mm/rev and depth of cut at 1 mm. Temperature on thecutting zone was recorded by infrared thermal camera.Secondary and backscatter electron detectors were used toinvestigate the deformed microstructure and to calculate theplastic strain. Experimental results showed the formation ofmicrocracks (build-up edge triggers) at the chip root stagnationzone of both samples. The austenite phase patterns wereevident against the cutting tool tip in the stagnation zone of thechip root fabricated at 65 m/min. The movement of thesepatterns caused the formation of the slip lines within thegrains. The backscatter diffraction maps showed the formationof special grain boundaries within the slip lines, workhardeninglayer and in the chip region. Strain measurementsin the microstructures of the chip roots fabricated at 94 and65 m/min showed high values of 6.5 and 5.7 (mm/mm) respectively.The finite element model was used to measure thestress, strain, temperature and chip morphology. Numericalresults were compared to the outcomes of the experimentalwork to validate the finite element model. The model validatingprocess showed good agreement between theexperimental and numerical results, and the error values werecalculated. For a 94- and 65-m/min cutting speeds, 7.5 and5.2% were the errors in the strain, 3 and 2.5% were the error inthe temperature and 4.7 and 6.8% were the error in the shearplane angles.