Zai You Yan¹, Qiang Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.75, No.3&4, pp. 205-222, 2011, DOI:10.3970/cmes.2011.075.205

Abstract In a boundary element method, the treatment of all the possible singular integrals is very important for the correctness and accuracy of the solutions. Generally, the directional derivative of a weakly singular integral is computed by an integral in the sense of Cauchy principal value if the directional derivative of the weakly singular integral kernel is strongly singular or in the sense of Hadamard finite part integral if it is hypersingular. In this paper, we try to discover how the strongly singular and hypersingular integrals are generated and propose an idea to avoid the appearance of such kind of strongly… More >

Y.C. Shiah¹, C. L. Tan², R.F. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.2, pp. 167-198, 2010, DOI:10.3970/cmes.2010.069.167

Abstract In this paper, fully explicit, algebraic expressions are derived for the first and second derivatives of the Green's function for the displacements in a three dimensional anisotropic, linear elastic body. These quantities are required in the direct formulation of the boundary element method (BEM) for determining the stresses at internal points in the body. To the authors' knowledge, similar quantities have never previously been presented in the literature because of their mathematical complexity. Although the BEM is a boundary solution numerical technique, solutions for the displacements and stresses at internal points are sometimes required for some engineering applications. To this… More >

Y.C. Shiah¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 95-108, 2011, DOI:10.3970/cmes.2011.078.095

Abstract By differentiating the Green function of Ting and Lee (1997) for 3D general anisotropic elastotatics in a spherical coordinate system as an intermediate step, and then using the chain rule, derivatives of up to the second order of this fundamental solution are obtained in exact, explicit, algebraic forms. No tensors of order higher than two are present in these derivatives, thereby allowing these quantities to be numerically evaluated quite expeditiously. These derivatives are required for the computation of the internal point displacements and stresses via Somigliana's identity in BEM analysis. Some examples are presented to demonstrate their successful implementation to… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 253-276, 2009, DOI:10.3970/cmes.2009.043.253

Abstract Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a₀ + a₁x + ... + a_nxⁿ expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a new polynomial interpolation… More >

R. Criado¹, J.E. Ortiz¹, V. Mantič¹, L.J. Gray^1,2, F. París¹

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 151-164, 2007, DOI:10.3970/cmes.2007.022.151

Abstract A numerical implementation of the Somigliana identity in displacements for the solution of 3D elastic problems in exponentially graded isotropic solids is presented. An expression for the fundamental solution in displacements, U_jl, was deduced by Martin et al. (Proc. R. Soc. Lond. A, 458, pp. 1931--1947, 2002). This expression was recently corrected and implemented in a Galerkin indirect 3D BEM code by Criado et al. (Int. J. Numer. Meth. Engng., 2008). Starting from this expression of U_jl, a new expression for the fundamental solution in tractions T_jl has been deduced in the present work. These quite complex expressions of the… More >

Yiorgos-Sokratis Smyrlis¹, Andreas Karageorghis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 535-550, 2003, DOI:10.3970/cmes.2003.004.535

Abstract In this study, we investigate the application of the Method of Fundamental Solutions for the solution of biharmonic Dirichlet problems on a disk. Modifications of the method for overcoming sources of inaccuracy are suggested. We also propose an efficient algorithm for the solution of the resulting systems which exploits the symmetries of the matrices involved. The techniques described in the paper are applied to standard test problems. More >

R. Djebali^1,2, M. El Ganaoui^2,3, H. Sammouda¹, R. Bennacer⁴

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.3, pp. 261-282, 2009, DOI:10.3970/fdmp.2009.005.261

Abstract The simplified thermal lattice Boltzmann model (STLBM) developed by Peng, Shu and Chew (2003) is used in this work to simulate low-Rayleigh-number natural convection in a heated rectangular cavity on a uniform grid. It is shown how by resorting to the double populations approach both hydrodynamic and thermal fields can be effectively simulated. Furthermore, a general benchmark is carried out to account for the effect of different parameters in relatively wide ranges. Results are compared with previous works available in the literature. More >

I. Ueno¹, T. Konisho², T. Kawase³, T. Watanabe⁴

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.1, pp. 21-26, 2008, DOI:10.3970/fdmp.2008.004.021

Abstract The present authors carried out an experimental study with a special interest upon the dynamics of the fluid in the vicinity of the boundary line of three phases; solid-liquid-gas interface, which is so-called `contact line.' The moving droplet on the solid substrate is accompanied with the movement of the boundary line of three phases; solid-liquid-gas interface, which is so-called macroscopic 'contact line.' Existing studies have indicated there is a thin liquid film known as 'precursor film' ahead the contact line of the droplet. In the present study the precursor film was detected by applying conventional ellipsometer, and its existing length… More >

Davood Kalantari¹

FDMP-Fluid Dynamics & Materials Processing, Vol.9, No.2, pp. 169-182, 2013, DOI:10.3970/fdmp.2013.009.169

Abstract In this study, a critical summary of existing spray/wall interaction models is given in synergy with a review of available experimental data. In particular, special attention is devoted to the limitations, difficulties and complexities of the most used approaches in the literatures. An attempt is also made to indicate the bottlenecks and criticalities which typically arise when investigators try to extend results obtained for isolated droplets to the more complex dynamics produced by spray impacts. More >

Gustav Amberg¹, Junichiro Shiomi²

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.1, pp. 81-96, 2005, DOI:10.3970/fdmp.2005.001.081

Abstract A few issues in materials science are reviewed with regard to the importance of fluid flows. The effect of convection on generic solidification problems is discussed. One relevant class of flows in melts is those driven by surface tension gradients. In welding this thermo- or solutocapillary flow will determine the penetration depth, and will depend very sensitively on the composition of the material, through the dependence of surface tension on temperature, presence of surfactants, etc. In crystal growth the convective motion in the melt may cause instabilities that are often undesired in practical processes. The unsteady flow structure can cause… More >