Mingxu Yi¹, Yiming Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 229-244, 2012, DOI:10.3970/cmes.2012.088.229

Abstract In this paper, Haar wavelet operational matrix method is proposed to solve a class of fractional partial differential equations. We derive the Haar wavelet operational matrix of fractional order integration. Meanwhile, the Haar wavelet operational matrix of fractional order differentiation is obtained. The operational matrix of fractional order differentiation is utilized to reduce the initial equation to a Sylvester equation. Some numerical examples are included to demonstrate the validity and applicability of the approach. More >

Taoufik Naffouti^1,2 and Ridha Djebali^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 211-228, 2012, DOI:10.3970/cmes.2012.088.211

Abstract This paper reports numerical results of natural convection flow evolving inside confined medium defined by two-dimensional square enclosure containing isothermal hot source placed asymmetrically at bottom wall. The sides-walls are isothermally cooled at a constant temperature; however the ceiling and the rest of bottom wall are insulated. The lattice Boltzmann method is used to solve the dimensionless governing equations with the associated boundary conditions. The flow is monitored by the Grashof number and the Prandtl number taken here 0.71. Numerical simulations are performed to study the effects of Grashof number ranging from 104 to 106, hot source length from 0.1… More >

Annamaria Mazzia¹, Giorgio Pini¹, Flavio Sartoretto²

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 183-210, 2012, DOI:10.3970/cmes.2012.088.183

Abstract Pure meshless techniques are promising methods for solving Partial Differential Equations (PDE). They alleviate difficulties both in designing discretization meshes, and in refining/coarsening, a task which is demanded e.g. in adaptive strategies. Meshless Local Petrov Galerkin (MLPG) methods are pure meshless techniques that receive increasing attention. Very recently, new methods, called Direct MLPG (DMLPG), have been proposed. They rely upon approximating PDE via the Generalized Moving Least Square method. DMLPG methods alleviate some difficulties of MLPG, e.g. numerical integration of tricky, non-polynomial factors, in weak forms. DMLPG techniques require lower computational costs respect to their MLPG counterparts. In this paper… More >

J. Wang¹, E. Andres, P. Simacek, S.G.Advani

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 155-182, 2012, DOI:10.3970/cmes.2012.088.155

Abstract In Liquid Composite Molding (LCM) processes, the process design requires an infusion and venting scheme which will saturate all the empty spaces between the fibers during mold filling resulting in a composite part without voids. However, the inherent material and process variability can change the filling patterns significantly which complicate this task. The objective of this work is to develop methodologies and tools to automate infusion process design and integrate it within the CAD design environment. The methodologies and algorithms developed examine the designed part geometry and material layups for ease of manufacturing with feasible infusion schemes by accounting for… More >

J.Y. Zhang¹, F.Z. Wang^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 141-154, 2012, DOI:10.3970/cmes.2012.088.141

Abstract The boundary knot method (BKM) is a kind of boundary-type meshless method, only boundary points are needed in the solution process. Since the BKM is mathematically simple and easy to implement, it is superior in dealing with Helmholtz problems with high wavenumbers and high dimensional problems. In this paper, we give an overview of the traditional BKM with collocation approach and provide three novel approaches for the BKM, as far as they are relevant for the other boundary-type techniques. The promising research directions are expected from an improved BKM aspect. More >

G. Trapani¹, T. Kipouros¹, A. M. Savill¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 107-140, 2012, DOI:10.3970/cmes.2012.088.107

Abstract This paper presents the development and the application of a multi-objective optimization framework for the design of two-dimensional multi-element high-lift airfoils. An innovative and efficient optimization algorithm, namely Multi-Objective Tabu Search (MOTS), has been selected as core of the framework. The flow-field around the multi-element configuration is simulated using the commercial computational fluid dynamics (cfd) suite Ansys cfx. Elements shape and deployment settings have been considered as design variables in the optimization of the Garteur A310 airfoil, as presented here. A validation and verification process of the cfd simulation for the Garteur airfoil is performed using available wind tunnel data.… More >

Nizami Gasilov¹, Sahin Emrah Amrahov² , Afet Golayoglu Fatullayev ˇ ¹, Halil Ibrahim Karaka¸s¹, Ömer Akın³

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 93-106, 2012, DOI:10.3970/cmes.2012.088.093

Abstract Uncertainties in some parameters of problems of Leontief input-output analysis lead naturally to fuzzy linear systems. In this work, we consider input-output model, where the technology matrix is crisp and the vector of final outputs is fuzzy. The model is expressed by a fuzzy linear system with crisp matrix and with fuzzy right-hand side vector. We apply a geometric method for solving the system. The method finds the solution in the form of a fuzzy set of vectors. The solution set is shown to be a parallelepiped in coordinate space and is expressed by an explicit formula. The features of… More >

Chih-Wen Chang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 65-92, 2012, DOI:10.3970/cmes.2012.088.065

Abstract We utilize a regularized integral equation scheme to resolve the three-dimensional backward advection-dispersion equation (BADE) for identifying the groundwater pollution source identification problems in this research. First, the Fourier series expansion method is employed to estimate the concentration field C(x, y, z, t) at any time t < T. Second, we contemplate a direct regularization by adding an extra term a(x, y, z, 0) to transform a second-kind Fredholm integral equation for C(x, y, z, 0). The termwise separable property of the kernel function permits us to acquire a closed-form regularized solution. In addition, a tactic to determine the regularization… More >

Hideki Fujii¹, Shinobu Yoshimura¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.1, pp. 49-64, 2012, DOI:10.3970/cmes.2012.088.049

Abstract Recently, global warming issues have been discussed all over the world. Of the total amount of CO₂ emitted in Japan, a transportation sector is responsible for 20%. In the transportation sector, 90% of the emission is due to road traffic. This amount must be reduced drastically to realize a low-carbon society. To do so, various measures have been discussed, and the effects of the measures must be estimated quantitatively. In conventional measurement methods, the amount of vehicle emission is simply calculated by multiplying travel distance or gasoline consumption by a specified emission coefficient. Such an approach neglects the effects of… More >

Chein-Shan Liu¹, Su-Ying Zhang², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.1, pp. 29-48, 2012, DOI:10.3970/cmes.2012.088.029

Abstract With a detailed investigation of n linear algebraic equations Bx=b, we find that the scaled residual dynamics for y∈Sⁿ⁻¹ is equipped with four structures: the Jordan dynamics, the rotation group SO(n), a generalized Hamiltonian formulation, as well as a metric bracket system. Therefore, it is the first time that we can compute the steplength used in the iterative method by a novel algorithm based on the Jordan structure. The algorithms preserving the length of y are developed as the structure preserving algorithms (SPAs), which can significantly accelerate the convergence speed and are robust enough against the noise in the numerical… More >