Mixed finite element method for coupled thermo‐hydro‐mechanical process in poro‐elasto‐plastic media at large strains

A mixed finite element for coupled thermo‐hydro‐mechanical (THM) analysis in unsaturated porous media is proposed. Displacements, strains, the net stresses for the solid phase; pressures, pressure gradients, Darcy velocities for pore water and pore air phases; temperature, temperature gradients, the total heat flux are interpolated as independent variables. The weak form of the governing equations of coupled THM problems in porous media within the element is given on the basis of the Hu–Washizu three‐filed variational principle. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of the THM constitutive model for unsaturated porous media based on the CAP model. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elasto‐plastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is utilized. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization and the softening behaviours caused by thermal and chemical effects. Copyright © 2005 John Wiley & Sons, Ltd.     

A thermo‐hydro‐mechanical model for multiphase geomaterials in dynamics with application to strain localization simulation

In this paper, the non‐isothermal elasto‐plastic behaviour of multiphase geomaterials in dynamics is investigated with a thermo‐hydro‐mechanical model of porous media. The supporting mathematical model is based on averaging procedures within the hybrid mixture theory. A computationally efficient reduced formulation of the macroscopic balance equations that neglects the relative acceleration of the fluids, and the convective terms is adopted. The modified effective stress state is limited by the Drucker–Prager yield surface. Small strains and dynamic loading conditions are assumed. The standard Galerkin procedure of the finite element method is applied to discretize the governing equations in space, while the generalized Newmark scheme is used for the time discretization. The final non‐linear set of equations is solved by the Newton method with a monolithic approach. Coupled dynamic analyses of strain localization in globally undrained samples of dense and medium dense sands are presented as examples. Vapour pressure below the saturation water pressure (cavitation) develops at localization in case of dense sands, as experimentally observed. A numerical study of the regularization properties of the finite element model is shown and discussed. A non‐isothermal case of incipient strain localization induced by temperature increase where evaporation takes place is also analysed. Copyright © 2016 John Wiley & Sons, Ltd.     

Fractional step methods for thermo‐mechanical damage analyses at transient elevated temperatures

In the interest of computational efficiency this paper describes the implementation of a coupled thermo‐damage constitutive model into a coupled time‐stepping analysis using fractional step methods. To begin it is demonstrated that a thermo‐damage model can be presented in a thermodynamic framework with the evolution equations satisfying the first and second laws of thermodynamics. The equations of evolution are partitioned in two ways, thus defining two fractional step methods: an isothermal method and an isentropic method. When implemented into a time‐stepping algorithm the isentropic method maintains a precise energy balance for the entire analysis where as the isothermal method can only provide an energy balance at the end of each thermal time step. In addition, a stability analysis shows that the isentropic analysis is unconditionally stable while a isothermal analysis is at best conditionally stable. Simulations of thermal fracture in a restrained specimen under heat show stable growth of damage to failure. Copyright © 2000 John Wiley & Sons, Ltd.     

Finite element formulation for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems

A finite element algorithm has been developed for the efficient analysis of smart composite structures with piezoelectric polymer sensors or/and actuators based on piezoelectro‐hygro‐thermo‐viscoelasticity. Variational principles for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems have also been proposed in this study. As illustrative studies, dynamic responses in laminated composite beams and plates with PVDF sensors and actuators are obtained as functions of time using the present finite element procedures. The voltage feedback control scheme is utilized. The proposed numerical method can be used for analysing problems in the design of smart structures as well as smart sensors and actuators. Copyright © 1999 John Wiley & Sons, Ltd.     

Object‐oriented finite element analysis of thermo‐hydro‐mechanical (THM) problems in porous media

The design, implementation and application of a concept for object‐oriented in finite element analysis of multi‐field problems is presented in this paper. The basic idea of this concept is that the underlying governing equations of porous media mechanics can be classified into different types of partial differential equations (PDEs). In principle, similar types of PDEs for diverse physical problems differ only in material coefficients. Local element matrices and vectors arising from the finite element discretization of the PDEs are categorized into several types, regardless of which physical problem they belong to (i.e. fluid flow, mass and heat transport or deformation processes). Element (ELE) objects are introduced to carry out the local assembly of the algebraic equations. The object‐orientation includes a strict encapsulation of geometrical (GEO), topological (MSH), process‐related (FEM) data and methods of element objects. Geometric entities of an element such as nodes, edges, faces and neighbours are abstracted into corresponding geometric element objects (ELE–GEO). The relationships among these geometric entities form the topology of element meshes (ELE–MSH). Finite element objects (ELE–FEM) are presented for the local element calculations, in which each classification type of the matrices and vectors is computed by a unique function. These element functions are able to deal with different element types (lines, triangles, quadrilaterals, tetrahedra, prisms, hexahedra) by automatically choosing the related element interpolation functions. For each process of a multi‐field problem, only a single instance of the finite element object is required. The element objects provide a flexible coding environment for multi‐field problems with different element types. Here, the C++ implementations of the objects are given and described in detail. The efficiency of the new element objects is demonstrated by several test cases dealing with thermo‐hydro‐mechanical (THM) coupled problems for geotechnical applications. Copyright © 2006 John Wiley & Sons, Ltd.     

A monolithic computational approach to thermo‐structure interaction

In the present work, a monolithic solution approach for thermo‐structure interaction problems motivated by the challenging application of the behaviour of rocket nozzles is proposed. Structural and thermal fields are independently discretised via finite elements. The resulting system of equations is solved via a monolithic thermo‐structure interaction scheme, which is constructed by a block Gauss–Seidel preconditioner in combination with algebraic multigrid methods. The proposed method is tested for four numerical examples, the second Danilovskaya problem, a simplified rocket nozzle configuration, an internally loaded hollow sphere, and a fully three‐dimensional nozzle configuration of a subscale thrust chamber. Good agreement of the numerical results with results from the literature is observed. Furthermore, it is shown that the monolithic solution algorithm can handle the complete range of the parameter spectrum, whereas partitioned algorithms are limited to a certain parameter range only. Moreover, the monolithic algorithm exhibits improved efficiency and robustness compared to partitioned algorithms. Copyright © 2013 John Wiley & Sons, Ltd.     

Thermo‐hydro‐mechanical modeling of fluid geological storage by Godunov‐mixed methods

The efficient solution to the coupled system of PDEs governing the mass and the energy balance in deformable porous media requires advanced numerical algorithms. A combination of mixed/Galerkin finite elements and finite volumes along with a staggered method are employed. Fluid flow and heat transfer are addressed iteratively by a fully coupled approach and the medium deformation by an explicitly coupled scheme, at each time step. Such formulation allows for stable numerical solutions, element‐wise conservative velocity fields and accurate prediction of sharp temperature convective fronts. The proposed model is experimented with in realistic applications of a deep aquifer fluid injection. Copyright © 2012 John Wiley & Sons, Ltd.     

A stabilized finite element formulation for the transonic small‐disturbance system

A stabilized finite element formulation for the transonic small‐disturbance system of equations is developed and used to solve a variety of problems in transonic aerodynamics. An adaptive mesh refinement technique and a common discontinuity capturing operator are used to resolve regions with large gradients in the velocity field. The scheme works well in both flow regimes, subsonic and supersonic, and captures shocks naturally. Agreement with available experimental observations and theoretical approximations is very good. Copyright © 2001 John Wiley & Sons, Ltd.     

A Galerkin/least‐squares finite element formulation for consolidation

A Galerkin/least‐squares (GLS) finite element formulation for problem of consolidation of fully saturated two‐phase media is presented. The elimination of spurious pressure oscillations appearing at the early stage of consolidation for standard Galerkin finite elements with equal interpolation order for both displacements and pressures is the goal of the approach. It will be shown that the least‐squares term, based exclusively on the residuum of the fluid flow continuity equation, added to the standard Galerkin formulation enhances its stability and can fully eliminate pressure oscillations. A reasonably simple framework designed for derivation of one‐dimensional as well as multi‐dimensional estimates of the stabilization factor is proposed and then verified. The formulation is validated on one‐dimensional and then on two‐dimensional, linear and non‐linear test problems. The effect of the fluid incompressibility as well as compressibility will be taken into account and investigated. Copyright © 2001 John Wiley & Sons Ltd.     

Eulerian formulation using stabilized finite element method for large deformation solid dynamics

This paper describes an Eulerian formulation for large deformation solid dynamics. In the present Eulerian approach, an advective equation is solved using the Stream‐Upwind/Petrov–Galerkin finite element method. The Eulerian finite element method is applied to path‐dependent solid analyses such as impact bar and ductile necking problems. These computational results using the Eulerian finite element method are compared with the results obtained from using the Lagrangian finite element method and an Eulerian formulation based on a finite difference method. Copyright © 2007 John Wiley & Sons, Ltd.     

Generalized finite element approaches for analysis of localized thermo‐structural effects

This work addresses computational modeling challenges associated with structures subjected to sharp, local heating, where complex temperature gradients in the materials cause three‐dimensional, localized, intense stress and strain variation. Because of the nature of the applied loadings, multiphysics analysis is necessary to accurately predict thermal and mechanical responses. Moreover, bridging spatial scales between localized heating and global responses of the structure is nontrivial. A large global structural model may be necessary to represent detailed geometry alone, and to capture local effects, the traditional approach of pre‐designing a mesh requires careful manual effort. These issues often lead to cumbersome and expensive global models for this class of problems. To address them, the authors introduce a generalized FEM (GFEM) approach for analyzing three‐dimensional solid, coupled physics problems exhibiting localized heating and corresponding thermomechanical effects. The capabilities of traditional hp‐adaptive FEM or GFEM as well as the GFEM with global–local enrichment functions are extended to one‐way coupled thermo‐structural problems, providing meshing flexibility at local and global scales while remaining competitive with traditional approaches. The methods are demonstrated on several example problems with localized thermal and mechanical solution features, and accuracy and (parallel) computational efficiency relative to traditional direct modeling approaches are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

A monolithic,mortar‐based interface coupling and solution scheme for finite element simulations of lithium‐ion cells

This work introduces a novel, mortar‐based coupling scheme for electrode‐electrolyte interfaces in 3‐dimensional finite element models for lithium‐ion cells and similar electrochemical systems. The coupling scheme incorporates the widely applied Butler‐Volmer charge transfer kinetics, but conceptually also works for other interface equations. Unlike conventional approaches, the coupling scheme allows flexible mesh generation for the electrode and electrolyte phases with nonmatching meshes at electrode‐electrolyte interfaces. As a result, the desired spatial mesh resolution in each phase and the resulting computational effort can be easily controlled, leading to improved efficiency. All governing equations are solved in a monolithic fashion as a holistic, unified system of linear equations for computational robustness and performance reasons. Consistency and optimal convergence behavior of the coupling scheme are demonstrated in elementary numerical tests, and the discharge of two different realistic lithium‐ion cells, each consisting of an anode, a cathode, and an electrolyte, is also simulated. One of the two cells involves about 1.35 million degrees of freedom and very complex microstructural geometries obtained from X‐ray tomography data. For validation purposes, characteristic numerical results from the literature are reproduced, and the coupling scheme is shown to require considerably fewer degrees of freedom than a standard discretization with matching interface meshes to achieve a similar level of accuracy.     

Efficient 3D‐finite element formulation for thin mechanical and piezoelectric structures

It is the aim to compute efficiently the deformations of the mechanical structures that are excited or damped by piezoelectric actuators. They are treated as 3D structures to have much flexibility. Special design of the finite element concept is required since the structures are thin walled and locking effects have to be avoided. Although the computations are performed in the framework of 3D elasticity, we use ideas from modern plate elements and mixed methods. Copyright © 2007 John Wiley & Sons, Ltd.     

A new path‐following constraint for strain‐softening finite element simulations

The application of strain‐softening constitutive relations to model the failure modes of real‐life structures is faced to numerical difficulties related to instabilities that appear as sharp snap‐backs of the structural response. A path‐following method has to complement the solution algorithm to achieve convergence despite these critical points. Because of the sharpness of the snap‐backs, it is believed essential that the path‐following constraint distinguish between a purely elastic unloading and a dissipative path. For that purpose, a new constraint based on the maximal value of the elastic predictor for the yield function is proposed. As it is highly non linear, a specific solution algorithm is required. The robustness of this constraint is illustrated by three applications: the study of crack propagations by means of a cohesive zone model, the failure of a structure submitted to nonlocal damage and the simulation of a nonlocal strain‐softening plastic specimen. Copyright © 2004 John Wiley & Sons, Ltd.     

Mixed finite element method for saturated poroelastoplastic media at large strains

A mixed finite element for hydro‐dynamic analysis in saturated porous media in the frame of the Biot theory is proposed. Displacements, effective stresses, strains for the solid phase and pressure, pressure gradients, and Darcy velocities for the fluid phase are interpolated as independent variables. The weak form of the governing equations of coupled hydro‐dynamic problems in saturated porous media within the element are given on the basis of the Hu–Washizu three‐field variational principle. In light of the stabilized one point quadrature super‐convergent element developed in solid continuum, the interpolation approximation modes for the primary unknowns and their spatial derivatives of the solid and the fluid phases within the element are assumed independently. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of pressure‐dependent non‐associated plasticity. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elastoplastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is used. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization due to strain softening in poroelastoplastic media subjected to dynamic loading at large strain. Copyright © 2003 John Wiley & Sons, Ltd.     

Assumed strain finite element methods for conserving temporal integrations in non‐linear solid dynamics

This paper presents a new class of assumed strain finite elements to use in combination with general energy‐momentum‐conserving time‐stepping algorithms so that these conservation properties in time are preserved by the fully discretized system in space and time. The case of interest corresponds to nearly incompressible material responses, in the fully non‐linear finite strain elastic and elastoplastic ranges. The new elements consider the classical scaling of the deformation gradient with an assumed Jacobian (its determinant) defined locally through a weighted averaging procedure at the element level. The key aspect of the newly proposed formulation is the definition of the associated linearized strain operator or B‐bar operator. The developments presented here start by identifying the conditions that this discrete operator must satisfy for the fully discrete system in time and space to inherit exactly the conservation laws of linear and angular momenta, and the conservation/dissipation law of energy for elastic and inelastic problems, respectively. Care is also taken of the preservation of the relative equilibria and the corresponding group motions associated with the momentum conservation laws, and characterized by purely rotational and translational motions superimposed to the equilibrium deformed configuration. With these developments at hand, a new general B‐bar operator is introduced that satisfies these conditions. The new operator not only accounts for the spatial interpolations (e.g. bilinear displacements with piece‐wise constant volume) but also depends on the discrete structure of the equations in time. The aforementioned conservation/dissipation properties of energy and momenta are then proven to hold rigorously for the final numerical schemes, unconditionally of the time step size and the material model (elastic or elastoplastic). Different finite elements are considered in this framework, including quadrilateral and triangular elements for plane problems and brick elements for three‐dimensional problems. Several representative numerical simulations are presented involving, in particular, the use of energy‐dissipating momentum‐conserving time‐stepping schemes recently developed by the author and co‐workers for general finite strain elastoplasticity in order to illustrate the properties of the new finite elements, including these conservation/dissipation properties in time and their locking‐free response in the quasi‐incompressible case. Copyright © 2007 John Wiley & Sons, Ltd.     

An evaluation of limited‐memory sparse linear solvers for thermo‐mechanical applications

We consider the performance of sparse linear solvers for problems that arise from thermo‐mechanical applications. Such problems have been solved using sparse direct schemes that enable robust solution at the expense of memory requirements that grow non‐linearly with the dimension of the coefficient matrix. In this paper, we consider a class of preconditioned iterative solvers as a limited‐memory alternative to direct solution schemes. However, such preconditioned iterative solvers typically exhibit complex trade‐offs between reliability and performance. We therefore characterize such trade‐offs for systems from thermo‐mechanical problems by considering several preconditioning schemes including multilevel methods and those based on sparse approximate inversion and incomplete matrix factorization. We provide an analysis of computational costs and memory requirements for model thermo‐mechanical problems, indicating that certain incomplete factorization schemes can achieve good performance. We also provide empirical evaluations that corroborate our analysis and indicate the relative effectiveness of different solution schemes. Our results indicate that our drop‐threshold incomplete Cholesky preconditioning is more robust, efficient and flexible than other popular preconditioning schemes. In addition, we propose preconditioner reuse to amortize preconditioner construction cost over a sequence of linear systems that arise from non‐linear solutions in a plastic regime. Copyright © 2007 John Wiley & Sons, Ltd.     

A new stabilized finite element method for reaction–diffusion problems: The source‐stabilized Petrov–Galerkin method

This paper proposes a new stabilized finite element method to solve singular diffusion problems described by the modified Helmholtz operator. The Galerkin method is known to produce spurious oscillations for low diffusion and various alternatives were proposed to improve the accuracy of the solution. The mostly used methods are the well‐known Galerkin least squares and Galerkin gradient least squares (GGLS). The GGLS method yields the exact nodal solution in the one‐dimensional case and for a uniform mesh. However, the behavior of the method deteriorates slightly in the multi‐dimensional case and for non‐uniform meshes. In this work we propose a new stabilized finite element method that leads to improved accuracy for multi‐dimensional problems. For the one‐dimensional case, the new method leads to the same results as the GGLS method and hence provides exact nodal solutions to the problem on uniform meshes. The proposed method is a Galerkin discretization used to solve a modified equation that includes a term depending on the gradient of the original partial differential equation. Copyright © 2008 John Wiley & Sons, Ltd.