Abstract: Railways in India remain primary mode of transport
having one of the largest networks in the world and catering to
billions of transits yearly. Catastrophic economic damage and loss
to life is encountered over the past few decades due to fire to
locomotives. Study of fire dynamics and fire propagation plays an
important role in evacuation planning and reducing losses. Simulation
based study of propagation of fire and soot inside an air conditioned
coach of Indian locomotive is done in this paper. Finite difference
based solver, Fire Dynamic Simulator (FDS) version 6 has been
used for analysis. A single air conditioned 3 tier coupe closed to
ambient surroundings by glass windows having occupancy for 8
people is the basic unit of the domain. A system of three such
coupes combined is taken to be fundamental unit for the entire
study to resemble effect to an entire coach. Analysis of flame and
soot contours and concentrations is done corresponding to variations
in heat release rate per unit volume (HRRPUA) of fire source,
variations in conditioned air velocity being circulated inside coupes
by vents and an alternate fire initiation and propagation mechanism
via ducts. Quantitative results of fractional area in top and front
view of the three coupes under fire and smoke are obtained using
MATLAB (IMT). Present simulations and its findings will be useful
for organizations like Commission of Railway Safety and others in
designing and implementing safety and evacuation measures.
Abstract: Propagation of fire through a non-air conditioned
railway compartment is studied by virtue of numerical simulations.
Simultaneous computational fire dynamics equations, such as
Navier-Stokes, lumped species continuity, overall mass and energy
conservation, and heat transfer are solved using finite volume based
(for radiation) and finite difference based (for all other equations)
solver, Fire Dynamics Simulator (FDS). A single coupe with an eight
berth occupancy is used to establish the numerical model, followed
by the selection of a three coupe system as the fundamental unit
of the locomotive compartment. Heat Release Rate Per Unit Area
(HRRPUA) of the initial fire is varied to consider a wide range of
compartmental fires. Parameters, such as air inlet velocity relative
to the locomotive at the windows, the level of interaction with the
ambiance and closure of middle berth are studied through a wide
range of numerical simulations. Almost all the loss of lives and
properties due to fire breakout can be attributed to the direct or
indirect exposure to flames or to the inhalation of toxic gases and
resultant suffocation due to smoke and soot. Therefore, the temporal
stature of fire and smoke are reported for each of the considered
cases which can be used in the present or extended form to develop
guidelines to be followed in case of a fire breakout.
Abstract: The spreadsheet engine is exploited via a non-conventional mechanism to enable novel worksheet solver functions for computational calculus. The solver functions bypass inherent restrictions on built-in math and user defined functions by taking variable formulas as a new type of argument while retaining purity and recursion properties. The enabling mechanism permits integration of numerical algorithms into worksheet functions for solving virtually any computational problem that can be modelled by formulas and variables. Several examples are presented for computing integrals, derivatives, and systems of deferential-algebraic equations. Incorporation of the worksheet solver functions with the ubiquitous spreadsheet extend the utility of the latter as a powerful tool for computational mathematics.
Abstract: In designing a low-energy-consuming buildings, the heat transfer through a large glass or wall becomes critical. Multiple layers of the window glasses and walls are employed for the high insulation. The gravity driven air flow between window glasses or wall layers is a natural heat convection phenomenon being a key of the heat transfer. For the first step of the natural heat transfer analysis, in this study the development and application of a finite volume method for the numerical computation of viscous incompressible flows is presented. It will become a part of the natural convection analysis with high-order scheme, multi-grid method, and dual-time step in the future. A finite volume method based on a fully-implicit second-order is used to discretize and solve the fluid flow on unstructured grids composed of arbitrary-shaped cells. The integrations of the governing equation are discretised in the finite volume manner using a collocated arrangement of variables. The convergence of the SIMPLE segregated algorithm for the solution of the coupled nonlinear algebraic equations is accelerated by using a sparse matrix solver such as BiCGSTAB. The method used in the present study is verified by applying it to some flows for which either the numerical solution is known or the solution can be obtained using another numerical technique available in the other researches. The accuracy of the method is assessed through the grid refinement.
Abstract: Modelica has many advantages and it is very useful in modeling and simulation especially for the multi-domain with a complex technical system. However, the big obstacle for a beginner is to understand the basic concept and to build a new system model for a real system. In order to understand how to solve the simple circuit model by hand translation and to get a better understanding of how modelica works, we provide a detailed explanation about solver ordering system in horizontal and vertical sorting and make some proposals for improvement. In this study, some difficulties in using modelica software with the original concept and the comparison with Finite Element Method (FEM) approach is discussed. We also present our textual modeling approach using FEM concept for acausal and causal model construction. Furthermore, simulation results are provided that demonstrate the comparison between using textual modeling with original coding in modelica and FEM concept.
Abstract: A computational study on bio-inspired NACA634-021 hydrofoils with leading-edge protuberances has been carried out to investigate their hydrodynamic flow control characteristics at a Reynolds number of 14,000 and different angles-of-attack. The numerical simulations were performed using ANSYS FLUENT and based on Reynolds-Averaged Navier-Stokes (RANS) solver mode incorporated with k-ω Shear Stress Transport (SST) turbulence model. The results obtained indicate varying flow phenomenon along the peaks and troughs over the span of the hydrofoils. Compared to the baseline hydrofoil with no leading-edge protuberances, the leading-edge modified hydrofoils tend to reduce flow separation extents along the peak regions. In contrast, there are increased flow separations in the trough regions of the hydrofoil with leading-edge protuberances. Interestingly, it was observed that dissimilar flow separation behaviour is produced along different peak- or trough-planes along the hydrofoil span, even though the troughs or peaks are physically similar at each interval for a particular hydrofoil. Significant interactions between adjacent flow structures produced by the leading-edge protuberances have also been observed. These flow interactions are believed to be responsible for the dissimilar flow separation behaviour along physically similar peak- or trough-planes.
Abstract: The agenda of showing the scheduled time for
performing certain tasks is known as timetabling. It is widely used in
many departments such as transportation, education, and production.
Some difficulties arise to ensure all tasks happen in the time and
place allocated. Therefore, many researchers invented various
programming models to solve the scheduling problems from several
fields. However, the studies in developing the general integer
programming model for many timetabling problems are still
questionable. Meanwhile, this thesis describes about creating a
general model which solves different types of timetabling problems
by considering the basic constraints. Initially, the common basic
constraints from five different fields are selected and analyzed. A
general basic integer programming model was created and then
verified by using the medium set of data obtained randomly which is
much similar to realistic data. The mathematical software, AIMMS
with CPLEX as a solver has been used to solve the model. The model
obtained is significant in solving many timetabling problems easily
since it is modifiable to all types of scheduling problems which have
same basic constraints.
Abstract: The need to save time and cost of soil testing at the
planning stage of road work has necessitated developing predictive
models. This study proposes a model for predicting the dry density of
lateritic soils stabilized with corn cob ash (CCA) and blended cement
- CCA. Lateritic soil was first stabilized with CCA at 1.5, 3.0, 4.5 and
6% of the weight of soil and then stabilized with the same
proportions as replacement for cement. Dry density, specific gravity,
maximum degree of saturation and moisture content were determined
for each stabilized soil specimen, following standard procedure.
Polynomial equations containing alpha and beta parameters for CCA
and blended CCA-cement were developed. Experimental values were
correlated with the values predicted from the Matlab curve fitting
tool, and the Solver function of Microsoft Excel 2010. The correlation
coefficient (R2) of 0.86 was obtained indicating that the model could
be accepted in predicting the maximum dry density of CCA stabilized
soils to facilitate quick decision making in roadworks.
Abstract: Operation enhancement in an air cooler depends on
rate of heat transfer, and pressure drop. In this paper for a given heat
duty, study of the effects of FPI (Fin Per Inch) and fin type (circular
and hexagonal fins) on heat transfer, and pressure drop in an air
cooler in Iran, Arvand petrochemical. A program in EES
(Engineering Equations Solver) software moreover, Aspen B-JAC
and HTFS+ softwares are used for this purpose to solve governing
equations. At first the simulated results obtained from this program is
compared to the experimental data for two cases of FPI. The effects
of FPI from 3 to 15 over heat transfer (Q) to pressure drop ratio
(Q/Δp ratio). This ratio is one of the main parameters in design, and
simulation heat exchangers. The results show that heat transfer (Q)
and pressure drop increase with increasing FPI steadily, and the Q/Δp
ratio increases to FPI=12 and then decreased gradually to FPI=15,
and Q/Δp ratio is maximum at FPI=12. The FPI value selection
between 8 and 12 obtained as a result to optimum heat transfer to
pressure drop ratio. Also by contrast, between circular and hexagonal
fins results, the Q/Δp ratio of hexagonal fins more than Q/Δp ratio of
circular fins for FPI between 8 and 12 (optimum FPI)
Abstract: Carefully scheduling the operations of pumps can be
resulted to significant energy savings. Schedules can be defined
either implicit, in terms of other elements of the network such as tank
levels, or explicit by specifying the time during which each pump is
on/off. In this study, two new explicit representations based on timecontrolled
triggers were analyzed, where the maximum number of
pump switches was established beforehand, and the schedule may
contain fewer switches than the maximum. The optimal operation of
pumping stations was determined using a Jumping Particle Swarm
Optimization (JPSO) algorithm to achieve the minimum energy cost.
The model integrates JPSO optimizer and EPANET hydraulic
network solver. The optimal pump operation schedule of VanZyl
water distribution system was determined using the proposed model
and compared with those from Genetic and Ant Colony algorithms.
The results indicate that the proposed model utilizing the JPSO
algorithm is a versatile management model for the operation of realworld
water distribution system.
Abstract: This paper presents a fully Lagrangian coupled
Fluid-Structure Interaction (FSI) solver for simulations of
fluid-structure interactions, which is based on the Moving Particle
Semi-implicit (MPS) method to solve the governing equations
corresponding to incompressible flows as well as elastic structures.
The developed solver is verified by reproducing the high velocity
impact loads of deformable thin wedges with three different materials
such as mild steel, aluminium and tin during water entry. The present
simulation results for aluminium are compared with analytical solution
derived from the hydrodynamic Wagner model and linear Wan’s
theory. And also, the impact pressure and strain on the water entry
wedge with three different materials, such as mild steel, aluminium
and tin, are simulated and the effects of hydro-elasticity are discussed.
Abstract: In this article, we deal with a variant of the classical
course timetabling problem that has a practical application in many
areas of education. In particular, in this paper we are interested in
high schools remedial courses. The purpose of such courses is to
provide under-prepared students with the skills necessary to succeed
in their studies. In particular, a student might be under prepared in
an entire course, or only in a part of it. The limited availability
of funds, as well as the limited amount of time and teachers at
disposal, often requires schools to choose which courses and/or which
teaching units to activate. Thus, schools need to model the training
offer and the related timetabling, with the goal of ensuring the
highest possible teaching quality, by meeting the above-mentioned
financial, time and resources constraints. Moreover, there are some
prerequisites between the teaching units that must be satisfied. We
first present a Mixed-Integer Programming (MIP) model to solve
this problem to optimality. However, the presence of many peculiar
constraints contributes inevitably in increasing the complexity of
the mathematical model. Thus, solving it through a general-purpose
solver may be performed for small instances only, while solving
real-life-sized instances of such model requires specific techniques
or heuristic approaches. For this purpose, we also propose a heuristic
approach, in which we make use of a fast constructive procedure
to obtain a feasible solution. To assess our exact and heuristic
approaches we perform extensive computational results on both
real-life instances (obtained from a high school in Lecce, Italy) and
randomly generated instances. Our tests show that the MIP model is
never solved to optimality, with an average optimality gap of 57%.
On the other hand, the heuristic algorithm is much faster (in about the
50% of the considered instances it converges in approximately half of
the time limit) and in many cases allows achieving an improvement
on the objective function value obtained by the MIP model. Such an
improvement ranges between 18% and 66%.
Abstract: Singular value decomposition based optimisation of
geometric design parameters of a 5-speed gearbox is studied. During
the optimisation, a four-degree-of freedom torsional vibration model
of the pinion gear-wheel gear system is obtained and the minimum
singular value of the transfer matrix is considered as the objective
functions. The computational cost of the associated singular value
problems is quite low for the objective function, because it is only
necessary to compute the largest and smallest singular values (μmax
and μmin) that can be achieved by using selective eigenvalue solvers;
the other singular values are not needed. The design parameters are
optimised under several constraints that include bending stress,
contact stress and constant distance between gear centres. Thus, by
optimising the geometric parameters of the gearbox such as, the
module, number of teeth and face width it is possible to obtain a
light-weight-gearbox structure. It is concluded that the all optimised
geometric design parameters also satisfy all constraints.
Abstract: The numerical simulation has made tremendous
advances in investigating the blood flow phenomenon through elastic
arteries. Such study can be useful in demonstrating the disease
progression and hemodynamics of cardiovascular diseases such as
atherosclerosis. In the present study, patient specific case diagnosed
with partially stenosed complete right ICA and normal left carotid
bifurcation without any atherosclerotic plaque formation is
considered. 3D patient specific carotid bifurcation model is generated
based on CT scan data using MIMICS-4.0 and numerical analysis is
performed using FSI solver in ANSYS-14.5. The blood flow is
assumed to be incompressible, homogenous and Newtonian, while
the artery wall is assumed to be linearly elastic. The two-way
sequentially coupled transient FSI analysis is performed using FSI
solver for three pulse cycles. The hemodynamic parameters such as
flow pattern, Wall Shear Stress, pressure contours and arterial wall
deformation are studied at the bifurcation and critical zones such as
stenosis. The variation in flow behavior is studied throughout the
pulse cycle. Also, the simulation results reveal that there is a
considerable increase in the flow behavior in stenosed carotid in
contrast to the normal carotid bifurcation system. The investigation
also demonstrates the disturbed flow pattern especially at the
bifurcation and stenosed zone elevating the hemodynamics,
particularly during peak systole and later part of the pulse cycle. The
results obtained agree well with the clinical observation and
demonstrates the potential of patient specific numerical studies in
prognosis of disease progression and plaque rupture.
Abstract: In this paper, we introduced a gradient-based inverse
solver to obtain the missing boundary conditions based on the
readings of internal thermocouples. The results show that the method
is very sensitive to measurement errors, and becomes unstable when
small time steps are used. The artificial neural networks are shown to
be capable of capturing the whole thermal history on the run-out
table, but are not very effective in restoring the detailed behavior of
the boundary conditions. Also, they behave poorly in nonlinear cases
and where the boundary condition profile is different.
GA and PSO are more effective in finding a detailed
representation of the time-varying boundary conditions, as well as in
nonlinear cases. However, their convergence takes longer. A
variation of the basic PSO, called CRPSO, showed the best
performance among the three versions. Also, PSO proved to be
effective in handling noisy data, especially when its performance
parameters were tuned. An increase in the self-confidence parameter
was also found to be effective, as it increased the global search
capabilities of the algorithm. RPSO was the most effective variation
in dealing with noise, closely followed by CRPSO. The latter
variation is recommended for inverse heat conduction problems, as it
combines the efficiency and effectiveness required by these
problems.
Abstract: The main aim of the presented experiments is to
improve behaviour of sandwich structures under dynamic loading,
such as crash or explosion. This paper describes experimental
investigation on the response of new advanced materials to low and
high velocity load. Blast wave energy absorbers were designed using
two types of porous lightweight raw particle materials based on
expanded glass and ceramics with dimensions of 0.5-1 mm,
combined with polymeric binder. The effect of binder amount on the
static and dynamic properties of designed materials was observed.
Prism shaped specimens were prepared and loaded to obtain physicomechanical
parameters – bulk density, compressive and flexural
strength under quasistatic load, the dynamic response was determined
using Split Hopkinson Pressure bar apparatus. Numerical
investigation of the material behaviour in sandwich structure was
performed using implicit/explicit solver LS-Dyna. As the last step,
the developed material was used as the interlayer of blast resistant
litter bin, and it´s functionality was verified by real field blast tests.
Abstract: In this paper we present the efficient parallel
implementation of elastoplastic problems based on the TFETI (Total
Finite Element Tearing and Interconnecting) domain decomposition
method. This approach allow us to use parallel solution and compute
this nonlinear problem on the supercomputers and decrease the
solution time and compute problems with millions of DOFs. In
our approach we consider an associated elastoplastic model with
the von Mises plastic criterion and the combination of linear
isotropic-kinematic hardening law. This model is discretized by
the implicit Euler method in time and by the finite element
method in space. We consider the system of nonlinear equations
with a strongly semismooth and strongly monotone operator. The
semismooth Newton method is applied to solve this nonlinear
system. Corresponding linearized problems arising in the Newton
iterations are solved in parallel by the above mentioned TFETI. The
implementation of this problem is realized in our in-house MatSol
packages developed in MatLab.
Abstract: Flow blockages referring to the increase in flow are
being considered as a vital equipment for marine current energy
conversion. However, the shape of these devices will result in
extracted energy under the operation. The present work investigates
the effect of two configurations of a grating, convergent and
divergent that located upstream, to the water flow velocity. The flow
characteristics are studied by Computational Fluid Dynamic
simulation by using the ANSYS Fluent solver for these specified
arrangements of the grating. The results indicate that distinguished
characteristics of flow velocity between “convergent” and
“divergent” grating placements is up to 10% in confined conditions.
Furthermore, the velocity in case of convergent grating is higher
than that of divergent grating.
Abstract: In order to study the aerodynamic performance of a
semi-flexible membrane wing, Fluid-Structure Interaction simulations
have been performed. The fluid problem has been modeled using
two different approaches which are the vortex panel method and the
numerical solution of the Navier-Stokes equations. Nonlinear analysis
of the structural problem is performed using the Finite Element
Method. Comparison between the two fluid solvers has been made.
Aerodynamic performance of the wing is discussed regarding its
lift and drag coefficients and they are compared with those of the
equivalent rigid wing.
Abstract: Diffusion stills have been effective in water
desalination. The present work represents a model of the distillation
process by using vertical single-effect diffusion stills. A semianalytical
model has been developed to model the process. A
software computer code using Engineering Equation Solver EES
software has been developed to solve the equations of the developed
model. An experimental setup has been constructed, and used for the
validation of the model. The model is also validated against former
literature results. The results obtained from the present experimental
test rig, and the data from the literature, have been compared with the
results of the code to find its best range of validity. In addition, a
parametric analysis of the system has been developed using the
model to determine the effect of operating conditions on the system's
performance. The dominant parameters that affect the productivity of
the still are the hot plate temperature that ranges from (55- 90°C) and
feed flow rate in range of (0.00694-0.0211 kg/m2-s).