Abstract: In the present study, a numerical analysis is carried
out to investigate unsteady MHD (magneto-hydrodynamic) flow and
heat transfer of a non-Newtonian second grade viscoelastic fluid
over an oscillatory stretching sheet. The flow is induced due to an
infinite elastic sheet which is stretched oscillatory (back and forth) in
its own plane. Effect of viscous dissipation and joule heating are
taken into account. The non-linear differential equations governing
the problem are transformed into system of non-dimensional
differential equations using similarity transformations. A newly
developed meshfree numerical technique Element free Galerkin
method (EFGM) is employed to solve the coupled non linear
differential equations. The results illustrating the effect of various
parameters like viscoelastic parameter, Hartman number, relative
frequency amplitude of the oscillatory sheet to the stretching rate and
Eckert number on velocity and temperature field are reported in
terms of graphs and tables. The present model finds its application in
polymer extrusion, drawing of plastic films and wires, glass, fiber
and paper production etc.
Abstract: In this paper we study some numerical methods to solve a model one-dimensional convection–diffusion equation. The semi-discretisation of the space variable results into a system of ordinary differential equations and the solution of the latter involves the evaluation of a matrix exponent. Since the calculation of this term is computationally expensive, we study some methods based on Krylov subspace and on Restrictive Taylor series approximation respectively. We also consider the Chebyshev Pseudospectral collocation method to do the spatial discretisation and we present the numerical solution obtained by these methods.
Abstract: One of the potential and effective ways of
storing thermal energy in buildings is the integration of brick with phase change materials (PCMs). This paper presents a two-dimensional model for simulating and analyzing of PCM
in order to minimize energy consumption in the buildings. The numerical approach has been used with the real weather data of a selected city of Iran (Tehran). Two kinds of brick integrated PCM are investigated and compared base on
outdoor weather conditions and the amount of energy
consumption. The results show a significant reduction in
maximum entering heat flux to building about 32.8%
depending on PCM quantity. The results are analyzed by
various temperature contour plots. The contour plots
illustrated the time dependent mechanism of entering heat flux for a brick integrated with PCM. Further analysis is developed to investigate the effect of PCM location on the inlet heat flux. The results demonstrated that to achieve maximum performance of PCM it is better to locate PCM near the outdoor.
Abstract: Meshing is the process of discretizing problem
domain into many sub domains before the numerical calculation can
be performed. One of the most popular meshes among many types of meshes is tetrahedral mesh, due to their flexibility to fit into almost
any domain shape. In both 2D and 3D domains, triangular and tetrahedral meshes can be generated by using Delaunay triangulation.
The quality of mesh is an important factor in performing any Computational Fluid Dynamics (CFD) simulations as the results is
highly affected by the mesh quality. Many efforts had been done in
order to improve the quality of the mesh. The paper describes a mesh
generation routine which has been developed capable of generating
high quality tetrahedral cells in arbitrary complex geometry. A few
test cases in CFD problems are used for testing the mesh generator.
The result of the mesh is compared with the one generated by a
commercial software. The results show that no sliver exists for the
meshes generated, and the overall quality is acceptable since the percentage of the bad tetrahedral is relatively small. The boundary
recovery was also successfully done where all the missing faces are
rebuilt.
Abstract: This paper presents a numerical study on surface heat
transfer characteristics of laminar air flows in parallel-plate dimpled
channels. The two-dimensional numerical model is provided by
commercial code FLUENT and the results are obtained for channels
with symmetrically opposing hemi-cylindrical cavities onto both
walls for Reynolds number ranging from 1000 to 2500. The influence
of variations in relative depth of dimples (the ratio of cavity depth to
the cavity curvature diameter), the number of them and the thermophysical
properties of channel walls on heat transfer enhancement is
studied. The results are evident for existence of an optimum value for
the relative depth of dimples in which the largest wall heat flux and
average Nusselt number can be achieved. In addition, the results of
conjugation simulation indicate that the overall influence of the ratio
of wall thermal conductivity to the one of the fluid on heat transfer
rate is not much significant and can be ignored.
Abstract: A numerical investigation of the effects of nanosecond
barrier discharge on the stability of a two-dimensional free shear layer
is performed. The computations are carried out using a compressible
Navier-Stokes algorithm coupled with a thermodynamic model of the
discharge. The results show that significant increases in the shear
layer-s momentum thickness and Reynolds stresses occur due to
actuation. Dependence on both frequency and amplitude of actuation
are considered, and a comparison is made of the computed growth
rates with those predicted by linear stability theory. Amplitude and
frequency ranges for the efficient promotion of shear-layer instabilities
are identified.
Abstract: This article presents a numerical study of the doublediffusive
mixed convection in a vertical channel filled with porous
medium by using non-equilibrium model. The flow is assumed
fully developed, uni-directional and steady state. The controlling
parameters are thermal Rayleigh number (RaT ), Darcy number (Da),
Forchheimer number (F), buoyancy ratio (N), inter phase heat transfer
coefficient (H), and porosity scaled thermal conductivity ratio
(γ). The Brinkman-extended non-Darcy model is considered. The
governing equations are solved by spectral collocation method. The
main emphasize is given on flow profiles as well as heat and solute
transfer rates, when two diffusive components in terms of buoyancy
ratio are in favor (against) of each other and solid matrix and fluid
are thermally non-equilibrium. The results show that, for aiding flow
(RaT = 1000), the heat transfer rate of fluid (Nuf ) increases upto a
certain value of H, beyond that decreases smoothly and converges
to a constant, whereas in case of opposing flow (RaT = -1000),
the result is same for N = 0 and 1. The variation of Nuf in (N,
Nuf )-plane shows sinusoidal pattern for RaT = -1000. For both cases
(aiding and opposing) the flow destabilize on increasing N by inviting
point of inflection or flow separation on the velocity profile. Overall,
the buoyancy force have significant impact on the non-Darcy mixed
convection under LTNE conditions.
Abstract: This paper focuses on the probabilistic numerical
solution of the problems in biomechanics and mining. Applications of
Simulation-Based Reliability Assessment (SBRA) Method are
presented in the solution of designing of the external fixators applied
in traumatology and orthopaedics (these fixators can be applied for
the treatment of open and unstable fractures etc.) and in the solution
of a hard rock (ore) disintegration process (i.e. the bit moves into the
ore and subsequently disintegrates it, the results are compared with
experiments, new design of excavation tool is proposed.
Abstract: A multiple-option analytical model for the evaluation of the energy performance and distribution of aerodynamic forces acting on a vertical-axis Darrieus wind turbine depending on both rotor architecture and operating conditions is presented. For this purpose, a numerical algorithm, capable of generating the desired rotor conformation depending on design geometric parameters, is coupled to a Single/Double-Disk Multiple-Streamtube Blade Element – Momentum code. Both single and double-disk configurations are analyzed and model predictions are compared to literature experimental data in order to test the capability of the code for predicting rotor performance. Effective airfoil characteristics based on local blade Reynolds number are obtained through interpolation of literature low-Reynolds airfoil databases. Some corrections are introduced inside the original model with the aim of simulating also the effects of blade dynamic stall, rotor streamtube expansion and blade finite aspect ratio, for which a new empirical relationship to better fit the experimental data is proposed. In order to predict also open field rotor operation, a freestream wind shear profile is implemented, reproducing the effect of atmospheric boundary layer.
Abstract: An optimal control problem for a mathematical model of efficiency of antiviral therapy in hepatitis B virus infections is considered. The aim of the study is to control the new viral production, block the new infection cells and maintain the number of uninfected cells in the given range. The optimal controls represent the efficiency of antiviral therapy in inhibiting viral production and preventing new infections. Defining the cost functional, the optimal control problem is converted into the constrained optimization problem and the first order optimality system is derived. For the numerical simulation, we propose the steepest descent algorithm based on the adjoint variable method. A computer program in MATLAB is developed for the numerical simulations.
Abstract: This paper deals with a periodic-review substitutable
inventory system for a finite and an infinite number of periods. Here
an upward substitution structure, a substitution of a more costly item
by a less costly one, is assumed, with two products. At the beginning
of each period, a stochastic demand comes for the first item only,
which is quality-wise better and hence costlier. Whenever an arriving
demand finds zero inventory of this product, a fraction of unsatisfied
customers goes for its substitutable second item. An optimal ordering
policy has been derived for each period. The results are illustrated
with numerical examples. A sensitivity analysis has been done to
examine how sensitive the optimal solution and the maximum profit
are to the values of the discount factor, when there is a large number
of periods.
Abstract: Based on the classical algorithm LSQR for solving (unconstrained) LS problem, an iterative method is proposed for the least-squares like-minimum-norm symmetric solution of AXB+CYD=E. As the application of this algorithm, an iterative method for the least-squares like-minimum-norm biymmetric solution of AXB=E is also obtained. Numerical results are reported that show the efficiency of the proposed methods.
Abstract: Appropriate ventilation in a classroom is helpful for
enhancing air exchange rate and student concentration. This study
focuses on the effects of fenestration in a four-story school building by
performing numerical simulation of a building when considering
indoor and outdoor environments simultaneously. The wind profile
function embedded in PHOENICS code was set as the inlet boundary
condition in a suburban environment. Sixteen fenestration
combinations were compared in a classroom containing thirty seats.
This study evaluates mean age of air (AGE) and airflow pattern of a
classroom on different floors. Considering both wind profile and
fenestration effects, the airflow on higher floors is channeled toward
the area near ceiling in a room and causes older mean age of air in the
breathing zone. The results in this study serve as a useful guide for
enhancing natural ventilation in a typical school building.
Abstract: This paper proposes improved delay-dependent stability conditions of the linear time-delay systems of neutral type. The proposed methods employ a suitable Lyapunov-Krasovskii’s functional and a new form of the augmented system. New delay-dependent stability criteria for the systems are established in terms of Linear matrix inequalities (LMIs) which can be easily solved by various effective optimization algorithms. Numerical examples showed that the proposed method is effective and can provide less conservative results.
Abstract: This paper proposes a power-controlled scheduling scheme for devices using a directional antenna in smart home. In the case of the home network using directional antenna, devices can concurrently transmit data in the same frequency band. Accordingly, the throughput increases compared to that of devices using omni-directional antenna in proportional to the number of concurrent transmissions. Also, the number of concurrent transmissions depends on the beamwidth of antenna, the number of devices operating in the network , transmission power, interference and so on. In particular, the less transmission power is used, the more concurrent transmissions occur due to small transmission range. In this paper, we considered sub-optimal scheduling scheme for throughput maximization and power consumption minimization. In the scheme, each device is equipped with a directional antenna. Various beamwidths, path loss components, and antenna radiation efficiencies are considered. Numerical results show that the proposed schemes outperform the scheduling scheme using directional antennas without power control.
Abstract: It is essential to have a uniform and calm flow field
for a settling tank to have high performance. In general, the
recirculation zones always occurred in sedimentation tanks. The
presence of these regions may have different effects. The nonuniformity
of the velocity field, the short-circuiting at the surface and
the motion of the jet at the bed of the tank that occurs because of the
recirculation in the sedimentation layer, are affected by the geometry
of the tank. There are some ways to decrease the size of these dead
zones, which would increase the performance. One of the ways is to
use a suitable baffle configuration. In this study, the presence of
baffle with different position has been investigated by a finite volume
method, with VOF (Volume of Fluid) model. Besides, the k-ε
turbulence model is used in the numerical calculations. The results
indicate that the best position of the baffle is obtained when the
volume of the recirculation region is minimized or is divided to
smaller part and the flow field trend to be uniform in the settling
zone.
Abstract: The performance of time-reversal MUSIC algorithm will be dramatically degrades in presence of strong noise and multiple scattering (i.e. when scatterers are close to each other). This is due to error in determining the number of scatterers. The present paper provides a new approach to alleviate such a problem using an information theoretic criterion referred as minimum description length (MDL). The merits of the novel approach are confirmed by the numerical examples. The results indicate the time-reversal MUSIC yields accurate estimate of the target locations with considerable noise and multiple scattering in the received signals.
Abstract: In this paper, the Gaussian type quadrature rules for fuzzy functions are discussed. The errors representation and convergence theorems are given. Moreover, four kinds of Gaussian type quadrature rules with error terms for approximate of fuzzy integrals are presented. The present paper complements the theoretical results of the paper by T. Allahviranloo and M. Otadi [T. Allahviranloo, M. Otadi, Gaussian quadratures for approximate of fuzzy integrals, Applied Mathematics and Computation 170 (2005) 874-885]. The obtained results are illustrated by solving some numerical examples.
Abstract: This work proposes a recursive weighted ELS
algorithm for system identification by applying numerically robust
orthogonal Householder transformations. The properties of the
proposed algorithm show it obtains acceptable results in a noisy
environment: fast convergence and asymptotically unbiased
estimates. Comparative analysis with others robust methods well
known from literature are also presented.
Abstract: In this work we study analytically and numerically the
performance of the mean heave motion of an OWC coupled with the
governing equation of the spreading ocean waves due to the wide
variation in an open parabolic channel with constant depth. This
paper considers that the ocean wave propagation is under the
assumption of a shallow flow condition. In order to verify the effect
of the waves in the OWC firstly we establish the analytical model in
a non-dimensional form based on the energy equation. The proposed
wave-power system has to aims: one is to perturb the ocean waves as
a consequence of the channel shape in order to concentrate the
maximum ocean wave amplitude in the neighborhood of the OWC
and the second is to determine the pressure and volume oscillation of
air inside the compression chamber.