A. P. S. Selvadurai¹, Wenjun Dong

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 41-54, 2006, DOI:10.3970/cmes.2006.012.041

Abstract A Fourier analysis conducted on both the spatial and the temporal discretizations of the governing partial differential equation shows that the Courant number as well as the time marching scheme have significant influences on the numerical behaviour of a Modified Least Squares (MLS) method for the solution of the advection equation. The variations of the amplification factor and the relative phase velocity with the Courant number and the dimensionless wave number indicate that when Courant number is equal to unity, the MLS method with the specified time-weighting and upwind function gives accurate results. This conclusion is confirmed by the numerical… More >

Gongsheng Li¹, De Yao², Hengyi Jiang³, Xianzheng Jia¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.2, pp. 83-108, 2011, DOI:10.3970/cmes.2011.074.083

Abstract This paper deals with an inverse problem of determining a time-depen -dent reaction coefficient arising from a disturbed soil-column infiltrating experiment based on measured breakthrough data. A purpose of doing such experiment is to simulate and study transport behaviors of contaminants when they vertically penetrating through the soils. Data compatibility of the inverse problem is discussed showing a sufficient condition to the solution's monotonicity and positivity with the help of an adjoint problem. Furthermore, an optimal perturbation regularization algorithm is applied to solve the inverse problem, and two typical numerical examples are presented to support the inversion algorithm. Finally, transport… More >

Ruben Avila¹, Zhidong Han², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.4, pp. 363-396, 2011, DOI:10.3970/cmes.2011.071.363

Abstract The two dimensional steady state Stokes equations are solved by using a novel MLPG-Mixed Finite Volume method, that is based on the independent meshless interpolations of the deviatoric velocity strain tensor, the volumetric velocity strain tensor, the velocity vector and the pressure. The pressure field directly obtained from this method does not suffer from the malady of checker-board patterns. Numerical simulations of the flow field, and trajectories of passive fluid elements in a new complex Stokes flow are also presented. The new flow geometry consists of three coaxial cylinders two of smaller diameter, that steadily rotate independently, inside a third… More >

N.Kringos¹, A.Scarpas¹, A.P.S. Selvadurai²

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.3, pp. 147-160, 2008, DOI:10.3970/cmes.2008.028.147

Abstract This paper presents a numerical approach for the modeling of water flow induced mastic erosion from a permeable asphalt mix and is part of an ongoing effort to model moisture-induced damage in asphalt mixes. Due to the complex composite structure of asphalt mixtures, moisture can infiltrate in various ways into the components and have an adverse effect on its mechanical performance. Depending on the gradation of the asphalt aggregates and the mixing procedure, asphalt structures with a variable permeability may result. Open-graded asphalt mixes are designed with a high interconnected air void content to serve as a drainage layer on… More >

Chein-Shan Liu¹, Chih-Wen Chang², Jiang-Ren Chang^2,3

CMC-Computers, Materials & Continua, Vol.9, No.2, pp. 111-136, 2009, DOI:10.3970/cmc.2009.009.111

Abstract In this paper, we take the advantage of an analytical method to solve the advection-dispersion equation (ADE) for identifying the contamination problems. First, the Fourier series expansion technique is employed to calculate the concentration field C(x, t) at any time t< T. Then, we consider a direct regularization by adding an extra term αC(x,0) on the final condition to carry off a second kind Fredholm integral equation. The termwise separable property of the kernel function permits us to transform itinto a two-point boundary value problem. The uniform convergence and error estimate of the regularized solution C^α(x,t) are provided and a… More >

Chih-Wen Chang¹, Chein-Shan Liu²

CMC-Computers, Materials & Continua, Vol.17, No.1, pp. 19-40, 2010, DOI:10.3970/cmc.2010.017.019

Abstract In this study, we employ a semi-analytical approach to solve a two-dimensional advection-dispersion equation (ADE) for identifying the contamination problems. First, the Fourier series expansion technique is used to calculate the concentration field C(x, y, t) at any time t < T. Then, we ponder a direct regularization by adding an extra termaC(x, y, 0) on the final time data C(x, y, T), to reach a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us obtaining a closed-form solution of the Fourier coefficients. A strategy to choose the regularization parameter is offered. The solver utilized in… More >