Abstract: In this paper, we present a fast and accurate numerical scheme for the solution of a Laplace equation with Dirichlet boundary conditions. The non-standard finite difference scheme (NSFD) is applied to construct the numerical solutions of a Laplace equation with two different Dirichlet boundary conditions. The solutions obtained using NSFD are compared with the solutions obtained using the standard finite difference scheme (SFD). The NSFD scheme is demonstrated to be reliable and efficient.
Abstract: In this paper, one dimensional advection diffusion
model is analyzed using finite difference method based on
Crank-Nicolson scheme. A practical problem of filter cake washing
of chemical engineering is analyzed. The model is converted into
dimensionless form. For the grid Ω × ω = [0, 1] × [0, T], the
Crank-Nicolson spatial derivative scheme is used in space domain
and forward difference scheme is used in time domain. The scheme is
found to be unconditionally convergent, stable, first order accurate in
time and second order accurate in space domain. For a test problem,
numerical results are compared with the analytical ones for different
values of parameter.
Abstract: A semi-implicit phase field method for droplet evolution is proposed. Using the phase field Cahn-Hilliard equation, we are able to track the interface in multiphase flow. The idea of a semi-implicit finite difference scheme is reviewed and employed to solve two nonlinear equations, including the Navier-Stokes and the Cahn-Hilliard equations. The use of a semi-implicit method allows us to have larger time steps compared to explicit schemes. The governing equations are coupled and then solved by a GMRES solver (generalized minimal residual method) using modified Gram-Schmidt orthogonalization. To show the validity of the method, we apply the method to the simulation of a rising droplet, a leaky dielectric drop and the coalescence of drops. The numerical solutions to the phase field model match well with existing solutions over a defined range of variables.
Abstract: In this numerical study, effects of using Al2O3-water
nanofluid on the rate of heat transfer have been investigated. Physical
model is a square enclosure with insulated top and bottom horizontal
walls, while the vertical walls are kept at different constant
temperatures. Two appropriate models are used to evaluate the
viscosity and thermal conductivity of nanofluid. The governing
stream-vorticity equations are solved using a second order central
finite difference scheme, coupled to the conservation of mass and
energy. The study has been carried out for the nanoparticle diameter
30, 60 and 90 nm and the solid volume fraction 0 to 0.04. Results are
presented by average Nusselt number and normalized Nusselt number
in different range of φ and D for mixed convection dominated
regime. It is found that different heat transfer rate is predicted when
the effect of nanoparticle diameter is taken into account.
Abstract: The formulated problem of optimization of the
technological process of water treatment for thermal power plants is
considered in this article. The problem is of multiparametric nature.
To optimize the process, namely, reduce the amount of waste water, a
new technology was developed to reuse such water. A mathematical
model of the technology of wastewater reuse was developed.
Optimization parameters were determined. The model consists of a
material balance equation, an equation describing the kinetics of ion
exchange for the non-equilibrium case and an equation for the ion
exchange isotherm. The material balance equation includes a
nonlinear term that depends on the kinetics of ion exchange. A direct
problem of calculating the impurity concentration at the outlet of the
water treatment plant was numerically solved. The direct problem
was approximated by an implicit point-to-point computation
difference scheme. The inverse problem was formulated as relates to
determination of the parameters of the mathematical model of the
water treatment plant operating in non-equilibrium conditions. The
formulated inverse problem was solved. Following the results of
calculation the time of start of the filter regeneration process was
determined, as well as the period of regeneration process and the
amount of regeneration and wash water. Multi-parameter
optimization of water treatment process for thermal power plants
allowed decreasing the amount of wastewater by 15%.
Abstract: In this paper, effects of using Alumina-water
nanofluid on the rate of heat transfer have been investigated
numerically. Physical model is a square enclosure with insulated top
and bottom horizontal walls, while the vertical walls are kept at
different constant temperatures. Two appropriate models are used to
evaluate the viscosity and thermal conductivity of nanofluid. The
governing stream-vorticity equations are solved using a second order
central finite difference scheme, coupled to the conservation of mass
and energy. The study has been carried out for the Richardson
number 0.1 to 10 and the solid volume fraction 0 to 0.04. Results are
presented by isotherms lines, average Nusselt number and normalized
Nusselt number in different range of φ and Ri for forced, combined
and natural convection dominated regime. It is found that higher heat
transfer rate is predicted when the effects of nanoparticle is taken into
account.
Abstract: The present analysis considers the steady stagnation point flow and heat transfer towards a permeable shrinking sheet in an upper-convected Maxwell (UCM) electrically conducting fluid, with a constant magnetic field applied in the transverse direction to flow and a local heat generation within the boundary layer, with a heat generation rate proportional to (T-T)p Using a similarity transformation, the governing system of partial differential equations is first transformed into a system of ordinary differential equations, which is then solved numerically using a finite-difference scheme known as the Keller-box method. Numerical results are obtained for the flow and thermal fields for various values of the stretching/shrinking parameter λ, the magnetic parameter M, the elastic parameter K, the Prandtl number Pr, the suction parameter s, the heat generation parameter Q, and the exponent p. The results indicate the existence of dual solutions for the shrinking sheet up to a critical value λc whose value depends on the value of M, K, and s. In the presence of internal heat absorption (Q
Abstract: A method for simulating flow around the solid bodies has been presented using hybrid meshfree and mesh-based schemes. The presented scheme optimizes the computational efficiency by combining the advantages of both meshfree and mesh-based methods. In this approach, a cloud of meshfree nodes has been used in the domain around the solid body. These meshfree nodes have the ability to efficiently adapt to complex geometrical shapes. In the rest of the domain, conventional Cartesian grid has been used beyond the meshfree cloud. Complex geometrical shapes can therefore be dealt efficiently by using meshfree nodal cloud and computational efficiency is maintained through the use of conventional mesh-based scheme on Cartesian grid in the larger part of the domain. Spatial discretization of meshfree nodes has been achieved through local radial basis functions in finite difference mode (RBF-FD). Conventional finite difference scheme has been used in the Cartesian ‘meshed’ domain. Accuracy tests of the hybrid scheme have been conducted to establish the order of accuracy. Numerical tests have been performed by simulating two dimensional steady and unsteady incompressible flows around cylindrical object. Steady flow cases have been run at Reynolds numbers of 10, 20 and 40 and unsteady flow problems have been studied at Reynolds numbers of 100 and 200. Flow Parameters including lift, drag, vortex shedding, and vorticity contours are calculated. Numerical results have been found to be in good agreement with computational and experimental results available in the literature.
Abstract: This paper presents the finite difference scheme and the numerical simulation of suspended string. The vibration solutions when the various external forces are taken into account are obtained and compared with the solutions without external force. In addition, we also investigate how the external forces and their powers and coefficients affect the amplitude of vibration.
Abstract: Discretization of spatial derivatives is an important
issue in meshfree methods especially when the derivative terms
contain non-linear coefficients. In this paper, various methods used
for discretization of second-order spatial derivatives are investigated
in the context of Smoothed Particle Hydrodynamics. Three popular
forms (i.e. "double summation", "second-order kernel derivation",
and "difference scheme") are studied using one-dimensional unsteady
heat conduction equation. To assess these schemes, transient response
to a step function initial condition is considered. Due to parabolic
nature of the heat equation, one can expect smooth and monotone
solutions. It is shown, however in this paper, that regardless of
the type of kernel function used and the size of smoothing radius,
the double summation discretization form leads to non-physical
oscillations which persist in the solution. Also, results show that when
a second-order kernel derivative is used, a high-order kernel function
shall be employed in such a way that the distance of inflection
point from origin in the kernel function be less than the nearest
particle distance. Otherwise, solutions may exhibit oscillations near
discontinuities unlike the "difference scheme" which unconditionally
produces monotone results.
Abstract: In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function ψ (x,y)and the function φ (x,y)as independent variables where for irrotational flow φ (x,y)can be recognized as the velocity potential function, for rotational flow φ (x,y)ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on the finite difference scheme on a uniform mesh is employed. The technique described is capable of tackling the so-called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct geometries. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists.
Abstract: A numerical analysis used to simulate the effects of wavy surfaces and thermal radiation on natural convection heat transfer boundary layer flow over an inclined wavy plate has been investigated. A simple coordinate transformation is employed to transform the complex wavy surface into a flat plate. The boundary layer equations and the boundary conditions are discretized by the finite difference scheme and solved numerically using the Gauss-Seidel algorithm with relaxation coefficient. Effects of the wavy geometry, the inclination angle of the wavy plate and the thermal radiation on the velocity profiles, temperature profiles and the local Nusselt number are presented and discussed in detail.
Abstract: Eight difference schemes and five limiters are applied to numerical computation of Riemann problem. The resolution of discontinuities of each scheme produced is compared. Numerical dissipation and its estimation are discussed. The result shows that the numerical dissipation of each scheme is vital to improve scheme-s accuracy and stability. MUSCL methodology is an effective approach to increase computational efficiency and resolution. Limiter should be selected appropriately by balancing compressive and diffusive performance.
Abstract: Several numerical schemes utilizing central difference
approximations have been developed to solve the Goursat problem.
However, in a recent years compact discretization methods which
leads to high-order finite difference schemes have been used since it
is capable of achieving better accuracy as well as preserving certain
features of the equation e.g. linearity. The basic idea of the new
scheme is to find the compact approximations to the derivative terms
by differentiating centrally the governing equations. Our primary
interest is to study the performance of the new scheme when applied
to two Goursat partial differential equations against the traditional
finite difference scheme.
Abstract: The paper deals with a mathematical model for fluid dynamic flows on road networks which is based on conservation laws. This nonlinear framework is based on the conservation of cars. We focus on traffic circle, which is a finite number of roads that meet at some junctions. The traffic circle with junctions having either one incoming and two outgoing or two incoming and one outgoing roads. We describe the numerical schemes with the particular boundary conditions used to produce approximated solutions of the problem.
Abstract: In the present analysis an unsteady laminar
forced convection water boundary layer flow is considered.
The fluid properties such as viscosity and Prandtl number are
taken as variables such that those are inversely proportional to
temperature. By using quasi-linearization technique the nonlinear
coupled partial differential equations are linearized and
the numerical solutions are obtained by using implicit finite
difference scheme with the appropriate selection of step sizes.
Non-similar solutions have been obtained from the starting
point of the stream-wise coordinate to the point where skin
friction value vanishes. The effect non-uniform mass transfer
along the surface of the cylinder through slot is studied on the
skin friction and heat transfer coefficients.
Abstract: A numerical solution of the initial boundary value
problem of the suspended string vibrating equation with the
particular nonlinear damping term based on the finite difference
scheme is presented in this paper. The investigation of how the
second and third power terms of the nonlinear term affect the
vibration characteristic. We compare the vibration amplitude as a
result of the third power nonlinear damping with the second power
obtained from previous report provided that the same initial shape
and initial velocities are assumed. The comparison results show that
the vibration amplitude is inversely proportional to the coefficient of
the damping term for the third power nonlinear damping case, while
the vibration amplitude is proportional to the coefficient of the
damping term in the second power nonlinear damping case.
Abstract: This paper examines the forced convection flow of
incompressible, electrically conducting viscous fluid past a sharp
wedge in the presence of heat generation or absorption with an
applied magnetic field. The system of partial differential equations
governing Falkner - Skan wedge flow and heat transfer is first
transformed into a system of ordinary differential equations using
similarity transformations which is later solved using an implicit
finite - difference scheme, along with quasilinearization technique.
Numerical computations are performed for air (Pr = 0.7) and
displayed graphically to illustrate the influence of pertinent physical
parameters on local skin friction and heat transfer coefficients and,
also on, velocity and temperature fields. It is observed that the
magnetic field increases both the coefficients of skin friction and heat
transfer. The effect of heat generation or absorption is found to be
very significant on heat transfer, but its effect on the skin friction is
negligible. Indeed, the occurrence of overshoot is noticed in the
temperature profiles during heat generation process, causing the
reversal in the direction of heat transfer.
Abstract: An efficient transient flow simulation for gas
pipelines and networks is presented. The proposed transient flow
simulation is based on the transfer function models and MATLABSimulink.
The equivalent transfer functions of the nonlinear
governing equations are derived for different types of the boundary
conditions. Next, a MATLAB-Simulink library is developed and
proposed considering any boundary condition type. To verify the
accuracy and the computational efficiency of the proposed
simulation, the results obtained are compared with those of the
conventional finite difference schemes (such as TVD, method of
lines, and other finite difference implicit and explicit schemes). The
effects of the flow inertia and the pipeline inclination are
incorporated in this simulation. It is shown that the proposed
simulation has a sufficient accuracy and it is computationally more
efficient than the other methods.
Abstract: This paper is a numerical investigation of a laminar
isothermal plane two dimensional wall jet. Special attention has been
paid to the effect of the inlet conditions at the nozzle exit on the
hydrodynamic and thermal characteristics of the flow. The
behaviour of various fluids evolving in both forced and mixed
convection regimes near a vertical plate plane is carried out. The
system of governing equations is solved with an implicit finite
difference scheme. For numerical stability we use a staggered non
uniform grid. The obtained results show that the effect of the Prandtl
number is significant in the plume region in which the jet flow is
governed by buoyant forces. Further for ascending X values, the
buoyancy forces become dominating, and a certain agreement
between the temperature profiles are observed, which shows that the
velocity profile has no longer influence on the wall temperature
evolution in this region. Fluids with low Prandtl number warm up
more importantly, because for such fluids the effect of heat diffusion
is higher.