Abstract: The aim of this paper is to investigate the
performance of the developed two point block method designed for
two processors for solving directly non stiff large systems of higher
order ordinary differential equations (ODEs). The method calculates
the numerical solution at two points simultaneously and produces
two new equally spaced solution values within a block and it is
possible to assign the computational tasks at each time step to a
single processor. The algorithm of the method was developed in C
language and the parallel computation was done on a parallel shared
memory environment. Numerical results are given to compare the
efficiency of the developed method to the sequential timing. For
large problems, the parallel implementation produced 1.95 speed-up
and 98% efficiency for the two processors.
Abstract: In this paper, a new approach is introduced to solve
Blasius equation using parameter identification of a nonlinear
function which is used as approximation function. Bees Algorithm
(BA) is applied in order to find the adjustable parameters of
approximation function regarding minimizing a fitness function
including these parameters (i.e. adjustable parameters). These
parameters are determined how the approximation function has to
satisfy the boundary conditions. In order to demonstrate the
presented method, the obtained results are compared with another
numerical method. Present method can be easily extended to solve a
wide range of problems.
Abstract: Einstein vacuum equations, that is a system of nonlinear
partial differential equations (PDEs) are derived from Weyl metric
by using relation between Einstein tensor and metric tensor. The
symmetries of Einstein vacuum equations for static axisymmetric
gravitational fields are obtained using the Lie classical method. We
have examined the optimal system of vector fields which is further
used to reduce nonlinear PDE to nonlinear ordinary differential
equation (ODE). Some exact solutions of Einstein vacuum equations
in general relativity are also obtained.
Abstract: The governing differential equations of laminated
plate utilizing trigonometric shear deformation theory are derived
using energy approach. The governing differential equations
discretized by different radial basis functions are used to predict the
free vibration behavior of symmetric laminated composite plates.
Effect of orthotropy and span to thickness ratio on frequency
parameter of simply supported laminated plate is presented.
Numerical results show the accuracy and good convergence of radial
basis functions.
Abstract: This paper presents investigation effects of a sharp edged gust on aeroelastic behavior and time-domain response of a typical section model using Jones approximate aerodynamics for pure plunging motion. Flutter analysis has been done by using p and p-k methods developed for presented finite-state aerodynamic model for a typical section model (airfoil). Introduction of gust analysis as a linear set of ordinary differential equations in a simplified procedure has been carried out by using transformation into an eigenvalue problem.
Abstract: This paper presents a new function expansion method for finding traveling wave solution of a non-linear equation and calls it the (G'/G)-expansion method. The shallow water wave equation is reduced to a non linear ordinary differential equation by using a simple transformation. As a result the traveling wave solutions of shallow water wave equation are expressed in three forms: hyperbolic solutions, trigonometric solutions and rational solutions.
Abstract: This work investigated the steady state and dynamic
simulation of a fixed bed industrial naphtha reforming reactors. The
performance of the reactor was investigated using a heterogeneous
model. For process simulation, the differential equations are solved
using the 4th order Runge-Kutta method .The models were validated
against measured process data of an existing naphtha reforming plant.
The results of simulation in terms of components yields and
temperature of the outlet were in good agreement with empirical data.
The simple model displays a useful tool for dynamic simulation,
optimization and control of naphtha reforming.
Abstract: The heart tissue is an excitable media. A Cellular
Automata is a type of model that can be used to model cardiac action
potential propagation. One of the advantages of this approach against
the methods based on differential equations is its high speed in large
scale simulations. Recent cellular automata models are not able to
avoid flat edges in the result patterns or have large neighborhoods. In
this paper, we present a new model to eliminate flat edges by
minimum number of neighbors.
Abstract: The aim of this paper is to investigate twodimensional unsteady flow of a viscous incompressible fluid about stagnation point on permeable stretching sheet in presence of time dependent free stream velocity. Fluid is considered in the influence of transverse magnetic field in the presence of radiation effect. Rosseland approximation is use to model the radiative heat transfer. Using time-dependent stream function, partial differential equations corresponding to the momentum and energy equations are converted into non-linear ordinary differential equations. Numerical solutions of these equations are obtained by using Runge-Kutta Fehlberg method with the help of Newton-Raphson shooting technique. In the present work the effect of unsteadiness parameter, magnetic field parameter, radiation parameter, stretching parameter and the Prandtl number on flow and heat transfer characteristics have been discussed. Skin-friction coefficient and Nusselt number at the sheet are computed and discussed. The results reported in the paper are in good agreement with published work in literature by other researchers.
Abstract: Aim. We have introduced the notion of order to multinormed spaces and countable union spaces and their duals. The topology of bounded convergence is assigned to the dual spaces. The aim of this paper is to develop the theory of ordered topological linear
spaces La,b, L(w, z), the dual spaces of ordered multinormed spaces
La,b, ordered countable union spaces L(w, z), with the topology of bounded convergence assigned to the dual spaces. We apply Laplace transformation to the ordered linear space of Laplace transformable
generalized functions. We ultimately aim at finding solutions to nonhomogeneous
nth order linear differential equations with constant
coefficients in terms of generalized functions and comparing different
solutions evolved out of different initial conditions.
Method. The above aim is achieved by
• Defining the spaces La,b, L(w, z).
• Assigning an order relation on these spaces by identifying a
positive cone on them and studying the properties of the cone.
• Defining an order relation on the dual spaces La,b, L(w, z) of La,b, L(w, z) and assigning a topology to these dual spaces which makes the order dual and the topological dual the same. • Defining the adjoint of a continuous map on these spaces
and studying its behaviour when the topology of bounded
convergence is assigned to the dual spaces.
• Applying the two-sided Laplace Transformation on the ordered
linear space of generalized functions W and studying some
properties of the transformation which are used in solving
differential equations.
Result. The above techniques are applied to solve non-homogeneous
n-th order linear differential equations with constant coefficients in
terms of generalized functions and to compare different solutions of the differential equation.
Abstract: The steady mixed convection boundary layer flow from
a vertical cone in a porous medium filled with a nanofluid is
numerically investigated using different types of nanoparticles as Cu
(copper), Al2O3 (alumina) and TiO2 (titania). The boundary value
problem is solved by using the shooting technique by reducing it
into an ordinary differential equation. Results of interest for the local
Nusselt number with various values of the constant mixed convection
parameter and nanoparticle volume fraction parameter are evaluated.
It is found that dual solutions exist for a certain range of mixed
convection parameter.
Abstract: Natural convection heat transfer from a heated
horizontal semi-circular cylinder (flat surface upward) has been
investigated for the following ranges of conditions; Grashof number,
and Prandtl number. The governing partial differential equations
(continuity, Navier-Stokes and energy equations) have been solved
numerically using a finite volume formulation. In addition, the role of
the type of the thermal boundary condition imposed at cylinder
surface, namely, constant wall temperature (CWT) and constant heat
flux (CHF) are explored. Natural convection heat transfer from a
heated horizontal semi-circular cylinder (flat surface upward) has
been investigated for the following ranges of conditions; Grashof
number, and Prandtl number, . The governing partial differential
equations (continuity, Navier-Stokes and energy equations) have
been solved numerically using a finite volume formulation. In
addition, the role of the type of the thermal boundary condition
imposed at cylinder surface, namely, constant wall temperature
(CWT) and constant heat flux (CHF) are explored. The resulting flow
and temperature fields are visualized in terms of the streamline and
isotherm patterns in the proximity of the cylinder. The flow remains
attached to the cylinder surface over the range of conditions spanned
here except that for and ; at these conditions, a separated flow
region is observed when the condition of the constant wall
temperature is prescribed on the surface of the cylinder. The heat
transfer characteristics are analyzed in terms of the local and average
Nusselt numbers. The maximum value of the local Nusselt number
always occurs at the corner points whereas it is found to be minimum
at the rear stagnation point on the flat surface. Overall, the average
Nusselt number increases with Grashof number and/ or Prandtl
number in accordance with the scaling considerations. The numerical
results are used to develop simple correlations as functions of
Grashof and Prandtl number thereby enabling the interpolation of the
present numerical results for the intermediate values of the Prandtl or
Grashof numbers for both thermal boundary conditions.
Abstract: In this paper, we study a new modified Novikov equation for its classical and nonclassical symmetries and use the symmetries to reduce it to a nonlinear ordinary differential equation (ODE). With the aid of solutions of the nonlinear ODE by using the modified (G/G)-expansion method proposed recently, multiple exact traveling wave solutions are obtained and the traveling wave solutions are expressed by the hyperbolic functions, trigonometric functions and rational functions.
Abstract: The group invariant solution for Prandtl-s boundary layer equations for an incompressible fluid governing the flow in radial free, wall and liquid jets having finite fluid velocity at the orifice are investigated. For each jet a symmetry is associated with the conserved vector that was used to derive the conserved quantity for the jet elsewhere. This symmetry is then used to construct the group invariant solution for the third-order partial differential equation for the stream function. The general form of the group invariant solution for radial jet flows is derived. The general form of group invariant solution and the general form of the similarity solution which was obtained elsewhere are the same.
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: For many industrial applications plate heat
exchangers are demonstrating a large superiority over the
other types of heat exchangers. The efficiency of such a
device depends on numerous factors the effect of which needs
to be analysed and accurately evaluated.
In this paper we present a theoretical analysis of a cocurrent
plate heat exchanger and the results of its numerical
simulation.
Knowing the hot and the cold fluid streams inlet temperatures,
the respective heat capacities mCp
and the value of the
overall heat transfer coefficient, a 1-D mathematical model
based on the steady flow energy balance for a differential
length of the device is developed resulting in a set of N first
order differential equations with boundary conditions where N
is the number of channels.For specific heat exchanger
geometry and operational parameters, the problem is
numerically solved using the shooting method.
The simulation allows the prediction of the temperature
map in the heat exchanger and hence, the evaluation of its
performances. A parametric analysis is performed to evaluate
the influence of the R-parameter on the e-NTU values. For
practical purposes effectiveness-NTU graphs are elaborated
for specific heat exchanger geometry and different operating
conditions.
Abstract: The objective of this paper is to analyse the
application of the Half-Sweep Gauss-Seidel (HSGS) method by using
the Half-sweep approximation equation based on central difference
(CD) and repeated trapezoidal (RT) formulas to solve linear fredholm
integro-differential equations of first order. The formulation and
implementation of the Full-Sweep Gauss-Seidel (FSGS) and Half-
Sweep Gauss-Seidel (HSGS) methods are also presented. The HSGS
method has been shown to rapid compared to the FSGS methods.
Some numerical tests were illustrated to show that the HSGS method
is superior to the FSGS method.
Abstract: An evolutionary computing technique for solving initial value problems in Ordinary Differential Equations is proposed in this paper. Neural network is used as a universal approximator while the adaptive parameters of neural networks are optimized by genetic algorithm. The solution is achieved on the continuous grid of time instead of discrete as in other numerical techniques. The comparison is carried out with classical numerical techniques and the solution is found with a uniform accuracy of MSE ≈ 10-9 .
Abstract: This paper presents an analytical method to solve
governing consolidation parabolic partial differential equation (PDE)
for inelastic porous Medium (soil) with consideration of variation of
equation coefficient under cyclic loading. Since under cyclic loads,
soil skeleton parameters change, this would introduce variable
coefficient of parabolic PDE. Classical theory would not rationalize
consolidation phenomenon in such condition. In this research, a
method based on time space mapping to a virtual time space along
with superimposing rule is employed to solve consolidation of
inelastic soils in cyclic condition. Changes of consolidation
coefficient applied in solution by modification of loading and
unloading duration by introducing virtual time. Mapping function is
calculated based on consolidation partial differential equation results.
Based on superimposing rule a set of continuous static loads in
specified times used instead of cyclic load. A set of laboratory
consolidation tests under cyclic load along with numerical
calculations were performed in order to verify the presented method.
Numerical solution and laboratory tests results showed accuracy of
presented method.
Abstract: An optimal mean-square fusion formulas with scalar
and matrix weights are presented. The relationship between them is
established. The fusion formulas are compared on the continuous-time
filtering problem. The basic differential equation for cross-covariance
of the local errors being the key quantity for distributed fusion is
derived. It is shown that the fusion filters are effective for multi-sensor
systems containing different types of sensors. An example
demonstrating the reasonable good accuracy of the proposed filters is
given.