Roman Chapko¹, B. Tomas Johansson², Oleh Sobeyko¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 57-76, 2010, DOI:10.3970/cmes.2010.062.057

Abstract We propose an iterative procedure for the inverse problem of determining the displacement vector on the boundary of a bounded planar inclusion given the displacement and stress fields on an infinite (planar) line-segment. At each iteration step mixed boundary value problems in an elastostatic half-plane containing the bounded inclusion are solved. For efficient numerical implementation of the procedure these mixed problems are reduced to integral equations over the bounded inclusion. Well-posedness and numerical solution of these boundary integral equations are presented, and a proof of convergence of the procedure for the inverse problem to the original solution is given. Numerical… More >

R. Avila¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.1, pp. 15-64, 2009, DOI:10.3970/cmes.2009.054.015

Abstract The Meshless Local Petrov Galerkin (MLPG) method is used to solve the non-steady two dimensional Navier-Stokes equations. Transient laminar flow field calculations have been carried out in domains wherein certain surfaces have: (i) a sliding motion, (ii) a harmonic motion, (iii) an undulatory movement, and (iv) a contraction-expansion movement. The weak form of the governing equations has been formulated in a Cartesian coordinate system and taking into account the primitive variables of the flow field. A fully implicit pressure correction approach, which requires at each time step an iterative process to solve in a sequential manner the equations which govern… More >

R. Avila¹, F.J. Solorio¹

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 177-202, 2009, DOI:10.3970/cmes.2009.044.177

Abstract Heat transfer processes involving phase change either, solidification or melting, appear frequently in nature and in industrial applications. In this paper the convective patterns that arise from a 2D shear driven annular flow (without and with melting), are presented. The convective annular flow with radial gravity can be considered as a simplified model of the atmospheric flow in the terrestrial equatorial plane (bounded by the warm surface of the Earth and the cold tropopause). The governing equations have been numerically solved by the Spectral Element Method. The numerical results reported in this paper, for the cases without melting (at two… More >

Sirajul Haq^1,2, Siraj-Ul-Islam³, Marjan Uddin^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 115-136, 2009, DOI:10.3970/cmes.2009.044.115

Abstract A mesh free method for the numerical solution of the nonlinear Schrodinger (NLS) and coupled nonlinear Schrodinger (CNLS) equation is implemented. The presented method uses a set of scattered nodes within the problem domain as well as on the boundaries of the domain along with approximating functions known as radial basis functions (RBFs). The set of scattered nodes do not form a mesh, means that no information of relationship between the nodes is needed. Error norms L₂, L_∞ are used to estimate accuracy of the method. Stability analysis of the method is given to demonstrate its practical applicability. More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.3, pp. 243-262, 2009, DOI:10.3970/cmes.2009.041.243

Abstract The Fictitious Time Integration Method (FTIM) previously developed by Liu and Atluri (2008a) is employed here to solve a system of ill-posed linear algebraic equations, which may result from the discretization of a first-kind linear Fredholm integral equation. We rationalize the mathematical foundation of the FTIM by relating it to the well-known filter theory. For the linear ordinary differential equations which are obtained through the FTIM (and which are equivalently used in FTIM to solve the ill-posed linear algebraic equations), we find that the fictitous time plays the role of a regularization parameter, and its filtering effect is better than… More >

Hossein M. Shodja^1,2,3, Alireza Hashemian^1,4

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

Abstract Gradient reproducing kernel particle method (GRKPM) is a meshless technique which incorporates the first gradients of the function into the reproducing equation of RKPM. Therefore, in two-dimensional space GRKPM introduces three types of shape functions rather than one. The robustness of GRKPM's shape functions is established by reconstruction of a third-order polynomial. To enforce the essential boundary conditions (EBCs), GRKPM's shape functions are modified by transformation technique. By utilizing the modified shape functions, the weak form of the nonlinear evolutionary Buckley-Leverett (BL) equation is discretized in space, rendering a system of nonlinear ordinary differential equations (ODEs). Subsequently, Gear's method is… More >

F. Bierbrauer, S.-P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 1-22, 2007, DOI:10.3970/cmes.2007.019.001

Abstract Recently the use of the one-field formulation in the numerical solution of the Navier-Stokes equations for two-phase incompressible flow has become a very attractive approach in CFD (computational fluid dynamics). While the presence of material discontinuities across fluid interfaces presents some difficulty, it is their combination with a non-solenoidal discontinuous initial velocity field, commonly occurring in the mathematical formulation, that has provided the greatest hindrance in the numerical solution. This paper presents three analytical solutions, the Bounded Creeping Flow, Solenoidal and Conserved Solenoidal Solutions, which are both continuous, incompressible, retain as much of the original mathematical formulation as possible and… More >

N. Mai-Duy, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 85-102, 2003, DOI:10.3970/cmes.2003.004.085

Abstract This paper reports a mesh-free Indirect Radial Basis Function Network method (IRBFN) using Thin Plate Splines (TPSs) for numerical solution of Differential Equations (DEs) in rectangular and curvilinear coordinates. The adjustable parameters required by the method are the number of centres, their positions and possibly the order of the TPS. The first and second order TPSs which are widely applied in numerical schemes for numerical solution of DEs are employed in this study. The advantage of the TPS over the multiquadric basis function is that the former, with a given order, does not contain the adjustable shape parameter (i.e. the… More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 365-372, 2001, DOI:10.3970/cmes.2001.002.365

Abstract The aim of this paper is to present a methodology for solving the plane elasticity problem using the Cell Method. It is shown that with the use of a parabolic interpolation in a vectorial problem, a convergence rate of 3.5 is obtained. Such a convergence rate compares with, or is even better than, the one obtained with FEM with the same interpolation – depending on the integration technique used by the FEM application. The accuracy of the solution is also comparable or better. More >

Igor Patlashenko¹, Dan Givoli²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 61-70, 2000, DOI:10.3970/cmes.2000.001.221

Abstract The method of Absorbing Boundary Conditions (ABCs) is considered for the numerical solution of a class of nonlinear exterior wave scattering problems. Recently, a scheme based on the exact nonlocal Dirichlet-to-Neumann (DtN) ABC has been proposed for such problems. Although this method is very accurate, it is also highly expensive computationally. In this paper, the nonlocal ABC is replaced by a low-order local ABC, which is obtained by localizing the DtN condition in a certain "optimal'' way. The performance of the new local scheme is compared to that of the nonlocal scheme via numerical experiments in two dimensions. More >