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 convection term are discretized using central differences,
stability problems occur when the grid spacing is chosen too coarse.
In this paper, we introduce the splitting upwind schemes for avoiding
stability problems and prove that it is consistent to the upwind scheme
with same accuracy. The splitting upwind schemes was adopted
to split the wave spectral action balance equation into four onedimensional
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-processor computer.
Abstract: The stability characteristics of water lubricated journal bearings having three axial grooves are obtained theoretically. In this lubricant (water) is fed under pressure from one end of the bearing, through the 3-axial grooves (groove angles may vary). These bearings can use the process fluid as the lubricant, as in the case of feed water pumps. The Reynolds equation is solved numerically by the finite difference method satisfying the boundary conditions. The stiffness and damping coefficient for various bearing number and eccentricity ratios, assuming linear pressure drop along the groove, shows that smaller groove angles better results.
Abstract: In this paper, the difference between the Alternating
Direction Method (ADM) and the Non-Splitting Method (NSM) is
investigated, while both methods applied to the simulations for 2-D
multimaterial radiation diffusion issues. Although the ADM have the
same accuracy orders with the NSM on the uniform meshes, the
accuracy of ADM will decrease on the distorted meshes or the
boundary of domain. Numerical experiments are carried out to
confirm the theoretical predication.
Abstract: A simple microstructure optical fiber design based on an octagonal cladding structure is presented for simultaneously controlling dispersion and leakage properties. The finite difference method with anisotropic perfectly matched boundary layer is used to investigate the guiding properties. It is demonstrated that octagonal photonic crystal fibers with four rings can assume negative ultra-flattened dispersion of -19 + 0.23 ps/nm/km in the wavelength range of 1.275 μm to 1.68 μm, nearly zero ultra-flattened dispersion of 0 ± 0.40 ps/nm/km in a 1.38 to 1.64 μm, and low confinement losses less than 10-3 dB/km in the entire band of interest.
Abstract: An important task in solving second order linear ordinary differential equations by the finite difference is to choose a suitable stepsize h. In this paper, by using the stochastic arithmetic, the CESTAC method and the CADNA library we present a procedure to estimate the optimal stepsize hopt, the stepsize which minimizes the global error consisting of truncation and round-off error.
Abstract: Smoke discharging is a main reason of air pollution
problem from industrial plants. The obstacle of a building has an
affect with the air pollutant discharge. In this research, a mathematical
model of the smoke dispersion from two sources and one source with
a structural obstacle is considered. The governing equation of the
model is an isothermal mass transfer model in a viscous fluid. The
finite element method is used to approximate the solutions of the
model. The triangular linear elements have been used for discretising
the domain, and time integration has been carried out by semi-implicit
finite difference method. The simulations of smoke dispersion in
cases of one chimney and two chimneys are presented. The maximum
calculated smoke concentration of both cases are compared. It is then
used to make the decision for smoke discharging and air pollutant
control problems on industrial area.
Abstract: In this article, our objective is the analysis of the resolution of non-linear differential systems by combining Newton and Continuation (N-C) method. The iterative numerical methods converge where the initial condition is chosen close to the exact solution. The question of choosing the initial condition is answered by N-C method.
Abstract: For fire safety purposes, the fire resistance and the
structural behavior of reinforced concrete members are assessed to
satisfy specific fire performance criteria. The available prescribed
provisions are based on standard fire load. Under various fire
scenarios, engineers are in need of both heat transfer analysis and
structural analysis. For heat transfer analysis, the study proposed a
modified finite difference method to evaluate the temperature profile
within a cross section. The research conducted is limited to concrete
sections exposed to a fire on their one side. The method is based on
the energy conservation principle and a pre-determined power
function of the temperature profile. The power value of 2.7 is found
to be a suitable value for concrete sections. The temperature profiles
of the proposed method are only slightly deviate from those of the
experiment, the FEM and the FDM for various fire loads such as
ASTM E 119, ASTM 1529, BS EN 1991-1-2 and 550 oC. The
proposed method is useful to avoid incontinence of the large matrix
system of the typical finite difference method to solve the
temperature profile. Furthermore, design engineers can simply apply
the proposed method in regular spreadsheet software.
Abstract: Few studies have been conducted on polymeric strip
and the behavior of soil retaining walls. This paper will present the
effect of frequency on the dynamic behavior of reinforced soil
retaining walls with polymeric strips. The frequency content
describes how the amplitude of a ground motion is distributed among
different frequencies. Since the frequency content of an earthquake
motion will strongly influence the effects of that motion, the
characterization of the motion cannot be completed without the
consideration of its frequency content. The maximum axial force of
reinforcements and horizontal displacement of the reinforced walls
are focused in this research. To clarify the dynamic behavior of
reinforced soil retaining walls with polymeric strips, a numerical
modeling using Finite Difference Method is benefited. As the results
indicate, the frequency of input base acceleration has an important
effect on the behavior of these structures. Because of resonant in the
system, where the frequency of the input dynamic load is equal to the
natural frequency of the system, the maximum horizontal
displacement and the maximum axial forces in polymeric strips is
occurred. Moreover, they were to increase the structure flexibility
because of the main advantages of polymeric strips; i.e. being simple
method of construction, having a homogeneous behavior with soils,
and possessing long durability, which are of great importance in
dynamic analysis.
Abstract: The aim of this paper is to study the internal
stabilization of the Bernoulli-Euler equation numerically. For this,
we consider a square plate subjected to a feedback/damping force
distributed only in a subdomain. An algorithm for obtaining an
approximate solution to this problem was proposed and implemented.
The numerical method used was the Finite Difference Method.
Numerical simulations were performed and showed the behavior of
the solution, confirming the theoretical results that have already been
proved in the literature. In addition, we studied the validation of the
numerical scheme proposed, followed by an analysis of the numerical
error; and we conducted a study on the decay of the energy associated.
Abstract: The present study investigates numerically the
phenomenon of vortex-shedding and its suppression in twodimensional
mixed convective flow past a square cylinder under the
joint influence of buoyancy and free-stream orientation with respect
to gravity. The numerical experiments have been conducted at a
fixed Reynolds number (Re) of 100 and Prandtl number (Pr) of 0.71,
while Richardson number (Ri) is varied from 0 to 1.6 and freestream
orientation, α, is kept in the range 0o≤ α ≤ 90o, with 0o
corresponding to an upward flow and 90o representing a cross-flow
scenario, respectively. The continuity, momentum and energy
equations, subject to Boussinesq approximation, are discretized using
a finite difference method and are solved by a semi-explicit pressure
correction scheme. The critical Richardson number, leading to the
suppression of the vortex-shedding (Ric), is estimated by using
Stuart-Landau theory at various free-stream orientations and the
neutral curve is obtained in the Ri-α plane. The neutral curve
exhibits an interesting non-monotonic behavior with Ric first
increasing with increasing values of α upto 45o and then decreasing
till 70o. Beyond 70o, the neutral curve again exhibits a sharp
increasing asymptotic trend with Ric approaching very large values
as α approaches 90o. The suppression of vortex shedding is not
observed at α = 90o (cross-flow). In the unsteady flow regime, the
Strouhal number (St) increases with the increase in Richardson
number.
Abstract: Artificial atoms are growing fields of interest due to their physical and optoelectronicapplications. The absorption spectra of the proposed artificial atom inpresence of Tera-Hertz field is investigated theoretically. We use the non-perturbativeFloquet theory and finite difference method to study the electronic structure of ArtificialAtom. The effect of static electric field on the energy levels of artificial atom is studied.The effect of orientation of static electric field on energy levels and diploe matrix elementsis also highlighted.
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: Subgrade moisture content varies with environmental and soil conditions and has significant influence on pavement performance. Therefore, it is important to establish realistic estimates of expected subgrade moisture contents to account for the effects of this variable on predicted pavement performance during the design stage properly. The initial boundary soil suction profile for a given pavement is a critical factor in determining expected moisture variations in the subgrade for given pavement and climatic and soil conditions. Several numerical models have been developed for predicting water and solute transport in saturated and unsaturated subgrade soils. Soil hydraulic properties are required for quantitatively describing water and chemical transport processes in soils by the numerical models. The required hydraulic properties are hydraulic conductivity, water diffusivity, and specific water capacity. The objective of this paper was to determine isothermal moisture profiles in a soil fill and predict the soil moisture movement above the ground water table using a simple one-dimensional finite difference model.
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: This paper adopted the hybrid differential transform approach for studying heat transfer problems in a gold/chromium thin film with an ultra-short-pulsed laser beam projecting on the gold side. The physical system, formulated based on the hyperbolic two-step heat transfer model, covers three characteristics: (i) coupling effects between the electron/lattice systems, (ii) thermal wave propagation in metals, and (iii) radiation effects along the interface. The differential transform method is used to transfer the governing equations in the time domain into the spectrum equations, which is further discretized in the space domain by the finite difference method. The results, obtained through a recursive process, show that the electron temperature in the gold film can rise up to several thousand degrees before its electron/lattice systems reach equilibrium at only several hundred degrees. The electron and lattice temperatures in the chromium film are much lower than those in the gold film.
Abstract: Water pollution assessment problems arise frequently
in environmental science. In this research, a finite difference method
for solving the one-dimensional steady convection-diffusion equation
with variable coefficients is proposed; it is then used to optimize
water treatment costs.
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: 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.