Search results

1 – 10 of over 1000
Article
Publication date: 15 May 2009

Evgeny Shavelzon and Dan Givoli

The interaction of a global model (GM) and a local (regional) model (LM) of heat flow is considered under the framework of so‐called “one‐way nesting”. In this framework, the GM…

Abstract

Purpose

The interaction of a global model (GM) and a local (regional) model (LM) of heat flow is considered under the framework of so‐called “one‐way nesting”. In this framework, the GM is constructed in a large domain with coarse discretization in space and time, while the LM is set in a small subdomain with fine discretization.

Design/methodology/approach

The GM is solved first, and its results are then used via some boundary transfer operator (BTO) on the GM–LM interface in order to solve the LM. Past experience in various fields of application has shown that one has to be careful in the choice of BTO to be used on the GM–LM interface, since this choice affects both the stability and accuracy of the computational scheme. Here the problem is first theoretically analyzed for the linear heat equation, and stable BTOs are identified. Then numerical experiments are performed with one‐way nesting in a two‐dimensional channel for heat flow with and without radiation emission and linear reaction, using four different BTOs.

Findings

Among other conclusions, it is shown that the “negative Robin” BTO is unstable, whereas the Dirichlet, Neumann and “positive Robin” BTO are all stable. It is also shown that in terms of accuracy, the Neumann and “positive Robin” BTOs should be preferred over the Dirichlet BTO.

Originality/value

This study may be the first step in analyzing BTO accuracy and stability for more general atmospheric systems.

Details

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

Keywords

Article
Publication date: 1 June 2000

A. Savini

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community…

1133

Abstract

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

Details

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

Keywords

Article
Publication date: 12 August 2020

Manki Cho

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…

Abstract

Purpose

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.

Design/methodology/approach

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.

Findings

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.

Originality/value

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.

Details

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

Keywords

Article
Publication date: 10 April 2009

C.C. Tsai

The purpose of this paper is to develop a meshless numerical method for three‐dimensional isotropic thermoelastic problems with arbitrary body forces.

Abstract

Purpose

The purpose of this paper is to develop a meshless numerical method for three‐dimensional isotropic thermoelastic problems with arbitrary body forces.

Design/methodology/approach

This paper combines the method of fundamental solutions (MFS) and the dual reciprocity method (DRM) as a meshless numerical method (MFS‐DRM) to solve three‐dimensional isotropic thermoelastic problems with arbitrary body forces. In the DRM, the arbitrarily distributed temperature and body force are approximated by polyharmonic splines with augmented polynomial basis, whose particular solutions and the corresponding tractions are reviewed and given explicitly. The MFS is then applied to solve the complementary solution. Numerical experiments of Dirchlet, Robin, and peanut‐shaped‐domain problems are carried out to validate the method.

Findings

In literature, it is commented that the Gaussian elimination can be used reliably to solve the MFS equations for non‐noisy boundary conditions. For noisy boundary conditions, the truncated singular value decomposition (TSVD) is more accurate than the Gaussian elimination. In this paper, it was found that the particular solutions obtained by the DRM act like noises and the use of TSVD improves the accuracy.

Originality/value

It is the first time that the MFS‐DRM is derived to solve three‐dimensional isotropic thermoelastic problems with arbitrary body forces.

Details

Engineering Computations, vol. 26 no. 3
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 1 July 2006

Daniel Ioan, Gabriela Ciuprina and Marius Radulescu

The paper has the purpose of proposing a new open boundary condition to be used in conjunction with the finite integration technique (FIT) for the modelling of passive on‐chip…

Abstract

Purpose

The paper has the purpose of proposing a new open boundary condition to be used in conjunction with the finite integration technique (FIT) for the modelling of passive on‐chip components.

Design/methodology/approach

This boundary condition is ensured by using a virtual layer that surrounds the computational domain.

Findings

The paper proves which are the optimal material properties of the equivalent layer of open boundary.

Practical implications

When modelling passive on‐chip components with FIT, the method proposed is more efficient than the strategic dual image technique.

Originality/value

The paper shows the advantage of this approach – that the analysis algorithm remains unchanged, while saving the field‐circuit compatibility properties, such as current conservation.

Details

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

Keywords

Article
Publication date: 1 February 1998

S. Alfonzetti, G. Borzì and N. Salerno

The Robin iteration procedure is a technique for the FEM computation of electromagnetic scattering fields in unbounded domains. It is based on the iterative improvement of the…

173

Abstract

The Robin iteration procedure is a technique for the FEM computation of electromagnetic scattering fields in unbounded domains. It is based on the iterative improvement of the known term of a non‐homogeneous Robin condition on a fictitious boundary enclosing the scatterer. In this paper it is shown that the procedure is equivalent to the application of the Richardson method to a reduced system and that the use of GMRES significantly reduces the computational effort.

Details

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

Keywords

Article
Publication date: 9 September 2013

Istvan Bardi, Kezhong Zhao, Rickard Petersson, John Silvestro and Nancy Lambert

– This paper aims to present a domain decomposition method to overcome the challenges posed by multi-domain, multi-scale high frequency problems.

Abstract

Purpose

This paper aims to present a domain decomposition method to overcome the challenges posed by multi-domain, multi-scale high frequency problems.

Design/methodology/approach

A hybrid finite element and boundary integral procedure is also presented that allows for domains to employ different solution methods in different subdomains.

Findings

By decomposing large electromagnetic regions into smaller domains, the finite element method can cope with the simulation of electrically large problems.

Practical implications

Real life examples demonstrate the accuracy and efficiency of the new method.

Originality/value

The Robin transmission condition (RTC) is applied to link the domains and preserve field continuity on interfaces.

Details

COMPEL: The International Journal for Computation and Mathematics in Electrical and Electronic Engineering, vol. 32 no. 5
Type: Research Article
ISSN: 0332-1649

Keywords

Article
Publication date: 5 May 2015

Martin Joseph Guillot and Steve C McCool

The purpose of this paper is to investigate the effect of numerical boundary condition implementation on local error and convergence in L2-norm of a finite volume discretization…

Abstract

Purpose

The purpose of this paper is to investigate the effect of numerical boundary condition implementation on local error and convergence in L2-norm of a finite volume discretization of the transient heat conduction equation subject to several boundary conditions, and for cases with volumetric heat generation, using both fully implicit and Crank-Nicolson time discretizations. The goal is to determine which combination of numerical boundary condition implementation and time discretization produces the most accurate solutions with the least computational effort.

Design/methodology/approach

The paper studies several benchmark cases including constant temperature, convective heating, constant heat flux, time-varying heat flux, and volumetric heating, and compares the convergence rates and local to analytical or semi-analytical solutions.

Findings

The Crank-Nicolson method coupled with second-order expression for the boundary derivatives produces the most accurate solutions on the coarsest meshes with the least computation times. The Crank-Nicolson method allows up to 16X larger time step for similar accuracy, with nearly negligible additional computational effort compared with the implicit method.

Practical implications

The findings can be used by researchers writing similar codes for quantitative guidance concerning the effect of various numerical boundary condition approximations for a large class of boundary condition types for two common time discretization methods.

Originality/value

The paper provides a comprehensive study of accuracy and convergence of the finite volume discretization for a wide range of benchmark cases and common time discretization methods.

Details

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

Keywords

Article
Publication date: 3 April 2018

Omar Abu Arqub

The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary

417

Abstract

Purpose

The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary conditions with parameters derivative arising in fluid flows, fluid dynamics, groundwater hydrology, conservation of energy, heat conduction and electric circuit.

Design/methodology/approach

The method provides appropriate representation of the solutions in convergent series formula with accurately computable components. This representation is given in the W(Ω) and H(Ω) inner product spaces, while the computation of the required grid points relies on the R(y,s) (x, t) and r(y,s) (x, t) reproducing kernel functions.

Findings

Numerical simulation with different order derivatives degree is done including linear and nonlinear terms that are acquired by interrupting the n-term of the exact solutions. Computational results showed that the proposed algorithm is competitive in terms of the quality of the solutions found and is very valid for solving such time-fractional models.

Research limitations/implications

Future work includes the application of the reproducing kernel algorithm to highly nonlinear time-fractional partial differential equations such as those arising in single and multiphase flows. The results will be published in forthcoming papers.

Practical implications

The study included a description of fundamental reproducing kernel algorithm and the concepts of convergence, and error behavior for the reproducing kernel algorithm solvers. Results obtained by the proposed algorithm are found to outperform in terms of accuracy, generality and applicability.

Social implications

Developing analytical and numerical methods for the solutions of time-fractional partial differential equations is a very important task owing to their practical interest.

Originality/value

This study, for the first time, presents reproducing kernel algorithm for obtaining the numerical solutions of some certain classes of Robin time-fractional partial differential equations. An efficient construction is provided to obtain the numerical solutions for the equations, along with an existence proof of the exact solutions based upon the reproducing kernel theory.

Article
Publication date: 7 September 2012

Salvatore Alfonzetti, Giuseppe Borzì and Nunzio Salerno

The purpose of this paper is to improve the accuracy of the integral equation of the hybrid FEM‐RBCI (Finite Element Method‐Robin Boundary Condition Iteration) method for the…

Abstract

Purpose

The purpose of this paper is to improve the accuracy of the integral equation of the hybrid FEM‐RBCI (Finite Element Method‐Robin Boundary Condition Iteration) method for the numerical solution of two‐dimensional electromagnetic (or acoustic) scattering problems.

Design/methodology/approach

This accuracy improvement is achieved by selecting the integration curve as straight segments lying in the middle of the triangular finite elements. An accuracy improvement is obtained as compared with selecting the integration curve as constituted by element sides.

Findings

The improved FEM‐RBCI method described in this paper leads to accuracies of the numerical results which are better than those obtained by selecting the integration curve by element sides.

Originality/value

The paper presents results for a simple two‐dimensional structure: a dielectric circular cylinder.

Details

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

Keywords

1 – 10 of over 1000