Chein-Shan Liu¹, Li-Wei Liu², Hong-Ki Hong²

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 1-18, 2007, DOI:10.3970/cmes.2007.017.001

Abstract In this paper we are concerned with the parameters identification of the inverse heat conduction problems governed by linear parabolic partial differential equations (PDEs). It is the first time that one can construct a closed-form estimation method for the inverse thermal problems of estimating the spatial-dependent thermophysical parameters. The key points hinge on an establishment of a one-step group preserving scheme (GPS) for the semi-discretization of PDEs, as well as a closed-form solution of the resulting algebraic equations. The new method, namely the Lie-group estimation method, has four advantages: it does not require any prior More >

Cheng-Hung Huang¹, Hsin-Hsien Wu²

CMC-Computers, Materials & Continua, Vol.5, No.3, pp. 197-212, 2007, DOI:10.3970/cmc.2007.005.197

Abstract The conjugate gradient method (CGM) is applied in an inverse fin design problem in estimating the optimum shapes for the non-Fourier spine and longitudinal fins based on the desired fin efficiency and fin volume at the specified time. One of the advantages in using CGM in the inverse design problem lies in that it can handle problems having a huge number of design parameters easily and converges very fast.
The validity of using CGM in solving the present inverse design problem is justified by performing the numerical experiments. Several test cases involving different design fin efficiency, More >

C.P. Providakis¹

Structural Durability & Health Monitoring, Vol.2, No.4, pp. 249-254, 2006, DOI:10.3970/sdhm.2006.002.249

Abstract This paper explores the concept of strain energy density rate in relation to the crack initiation in fracture analysis problems arising in creeping cracked structural components. The analysis of the components is performed by using the boundary element methodology in association with the employment of singular boundary elements for the modeling of the crack tip region. The deformation of the material is assumed to be described by an elastic power law creep model. The strain energy density rate theory is applied to determine the direction of the crack initiation for a center cracked plate in More >

R. Springhetti¹, G. Novati¹, M. Margonari²

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.1, pp. 67-80, 2006, DOI:10.3970/cmes.2006.013.067

Abstract With reference to potential and elastostatic problems, a BEM-FEM coupling procedure, based on the symmetric Galerkin version of the BEM, is developed; the continuity conditions at the interface of the BE and FE subdomains are enforced in weak form; the global linear system is characterized by a symmetric coefficient matrix. The procedure is numerically tested with reference first to 2D potential problems and successively to 3D elastoplastic problems (with plastic strains confined to the FE subdomain). More >

N. Mai-Duy^1,2, T. Tran-Cong², R.I. Tanner³

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.3, pp. 177-186, 2006, DOI:10.3970/cmes.2006.016.177

Abstract This paper presents a new high-order time-kernel boundary-integral-equation method (BIEM) for numerically solving transient problems governed by the Burgers equation. Instead of using high-order Lagrange polynomials such as quadratic and quartic interpolation functions, the proposed method employs integrated radial-basis-function networks (IRBFNs) to represent the unknown functions in boundary and volume integrals. Numerical implementations of ordinary and double integrals involving time in the presence of IRBFNs are discussed in detail. The proposed method is verified through the solution of diffusion and convection-diffusion problems. A comparison of the present results and those obtained by low-order BIEMs and More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², C.L. Tan³

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 57-68, 2006, DOI:10.3970/cmes.2006.016.057

Abstract The Meshless Local Petrov-Galerkin (MLPG) method for linear transient coupled thermoelastic analysis is presented. Orthotropic material properties are considered here. A Heaviside step function as the test functions is applied in the weak-form to derive local integral equations for solving two-dimensional (2-D) problems. In transient coupled thermoelasticity an inertial term appears in the equations of motion. The second governing equation derived from the energy balance in coupled thermoelasticity has a diffusive character. To eliminate the time-dependence in these equations, the Laplace-transform technique is applied to both of them. Local integral equations are written on small More >

K.H. Chen¹, J.T. Chen², J.H. Kao³

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 27-40, 2006, DOI:10.3970/cmes.2006.016.027

Abstract In this paper, we employ the regularized meshless method (RMM) to search for eigenfrequency of two-dimension acoustics with multiply-connected domain. The solution is represented by using the double layer potentials. The source points can be located on the physical boundary not alike method of fundamental solutions (MFS) after using the proposed technique to regularize the singularity and hypersingularity of the kernel functions. The troublesome singularity in the MFS methods is desingularized and the diagonal terms of influence matrices are determined by employing the subtracting and adding-back technique. Spurious eigenvalues are filtered out by using singular More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 179-196, 2006, DOI:10.3970/cmes.2006.015.179

Abstract This paper studies the numerical computations of the second-order singularly perturbed boundary value problems (SPBVPs). In order to depress the singularity we consider a coordinate transformation from the x-domain to the t-domain. The relation between singularity and stiffness is demonstrated, of which the coordinate transformation parameter λ plays a key role to balance these two tendencies. Then we construct a very effective Lie-group shooting method to search the missing initial condition through a weighting factor r ∈ (0,1) in the t-domain formulation. For stabilizing the new method we also introduce two new systems by a translation of More >

A. Wacher¹, D. Givoli²

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 147-164, 2006, DOI:10.3970/cmes.2006.015.147

Abstract The recently proposed String Gradient Weighted Moving Finite Element (SGWMFE) method is extended to include remeshing and refining. The method simultaneously determines, at each time step, the solution of the governing partial differential equations and an optimal location of the finite element nodes. It has previously been applied to the nonlinear time-dependent two-dimensional shallow water equations, under the demanding conditions of large Coriolis forces, inducing large mesh and field rotation. Such effects are of major importance in geophysical fluid dynamics applications. Two deficiencies of the original SGWMFE method are (1) possible tangling of the mesh… More >

S. Chantasiriwan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 137-146, 2006, DOI:10.3970/cmes.2006.015.137

Abstract The multiquadric collocation method is the collocation method based on radial basis function known as multiquadrics. It has been successfully used to solve several linear and nonlinear problems. Although fluid flow problems are among problems previously solved by this method, there is still an outstanding issue regarding the influence of the free parameter of multiquadrics (or the shape parameter) on the performance of the method. This paper provides additional results of using the multiquadric collocation method to solve the lid-driven cavity flow problem. The method is used to solve the problem in the stream function-vorticity More >