Liming Xiong¹, Youping Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 203-214, 2008, DOI:10.3970/cmes.2008.024.203

Abstract This paper presents the application of an atomistic field theory (AFT) in modeling and simulation of boron- , boron/nitrogen and silicon/nitrogen-doped nanocrystalline silicon carbide (B-, BN-, SiN-SiC). Intergranular glassy films (IGFs) and nano-sized pores have been obtained in triple junctions of the grains in nanocrystalline SiC (nc-SiC). Residual tensile stress in the SiC grains and compressive stress in the grain boundaries (GBs) are observed. Under uniaxial tension, the constitutive responses of nanocrystalline SiC were reproduced from the simulations. It is found that the mechanical properties of nanocrystalline SiC are strongly dependent on the compositions of More >

Hanqing Jiang¹, Junqiang Lu², Min-Feng Yu², Yonggang Huang³

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 95-102, 2008, DOI:10.3970/cmes.2008.024.095

Abstract The periodic lattice registry of multi-walled carbon nanotubes (MWCNTs) have been exploited for the possibilities of development of nanodevices. This paper studied the telescoping behaviors of double-walled carbon nanotubes (DWCNTs) by atomic-scale finite element and tight-bind Green function methods. It was found that telescoping a DWCNT (e.g., (6,3)/(12,6)) will induce a rotational motion of the inner CNT that has a chirl angle θ (0° < θ < 30°). This telescoping-induced rotational motion does not exist for armchair and zigzag DWCNTs due to the symmetry of CNTs. The rotational angle is completely determined by the chirality of More >

S.Y. Wang^1,2, K.M. Lim^2,3, B.C. Khoo^2,3, M.Y. Wang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.1, pp. 21-50, 2008, DOI:10.3970/cmes.2008.024.021

Abstract In topology optimization, filtering techniques have become quite popular in practice. In this paper, an accurate and efficient hybrid sensitivity filtering approach based on the traditional and bilateral sensitivity filtering approaches is proposed. In the present hybrid approach, the traditional sensitivity filter is applied to a sub-domain where numerical instabilities are likely to occur to overcome the numerical instabilites robustly. Filtering on mesh-independent holes identified by an image-processing-based technique is prohibited to reduce the computational cost. The bilateral approach is employed for the corresponding nearest neighboring elements of the mesh-independent holes to drive the 0-1 More >

H. Nguyen-Van¹, N. Mai-Duy², T. Tran-Cong³

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 209-222, 2008, DOI:10.3970/cmes.2008.023.209

Abstract This paper reports a study of linear elastic analysis of two-dimensional piezoelectric structures using a smoothed four-node piezoelectric element. The element is built by incorporating the strain smoothing method of mesh-free conforming nodal integration into the standard four-node quadrilateral piezoelectric finite element. The approximations of mechanical strains and electric potential fields are normalized using a constant smoothing function. This allows the field gradients to be directly computed from shape functions. No mapping or coordinate transformation is necessary so that the element can be used in arbitrary shapes. Through several examples, the simplicity, efficiency and reliability More >

Nagalinga Rajan, Soumyendu Raha¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 91-116, 2008, DOI:10.3970/cmes.2008.023.091

Abstract The article describes a numerical method for time domain integration of noisy dynamical systems originating from engineering applications. The models are second order stochastic differential equations (SDE). The stochastic process forcing the dynamics is treated mainly as multiplicative noise involving a Wiener Process in the Itô sense. The developed numerical integration method is a drift implicit strong order 2.0 method. The method has user-selectable numerical dissipation properties that can be useful in dealing with both multiplicative noise and stiffness in a computationally efficient way. A generalized analysis of the method including the multiplicative noise is More >

Xiang Jiawei¹, Chen Xuefeng², Yang Lianfa³, He Zhengjia⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 1-12, 2008, DOI:10.3970/cmes.2008.023.001

Abstract A class of flat shell elements is constructed by using the scaling functions of two-dimensional tensor product B-spline wavelet on the interval (BSWI). Unlike the process of direct wavelets adding in the wavelet Galerkin method, the element displacement field represented by the coefficients of wavelets expansions was transformed from wavelet space into physical space via the constructed two-dimensional transformation matrix. Then, the BSWI flat shell element is constructed by the assembly of BSWI plane elastomechanics and Mindlin plate elements. Because of the good character of BSWI scaling functions, the BSWI flat shell element combine the More >

G. Accary¹, S. Meradji², D. Morvan², D. Fougère²

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.4, pp. 263-270, 2008, DOI:10.3970/fdmp.2008.004.263

Abstract In the literature, only few references have dealt with mixed-convection flows in the low Mach number approximation. For this reason, in the present study we propose to extend the standard 3D benchmark for mixed convection in a rectangular channel heated from below (Medale and Nicolas, 2005) to the case of large temperature variations (for which the Boussinesq approximation is no longer valid). The Navier-Stokes equations, obtained under the assumption of a low Mach number flow, are solved using a finite volume method. The results, corresponding to the steady-state case of the benchmark, lead to the More >

R. Henda¹, W. Quesnel², Z. Saghir³

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.4, pp. 237-244, 2008, DOI:10.3970/fdmp.2008.004.237

Abstract The transient thermal behavior of a two-dimensional circulating porous bed is analytically investigated. A one-energy equation model, representing both the gas and solid phases via a unified temperature, is employed to describe the thermal behavior of the circulating bed. The latter is essentially a tube and shell heat exchanger commonly used in technologically important applications. The model equation is transformed into a simpler set of partial differential equations using an analytical procedure. The analytical solution, based on the method of separation of variables and the principle of superposition, is formulated for the calculation of the More >

Munther Issa K,ah¹, Jean-Luc Meunier²

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.4, pp. 231-236, 2008, DOI:10.3970/fdmp.2008.004.231

Abstract Multi-walled carbon nanotubes (MWCNTs) have been produced by chemical vapor deposition using acetylene as a hydrocarbon source at 700 \textdegree C. The Monel substrates coated with MWCNTs can be used later as cathodes in the Physical Vapor Deposition (PVD) systems to produce CNTs embedded in diamond-like carbon (DLC) film. This new method of generating MWCNTs on Monel substrates without using any catalyst is very simple and uses very little power compared to other CNTs generating techniques. Embedded CNTs in a crystalline carbon coat was observed on mirror-like polished Monel substrates at 700 \textdegree C but More >

L. Bennamoun¹, A. Belhamri²

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.4, pp. 221-230, 2008, DOI:10.3970/fdmp.2008.004.221

Abstract This work focuses on tyipical heat and mass transfer phenomena during the processing of products in the context of the packed-bed drying method (products arranged in thick layers into dryers working in forced convection mode). The dryers are modeled as porous media at the macroscopic level. The simulations are carried out using the mass, momentum (written in the framework of the Darcy's law approximation) and energy equations applied for the different components. A diffusion model based on Fick's law is also used to take into account the drying kinetics. This approach allows monitoring of the More >