Abstract: The finite-difference time-domain (FDTD) method is
one of the most widely used computational methods in
electromagnetic. This paper describes the design of two-dimensional
(2D) FDTD simulation software for transverse magnetic (TM)
polarization using Berenger's split-field perfectly matched layer
(PML) formulation. The software is developed using Matlab
programming language. Numerical examples validate the software.
Abstract: Chemical reaction and diffusion are important phenomena in quantitative neurobiology and biophysics. The knowledge of the dynamics of calcium Ca2+ is very important in cellular physiology because Ca2+ binds to many proteins and regulates their activity and interactions Calcium waves propagate inside cells due to a regenerative mechanism known as calcium-induced calcium release. Buffer-mediated calcium diffusion in the cytosol plays a crucial role in the process. A mathematical model has been developed for calcium waves by assuming the buffers are in equilibrium with calcium i.e., the rapid buffering approximation for a one dimensional unsteady state case. This model incorporates important physical and physiological parameters like dissociation rate, diffusion rate, total buffer concentration and influx. The finite difference method has been employed to predict [Ca2+] and buffer concentration time course regardless of the calcium influx. The comparative studies of the effect of the rapid buffered diffusion and kinetic parameters of the model on the concentration time course have been performed.
Abstract: Double-diffusive natural convection in an open top
square cavity and heated from the side is studied numerically.
Constant temperatures and concentration are imposed along the right
and left walls while the heat balance at the surface is assumed to obey
Newton-s law of cooling. The finite difference method is used to
solve the dimensionless governing equations. The numerical results
are reported for the effect of Marangoni number, Biot number and
Prandtl number on the contours of streamlines, temperature and
concentration. The predicted results for the average Nusselt number
and Sherwood number are presented for various parametric
conditions. The parameters involved are as follows; the thermal
Marangoni number, 0 ≤ MaT ≤1000 , the solutal Marangoni number,
0 1000 c ≤ Ma ≤ , the Biot number, 0 ≤ Bi ≤ 6 , Grashof number,
5 Gr = 10 and aspect ratio 1. The study focused on both flows; thermal
dominated, N = 0.8 , and compositional dominated, N = 1.3 .
Abstract: A new numerical method for solving the twodimensional,
steady, incompressible, viscous flow equations on a
Curvilinear staggered grid is presented in this paper. The proposed
methodology is finite difference based, but essentially takes
advantage of the best features of two well-established numerical
formulations, the finite difference and finite volume methods. Some
weaknesses of the finite difference approach are removed by
exploiting the strengths of the finite volume method. In particular,
the issue of velocity-pressure coupling is dealt with in the proposed
finite difference formulation by developing a pressure correction
equation in a manner similar to the SIMPLE approach commonly
used in finite volume formulations. However, since this is purely a
finite difference formulation, numerical approximation of fluxes is
not required. Results obtained from the present method are based on
the first-order upwind scheme for the convective terms, but the
methodology can easily be modified to accommodate higher order
differencing schemes.
Abstract: In this paper, numerical solution for the generalized Rosenau-Burgers equation is considered and Crank-Nicolson finite difference scheme is proposed. Existence of the solutions for the difference scheme has been shown. Stability, convergence and priori error estimate of the scheme are proved. Numerical results demonstrate that the scheme is efficient and reliable.
Abstract: Sampling and analysis of leachate from Bhalaswa
landfill and groundwater samples from nearby locations, clearly
indicated the likely contamination of groundwater due to landfill
leachate. The results of simulation studies carried out for the
migration of Chloride from landfill shows that the simulation results
are in consonance with the observed concentration of Chloride in the
vicinity of landfill facility. The solid waste disposal system presently
being practiced in Delhi consists of mere dumping of wastes
generated, at three locations Bhalaswa, Ghazipur, and Okhla without
any regard to proper care for the protection of surrounding
environment. Bhalaswa landfill site in Delhi, which is being operated
as a dump site, is expected to become cause of serious groundwater
pollution in its vicinity. The leachate from Bhalaswa landfill was
found to be having a high concentration of chlorides, as well as DOC,
COD. The present study was undertaken to determine the likely
concentrations of principle contaminants in the groundwater over a
period of time due to the discharge of such contaminants from
landfill leachates to the underlying groundwater. The observed
concentration of chlorides in the groundwater within 75m of the
radius of landfill facility was found to be in consonance with the
simulated concentration of chloride in groundwater considering one
dimensional transport model, with finite mass of contaminant source.
Governing equation of contaminant transport involving advection and
diffusion-dispersion was solved in matlab7.0 using finite difference
method.
Abstract: In this paper, collocation based cubic B-spline and
extended cubic uniform B-spline method are considered for
solving one-dimensional heat equation with a nonlocal initial
condition. Finite difference and θ-weighted scheme is used for
time and space discretization respectively. The stability of the
method is analyzed by the Von Neumann method. Accuracy of
the methods is illustrated with an example. The numerical results
are obtained and compared with the analytical solutions.
Abstract: In the present paper, an improved initial value
numerical technique is presented to analyze the free vibration of
symmetrically laminated rectangular plate. A combination of the
initial value method (IV) and the finite differences (FD) devices is
utilized to develop the present (IVFD) technique. The achieved
technique is applied to the equation of motion of vibrating laminated
rectangular plate under various types of boundary conditions. Three
common types of laminated symmetrically cross-ply, orthotropic and
isotropic plates are analyzed here. The convergence and accuracy of
the presented Initial Value-Finite Differences (IVFD) technique have
been examined. Also, the merits and validity of improved technique
are satisfied via comparing the obtained results with those available
in literature indicating good agreements.
Abstract: A numerical study has been carried out to investigate
the heat transfer by natural convection of nanofluid taking Cu as
nanoparticles and the water as based fluid in a three dimensional
annulus enclosure filled with porous media (silica sand) between two
horizontal concentric cylinders with 12 annular fins of 2.4mm
thickness attached to the inner cylinder under steady state conditions.
The governing equations which used are continuity, momentum and
energy equations under an assumptions used Darcy law and
Boussinesq-s approximation which are transformed to dimensionless
equations. The finite difference approach is used to obtain all the
computational results using the MATLAB-7. The parameters affected
on the system are modified Rayleigh number (10 ≤Ra*≤ 1000), fin
length Hf (3, 7 and 11mm), radius ratio Rr (0.293, 0.365 and 0.435)
and the volume fraction(0 ≤ ¤ò ≤ 0 .35). It was found that the
average Nusselt number depends on (Ra*, Hf, Rr and φ). The results
show that, increasing of fin length decreases the heat transfer rate and
for low values of Ra*, decreasing Rr cause to decrease Nu while for
Ra*
greater than 100, decreasing Rr cause to increase Nu and adding
Cu nanoparticles with 0.35 volume fraction cause 27.9%
enhancement in heat transfer. A correlation for Nu in terms of Ra*,
Hf and φ, has been developed for inner hot cylinder.
Abstract: The spectral action balance equation is an equation that
used to simulate short-crested wind-generated waves in shallow water
areas such as coastal regions and inland waters. This equation consists
of two spatial dimensions, wave direction, and wave frequency which
can be solved by finite difference method. When this equation with
dominating propagation velocity terms are discretized using central
differences, stability problems occur when the grid spacing is chosen
too coarse. In this paper, we introduce the splitting modified donorcell
scheme for avoiding stability problems and prove that it is
consistent to the modified donor-cell scheme with same accuracy. The
splitting modified donor-cell scheme was adopted to split the wave
spectral action balance equation into four one-dimensional problems,
which for each small problem obtains the independently tridiagonal
linear systems. For each smaller system can be solved by direct or
iterative methods at the same time which is very fast when performed
by a multi-cores computer.
Abstract: This paper uses quasi-steady molecular statics model
and diamond tool to carry out simulation temperature rise of nanoscale
orthogonal cutting single-crystal silicon. It further qualitatively
analyzes temperature field of silicon workpiece without considering
heat transfer and considering heat transfer. This paper supposes that
the temperature rise of workpiece is mainly caused by two heat sources:
plastic deformation heat and friction heat. Then, this paper develops a
theoretical model about production of the plastic deformation heat and
friction heat during nanoscale orthogonal cutting. After the increased
temperature produced by these two heat sources are added up, the
acquired total temperature rise at each atom of the workpiece is
substituted in heat transfer finite difference equation to carry out heat
transfer and calculates the temperature field in each step and makes
related analysis.
Abstract: The present paper develops and validates a numerical procedure for the calculation of turbulent combustive flow in converging and diverging ducts and throuh simulation of the heat transfer processes, the amount of production and spread of Nox pollutant has been measured. A marching integration solution procedure employing the TDMA is used to solve the discretized equations. The turbulence model is the Prandtl Mixing Length method. Modeling the combustion process is done by the use of Arrhenius and Eddy Dissipation method. Thermal mechanism has been utilized for modeling the process of forming the nitrogen oxides. Finite difference method and Genmix numerical code are used for numerical solution of equations. Our results indicate the important influence of the limiting diverging angle of diffuser on the coefficient of recovering of pressure. Moreover, due to the intense dependence of Nox pollutant to the maximum temperature in the domain with this feature, the Nox pollutant amount is also in maximum level.
Abstract: Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modelling reaction diffusion for such branches of sciences. In this paper, a numerical solution for the one-dimensional hyperbolic telegraph equation by using the collocation method using the septic splines is proposed. The scheme works in a similar fashion as finite difference methods. Test problems are used to validate our scheme by calculate L2-norm and L∞-norm. The accuracy of the presented method is demonstrated by two test problems. The numerical results are found to be in good agreement with the exact solutions.
Abstract: Heating is inevitable in any bearing operation. This
leads to not only the thinning of the lubricant but also could lead to a
thermal deformation of the bearing. The present work is an attempt to
analyze the influence of thermal deformation on the thermohydrodynamic
lubrication of infinitely long tilted pad slider rough
bearings. As a consequence of heating the slider is deformed and is
assumed to take a parabolic shape. Also the asperities expand leading
to smaller effective film thickness. Two different types of surface
roughness are considered: longitudinal roughness and transverse
roughness. Christensen-s stochastic approach is used to derive the
Reynolds-type equations. Density and viscosity are considered to be
temperature dependent. The modified Reynolds equation, momentum
equation, continuity equation and energy equation are decoupled and
solved using finite difference method to yield various bearing
characteristics. From the numerical simulations it is observed that the
performance of the bearing is significantly affected by the thermal
distortion of the slider and asperities and even the parallel sliders
seem to carry some load.
Abstract: In this paper, we have combined some spatial derivatives with the optimised time derivative proposed by Tam and Webb in order to approximate the linear advection equation which is given by = 0. Ôêé Ôêé + Ôêé Ôêé x f t u These spatial derivatives are as follows: a standard 7-point 6 th -order central difference scheme (ST7), a standard 9-point 8 th -order central difference scheme (ST9) and optimised schemes designed by Tam and Webb, Lockard et al., Zingg et al., Zhuang and Chen, Bogey and Bailly. Thus, these seven different spatial derivatives have been coupled with the optimised time derivative to obtain seven different finite-difference schemes to approximate the linear advection equation. We have analysed the variation of the modified wavenumber and group velocity, both with respect to the exact wavenumber for each spatial derivative. The problems considered are the 1-D propagation of a Boxcar function, propagation of an initial disturbance consisting of a sine and Gaussian function and the propagation of a Gaussian profile. It is known that the choice of the cfl number affects the quality of results in terms of dissipation and dispersion characteristics. Based on the numerical experiments solved and numerical methods used to approximate the linear advection equation, it is observed in this work, that the quality of results is dependent on the choice of the cfl number, even for optimised numerical methods. The errors from the numerical results have been quantified into dispersion and dissipation using a technique devised by Takacs. Also, the quantity, Exponential Error for Low Dispersion and Low Dissipation, eeldld has been computed from the numerical results. Moreover, based on this work, it has been found that when the quantity, eeldld can be used as a measure of the total error. In particular, the total error is a minimum when the eeldld is a minimum.
Abstract: A steady two-dimensional magnetohydrodynamics
flow and heat transfer over a stretching vertical sheet influenced by
radiation and porosity is studied. The governing boundary layer
equations of partial differential equations are reduced to a system of
ordinary differential equations using similarity transformation. The
system is solved numerically by using a finite difference scheme
known as the Keller-box method for some values of parameters,
namely the radiation parameter N, magnetic parameter M, buoyancy
parameter l , Prandtl number Pr and permeability parameter K. The
effects of the parameters on the heat transfer characteristics are
analyzed and discussed. It is found that both the skin friction
coefficient and the local Nusselt number decrease as the magnetic
parameter M and permeability parameter K increase. Heat transfer
rate at the surface decreases as the radiation parameter increases.
Abstract: A numerical study on the heat transfer in the thermal
barrier coatings and the substrates of a parallel-plate enclosure is
carried out. Some of the thermal barrier coatings, such as ceramics, are
semitransparent and are of interest for high-temperature applications
where radiation effects are significant. The radiative transfer equations
and the energy equations are solved by using the discrete ordinates
method and the finite difference method. Illustrative results are
presented for temperature distributions in the coatings and the opaque
walls under various heating conditions. The results show that the
temperature distribution is more uniform in the interior portion of each
coating away from its boundary for the case with a larger average of
varying refractive index and a positive gradient of refractive index
enhances radiative transfer to the substrates.
Abstract: In this paper, growth and collapse of a vapour bubble
generated due to a local energy input inside a rigid cylinder and in
the absence of buoyancy forces is investigated using Boundary
Integral Equation Method and Finite Difference Method .The fluid is
treated as potential flow and Boundary Integral Equation Method is
used to solve Laplace-s equation for velocity potential. Different
ratios of the diameter of the rigid cylinder to the maximum radius of
the bubble are considered. Results show that during the collapse
phase of the bubble inside a vertical rigid cylinder, two liquid micro
jets are developed on the top and bottom sides of the vapour bubble
and are directed inward. It is found that by increasing the ratio of the
cylinder diameter to the maximum radius of the bubble, the rate of
the growth and collapse phases of the bubble increases and the life
time of the bubble decreases.
Abstract: Simultaneous transient conduction and radiation heat
transfer with heat generation is investigated. Analysis is carried out
for both steady and unsteady situations. two-dimensional gray
cylindrical enclosure with an absorbing, emitting, and isotropically
scattering medium is considered. Enclosure boundaries are assumed
at specified temperatures. The heat generation rate is considered
uniform and constant throughout the medium. The lattice Boltzmann
method (LBM) was used to solve the energy equation of a transient
conduction-radiation heat transfer problem. The control volume finite
element method (CVFEM) was used to compute the radiative
information. To study the compatibility of the LBM for the energy
equation and the CVFEM for the radiative transfer equation, transient
conduction and radiation heat transfer problems in 2-D cylindrical
geometries were considered. In order to establish the suitability of the
LBM, the energy equation of the present problem was also solved
using the the finite difference method (FDM) of the computational
fluid dynamics. The CVFEM used in the radiative heat transfer was
employed to compute the radiative information required for the
solution of the energy equation using the LBM or the FDM (of the
CFD). To study the compatibility and suitability of the LBM for the
solution of energy equation and the CVFEM for the radiative
information, results were analyzed for the effects of various
parameters such as the boundary emissivity. The results of the LBMCVFEM
combination were found to be in excellent agreement with
the FDM-CVFEM combination. The number of iterations and the
steady state temperature in both of the combinations were found
comparable. Results are found for situations with and without heat
generation. Heat generation is found to have significant bearing on
temperature distribution.
Abstract: In this work the numerical simulation of transient heat
transfer in a cylindrical probe is done. An experiment was conducted
introducing a steel cylinder in a heating chamber and registering its
surface temperature along the time during one hour. In parallel, a
mathematical model was solved for one dimension transient heat
transfer in cylindrical coordinates, considering the boundary
conditions of the test. The model was solved using finite difference
method, because the thermal conductivity in the cylindrical steel bar
and the convection heat transfer coefficient used in the model are
considered temperature dependant functions, and both conditions
prevent the use of the analytical solution. The comparison between
theoretical and experimental results showed the average deviation is
below 2%. It was concluded that numerical methods are useful in
order to solve engineering complex problems. For constant k and h,
the experimental methodology used here can be used as a tool for
teaching heat transfer in mechanical engineering, using mathematical
simplified models with analytical solutions.