Search results

1 – 10 of 114
To view the access options for this content please click here
Article
Publication date: 4 June 2020

Tiago Oliveira, Wilber Vélez and Artur Portela

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.

Abstract

Purpose

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.

Design/methodology/approach

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.

Findings

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.

Originality/value

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.

Details

Multidiscipline Modeling in Materials and Structures, vol. 16 no. 5
Type: Research Article
ISSN: 1573-6105

Keywords

To view the access options for this content please click here
Article
Publication date: 15 November 2011

Emre Erkmen and M.A. Bradford

The purpose of this paper is to develop a computational technique to couple finite element and meshfree methods for locking‐free analysis of shear deformable beams and…

Abstract

Purpose

The purpose of this paper is to develop a computational technique to couple finite element and meshfree methods for locking‐free analysis of shear deformable beams and plates, and to impose the boundary conditions directly when the matching field approach is adopted in the meshfree region.

Design/methodology/approach

Matching field approach eliminates shear‐locking which may occur due to inconsistencies in the approximations of the transverse displacement and rotation fields in shear‐deformable beams and plates. Continuous blending method is modified in order to be able to satisfy the constraint conditions of the matching field strategy.

Findings

For both transverse displacement and rotation fields, the developed technique produces approximation functions that satisfy the Kronecker delta property at the required nodes of the meshfree region when the matching field approach is adopted.

Originality/value

This approach allows for direct assembly of the stiffness matrices that are built for separate finite element and meshfree regions when the matching field approach is adopted. The boundary conditions can be directly applied, and the reaction forces can also be calculated directly from the structural stiffness matrix by using the developed technique.

Details

Engineering Computations, vol. 28 no. 8
Type: Research Article
ISSN: 0264-4401

Keywords

To view the access options for this content please click here
Article
Publication date: 20 January 2021

Ram Jiwari and Alf Gerisch

This paper aims to develop a meshfree algorithm based on local radial basis functions (RBFs) combined with the differential quadrature (DQ) method to provide numerical…

Abstract

Purpose

This paper aims to develop a meshfree algorithm based on local radial basis functions (RBFs) combined with the differential quadrature (DQ) method to provide numerical approximations of the solutions of time-dependent, nonlinear and spatially one-dimensional reaction-diffusion systems and to capture their evolving patterns. The combination of local RBFs and the DQ method is applied to discretize the system in space; implicit multistep methods are subsequently used to discretize in time.

Design/methodology/approach

In a method of lines setting, a meshless method for their discretization in space is proposed. This discretization is based on a DQ approach, and RBFs are used as test functions. A local approach is followed where only selected RBFs feature in the computation of a particular DQ weight.

Findings

The proposed method is applied on four reaction-diffusion models: Huxley’s equation, a linear reaction-diffusion system, the Gray–Scott model and the two-dimensional Brusselator model. The method captured the various patterns of the models similar to available in literature. The method shows second order of convergence in space variables and works reliably and efficiently for the problems.

Originality/value

The originality lies in the following facts: A meshless method is proposed for reaction-diffusion models based on local RBFs; the proposed scheme is able to capture patterns of the models for big time T; the scheme has second order of convergence in both time and space variables and Nuemann boundary conditions are easy to implement in this scheme.

Details

Engineering Computations, vol. ahead-of-print no. ahead-of-print
Type: Research Article
ISSN: 0264-4401

Keywords

To view the access options for this content please click here
Article
Publication date: 1 May 2020

Rajul Garg, Harishchandra Thakur and Brajesh Tripathi

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…

Abstract

Purpose

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.

Design/methodology/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.

Findings

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.

Research limitations/implications

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.

Practical implications

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.

Originality/value

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.

Details

Engineering Computations, vol. 37 no. 8
Type: Research Article
ISSN: 0264-4401

Keywords

To view the access options for this content please click here
Article
Publication date: 13 May 2019

Ram Jiwari, Sanjay Kumar and R.C. Mittal

The purpose of this paper is to develop two meshfree algorithms based on multiquadric radial basis functions (RBFs) and differential quadrature (DQ) technique for…

Abstract

Purpose

The purpose of this paper is to develop two meshfree algorithms based on multiquadric radial basis functions (RBFs) and differential quadrature (DQ) technique for numerical simulation and to capture the shocks behavior of Burgers’ type problems.

Design/methodology/approach

The algorithms convert the problems into a system of ordinary differential equations which are solved by the Runge–Kutta method.

Findings

Two meshfree algorithms are developed and their stability is discussed. Numerical experiment is done to check the efficiency of the algorithms, and some shock behaviors of the problems are presented. The proposed algorithms are found to be accurate, simple and fast.

Originality/value

The present algorithms LRBF-DQM and GRBF-DQM are based on radial basis functions, which are new for Burgers’ type problems. It is concluded from the numerical experiments that LRBF-DQM is better than GRBF-DQM. The algorithms give better results than available literature.

To view the access options for this content please click here
Article
Publication date: 13 November 2007

Gleber Nelson Marques, José Márcio Machado, Sérgio Luis Lopes Verardi, Stephan Stephany and Airam Jonatas Preto

This paper proposes an interpolating approach of the element‐free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value…

Abstract

Purpose

This paper proposes an interpolating approach of the element‐free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non‐homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.

Design/methodology/approach

The authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.

Findings

The local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM‐oriented programs can, thus, be easily extended to provide EFGM approximations.

Research limitations/implications

The robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high‐material discontinuity.

Practical implications

The proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.

Originality/value

This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well as of a penalty strategy. The resulting formulation provided accurate results including in the case of high‐material discontinuity.

Details

COMPEL - The international journal for computation and mathematics in electrical and electronic engineering, vol. 26 no. 5
Type: Research Article
ISSN: 0332-1649

Keywords

To view the access options for this content please click here
Article
Publication date: 7 October 2019

Wenan Wu and Hong Zheng

This study aims to introduce the hybrid finite element (FE) – meshfree method and multiscale variational principle into the traditional mixed FE formulation, leading to a…

Abstract

Purpose

This study aims to introduce the hybrid finite element (FE) – meshfree method and multiscale variational principle into the traditional mixed FE formulation, leading to a stable mixed formulation for incompressible linear elasticity which circumvents the need to satisfy inf-sup condition.

Design/methodology/approach

Using the hybrid FE–meshfree method, the displacement and pressure are interpolated conveniently with the same order so that a continuous pressure field can be obtained with low-order elements. The multiscale variational principle is then introduced into the Galerkin form to obtain stable and convergent results.

Findings

The present method is capable of overcoming volume locking and does not exhibit unphysical oscillations near the incompressible limit. Moreover, there are no extra unknowns introduced in the present method because the fine-scale unknowns are eliminated using the static condensation technique, and there is no need to evaluate any user-defined stability parameter as the classical stabilization methods do. The shape functions constructed in the present model possess continuous derivatives at nodes, which gives a continuous and more precise stress field with no need of an additional smooth process. The shape functions in the present model also possess the Kronecker delta property, so that it is convenient to impose essential boundary conditions.

Originality/value

The proposed model can be implemented easily. Its convergence rates and accuracy in displacement, energy and pressure are even comparable to those of second-order mixed elements.

To view the access options for this content please click here
Article
Publication date: 15 June 2015

Chun Feng, Shi-hai Li and Eugenio Onate

Continuum-based discrete element method is an explicit numerical method, which is a combination of block discrete element method (DEM) and FEM. When simulating large…

Abstract

Purpose

Continuum-based discrete element method is an explicit numerical method, which is a combination of block discrete element method (DEM) and FEM. When simulating large deformation problems, such as cutting, blasting, water-like material flowing, the distortion of elements will lead to no convergence of the numerical system. To solve the convergence problem, a particle contact-based meshfree method (PCMM) is introduced in. The paper aims to discuss this issue.

Design/methodology/approach

PCMM is based on traditional particle DEM, and use particle contacts to generate triangular elements. If three particles are contact with each other, the element will be created. Once elements are created, the macroscopic constitutive law could be introduced in. When large deformation of element occurs, the contact relationship between particles will be changed. Those elements that do not meet the contact condition will be deleted, and new elements that coincide with the relationship will be generated. By the deletion and creation of elements, the convergence problem induced by element distortion will be eliminated. To solve FEM and PCMM coupled problems, a point-edge contact model is introduced in, and normal and tangential springs are adopted to transfer the contact force between particles and blocks.

Findings

According to the deletion and recreation of elements based on particle contacts, PCMM could simulate large deformation problems. Some numerical cases (i.e. elastic field testing, uniaxial compression analysis and wave propagation simulation) show the accuracy of PCMM, and others (i.e. soil cutting, contact burst and water-like material flowing) show the rationality of PCMM.

Originality/value

In traditional particle DEM, contact relationships are used to calculate contact forces. But in PCMM, contact relationships are adopted to generate elements. Compared to other meshfree methods, in PCMM, the element automatic deletion and recreation technique is used to solve large deformation problems.

Details

Engineering Computations, vol. 32 no. 4
Type: Research Article
ISSN: 0264-4401

Keywords

To view the access options for this content please click here
Article
Publication date: 13 May 2019

Rituraj Singh and Krishna Mohan Singh

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…

Abstract

Purpose

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.

Design/methodology/approach

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.

Findings

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.

Originality/value

The MLPG method based on the stabilised MLS scheme.

Details

Engineering Computations, vol. 36 no. 4
Type: Research Article
ISSN: 0264-4401

Keywords

To view the access options for this content please click here
Article
Publication date: 29 May 2019

Prashant Dineshbhai Vyas, Harish C. Thakur and Veera P. Darji

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

Abstract

Purpose

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

Design/methodology/approach

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

Findings

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.

Research limitations/implications

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.

Originality/value

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

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 30 no. 6
Type: Research Article
ISSN: 0961-5539

Keywords

1 – 10 of 114