Two-grid methods for semi-linear elliptic interface problems by immersed finite element methods

In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.

P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.

Solving Interface Problems of the Helmholtz Equation by Immersed Finite Element Methods

This article reports our explorations for solving interface problems of the Helmholtz equation by immersed finite elements (IFE) on interface independent meshes. Two IFE methods are investigated: the partially penalized IFE (PPIFE) and discontinuous Galerkin IFE (DGIFE) methods. Optimal convergence rates are observed for these IFE methods once the mesh size is smaller than the optimal mesh size which is mainly dictated by the wave number. Numerical experiments also suggest that higher degree IFE methods are advantageous because of their larger optimal mesh size and higher convergence rates.

Study of adaptive finite element techniques for contact problem in elastic bodies

The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.

Numerical Exploration of Asymmetrical Impact Dynamics: Unveiling Nonlinearities in Collision Problems and Resilience of Reinforced Concrete Structures

This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.

Mathematical Description and Finite Element Equation of 3D Coupled Thermo-elastic Contact Problem

Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.

A Stable FE-FD Method for Anisotropic Parabolic PDEs with Moving Interfaces

In this paper,a new finite element and finite difference(FE-FD)method has been developed for anisotropic parabolic interface problems with a known moving interface using Cartesian meshes.In the spatial discretization,the standard P,FE discretization is applied so that the part of the coefficient matrix is symmetric positive definite,while near the interface,the maximum principle preserving immersed interface discretization is applied.In the time discretization,a modified Crank-Nicolson discretization is employed so that the hybrid FE-FD is stable and second order accurate.Correction terms are needed when the interface crosses grid lines.The moving interface is represented by the zero level set of a Lipschitz continuous function.Numerical experiments presented in this paper confirm second orderconvergence.

Arc-length technique for nonlinear finite element analysis

Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.

RANDOM MICROSTRUCTURE FINITE ELEMENT METHOD FOR EFFECTIVE NONLINEAR PROPERTIES OF COMPOSITE MATERIALS

Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.

A Simple FEM for Solving Two-Dimensional Diffusion Equation with Nonlinear Interface Jump Conditions

In this paper,we propose a numerical method for solving parabolic interface problems with nonhomogeneous flux jump condition and nonlinear jump condition.The main idea is to use traditional finite element method on semi-Cartesian mesh coupled with Newton’s method to handle nonlinearity.It is easy to implement even though variable coefficients are used in the jump condition instead of constant in previous work for elliptic interface problem.Numerical experiments show that our method is about second order accurate in the L1 norm.

Formulations of a hydromechanical interface element

A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.

Geotechnical particle finite element method for modeling of soilstructure interaction under large deformation conditions

The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems that are beyond the capabilities of classical finite element method(FEM).In PFEM,the computational domain is reconfigured for optimal solution by frequent remeshing and boundary updating.PFEM inherits many concepts,such as a Lagrangian description of continuum,from classic geomechanical FEM.This familiarity with more popular numerical methods facilitates learning and application.This work focuses on G-PFEM,a code specifically developed for the use of PFEM in geotechnical problems.The article has two purposes.The first is to give the reader an overview of the capabilities and main features of the current version of the G-PFEM and the second is to illustrate some of the newer developments of the code.G-PFEM can solve coupled hydro-mechanical static and dynamic problems involving the interaction of solid and/or deformable bodies.Realistic constitutive models for geomaterials are available,including features,such as structure and destructuration,which result in brittle response.The solutions are robust,solidly underpinned by numerical technology including mixedfield formulations,robust and mesh-independent integration of elastoplastic constitutive models and a rigorous and flexible treatment of contact interactions.The novel features presented in this work include the contact domain technique,a natural way to capture contact interactions and impose contact constraints between different continuum bodies,as well as a new simplified formulation for dynamic impact problems.The code performance is showcased by the simulation of several soil-structure interaction problems selected to highlight the novel code features:a rigid footing insertion in soft rock,pipeline insertion and subsequent lateral displacement on over-consolidated clay,screw-pile pull-out and the dynamic impact of a free-falling spherical penetrometer into clay.

Finite Element Analysis of the Ramberg-Osgood Bar

In this work, we present a priori error estimates of finite element approximations of the solution for the equilibrium equation of an axially loaded Ramberg-Osgood bar. The existence and uniqueness of the solution to the associated nonlinear two point boundary value problem is established and used as a foundation for the finite element analysis.

EFFECTIVE NUMERICAL METHODS FOR ELASTO-PLASTIC CONTACT PROBLEMS WITH FRICTION

Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.

OPTIMAL QUADRATIC NITSCHE EXTENDED FINITE ELEMENT METHOD FOR INTERFACE PROBLEM OF DIFFUSION EQUATION

In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.

A Well-Conditioned, Nonconforming Nitsche’s Extended Finite Element Method for Elliptic Interface Problems

In this paper,we introduce a nonconforming Nitsche’s extended finite element method(NXFEM)for elliptic interface problems on unfitted triangulation elements.The solution on each side of the interface is separately expanded in the standard nonconforming piecewise linear polynomials with the edge averages as degrees of freedom.The jump conditions on the interface and the discontinuities on the cut edges(the segment of edges cut by the interface)are weakly enforced by the Nitsche’s approach.In the method,the harmonic weighted fluxes are used and the extra stabilization terms on the interface edges and cut edges are added to guarantee the stability and the well conditioning.We prove that the convergence order of the errors in energy and L 2 norms are optimal.Moreover,the errors are independent of the position of the interface relative to the mesh and the ratio of the discontinuous coefficients.Furthermore,we prove that the condition number of the system matrix is independent of the interface position.Numerical examples are given to confirm the theoretical results.

THE ADAPTIVE IMMERSED INTERFACE FINITE ELEMENT METHOD FOR ELASTICITY INTERFACE PROBLEMS

In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a residual-based a posteriori error estimate which is λ- independent multiplicative constants; the Lame constant λ steers the incompressibility. The error estimators are then implemented and tested with promising numerical results which will show the competitive behavior of the adaptive algorithm.

SOLVING 2-D NONLINEAR COUPLED HEAT AND MOISTURE TRANSFER PROBLEMS VIA A SELF-ADAPTIVE PRECISE ALGORITHM IN TIME DOMAIN

A self-adaptive precise algorithm in the time domain was employed to solve 2-D nonlinear coupled heat and moisture transfer problems. By expanding variables at a discretized time interval, the variations of variables can be described more precisely,and a nonlinear coupled initial and boundary value problem was converted into a series of recurrent linear boundary value problems which are solved by FE technique. In the computation, no additional assumption and the nonlinear iteration are required, and a criterion for self-adaptive computation is proposed to maintain sufficient computing accuracy for the change sizes of time steps. In the numerical comparison, the variations of material properties with temperature, moisture content, and both temperature and moisture content are taken into account, respectively. Satisfactory results have been obtained, indicating that the proposed approach is capable of dealing with complex nonlinear problems.