J.H. Wu¹, L.Y. Zhu², J.T. Tang³

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.3, pp. 225-240, 2012, DOI:10.3970/cmes.2012.086.225

Abstract Magnetic induction hyperthermia is one of hopeful methods for tumor therapy. In this method, several ferromagnetic seeds are needed to be implanted into the tumor. The seeds would produce energy, and cause the nearby tumor to die. Temperature prediction is significant before treatment. In addition, in clinical treatment, the tumor temperature has to be monitored in realtime. However, using as few thermometers as possible is the basic principle. Fortunately, the numerical simulation can contribute to realtime measurement. The seed temperature is modeled based on the Haider's method, which is treated as the thermal boundary in numerical simulation. We employ the… More >

J.E. Lee¹and S.J. Kim²

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.2, pp. 93-112, 2012, DOI:10.3970/cmes.2012.086.093

Abstract As the option trading nowadays has become popular, it is important to simulate efficiently large amounts of option pricings. The purpose of this paper is to show valuations of large amount of options, using network distribute computing resources. We valuated 108 options simultaneously on the self-made cluster computer system which is very inexpensive, compared to the supercomputer or the GPU adopting system. For the numerical valuations of options, we developed the option pricing software to solve the Black-Scholes partial differential equation by the finite element method. This yielded accurate values of options and the Greeks with reasonable computational times. This… More >

C.-D. Tran¹, N. Mai-Duy^1,1, K. Le-Cao¹, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.6, pp. 499-520, 2012, DOI:10.3970/cmes.2012.085.499

Abstract A numerical computation of continuum-microscopic model for visco-elastic flows based on the Integrated Radial Basis Function (IRBF) Control Volume and the Stochastic Simulation Techniques (SST) is reported in this paper. The macroscopic flow equations are closed by a stochastic equation for the extra stress at the microscopic level. The former are discretised by a 1D-IRBF-CV method while the latter is integrated with Euler explicit or Predictor-Corrector schemes. Modelling is very efficient as it is based on Cartesian grid, while the integrated RBF approach enhances both the stability of the procedure and the accuracy of the solution. The proposed method is… More >

Jiaoyan Li¹, Xianqiao Wang², James D. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.5, pp. 463-480, 2012, DOI:10.3970/cmes.2012.085.463

Abstract This paper presents a novel multiple time scale algorithm integrated with the concurrent atomic/atom-based continuum modeling, which involves molecular dynamic (MD) simulation and coarse-grained molecular dynamic (CG-MD) simulation. To capture the key features of the solution region while still considering the computational efficiency, we decompose it into two sub-regions in space and utilize the central difference method with different time steps for different sub-regions to march on in time. Usually, the solution region contains a critical field and a non-critical far field. For the critical field (named atomic region) modeled by MD simulation, a relatively small time step is used… More >

V.A. Buryachenko^1,2, T.L. Jackson^2,3, G. Amadio³

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.5, pp. 379-416, 2012, DOI:10.3970/cmes.2012.085.379

Abstract We consider a composite medium, which consists of a homogeneous matrix containing a statistically homogeneous set of multimodal spherical inclusions. This model is used to represent the morphology of heterogeneous solid propellants (HSP) that are widely used in the rocket industry. The Lubachevsky-Stillinger algorithm is used to generate morphological models of HSP with large polydisperse packs of spherical inclusions. We modify the algorithm by proposing a random shaking procedure that leads to the stabilization of a statistical distribution of the simulated structure that is homogeneous, highly mixed, and protocol independent (in sense that the statistical parameters estimated do not depend… More >

Pengfei Yang¹, Yan Liu¹, Xiong Zhang^1,2, Xu Zhou³, Yuli Zhao³

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 207-238, 2012, DOI:10.3970/cmes.2012.085.207

Abstract The material point method is extended to the simulations of fragmentation driven by detonation. A crack modeling scheme based on contact algorithm with material failure process is developed to study the dynamic crack propagation in plastic media. When considering microscopic damage of material, the plastic behavior is described by Gurson model with randomly-distributed initial void of material points. Gurson model can degenerate to J2 plastic theory while the microscopic void is ignored, in which situation the Weibull random failure scheme will be used. Meanwhile, a background-grid-based searching method is proposed to capture the statistical feature of the fragmentation. The scaling… More >

V.G.Yakhno¹, D. Ozdek²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 157-176, 2012, DOI:10.3970/cmes.2012.085.157

Abstract The present paper describes computation of the time-dependent Green's function for the equations of longitudinal vibration of a multi-step rod with a piecewise constant varying cross-section. This computation is based on generalization of the Fourier series expansion method. The time-dependent Green's function is derived in the form of the Fourier series. The basic functions of this series are eigenfunctions of an ordinary differential equation with boundary and matching conditions. Constructing the eigenvalues and eigenfunctions of this differential equation and then derivation of the Fourier coefficients of the Green's function are main steps of the method. Computational experiments confirm the robustness… More >

K. Ibrahim¹, W.A. El-Askary^1,2, A. Balabel¹, I.M. Sakr¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 79-104, 2012, DOI:10.3970/cmes.2012.085.079

Abstract In the present paper, a numerical code has been developed with different turbulence models aiming at simulating turbulent bubbly flows in vertical circular pipes. The mass and momentum conservation equations are used to describe the motion of both phases (water/air). Because of the averaging process additional models are needed for the inter-phase momentum transfer and turbulence quantities for closure. The continuous phase (water) turbulence is represented using different turbulence models namely: two-equation k-ε, extended k-ε and shear-stress transport (SST) k-ω turbulence models which contains additional term to account for the effect of the dispersed phase (air) on the continuous phase… More >

W. A. El-Askary, A. Abdel-Fattah

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.4, pp. 425-458, 2012, DOI:10.3970/cmes.2012.083.425

Abstract In the present paper, experimental and numerical studies of heat transfer and the frictional head loss of turbulent flow in a duct with a heated upper surface are performed. Four different arrangements are considered (case 1: without baffles, case 2: one perforated baffle on the upper wall and one solid baffle on the lower wall, case 3: one perforated baffle on the upper wall and one perforated baffle on the lower wall and case 4: two perforated baffles on the upper wall). A numerical code developed by the present authors is simultaneously presented including four different turbulence models; namely: the… More >

M. Chatterjee, A.K. Mahendra, A.Sanyal, G. Gouthaman

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.4, pp. 385-402, 2012, DOI:10.3970/cmes.2012.083.385

Abstract In this paper, we describe an Incompressible Navier-Stokes (INS) sol -ver using mesh less least square based discretisation on arbitrary distribution of points. The method uses modified Artificial Compressibility Method (ACM) with least square based discretisation. The Solver operates on an arbitrary distribution of points and uses a novel least squares based method that replaces the normal equations approach. This method generates the non-symmetric cross-product matrix by suitable selection of sub stencils such that the matrix is diagonally dominant and well conditioned. The INS solver has been validated with results available in literature for standard test cases. More >