Haruhiko Kohno, Takahiko Tanahashi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 155-170, 2001, DOI:10.3970/cmes.2001.002.155

Abstract Three-dimensional (3D) unsteady numerical simulations are carried out by means of the finite element method (FEM) with the generalized simplified marker and cell (GSMAC) method in silicon melt with a non-deformable free surface with Prandtl number Pr = 1.8534 × 10^-2, Marangoni number Ma = 0.0 - 6.2067 × 10², Grashof number Gr = 7.1104 × 10⁶, and the aspect ratio As = 1.0 in the Czochralski (CZ) method. The flow state becomes unstable earlier by increasing the absolute value of the thermal coefficient of surface tension in the range of σ_T =0.0 - 1.5 × 10^-5N/mK. Although… More >

K. Nishimura¹, N. Miyazaki²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 143-154, 2001, DOI:10.3970/cmes.2001.002.143

Abstract In this paper, we present a classical molecular dynamics algorithm and its implementation on Cray C90 and Fujitsu VPP700. The characters of this algorithm consist in a grid based on the block division of the atomic system and a neighbor list based on the use of a short range potential. The computer program is used for large scale simulations on a Cray C90 and a 32-node VPP700, and measurements of computational performance are reported. Then, we examine the interaction between a crack propagating and a tilt grain boundary under uniaxial tension using this computer program.… More >

H. Lin, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 117-142, 2001, DOI:10.3970/cmes.2001.002.117

Abstract The truly Meshless Local Petrov-Galerkin (MLPG) method is extended to solve the incompressible Navier-Stokes equations. The local weak form is modified in a very careful way so as to ovecome the so-called Babus^~ka-Brezzi conditions. In addition, The upwinding scheme as developed in Lin and Atluri (2000a) and Lin and Atluri (2000b) is used to stabilize the convection operator in the streamline direction. Numerical results for benchmark problems show that the MLPG method is very promising to solve the convection dominated fluid mechanics problems. More >

Xiang-YangLi¹, Shang-Hua Teng², Peng-Jun Wan³

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 97-116, 2001, DOI:10.3970/cmes.2001.002.097

Abstract To conduct numerical simulations by finite element methods, we often need to generate a high quality mesh, yet with a smaller number of elements. Moreover, the size of each of the elements in the mesh should be approximately equal to a given size requirement. Li et al. recently proposed a new method, named biting, which combines the strengths of advancing front and sphere packing. It generates high quality meshes with a theoretical guarantee. In this paper, we show that biting squares instead of circles not only generates high quality meshes but also has the following advantages. It More >

M.A. Golberg¹, C.S. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 87-96, 2001, DOI:10.3970/cmes.2001.002.087

Abstract The purpose of this paper is to develop a highly efficient mesh-free method for solving nonlinear diffusion-reaction equations in R^d, d=2, 3. Using various time difference schemes, a given time-dependent problem can be reduced to solving a series of inhomogeneous Helmholtz-type equations. The solution of these problems can then be further reduced to evaluating particular solutions and the solution of related homogeneous equations. Recently, radial basis functions have been successfully implemented to evaluate particular solutions for Possion-type equations. A more general approach has been developed in extending this capability to obtain particular solutions for Helmholtz-type equations More >

Ru-Min Chao, Yen-Ji Chen, F.C. Lin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 73-86, 2001, DOI:10.3970/cmes.2001.002.073

Abstract The two-dimensional elasticity problem of an isotropic material, containing a centered-crack with unknown boundary traction is studied by the inverse procedure. The dual boundary integral equations are used to analyze the problem. While solving the ill-posed inverse problem, both of the conjugate gradient method and the regularization method are used. A scaling factor depending upon the material constant μ is introduced into the sensitivity matrix in order to keep the order of magnitude the same throughout the formulation. The result by using the displacement measurement will be compared with those by stress measurement, and an extensive More >

A. Le van¹, P. Le Grognec¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 63-72, 2001, DOI:10.3970/cmes.2001.002.063

Abstract Necking is a bifurcation phenomenon observed in round bars under tensile loading and has been investigated in numbers of papers. In the present work, it is modeled within the framework of finite rate-independent plasticity. The theory is based on thermodynamic foundations developed for standard materials and results in a total Lagrangian formulation for finite plasticity, where the total strain is decomposed additively according to [Green and Nagdhi 1965)] and the hardening is characterized by a nonlinear isotropic hardening law of the saturation type.
The discretization and consistent linearization of the elastic-plastic equation set using the More >

António Tadeu, Andreia Pereira, Luís Godinho¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 49-62, 2001, DOI:10.3970/cmes.2001.002.049

Abstract The Boundary Element Method (BEM) is used to compute the three-dimensional variation pressure field generated by a point pressure source inside a flat waveguide channel filled with a homogeneous fluid, in the presence of infinite rigid circular pipelines. The problem is solved in the frequency domain, using boundary elements to model the pipeline and an appropriate Green's function to simulate the free surface and the rigid floor of the channel. Because of the 2 ---1/2 ---D geometry of the problem, the separation of variables has been used, and the solution at each frequency is expressed in… More >

A. Mandolini¹, V. Minutolo¹, E. Ruocco¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 39-48, 2001, DOI:10.3970/cmes.2001.002.039

Abstract Many sloping areas in the world are affected by slow movements. If they are occupied by settlements or are crossed by roads, pipelines or other infrastructures, a correct evaluation of future displacements is crucial for land management and sometimes for men safety. It is widely recognized that rainfall is the main triggering factor, producing an intermittent and delayed recharge of the groundwater; as a consequence, the displacement rate is cyclic, following a seasonal trend. In Italy this problem is particularly relevant since many exploited sloping areas are affected by slowly moving landslides that interact with More >

A. El-Zafrany¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 27-38, 2001, DOI:10.3970/cmes.2001.002.027

Abstract In a recent publication, the author has introduced boundary integral equations for thick plate bending problems, for cases with generalized types of loading. Internal bending moments and shear forces, required for stress analysis, were calculated by means of a finite difference procedure, which requires fine boundary element meshes to achieve an acceptable degree of accuracy. In this paper, boundary integral equations for internal bending moments and shear forces are presented for thick Reissner plates in bending. Domain loading terms in those boundary integral equations have also been simplified for a variety of loading types including More >