Mechanical-electro-magnetic coupling enriched finite element method for static analysis of magneto-electro-elastic composites

An enriched finite element method with interpolation cover functions for acoustic analysis in high frequencies

In this work, the functionally graded magneto-electro-elastic (FG-MEE) structures were used to study mechanical performance under a thermal environment using the thermo-mechanical-electro-magnetic enriched finite element method (TMEM-EFEM). Through enhancing the traditional solution space of the finite element method (FEM), additional degrees of freedom are introduced without modifying the mesh structure. Therefore, as a class of higher-order finite element methods, field variables with high gradients can be captured without a relatively finer mesh. The accuracy of TMEM-EFEM is improved by using interpolation covers. The feasibility of TMEM-EFEM is proven by numerical examples. The proposed TMEM-EFEM is proven to be accurate under the relatively coarse mesh. Through numerical examples, the high accuracy of TMEM-EFEM is demonstrated. Therefore, TMEM-EFEM can be a good alternative for the statics study of FG-MEE structures.

ABSTRACTThis paper proposes an enriched equilibrium finite element method (EFEM) for addressing crack problems within the framework of linear elastic fracture mechanics and complementary energy principle. The key ingredient of enriched EFEM is to construct the equilibrated stress field using equilibrated tractions along with asymptotic enrichment. Specifically, an auxiliary field is employed to capture the asymptotic singular behavior of the near‐tip stress field, while the traction‐based equilibrium finite element accounts for the residual terms. The remarkable property of the proposed enriched EFEM is that the stress intensity factors (SIFs) and the stress fields are directly obtained with high accuracy, even on a very coarse mesh. Moreover, in combination with the compatible solution acquired by the enriched finite element method (FEM), strict estimation of the discretization error is obtained for crack problems by dual analysis. Numerical examples are finally conducted to verify the performance of the proposed method on strict error estimation and accurate computation of the SIFs.

Elastodynamic problems are investigated in this work by employing the enriched finite element method (EFEM) with various enrichment functions. By performing the dispersion analysis, it is confirmed that for elastodynamic analysis, the amount of numerical dispersion, which is closely related to the numerical error from the space domain discretization, can be suppressed to a very low level when quadric polynomial bases are employed to construct the local enrichment functions, while the amount of numerical dispersion from the EFEM with other types of enrichment functions (linear polynomial bases or first order of trigonometric functions) is relatively large. Consequently, the present EFEM with a quadric polynomial enrichment function shows more powerful capacities in elastodynamic analysis than the other considered numerical techniques. More importantly, the attractive monotonic convergence property can be broadly realized by the present approach with the typical two-step Bathe temporal discretization technique. Three representative numerical experiments are conducted in this work to verify the abilities of the present approach in elastodynamic analysis.

This paper presents a novel numerical method for simulating the transport and dispersion of pollutants in the Mediterranean sea. The governing mathematical equations consist of a barotropic ocean model with friction terms, bathymetric forces, Coriolis and wind stresses coupled to an advection–diffusion equation with anisotropic dispersion tensor and source terms. The proposed numerical solver uses a multilevel adaptive semi-Lagrangian finite element method that combines various techniques, including the modified method of characteristics, finite element discretization, coupled projection scheme based on a rotational pressure correction algorithm, and an adaptive L2-projection. The approach employs the gradient of the concentration as an error indicator for enrichment adaptations and increasing the number of quadrature points where needed without refining the mesh. The method is shown to provide accurate and efficient simulations for pollution transport in the Mediterranean sea. The proposed approach distinguishes itself from the well-established adaptive finite element methods for incompressible viscous flows by retaining the same structure and dimension of linear systems during the adaptation process.

Extended finite element modeling of deformable porous media with arbitrary interfaces

Enriched finite element methods have gained traction in recent years for modeling problems with material interfaces and cracks. By means of enrichment functions that incorporate a priori behavior about the solution, these methods decouple the finite element (FE) discretization from the geometric configuration of such discontinuities. Taking advantage of this greater flexibility, recent studies have proposed the adoption of Non-Uniform Rational B-Splines (NURBS) to preserve the interfaces' exact geometries throughout the analysis. In this article, we investigate NURBS-based geometries in the context of the Discontinuity-Enriched Finite Element Method (DE-FEM) based on linear field approximations. While optimal convergence is retained for problems with weak discontinuities without singularities, representing exact geometry via NURBS does not yield noticeable improvements when extracting stress intensity factors of cracked specimens. For low-order elements, we conclude that the benefits of exact geometry representation do not outweigh the increased complexity in formulation and implementation. The choice of linear FEs hinders the accuracy of the proposed formulation, suggesting that its full potential may only be unleashed by increasing the field representation order.

An enriched finite element method is presented to numerically solve the eigenvalue problem on electromagnetic waveguides governed by the Helmholtz equation. In this work, a highly efficient, simple and precise higher order subparametric method was developed using a 2D automated mesh generator performed with JuliaFEM. The transcendence computerized discretization code in Julia is developed for the present work. For curved waveguide structures, meshes with one side and curving higher orders are proposed with triangular elements with parabolic arcs. The technique is shown for distinct waveguide constructions, and the results are compared with the strongest numerical or analytical results available. The results demonstrate that the proposed methodology is effective and accurate for generating finite element simulations for complex structures with black holes and irregular topology due to no curvature loss. This article presents a finding cutoff frequency performed with JuliaFEM—an open-source program. Analysis results produced by commercial software are considered for the comparison and show that the calculation results between the two programs do not differ significantly. This procedure can be used to achieve the most effective transmission of energy for electromagnetic applications.

A weakly singular symmetric Galerkin boundary element method (SGBEM) is coupled with the standard finite element method (FEM) in order to establish an accurate and efficient numerical technique for analysis of fractures in three‐dimensional, anisotropic, linearly elastic media. In the strategy, the weakly singular SGBEM developed by Rungamornrat and Mear (textitInt. J. Solids Struct. 2008; 45:1283–1301; Comput. Methods Appl. Mech. Engrg 2008; 197:4319–4332) is utilized to model a small‐scale region containing the crack while the (possibly large and complex) compliment region is treated by the FEM. The coupled technique exploits the positive features of both methods; the SGBEM proves to be a convenient and highly accurate method for obtaining mixed‐mode stress intensity factors along the crack front, whereas the FEM is very efficient for modeling large‐scale problems in the absence of cracks. An important aspect of the formulation and implementation of the technique is that continuity of displacement and traction across the interface between the SGBEM and FEM regions is enforced in a weak sense. This allows the two regions (one modeled by the SGBEM and the other by the FEM) to be discretized independently without the need for the resulting meshes to conform on the interface separating the regions, and this flexibility in the discretization process leads to a significant reduction in the modeling effort. To demonstrate the utility and accuracy of the technique, several boundary value problems involving both embedded and surface breaking cracks are treated, and it is shown that the coupled technique yields highly accurate stress intensity factors that exhibit only a slight dependence upon mesh refinement. Copyright © 2010 John Wiley & Sons, Ltd.

The workshop Nonstandard Finite Element Methods , organized by Susanne C. Brenner (Baton Rouge), Carsten Carstensen (Berlin) and Peter Monk (Newark) was held August 10 – 16, 2008. This meeting was well attended with over 40 participants with broad geographic representation. Although Courant is often credited with the discovery of finite elements [2, 3], the first practical use of such methods in engineering seems to date back to the work of Argyris and Clough et. al. (apparently Clough coined the name “finite elements”) in the 1950s [1, 4]. By the 1970s the finite element method was firmly established in engineering practice and the basic theory of conforming elements for elliptic problems was well understood [2]. Engineers and mathematicians have since expanded the use of finite elements to a wide variety of new applications and have improved the theoretical underpinnings of the method. In addition, efficient computational algorithms have been developed, such as the multigrid method, that allow for the application of finite elements to sophisticated three dimensional problems. Within current research, we can distinguish classical or standard finite element methods, and newer nonstandard finite elements The former include conforming elements applied to “standard” Galerkin discretizations of problems. Usually accuracy is obtained by mesh refinement. As the range of applications has increased, and to alleviate perceived shortcomings in standard conforming elements, numerical analysts and engineers have sought more general finite element methods based, for example, on mixed formulations, and non-conforming or even discontinuous Galerkin schemes. These methods may be superior to standard finite elements in particular applications because of enhancements to stability, robustness or conservation properties. We refer to them as non-standard finite element methods. This Oberwolfach workshop was devoted to non-standard finite element methods and their analyses. Inevitability such a snapshot of the subject is biased, but we have tried to include non-standard methods in the broadest sense including mixed, discontinuous Galerkin, generalized, partition of unity and mortar FEM methods. All these schemes have in common that stability and convergence is not obvious and requires mathematical analysis. This is even more true for developing fast solvers, a posteriori error estimation and adaptive mesh design. Through the week there were 29 presentations on various non-standard topics. These provided an overview of current research directions, new developments and open problems.

Fluorescence diffuse optical tomography is a powerful tool for the investigation of molecular events in studies for new therapeutic developments. Here, the emphasis is put on the mathematical problem of tomography, which can be formulated in terms of an estimation of physical parameters appearing as a set of Partial Differential Equations (PDEs). The standard polynomial Finite Element Method (FEM) is a method of choice to solve the diffusion equation because it has no restriction in terms of neither the geometry nor the homogeneity of the system, but it is time consuming. In order to speed up computation time, this paper proposes an alternative numerical model, describing the diffusion operator in orthonormal basis of compactly supported wavelets. The discretization of the PDEs yields to matrices which are easily computed from derivative wavelet product integrals. Due to the shape of the wavelet basis, the studied domain is included in a regular fictitious domain. A validation study and a comparison with the standard FEM are conducted on synthetic data.

A new numerical approach for fluid-structure interaction of viscous fluid flow and flexible structures of negligible thickness (e.g. membranes, plates) on a topologically fixed fluid discretization is presented. The linear elastic and geometrically nonlinear structure is embedded [1] into the flow field described by the incompressible Navier-Stokes equations using locally enriched space-time (EST) finite elements [2], [3] to consider resulting strong and weak discontinuities in the fluid field appropriately [4]. Since the formulation of fluid, structure and coupling conditions uniformly uses velocities as unknowns and the integration of governing equations is perfomed on the deformed space-time mesh, the simulation of the strongly coupled physical domains becomes very comfortable and results in a monolithic system [5]. Numerical examples of fluid-stucture interaction processes show the ability of the developed numerical method to describe those situations of coupled systems for that common mesh moving strategies are not useful applicable and would require the introduction of time-consuming remeshing algorithms.

We present a framework for coupled multiphysics in computational fluid dynamics, targeting massively parallel systems. Our strategy is based on general problem formulations in the form of partial differential equations and the finite element method, which open for automation, and optimization of a set of fundamental algorithms. We describe these algorithms, including finite element matrix assembly, adaptive mesh refinement and mesh smoothing; and multiphysics coupling methodologies such as unified continuum fluid-structure interaction (FSI), and aeroacoustics by coupled acoustic analogies. The framework is implemented as FEniCS open source software components, optimized for massively parallel computing. Examples of applications are presented, including simulation of aeroacoustic noise generated by an airplane landing gear, simulation of the blood flow in the human heart, and simulation of the human voice organ.

It has been well known that exact closed-form solutions are not available for non-Levy-type plates. Thus, more accurate and efficient computational methods have been required for the plates subjected to arbitrary boundary conditions. This paper presents a frequency-domain spectral element model for the rectangular finite plate element. The spectral element model is developed by using two methods in combination: (1) the boundary splitting and (2) the super spectral element method in which the Kantorovich method-based finite strip element method and the frequency-domain waveguide method are utilized. The present spectral element model has nodes on four edges of the finite plate element, but no nodes inside. This can reduce the total number of degrees of freedom a lot to improve the computational efficiency significantly, when compared with the standard finite element method (FEM). The high solution accuracy and computational efficiency of the present spectral element model are evaluated by the comparison with exact solutions and the solutions by the standard FEM.

An enriched finite element method for the multi‐dimensional Stefan problems is presented. In this method the standard finite element basis is enriched with a discontinuity in the derivative of the temperature normal to the interface. The approximation can then represent the phase interface and the associated discontinuity in the temperature gradient within an element. The phase interface can be evolved without re‐meshing or the use of artificial heat capacity techniques. The interface is described by a level set function that is updated by a stabilized finite element scheme. Several examples are solved by the proposed method to demonstrate the accuracy and robustness of the method. Copyright © 2001 John Wiley & Sons, Ltd.