Leonid V. Zhigilei¹, Avinash M. Dongare¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.5, pp. 539-556, 2002, DOI:10.3970/cmes.2002.003.539

Abstract Computational modeling has a potential of making an important contribution to the advancement of laser-driven methods in nanotechnology. In this paper we discuss two computational schemes developed for simulation of laser coupling to organic materials and metals and present a multiscale model for laser ablation and cluster deposition of nanostructured materials. In the multiscale model the initial stage of laser ablation is reproduced by the classical molecular dynamics (MD) method. For organic materials, the breathing sphere model is used to simulate the primary laser excitations and the vibrational relaxation of excited molecules. For metals, the two temperature model coupled to… More >

L. Chen¹ J. H. Lee¹, C.-f. Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.3, pp. 199-224, 2012, DOI:10.3970/cmes.2012.086.199

Abstract This paper presents a numerical procedure to model surface tension using the Generalized Interpolation Material Point (GIMP) method which employs a background mesh in solving the equations of motion. The force due to surface tension is formulated at the mesh grid points by using the continuum surface force (CSF) model and then added to the equations of motion at each grid point. In GIMP, we use the grid mass as the color function in CSF and apply a moving average smoothing scheme to the grid mass to improve the accuracy in calculating the surface interface. The algorithm, named as GIMP-CSF,… More >

S. Gerace¹, K. Erhart¹, E. Divo^1,2, A. Kassab¹

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.1, pp. 35-68, 2011, DOI:10.3970/cmes.2011.081.035

Abstract The focus of this work is to demonstrate a novel approach to true CFD automation based on an adaptive Cartesian point distribution process coupled with a Meshless flow solution algorithm. As Meshless method solutions require only an underlying nodal distribution, this approach works well even for complex flow geometries with non-aligned domain boundaries. Through the addition of a so-called shadow layer of body-fitted nodes, application of boundary conditions is simplified considerably, eliminating the stair-casing issues of typical Cartesian-based techniques. This paper describes the approach taken to automatically generate the Meshless nodal distribution, along with the details of an automatic local… More >

H.F.Qiang¹, F.Z.Chen¹, W.R. Gao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 239-262, 2011, DOI:10.3970/cmes.2011.077.239

Abstract Based on smoothed particle hydrodynamics (SPH) method with surface tension proposed by Morris, this paper is intended to modify equations for surface tension by modifying normal and curvature with corrective smoothing particle method (CSPM). Compared with the continuum surface force (CSF) model for surface tension employed in the traditional SPH method, the accuracy in the present paper is much higher in terms of handling the problems with large deformation and surface tension. The reason is that in the traditional SPH method the deficiency of particles is near the boundary and sharp-angled areas, and it causes gross errors of curvature calculation.… More >

A. Gakwaya¹, H. Sharifi², M. Guillot¹, M. Souli³, F. Erchiqui⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.3, pp. 209-266, 2011, DOI:10.3970/cmes.2011.073.209

Abstract A mechanical computer aided design and engineering system can be used to reduce the design-to-manufacture cycle time in metal forming process. Such a system could be built upon a solid modeling geometry engine and an efficient finite element (FE) solver. The maintenance of a high-quality mesh throughout the analysis is an essential feature of an efficient finite element simulation of large strain metal forming problems. In this paper, a mesh adaptation technique employing the Arbitrary Lagrangian-Eulerian formulation (ALE) is applied to some industrial metal forming problems. An ACIS boundary representation of the solid model is employed. This type of representation… More >

Sunilkumar N¹, D Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.3, pp. 249-284, 2010, DOI:10.3970/cmes.2010.066.249

Abstract We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and 1D NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of… More >

M.K. Kang¹, E.H. Kim¹, M.S. Rim¹, I. Lee^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.64, No.3, pp. 227-250, 2010, DOI:10.3970/cmes.2010.064.227

Abstract An engineering model for predicting the behavior of shape memory alloy (SMA) wire is presented in this study. Piecewise linear relations between stress and strain at a given temperature are assumed and the mixture rule of Reuss bounds is applied to get the elastic modulus of the SMAs in the mixed phase. Critical stresses and strains of the start and finish of the phase transformation are calculated at a given temperature by means of a linear constitutive equation and a stress-temperature diagram. Transformation conditions based on the critical stresses are translated in terms of critical strains. Martensite volume fraction and… More >

Y.C. Chen¹, Chyanbin Hwu²

CMES-Computer Modeling in Engineering & Sciences, Vol.57, No.1, pp. 31-50, 2010, DOI:10.3970/cmes.2010.057.031

Abstract The Green's function for anisotropic bimaterials has been investigated around three decades ago. Since the mathematical formulation of piezoelectric elasticity can be organized into the same form as that of anisotropic elasticity by just expanding the dimension of the corresponding matrix to include the piezoelectric effects, the extension of the Green's function to piezoelectric bimaterials can be obtained immediately through the associated anisotropic bimaterials. In this paper, the Green's function for the bimaterials bonded together with one anisotropic material and one piezoelectric material is derived by applying Stroh's complex variable formalism with the aid of analytical continuation method. For this… More >

L.M.J.S. Dinis¹, R.M. Natal Jorge², J. Belinha³

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 1-34, 2009, DOI:10.3970/cmes.2009.044.001

Abstract The Natural Neighbour Radial Point Interpolation Method (NNRPIM) is extended to the large deformation analysis of non-linear elastic structures. The nodal connectivity in the NNRPIM is enforced using the Natural Neighbour concept. After the Voronoï diagram construction of the unstructured nodal mesh, which discretize the problem domain, small cells are created, the "influence-cells". These cells are in fact influence-domains entirely nodal dependent. The Delaunay triangles are used to create a node-depending background mesh used in the numerical integration of the NNRPIM interpolation functions. The NNRPIM interpolation functions, used in the Galerkin weak form, are constructed with the Radial Point Interpolators.… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 253-276, 2009, DOI:10.3970/cmes.2009.043.253

Abstract Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a₀ + a₁x + ... + a_nxⁿ expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a new polynomial interpolation… More >