Ma Tianbao¹, Wang Cheng, Ning Jianguo

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 155-178, 2008, DOI:10.3970/cmes.2008.033.155

Abstract A multi-material Eulerian hydrodynamic numerical method and hydrocode that can effectively simulate explosion problems in engineering practice were developed in this study. A modified Youngs' interface reconstruction algorithm was proposed for mixed cells, in which the material's volume fractions of the surrounding cells are not only used to reconstruct the material interface but also adopted to determine the transport order of the material. The algorithm developed herein was validated by the modeling of several tests, such as objects with different shapes moving in translational, rotating and shear flow field in two dimensional Descartes coordinates and axis-symmetric cylindrical coordinates. Results show… More >

Marco La Mantia¹, Peter Dabnichki^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 131-154, 2008, DOI:10.3970/cmes.2008.033.131

Abstract A potential flow based boundary element method was devised to obtain the hydrodynamic forces acting on oscillating wings. A new formulation of the unsteady Kutta condition, postulating a finite pressure difference at the trailing edge of the flapping wing and proposed earlier by the authors, is implemented in the numerical procedure. A comparison with published experimental data (Read et al., 2003) is carried out and the three-dimensional computational results showed good agreement, especially if compared with a similar two-dimensional numerical approach (La Mantia and Dabnichki, 2008) and the potential analytical model of Garrick (1936). The need of considering the differences… More >

X.G. Yuan^1,2, H.W. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 175-184, 2008, DOI:10.3970/cmes.2008.032.175

Abstract In this paper, a dynamic problem that describes void formation and motion in an incompressible hyperelastic solid sphere composed of a transversely isotropic Valanis-Landel material is examined, where the sphere is subjected to a class of periodic step tensile loads on its surface. A motion equation of void is derived. On analyzing the dynamical properties of the motion equation and examining the effect of material anisotropy on void formation and motion in the sphere, we obtain some new and interesting results. Firstly, under a constant surface tensile load, it is proved that a void would form in the sphere as… More >

J. Sladek¹, V. Sladek¹, C.L. Tan², S.N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 161-174, 2008, DOI:10.3970/cmes.2008.032.161

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of steady-state and transient heat conduction problems in a continuously non-homogeneous anisotropic medium. The Laplace transform is used to treat the time dependence of the variables for transient problems. The analyzed domain is covered by small subdomains with a simple geometry. A weak formulation for the set of governing equations is transformed into local integral equations on local subdomains by using a unit test function. Nodal points are randomly distributed in the 3D analyzed domain and each node is surrounded by a spherical subdomain to which… More >

Vikas Namdeo^1,2, C S Manohar^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 123-160, 2008, DOI:10.3970/cmes.2008.032.123

Abstract The problem of identification of nonlinear system parameters from measured time histories of response under known excitations is considered. Solutions to this problem are obtained by using the force state mapping technique with two alternative functional representation schemes. These schemes are based on the application of reproducing kernel particle method (RKPM) and kriging techniques to fit the force state map. The RKPM has the capability to reproduce exactly polynomials of specified order at any point in a given domain. The kriging based methods represent the function under study as a random field and the parameters describing this field are optimally… More >

K. Lee¹, H. Park², S.W. Lee³

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 113-122, 2008, DOI:10.3970/cmes.2008.032.113

Abstract A post-processing scheme is presented to accurately determine transverse shear and normal stresses in composite panels undergoing geometrically nonlinear deformation under dynamic loading conditions. Transverse stresses are assumed through thickness at a point of interest and are recovered via a one-dimensional finite element method. The finite element method is based on the least square functional of the error in the equilibrium equation along the thickness direction and utilizes the in-plane stresses and resultant transverse shear forces per unit length obtained by a shell element analysis. Numerical results demonstrate that, with minimal addition of computational efforts, the present post-processing approach can… More >

N. Sageresan, R. Drathi

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 103-112, 2008, DOI:10.3970/cmes.2008.032.103

Abstract Crack propagation in concrete is computed with a simplified meshless method. The material is elastic of Neo-Hookian type until fracture. Then a discrete cohesive crack method is used. In the cohesive crack method, cohesive segments are introduced at the meshless nodes. No representation of the crack surface is needed. The method is well-suited for concrete since concrete develops many cracks. Mesh independent results are obtained due to the cohesive model that takes into account the correct energy dissipation during crack opening. We show the accuracy of our method by comparison to experimental data. More >

Shishir Sinha¹, Praveen Ch.

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 85-102, 2008, DOI:10.3970/cmes.2008.032.085

Abstract Fluid catalytic cracking (FCC) unit plays most important role in the economy of a modern refinery that it is use for value addition to the refinery products. Because of the importance of FCC unit in refining, considerable effort has been done on the modeling of this unit for better understanding and improved productivity. The process is characterized by complex interactions among feed quality, catalyst properties, unit hardware parameters and process conditions. \newline The traditional and global approach of cracking kinetics is lumping. In the present paper, five lump kinetic scheme is considered, where gas oil crack to give lighter fractions… More >

N.B.Petrovskaya ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 69-84, 2008, DOI:10.3970/cmes.2008.032.069

Abstract Discontinuous weighted least--squares (DWLS) approximation is a modification of a standard weighted least-squares approach that nowadays is intensively exploited in computational aerodynamics. A DWLS method is often employed to approximate a solution function over an unstructured computational grid that results in an irregular local support for the approximation. While the properties of a weighted least-squares reconstruction are well known for regular geometries, the approximation over a non-uniform grid is not a well researched area so far. In our paper we demonstrate the difficulties related to the performance of a DWLS method on distorted grids and outline a new approach based… More >

Delfim Soares Jr. ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 57-68, 2008, DOI:10.3970/cmes.2008.032.057

Abstract In this work, a time-domain finite element -- boundary element iterative coupling technique is presented in order to analyse electromagnetic scattering from two-dimensional inhomogeneous objects. In the iterative coupling algorithm, the domain of the original problem is subdivided into sub-domains and each sub-domain is analysed independently (as an uncoupled model) taking into account successive renewals of variables at common interfaces. In order to improve the effectiveness of the iterative coupling approach, the evaluation of an optimised relaxation parameter is discussed, taking into account the minimisation of a square error functional. The algorithm that arises is efficient, accurate and flexible. At… More >