Srinivasan Natesan¹, S. Gowrisankar²

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.4, pp. 245-268, 2012, DOI:10.3970/cmes.2012.088.245

Abstract In this article, we propose a parameter-uniform computational technique to solve singularly perturbed parabolic initial-boundary-value problems exhibiting parabolic layers. The domain is discretized with a uniform mesh on the time direction and a nonuniform mesh obtained via equidistribution of a monitor function for the spatial variable. The numerical scheme consists of the implicit-Euler scheme for the time derivative and the classical central difference scheme for the spatial derivative. Truncation error, and stability analysis are carried out. Error estimates are derived, and numerical examples are presented. More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 179-196, 2006, DOI:10.3970/cmes.2006.015.179

Abstract This paper studies the numerical computations of the second-order singularly perturbed boundary value problems (SPBVPs). In order to depress the singularity we consider a coordinate transformation from the x-domain to the t-domain. The relation between singularity and stiffness is demonstrated, of which the coordinate transformation parameter λ plays a key role to balance these two tendencies. Then we construct a very effective Lie-group shooting method to search the missing initial condition through a weighting factor r ∈ (0,1) in the t-domain formulation. For stabilizing the new method we also introduce two new systems by a translation of the dependent variable.… More >

Wei-Chung Tien¹, Yue-Tzu Yang¹, Cha’o-Kuang Chen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.5, pp. 447-462, 2012, DOI:10.3970/cmes.2012.085.447

Abstract Both buoyancy effects of thermal and mass diffusion in the natural convection flow about a vertical plate are considered in this paper. The non-linear coupled differential governing equations for velocity, temperature and concentration fields are solved by using the homotopy analysis method. Without the need of iteration, the obtained solution is in the form of an infinite power series which indicates those series have high accuracy when comparing it with other-generated by the traditional method. The impact of the Prandtl number, Schmidt number and the buoyancy parameter on the flow are widely discussed in detail. More >

W. A. El-Askary¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.3&4, pp. 203-232, 2011, DOI:10.3970/cmes.2011.074.203

Abstract Large eddy simulation (LES) is a viable and powerful tool to analyze unsteady three- dimensional turbulent flows. In this paper, the method of LES is used to compute a plane turbulent supersonic boundary layer subjected to different pressure gradients. The pressure gradients are generated by allowing the flow to pass in the vicinity of an expansion-compression ramp (inclined backward-facing step with leeward-face angle of 25 degrees) for an upstream Mach number of 2.9. The inflow boundary condition is the main problem for all turbulent wall-bounded flows. An approach to solve this problem is to extract instantaneous velocities, temperature and density… More >

Chein-Shan Liu¹, Weichung Yeih², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.3, pp. 301-324, 2010, DOI:10.3970/cmes.2010.059.301

Abstract When problems in engineering and science are discretized, algebraic equations appear naturally. In a recent paper by Liu and Atluri, non-linear algebraic equations (NAEs) were transformed into a nonlinear system of ODEs, which were then integrated by a method labelled as the Fictitious Time Integration Method (FTIM). In this paper, the FTIM is enhanced, by using the concept of arepellorin the theory ofnonlinear dynamical systems, to situations where multiple-solutions exist. We label this enhanced method as MSFTIM. MSFTIM is applied and illustrated in this paper through solving boundary-layer problems, boundary-value problems, and eigenvalue problems with multiple solutions. More >

Y.M. Zhang¹, Y. Gu¹, J.T. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.45, No.3, pp. 227-248, 2009, DOI:10.3970/cmes.2009.045.227

Abstract The accurate evaluation of nearly singular integrals is one of the major concerned problems in the boundary element method (BEM). Although the current methods have achieved great progress, it is often possible only for problems defined in the simplest geometrical domains when the nearly singular integrals need to be calculated. However, engineering processes occur mostly in complex geometrical domains, and always, involve nonlinearities of the unknown variables and its derivatives. Therefore, effective methods of dealing with nearly singular integrals for such practical problems are necessary and need to be further investigated. In this paper, a general strategy based on a… More >

Briti Sundar Deb¹, Srinivasan Natesan²

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.2, pp. 179-200, 2008, DOI:10.3970/cmes.2008.038.179

Abstract This paper presents an almost second--order uniformly convergent Richardson extrapolation method for convection- dominated coupled system of boundary value problems. First, we solve the system by using the classical finite difference scheme on the layer resolving Shishkin mesh, and then we construct the Richardson approximation solution using the solutions obtained on N and 2N mesh intervals. Second-order parameter--uniform error estimate is derived. The proposed method is applied to a test example for verification of the theoretical results for the case ε ≤ N⁻¹. More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.1, pp. 1-16, 2008, DOI:10.3970/cmes.2008.032.001

Abstract In this paper, we propose a new method to tackle of two famous boundary layer equations in fluid mechanics, namely, the Falkner-Skan and the Blasius equations. We can employ this method to find unknown initial conditions. The pivotal point is based on the erection of a one-step Lie group element$\mathbf {G}(T)$ and the formation of a generalized mid-point Lie group element$\mathbf {G}(r)$. Then, by imposing$\mathbf {G}(T) = \mathbf {G}(r)$ we can seek the missing initial conditions through a minimum discrepancy from the target in terms of a weighting factor$r \in (0, 1)$. Numerical examples are worked out to persuade that… More >

Hayder A. Abdulbari ^1,2, Hassan D. Mahammed¹, Z. Hassan, Wafaa K. Mahmood³

FDMP-Fluid Dynamics & Materials Processing, Vol.12, No.2, pp. 86-101, 2016, DOI:10.3970/fdmp.2016.012.086

Abstract An experimental investigation was carried out to study the response of a turbulent pressure drop fluctuations to longitudinal groove riblets, involved two configurations being triangular and spaced triangular grooves with height 600, 800, 1000 μm and peak to peak spacing 1000 μm and 2000 μm respectively. Experiments were therefore performed at free stream velocity up to 0.44 m/sec, which were corresponding to Reynolds number (Re) 53000. The development of the obtained turbulent layer downstream of the grooves was then compared with the results from the corresponding smooth-wall case. To conclude, the effect of the spaced triangular riblets on the turbulent… More >

S. Aggoune¹, C. Abid², E.H. Amara^1,3

FDMP-Fluid Dynamics & Materials Processing, Vol.11, No.2, pp. 115-125, 2015, DOI:10.3970/fdmp.2015.011.115

Abstract This paper focuses on the vortex formation effect during the application of a laser-fusion cutting technique. This industrial technique is typically associated with the ejection of a film of molten stainless steel blown off by a subsonic laminar jet of nitrogen gas used to assist the process. Without taking into account the transverse movement of the workpiece, we consider a 4 mm thick stainless steel plate. The resulting molten metal flow is assumed to be laminar, steady, viscous and incompressible. The numerical results reveal vortex structures adjacent to the walls at the entrance of the kerf, and a pair of… More >