E. Shaik^∗, K.A. Hoffmann^∗, J-F. Dietiker^∗

Molecular & Cellular Biomechanics, Vol.4, No.1, pp. 41-54, 2007, DOI:10.3970/mcb.2007.004.041

Abstract A potential interaction between the local hemodynamics and the artery wall response has been suggested for vascular graft failure by intimal hyperplasia (IH). Among the various hemodynamic factors, wall shear stress has been implicated as the primary factor responsible for the development of IH. In order to explore the role of hemodynamics in the formation of IH in end-to-side anastomosis, computational fluid dynamics is employed. To validate the numerical simulations, comparisons with existing experimental data are performed for both steady and pulsatile flows. Generally, good agreement is observed with the velocity profiles whereas some discrepancies… More >

K. Han ¹, Y. T. Feng ¹, D. R. J. Owen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.2, pp. 87-100, 2007, DOI:10.3970/cmes.2007.018.087

Abstract Numerical procedures are introduced for simulations of irregular particle transport in turbulent flows using the coupled lattice Boltzmann method (LBM) and the discrete element method (DEM). The fluid field is solved by the extended LBM with the incorporation of the Smagorinsky turbulence approach, while particle interaction is modeled by the DEM. The hydrodynamic interactions between fluid and particles are realised through an immersed boundary condition, which gives rise to a coupled solution strategy to model the fluid-particle system under consideration. Main computational aspects comprise the lattice Boltzmann formulation for the solution of fluid flows; the More >

Joseph G. Lawrence, Arunan Nadarajah¹

FDMP-Fluid Dynamics & Materials Processing, Vol.3, No.3, pp. 247-254, 2007, DOI:10.3970/fdmp.2007.003.247

Abstract A simplified 2D pseudo steady state model was developed for an atmospheric chemical vapor deposition (CVD) process on glass. This is used to study the feasibility of converting a continuous coating process to one with discrete glass plates with a gap between them. A preliminary estimate employing mass transfer correlations suggested that there would be significant concentration variations due to the gap between the plates. More detailed studies were done by solving the model numerically employing a finite difference scheme with a vorticity-stream function formulation, and employing the commercial computational fluid dynamics program FIDAP which… More >

Evgeny V. Votyakov¹, Egbert A. Zienicke¹

FDMP-Fluid Dynamics & Materials Processing, Vol.3, No.2, pp. 97-114, 2007, DOI:10.3970/fdmp.2007.003.097

Abstract We simulated numerically the laminar flow in the geometry and the magnetic field of the experimental channel used in [Andreev, Kolesnikov, and Thess (2006)]. This provides detailed information about the electric potential distribution for the laminar regime (numerical simulation) and in the turbulent regime as well (experiment). As follows from comparison of simulated and experimental results, the flow under the magnet is determined by the interaction parameter N = Ha² / Re representing the ratio between magnetic force, determined by the Hartmann number Ha, and inertial force, determined by the Reynolds number Re. We compared two variants: (i)(Re,N)=(2000,18.6)… More >

Haomiao Zhou^1,2, Youhe Zhou^1,3, Xiaojing Zheng¹

CMC-Computers, Materials & Continua, Vol.6, No.3, pp. 201-212, 2007, DOI:10.3970/cmc.2007.006.201

Abstract This paper presents some simulation results of nonlinear dynamic responses for a laminated composite beam embedded by actuators of the giant magnetostrictive material (Terfenol-D) subjected to external magnetic fields, where the giant magnetostrictive materials utilizing the realignment of magnetic moments in response to applied magnetic fields generate nonlinear strains and forces significantly larger than those generated by other smart materials. To utilize the full potential application of the materials in the function and safety designs, e.g., active control of vibrations, the analysis of dynamic responses is requested in the designs as accurately as possible on… More >

Y.R. Li¹, L. Peng¹, W.Y. Shi¹, N. Imaishi²

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.3, pp. 153-166, 2006, DOI:10.3970/fdmp.2006.002.153

Abstract The convective instabilities in semiconduc-tor or oxide melts, significantly affect the quality oflarge crystals grown from the melts by the Czochralskimethod. This paper reviewsour recent numerical stud-ies of thermal convection in annular pools of low-Pr sili-con melt and moderate-Pr silicone oil. The mechanismsof the convective instability are discussed and the criticalconditions for the onset of three-dimensional flow are de-termined. The results show that the hydrothermal wave,characterized by curved spokes, is dominant in a shallowthin pool. In a thick pool of the low-Pr silicon melt, thereappears a standing wave type of oscillatory longitudinalrolls, which moves in the azimuthal More >

H. Kawai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 111-120, 2006, DOI:10.3970/cmes.2006.011.111

Abstract ADVENTURE AutoGL (pronounced as ‘Ote- ga-lu’) is a graphics and GUI library, dedicated for simulation-based research and development. It is designed for the simulation users to develop their own data viewers and editors. Currently, the library is used among many researchers and simulation users, mainly in universities and national research centers. Its functionalities and supported platforms are explained. AutoGL applications of various kinds of simulation methods are demonstrated also. More >

T. Fujimoto¹ and T. Nishioka¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.2, pp. 91-102, 2006, DOI:10.3970/cmes.2006.011.091

Abstract In the dynamic fracture of metallic material, some cracks propagate with the incidence of plastic deformation, and distinct plastic strain remains near the post-propagation area. In order to elucidate these dynamic nonlinear fracture processes, the moving finite element method is developed for nonlinear crack propagation. The T* integral is used as the parameter to estimate crack tip condition. First, the effect of material viscosity and crack propagation velocity have been discussed based on the numerical results for fracture under pure mode I high speed loading. Under mixed mode loading, numerical simulations for fracture path prediction More >

Lingfei Gao¹, Xiaoping Zheng^1,2, Zhenhan Yao¹

CMC-Computers, Materials & Continua, Vol.3, No.1, pp. 25-36, 2006, DOI:10.3970/cmc.2007.003.025

Abstract A general numerical approach is developed to model the elastic behaviours and failure processes of heterogeneous materials. The heterogeneous material body is assumed composed of a large number of convex polygon lattices with different phases. These phases are locally isotropic and elastic-brittle with the different lattices displaying variable material parameters and a Weibull-type statistical distribution. When the effective strain exceeds a local fracture criterion, the full lattice exhibits failure uniformly, and this is modelled by assuming a very small Young modulus value. An auto-select loading method is employed to model the failure process. The proposed More >

J.M.M. Sousa¹

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 211-220, 2005, DOI:10.3970/cmes.2005.009.211

Abstract Time-dependent numerical simulations of the flow through a spherical bulge in a 90-degree asymmetrical bend have been performed for Reynolds numbers in the range 100-400. The present results have demonstrated that the flow reaches asymptotically steady, symmetrical solutions for Reynolds numbers up to 300, whereas a value of 400 for this parameter leads to unsteadiness. The computed flow behavior at this higher Reynolds number has shown to be characterized by an intermittent transition between small-amplitude, irregular oscillations and large-amplitude bursts occurring at a low frequency. In addition, the unsteady flow was asymmetrical and exhibited swirl More >