Yoshihiro OCHIAI¹

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.2, pp. 83-102, 2011, DOI:10.3970/cmes.2011.079.083

Abstract In general, internal cells are required to solve unsteady thermo-elasto -plastic problems using a conventional boundary element method (BEM). However, in this case, the merit of BEM, which is the easy preparation of data, is lost. The conventional multiple-reciprocity boundary element method (MRBEM) cannot be used to solve thermo-elastoplastic problems, because the distribution of initial stress cannot be determined analytically. In this paper, it is shown that two-dimensional unsteady thermo-elastoplastic problems can be solved without the use of internal cells by using the triple-reciprocity BEM and a thin plate spline. The initial stress formulation is More >

D.M.S. Albuquerque¹, J.M.C. Pereira¹, J.C.F. Pereira^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.4, pp. 337-366, 2011, DOI:10.3970/cmes.2011.072.337

Abstract A numerical study of fluid-structure interaction is presented for the analysis of viscous flow over a resonant membrane airfoil in hovering flight. A flexible membrane moving with a prescribed stroke period was naturally excited to enter into 1st, 2nd and 3rd mode of vibration according to the selected membrane tension. The Navier-Stokes equations were discretized on a moving body unstructured grid using the finite volume method. The instantaneous membrane position was predicted by the 1D unsteady membrane equation with input from the acting fluid flow forces. Following initial validation against reported rigid airfoils predictions, the… More >

Anant Diwakar¹, D. N.Srinath¹, Sanjay Mittal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.1, pp. 61-90, 2010, DOI:10.3970/cmes.2010.069.061

Abstract Aerodynamic shape optimization of airfoils is carried out for two values of Reynolds numbers: 10³ and 10⁴, for an angle of attack of 5^o. The objective functions used are (a) maximization of lift (b) minimization of drag and (c) minimization of drag to lift ratio. The surface of the airfoil is parametrized by a 4^th order non-uniform rational B-Spline (NURBS) curve with 61 control points. Unlike the efforts in the past, the relatively large number of control points used in this study offer a rich design shape with the possibility of local bumps and valleys on the… More >

Peter Lucas¹, Alexander H. van Zuijlen¹, Hester Bijl¹

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.1, pp. 79-106, 2010, DOI:10.3970/cmes.2010.059.079

Abstract Despite the advances in computer power and numerical algorithms over the last decades, solutions to unsteady flow problems remain computing time intensive.
In previous work [Lucas, P.,Bijl, H., and Zuijlen, A.H. van(2010)], we have shown that a Jacobian-free Newton-Krylov (JFNK) algorithm, preconditioned with an approximate factorization of the Jacobian which approximately matches the target residual operator, enables a speed up of a factor of 10 compared to nonlinear multigrid (NMG) for two-dimensional, large Reynolds number, unsteady flow computations. Furthermore, in [Lucas, P., Zuijlen, A.H. van, and Bijl, H. (2010)] we show that this algorithm also greatly… More >

Sujoy Mukherjee¹, Ranjan Ganguli^1,2

CMC-Computers, Materials & Continua, Vol.19, No.2, pp. 105-134, 2010, DOI:10.3970/cmc.2010.019.105

Abstract In this study, variational principle is used for dynamic modeling of an Ionic Polymer Metal Composite (IPMC) flapping wing. The IPMC is an Electro-active Polymer (EAP) which is emerging as a useful smart material for `artificial muscle' applications. Dynamic characteristics of IPMC flapping wings having the same size as the actual wings of three different dragonfly species Aeshna Multicolor, Anax Parthenope Julius and Sympetrum Frequens are analyzed using numerical simulations. An unsteady aerodynamic model is used to obtain the aerodynamic forces. A comparative study of the performances of three IPMC flapping wings is conducted. Among More >

Alfredo Nicolás¹, Blanca Bermúdez², Elsa Báez³

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.1, pp. 83-106, 2009, DOI:10.3970/cmes.2009.048.083

Abstract Numerical results of natural convection flows in two-dimensional cavities, filled with air, are presented to study the effects on the characteristics of the flows as some parameters vary: the Rayleigh number Ra and the aspect ratio A of the cavity. This kind of thermal flows may be modeled by the unsteady Boussinesq approximation in stream function-vorticity variables. The results are obtained with a simple numerical scheme, previously reported for isothermal/mixed convection flows, based mainly on a fixed point iterative process applied to the non-linear elliptic system that results after time discretization. The evolution of the flows,… More >

Edoardo Bucchignani¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.1, pp. 37-60, 2009, DOI:10.3970/fdmp.2009.005.037

Abstract This article describes an implicit method for the solution of time dependent Navier-Stokes equations written in terms of vorticity and velocity. The field equations are discretized using a finite volume technique over quadrilateral meshes.
The numerical code has been applied to the classical window cavity test, employing a fine stretched non-uniform grid, in order to provide an accurate steady solution for a high value of the Rayleigh number (10⁸). It has also been performed a simulation for a value of Rayleigh larger than the critical value, in order to show the capabilities of the proposed More >

S.P. Hu¹, D.L. Young^1,2, C.M. Fan¹

The International Conference on Computational & Experimental Engineering and Sciences, Vol.6, No.1, pp. 25-50, 2008, DOI:10.3970/icces.2008.006.025

Abstract In this paper, a novel numerical scheme called (FDMFS), which combines the finite difference method (FDM) and the method of fundamental solutions (MFS), is proposed to simulate the nonhomogeneous diffusion problem with an unsteady forcing function. Most meshless methods are confined to the investigations of nonhomogeneous diffusion equations with steady forcing functions due to the difficulty to find an unsteady particular solution. Therefore, we proposed a FDM with Cartesian grid to handle the unsteady nonhomogeneous term of the equations. The numerical solution in FDMFS is decomposed into a particular solution and a homogeneous solution. The… More >

Marco La Mantia¹, Peter Dabnichki^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 131-154, 2008, DOI:10.3970/cmes.2008.033.131

Abstract A potential flow based boundary element method was devised to obtain the hydrodynamic forces acting on oscillating wings. A new formulation of the unsteady Kutta condition, postulating a finite pressure difference at the trailing edge of the flapping wing and proposed earlier by the authors, is implemented in the numerical procedure. A comparison with published experimental data (Read et al., 2003) is carried out and the three-dimensional computational results showed good agreement, especially if compared with a similar two-dimensional numerical approach (La Mantia and Dabnichki, 2008) and the potential analytical model of Garrick (1936). The More >

S.P. Hu¹, D.L. Young², C.M. Fan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.1, pp. 1-20, 2008, DOI:10.3970/cmes.2008.024.001

Abstract In this paper, a novel numerical scheme called (FDMFS), which combines the finite difference method (FDM) and the method of fundamental solutions (MFS), is proposed to simulate the nonhomogeneous diffusion problem with an unsteady forcing function. Most meshless methods are confined to the investigations of nonhomogeneous diffusion equations with steady forcing functions due to the difficulty to find an unsteady particular solution. Therefore, we proposed a FDM with Cartesian grid to handle the unsteady nonhomogeneous term of the equations. The numerical solution in FDMFS is decomposed into a particular solution and a homogeneous solution. The… More >