Abstract:Purpose - The purpose of this paper is to propose an efficient iterative method for large-scale finite element equations of bad numerical stability arising from deformation analysis with multi-point constraint using Lagrange multiplier method. Design/methodology/approach - In this paper, taking warpage analysis of polymer injection molding based on surface model as an example, the performance of several popular Krylov subspace methods, including conjugate gradient, BiCGSTAB and generalized minimal residual (GMRES), with diffident Incomplete LU (ILU)-type preconditions is investigated and compared. For controlling memory usage, GMRES(m) is also considered. And the ordering technique, commonly used in the direct method, is introduced into the presented iterative method to improve the preconditioner. Findings - It is found that the proposed preconditioned GMRES method is robust and effective for solving problems considered in this paper, and approximate minimum degree (AMD) ordering is most beneficial for the reduction of fill-ins in the ILU preconditioner and acceleration of the convergence, especially for relatively accurate ILU-type preconditioning. And because of concerns about memory usage, GMRES(m) is a good choice if necessary. Originality/value - In this paper, for overcoming difficulties of bad numerical stability resulting from Lagrange multiplier method, together with increasing scale of problems in engineering applications and limited hardware conditions of computer, a stable and efficient preconditioned iterative method is proposed for practical purpose. Before the preconditioning, AMD reordering, commonly used in the direct method, is introduced to improve the preconditioner. The numerical experiments show the good performance of the proposed iterative method for practical cases, which is implemented in in-house and commercial codes on PC.

Acceleration Technique Using Krylov Subspace Methods for 2d Arbitrary Geometry Characteristics Solver

An acceleration technique for 2D MOC based on Krylov subspace and domain decomposition methods

Two-Dimensional Matrix Characteristics Method Based on Krylov Subspace and Domain Decomposition Theories

Acceleration and Parallelization of Arbitrary Trajectory Based Three-Dimensional Method of Characteristics

Comparative study of inner-outer Krylov solvers for linear systems in structured and high-order unstructured CFD problems

GMRES Implicit Algorithm Based on 2D Unstructured Meshes for Solving Euler Equations

The research of discontinuous galerkin p-multigrid solver

The Simpler GMRES Method Combined with Finite Volume Method for Simulating Viscoelastic Flows on Triangular Grid.

An efficient preconditioned Krylov subspace method for large-scale finite element equations with MPC using Lagrange multiplier method

Application of Discontinuous Galerkin Method in Solving Euler's Equations

Generalized coarse-mesh finite difference acceleration for the method of characteristics

Anderson Acceleration for Geometry Optimization and Physics Simulation

A dynamically deflated GMRES adjoint solver for aerodynamic shape optimization

Effectively Applying High-Order Discontinuous Galerkin Method(DGM) to Solving 3-D Euler Equations on Unstructured Grids

A novel family of iterative solvers for method of moments discretizations of maxwell's equations

Achieving high performance and portable parallel GMRES algorithm for compressible flow simulations on unstructured grids

Accurate and Robust Cfd Algorithms Applied to 3d Arbitrary Polyhedral Grids

An energy-efficient GMRES–multigrid solver for space-time finite element computation of dynamic poroelasticity

Fast block Krylov subspace methods for solving sequences of dense MoM linear systems with multiple right-hand sides

Implementation of a Newton−Krylov Algorithm in the Open-source Solver of PHengLEI

Accelerated Direct Solver With Multiscale Compressed Block Decomposition and Multilevel Characteristic Basis Function Method