Qiang Li¹, Jie Ouyang¹, Guorong Wu¹, Xiaoyang Xu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.82, No.3&4, pp. 215-232, 2011, DOI:10.32604/cmes.2011.082.215

Abstract The governing equations of two-phase flows including gas and polymer melt are presented, which are solved using finite volume and domain extension methods with SIMPLEC technology. The melt filling and primary gas penetration in gas-assisted injection molding (GAIM) process are simulated, where the Cross-viscosity model is employed to describe the melt rheological behavior, and the CLSVOF(coupled Level Set and Volume of fluid) method is employed to capture the moving interfaces. In order to test and verify the coupled methods, melt filling in a rectangular plate with an insert is simulated, and the numerical results are More >

M.A. Argenta¹, A.P. Gebert², E.S. Filho³, B.A. Felizari⁴, M.B. Hecke⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.3&4, pp. 159-182, 2011, DOI:10.3970/cmes.2011.079.159

Abstract Various methods in the literature proposesequations to calculate the stiffness as a function of density of bone tissue such as apparent density and ash density among others[Helgason, Perilli, Schileo, Taddei, Brynjolfsson and Viceconti, 2008]. Other ones present a value of an equivalent elasticity modulus, obtained by statistical adjustments of curves generated through mechanical compression tests over various specimens[Chevalier, Pahr, Allmer, Charlebois and Zysset, 2007; Cuppone, Seedhom, Berry and Ostell, 2004]. Bone tissue is a material withdifferent behaviors according to the scale of observation. It has a complex composite hierarchical structure, which is responsible for assign More >

S. Abbasbandy^1,2, V. Sladek³, A. Shirzadi¹, J. Sladek³

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.1, pp. 15-38, 2011, DOI:10.3970/cmes.2011.071.015

Abstract This paper deals with the development of a new method for solution of initial-boundary value problems governed by a couple of nonlinear diffusion equations occurring in the theory of self-organization in non-equilibrium systems. The time dependence is treated by linear interpolation using the finite difference method and the semi-discrete partial differential equations are considered in a weak sense by using the local integral equation method with approximating 2-d spatial variations of the field variables by the Moving Least Squares. The evaluation techniques are discussed and the applicability of the presented method is demonstrated on two More >

T. Sanada¹, K. Ando², T. Colonius²

FDMP-Fluid Dynamics & Materials Processing, Vol.7, No.4, pp. 329-340, 2011, DOI:10.3970/fdmp.2011.007.329

Abstract When a droplet impacts a solid surface at high speed, the contact periphery expands very quickly and liquid compressibility plays an important role in the initial dynamics and the formation of lateral jets. The high speed impact results in high pressures that can account for the surface erosion. In this study, we numerically investigated a high speed droplet impacts on a solid wall. The multicomponent Euler equations with the stiffened equation of state are computed using a FV-WENO scheme with an HLLC Riemann solver that accurately captures shocks and interfaces. In order to compare the More >

S. Homma¹, M. Yokotsuka¹, T. Tanaka¹, K. Moriguchi¹, J. Koga¹, G. Tryggvason²

FDMP-Fluid Dynamics & Materials Processing, Vol.7, No.3, pp. 231-240, 2011, DOI:10.3970/fdmp.2011.007.231

Abstract We develop a three-fluid front-tracking method in order to simulate the motion of an axisymmetry compound droplet, which consists of three immiscible fluids separated by two different interfaces. The two interfaces of the compound droplet are represented by two different sets of the front-tracking elements immersed on the Eulerian grid mesh, where the velocities and the pressure are calculated. The density and viscosity profiles with jumps at the interfaces are successfully determined from the location and the connection information of the front-tracking elements. The motion of a compound droplet is simulated on axisymmetric cylindrical coordinates. More >

Ahmed Meskini, Mostafa Najam, Mustapha El Alami

FDMP-Fluid Dynamics & Materials Processing, Vol.7, No.1, pp. 97-110, 2011, DOI:10.3970/fdmp.2011.007.097

Abstract A numerical mixed convection investigation was carried out to study the enhancement of heat transfer in a square cavity with identical heated rectangular blocks adjacent to its upper wall, and submitted to a vertical jet of fresh air from below. The configuration so defined is an inverted "T"-shaped cavity presenting symmetry with respect to a vertical axis passing by the middle of the openings. The governing equations have been solved using the finite difference method. The parameters of this study are: Rayleigh number 10⁴ ≤ Ra ≤ 10⁶, Reynolds number 1 ≤ Re ≤ 1000, the opening More >

A. Kidess¹, N. Armour¹, S. Dost^1,2

FDMP-Fluid Dynamics & Materials Processing, Vol.7, No.1, pp. 29-56, 2011, DOI:10.3970/fdmp.2011.007.029

Abstract When a metallic liquid is subject to strong magnetic body forces, the issues of convergence and numerical stability may arise in numerical simulations. Handling of magnetic body force terms needs care. In this work we have studied two open codes and discussed the related issues. Magnetic force and mass transport terms were added to these codes. Handling the stability issues was discussed. The developed systems were validated by two benchmark cases. Then, the dissolution process of silicon into the germanium melt was selected as an application. The objective was the numerical study of the dissolution More >

Yong Kou¹, Ke Jin¹, Xiaojing Zheng^1,2

CMC-Computers, Materials & Continua, Vol.22, No.2, pp. 129-146, 2011, DOI:10.3970/cmc.2011.022.129

Abstract Magnetic shielding needs to be employed to ensure proper operation of some electronic equipment which are sensitive to external magnetic interference, such as cryogenic valves located inside the ITER feeder cubicles. This paper is concerned with the shielding efficiency of the magnetic shielding enclosures. A 3-D theoretical model for Fe-Ni alloy magnetic shielding enclosures based on finite element method (FEM) is obtained with the nonlinear law of magnetization. The influence of shielding materials, enclosure configurations, single or multi- layer designs, and apertures on the shielding efficiency is investigated. It is shown that the proposed model More >

Cheng Wang¹, Jianguo Ning¹, Tianbao Ma¹

CMC-Computers, Materials & Continua, Vol.22, No.1, pp. 73-96, 2011, DOI:10.3970/cmc.2011.022.073

Abstract In this paper, we report high resolution simulations using a fifth-order weighted essentially non-oscillatory (WENO) scheme with a third-order TVD Runge-Kutta time stepping method to examine the features of the detonation for gas and condensed explosives. A two-stage chemical reaction model and an ignition and growth model are employed to describe the chemical reaction process for gas and condensed explosives. Based on the Steger-Warming vector flux splitting method, a splitting method is employed when the vector flux does not satisfy the homogeneity property for simulating detonation wave propagation for condensed explosives. The sensibility of flame More >

C.S. Lin, T.C. Chen, C.C. Yu ,M.E. Wu, Y.H. Tu¹

The International Conference on Computational & Experimental Engineering and Sciences, Vol.15, No.3, pp. 103-110, 2010, DOI:10.3970/icces.2010.015.103

Abstract Most old apartments in Taiwan lack fire prevention equipment, making fire awareness and escape difficult, as well as timely fire fighting, which leads to increased death tolls from fire incidents. This research utilizes Fire Dynamics Simulator (FDS) software to analyze and simulate the fire accident that occurred in a single old-style five-story apartment on Siu-Lang Road, Chung-Ho City. In this event, many occupant vehicles were parked at the front door of the apartment building or in nearby parking lanes. The fire engine can only drive in after vehicles were cleared from the fire area, a More >