S. Rajagopal¹, Ranjan Ganguli²

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.1, pp. 1-34, 2011, DOI:10.3970/cmes.2011.081.001

Abstract The preliminary wing design of a low speed, long endurance UAV is formulated as a two step optimization problem. The first step performs a single objective aerodynamic optimization and the second step involves a coupled dual objective aerodynamic and structural optimization. During the first step, airfoil geometry is optimized to get maximum endurance parameter at a 2D level with maximum thickness to chord ratio and maximum camber as design variables. Leading edge curvature, trailing edge radius, zero lift drag coefficient and zero lift moment coefficient are taken as constraints. Once the airfoil geometry is finalized,… More >

Chein-Shan Liu^1,2, Satya N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 275-298, 2011, DOI:10.3970/cmes.2011.080.275

Abstract To solve an ill-conditioned system of linear algebraic equations (LAEs): Bx - b = 0, we define an invariant-manifold in terms of r := Bx - b, and a monotonically increasing function Q(t) of a time-like variable t. Using this, we derive an evolution equation for dx / dt, which is a system of Nonlinear Ordinary Differential Equations (NODEs) for x in terms of t. Using the concept of discrete dynamics evolving on the invariant manifold, we arrive at a purely iterative algorithm for solving x, which we label as an Optimal Iterative Algorithm (OIA) involving an Optimal Descent Vector More >

V. Kozulić¹, B. Gotovac¹

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 251-274, 2011, DOI:10.3970/cmes.2011.080.251

Abstract The numerical model for the elasto-plastic analysis of prismatic bars subjected to torsion is developed. The functions implemented in this model are Fup basis functions which belong to the class of atomic functions. The collocation method is used to form a system of equations in which the differential equation of the problem is satisfied in collocation points of closed domain, while boundary conditions are satisfied exactly at the domain boundary. The propagation of plastic zones in the cross-section is monitored by applying the incremental-iterative procedure until failure. An approximate solution of arbitrary accuracy is attained More >

V.G.Yakhno¹, T.M. Yakhno²

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 233-250, 2011, DOI:10.3970/cmes.2011.080.233

Abstract The time-dependent Maxwell's equations in non-dispersive homogeneous bi-anisotropic materials are considered in the paper. These equations are written as a symmetric hyperbolic system. A new method of the computation of the electric and magnetic fields arising from electric current is suggested in the paper. This method consists of the following. The Maxwell's equations are written in terms of the Fourier transform with respect to the space variables. The Fourier image of the obtained system is a system of ordinary differential equations whose coefficients depend on the 3D Fourier parameter. The formula for the solution of More >

Gregory Baker¹, Huaijian Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 219-232, 2011, DOI:10.3970/cmes.2011.080.219

Abstract Boundary integral methods have proved very useful in the simulation of free surface motion, in part, because only information at the surface is necessary to track its motion. However, the velocity of the surface must be calculated quite accurately, and the error must be reasonably smooth, otherwise the surface buckles as numerical inaccuracies grow, leading to a failure in the simulation. For two-dimensional motion, the surface is just a curve and the boundary integrals are simple poles that may be removed, allowing spectrally accurate numerical integration. For three-dimensional motion, the singularity in the integrand, although More >

F.Z. Wang ^1,2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 201-218, 2011, DOI:10.3970/cmes.2011.080.201

Abstract In this paper, we consider the boundary particle method (BPM) which is excellent in solving inhomogeneous partial differential equations in terms of solution accuracy and simplicity. In order to investigate the applicability of the BPM, we examine the relationship between its solution accuracy and the effective condition number. We show that the effective condition number, which estimates system stability with the right-hand side vector taken into account, is inversely proportional to the root mean square error in the numerical approximation. Moreover, for noisy-boundary cases, we find that the BPM can not yield reasonable results, for More >

Sushil Mohan Ratnaker¹, Sanjay Mittal^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 179-200, 2011, DOI:10.3970/cmes.2011.080.179

Abstract Flow past a circular cylinder looses stability at Re ~ 47 via Hopf bifurcation. The eigenmode responsible for the instability leads to the von Kármán vortex shedding. In this work the linear stability of the flow to other modes, near the critical Re, is investigated. In particular, the study explores the possibility of modes other than the Kármán mode having the largest growth rate for Re < Re_cr. To this extent, global linear stability analysis (LSA) of the steady flow past a circular cylinder is carried out for Re = 45 and 48. In addition to the Kármán… More >

D.-A. An-Vo¹, N. Mai-Duy¹, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.2, pp. 141-178, 2011, DOI:10.3970/cmes.2011.080.141

Abstract In this paper, 2-node integrated radial basis function elements (IRBFEs) [CMES, vol.72, no.4, pp.299-334, 2011] are further developed for the simulation of incompressible viscous flows in two dimensions. Emphasis is placed on (i) the incorporation of C²-continuous 2-node IRBFEs into the subregion and point collocation frameworks for the discretisation of the stream function-vorticity formulation on Cartesian grids; and (ii) the development of high order upwind schemes based on 2-node IRBFEs for the case of convection-dominant flows. High levels of accuracy and efficiency of the present methods are demonstrated by solutions of several benchmark problems defined on More >

M.M. Kaminski

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.2, pp. 113-140, 2011, DOI:10.3970/cmes.2011.080.113

Abstract The main aim of this work is to demonstrate a solution to the transient problems for the statistically homogeneous media with random physical parameters. This is done with the use of the stochastic perturbation technique based on the general order Taylor series expansions and the additionally modified implementation of the Finite Element Method. Now, both the Direct Differentiation Method as well as the Response Function Method are employed to form and solve up to the nth order state equations. Computational implementation of both approaches is illustrated using two examples - by determination of the probabilistic More >

Hui-Fen Guo^1,2, Bin-Gang Xu^1,3, Sheng-Yan Li¹, Chong-Wen Yu²

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.2, pp. 87-112, 2011, DOI:10.3970/cmes.2011.080.087

Abstract Three-dimensional transient simulation is presented for swirling recirculating flow in a nozzle with a slotted-tube (different grooves) and the effect of the groove number is also investigated. The numerical results on the streamline angles are validated by experimental visualization using the surface oil flow technology. In the downstream center of the injectors, the vortex breakdown experiences a transition from bubble- to spiral- breakdown as time is increased. For all cases under study, as the sizes of two recirculation zones near the injector upstream wall and the step retain almost constant, the spiral breakdown shows a More >