Abstract: Flow movement in unsaturated soil can be expressed
by a partial differential equation, named Richards equation. The
objective of this study is the finding of an appropriate implicit
numerical solution for head based Richards equation. Some of the
well known finite difference schemes (fully implicit, Crank Nicolson
and Runge-Kutta) have been utilized in this study. In addition, the
effects of different approximations of moisture capacity function,
convergence criteria and time stepping methods were evaluated. Two
different infiltration problems were solved to investigate the
performance of different schemes. These problems include of vertical
water flow in a wet and very dry soils. The numerical solutions of
two problems were compared using four evaluation criteria and the
results of comparisons showed that fully implicit scheme is better
than the other schemes. In addition, utilizing of standard chord slope
method for approximation of moisture capacity function, automatic
time stepping method and difference between two successive
iterations as convergence criterion in the fully implicit scheme can
lead to better and more reliable results for simulation of fluid
movement in different unsaturated soils.
Abstract: The three-dimensional incompressible flow past a
rectangular open cavity is investigated, where the aspect ratio of the
cavity is considered as 4. The principle objective is to use large-eddy
simulation to resolve and control the large-scale structures, which are
largely responsible for flow oscillations in a cavity. The flow past an
open cavity is very common in aerospace applications and can be a
cause of acoustic source due to hydrodynamic instability of the shear
layer and its interactions with the downstream edge. The unsteady
Navier-stokes equations have been solved on a staggered mesh using
a symmetry-preserving central difference scheme. Synthetic jet has
been used as an active control to suppress the cavity oscillations in
wake mode for a Reynolds number of ReD = 3360. The effect of
synthetic jet has been studied by varying the jet amplitude and
frequency, which is placed at the upstream wall of the cavity. The
study indicates that there exits a frequency band, which is larger than
a critical value, is effective in attenuating cavity oscillations when
blowing ratio is more than 1.0.
Abstract: In this article, we aim to discuss the formulation of two explicit group iterative finite difference methods for time-dependent two dimensional Burger-s problem on a variable mesh. For the non-linear problems, the discretization leads to a non-linear system whose Jacobian is a tridiagonal matrix. We discuss the Newton-s explicit group iterative methods for a general Burger-s equation. The proposed explicit group methods are derived from the standard point and rotated point Crank-Nicolson finite difference schemes. Their computational complexity analysis is discussed. Numerical results are given to justify the feasibility of these two proposed iterative methods.
Abstract: Numerical investigation of the characteristics of an 80°
delta wing in combined force-pitch and free-roll is presented. The
implicit, upwind, flux-difference splitting, finite volume scheme and
the second-order-accurate finite difference scheme are employed to
solve the flow governing equations and Euler rigid-body dynamics
equations, respectively. The characteristics of the delta wing in
combined free-roll and large amplitude force-pitch is obtained
numerically and shows a well agreement with experimental data
qualitatively. The motion in combined force-pitch and free-roll
significantly reduces the lift force and transverse stabilities of the delta
wing, which is closely related to the flying safety. Investigations on
sensitive factors indicate that the roll-axis moment of inertia and the
structural damping have great influence on the frequency and
amplitude, respectively. Moreover, the turbulence model is considered
as an influencing factor in the investigation.
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: 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: 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: In this paper, we research the standard 13-point difference schemes for solving the biharmonic equation. Heuristic method is applied to judging the stability of multi-level difference schemes of the biharmonic equation. It is showed that the standard 13-point difference schemes are stable.
Abstract: The present paper considers the steady free convection
boundary layer flow of a viscoelastic fluid on solid sphere with
Newtonian heating. The boundary layer equations are an order higher
than those for the Newtonian (viscous) fluid and the adherence
boundary conditions are insufficient to determine the solution of
these equations completely. Thus, the augmentation an extra
boundary condition is needed to perform the numerical
computational. The governing boundary layer equations are first
transformed into non-dimensional form by using special
dimensionless group and then solved by using an implicit finite
difference scheme. The results are displayed graphically to illustrate
the influence of viscoelastic K and Prandtl Number Pr parameters on
skin friction, heat transfer, velocity profiles and temperature profiles.
Present results are compared with the published papers and are found
to concur very well.
Abstract: An analysis is performed to study the influence of nonuniform double slot suction on a steady laminar boundary layer flow over a rotating sphere when fluid properties such as viscosity and Prandtl number are inverse linear functions of temperature. Nonsimilar solutions have been obtained from the starting point of the streamwise co-ordinate to the exact point of separation. The difficulties arising at the starting point of the streamwise co-ordinate, at the edges of the slot and at the point of separation have been overcome by applying an implicit finite difference scheme in combination with the quasi-linearization technique and an appropriate selection of the finer step sizes along the stream-wise direction. The present investigation shows that the point of ordinary separation can be delayed by nonuniform double slot suction if the mass transfer rate is increased and also if the slots are positioned further downstream. In addition, the investigation reveals that double slot suction is found to be more effective compared to a single slot suction in delaying ordinary separation. As rotation parameter increase the point of separation moves upstream direction.