Search results

1 – 10 of over 4000
To view the access options for this content please click here
Article
Publication date: 8 March 2011

Jianhua Dai, Helder Pinheiro, Jonathan P. Webb and Igor Tsukerman

The purpose of this paper is to extend the generalized finitedifference calculus of flexible local approximation methods (FLAME) to problems where local analytical…

Abstract

Purpose

The purpose of this paper is to extend the generalized finitedifference calculus of flexible local approximation methods (FLAME) to problems where local analytical solutions are unavailable.

Design/methodology/approach

FLAME uses accurate local approximations of the solution to generate difference schemes with small consistency errors. When local analytical approximations are too complicated, semi‐analytical or numerical ones can be used instead. In the paper, this strategy is applied to electrostatic multi‐particle simulations and to electromagnetic wave propagation and scattering. The FLAME basis is constructed by solving small local finite‐element problems or, alternatively, by a local multipole‐multicenter expansion. As yet another alternative, adaptive FLAME is applied to problems of wave propagation in electromagnetic (photonic) crystals.

Findings

Numerical examples demonstrate the high rate of convergence of new five‐ and nine‐point schemes in 2D and seven‐ and 19‐point schemes in 3D. The accuracy of FLAME is much higher than that of the standard FD scheme. This paves the way for solving problems with a large number of particles on relatively coarse grids. FLAME with numerical bases has particular advantages for the multi‐particle model of a random or quasi‐random medium.

Research limitations/implications

Irregular stencils produced by local refinement may adversely affect the accuracy. This drawback could be rectified by least squares FLAME, where the number of stencil nodes can be much greater than the number of basis functions, making the method more robust and less sensitive to the irregularities of the stencils.

Originality/value

Previous applications of FLAME were limited to purely analytical basis functions. The present paper shows that numerical bases can be successfully used in FLAME when analytical ones are not available.

Details

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

Keywords

To view the access options for this content please click here
Article
Publication date: 30 September 2014

Seyed Mahmoud Hosseini

The purpose of this paper is to propose a hybrid mesh-free method based on generalized finite difference (GFD) and Newmark finite difference methods to study the elastic…

Downloads
103

Abstract

Purpose

The purpose of this paper is to propose a hybrid mesh-free method based on generalized finite difference (GFD) and Newmark finite difference methods to study the elastic wave propagation in functionally graded nanocomposite reinforced by carbon nanotubes (FGNRCN). The presented hybrid mesh-free method is applied for a thick hollow cylinder, which is made of FGNRCN and excited by various mechanical shock loadings.

Design/methodology/approach

The FG nanocomposite cylinder is assumed to be under shock loading. The elastic wave propagation is obtained and studied for various nonlinear grading patterns and distributions of the aligned carbon nanotubes. The distribution of carbon naotubes in FG nanocomposite are considered to vary as nonlinear function of radius, which varies with various nonlinear grading patterns continuously through radial direction. The effective material properties of functionally graded carbon nanotube are estimated using a micro-mechanical model.

Findings

The mechanical shock analysis of FGNRCN thick hollow cylinder is carried out and the dynamic behavior of displacement field and the time history of radial displacement are obtained for various grading patterns. An effective hybrid mesh-free method based on GFD and Newmark finite difference methods is presented to calculate the average velocity of elastic wave propagation in FGNRCN. The average velocity of elastic wave propagation is obtained for various grading patterns and various kinds of volume fraction. The effects of some parameters on average velocity of elastic wave propagation are obtained and studied in detail.

Originality/value

The calculation of elastic radial wave propagation in a FGNRCN thick hollow cylinder is presented using a hybrid mesh-free method. The effects of some parameters on wave propagation such as various grading patterns of distribution of carbon nanotubes are studied in details.

Details

Engineering Computations, vol. 31 no. 7
Type: Research Article
ISSN: 0264-4401

Keywords

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

L. De Biase, F. Feraudi and V. Pennati

A new finite volume (FV) method is proposed for the solution ofconvection‐diffusion equations defined on 2D convex domains of general shape.The domain is approximated by a…

Abstract

A new finite volume (FV) method is proposed for the solution of convection‐diffusion equations defined on 2D convex domains of general shape. The domain is approximated by a polygonal region; a structured non‐uniform mesh is defined; the domain is partitioned in control volumes. The conservative form of the problem is solved by imposing the law to be verified on each control volume. The dependent variable is approximated to the second order by means of a quadratic profile. When, for the hyperbolic equation, discontinuities are present, or when the gradient of the solution is very high, a cubic profile is defined in such a way that it enjoys unidirectional monotonicity. Numerical results are given.

Details

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

Keywords

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

Ewa Majchrzak and Bohdan Mochnacki

The purpose of this paper is the application of the finite difference method (FDM) for numerical modeling of the microscale heat transfer processes occurring in the domain…

Abstract

Purpose

The purpose of this paper is the application of the finite difference method (FDM) for numerical modeling of the microscale heat transfer processes occurring in the domain of thin metal film subjected to a laser pulse. The problem discussed is described by the different variants of the second-order dual-phase-lag equation (DPLE). The laser action is taken into account by the introduction of internal volumetric heat source to the governing equation. The capacity of the source is dependent on the geometrical co-ordinates and duration of the laser beam. The modified forms of DPLE presented in the paper, resulting from the certain substitutions introduced to the basic equation.

Design/methodology/approach

At the stage of numerical computations, the different variants of the FDM are applied. Both the explicit and implicit FDM schemes are used. The formula determining the capacity of the internal heat source suggests the formulation of the task discussed using the cylindrical co-ordinate system. The in-house programs realizing the numerical computations concern the axially-symmetrical tasks. In this paper, the metal films made of the nickel and gold are considered.

Findings

The algorithms presented make possible to analyze the heating/cooling processes occurring in the domain of metal film having a thickness Z for the different laser parameters (laser intensity, characteristic time of laser pulse and laser beam radius) and the different materials (optical penetration depth, reflectivity of irradiated surface, lag times, thermal conductivity and volumetric specific heat).

Research limitations/implications

Not for all metals, one can find information on lag times. In the literature, analytical formulas can be found to calculate these values, but they are strongly approximated. It should be pointed out that there are some limitations concerning the delay times of material considered, which assure the physical correctness of the second-order DPLE.

Originality/value

The FDM algorithm concerns the three-dimensional cylindrical domain while a large majority of the second-order DPLE numerical solutions have been obtained for the one-dimensional tasks. Both the implicit and explicit numerical schemes are proposed and the testing computations confirm the correctness and effectiveness of the algorithms presented.

Details

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

Keywords

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

A. A. Alanazi, Sultan Z. Alamri, S. Shafie and Shazirawati Mohd Puzi

The purpose of this paper is to obtain the nonlinear Schrodinger equation (NLSE) numerical solutions in the presence of the first-order chromatic dispersion using a…

Abstract

Purpose

The purpose of this paper is to obtain the nonlinear Schrodinger equation (NLSE) numerical solutions in the presence of the first-order chromatic dispersion using a second-order, unconditionally stable, implicit finite difference method. In addition, stability and accuracy are proved for the resulting scheme.

Design/methodology/approach

The conserved quantities such as mass, momentum and energy are calculated for the system governed by the NLSE. Moreover, the robustness of the scheme is confirmed by conducting various numerical tests using the Crank-Nicolson method on different cases of solitons to discuss the effects of the factor considered on solitons properties and on conserved quantities.

Findings

The Crank-Nicolson scheme has been derived to solve the NLSE for optical fibers in the presence of the wave packet drift effects. It has been founded that the numerical scheme is second-order in time and space and unconditionally stable by using von-Neumann stability analysis. The effect of the parameters considered in the study is displayed in the case of one, two and three solitons. It was noted that the reliance of NLSE numeric solutions properties on coefficients of wave packets drift, dispersions and Kerr nonlinearity play an important control not only the stable and unstable regime but also the energy, momentum conservation laws. Accordingly, by comparing our numerical results in this study with the previous work, it was recognized that the obtained results are the generalized formularization of these work. Also, it was distinguished that our new data are regarding to the new communications modes that depend on the dispersion, wave packets drift and nonlinearity coefficients.

Originality/value

The present study uses the first-order chromatic. Also, it highlights the relationship between the parameters of dispersion, nonlinearity and optical wave properties. The study further reports the effect of wave packet drift, dispersions and Kerr nonlinearity play an important control not only the stable and unstable regime but also the energy, momentum conservation laws.

Details

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

Keywords

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

M.J. Huntul, Mohammad Tamsir and Abdullah Ahmadini

The paper aims to numerically solve the inverse problem of determining the time-dependent potential coefficient along with the temperature in a higher-order…

Abstract

Purpose

The paper aims to numerically solve the inverse problem of determining the time-dependent potential coefficient along with the temperature in a higher-order Boussinesq-Love equation (BLE) with initial and Neumann boundary conditions supplemented by boundary data, for the first time.

Design/methodology/approach

From the literature, the authors already know that this inverse problem has a unique solution. However, the problem is still ill-posed by being unstable to noise in the input data. For the numerical realization, the authors apply the generalized finite difference method (GFDM) for solving the BLE along with the Tikhonov regularization to find stable and accurate numerical solutions. The regularized nonlinear minimization is performed using the MATLAB subroutine lsqnonlin. The stability analysis of solution of the BLE is proved using the von Neumann method.

Findings

The present numerical results demonstrate that obtained solutions are stable and accurate.

Practical implications

Since noisy data are inverted, the study models real situations in which practical measurements are inherently contaminated with noise.

Originality/value

The knowledge of this physical property coefficient is very important in various areas of human activity such as seismology, mineral exploration, biology, medicine, quality control of industrial products, etc. The originality lies in the insight gained by performing the numerical simulations of inversion to find the potential co-efficient on time in the BLE from noisy measurement.

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

S. BRANDON and J.J. DERBY

A finite element method for the analysis of combined radiative and conductive heat transport in a finite axisymmetric configuration is presented. The appropriate…

Abstract

A finite element method for the analysis of combined radiative and conductive heat transport in a finite axisymmetric configuration is presented. The appropriate integro‐differential governing equations for a grey and non‐scattering medium with grey and diffuse walls are developed and solved for several model problems. We consider axisymmetric, cylindrical geometries with top and bottom boundaries of arbitrary convex shape. The method is accurate for media of any optical thickness and is capable of handling a wide array of axisymmetric geometries and boundary conditions. Several techniques are presented to reduce computational overhead, such as employing a Swartz‐Wendroff approximation and cut‐off criteria for evaluating radiation integrals. The method is successfully tested against several cases from the literature and is applied to some additional example problems to demonstrate its versatility. Solution of a free‐boundary, combined‐mode heat transfer problem representing the solidification of a semitransparent material, the Bridgman growth of an yttrium aluminium garnet (YAG) crystal, demonstrates the utility of this method for analysis of a complex materials processing system. The method is suitable for application to other research areas, such as the study of glass processing and the design of combustion furnace systems.

Details

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

Keywords

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

Suvranu De and Klaus‐Jürgen Bathe

Computational efficiency and reliability are clearly the most important requirements for the success of a meshless numerical technique. While the basic ideas of meshless…

Abstract

Computational efficiency and reliability are clearly the most important requirements for the success of a meshless numerical technique. While the basic ideas of meshless techniques are simple and well understood, an effective meshless method is very difficult to develop. The efficiency depends on the proper choice of the interpolation scheme, numerical integration procedures and techniques of imposing the boundary conditions. These issues in the context of the method of finite spheres are discussed.

Details

Engineering Computations, vol. 18 no. 1/2
Type: Research Article
ISSN: 0264-4401

Keywords

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

Mas Irfan Purbawanto Hidayat, Bambang Ariwahjoedi and Setyamartana Parman

The purpose of this paper is to present a new approach of meshless local B-spline based finite difference (FD) method for solving two dimensional transient heat conduction…

Downloads
230

Abstract

Purpose

The purpose of this paper is to present a new approach of meshless local B-spline based finite difference (FD) method for solving two dimensional transient heat conduction problems.

Design/methodology/approach

In the present method, any governing equations are discretized by B-spline approximation which is implemented in the spirit of FD technique using a local B-spline collocation scheme. The key aspect of the method is that any derivative is stated as neighbouring nodal values based on B-spline interpolants. The set of neighbouring nodes are allowed to be randomly distributed thus enhanced flexibility in the numerical simulation can be obtained. The method requires no mesh connectivity at all for either field variable approximation or integration. Time integration is performed by using the Crank-Nicolson implicit time stepping technique.

Findings

Several heat conduction problems in complex domains which represent for extended surfaces in industrial applications are examined to demonstrate the effectiveness of the present approach. Comparison of the obtained results with solutions from other numerical method available in literature is given. Excellent agreement with reference numerical method has been found.

Research limitations/implications

The method is presented for 2D problems. Nevertheless, it would be also applicable for 3D problems.

Practical implications

A transient two dimensional heat conduction in complex domains which represent for extended surfaces in industrial applications is presented.

Originality/value

The presented new meshless local method is simple and accurate, while it is also suitable for analysis in domains of arbitrary geometries.

Details

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

Keywords

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

Soheil Bazazzadeh, Arman Shojaei, Mirco Zaccariotto and Ugo Galvanetto

The purpose of this paper is to apply the Peridynamic differential operator (PDDO) to incompressible inviscid fluid flow with moving boundaries. Based on the potential…

Abstract

Purpose

The purpose of this paper is to apply the Peridynamic differential operator (PDDO) to incompressible inviscid fluid flow with moving boundaries. Based on the potential flow theory, a Lagrangian formulation is used to cope with non-linear free-surface waves of sloshing water in 2D and 3D rectangular and square tanks.

Design/methodology/approach

In fact, PDDO recasts the local differentiation operator through a nonlocal integration scheme. This makes the method capable of determining the derivatives of a field variable, more precisely than direct differentiation, when jump discontinuities or gradient singularities come into the picture. The issue of gradient singularity can be found in tanks containing vertical/horizontal baffles.

Findings

The application of PDDO helps to obtain the velocity field with a high accuracy at each time step that leads to a suitable geometry updating for the procedure. Domain/boundary nodes are updated by using a second-order finite difference time algorithm. The method is applied to the solution of different examples including tanks with baffles. The accuracy of the method is scrutinized by comparing the numerical results with analytical, numerical and experimental results available in the literature.

Originality/value

Based on the investigations, PDDO can be considered a reliable and suitable approach to cope with sloshing problems in tanks. The paper paves the way to apply the method for a wider range of problems such as compressible fluid flow.

Details

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

Keywords

1 – 10 of over 4000