S.P. Hu¹, D.L. Young², C.M. Fan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.1, pp. 1-20, 2008, DOI:10.3970/cmes.2008.024.001

Abstract In this paper, a novel numerical scheme called (FDMFS), which combines the finite difference method (FDM) and the method of fundamental solutions (MFS), is proposed to simulate the nonhomogeneous diffusion problem with an unsteady forcing function. Most meshless methods are confined to the investigations of nonhomogeneous diffusion equations with steady forcing functions due to the difficulty to find an unsteady particular solution. Therefore, we proposed a FDM with Cartesian grid to handle the unsteady nonhomogeneous term of the equations. The numerical solution in FDMFS is decomposed into a particular solution and a homogeneous solution. The particular solution is constructed using… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.023.075

Abstract In the MLPG_R (Meshless Local Petrove-Galerkin based on Rankine source solution) method, one needs a meshless interpolation scheme for an unknown function to discretise the governing equation. The MLS (moving least square) method has been used for this purpose so far. The MLS method requires inverse of matrix or solution of a linear algebraic system and so is quite time-consuming. In this paper, a new scheme, called simplified finite difference interpolation (SFDI), is devised. This scheme is generally as accurate as the MLS method but does not need matrix inverse and consume less CPU time to evaluate. Although this scheme… More >

A. J. Davies¹, D. Crann¹, S. J. Kane¹, C-H. Lai²

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.2, pp. 79-86, 2007, DOI:10.3970/cmes.2007.018.079

Abstract The solution process for diffusion problems usually involves the time development separately from the space solution. A finite difference algorithm in time requires a sequential time development in which all previous values must be determined prior to the current value. The Stehfest Laplace transform algorithm, however, allows time solutions without the knowledge of prior values. It is of interest to be able to develop a time-domain decomposition suitable for implementation in a parallel environment. One such possibility is to use the Laplace transform to develop coarse-grained solutions which act as the initial values for a set of fine-grained solutions. The… More >

D. Soares Jr.^1,2, W.J. Mansur¹, D.L. Lima³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 19-34, 2007, DOI:10.3970/cmes.2007.017.019

Abstract The present paper discussion is concerned with the development of robust and efficient algorithms to model propagation of interacting acoustic and elastic waves. The paper considers acoustic-elastic, acoustic-acoustic and elastic-elastic partitioned analyses of coupled systems; however, the focus here is the acoustic-elastic coupling considering finite elements and the acoustic-acoustic coupling considering finite elements and finite differences (other coupling procedures can be implemented analogously). One important feature of the algorithms presented is that they allow considering different time-steps for different sub-domains; so it is possible to substantially improve efficiency, accuracy and stability of the central difference time integration algorithm employed here.… More >

Sacia Kachi^1,*, Fatima-zohra Bensouici¹, Nawel Ferroudj¹, Saadoun Boudebous²

FDMP-Fluid Dynamics & Materials Processing, Vol.15, No.2, pp. 89-105, 2019, DOI:10.32604/fdmp.2019.00263

Abstract The objective of the present study is to analyze the laminar mixed convection in a square cavity with moving cooled vertical sidewalls. A constant flux heat source with relative length l is placed in the center of the lower wall while all the other horizontal sides of the cavity are considered adiabatic. The numerical method is based on a finite difference technique where the spatial partial derivatives appearing in the governing equations are discretized using a high order scheme, and time advance is dealt with by a fourth order Runge Kutta method. The Richardson number (Ri), which represents the relative… More >

A. Khechiba¹, Y. Benakcha¹, A. Ghezal¹, P. Spetiri²

FDMP-Fluid Dynamics & Materials Processing, Vol.14, No.2, pp. 137-154, 2018, DOI: 10.3970/fdmp.2018.04054

Abstract This work investigates the dynamic behavior of a pulsatile flow electrically conducting through porous medium in a cylindrical conduit under the influence of a magnetic field. The imposed magnetic field is assumed to be uniform and constant. An exact solution of the equations governing magneto hydro-dynamics (MHD) flow in a conduit has been obtained in the form of Bessel functions. The analytical study has been used to establish an expression between the Hartmann number, Darcy number and the stress coefficient. The numerical method is based on an implicit finite difference time marching scheme using the Thomas algorithm and Gauss Seidel… More >

Anupam Bh,ari¹, Vipin Kumar²

FDMP-Fluid Dynamics & Materials Processing, Vol.10, No.3, pp. 359-375, 2014, DOI:10.3970/fdmp.2014.010.359

Abstract The purpose of this paper is to study the flow characteristics of a ferrofluid revolving through a porous medium with a magnetic-field-dependent viscosity in the presence of a stationary disk. A Finite Difference Method is employed to discretize the set of nonlinear coupled differential equations involved in the problem. The discretized nonlinear equations, in turn, are solved by a Newton method (using MATLAB) taking the initial guess with the help of a PDE Solver. Results displayed in graphical form are used to assess the effect of the variable viscosity and porosity parameters on the velocity components. The displacement thickness of… More >

Paras Ram^1,2, Vikas Kumar³

FDMP-Fluid Dynamics & Materials Processing, Vol.10, No.2, pp. 179-196, 2014, DOI:10.3970/fdmp.2014.010.179

Abstract The present study is carried out to examine the effects of temperature dependent variable viscosity on the three dimensional steady axi-symmetric Ferrohydrodynamic (FHD) boundary layer flow of an incompressible electrically nonconducting magnetic fluid in the presence of a rotating disk. The disk is subjected to an externally applied magnetic field and is maintained at a uniform temperature. The nonlinear coupled partial differential equations governing the boundary layer flow are non dimensionalized using similarity transformations and are reduced to a system of coupled ordinary differential equations. To study the effects of temperature dependent viscosity on velocity profiles and temperature distribution within… More >

Cosmina S. Hogea¹, Bruce T. Murray¹, James A. Sethian^2,3

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.2, pp. 109-130, 2005, DOI:10.3970/fdmp.2005.001.109

Abstract A computational framework for simulating growth and transport in biological materials based on continuum models is proposed. The advantages of the finite difference methodology employed are generality and relative simplicity of implementation. The Cartesian mesh/level set method developed here provides a computational tool for the investigation of a host of transport-based tissue/tumor growth models, that are posed as free or moving boundary problems and may exhibit complicated boundary evolution including topological changes. The methodology is tested here on a widely studied "incompressible flow" type tumor growth model with a numerical implementation in two dimensions; comparisons with results obtained from a… More >

Chia-Ming Fan^1,2, Chein-Shan Liu³, Weichung Yeih¹, Hsin-Fang Chan¹

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 67-86, 2010, DOI:10.3970/cmc.2010.015.067

Abstract In this study, the nonlinear obstacle problems, which are also known as the nonlinear free boundary problems, are analyzed by the scalar homotopy method (SHM) and the finite difference method. The one- and two-dimensional nonlinear obstacle problems, formulated as the nonlinear complementarity problems (NCPs), are discretized by the finite difference method and form a system of nonlinear algebraic equations (NAEs) with the aid of Fischer-Burmeister NCP-function. Additionally, the system of NAEs is solved by the SHM, which is globally convergent and can get rid of calculating the inverse of Jacobian matrix. In SHM, by introducing a scalar homotopy function and… More >