Finite-element EM modelling on hexahedral grids with an FD solver as a pre-conditioner

2020 ◽  
Vol 223 (2) ◽  
pp. 840-850
Author(s):  
Nikolay Yavich ◽  
Michael S Zhdanov

SUMMARY The finite-element (FE) method is one of the most powerful numerical techniques for modelling 3-D electromagnetic fields. At the same time, there still exists the problem of efficient and economical solution of the respective system of FE equations in the frequency domain. In this paper, we concentrate on modelling with adapted hexahedral or logically rectangular grids. These grids are easy to generate, yet they are flexible enough to incorporate real topography and seismic horizons. The goal of this work is to show how a finite-difference (FD) solver can be used as a pre-conditioner for hexahedral FE modelling. Applying the lowest order Nédélec elements, we present a novel pre-conditioned iterative solver for the arising system of linear equations that combines an FD solver and simple smoothing procedure. The particular FD solver that we use relies on the implicit factorization of the horizontally layered earth matrix. We assessed runtime and accuracy of the presented approach on synthetic and real resistivity models (topography of the Black Sea continental slope). We further compared performance of our program versus publicly available Mare2DEM, ModEM and MUMPS programs/libraries. Our examples involve plane-wave and controlled source modelling. The numerical examples demonstrate that the presented approach is fast and robust for models with moderate contrast, supports highly deformed cells, and is quite memory-economical.

Author(s):  
S. Chakraverty ◽  
Diptiranjan Behera

This chapter presents the static and dynamic analysis of structures with uncertain parameters using fuzzy finite element method. Uncertainties presents in the parameters are modelled through convex normalised fuzzy sets. Fuzzy finite element method converts the structures into fuzzy system of linear equations and fuzzy eigenvalue problem for static and dynamic problems respectively. As such method to solve fuzzy system of linear equations, fully fuzzy system of linear equations and fuzzy eigenvalue problems are presented. These methods are applied to various structural problems to find out the fuzzy static and dynamic responses of the structures. Also the chapter analyses the numerical solution of uncertain fractionally damped spring-mass system. Uncertainties considered in the initial condition of the system.


Insects ◽  
2020 ◽  
Vol 11 (8) ◽  
pp. 546
Author(s):  
Shahab Eshghi ◽  
Vahid Nooraeefar ◽  
Abolfazl Darvizeh ◽  
Stanislav N. Gorb ◽  
Hamed Rajabi

The finite element (FE) method is one of the most widely used numerical techniques for the simulation of the mechanical behavior of engineering and biological objects. Although very efficient, the use of the FE method relies on the development of accurate models of the objects under consideration. The development of detailed FE models of often complex-shaped objects, however, can be a time-consuming and error-prone procedure in practice. Hence, many researchers aim to reach a compromise between the simplicity and accuracy of their developed models. In this study, we adapted Distmesh2D, a popular meshing tool, to develop a powerful application for the modeling of geometrically complex objects, such as insect wings. The use of the burning algorithm (BA) in digital image processing (DIP) enabled our method to automatically detect an arbitrary domain and its subdomains in a given image. This algorithm, in combination with the mesh generator Distmesh2D, was used to develop detailed FE models of both planar and out-of-plane (i.e., three-dimensionally corrugated) domains containing discontinuities and consisting of numerous subdomains. To easily implement the method, we developed an application using the Matlab App Designer. This application, called WingMesh, was particularly designed and applied for rapid numerical modeling of complicated insect wings but is also applicable for modeling purposes in the earth, engineering, mathematical, and physical sciences.


Fuzzy Systems ◽  
2017 ◽  
pp. 1-30 ◽  
Author(s):  
S. Chakraverty ◽  
Diptiranjan Behera

This chapter presents the static and dynamic analysis of structures with uncertain parameters using fuzzy finite element method. Uncertainties presents in the parameters are modelled through convex normalised fuzzy sets. Fuzzy finite element method converts the structures into fuzzy system of linear equations and fuzzy eigenvalue problem for static and dynamic problems respectively. As such method to solve fuzzy system of linear equations, fully fuzzy system of linear equations and fuzzy eigenvalue problems are presented. These methods are applied to various structural problems to find out the fuzzy static and dynamic responses of the structures. Also the chapter analyses the numerical solution of uncertain fractionally damped spring-mass system. Uncertainties considered in the initial condition of the system.


2020 ◽  
Vol 28 (4) ◽  
pp. 247-261
Author(s):  
Chunjae Park

AbstractRecently, the P1-nonconforming finite element space over square meshes has been proved stable to solve Stokes equations with the piecewise constant space for velocity and pressure, respectively. In this paper, we will introduce its locally divergence-free subspace to solve the elliptic problem for the velocity only decoupled from the Stokes equation. The concerning system of linear equations is much smaller compared to the Stokes equations. Furthermore, it is split into two smaller ones. After solving the velocity first, the pressure in the Stokes problem can be obtained by an explicit method very rapidly.


2020 ◽  
Vol 14 (3) ◽  
pp. 7141-7151 ◽  
Author(s):  
R. Omar ◽  
M. N. Abdul Rani ◽  
M. A. Yunus

Efficient and accurate finite element (FE) modelling of bolted joints is essential for increasing confidence in the investigation of structural vibrations. However, modelling of bolted joints for the investigation is often found to be very challenging. This paper proposes an appropriate FE representation of bolted joints for the prediction of the dynamic behaviour of a bolted joint structure. Two different FE models of the bolted joint structure with two different FE element connectors, which are CBEAM and CBUSH, representing the bolted joints are developed. Modal updating is used to correlate the two FE models with the experimental model. The dynamic behaviour of the two FE models is compared with experimental modal analysis to evaluate and determine the most appropriate FE model of the bolted joint structure. The comparison reveals that the CBUSH element connectors based FE model has a greater capability in representing the bolted joints with 86 percent accuracy and greater efficiency in updating the model parameters. The proposed modelling technique will be useful in the modelling of a complex structure with a large number of bolted joints.


Author(s):  
David Ek ◽  
Anders Forsgren

AbstractThe focus in this paper is interior-point methods for bound-constrained nonlinear optimization, where the system of nonlinear equations that arise are solved with Newton’s method. There is a trade-off between solving Newton systems directly, which give high quality solutions, and solving many approximate Newton systems which are computationally less expensive but give lower quality solutions. We propose partial and full approximate solutions to the Newton systems. The specific approximate solution depends on estimates of the active and inactive constraints at the solution. These sets are at each iteration estimated by basic heuristics. The partial approximate solutions are computationally inexpensive, whereas a system of linear equations needs to be solved for the full approximate solution. The size of the system is determined by the estimate of the inactive constraints at the solution. In addition, we motivate and suggest two Newton-like approaches which are based on an intermediate step that consists of the partial approximate solutions. The theoretical setting is introduced and asymptotic error bounds are given. We also give numerical results to investigate the performance of the approximate solutions within and beyond the theoretical framework.


2019 ◽  
Vol 2019 (1) ◽  
Author(s):  
A. Khalid ◽  
M. N. Naeem ◽  
P. Agarwal ◽  
A. Ghaffar ◽  
Z. Ullah ◽  
...  

AbstractIn the current paper, authors proposed a computational model based on the cubic B-spline method to solve linear 6th order BVPs arising in astrophysics. The prescribed method transforms the boundary problem to a system of linear equations. The algorithm we are going to develop in this paper is not only simply the approximation solution of the 6th order BVPs using cubic B-spline, but it also describes the estimated derivatives of 1st order to 6th order of the analytic solution at the same time. This novel technique has lesser computational cost than numerous other techniques and is second order convergent. To show the efficiency of the proposed method, four numerical examples have been tested. The results are described using error tables and graphs and are compared with the results existing in the literature.


2014 ◽  
Vol 2014 ◽  
pp. 1-8 ◽  
Author(s):  
Mohsen Alipour ◽  
Dumitru Baleanu ◽  
Fereshteh Babaei

We introduce a new combination of Bernstein polynomials (BPs) and Block-Pulse functions (BPFs) on the interval [0, 1]. These functions are suitable for finding an approximate solution of the second kind integral equation. We call this method Hybrid Bernstein Block-Pulse Functions Method (HBBPFM). This method is very simple such that an integral equation is reduced to a system of linear equations. On the other hand, convergence analysis for this method is discussed. The method is computationally very simple and attractive so that numerical examples illustrate the efficiency and accuracy of this method.


Sign in / Sign up

Export Citation Format

Share Document