Search results

This study aims to overcome the involved challenging issues and provide high-precision eigensolutions. General eigenproblems in the system of ordinary differential equations (ODEs) serve as mathematical models for vector Sturm-Liouville (SL) and free vibration problems. High-precision eigenvalue and eigenfunction solutions are crucial bases for the reliable dynamic analysis of structures. However, solutions that meet the error tolerances specified are difficult to obtain for issues such as coefficients of variable matrices, coincident and adjacent approximate eigenvalues, continuous orders of eigenpairs and varying boundary conditions.

This study presents an h-version adaptive finite element method based on the superconvergent patch recovery displacement method for eigenproblems in system of second-order ODEs. The high-order shape function interpolation technique is further introduced to acquire superconvergent solution of eigenfunction, and superconvergent solution of eigenvalue is obtained by computing the Rayleigh quotient. Superconvergent solution of eigenfunction is used to estimate the error of finite element solution in the energy norm. The mesh is then, subdivided to generate an improved mesh, based on the error.

Representative eigenproblems examples, containing typical vector SL and free vibration of beams problems involved the aforementioned challenging issues, are selected to evaluate the accuracy and reliability of the proposed method. Non-uniform refined meshes are established to suit eigenfunctions change, and numerical solutions satisfy the pre-specified error tolerance.

The proposed combination of methodologies described in the paper, leads to a powerful h-version mesh refinement algorithm for eigenproblems in system of second-order ODEs, that can be extended to other classes of applications in damage detection of multiple cracks in structures based on the high-precision eigensolutions.

The purpose of this paper is to propose the generalized integral transform technique (GITT) to the solution of convection-diffusion problems with nonlinear boundary conditions by employing the corresponding nonlinear eigenvalue problem in the construction of the expansion basis.

The original nonlinear boundary condition coefficients in the problem formulation are all incorporated into the adopted eigenvalue problem, which may be itself integral transformed through a representative linear auxiliary problem, yielding a nonlinear algebraic eigenvalue problem for the associated eigenvalues and eigenvectors, to be solved along with the transformed ordinary differential system. The nonlinear eigenvalues computation may also be accomplished by rewriting the corresponding transcendental equation as an ordinary differential system for the eigenvalues, which is then simultaneously solved with the transformed potentials.

An application on one-dimensional transient diffusion with nonlinear boundary condition coefficients is selected for illustrating some important computational aspects and the convergence behavior of the proposed eigenfunction expansions. For comparison purposes, an alternative solution with a linear eigenvalue problem basis is also presented and implemented.

This novel approach can be further extended to various classes of nonlinear convection-diffusion problems, either already solved by the GITT with a linear coefficients basis, or new challenging applications with more involved nonlinearities.

The purpose of this paper is to present a numerically adaptive finite element (FE) method for accurate, efficient and reliable eigensolutions of regular second- and fourth-order Sturm–Liouville (SL) problems with variable coefficients.

After the conventional FE solution for an eigenpair (i.e. eigenvalue and eigenfunction) of a particular order has been obtained on a given mesh, a novel strategy is introduced, in which the FE solution of the eigenproblem is equivalently viewed as the FE solution of an associated linear problem. This strategy allows the element energy projection (EEP) technique for linear problems to calculate the super-convergent FE solutions for eigenfunctions anywhere on any element. These EEP super-convergent solutions are used to estimate the FE solution errors and to guide mesh refinements, until the accuracy matches user-preset error tolerance on both eigenvalues and eigenfunctions.

Numerical results for a number of representative and challenging SL problems are presented to demonstrate the effectiveness, efficiency, accuracy and reliability of the proposed method.

The method is limited to regular SL problems, but it can also solve some singular SL problems in an indirect way.

Comprehensive utilization of the EEP technique yields a simple, efficient and reliable adaptive FE procedure that finds sufficiently fine meshes for preset error tolerances on eigenvalues and eigenfunctions to be achieved, even on problems which proved troublesome to competing methods. The method can readily be extended to vector SL problems.

A general solution for the small deflexions of thin plates of slowly varying thickness under lateral loading in the form of an influence function is briefly presented. It is known that the influence function may be represented as an infinite series in terms of the eigenfunctions and eigenvalues associated with a homogeneous form of the plate differential equation. It is suggested that the series may give an acceptable approximation to the influence function when summed over a small number of terms when also the eigenfunctions and eigenvalues involved are deduced by an approximate procedure of the Rayleigh‐Ritz type. In order to test this assertion a numerical example is given for a uniform canti‐lever plate and the results are compared with experiment and with similar results deduced by an alternative theoretical procedure. Thus the calculation of a sufficient number of approximate normal vibration modes and frequencies for the plate as normally required for aeroelastic investigations may in this way be made to serve as the basis for a complete analysis of the plate. A simple approximate allowance for shear deflexion of the plate is presented and illustrated.

This paper aims to present a meshless technique to find the Green’s functions for solutions of Laplacian boundary value problems on rectangular domains. This paper also investigates a theoretical basis for the Steklov series expansion methods to reduce and estimate the error of numerical approaches for the boundary correction kernel of the Laplace operator.

The main interest is how the Green's functions differ from the fundamental solution of the Laplace operator. Steklov expansion methods for finding the correction term are supported by the analysis that bases of the class of all finite harmonic functions can be formed using harmonic Steklov eigenfunctions. These functions construct a basis of the space of solutions of harmonic boundary value problems and their boundary traces generate an orthogonal basis of the trace space of solutions on the boundary.

The main conclusion is that the boundary correction term for the Green's functions is well-approximated by Steklov expansions with a few Steklov eigenfunctions. The error estimates for the Steklov approximations of the boundary correction term involved in Dirichlet or Robin boundary value problems are found. They appear to provide very good approximations in the interior of the region and become quite oscillatory close to the boundary.

This paper concentrates to document the first attempt to find the Green's function for various harmonic boundary value problems with the explicit Steklov eigenfunctions without concerns regarding discretizations when the region is a rectangle.

The purpose of this paper is to present a computationally efficient approach for the stochastic modeling of an inhomogeneous reluctivity of magnetic materials. These materials can be part of electrical machines such as a single-phase transformer (a benchmark example that is considered in this paper). The approach is based on the Karhunen–Loève expansion (KLE). The stochastic model is further used to study the statistics of the self-inductance of the primary coil as a quantity of interest (QoI).

The computation of the KLE requires solving a generalized eigenvalue problem with dense matrices. The eigenvalues and the eigenfunction are computed by using the Lanczos method that needs only matrix vector multiplications. The complexity of performing matrix vector multiplications with dense matrices is reduced by using hierarchical matrices.

The suggested approach is used to study the impact of the spatial variability in the magnetic reluctivity on the QoI. The statistics of this parameter are influenced by the correlation lengths of the random reluctivity. Both, the mean value and the standard deviation increase as the correlation length of the random reluctivity increases.

The KLE, computed by using hierarchical matrices, is used for uncertainty quantification of low frequency electrical machines as a computationally efficient approach in terms of memory requirement, as well as computation time.

The purpose of this study is to propose the generalised integral transform technique to investigate the natural convection behaviour in a vertical cylinder under different boundary conditions, adiabatic and isothermal walls and various aspect ratios.

GITT was used to investigate the steady-state natural convection behaviour in a vertical cylinder with internal uniformed heat generation. The governing equations of natural convection were transferred to a set of ordinary differential equations by using the GITT methodology. The coefficients of the ODEs were determined by the integration of the eigenfunction of the auxiliary eigenvalue problems in the present natural convection problem. The ordinary differential equations were solved numerically by using the DBVPFD subroutine from the IMSL numerical library. The convergence was achieved reasonably by using low truncation orders.

GITT is a powerful computational tool to explain the convection phenomena in the cylindrical cavity. The convergence analysis shows that the hybrid analytical–numerical technique (GITT) has a good convergence performance in relatively low truncation orders in the stream-function and temperature fields. The effect of the Rayleigh number and aspect ratio on the natural convection behaviour under adiabatic and isothermal boundary conditions has been discussed in detail.

The present hybrid analytical–numerical methodology can be extended to solve various convection problems with more involved nonlinearities. It exhibits potential application to solve the convection problem in the nuclear, oil and gas industries.

The purpose of this study is to develop an efficient model to solve the electromagnetic forward problem using a novel semi-analytical approach to compute the electromagnetic fields because of the presence of a scatterer.

The proposed model involves a novel formulation of a complete orthonormal set of radiating/nonradiating polarization currents. Furthermore, an integral equation-based representation is derived, and the appropriate boundary conditions are imposed to get the scattered electromagnetic field. An error term is introduced to evaluate the obtained results.

The proposed model was tested using several examples at different frequencies. The results of this study show that the novel representation exhibits fast convergence behavior and achieves highly accurate results, when compared to the results provided by the transmission line method.

The derived formulations presented in this study are significant in the electromagnetic forward modelling field because of the meaningful physical representation they provide. This is an important aspect that leads to precise calculation of electromagnetic fields for various applications.

Different possibilities for the enhancement of convergence rates in eigenfunction expansions are investigated in the realm of integral transform solutions for partial differential equations. A representative parabolic problem is chosen to illustrate two schemes and their combinations; a filtering technique and an integral balance approach. Numerical results are presented to confirm the relative merits in each proposed procedure.

This study aims to develop an adaptive finite element method for structural eigenproblems of cracked Euler–Bernoulli beams via the superconvergent patch recovery displacement technique. This research comprises the numerical algorithm and experimental results for free vibration problems (forward eigenproblems) and damage detection problems (inverse eigenproblems).

The weakened properties analogy is used to describe cracks in this model. The adaptive strategy proposed in this paper provides accurate, efficient and reliable eigensolutions of frequency and mode (i.e. eigenpairs as eigenvalue and eigenfunction) for Euler–Bernoulli beams with multiple cracks. Based on the frequency measurement method for damage detection, using the difference between the actual and computed frequencies of cracked beams, the inverse eigenproblems are solved iteratively for identifying the residuals of locations and sizes of the cracks by the Newton–Raphson iteration technique. In the crack detection, the estimated residuals are added to obtain reliable results, which is an iteration process that will be expedited by more accurate frequency solutions based on the proposed method for free vibration problems.

Numerical results are presented for free vibration problems and damage detection problems of representative non-uniform and geometrically stepped Euler–Bernoulli beams with multiple cracks to demonstrate the effectiveness, efficiency, accuracy and reliability of the proposed method.

The proposed combination of methodologies described in the paper leads to a very powerful approach for free vibration and damage detection of beams with cracks, introducing the mesh refinement, that can be extended to deal with the damage detection of frame structures.