Boundary Element Method Analysis for the Bioheat Transfer Equation

1992 ◽  
Vol 114 (3) ◽  
pp. 358-365 ◽  
Author(s):  
Cho Lik Chan
Keyword(s):  
Steady State ◽  
Boundary Element Method ◽  
Boundary Element ◽  
Numerical Solutions ◽  
Conjugate Problem ◽  
Bioheat Transfer ◽  
Bioheat Equation ◽  
Method Analysis ◽  
Good Agreement ◽  
Element Method

In this paper, the boundary element method (BEM) approach is applied to solve the Pennes (1948) bioheat equation. The objective is to develop the BEM formulation and demonstrate its feasibility. The basic BEM formulations for the transient and steady-state cases are first presented. To demonstrate the usefulness of the BEM approach, numerical solutions for 2-D steady-state problems are obtained and compared to analytical solutions. Further, the BEM formulation is applied to model a conjugate problem for an artery imbedded in a perfused heated tissue. Analytical solution is possible when the conduction in the x-direction is negligible. The BEM and analytical results have very good agreement.

Numerical modelling of convection-diffusion problems with first-order chemical reaction using the dual reciprocity boundary element method

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Salam Adel Al-Bayati ◽  
Luiz C. Wrobel
Keyword(s):  
Steady State ◽  
Boundary Element Method ◽  
Boundary Element ◽  
Content Type ◽  
Convection Diffusion ◽  
First Order ◽  
Order Chemical Reaction ◽  
First Order Chemical Reaction ◽  
Diffusion Problems ◽  
Element Method

Purpose The purpose of this paper is to describe an extension of the boundary element method (BEM) and the dual reciprocity boundary element method (DRBEM) formulations developed for one- and two-dimensional steady-state problems, to analyse transient convection–diffusion problems associated with first-order chemical reaction. Design/methodology/approach The mathematical modelling has used a dual reciprocity approximation to transform the domain integrals arising in the transient equation into equivalent boundary integrals. The integral representation formula for the corresponding problem is obtained from the Green’s second identity, using the fundamental solution of the corresponding steady-state equation with constant coefficients. The finite difference method is used to simulate the time evolution procedure for solving the resulting system of equations. Three different radial basis functions have been successfully implemented to increase the accuracy of the solution and improving the rate of convergence. Findings The numerical results obtained demonstrate the excellent agreement with the analytical solutions to establish the validity of the proposed approach and to confirm its efficiency. Originality/value Finally, the proposed BEM and DRBEM numerical solutions have not displayed any artificial diffusion, oscillatory behaviour or damping of the wave front, as appears in other different numerical methods.

An Improved Assembling Algorithm in Boundary Elements With Galerkin Weighting Applied to Three-Dimensional Stokes Flows

2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Sofia Sarraf ◽  
Ezequiel López ◽  
Laura Battaglia ◽  
Gustavo Ríos Rodríguez ◽  
Jorge D'Elía
Keyword(s):  
Boundary Element Method ◽  
Linear System ◽  
Boundary Element ◽  
Numerical Solutions ◽  
Three Dimensional ◽  
Boundary Integral ◽  
Computational Time ◽  
Creeping Flows ◽  
Order Of Magnitude ◽  
Element Method

In the boundary element method (BEM), the Galerkin weighting technique allows to obtain numerical solutions of a boundary integral equation (BIE), giving the Galerkin boundary element method (GBEM). In three-dimensional (3D) spatial domains, the nested double surface integration of GBEM leads to a significantly larger computational time for assembling the linear system than with the standard collocation method. In practice, the computational time is roughly an order of magnitude larger, thus limiting the use of GBEM in 3D engineering problems. The standard approach for reducing the computational time of the linear system assembling is to skip integrations whenever possible. In this work, a modified assembling algorithm for the element matrices in GBEM is proposed for solving integral kernels that depend on the exterior unit normal. This algorithm is based on kernels symmetries at the element level and not on the flow nor in the mesh. It is applied to a BIE that models external creeping flows around 3D closed bodies using second-order kernels, and it is implemented using OpenMP. For these BIEs, the modified algorithm is on average 32% faster than the original one.

Sign in / Sign up

Export Citation Format

Share Document