Guoqing Yang¹, Baotong Li^2,3, Yang Wang², Jun Hong²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 251-279, 2014, DOI:10.3970/cmes.2014.103.251

Abstract Mechanical behaviors arising at the contact interface largely depend on its surface topographies, particularly when it comes to rough surfaces. A numerical simulation based on an appropriate characterization of rough surfaces especially in terms of three dimensional can be of great significance when it comes to capturing the deformation patterns of micro-scale contacts. In this paper, a simple and practical scheme is developed to generate 3D rough surfaces and to analyze and evaluate the contact characteristics. Firstly amplitude and spatial statistical characterizations of asperities are introduced to avert from the redundancy of topography data caused by traditional measuring methods. A… More >

K. Kakuda¹, Y. Hayashi¹, J. Toyotani¹

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 229-249, 2014, DOI:10.3970/cmes.2014.103.229

Abstract In this paper, we present the application of the particle-based simulations to complicated fluid flow problem with free surfaces. The particle approach is based on the MPS (Moving Particle Simulation) method using hyperbolic-type weighting function to stabilize the spurious oscillatory solutions for solving the Poisson equation with respect to the pressure fields. The hyperbolic-type weighting function is constructed by differentiating the characteristic function based on neural network framework. The weighting function proposed herein is collaterally applied to the kernel function in the SPH-framework. Numerical results demonstrate the workability and validity of the present MPS approach through the dambreaking flow problem. More >

J.R. Zabadal¹, B. Bodmann¹, V. G. Ribeiro², A. Silveira², S. Silveira²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 215-227, 2014, DOI:10.3970/cmes.2014.103.215

Abstract This work presents a new analytical method to transform exact solutions of linear diffusion equations into exact ones for nonlinear advection-diffusion models. The proposed formulation, based on Bäcklund transformations, is employed to obtain velocity fields for the unsteady two-dimensional Helmholtz equation, starting from analytical solutions of a heat conduction type model. More >

Shueei-Muh Lin¹, Min-Jun Teng²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.3, pp. 189-213, 2014, DOI:10.3970/cmes.2014.103.189

Abstract The control and separation of a scanning device moving along an arbitrary trajectory on an elastic plate is investigated. The system is a moving mass problem and is difficult to analyze directly. A semi-analytical method for the movingmass model is presented here. Without vibration control, the separation of a vehicle from a plate is likely to happen. The mechanism of separation of a vehicle from a plate is studied. Moreover, the effects of several parameters on vibration separation and the critical speed of system are studied. An effective control methodology is proposed for suppressing vibration and separation This model is… More >

Haitao Liao

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.3, pp. 155-187, 2014, DOI:10.3970/cmes.2014.103.155

Abstract The improved constrained optimization harmonic balance method(COHBM) is presented to solve the Duffing oscillator with time delayed term. Within the framework of the proposed method, the analytical gradients of the objective function and nonlinear quality constraints with respect to optimization variables are formulated and the sensitivity information of the Fourier coefficients can also obtained. The general formulas of the geometrically nonlinear and time delayed terms are analytically derived, which makes the calculations of nonlinear differential equations in the frequency domain easily. A stability analysis method based on the analytical formulation of the nonlinear equality constraints is presented for the nonlinear… More >

M.D. Campos¹, E.C. Romão²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.3, pp. 139-154, 2014, DOI:10.3970/cmes.2014.103.139

Abstract The objective of this paper aims to present a numerical solution of high accuracy and low computational cost for the three-dimensional Burgers equations. It is a well-known problem and studied the form for one and two-dimensional, but still little explored numerically for three-dimensional problems. Here, by using the High-Order Finite Difference Method for spatial discretization, the Crank-Nicolson method for time discretization and an efficient linearization technique with low computational cost, two numerical applications are used to validate the proposed formulation. In order to analyze the numerical error of the proposed formulation, an unpublished exact solution was used. More >

J. E. Sebold¹, L. A. Lacerda², J. A. M. Carrer³

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.2, pp. 111-137, 2014, DOI:10.3970/cmes.2014.103.111

Abstract This work is concerned with the development of the so-called Whitney and Nédélec edge finite element method for the solution of the time-harmonic Maxwell equations. Initially, the second order time harmonic Maxwell systems, as well as their variational formulation, are presented. In the sequence, Whitney and Nédélec element spaces, whose functions present continuous tangential components along the interface are built of adjacent elements. Then, numerical experiments validate the performance of Whitney and Nédélec first order elements in a two-dimensional domain. The discrete dispersion relation for the elements shows that the numerical phase velocity can be used as an error estimator.… More >

A. Roychowdhury^1,2, A. Nandy¹, C.S. Jog¹, R. Pratap^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.2, pp. 91-110, 2014, DOI:10.3970/cmes.2014.103.091

Abstract We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on… More >

S. Chakraverty^1,2, Smita Tapaswini¹

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.2, pp. 71-90, 2014, DOI:10.3970/cmes.2014.103.071

Abstract Fractional Fornberg-Whitham equation has a vast application in physics. There exist various investigations for the above problem by considering the variables and parameters as crisp/exact. In practice, we may not have these parameters exactly but those may be known in some uncertain form. In the present paper, these uncertainties are taken as interval/fuzzy and the authors proposed here a new method viz. that of the double parametric form of fuzzy numbers to handle the uncertain fractional Fornberg-Whitham equation. Using the single parametric form of fuzzy numbers, original fuzzy fractional Fornberg-Whitham equation is converted first to an interval based fuzzy differential… More >

V. C. Loukopoulos¹, G. C. Bourantas²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.1, pp. 49-70, 2014, DOI:10.3970/cmes.2014.103.049

Abstract A numerical solution of the linear and nonlinear time-dependent Schrödinger equation is obtained, using the strong form MLPG Collocation method. Schrödinger equation is replaced by a system of coupled partial differential equations in terms of particle density and velocity potential, by separating the real and imaginary parts of a general solution, called a quantum hydrodynamic (QHD) equation, which is formally analogous to the equations of irrotational motion in a classical fluid. The approximation of the field variables is obtained with the Moving Least Squares (MLS) approximation and the implicit Crank-Nicolson scheme is used for time discretization. For the two-dimensional nonlinear… More >