Chih-Wen Chang^1,2, Chein-Shan Liu³

CMC-Computers, Materials & Continua, Vol.34, No.2, pp. 143-175, 2013, DOI:10.3970/cmc.2013.034.143

Abstract The nonlinear Poisson problems in heat diffusion governed by elliptic type partial differential equations are solved by a modified globally optimal iterative algorithm (MGOIA). The MGOIA is a purely iterative method for searching the solution vector x without using the invert of the Jacobian matrix D. Moreover, we reveal the weighting parameter αc in the best descent vector w = αcE + DTE and derive the convergence rate and find a criterion of the parameter γ. When utilizing αc and γ, we can further accelerate the convergence speed several times. Several numerical experiments are carefully discussed and validated the proposed… More >

Cheng-Te Chi¹, Ming-Hsiao Lee², Wen-Hwa Chen^1,3

CMC-Computers, Materials & Continua, Vol.24, No.3, pp. 239-256, 2011, DOI:10.3970/cmc.2011.024.239

Abstract A novel three-dimensional adaptive element-free Galerkin method (EFGM) based on a uniform background grid is proposed to cope with the problems with extremely large deformation. On the basis of this uniform background grid, an interior adaptive strategy through an error estimation within the analysis domain is developed. By this interior adaptive scheme, additional adaptive nodes are inserted in those regions where the solution accuracy needs to be improved. As opposed to the fixed uniform background grid, these inserted nodes can move along with deformation to describe the particular local deformation of the structure. In addition, a triangular surface technique is… More >

Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.24, No.3, pp. 209-238, 2011, DOI:10.3970/cmc.2011.024.209

Abstract In this study, we employ a semi-analytical scheme to resolve the three-dimensional backward heat conduction problem (BHCP) by utilizing a quasi-bound -ary concept. First, the Fourier series expansion method is used to estimate the temperature field u(x, y, z, t) at any time t < T. Second, we ponder a direct regularization by adding an extra term a(x, y, z, 0) to transform a second-kind Fredholm integral equation for u(x, y, z, 0). The termwise separable property of the kernel function allows us to acquire a closed-form regularized solution. In addition, a tactic to determine the regularization parameter is recommended.… More >

Hsin-Fang Chan¹, Chia-Ming Fan^1,2, Weichung Yeih¹

CMC-Computers, Materials & Continua, Vol.24, No.2, pp. 125-142, 2011, DOI:10.3970/cmc.2011.024.125

Abstract The inverse boundary optimization problem, governed by the Helmholtz equation, is analyzed by the Trefftz method (TM) and the exponentially convergent scalar homotopy algorithm (ECSHA). In the inverse boundary optimization problem, the position for part of boundary with given boundary condition is unknown, and the position for the rest of boundary with additionally specified boundary conditions is given. Therefore, it is very difficult to handle the boundary optimization problem by any numerical scheme. In order to stably solve the boundary optimization problem, the TM, one kind of boundary-type meshless methods, is adopted in this study, since it can avoid the… More >

Chein-Shan Liu^1,2, Chung-Lun Kuo³, Dongjie Liu⁴

CMC-Computers, Materials & Continua, Vol.24, No.2, pp. 105-124, 2011, DOI:10.3970/cmc.2011.024.105

Abstract The inverse Cauchy problems for elliptic equations, such as the Laplace equation, the Poisson equation, the Helmholtz equation and the modified Helmholtz equation, defined in annular domains are investigated. The outer boundary of the annulus is imposed by overspecified boundary data, and we seek unknown data on the inner boundary through the numerical solution by a spring-damping regularization method and its Lie-group shooting method (LGSM). Several numerical examples are examined to show that the LGSM can overcome the ill-posed behavior of inverse Cauchy problem against the disturbance from random noise, and the computational cost is very cheap. More >

Gregor Kosec¹, Miha Založnik², Božidar Šarler¹, Hervé Combeau²

CMC-Computers, Materials & Continua, Vol.22, No.2, pp. 169-196, 2011, DOI:10.3970/cmc.2011.022.169

Abstract The simulation of macrosegregation as a consequence of solidification of a binary Al-4.5%Cu alloy in a 2-dimensional rectangular enclosure is tackled in the present paper. Coupled volume-averaged governing equations for mass, energy, momentum and species transfer are considered. The phase properties are resolved from the Lever solidification rule, the mushy zone is modeled by the Darcy law and the liquid phase is assumed to behave like an incompressible Newtonian fluid. Double diffusive effects in the melt are modeled by the thermal and solutal Boussinesq hypothesis. The physical model is solved by the novel Local Radial Basis Function Collocation Method (LRBFCM).… More >

Wenqing Wang¹, Thomas Schnicke², Olaf Kolditz³

CMC-Computers, Materials & Continua, Vol.21, No.3, pp. 217-236, 2011, DOI:10.3970/cmc.2011.021.217

Abstract This work focuses on parallel finite element simulation of thermal hydraulic and mechanical (THM) coupled processes in porous media, which is a common phenomenon in geological applications such as nuclear waste repository and CO2 storage facilities. The Galerkin finite element method is applied to solve the derived partial differential equations. To deal with the coupling terms among the equations, the momentum equation is solved individually in a monolithic manner, and moreover their solving processes are incorporated into the solving processes of nonisothermal hydraulic equation and heat transport equation in a staggered manner. The computation task arising from the present method… More >

Zhao Huiming¹, Dong Zhengzhu¹, Chen Jiarui¹, Yang Min¹

CMC-Computers, Materials & Continua, Vol.21, No.3, pp. 209-216, 2011, DOI:10.3970/cmc.2011.021.209

Abstract The advantages of coupling of a natural boundary element method and a finite element method are introduced. Then we discuss the principle of the direct coupling of NBEM and FEM and its implementation. The comparison of the results between the direct coupling method and FEM proves that the direct coupling method is simple, feasible and valid in practice. More >

Zhuo-Jia Fu^1,2, Wen Chen¹, Qing-Hua Qin^2,3

CMC-Computers, Materials & Continua, Vol.21, No.3, pp. 187-208, 2011, DOI:10.3970/cmc.2011.021.187

Abstract This paper presents a hybrid finite element model (FEM) with a new type of general solution as interior trial functions, named as HGS-FEM. A variational functional corresponding to the proposed general solution is then constructed for deriving the element stiffness matrix of the proposed element model and the corresponding existence of extremum is verified. Then the assumed intra-element potential field is constructed by a linear combination of novel general solutions at the points on the element boundary under consideration. Furthermore, the independent frame field is introduced to guarantee the intra-element continuity. The present scheme inherits the advantages of hybrid Trefftz… More >

Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.21, No.2, pp. 87-106, 2011, DOI:10.3970/cmc.2011.021.087

Abstract We address a new numerical approach to deal with these multi-dimensional backward wave problems (BWPs) in this study. A fictitious time τ is utilized to transform the dependent variable u(x, y, z, t) into a new one by (1+τ)u(x, y, z, t)=: v(x, y, z, t, τ), such that the original wave equation is written as a new hyperbolic type partial differential equation in the space of (x, y, z, t, τ). Besides, a fictitious viscous damping coefficient can be employed to strengthen the stability of numerical integration of the discretized equations by using a group preserving scheme. Several numerical… More >