Abstract

In this paper a stabilized finite element method to deal with incompressibility in solid mechanics is presented. A mixed formulation involving pressure and displacement fields is used and a continuous linear interpolation is considered for both fields. To overcome the Babuška–Brezzi condition, a stabilization technique based on the orthogonal sub-scale method is introduced. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate for real industrial applications. Results are compared with standard Galerkin and Q1P0 mixed formulations for nearly incompressible problems in the context of linear elasticity.

Abstract

The use of stabilization methods is becoming an increasingly well-accepted technique due to their success in dealing with numerous numerical pathologies that arise in a variety of applications in computational mechanics.

In this paper a multiscale finite element method technique to deal with pressure stabilization of nearly incompressibility problems in nonlinear solid mechanics at finite deformations is presented. A J2-flow theory plasticity model at finite deformations is considered. A mixed formulation involving pressure and displacement fields is used as starting point. Within the finite element discretization setting, continuous linear interpolation for both fields is considered. To overcome the Babuška–Brezzi stability condition, a multiscale stabilization method based on the orthogonal subgrid scale (OSGS) technique is introduced. A suitable nonlinear expression of the stabilization parameter is proposed. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate and, therefore, very convenient for practical industrial applications.

Numerical results obtained using the OSGS stabilization technique are compared with results provided by the P1 standard Galerkin displacements linear triangular/tetrahedral element, P1/P1 standard mixed linear displacements/linear pressure triangular/tetrahedral element and Q1/P0 mixed bilinear/trilinear displacements/constant pressure quadrilateral/hexahedral element for 2D/3D nearly incompressible problems in the context of a nonlinear finite deformation J2 plasticity model.

Abstract

This paper deals with the computational modeling and sub-grid scale stabilization of incompressibility and convection in the numerical simulation of the material flow around the probe tool in a friction stir welding (FSW) process. Within the paradigmatic framework of the multiscale stabilization methods, suitable pressure and convective derivative of the temperature sub-grid scale stabilized coupled thermomechanical formulations have been developed using an Eulerian description. Norton-Hoff and Sheppard-Wright thermo-rigid-viscoplastic constitutive material models have been considered. Constitutive equations for the sub-grid scale models have been proposed and an approximation of the sub-grid scale variables has been given. In particular, algebraic sub-grid scale (ASGS) and orthogonal sub-grid scale (OSGS) methods for mixed velocity, pressure and temperature P1/P1/P1 linear elements have been considered. Furthermore, it has been shown that well known classical stabilized formulations, such as the Galerkin least-squares (GLS) for incompressible (or quasi-incompressible) problems or the Streamline Upwind/Petrov-Galerkin (SUPG) method for convection dominant problems, can be recovered as particular cases of the multiscale stabilization framework considered. Using a product formula algorithm for the solution of the coupled thermomechanical problem, the resulting algebraic system of equations has been solved using a staggered procedure in which a mechanical problem, defined by the linear momentum balance equation, under quasi-static conditions, and the incompressibility equation, is solved first at constant temperature. Then a thermal problem, defined by the energy balance equation, is solved keeping constant the mechanical variables, i.e. velocity and pressure. The computational model has been implemented in an enhanced version of the finite element software COMET, developed by the authors at the International Center for Numerical Methods in Engineering (CIMNE). Two numerical examples have been considered. The first one deals with the numerical simulation of a coupled thermomechanical flow in a 2D rectangular domain. Steady-state and transient conditions have been considered. The goal of this numerical example has been the comparison between different sub-grid scale stabilization methods for the velocity and temperature equations. In particular, using a GLS stabilization method for the pressure equation, a comparison between SUPG and OSGS convective stabilization methods has been performed. Additionally, using a SUPG stabilization method for the temperature equation, a comparison between GLS and OSGS pressure stabilization methods has been done. The second example deals with the 3D numerical simulation of a representative FSW process. Numerical results obtained have been compared with experimental results available in the literature. A good agreement on the temperature distribution has been obtained and predicted peak temperatures compare well, both in value and position, with the experimental results available.

Abstract

The use of stabilization methods is becoming an increasingly well-accepted

technique due to their success in dealing with numerous numerical pathologies

that arise in a variety of applications in computational mechanics. In this monograph a multiscale finite element method technique to deal with pressure stabilization of nearly incompressibility problems in nonlinear solid mechanics at small and finite deformations J2 plasticity is presented. A mixed formulation involving pressure and displacement fields is used as starting point. Within the finite element discretization setting, continuous linear interpolation for both fields is considered. To overcome the Babuˇska-Brezzi stability condition, a multiscale stabilization method based on the Orthogonal Subgrid Scale (OSGS) technique is introduced. Suitable nonlinear expression of the stabilization parameters are proposed. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate and, therefore, very convenient for practical industrial applications. Numerical results obtained using the OSGS stabilization technique are compared with results provided by the P1 standard Galerkin displacements linear triangular/tehrahedral element, P1/P1 standard mixed linear displacements/ linear pressure triangular/tetrahedral element and Q1/P0 mixed bilinear/ trilinear displacements/constant pressure quadrilateral/hexahedral element for 2D/3D nearly incompressible problems in the context of nonlinear small and finite deformation J2 plasticity models.

Abstract

In this paper we present a stabilized finite element formulation to solve the Oseen equations as a model problem involving both convection effects and the incompressibility restriction. The need for stabilization techniques to solve this problem arises because of the restriction in the possible choices for the velocity and pressure spaces dictated by the inf–sup condition, as well as the instabilities encountered when convection is dominant. Both can be overcome by resorting from the standard Galerkin method to a stabilized formulation. The one presented here is based on the subgrid scale concept, in which unresolvable scales of the continuous solution are approximately accounted for. In particular, the approach developed herein is based on the assumption that unresolved subscales are orthogonal to the finite element space. It is shown that this formulation is stable and optimally convergent for an adequate choice of the algorithmic parameters on which the method depends.

Abstract

Purpose

This paper aims to present a finite element formulation to approximate systems of reaction–diffusion–advection equations, focusing on cases with nonlinear reaction. The formulation is based on the orthogonal sub-grid scale approach, with some simplifications that allow one to stabilize only the convective term, which is the source of potential instabilities. The space approximation is combined with finite difference time integration and a Newton–Raphson linearization of the reactive term. Some numerical examples show the accuracy of the resulting formulation. Applications using classical nonlinear reaction models in population dynamics are also provided, showing the robustness of the approach proposed.

Design/methodology/approach

A stabilized finite element method for advection–diffusion–reaction equations to the problem on nonlinear reaction is adapted. The formulation designed has been implemented in a computer code. Numerical examples are run to show the accuracy and robustness of the formulation.

Findings

The stabilized finite element method from which the authors depart can be adapted to problems with nonlinear reaction. The resulting method is very robust and accurate. The framework developed is applicable to several problems of interest by themselves, such as the predator–prey model.

Originality/value

A stabilized finite element method to problems with nonlinear reaction has been extended. Original contributions are the design of the stabilization parameters and the linearization of the problem. The application examples, apart from demonstrating the validity of the numerical model, help to get insight in the system of nonlinear equations being solved.

Abstract

The thermally coupled incompressible inductionless magnetohydrodynamics (MHD) problem models the flow of an electrically charged fluid under the influence of an external electromagnetic field with thermal coupling. This system of partial differential equations is strongly coupled and highly nonlinear for real cases of interest. Therefore, fully implicit time integration schemes are very desirable in order to capture the different physical scales of the problem at hand. However, solving the multiphysics linear systems of equations resulting from such algorithms is a very challenging task which requires efficient and scalable preconditioners. In this work, a new family of recursive block LU preconditioners is designed and tested for solving the thermally coupled inductionless MHD equations. These preconditioners are obtained after splitting the fully coupled matrix into one-physics problems for every variable (velocity, pressure, current density, electric potential and temperature) that can be optimally solved, e.g., using preconditioned domain decomposition algorithms. The main idea is to arrange the original matrix into an (arbitrary) 2 x 2 block matrix, and consider an LU preconditioner obtained by approximating the corresponding Schur complement. For every one of the diagonal blocks in the LU preconditioner, if it involves more than one type of unknowns, we proceed the same way in a recursive fashion. This approach is stated in an abstract way, and can be straightforwardly applied to other multiphysics problems. Further, we precisely explain a flexible and general software design for the code implementation of this type of preconditioners.

Abstract

We present three new stabilized finite element (FE) based Petrov–Galerkin methods for the convection–diffusion–reaction (CDR), the Helmholtz and the Stokes problems, respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework.

It was found that the use of some standard practices (e. g.M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not appropriate when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the CDR problem. The structure
of the method in 1d is identical to the consistent approximate upwind (CAU) Petrov–Galerkin method [68] except for the definitions of the stabilization parameters.

Such a structure may also be attained via the Finite Calculus (FIC) procedure [141] by an appropriate definition of the characteristic length. The prefix high-resolution is used here in the sense popularized by Harten, i. e.second order accuracy for smooth/regular regimes and good shock-capturing in non-regular regimes. The design procedure in 1d embarks on the problem of circumventing the Gibbs phenomenon observed in L2 projections. Next, we study the conditions on the stabilization parameters to circumvent the global oscillations due to the convective term. A conjuncture of the two results is made to deal with the problem at hand that is usually plagued by Gibbs, global and dispersive oscillations in the numerical solution. A multi dimensional extension of the HRPG method using multi-linear block finite elements is also presented.
Next, we propose a higher-order compact scheme (involving two parameters) on structured meshes for the Helmholtz equation. Making the parameters equal, we recover the alpha-interpolation of the Galerkin finite element method (FEM) and the classical central finite difference method. In 1d this scheme is identical to the alphainterpolation method [140] and in 2d choosing the value 0.5 for both the parameters, we recover the generalized fourth-order compact Padé approximation [81, 168] (therein using the parameter
gamma = 2). We follow [10] for the analysis of this scheme and its performance on square meshes is compared with that of the quasi-stabilized FEM [10]. Generic expressions for the parameters are given that guarantees a dispersion accuracy of sixth-order should the parameters be distinct and fourth-order should they be equal. In the later case, an expression for the parameter is given that minimizes the maximum relative phase error in 2d. A Petrov–Galerkin formulation that yields the aforesaid scheme on structured meshes is also presented. Convergence studies of the
error in the L2 norm, the H1 semi-norm and the l1 Euclidean norm is done and the pollution effect is found to be small.

Finally, we present a collection of stabilized FE methods derived via first-order and second-order FIC procedures for the Stokes problem. It is shown that several well known existing stabilized FE methods such as the penalty technique, the Galerkin Least Square (GLS) method [93], the Pressure Gradient Projection (PGP) method [35] and the orthogonal sub-scales (OSS) method [34] are recovered from the general
iii residual-based FIC stabilized form. A new family of Pressure Laplacian Stabilization (PLS) FE methods with consistent nonlinear forms of the stabilization parameters are derived. The distinct feature of the family of PLS methods is that they are nonlinear and residual-based, i. e.the stabilization terms depend on the discrete residuals of the momentum and/or the incompressibility equations. The advantages and disadvantages of these stabilization techniques are discussed and several examples of application are presented.

Abstract

In this paper we present a stabilized finite element formulation to solve the Oseen equations as a model problem involving both convection effects and the incompressibility restriction. The need for stabilization techniques to solve this problem arises because of the restriction in the possible choices for the velocity and pressure spaces dictated by the inf-sup condition, as well as the instabilities encountered when convection is dominant. Both can be overcome by resorting from the standard Galerkin method to a stabilized formulation. The one presented here is based on the subgrid scale concept, in which unresolvable scales of the continuous solution are approximately accounted for. In particular, the approach developed herein is based on the assumption that unresolved subscales are orthogonal to the finite element space. The motivation of the method is fully described. It is also shown that this formulation is table and optimally convergent for an adequate choice of the algorithmic on which the method depends.

Abstract

In this work we present a stabilized finite element method for the stationary magneto-hydrodynamic equations based on a simple algebraic version of the subgrid scale variational concept. The linearization that yields a well posed linear problem is first identified, and for this linear problem the stabilization method is designed. The key point is the correct behavior of the stabilization parameters on which the formulation depends. It is shown that their expression can be obtained only on the basis of having a correct error estimate. For the stabilization parameters chosen, a stability estimate is proved in detail, as well as the convergence of the numerical solution to the continuous one. The method is then extended to nonlinear problems and its performance checked through numerical experiments.