Search results

This paper aims to propose a new boundary condition and a web-spline basis of finite element space approximation to remedy the problems of constraints due to homogeneous and non-homogeneous; Dirichlet boundary conditions. This paper considered the two-dimensional linear elasticity equation of Navier–Lamé with the condition CAB. The latter allows to have a total insertion of the essential boundary condition in the linear system obtained; without using a numerical method as Lagrange multiplier. This study have developed mixed finite element; method using the B-splines Web-spline space. These provide an exact implementation of the homogeneous; Dirichlet boundary conditions, which removes the constraints caused by the standard; conditions. This paper showed the existence and the uniqueness of the weak solution, as well as the convergence of the numerical solution for the quadratic case are proved. The weighted extended B-spline; approach have become a much more workmanlike solution.

In this paper, this study used the implementation of weighted finite element methods to solve the Navier–Lamé system with a new boundary condition CA, B (Koubaiti et al., 2020), that generalises the well-known basis, especially the Dirichlet and the Neumann conditions. The novel proposed boundary condition permits to use a single Matlab code, which summarises all kind of boundary conditions encountered in the system. By using this model is possible to save time and programming recourses while reap several programs in a single directory.

The results have shown that the Web-spline-based quadratic-linear finite elements satisfy the inf–sup condition, which is necessary for existence and uniqueness of the solution. It was demonstrated by the existence of the discrete solution. A full convergence was established using the numerical solution for the quadratic case. Due to limited regularity of the Navier–Lamé problem, it will not change by increasing the degree of the Web-spline. The computed relative errors and their rates indicate that they are of order 1/H. Thus, it was provided their theoretical validity for the numerical solution stability. The advantage of this problem that uses the CA, B boundary condition is associated to reduce Matlab programming complexity.

The mixed finite element method is a robust technique to solve difficult challenges from engineering and physical sciences using the partial differential equations. Some of the important applications include structural mechanics, fluid flow, thermodynamics and electromagnetic fields (Zienkiewicz and Taylor, 2000) that are mainly based on the approximation of Lagrange. However, this type of approximation has experienced a great restriction in the level of domain modelling, especially in the case of complicated boundaries such as that in the form of curvilinear graphs. Recently, the research community tried to develop a new way of approximation based on the so-called B-spline that seems to have superior results in solving the engineering problems.

This paper aims to proposes and analyzes a novel recovery-based posteriori error estimator for the stationary natural-convection problem based on penalized finite element method.

The optimal error estimates of the penalty FEM are established by using the lower-order finite element pair P₁-P₀-P₁ which does not satisfy the discrete inf-sup condition. Besides, a new recovery type posteriori estimator in view of the gradient recovery and superconvergent theory to deal with the discontinuity of the gradient of numerical solution.

The stability, accuracy and efficiency of the proposed method are confirmed by several numerical investigations.

The provided reliability and efficiency analysis is shown that the true error can be effectively bounded by the recovery-based error estimator.

The use of the prominent finite difference time‐domain (FDTD) method for the time‐domain solution of electromagnetic wave propagation past devices with small geometrical details can require very fine grids and can lead to unmanageable computational time and storage. The purpose of this paper is to extend the analysis of a discontinuous Galerkin time‐domain (DGTD) method (able to handle possibly non‐conforming locally refined grids, based on portions of Cartesian grids) and investigate the use of perfectly matched layer regions and the coupling with a fictitious domain approach. The use of a DGTD method with a locally refined, non‐conforming mesh can help focusing on these small details. In this paper, the adaptation to the DGTD method of the fictitious domain approach initially developed for the FDTD is considered, in order to avoid the use of a volume mesh fitting the geometry near the details.

Based on a DGTD method, a fictitious domain approach is developed to deal with complex and small geometrical details.

The fictitious domain approach is a very interesting complement to the FDTD method, since it makes it possible to handle complex geometries. However, the fictitious domain approach requires small volume elements, thus making the use of the FDTD on wide, regular, fine grids often unmanageable. The DGTD method has the ability to handle easily locally refined grids and the paper shows it can be coupled to a fictitious domain approach.

Although the stability and dispersion analysis of the DGTD method is complete, the theoretical analysis of the fictitious domain approach in the DGTD context is not. It is a subject of further investigation (which could provide important insights for potential improvements).

This is believed to be the first time a DGTD method is coupled with a fictitious domain approach.

This paper aims to propose a characteristic stabilized finite element method for non-stationary conduction-convection problems.

To avoid difficulty caused by the trilinear term, the authors use the characteristic method to deal with the time derivative term and the advection term. The space discretization adopts the low-order triples (i.e. P₁-P₁-P₁ and P₁-P₀-P₁ triples). As low-order triples do not satisfy inf-sup condition, the authors use the stability technique to overcome this flaw.

The stability and the convergence analysis shows that the method is stable and has optimal-order error estimates.

Numerical experiments confirm the theoretical analysis and illustrate that the authors’ method is highly effective and reliable, and consumes less CPU time.

The purpose of this paper is to propose three iterative finite element methods for equations of thermally coupled incompressible magneto-hydrodynamics (MHD) on 2D/3D bounded domain. The detailed theoretical analysis and some numerical results are presented. The main results show that the Stokes iterative method has the strictest restrictions on the physical parameters, and the Newton’s iterative method has the higher accuracy and the Oseen iterative method is stable unconditionally.

Three iterative finite element methods have been designed for the thermally coupled incompressible MHD flow on 2D/3D bounded domain. The Oseen iterative scheme includes solving a linearized steady MHD and Oseen equations; unconditional stability and optimal error estimates of numerical approximations at each iterative step are established under the uniqueness condition. Stability and convergence of numerical solutions in Newton and Stokes’ iterative schemes are also analyzed under some strong uniqueness conditions.

This work was supported by the NSF of China (No. 11971152).

This paper presents the best choice for solving the steady thermally coupled MHD equations with different physical parameters.

The purpose of this paper is to design a parallel finite element (FE) algorithm based on fully overlapping domain decomposition for solving the nonstationary incompressible magnetohydrodynamics (MHD).

The fully discrete Euler implicit/explicit FE subproblems, which are defined in the whole domain with vast majority of the degrees of freedom associated with the particular subdomain, are solved in parallel. In each subproblem, the linear term is treated by implicit scheme and the nonlinear term is solved by explicit one.

For the algorithm, the almost unconditional convergence with optimal orders is validated by numerical tests. Some interesting phenomena are presented.

The proposed algorithm is effective, easy to realize with low communication costs and preferred for solving the strong nonlinear MHD system.

The purpose of this paper is to propose a new stabilized finite volume element method for the Navier-Stokes problem.

This new method is based on the multiscale enrichment and uses the lowest equal order finite element pairs P₁/P₁.

The stability and convergence of the optimal order in H¹-norm for velocity and L²-norm for pressure are obtained.

Using a dual problem for the Navier-Stokes problem, the convergence of the optimal order in L²-norm for the velocity is obtained. Finally, numerical example confirms the theory analysis and validates the effectiveness of this new method.

This paper aims to propose a new highly efficient iterative method based on classical Oseen iteration for the natural convection equations.

First, the authors solve the problem by the Oseen iterative scheme based on finite element method, then use the error correction strategy to control the error arising.

The new iterative method not only retains the advantage of the Oseen scheme but also saves computational time and iterative step for solving the considered problem.

In this work, the authors introduce a new iterative method to solve the natural convection equations. The new algorithm consists of the Oseen scheme and the error correction which can control the errors from the iterative step arising for solving the nonlinear problem. Comparing with the classical iterative method, the new scheme requires less iterations and is also capable of solving the natural convection problem at higher Rayleigh number.

Engineers have developed robust and efficient incompressible finite element formulations using tools such as the Patch Test and the counting of constraints/variables, the first one aimed at the development of consistent elements and the second one aimed at the development of non‐locking and stable elements. The mentioned tools are rooted in the physics of the continuum mechanics problem. Mathematicians, on the other side, developed complex and powerful tools to examine the convergence of finite element formulations, such as the inf‐sup condition, these methods are based on the properties of the elliptical PDEs that constitute the mathematical model of the continuum mechanics problem. In this paper we intend to understand the inf‐sup condition from an engineering perspective, so as to be able to incorporate it into the package of tools used in the development of finite element formulations.

This study aims to introduce the hybrid finite element (FE) – meshfree method and multiscale variational principle into the traditional mixed FE formulation, leading to a stable mixed formulation for incompressible linear elasticity which circumvents the need to satisfy inf-sup condition.

Using the hybrid FE–meshfree method, the displacement and pressure are interpolated conveniently with the same order so that a continuous pressure field can be obtained with low-order elements. The multiscale variational principle is then introduced into the Galerkin form to obtain stable and convergent results.

The present method is capable of overcoming volume locking and does not exhibit unphysical oscillations near the incompressible limit. Moreover, there are no extra unknowns introduced in the present method because the fine-scale unknowns are eliminated using the static condensation technique, and there is no need to evaluate any user-defined stability parameter as the classical stabilization methods do. The shape functions constructed in the present model possess continuous derivatives at nodes, which gives a continuous and more precise stress field with no need of an additional smooth process. The shape functions in the present model also possess the Kronecker delta property, so that it is convenient to impose essential boundary conditions.

The proposed model can be implemented easily. Its convergence rates and accuracy in displacement, energy and pressure are even comparable to those of second-order mixed elements.