In this paper, the reciprocal work theorem for viscous fluid flow is established for Newtonian fluids and, based on this theorem, a set of boundary-domain integral equations is derived from the continuity and momentum equations for two-dimensional viscous flow The complex-variable technique is used to compute velocity gradients in the use of the continuity equation. The primary variables involved in these integral equations are velocity, traction and pressure. Although the numerical implementation is only focused on steady incompressible flow, these equations are applicable to solving steady, unsteady, compressible and incompressible problems. In this method, the pressure can be expressed in terms of velocity and traction such that the final system of equations entering the iteration procedure only involves velocity and traction as unknowns. Two commonly cited numerical examples are presented to validate the derived equations.

In this paper, a new set of boundary-domain integral equations is derived from the continuity and momentum equations for three-dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex-variable technique is used to compute the divergence of velocity for internal points, while the traction-recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations.

This paper presents a new inverse analysis approach for identifying material properties and unknown geometries for multi-region problems using the Boundary Element Method (BEM). In this approach, the material properties and coordinates of an unknown region boundary are taken as the optimization variables, and the sensitivity coefficients are computed by the Complex-Variable-Differentiation Method (CVDM). Due to the use of CVDM, the sensitivity coefficients can be accurately determined in a way that is as simple to use as the Finite Difference Method (FDM) and an inverse analysis for a complex composite structure can be easily performed through a similar procedure to the direct computation. Although basic integral equations are presented for heat conduction problems, the application of the proposed algorithm to other problems, such as elastic problems, is straightforward. Two numerical examples are given to demonstrate the potential of the proposed approach.

In this paper, a set of boundary integrals are derived based on a radial integration technique to accurately evaluate two dimensional (2D) and three dimensional (3D), regular and singular domain integrals. A self-contained Fortran code is listed and described for numerical implementation of these boundary integrals. The main feature of the theory is that only the boundary of the integration domain needs to be discretized into elements. This feature can not only save considerable efforts in discretizing the integration domain into internal cells (as in the conventional method), but also make computational results for singular domain integrals more accurate since the integrals have been regularized. Some examples are provided to verify the correctness of the presented formulations and the included code.

unctionally graded materials (FGMs). A meshless boundary element method (BEM) is developed for this purpose. Young’s modulus of the FGMs is assumed to have an exponential variation, while Poisson’s ratio is taken as constant. Since no simple fundamental solutions are available for general FGMs, fundamental solutions for homogeneous, isotropic and linear elastic solids are used in the present BEM, which contains a domain-integral due to the material non-homogeneity. Normalized displacements are introduced to avoid displacement gradients in the domain-integral. The domain-integral is transformed into a boundary integral along the global boundary by using the radial integration method (RIM). To approximate the normalized displacements arising in the domain-integral, basis functions consisting of radial basis functions and polynomials in terms of global coordinates are applied. Numerical results are presented and discussed to show the accuracy and the efficiency of the present meshless BEM.