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 >

Selden B. Crary¹, Peter Cousseau², David Armstrong¹, David M. Woodcock³, Eva H. Mok¹, Olivier Dubochet⁴, Philippe Lerch⁴, Philippe Renaud²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.1, pp. 127-140, 2000, DOI:10.3970/cmes.2000.001.127

Abstract We present a new and unique software capability for finding statistical optimal designs of deterministic experiments on continuous cuboidal regions. The objective function for the design optimization is the minimization of the expected integrated mean squared error of prediction of the metamodel that will be found, subsequent to the running of the computer simulations, using the best linear unbiased predictor (BLUP). The assumed response-model function includes an unknown, stochastic term, Z. We prove that this criterion, which we name I_Z-optimality, is equivalent to I-optimality for non-deterministic experiments, in the limit of zero correlations among the Z's for More >

Wenjing Ye¹, Subrata Mukherjee²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.1, pp. 111-120, 2000, DOI:10.3970/cmes.2000.001.111

Abstract Polynomial driving-force comb drives are designed using numerical simulation. The electrode shapes are obtained using the indirect boundary element method. Variable gap comb drives that produce combinations of linear, quadratic, and cubic driving-force profiles are synthesized. This inverse problem is solved by an optimization procedure. Sensitivity analysis is carried out by the direct differentiation approach (DDA) in order to compute design sensitivity coefficients (DSCs) of force profiles with respect to parameters that define the shapes of the fingers of a comb drive. The DSCs are then used to drive iterative optimization procedures. Designs of variable More >