Unstructured finite element method for the solution of the Boussinesq problem in three dimensions

Author(s):  
C.A. Smethurst ◽  
D.J. Silvester ◽  
M.D. Mihajlović
2017 ◽  
Vol 22 (1) ◽  
pp. 133-156 ◽  
Author(s):  
Yu Du ◽  
Zhimin Zhang

AbstractWe study the error analysis of the weak Galerkin finite element method in [24, 38] (WG-FEM) for the Helmholtz problem with large wave number in two and three dimensions. Using a modified duality argument proposed by Zhu and Wu, we obtain the pre-asymptotic error estimates of the WG-FEM. In particular, the error estimates with explicit dependence on the wave numberkare derived. This shows that the pollution error in the brokenH1-norm is bounded byunder mesh conditionk7/2h2≤C0or (kh)2+k(kh)p+1≤C0, which coincides with the phase error of the finite element method obtained by existent dispersion analyses. Herehis the mesh size,pis the order of the approximation space andC0is a constant independent ofkandh. Furthermore, numerical tests are provided to verify the theoretical findings and to illustrate the great capability of the WG-FEM in reducing the pollution effect.


2021 ◽  
Author(s):  
Alejandro Ortega Laborin ◽  
Yann MALECOT ◽  
Emmanuel ROUBIN ◽  
Laurent DAUDEVILLE

This paper studies the use of the Embedded Finite Element Method (E-FEM) for the numerical modelling of triaxial fracture processes in non-homogeneous quasi-brittle materials. The E-FEM framework used in this study combines two kinematics enhancements: a weak discontinuity allowing the model to account for material heterogeneities and a strong discontinuity allowing the model to represent local fractures. The strong discontinuity features enriched fracture kinematics that allow the modelling of all typical fracture modes in three dimensions. A brief review is done of past work using similar enriched finite element frameworks to approach this problem. The work continues by establishing the theoretical basis of each kind of discontinuity formulation and their superposition through the Hu-Washizu variational principle. Afterwards, two groups of simulations have been done for discussing the performance of this combined E-FEM model: homogeneous simulations and simple heterogeneous simulations. Simple homogeneous material simulations aim to test the capabilities of the strong discontinuity model featuring full 3-D kinematics. Simple heterogeneous simulations show numerical applications of the model to the problem of a single spherical inclusion embedded into a homogeneous matrix. Comparisons will be made with another E-FEM model considering a single local fracture mode approach to discuss the differences on the representation of fracture physics under all explored conditions. A concluding statement is made on the benefits and complications identified for the E-FEM framework in this kind of applications.


Author(s):  
Karl Hollaus

Purpose The simulation of eddy currents in laminated iron cores by the finite element method (FEM) is of great interest in the design of electrical devices. Modeling each laminate by finite elements leads to extremely large nonlinear systems of equations impossible to solve with present computer resources reasonably. The purpose of this study is to show that the multiscale finite element method (MSFEM) overcomes this difficulty. Design/methodology/approach A new MSFEM approach for eddy currents of laminated nonlinear iron cores in three dimensions based on the magnetic vector potential is presented. How to construct the MSFEM approach in principal is shown. The MSFEM with the Biot–Savart field in the frequency domain, a higher-order approach, the time stepping method and with the harmonic balance method are introduced and studied. Findings Various simulations demonstrate the feasibility, efficiency and versatility of the new MSFEM. Originality/value The novel MSFEM solves true three-dimensional eddy current problems in laminated iron cores taking into account of the edge effect.


2020 ◽  
Vol 54 (5) ◽  
pp. 1525-1568 ◽  
Author(s):  
Eligio Colmenares ◽  
Gabriel N. Gatica ◽  
Sebastián Moraga

In this paper we propose and analyze, utilizing mainly tools and abstract results from Banach spaces rather than from Hilbert ones, a new fully-mixed finite element method for the stationary Boussinesq problem with temperature-dependent viscosity. More precisely, following an idea that has already been applied to the Navier–Stokes equations and to the fluid part only of our model of interest, we first incorporate the velocity gradient and the associated Bernoulli stress tensor as auxiliary unknowns. Additionally, and differently from earlier works in which either the primal or the classical dual-mixed method is employed for the heat equation, we consider here an analogue of the approach for the fluid, which consists of introducing as further variables the gradient of temperature and a vector version of the Bernoulli tensor. The resulting mixed variational formulation, which involves the aforementioned four unknowns together with the original variables given by the velocity and temperature of the fluid, is then reformulated as a fixed point equation. Next, we utilize the well-known Banach and Brouwer theorems, combined with the application of the Babuška-Brezzi theory to each independent equation, to prove, under suitable small data assumptions, the existence of a unique solution to the continuous scheme, and the existence of solution to the associated Galerkin system for a feasible choice of the corresponding finite element subspaces. Finally, we derive optimal a priori error estimates and provide several numerical results illustrating the performance of the fully-mixed scheme and confirming the theoretical rates of convergence.


2013 ◽  
Vol 444-445 ◽  
pp. 671-675
Author(s):  
Jian Ming Zhang ◽  
Yong He

In recent three decades, the finite element method (FEM) has rapidly developed as an important numerical method and used widely to solve large-scale scientific and engineering problems. In the fields of structural mechanics such as civil engineering , automobile industry and aerospace industry, the finite element method has successfully solved many engineering practical problems, and it has penetrated almost every field of today's sciences and engineering, such as material science, electricmagnetic fields, fluid dynamics, biology, etc. In this paper, we will overview and summarize the development of the p and h-p version finite element method, and introduce some recent new development and our newest research results of the p and h-p version finite element method with quasi-uniform meshes in three dimensions for elliptic problems.


Sign in / Sign up

Export Citation Format

Share Document