Search results

In the present work, we focus on developing an in-house parallel meshless local Petrov-Galerkin (MLPG) code for the analysis of heat conduction in two-dimensional and three-dimensional regular as well as complex geometries.

The parallel MLPG code has been implemented using open multi-processing (OpenMP) application programming interface (API) on the shared memory multicore CPU architecture. Numerical simulations have been performed to find the critical regions of the serial code, and an OpenMP-based parallel MLPG code is developed, considering the critical regions of the sequential code.

Based on performance parameters such as speed-up and parallel efficiency, the credibility of the parallelization procedure has been established. Maximum speed-up and parallel efficiency are 10.94 and 0.92 for regular three-dimensional geometry (343,000 nodes). Results demonstrate the suitability of parallelization for larger nodes as parallel efficiency and speed-up are more for the larger nodes.

Few attempts have been made in parallel implementation of the MLPG method for solving large-scale industrial problems. Although the literature suggests that message-passing interface (MPI) based parallel MLPG codes have been developed, the OpenMP model has rarely been touched. This work is an attempt at the development of OpenMP-based parallel MLPG code for the very first time.

The study aims to highlight the behaviour of one-dimensional and two-dimensional fin models under the natural room conditions, considering the different values of dimensionless Biot number (Bi). The effect of convection and radiation on the heat transfer process has also been demonstrated using the meshless local Petrov–Galerkin (MLPG) approach.

It is true that MLPG method is time-consuming and expensive in terms of man-hours, as it is in the developing stage, but with the advent of computationally fast new-generation computers, there is a big possibility of the development of MLPG software, which will not only reduce the computational time and cost but also enhance the accuracy and precision in the results. Bi values of 0.01 and 0.10 have been taken for the experimental investigation of one-dimensional and two-dimensional rectangular fin models. The numerical simulation results obtained by the analytical method, benchmark numerical method and the MLPG method for both the models have been compared with that of the experimental investigation results for validation and found to be in good agreement. Performance of the fin has also been demonstrated.

The experimental and numerical investigations have been conducted for one-dimensional and two-dimensional linear and nonlinear fin models of rectangular shape. MLPG is used as a potential numerical method. Effect of radiation is also, implemented successfully. Results are found to be in good agreement with analytical solution, when one-dimensional steady problem is solved; however, two-dimensional results obtained by the MLPG method are compared with that of the finite element method and found that the proposed method is as accurate as the established method. It is also found that for higher Bi, the one-dimensional model is not appropriate, as it does not demonstrate the appreciated error; hence, a two-dimensional model is required to predict the performance of a fin. Radiative fin illustrates more heat transfer than the pure convective fin. The performance parameters show that as the Bi increases, the performance of fin decreases because of high thermal resistance.

Though, best of the efforts have been put to showcase the behaviour of one-dimensional and two-dimensional fins under nonlinear conditions, at different Bi values, yet lot more is to be demonstrated. Nonlinearity, in the present paper, is exhibited by using the thermal and material properties as the function of temperature, but can be further demonstrated with their dependency on the area. Additionally, this paper can be made more elaborative by extending the research for transient problems, with different fin profiles. Natural convection model is adopted in the present study but it can also be studied by using forced convection model.

Fins are the most commonly used medium to enhance heat transfer from a hot primary surface. Heat transfer in its natural condition is nonlinear and hence been demonstrated. The outcome is practically viable, as it is applicable at large to the broad areas like automobile, aerospace and electronic and electrical devices.

As per the literature survey, lot of work has been done on fins using different numerical methods; but to the best of authors’ knowledge, this study is first in the area of nonlinear heat transfer of fins using dimensionless Bi by the truly meshfree MLPG method.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

The purpose of this paper is to assess the performance of the stabilised moving least squares (MLS) scheme in the meshless local Petrov–Galerkin (MLPG) method for heat conduction method.

In the current work, the authors extend the stabilised MLS approach to the MLPG method for heat conduction problem. Its performance has been compared with the MLPG method based on the standard MLS and local coordinate MLS. The patch tests of MLS and modified MLS schemes have been presented along with the one- and two-dimensional examples for MLPG method of the heat conduction problem.

In the stabilised MLS, the condition number of moment matrix is independent of the nodal spacing and it is nearly constant in the global domain for all grid sizes. The shifted polynomials based MLS and stabilised MLS approaches are more robust than the standard MLS scheme in the MLPG method analysis of heat conduction problems.

The MLPG method based on the stabilised MLS scheme.

The work presents a novel implementation of the generalized minimum residual (GMRES) solver in conjunction with the interpolating meshless local Petrov–Galerkin (MLPG) method to solve steady-state heat conduction in 2-D as well as in 3-D domains.

The restarted version of the GMRES solver (with and without preconditioner) is applied to solve an asymmetric system of equations, arising due to the interpolating MLPG formulation. Its performance is compared with the biconjugate gradient stabilized (BiCGSTAB) solver on the basis of computation time and convergence behaviour. Jacobi and successive over-relaxation (SOR) methods are used as the preconditioners in both the solvers.

The results show that the GMRES solver outperforms the BiCGSTAB solver in terms of smoothness of convergence behaviour, while performs slightly better than the BiCGSTAB method in terms of Central processing Unit (CPU) time.

MLPG formulation leads to a non-symmetric system of algebraic equations. Iterative methods such as GMRES and BiCGSTAB methods are required for its solution for large-scale problems. This work presents the use of GMRES solver with the MLPG method for the very first time.

The purpose of this paper is to solve both eigenvalue and boundary value problems coming from the field of quantum mechanics through the application of meshless methods, particularly the one known as meshless local Petrov‐Galerkin (MLPG).

Regarding eigenvalue problems, the authors show how to apply MLPG to the time‐independent Schrödinger equation stated in three dimensions. Through a special procedure, the numerical integration of weak forms is carried out only for internal nodes. The boundary conditions are enforced through a collocation method. The final result is a generalized eigenvalue problem, which is readily solved for the energy levels. An example of boundary value problem is described by the time‐dependent nonlinear Schrödinger equation. The weak forms are again stated only for internal nodes, whereas the same collocation scheme is employed to enforce the boundary conditions. The nonlinearity is dealt with by a predictor‐corrector scheme.

Results show that the combination of MLPG and a collocation scheme works very well. The numerical results are compared to those provided by analytical solutions, exhibiting good agreement.

The flexibility of MLPG is made explicit. There are different ways to obtain the weak forms, and the boundary conditions can be enforced through a number of ways, the collocation scheme being just one of them. The shape functions used to approximate the solution can incorporate modifications that reflect some feature of the problem, like periodic boundary conditions. The value of this work resides in the fact that problems from other areas of electromagnetism can be attacked by the very same ideas herein described.

This paper aims to study nonlinear heat transfer through a longitudinal fin of three different profiles.

A truly meshfree method is used to undertake a nonlinear analysis to predict temperature distribution and heat-transfer rate.

A longitudinal fin of three different profiles, such as rectangular, triangular and concave parabolic, are analyzed. Temperature variation, along with the fin length and rate of heat transfer in steady state, under convective and convective-radiative environments has been demonstrated and explained. Moving least square (MLS) approximants are used to approximate the unknown function of temperature T(x) with Th(x). Essential boundary conditions are imposed using the penalty method. An iterative predictor–corrector scheme is used to handle nonlinearity.

Modelling fin in a convective-radiative environment removes the assumption of no radiation condition. It also allows to vary convective heat-transfer coefficient and predict the closer values to the real problems for the corresponding fin surfaces.

The meshless local Petrov–Galerkin method can solve nonlinear fin problems and predict an accurate solution.

The problems of shear flexible beams are analyzed by the MLPG method based on a locking‐free weak formulation. In order for the weak formulation to be locking‐free, the numerical characteristics of the variational functional for a shear flexible beam, in the thin beam limit, are discussed. Based on these discussions a locking‐free local symmetric weak form is derived by changing the set of two dependent variables in governing equations from that of transverse displacement and total rotation to the set of transverse displacement and transverse shear strain. For the interpolation of the chosen set of dependent variables (i.e. transverse displacement and transverse shear strain) in the locking‐free local symmetric weak form, the recently proposed generalized moving least squares (GMLS) interpolation scheme is utilized, in order to introduce the derivative of the transverse displacement as an additional nodal degree of freedom, independent of the nodal transverse displacement. Through numerical examples, convergence tests are performed. To identify the locking‐free nature of the proposed method, problems of shear flexible beams in the thick beam limit and in the thin beam limit are analyzed, and the numerical results are compared with analytical solutions. The potential of using the truly meshless local Petrov‐Galerkin (MLPG) method is established as a new paradigm in totally locking‐free computational analyses of shear flexible plates and shells.

The purpose of this paper is to present a micromechanical model based on a new truly local meshless method for analysis of heat transfer in composite materials.

The presented meshless method is based on the integral form of energy equation in the sub‐particles in the material. In the presented meshless method due to elimination of domain integration the computational efforts are decreased substantially.

Numerical results are presented for temperature distribution, heat flux and thermal conductivity. Numerical results show that the presented meshless method is simple, effective, accurate and less costly method in micromechanical modeling of heat conduction in heterogeneous materials.

A small area of the composite system called representative volume element is considered as the solution domain. The fully bonded fiber‐matrix interface is considered and contact thermal resistant is neglected from the fiber matrix interface and so the continuity of temperature and reciprocity of heat flux is satisfied in the fiber‐matrix interface.

For the first time a new truly local meshless method based on the integral form of energy equation for the sub‐particles in the materials is presented for analysis of heat transfer in composite materials.

This paper is concerned with new formulations of local meshfree and finite element numerical methods, for the solution of two-dimensional problems in linear elasticity.

In the local domain, assigned to each node of a discretization, the work theorem establishes an energy relationship between a statically admissible stress field and an independent kinematically admissible strain field. This relationship, derived as a weighted residual weak form, is expressed as an integral local form. Based on the independence of the stress and strain fields, this local form of the work theorem is kinematically formulated with a simple rigid-body displacement to be applied by local meshfree and finite element numerical methods. The main feature of this paper is the use of a linearly integrated local form that implements a quite simple algorithm with no further integration required.

The reduced integration, performed by this linearly integrated formulation, plays a key role in the behavior of local numerical methods, since it implies a reduction of the nodal stiffness which, in turn, leads to an increase of the solution accuracy and, which is most important, presents no instabilities, unlike nodal integration methods without stabilization. As a consequence of using such a convenient linearly integrated local form, the derived meshfree and finite element numerical methods become fast and accurate, which is a feature of paramount importance, as far as computational efficiency of numerical methods is concerned. Three benchmark problems were analyzed with these techniques, in order to assess the accuracy and efficiency of the new integrated local formulations of meshfree and finite element numerical methods. The results obtained in this work are in perfect agreement with those of the available analytical solutions and, furthermore, outperform the computational efficiency of other methods. Thus, the accuracy and efficiency of the local numerical methods presented in this paper make this a very reliable and robust formulation.

Presentation of a new local mesh-free numerical method. The method, linearly integrated along the boundary of the local domain, implements an algorithm with no further integration required. The method is absolutely reliable, with remarkably-accurate results. The method is quite robust, with extremely-fast computations.