Po-Jen Su¹, Cha’o-Kung Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.2, pp. 135-151, 2013, DOI:10.3970/cmes.2013.091.135

Abstract The residual correction method is used to predict the temperature distribution of non-Fourier heat transfer with time-dependent boundary condition. The approximate solution of temperature field is obtained by the residual correction method based on the maximum principle in combination with the finite difference method, making it easier and faster to obtain upper and lower approximations of exact solutions, and even can provide clear definitions of the maximum error bounds of the approximate solutions. The proposed method is found to be an effective numerical method with satisfactory accuracy. More >

Souli, M.¹, Aquelet, N.², Al-Bahkali, E.³, Moatamedi, M.⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.2, pp. 119-133, 2013, DOI:10.3970/cmes.2013.091.119

Abstract To simulate accurately a pressure wave propagation problem, a fine mesh is required in order to capture peak pressures accurately. This may require a very large size problem with several millions of elements. To reduce CPU time and prevent high mesh distortion, a two-dimensional problem for blast ignition and pressure propagation is performed first on a fixed Eulerian mesh. When the pressure wave gets closer to the structure, a three dimensional ALE simulation follows, where the fluid mesh and structure mesh at the fluid structure interface are coincident. The three dimensional problem is performed after mapping history variables from the… More >

J.H. Lv¹, Y. Miao^1,2, H.P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.2, pp. 101-117, 2013, DOI:10.3970/cmes.2013.091.101

Abstract The accurate and efficient evaluation of nearly singular integrals is one of the major concerned problems in the implementation of the boundary element method (BEM). Among the various commonly used nonlinear transformation methods, the distance transformation technique seems to be a promising method to deal with various orders of nearly singular integrals both in potential and elasticity problems. In this paper, some drawbacks of the conventional distance transformation, such as the sensitivity to the position of projection point, are investigated by numerical tests. A general distance transformation technique is developed to circumvent these drawbacks, which is aimed to remove or… More >

Banghua Yao^1,2, Jianping Wei¹, Dengke Wang¹, Dan Ma^2,3, Zhanqing Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.2, pp. 81-100, 2013, DOI:10.3970/cmes.2013.091.081

Abstract Presently, there is an increasing number of water outburst accidents in China as mining activity continues to develop to deeper ground. In these accidents, water outburst caused by karst collapse columns often results in serious damage, involving both the loss of lives and significant economic loss. Therefore, it is of utmost importance to study the seepage property and water outburst mechanism of karst collapse columns. In this paper, based on the seepage theory and the groundwater dynamic theory of porous media, a fluid-solid coupling model for karst collapse columns was built and then imported into COMSOL Multiphysics to be solved,… More >

Xinhua Liu^1,2, Yongzhi Liu¹, Hao Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.1, pp. 65-80, 2013, DOI:10.3970/cmes.2013.091.065

Abstract In order to study the rheological characteristics of magnetorheological fluids, a simulation approach through integration of Monte Carlo method and GPU accelerated technology was proposed and the three-dimensional micro-structure of magnetic particles in different strength magnetic fields were simulated. The Monte Carlo method to magnetic particles of magnetorheological fluids and its key steps such as particle modeling, magnetic energy equations calculating and system state updating were elaborated. Moreover, GPU accelerated technology was applied to the simulation of magnetorheological fluids to reduce computational time and a flowchart for the proposed approach was designed. Finally, a physics experiment was carried out and… More >

Byung-Hyuk Lee¹, Se-min Jeong², Sung-Chul Hwang², Jong-Chun Park³, Moo-Hyun Kim⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.1, pp. 43-63, 2013, DOI:10.3970/cmes.2013.091.043

Abstract The violent free-surface motions interacting with floating vessels containing inner liquid tanks are investigated by using the newly developed Moving Particle Semi-implicit (MPS) method for 2-dimensional incompressible flow simulation. In the present numerical examples, many efficient and robust algorithms have been developed and applied to improve the overall quality and efficiency in solving various highly nonlinear free-surface problems and evaluating impact pressures compared to the original MPS method proposed by Koshizuka and Oka (1996). For illustration, the improved MPS method is applied to the simulation of nonlinear floating-body motions, violent sloshing motions and corresponding impact loads, and vessel motions with… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.1, pp. 17-42, 2013, DOI:10.3970/cmes.2013.091.017

Abstract In the present paper, we propose a novel two-side equilibration method to properly reduce the condition number of a given non-singular matrix only through a few operations. Then, two different conditioners together with the conjugate gradient method (CGM) are developed, which can overcome the defect of CGM, being not vulnerable to noisy disturbance exerted on an ill-posed linear system. The twoside CGM (TSCGM) and the pre-conditioning CGM (PrCGM) are convergent fast and accurate in solving linear inverse problems and the linear Hilbert problem under a large random noise. More >

D. Bergmair^1,2, S. J. Metz¹, H. C. de Lange², A. A. van Steenhoven²

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.1, pp. 1-15, 2013, DOI:10.3970/cmes.2013.091.001

Abstract A method to predict the permeation of water vapor, present in a laminar flowing humid carrier gas, through a hollow fiber membrane is presented. The method uses simulation particles that move like molecules, according to the kinetic gas theory, but carry the physical properties of an ensemble of molecules which they statistically represent. With this approach an ideal operational window for membrane modules can be found and parameters tested for, can be varied over orders of magnitude. The results show that the right dimensioning is essential for the efficient use of the membrane area. More >

M.R. Hematiyan^1,2, A. Khosravifard¹, Y.C. Shiah³, C.L. Tan⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 55-76, 2012, DOI:10.3970/cmes.2012.087.055

Abstract An inverse technique, based on the boundary element method (BEM) and elastostatic experiments for identification of elastic constants of orthotropic and general anisotropic 2D bodies is presented. Displacement measurements at several points on the boundary of the body, obtained by a few known load cases are used in the inverse analysis to find the unknown elastic constants of the body. Using data from more than one elastostatic experiment results in a more accurate and stable solution for the identification problem. In the inverse analysis, sensitivities of displacements of only boundary points with respect to the elastic constants are needed. Therefore,… More >

Yoshihiro Ochiai¹, Vladimir Sladek², Jan Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 41-54, 2012, DOI:10.3970/cmes.2012.087.041

Abstract The conventional boundary element method (BEM) requires a domain integral in unsteady thermal stress analysis with heat generation or an initial temperature distribution. In this paper it is shown that the three-dimensional unsteady thermal stress problem can be solved effectively using the triple-reciprocity boundary element method without internal cells. In this method, the distributions of heat generation and initial temperature are interpolated using integral equations and time-dependent fundamental solutions are used. A new computer program was developed and applied to solving several problems. More >