Search results
1 – 10 of 632Francisco Bernal and Manuel Kindelan
The Motz problem can be considered as a benchmark problem for testing the performance of numerical methods in the solution of elliptic problems with boundary singularities. The…
Abstract
Purpose
The Motz problem can be considered as a benchmark problem for testing the performance of numerical methods in the solution of elliptic problems with boundary singularities. The purpose of this paper is to address the solution of the Motz problem using the radial basis function (RBF) method, which is a truly meshfree scheme.
Design/methodology/approach
Both the global RBF collocation method (also known as Kansa's method) and the recently proposed local RBF‐based differential quadrature (LRBFDQ) method are considered. In both cases, it is shown that the accuracy of the solution can be significantly increased by using special functions which capture the behavior of the singularity. In the case of global collocation, the functional space spanned by the RBF is enlarged by adding singular functions which capture the behavior of the local singular solution. In the case of local collocation, the problem is modified appropriately in order to eliminate the singularities from the formulation.
Findings
The paper shows that the exponential convergence both with increasing resolution and increasing shape parameter, which is typical of the RBF method, is lost in problems containing singularities. The accuracy of the solution can be increased by collocation of the partial differential equation (PDE) at boundary nodes. However, in order to restore the exponential convergence of the RBF method, it is necessary to use special functions which capture the behavior of the solution near the discontinuity.
Practical implications
The paper uses Motz's problem as a prototype for problems described by elliptic partial differential equations with boundary singularities. However, the results obtained in the paper are applicable to a wide range of problems containing boundaries with conditions which change from Dirichlet to Neumann, thus leading to singularities in the first derivatives.
Originality/value
The paper shows that both the global RBF collocation method and the LRBFDQ method, are truly meshless methods which can be very useful for the solution of elliptic problems with boundary singularities. In particular, when complemented with special functions that capture the behavior of the solution near the discontinuity, the method exhibits exponential convergence both with resolution and with shape parameter.
Details
Keywords
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.
Abstract
Purpose
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.
Design/methodology/approach
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.
Findings
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.
Originality/value
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.
Details
Keywords
S. D'Heedene, K. Amaratunga and J. Castrillón‐Candás
This paper presents a novel framework for solving elliptic partial differential equations (PDEs) over irregularly spaced meshes on bounded domains.
Abstract
Purpose
This paper presents a novel framework for solving elliptic partial differential equations (PDEs) over irregularly spaced meshes on bounded domains.
Design/methodology/approach
Second‐generation wavelet construction gives rise to a powerful generalization of the traditional hierarchical basis (HB) finite element method (FEM). A framework based on piecewise polynomial Lagrangian multiwavelets is used to generate customized multiresolution bases that have not only HB properties but also additional qualities.
Findings
For the 1D Poisson problem, we propose – for any given order of approximation – a compact closed‐form wavelet basis that block‐diagonalizes the stiffness matrix. With this wavelet choice, all coupling between the coarse scale and detail scales in the matrix is eliminated. In contrast, traditional higher‐order (n>1) HB do not exhibit this property. We also achieve full scale‐decoupling for the 2D Poisson problem on an irregular mesh. No traditional HB has this quality in 2D.
Research limitations/implications
Similar techniques may be applied to scale‐decouple the multiresolution finite element (FE) matrices associated with more general elliptic PDEs.
Practical implications
By decoupling scales in the FE matrix, the wavelet formulation lends itself particularly well to adaptive refinement schemes.
Originality/value
The paper explains second‐generation wavelet construction in a Lagrangian FE context. For 1D higher‐order and 2D first‐order bases, we propose a particular choice of wavelet, customized to the Poisson problem. The approach generalizes to other elliptic PDE problems.
Details
Keywords
Discusses the 27 papers in ISEF 1999 Proceedings on the subject of electromagnetisms. States the groups of papers cover such subjects within the discipline as: induction machines;…
Abstract
Discusses the 27 papers in ISEF 1999 Proceedings on the subject of electromagnetisms. States the groups of papers cover such subjects within the discipline as: induction machines; reluctance motors; PM motors; transformers and reactors; and special problems and applications. Debates all of these in great detail and itemizes each with greater in‐depth discussion of the various technical applications and areas. Concludes that the recommendations made should be adhered to.
Details
Keywords
Zohreh Heydari, Gholamreza Shobeyri and Seyed Hossein Ghoreishi Najafabadi
This paper aims to examine the accuracy of several higher-order incompressible smoothed particle hydrodynamics (ISPH) Laplacian models and compared with the classic model (Shao…
Abstract
Purpose
This paper aims to examine the accuracy of several higher-order incompressible smoothed particle hydrodynamics (ISPH) Laplacian models and compared with the classic model (Shao and Lo, 2003).
Design/methodology/approach
The numerical errors in solving two-dimensional elliptic partial differential equations using the Laplacian models are investigated for regular and highly irregular node distributions over a unit square computational domain.
Findings
The numerical results show that one of the Laplacian models, which is newly developed by one of the authors (Shobeyri, 2019) can get the smallest errors for various used node distributions.
Originality/value
The newly proposed model is formulated by the hybrid of the standard ISPH Laplacian model combined with Taylor expansion and moving least squares method. The superiority of the proposed model is significant when multi-resolution irregular node distributions commonly seen in adaptive refinement strategies used to save computational cost are applied.
Details
Keywords
Gives a bibliographical review of the error estimates and adaptive finite element methods from the theoretical as well as the application point of view. The bibliography at the…
Abstract
Gives a bibliographical review of the error estimates and adaptive finite element methods from the theoretical as well as the application point of view. The bibliography at the end contains 2,177 references to papers, conference proceedings and theses/dissertations dealing with the subjects that were published in 1990‐2000.
Details
Keywords
Thierry Maré, Nicolas Galanis, Sylvie Prétot and Jacques Miriel
To determine the axial evolution of the hydrodynamic and the thermal fields for mixed convection in inclined tubes and to investigate the presence of flow reversal.
Abstract
Purpose
To determine the axial evolution of the hydrodynamic and the thermal fields for mixed convection in inclined tubes and to investigate the presence of flow reversal.
Design/methodology/approach
The elliptical, coupled, steady state, three‐dimensional governing partial differential equations for heated ascending laminar mixed convection in an inclined isothermal tube were solved numerically using a finite volume staggered grid approach.
Findings
The axial evolution of the velocity profiles and fluid temperatures show that upstream diffusion has an important effect near the inlet of the heating region. As a result, both the wall shear stress and the Nusselt number are affected upstream of the heating zone. Flow reversal occurs of GF≥9 × 105. The shape and size of the region with negative velocities depends strongly on the value of the Grashof number. The effect of the Grashof number on the axial evolution of the wall shear stress and the Nusselt number is shown to be very important in the region of developing flow.
Research limitations/implications
The results have been calculated for one Reynolds number (Re=100), a single fluid (air) and one tube inclination (45°).
Practical implications
Further results of this type can be mapped and would be useful for heat exchanger design.
Originality/value
This is the first time that flow reversal has been calculated numerically for inclined tubes. Most previous studies concern horizontal or vertical tubes and use axially parabolic equations which are easier to solve but can not calculate the flow field in the region with backflow.
Details
Keywords
This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics…
Abstract
This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.
Details
Keywords
Gopal Priyadarshi and B.V. Rathish Kumar
In the past few years, Haar wavelet-based numerical methods have been applied successfully to solve linear and nonlinear partial differential equations. This study aims to propose…
Abstract
Purpose
In the past few years, Haar wavelet-based numerical methods have been applied successfully to solve linear and nonlinear partial differential equations. This study aims to propose a wavelet collocation method based on Haar wavelets to identify a parameter in parabolic partial differential equations (PDEs). As Haar wavelet is defined in a very simple way, implementation of the Haar wavelet method becomes easier than the other numerical methods such as finite element method and spectral method. The computational time taken by this method is very less because Haar matrices and Haar integral matrices are stored once and used for each iteration. In the case of Haar wavelet method, Dirichlet boundary conditions are incorporated automatically. Apart from this property, Haar wavelets are compactly supported orthonormal functions. These properties lead to a huge reduction in the computational cost of the method.
Design/methodology/approach
The aim of this paper is to reconstruct the source control parameter arises in quasilinear parabolic partial differential equation using Haar wavelet-based numerical method. Haar wavelets possess various properties, for example, compact support, orthonormality and closed form expression. The main difficulty with the Haar wavelet is its discontinuity. Therefore, this paper cannot directly use the Haar wavelet to solve partial differential equations. To handle this difficulty, this paper represents the highest-order derivative in terms of Haar wavelet series and using successive integration this study obtains the required term appearing in the problem. Taylor series expansion is used to obtain the second-order partial derivatives at collocation points.
Findings
An efficient and accurate numerical method based on Haar wavelet has been proposed for parameter identification in quasilinear parabolic partial differential equations. Numerical results are obtained from the proposed method and compared with the existing results obtained from various finite difference methods including Saulyev method. It is shown that the proposed method is superior than the conventional finite difference methods including Saulyev method in terms of accuracy and CPU time. Convergence analysis is presented to show the accuracy of the proposed method. An efficient algorithm is proposed to find the wavelet coefficients at target time.
Originality/value
The outcome of the paper would have a valuable role in the scientific community for several reasons. In the current scenario, the parabolic inverse problem has emerged as very important problem because of its application in many diverse fields such as tomography, chemical diffusion, thermoelectricity and control theory. In this paper, higher-order derivative is represented in terms of Haar wavelet series. In other words, we represent the solution in multiscale framework. This would enable us to understand the solution at various resolution levels. In the case of Haar wavelet, this paper can achieve a very good accuracy at very less resolution levels, which ultimately leads to huge reduction in the computational cost.
Details