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… 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 More >

Patrícia C.T. Gonçalves^1,2, João Manuel R.S. Tavares^1,2, R.M. Natal Jorge^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.1, pp. 45-56, 2008, DOI:10.3970/cmes.2008.032.045

Abstract The main goals of the present work are to automatically extract the contour of an object and to simulate its deformation using a physical approach. In this work, to segment an object represented in an image, an initial contour is manually defined for it that will then automatically evolve until it reaches the border of the desired object. In this approach, the contour is modelled by a physical formulation using the finite element method, and its temporal evolution to the desired final contour is driven by internal and external forces. The internal forces are defined… More >

Chein-Shan Liu¹, Chih-Wen Chang², Jiang-Ren Chang^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.1, pp. 1-16, 2008, DOI:10.3970/cmes.2008.032.001

Abstract In this paper, we propose a new method to tackle of two famous boundary layer equations in fluid mechanics, namely, the Falkner-Skan and the Blasius equations. We can employ this method to find unknown initial conditions. The pivotal point is based on the erection of a one-step Lie group element$\mathbf {G}(T)$ and the formation of a generalized mid-point Lie group element$\mathbf {G}(r)$. Then, by imposing$\mathbf {G}(T) = \mathbf {G}(r)$ we can seek the missing initial conditions through a minimum discrepancy from the target in terms of a weighting factor$r \in (0, 1)$. Numerical examples are More >

Y.Y Zhang, L. Chen

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 189-200, 2008, DOI:10.3970/cmes.2008.031.189

Abstract A simplified meshless method for dynamic crack growth is presented. The method uses an extrinsic enrichment based on a local partition of unity concept. The crack is represented by a set of crack segments. The crack segments are required to pass through the entire domain of influence of node. They are introduced when the maximum principal stress exceeds the uniaxial tensile strength. The crack segments are allowed to rotate in order to avoid too stiff system responses. The major advantage of our method is that it does not require algorithms to track the crack path. More >

L. Mai-Cao¹, T. Tran-Cong²

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 157-188, 2008, DOI:10.3970/cmes.2008.031.157

Abstract This paper presents a new meshless numerical approach to solving a special class of moving interface problems known as the passive transport where an ambient flow characterized by its velocity field causes the interfaces to move and deform without any influences back on the flow. In the present approach, the moving interface is captured by the level set method at all time as the zero contour of a smooth function known as the level set function whereas one of the two new meshless schemes, namely the SL-IRBFN based on the semi-Lagrangian method and the Taylor-IRBFN More >

L. Codecasa¹, R. Specogna², F. Trevisan³

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 129-144, 2008, DOI:10.3970/cmes.2008.031.129

Abstract In the paper we introduce a methodology to construct discrete constitutive matrices relating magnetic fluxes with magneto motive forces (reluctance matrix) and electro motive forces with currents (conductance matrix) needed for discretizing eddy current problems over hexahedral primal grids by means of the Finite Integration Technique (FIT) and the Cell Method (CM). We prove that, unlike the mass matrices of Finite Elements, the proposed matrices ensure both the stability and the consistency of the discrete equations introduced in FIT and CM. More >

M. Steffen¹, P.C. Wallstedt², J.E. Guilkey^2,3, R.M. Kirby¹, M. Berzins¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 107-128, 2008, DOI:10.3970/cmes.2008.031.107

Abstract The Material Point Method (MPM) has shown itself to be a powerful tool in the simulation of large deformation problems, especially those involving complex geometries and contact where typical finite element type methods frequently fail. While these large complex problems lead to some impressive simulations and solutions, there has been a lack of basic analysis characterizing the errors present in the method, even on the simplest of problems. The large number of choices one has when implementing the method, such as the choice of basis functions and boundary treatments, further complicates this error analysis.\newline In More >

Francisco. P. M. Oliveira¹, João Manuel R. S. Tavares¹

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

Abstract This paper presents a new assignment algorithm with order restriction. Our optimization algorithm was developed using dynamic programming. It was implemented and tested to determine the best global matching that preserves the order of the points that define two contours to be matched. In the experimental tests done, we used the affinity matrix obtained via the method proposed by Shapiro, based on geometric modeling and modal matching. \newline The proposed algorithm revealed an optimum performance, when compared with classic assignment algorithms: Hungarian Method, Simplex for Flow Problems and LAPm. Indeed, the quality of the matching More >

P.H. Wen¹, M.H. Aliabadi², Y.W. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 133-148, 2008, DOI:10.3970/cmes.2008.030.133

Abstract Based on the variation of potential energy, the element-free Galerkin method (MFGM) has been investigated for structures with crack on the basis of radial base function interpolation. An enriched radial base function is introduced to capture the singularities of stress at the crack tips. The advantages of the finite element method are remained in this method and there is a significant improvement of accuracy, particularly for the crack problems of fracture mechanics. The applications of the element-free Galerkin method with enriched radial base function to two-dimensional fracture mechanics in functionally graded materials have been presented More >