No Arabic abstract
This paper presents a steady-state and transient heat conduction analysis framework using the polygonal scaled boundary finite element method (PSBFEM) with polygon/quadtree meshes. The PSBFEM is implemented with commercial finite element code Abaqus by the User Element Sub-routine (UEL) feature. The detailed implementation of the framework, defining the UEL element, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from heat conduction are solved to validate the proposed implementation. Results show that the PSBFEM is reliable and accurate for solving heat conduction problems. Not only can the proposed implementation help engineering practitioners analyze the heat conduction problem using polygonal mesh in Abaqus, but it also provides a reference for developing the UEL to solve other problems using the scaled boundary finite element method.
The scaled boundary finite element method (SBFEM) is a semi-analytical computational scheme, which is based on the characteristics of the finite element method (FEM) and combines the advantages of the boundary element method (BEM). This paper integrates the scaled boundary finite element method (SBFEM) and the polygonal mesh technique into a new approach to solving the steady-state and transient seepage problems. The proposed method is implemented in Abaqus using a user-defined element (UEL). The detailed implementations of the procedure, defining the UEL element, updating the RHS and AMATRX, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from seepage are solved to validate the proposed implementation. Results show that the polygonal element of PS-SBFEM has a higher accuracy rate than the standard FEM element in the same element size. For the transient problems, the results between PS-SBFEM and the FEM are in excellent agreement. Furthermore, the PS-SBFEM with quadtree meshes shows a good effect for solving complex geometric in the seepage problem. Hence, the proposed method is robust accurate for solving the steady-state and transient seepage problems. The developed UEL source code and the associated input files can be downloaded from https://github.com/yangyLab/PS-SBFEM.
In this paper, we examine the effectiveness of classic multiscale finite element method (MsFEM) (Hou and Wu, 1997; Hou et al., 1999) for mixed Dirichlet-Neumann, Robin and hemivariational inequality boundary problems. Constructing so-called boundary correctors is a common technique in existing methods to prove the convergence rate of MsFEM, while we think not reflects the essence of those problems. Instead, we focus on the first-order expansion structure. Through recently developed estimations in homogenization theory, our convergence rate is provided with milder assumptions and in neat forms.
The locally modified finite element method, which is introduced in [Frei, Richter: SINUM 52(2014), p. 2315-2334] is a simple fitted finite element method that is able to resolve weak discontinuities in interface problems. The method is based on a fixed structured coarse mesh, which is then refined into sub-elements to resolve an interior interface. In this work, we extend the locally modified finite element method to second order using an isoparametric approach in the interface elements. Thereby we need to take care that the resulting curved edges do not lead to degenerate sub-elements. We prove optimal a priori error estimates in the $L^2$-norm and in a modified energy norm, as well as a reduced convergence order of ${cal O}(h^{3/2})$ in the standard $H^1$-norm. Finally, we present numerical examples to substantiate the theoretical findings.
In this paper we propose a finite element method for solving elliptic equations with the observational Dirichlet boundary data which may subject to random noises. The method is based on the weak formulation of Lagrangian multiplier. We show the convergence of the random finite element error in expectation and, when the noise is sub-Gaussian, in the Orlicz 2- norm which implies the probability that the finite element error estimates are violated decays exponentially. Numerical examples are included.
In this paper we proposed two new quasi-boundary value methods for regularizing the ill-posed backward heat conduction problems. With a standard finite difference discretization in space and time, the obtained all-at-once nonsymmetric sparse linear systems have the desired block $omega$-circulant structure, which can be utilized to design an efficient parallel-in-time (PinT) direct solver that built upon an explicit FFT-based diagonalization of the time discretization matrix. Convergence analysis is presented to justify the optimal choice of the regularization parameter. Numerical examples are reported to validate our analysis and illustrate the superior computational efficiency of our proposed PinT methods.