Hrvoje Gotovac¹, Vedrana Kozulić², Blaž Gotovac¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 173-198, 2010, DOI:10.3970/cmc.2010.015.173

Abstract The space-time Adaptive Fup Collocation Method (AFCM) for solving boundary-initial value problems is presented. To solve the one-dimensional initial boundary value problem, we convert the problem into a two-dimensional boundary value problem. This quasi-boundary value problem is then solved simultaneously in the space-time domain with a collocation technique and by using atomic Fup basis functions. The proposed method is a generally meshless methodology because it requires only the addition of collocation points and basis functions over the domain, instead of the classical domain discretization and numerical integration. The grid is adapted progressively by setting the threshold as a direct measure… More >

Y.B. Zhang¹, A. Godfrey², D. Juul Jensen¹

CMC-Computers, Materials & Continua, Vol.14, No.3, pp. 197-208, 2009, DOI:10.3970/cmc.2009.014.197

Abstract Two methods to quantify protrusions/retrusions and to estimate local boundary curvature from sample plane sections are proposed. The methods are used to evaluate the driving force due to curvature of the protrusions/retrusions for partially recrystallized pure nickel cold rolled to 96% reduction in thickness. The results reveal that the values calculated by both these methods are reasonable when compared with the stored energy measured by differential scanning calorimetry. The relationship between protrusions and the average stored energy density in the deformed matrix is also investigated for partially recrystallized pure aluminum cold rolled to 50%. The results show that the local… More >

B. Tomas Johansson¹, Liviu Marin²

CMC-Computers, Materials & Continua, Vol.13, No.2, pp. 153-190, 2009, DOI:10.3970/cmc.2009.013.153

Abstract We propose two algorithms involving the relaxation of either the given Dirichlet data or the prescribed Neumann data on the over-specified boundary, in the case of the alternating iterative algorithm of Kozlov, Maz'ya and Fomin(1991) applied to Cauchy problems for the modified Helmholtz equation. A convergence proof of these relaxation methods is given, along with a stopping criterion. The numerical results obtained using these procedures, in conjunction with the boundary element method (BEM), show the numerical stability, convergence, consistency and computational efficiency of the proposed methods. More >

D.P.Mondal¹, S.Das¹, Anshul Badkul¹, Nidhi Jha¹

CMC-Computers, Materials & Continua, Vol.13, No.2, pp. 115-134, 2009, DOI:10.3970/cmc.2009.013.115

Abstract The cooling effect by evaporative liquid is modeled by considering that heat is transferred from the system to the surrounding due to evaporation of liquid through the pores present in the medium. The variation of cooling rate with cell size, volume fraction of pores and physical conditions has been analyzed. The model demonstrates that it increases with increase in thickness of the foam slab and with increase in velocity of air. It is also observed that cooling effect decreases with decrease in volume fraction of porosity and with increase in relative density, cell size, thermal conductivity and relative humidity. More >

Shu-Feng Lee¹, Che-Wun Hong^1,2

CMC-Computers, Materials & Continua, Vol.12, No.3, pp. 223-236, 2009, DOI:10.3970/cmc.2009.012.223

Abstract Solid oxides, such as ceria (CeO₂) doped with cations of lower valance, are potential electrolytes for future solid oxide fuel cells. This is due to the theoretically high ionic conductivity at low operation temperature. This paper investigates the feasibility of two potential electrolytes which are samarium-doped ceria (SDC) and gadolinium-doped ceria (GDC) to replace the traditional yttria-stablized zirconia (YSZ). Molecular simulation techniques were employed to study the influence of different dopant concentrations at different operation temperatures on the ionic conductivity from the atomistic perspective. Simulation results show that the optimized ionic conductivity occurs at 11.11mol% concentration using both dopants of… More >

Aimin Jiang^1,2, Yili Wu²

CMC-Computers, Materials & Continua, Vol.12, No.2, pp. 101-120, 2009, DOI:10.3970/cmc.2009.012.101

Abstract This paper presents a development of the boundary contour method (BCM) for piezoelectric media. Firstly, the divergence-free of the integrand of the piezoelectric boundary element method is proved. Secondly, the boundary contour method formulations are obtained by introducing quadratic shape functions and Green's functions (Computer Methods in Applied Mechanics and Engineering1998;158: 65-80) for piezoelectric media and using the rigid body motion solution to regularize the BCM and avoid computation of the corner tensor. The BCM is applied to the problem of piezoelectric media. Finally, numerical solutions for illustrative examples are compared with exact ones. The numerical results of the BCM… More >

LiviuMarin ¹

CMC-Computers, Materials & Continua, Vol.12, No.1, pp. 71-100, 2009, DOI:10.3970/cmc.2009.012.071

Abstract In this paper, the alternating iterative algorithm originally proposed by Kozlov, Maz'ya and Fomin (1991) is numerically implemented for the Cauchy problem in anisotropic heat conduction using a meshless method. Every iteration of the numerical procedure consists of two mixed, well-posed and direct problems which are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise… More >

Amiya K. Samanta¹, Somnath Ghosh²

CMC-Computers, Materials & Continua, Vol.12, No.1, pp. 1-38, 2009, DOI:10.3970/cmc.2009.012.001

Abstract This paper presents a formulation of hypo-elasticity based RC beam model with bond-slip. Details of the constitutive model and analysis method used are provided. A procedure has been described to carry out three-dimensional analysis considering both geometrical as well as material nonlinearity for a simply supported RC beam employing finite element technique, which uses 8-noded isoparametric hexahedral element HCiS18. Enhanced assumed strain (EAS) formulation has been utilized to predict load-deformation and internal stresses both in the elastic as well as nonlinear regime. It models the composite behaviour of concrete and reinforcements in rigid /perfect bond situation and their mutual interaction… More >

D. S. Liu¹, C. Y. Tsai¹, S. R. Lyu²

CMC-Computers, Materials & Continua, Vol.11, No.2, pp. 147-164, 2009, DOI:10.3970/cmc.2009.011.147

Abstract This study presents a novel numerical method for extracting the tempe -rature-dependent mechanical properties of the gold and aluminum thin-films. In the proposed approach, molecular dynamics (MD) simulations are performed to establish the load-displacement response of the thin substrate nanoindented at temperatures ranging from 300-900 K. A simple but effective procedure involving genetic algorithm (GA) and finite element method (FEM) is implemented to extract the material constants of the gold and aluminum substrates. The material constants are then used to construct the corresponding stress-strain curve, from which the elastic modulus, yield stress and the tangent modulus of the thin film… More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.11, No.1, pp. 15-32, 2009, DOI:10.3970/cmc.2009.011.015

Abstract Motivated by the evolutionary and dissipative properties of parabolic type partial differential equation (PDE), Liu (2008a) has proposed a natural and mathematically equivalent approach by transforming the quasilinear elliptic PDE into a parabolic one. However, the above paper only considered a rectangular domain in the plane, and did not treat the difficulty arisen from the quasilinear PDE defined in an arbitrary plane domain. In this paper we propose a new technique of internal and boundary residuals in a fictitious rectangular domain, which are driving forces for the ordinary differential equations based on the Fictitious Time Integration Method (FTIM). Several numerical… More >