Yiming Chen, Mingxu Yi, Chen Chen, Chunxiao Yu

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.6, pp. 639-654, 2012, DOI:10.3970/cmes.2012.083.639

Abstract In this paper, Bernstein polynomials method is proposed for the numerical solution of a class of space-time fractional convection-diffusion equation with variable coefficients. This method combines the definition of fractional derivatives with some properties of Bernstein polynomials and are dispersed the coefficients efficaciously. The main characteristic behind this method is that the original problem is translated into a Sylvester equation. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples show that the method is effective. 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 >

Han-Taw Chen¹, Li-Shie Liu¹, Shin-Ku Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.2, pp. 155-178, 2010, DOI:10.3970/cmes.2010.065.155

Abstract The finite difference method in conjunction with the least-squares scheme and experimental measured temperatures is proposed to solve a two-dimensional steady-state inverse heat conduction problem in order to predict the natural-convection heat transfer coefficient under the isothermal situation h^−iso from a three fin array mounted on a horizontal plate and fin efficiency η_f for various values of the fin spacing and fin height. The measured fin temperatures and ambient temperature are obtained from the present experimental apparatus conducted in a small wind tunnel. The heat transfer coefficient on a fin is non-uniform for the present problem, and its functional… More >

Sung-Eun Kim¹, Shin Hyung Rhee²

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.3, pp. 291-310, 2010, DOI:10.3970/cmes.2010.062.291

Abstract Turbulent flow past a finite wing has been computed to assess the fidelity of modern computational fluid dynamics in predicting tip vortex flows. The efficacy of a feature-adaptive local mesh refinement to resolve the steep gradients in the flow field near the tip vortex is demonstrated. The impact of turbulence modeling is evaluated using several popular eddy viscosity models and a Reynolds stress transport model. The results indicate that the combination of a computational mesh with an adequate resolution, high-order spatial discretization scheme along with the use of advanced turbulence models can predict tip vortex flows with acceptable accuracy. More >

Bart Bergen¹, Bert Van Genechten¹, Dirk Vandepitte¹, Wim Desmet¹

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.2, pp. 155-176, 2010, DOI:10.3970/cmes.2010.061.155

Abstract The Wave Based Method (WBM) is a numerical prediction technique for Helmholtz problems. It is an indirect Trefftz method using wave functions, which satisfy the Helmholtz equation, for the description of the dynamic variables. In this way, it avoids both the large systems and the pollution errors that jeopardize accurate element-based predictions in the mid-frequency range. The enhanced computational efficiency of the WBM as compared to the element-based methods has been proven for the analysis of both three-dimensional bounded and two-dimensional unbounded problems. This paper presents an extension of the WBM to the application of three-dimensional acoustic scattering and radiation… More >

Paweł Dłotko¹, Ruben Specogna²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.3, pp. 247-278, 2010, DOI:10.3970/cmes.2010.060.247

Abstract The systematic potential design is of high importance in computational electromagnetics. For example, it is well known that when the efficient eddy-current formulations based on a magnetic scalar potential are employed in problems which involve conductive regions with holes, the so-calledthick cutsare needed to make the boundary value problem well defined. Therefore, a considerable effort has been invested over the past twenty-five years to develop fast and general algorithms to compute thick cuts automatically. Nevertheless, none of the approaches proposed in literature meet all the requirements of being automatic, computationally efficient and general. In this paper, an automatic, computationally efficient… More >

Han Taw Chen¹, Shao Lun Sun¹, Hui Chen Huang¹, Sen Yung Lee^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.2, pp. 107-126, 2010, DOI:10.3970/cmes.2010.059.107

Abstract A new solution method is proposed to develop the analytic closed form solution for the one dimensional heat conduction with one mixed type boundary condition and general time dependent heat convection coefficient for the first time. The solution method is the combination of an extension of the shifting function method developed by Lee and his colleagues and a series expansion. It is shown that the solution is simple and accurate. The convergence of the present analysis is very fast. One can find that when the dimensionless Fourier number is greater than 0.2, the error for the one term approximation solution… More >

Arif Amirov¹, Zekeriya Ustaoglu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.3, pp. 283-300, 2009, DOI:10.3970/cmes.2009.054.283

Abstract We present the solvability of a two space dimensional coefficient inverse problem for a transport-like equation and investigate the approximate solution of this problem with the use of centered difference formulas and a symbolic approximation method. Since this inverse problem is overdetermined, which is the main difficulty in studying of its solvability, it is replaced by a related determined one by using some extension of the class of unknown functions. More >

Dan Gordon¹, Rachel Gordon²

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.1, pp. 23-46, 2009, DOI:10.3970/cmes.2009.053.023

Abstract Linear systems with very large off-diagonal elements and discontinuous coefficients (LODC systems) arise in some modeling cases, such as those involving heterogeneous media. Such problems are usually solved by domain decomposition methods, but these can be difficult to implement on unstructured grids or when the boundaries between subdomains have a complicated geometry. Gordon and Gordon have shown that Björck and Elfving's (sequential) CGMN algorithm and their own block-parallel CARP-CG are very robust and efficient on strongly convection dominated cases (but without discontinuous coefficients). They have also shown that scaling the equations by dividing each equation by the L2-norm of its… More >

A.R. da Silva¹, A. Silveira-Neto^2,3, D.A. Rade^2,4, R.Francis⁴, E.A. Santos⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.3, pp. 285-304, 2009, DOI:10.3970/cmes.2009.050.285

Abstract In the context of computational fluid dynamics a numerical investigation of incompressible flow around fixed pairs of rigid circular cylinders was carried out. The two-dimensional filtered Navier-Stokes equations with the Smagorinsky sub-grid scale model were solved using a Cartesian non-uniform grid. The immersed Boundary Method with the Virtual Physical Model was used in order to model the presence of two circular cylinders embedded in the flow. The fractional time step method was used to couple pressure and velocity fields. The simulations were carried out for Reynolds number equal to 72,000 for pitch ratio equal to 2 and different arrangements regarding… More >