Search results

Piezoelectric actuators and sensors are an invaluable part of lightweight designs for several reasons. They can either be used in noise cancellation devices as thin‐walled structures are prone to acoustic emissions, or in shape control approaches to suppress unwanted vibrations. Also in Lamb wave based health monitoring systems piezoelectric patches are applied to excite and to receive ultrasonic waves. The purpose of this paper is to develop a higher order finite element with piezoelectric capabilities in order to simulate smart structures efficiently.

In the paper the development of a new fully three‐dimensional piezoelectric hexahedral finite element based on the p‐version of the finite element method (FEM) is presented. Hierarchic Legendre polynomials in combination with an anisotropic ansatz space are utilized to derive an electro‐mechanically coupled element. This results in a reduced numerical effort. The suitability of the proposed element is demonstrated using various static and dynamic test examples.

In the current contribution it is shown that higher order coupled‐field finite elements hold several advantages for smart structure applications. All numerical examples have been found to agree well with previously published results. Furthermore, it is demonstrated that accurate results can be obtained with far fewer degrees of freedom compared to conventional low order finite element approaches. Thus, the proposed finite element can lead to a significant reduction in the overall numerical costs.

To the best of the author's knowledge, no piezoelectric finite element based on the hierarchical‐finite‐element‐method has yet been published in the literature. Thus, the proposed finite element is a step towards a holistic numerical treatment of structural health monitoring (SHM) related problems using p‐version finite elements.

The paper presents the principal elements of automatic adaptivity built in our 2D software for monolithic solution of multiphysics problems based on a fully adaptive finite element method of higher order of accuracy. The adaptive techniques are illustrated by appropriate examples.

Presented are algorithms for realization of the h‐adaptivity, p‐adaptivity, hp‐adaptivity, creation of curvilinear elements for modelling general boundaries and interfaces. Indicated also is the possibility of combining triangular and quadrilateral elements (both classical and curved).

The presented higher‐order adaptive processes are reliable, robust and lead to a substantial reduction of the degrees of freedom in comparison with the techniques used in low‐order finite element methods. They allow solving examples that are by classical approaches either unsolvable or solvable at a cost of high memory and time of computation.

The adaptive processes described in the paper are still limited to 2D computations. Their computer implementation is highly nontrivial (every physical field in a multiphysics task is generally solved on a different mesh satisfying its specific features) and in 3D the number of possible adaptive steps is many times higher.

The described adaptive techniques may represent a powerful tool for the monolithic solution of complex multiphysics problems.

The presented higher‐order adaptive approach of solution is shown to provide better results than the schemes implemented in professional codes based on low‐order finite element methods. Obtaining the results, moreover, requires less time and computer memory.

The purpose of this paper is to present a methodology of high‐precision finite element modeling of induction heating of rotating nonferromagnetic cylindrical billets in static magnetic field produced by appropriately arranged permanent magnets.

The mathematical model consisting of two partial differential equations describing the distribution of the magnetic and temperature fields are solved by a fully adaptive higher‐order finite element method in the monolithic formulation and selected results are validated experimentally.

The method of solution realized by own code is very fast, robust and exhibits much more powerful features when compared with classical low‐order numerical methods implemented in existing commercial codes.

For sufficiently long arrangements the method provides good results even for 2D model. The principal limitation consists in problems with determining correct boundary conditions for the temperature field (generalized coefficient of convective heat transfer as a function of the temperature and revolutions).

The methodology can successfully be used for design of devices for induction heating of cylindrical nonmagnetic bodies by rotation and determination of their operation parameters.

The paper is a presentation of the fully adaptive higher‐order finite element and its utilization for a monolithic numerical solution of a relatively complicated coupled problem.

The purpose of this paper is to propose a solution strategy for both accurate and efficient simulation of nonlinear magnetostatic problems in thin structures using higher order finite element methods. Special interest is put in the investigation of the step-lap joints of transformer cores, with a focus on the spatial resolution of the field quantities.

The usage of hierarchical finite elements of higher order makes it possible to adapt the local accuracy in different spatial directions in thin steel sheets. Due to explicit representation of gradients in the basis functions, a simple Schwarz-type block preconditioner with a conjugate gradient solver can efficiently solve the arising algebraic system. By adapting the block size automatically according to the aspect ratio, deterioration of convergence in case of thin elements can be prevented. The resulting Newton scheme is accelerated utilizing the hierarchical splitting in a two-level scheme, where an initial guess is computed on a coarse sub-space.

Compared to an isotropic choice of polynomial order for the basis functions, significant runtime and memory can be saved in the simulation of thin structures without losing accuracy. The iterative solution scheme proves to be robust with respect to the polynomial order, even for aspect ratios of 1:1000 and anisotropies in two directions. An additional saving in runtime and Newton iterations can be achieved by solving the nonlinear problem initially on the lowest order basis functions only and projecting the solution to the complete space as starting value, analogous to a full multigrid scheme.

Within the presented solution strategy, especially the anisotropic block preconditioner and the accelerated Newton scheme based on the two-level splitting constitute a novel contribution. They provide building blocks, which can be utilized for other types of magnetic field problems like transient nonlinear problems or hysteresis modeling as well.

The purpose of this paper is to present a model of induction heating of aluminium billets rotating in a static magnetic field generated by permanent magnets. The model is solved by the authors' own software and the results are verified experimentally.

The mathematical model of the problem given by two partial differential equations describing the distribution of the magnetic and temperature fields in the system is solved by a fully adaptive higher‐order finite element method in the hard‐coupled formulation. All material nonlinearities are taken into account.

The method of solution realized by the code is reliable and works faster in comparison with the existing low‐order finite element codes.

The method works for 2D arrangements with an extremely high accuracy. Its limitations consist mainly in problems of determining the coefficients of convection and radiation for temperature field in the system (respecting both temperature and revolutions).

The methodology can successfully be used for design of devices for induction heating of cylindrical nonmagnetic bodies by rotation and anticipation of their operation parameters.

The paper presents a fully adaptive higher‐order finite element and its utilization for a hard‐coupled numerical solution of the problem of induction heating.

The purpose of this paper is to present an application of the Generalized Finite Element Method (GFEM) for modal analysis of 2D wave equation.

The GFEM can be viewed as an extension of the standard Finite Element Method (FEM) that allows non-polynomial enrichment of the approximation space. In this paper the authors enrich the approximation space with sine e cosine functions, since these functions frequently appear in the analytical solution of the problem under study. The results are compared with the ones obtained with the polynomial FEM using higher order elements.

The results indicate that the proposed approach is able to obtain more accurate results for higher vibration modes than standard polynomial FEM.

The examples studied in this paper indicate a strong potential of the GFEM for the approximation of higher vibration modes of structures, analysis of structures subject to high frequency excitations and other problems that concern high frequency oscillatory phenomena.

This paper aims to present a high-order algorithm implemented with the modal spectral element method and simulations of three-dimensional thermal convective flows by using the full viscous dissipation function in the energy equation. Three benchmark problems were solved to validate the algorithm with exact or theoretical solutions. The heated rotating sphere at different temperatures inside a cold planar Poiseuille flow was simulated parametrically at varied angular velocities with positive and negative rotations.

The fourth-order stiffly stable schemes were implemented and tested for time integration. To provide the hp-refinement and spatial resolution enhancement, a modal spectral element method using hierarchical basis functions was used to solve governing equations in a three-dimensional space.

It was found that the direction of rotation of the heated sphere has totally different effects on drag, lateral force and torque evaluated on surfaces of the sphere and walls. It was further concluded that the angular velocity of the heated sphere has more influence on the wall normal velocity gradient than on the wall normal temperature gradients and therefore, more influence on the viscous dissipation than on the thermal dissipation.

This paper concerns incompressible fluid flow at constant properties with up to medium temperature variations in the absence of thermal radiation and ignoring the pressure work.

This paper contributes a viable high-order algorithm in time and space for modeling convective heat transfer involving an internal heated rotating sphere with the effect of viscous heating.

Results of this paper could provide reference for related topics such as enhanced heat transfer forced convection involving rotating spheres and viscous thermal effect.

The merits include resolving viscous dissipation and thermal diffusion in stationary and rotating boundary layers with both h- and p-type refinements, visualizing the viscous heating effect with the full viscous dissipation function in the energy equation and modeling the forced advection around a rotating sphere with varied positive and negative angular velocities subject to a shear flow.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional time‐harmonic electromagnetic wave propagation problems.

The proposed HDG method for the discretization of the two‐dimensional transverse magnetic Maxwell equations relies on an arbitrary high order nodal interpolation of the electromagnetic field components and is formulated on triangular meshes. In the HDG method, an additional hybrid variable is introduced on the faces of the elements, with which the element‐wise (local) solutions can be defined. A so‐called conservativity condition is imposed on the numerical flux, which can be defined in terms of the hybrid variable, at the interface between neighbouring elements. The linear system of equations for the unknowns associated with the hybrid variable is solved here using a multifrontal sparse LU method. The formulation is given, and the relationship between the considered HDG method and a standard upwind flux‐based DG method is also examined.

The approximate solutions for both electric and magnetic fields converge with the optimal order of p+1 in L₂ norm, when the interpolation order on every element and every interface is p and the sought solution is sufficiently regular. The presented numerical results show the effectiveness of the proposed HDG method, especially when compared with a classical upwind flux‐based DG method.

The work described here is a demonstration of the viability of a HDG formulation for solving the time‐harmonic Maxwell equations through a detailed numerical assessment of accuracy properties and computational performances.

This study seeks to focus on the annular flow between rectangular and equilateral‐triangular ducts under all possible arrangements. The aim of this work is to obtain accurate prediction of the friction factor of this flow using high‐order finite element method.

Steady and fully developed laminar flow of incompressible Newtonian fluid in an annulus of variable cross‐sectional geometry is investigated numerically. Accurate prediction of the friction factor of this flow was obtained using high‐order finite element method.

The results were in agreement with already published findings in the literature. It was found that a higher annular area ratio will lead to a monotonic increase in fRe value in the case of regular annuli, and will lead to an increase followed by a decrease in fRe value in the case of irregular annuli. Also, it was, found that irregular annuli have lower fRe value than regular annuli, and that the square‐in‐triangle case has the lowest fRe value, whereas the square‐in‐square case has the highest fRe value.

Accurate prediction of the friction factor of the laminar flow in irregular annuli was obtained. Also, the obtained results can be utilized to optimize the annular geometries under consideration. In addition, the obtained results can lead to the design of more efficient heat exchangers.