Jinsheng Wang¹, Liqing Liu², Yiming Chen², Lechun Liu², Dayan Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.104, No.1, pp. 69-85, 2015, DOI:10.3970/cmes.2015.104.069

Abstract The aim of this paper is to seek the numerical solution of a class of variable order fractional integral-differential equation in terms of Bernstein polynomials. The fractional derivative is described in the Caputo sense. Four kinds of operational matrixes of Bernstein polynomials are introduced and are utilized to reduce the initial equation to the solution of algebraic equations after dispersing the variable. By solving the algebraic equations, the numerical solutions are acquired. The method in general is easy to implement and yields good results. Numerical examples are provided to demonstrate the validity and applicability of the method. More >

P. Redou¹, L. Gaubert¹, G. Desmeulles¹, P-A. Béal², C. Le Gal², V. Rodin³

CMES-Computer Modeling in Engineering & Sciences, Vol.70, No.1, pp. 11-40, 2010, DOI:10.3970/cmes.2010.070.011

Abstract Multi Interaction Systems, used in the context of Virtual Reality, are dedicated to real-time interactive simulations. They open the way to the in virtuo experimentation, especially useful in the domain of biochemical kinetics. To this purpose, chaotic and asynchronous scheduling of autonomous processes is based upon desynchronization of phenomena involved in the system. It permits interactivity, especially the capability to add or remove phenomena in the course of a simulation. It provides methods of resolution of ordinary differential systems and partial derivative equations. Proofs of convergence for these methods have been established, but the problem of absolute stability, although it… More >

A. Brancati¹, M. H. Aliabadi¹, A. Milazzo²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 1-24, 2011, DOI:10.3970/cmes.2011.078.001

Abstract This paper presents an improved adaptive cross approximation (ACA) approach developed in conjunction with the Hierarchical format matrix and the GMRES solver. A novel scheme to generate the cluster tree (based upon preliminary considerations of the prescribed boundary conditions) and an improved ACA algorithm (approximating the system matrix for mixed Robin conditions) are described. The asymptotic smoothness property of a kernel generated by a linear combination of two asymptotic smooth kernels is demonstrated. Numerical results show the new approach to be up to 50% faster than the conventional ACA approach. More >

Roland Martin¹, Carlos Couder-Castaneda¹

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 47-78, 2010, DOI:10.3970/cmes.2010.063.047

Abstract We develop an unsplit convolutional perfectly matched layer (CPML) technique to absorb efficiently compressible viscous flows and their related supersonic or subsonic regimes at the outer boundary of a distorted computational domain. More particularly subsonic outgoing flows or subsonic wall-boundary layers close to the PML are well absorbed, which is difficult to obtain without creating numerical instabilities over long time periods. This new PML (CPML) introduces the calculation of auxiliary memory variables at each time step and allows an unsplit formulation of the PML. Damping functions involving a high shift in the frequency domain allow a much better absorption of… More >

R. Martin¹, D. Komatitsch^1,2, S. D. Gedney³, E. Bruthiaux^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.1, pp. 17-42, 2010, DOI:10.3970/cmes.2010.056.017

Abstract Unsplit convolutional perfectly matched layers (CPML) for the velocity and stress formulation of the seismic wave equation are classically computed based on a second-order finite-difference time scheme. However it is often of interest to increase the order of the time-stepping scheme in order to increase the accuracy of the algorithm. This is important for instance in the case of very long simulations. We study how to define and implement a new unsplit non-convolutional PML called the Auxiliary Differential Equation PML (ADE-PML), based on a high-order Runge-Kutta time-stepping scheme and optimized at grazing incidence. We demonstrate that when a second-order time-stepping… More >

C. J. F. Solano, A. Costales, E. Francisco, A. Martín Pendás, Miguel A. Blanco¹, K.-C. Lau, H. He, Ravindra Pandey²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 143-156, 2008, DOI:10.3970/cmes.2008.024.143

Abstract The buckling of hexagonal layers in bulk and nanostructures of AlN is analyzed in the framework of atomistic and first principles techniques. At ambient conditions, the wurtzite structure (B4) of AlN consists of buckled hexagons. On the other hand, a non-buckled B_k structure is found to be metastable at zero pressure, being favored at higher pressures. It is suggested that the energy ordering of B4 and B_k may change in finite systems; an assertion tested in this study by considering finite slabs, nanobelts, and nanorings, and comparing the results with the previous studies on small clusters, and periodic nanostructures. We… More >

Shamsh Tabrez, Mira Mitra, S. Gopalakrishnan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.1, pp. 77-90, 2007, DOI:10.3970/cmes.2007.022.077

Abstract In this paper, wave propagation is studied in degraded composite beam due to moisture absorption. The obtained wave responses are then used for diagnosis of the degraded zone. Moisture absorption causes an irreversible hygrothermal deterioration of the material. The change in temperature and moisture absorption changes the mechanical properties. Thus this affects the structure in dimensional stability as well as material degradation due to reduction in mechanical properties. Here, the composite beam is modeled as Timoshenko beam using wavelet based spectral finite element (WSFE) method. The WSFE technique is especially tailored for simulation of wave propagation. It involves Daubechies scaling… More >

S.-P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 227-236, 2001, DOI:10.3970/cmes.2001.002.227

Abstract Problems defined on infinite domains must be treated on a finite computational domain. The treatment of the artificially placed boundaries (usually referred to as open boundaries) of such domain truncations can be quite subtle; an over truncation would normally result in large, undesirable reflection of signals back to the computational domain whereas an under truncation would imply an injudicious use of computational resources. In particular, problems occur when strongly nonlinear free surface waves generated in a numerical wave tank are passing through such an open boundary.
In this paper, some recent numerical test results of an innovative treatment of… 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 in two dimensions. More >

Mohamed Ammar Abbassi¹, Bouchmel Mliki¹, Ridha Djebali^1,2

FDMP-Fluid Dynamics & Materials Processing, Vol.13, No.2, pp. 59-83, 2017, DOI:10.3970/fdmp.2017.013.059

Abstract This article reports a numerical study of nanoparticle Brownian motion and magnetic field effects by natural convection in a nanofluid-filled open cavity with non uniform boundary condition. Lattice Boltzmann Method (LBM) is used to simulate nanofluid flow and heat transfer. The effective thermal conductivity and viscosity of nanofluid are calculated by KKL (Koo-Kleinstreuer-Li) correlation. In this model effect of Brownian motion on the effective thermal conductivity and effective viscosity is considered and examined. Simulations have been carried out for the pertinent parameters in the following ranges: Rayleigh number (Ra=10³−10⁶), Hartmann number (Ha=0-60), nanoparticle volume concentration (Φ=0–0.04) and heat generation or… More >