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 >

M.Ganesh¹, D. Sheen²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 183-194, 2001, DOI:10.3970/cmes.2001.002.183

Abstract A parallel numerical algorithm is introduced and analyzed for solving inhomogeneous initial-boundary value parabolic problems. The scheme is based on the method recently introduced in Sheen, Sloan, and Thomée (2000) for homogeneous problems. We give a method based on a suitable choice of multiple parameters. Our scheme allows one to compute solutions in a wide range of time. Instead of using a standard time-marching method, which is not easily parallelizable, we take the Laplace transform in time of the parabolic problems. The resulting elliptic problems can be solved in parallel. Solutions are then computed by More >

S.T. Lie¹, Guoyou Yu, Z. Zhao²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 171-182, 2001, DOI:10.3970/cmes.2001.002.171

Abstract The BEM/FEM coupling procedure is applied to 2-D time domain structural-acoustic interaction problems. The acoustic domain for fluid or air is modeled by BEM scheme that is suitable for both finite and infinite domains, while the structure is modeled by FEM scheme. The input impact, which can be either plane waves or non-plane waves, can either be forces acting directly on the structural-acoustic system or be explosion sources. The far field or near field explosion sources which are difficult to be simulated by finite element modeling, can be simulated exactly by boundary element modeling as More >

W.W. Liou¹, Y.C. Fang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.4, pp. 119-128, 2000, DOI:10.3970/cmes.2000.001.571

Abstract A simple implicit treatment for the low speed inflow and outflow boundary conditions for the direct simulation Monte Carlo (DSMC) of the flows in microelectromechanical systems (MEMS) is proposed. The local mean flow velocity, temperature, and number density near the subsonic boundaries were used to determine the number of molecules entering the computational domain and their corresponding velocities at every sample average step. The proposed boundary conditions were validated against micro-Poiseuille flows and micro-Couette flows. The results were compared with analytical solutions derived from the Navier-Stokes equations using first-order and second order slip-boundary conditions. The More >

John L. Junkins¹, Kamesh Subbarao², Ajay Verma³

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.4, pp. 99-118, 2000, DOI:10.3970/cmes.2000.001.551

Abstract Model reference adaptive control formulations are presented that rigorously impose the dynamical structure of the state space descriptions of several distinct large classes of dynamical systems. Of particular interest, the formulations enable the imposition of exact kinematic differential equation constraints upon the adaptation process that compensates for model errors and disturbances at the acceleration level. Other adaptive control formulations are tailored for redundantly actuated and constrained dynamical systems. The utility of the resulting structured adaptive control formulations is studied by considering examples from nonlinear oscillations, aircraft control, spacecraft control, and cooperative robotic system control. The More >

Changyu Hwang¹, Philippe H. Geubelle²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.4, pp. 45-56, 2000, DOI:10.3970/cmes.2000.001.497

Abstract This paper presents the formulation and numerical implementation of a spectral scheme specially developed to simulate dynamic fracture events in unidirectional and cross-ply fiber-reinforced composites. The formulation is based on the spectral representation of the transversely isotropic elastodynamic relations between the traction stresses along the fracture plane and the resulting displacements. Example problems involving stationary or dynamically propagating cracks in fiber-reinforced composites are investigated and compared with reference solutions available in the literature and/or experimental observations. More >

José A. González, Ramón Abascal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 141-150, 2000, DOI:10.3970/cmes.2000.001.443

Abstract In this work an approach to the two-dimensional steady-state rolling contact problem, with and without force transmission, is presented. The problem is solved by the combination of the Boundary Element Method with a formulation of the variational inequalities that govern the problem in the contact area, producing finally a mathematical programming problem. This formulation avoids the direct use of the contact constrains, but it drives to the minimisation of a non-differentiable function, being necessary the use of an specific numerical tool as the modified Newton's method. More >

R. Mustata¹, S. D. Harris², L. Elliott¹, D. Lesnic¹, D. B. Ingham¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 107-116, 2000, DOI:10.3970/cmes.2000.001.409

Abstract An inverse boundary element method is developed to characterise the components of the hydraulic conductivity tensor K of anisotropic materials. Surface measurements at exposed boundaries serve as additional input to a Genetic Algorithm (GA) using a modified least squares functional that minimises the difference between observed and BEM-predicted boundary pressure and/or hydraulic flux measurements under current hydraulic conductivity tensor component estimates. 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 More >

H. Lin, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 45-60, 2000, DOI:10.3970/cmes.2000.001.205

Abstract Due to the very general nature of the Meshless Local Petrov-Galerkin (MLPG) method, it is very easy and natural to introduce the upwinding concept (even in multi-dimensional cases) in the MLPG method, in order to deal with convection-dominated flows. In this paper, several upwinding schemes are proposed, and applied to solve steady convection-diffusion problems, in one and two dimensions. Even for very high Peclet number flows, the MLPG method, with upwinding, gives very good results. It shows that the MLPG method is very promising to solve the convection-dominated flow problems, and fluid mechanics problems. More >