13
13
net/publication/382264625
CITATIONS READS
0 197
2 authors, including:
SEE PROFILE
All content following this page was uploaded by Shehroze Tahir Khan on 15 July 2024.
1
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
2
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Preface
The 13th International Mechanical Engineering Conference was held on 6th & 7th March 2024, an annual
conference collaboratively organized by Institution of Engineers Pakistan (IEP), NED University of
Engineering & Technology, Karachi and Project GetInnovative4Impact. It aims to bring together people
who are working on the Industry 4.0 and related domains.
This year, the IMEC conference was structured around 8 tracks on various thematic topics, including
IMEC-24 received 30 papers from all over the world. All submissions were peer-reviewed and also
evaluated during the presentation of authors on conference days. We are very pleased to report that the
quality of the submissions this year turned out to be very high. A total of 21 papers were accepted as full
papers which is around 70% full paper acceptance rate.
We acknowledge the invaluable assistance of the track chairs and the technical review committee members.
A complete list of all committee members can be found on the previous page. Most reviewers provided
detailed and constructive comments which were valuable for the authors to continue improving their papers,
even if their submissions were not selected for the conference.
Given the high-quality works done by authors, reviewers, and track chairs, we are confident that the IMEC
2024 proceedings capture the current state-of-the-art research in the learning technology field and will have
a significant impact on the research community in the longer term.
3
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Conference Proceedings
International Mechanical Engineering
Conference 2024
Table of Contents
S. No. Title of Paper / Authors Page No.
1. DESIGN AND SIMULATION OF PHOTOVOLTAIC POWER PARK FOR EVACUATING 7
SINDH SOLAR POTENTIAL USING HVDC TRANSMISSION SYSTEM
4
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Abdul Qadir, Shehroze Tahir Khan, Murtaza Mehdi, Meraj Ali Shah
15. MONITORING OF WELDING DIGITIZATION USING VISION-BASED 166
CLASSIFICATION: A COMPARATIVE ANALYSIS OF CNN AND RESNET-BASED
APPROACH
5
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Syed Abdullah Ali, Furqan, Hasnain, Abdul Rauf, Izhan, Muhammad Hateem
Arif
20. INVERSE KINEMATICS OF ROBOTIC ARM FOR WELDING DIGITIZATION IN THE 220
ARENA OF INDUSTRY 4.0
Shehroze Tahir Khan, Murtuza Mehdi, Tariq Jamil and Abdul Qadir
6
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (M.M. Khan)
ABSTRACT
Mega-sized power plants that use both conventional and unconventional energy sources and are
situated in remote locations are producing electricity. The electricity industry in Pakistan has
several issues due to outdated and inadequate domestic power plants as well as unstable, costly,
and insecure power generation. While finding adequate land to establish a Photovoltaic power Park
(PPP) is challenging, on-grid solar photovoltaic power parks (OGSPPPs) are a key step towards
addressing the energy shortfall in urban load centers. Despite this, PPPs are the most simple and
straightforward technology to deploy. The strategic challenge, though, is in determining the solar
energy potential of Pakistan's Sindh and Punjab Provinces. Sindh province is perfect for large-
scale solar installation because of its high solar irradiance and lack of fog throughout the year, the
yearly average global horizontal solar potential of Sindh is 1978.9 kWh/m2 and diffuse solar
irradiation of 871.7 kWh/m2, greater than Punjab. The ±660 KV HVDC transmission line in Sindh
is designed to transfer large amounts of power from the various Sindhi areas, which generate from
various power plants, to the load centers located in North Punjab. The findings indicate that the
PVSyst software records an annual average performance ratio of 84.8% for the Jamshoro area of
Sindh. The paper's findings show that the execution of this PPP project involves adding more
energy to load centers in order to combat energy shortages and increases in per-unit costs. This
can lead to both energy security and national economic stability.
Keywords: Photovoltaic Power Park (PPP); Solar irradiation; Power Evacuating; Sindh solar
Potential; Solar irradiance; Power Generation; On grid Solar Photovoltaic Power Park (OGSPPP);
Load centers
1. INTRODUCTION
The need for energy is growing daily as time for energy production, such as coal, natural
goes on. Nearly all homes and businesses gas, and oil, results in rising greenhouse gas
require electricity, which puts a heavy emissions that have a negative effect on the
company. The use of traditional fossil fuels environment and drive global warming [1,2].
7
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Cleaner energy sources including sun, wind, horizontal irradiance is decreasing and peaks
biomass, magneto hydrodynamics (MHD), in the Himalayas. While other portions of
and hydropower must be used in order to Pakistan also have significant solar potential,
replace traditional power generating methods with about 1500 kWh/m2 per year,
with new ones. These resources are encompassing more than 90% of Pakistan's
renewable since nature replenishes them total area, Figure 1 [5] shows that the
maximum solar radiation in Sindh as a whole
constantly [3]. Solar energy is among the
is slightly over 2300 kWh/m2 per year. This
most practical answers to the global energy
energy can be captured efficiently enough to
crisis. Due to the expensive impossibility of meet all of the country's energy
connecting power lines to remote and rural requirements. With Pakistan's population
areas, PPP systems have emerged as a expanding at an exponential rate and
dependable global energy generating source conventional energy supplies like gas,
[4]. furnace oil, and coal running out quickly,
photovoltaic energy is the only dependable
substitute for traditional energy sources. In
this context, a number of actions have
already been taken. To maximize the benefit
of the solar resource, system parameters must
be sized and simulated for an efficient energy
production. The energy production and
necessary PV system size may be assessed
using the modeling tool PVSyst [6,7].
Numerous studies have used PVSyst to
assess the system's performance. The
efficacy of a grid-connected, 190kWp solar
power plant was investigated by the authors
[9]. PVSyst was used to analyze the electrical
power produced by solar arrays and
photovoltaic systems in grid-tied systems in
Berlin and Kathmandu, together with all
conversion losses and losses from the
inverter system. Regarding PVSyst simulator
usage PV system-based power generating.
The authors used PVSyst to model a grid tie
system in order to determine if building a 1
MW solar photovoltaic (PV) power station
would be feasible. They achieved this by
Figure 1. Pakistan's yearly global
simulating the southern Tamandua area of
horizontal irradiance (GHI) in kWh/m2.
India and comparing performance analyses
In Pakistan's Sindh and Baluchistan, where from various geographic locations [8-11].
horizontal surface global irradiations are The objective of this research is to utilize
highest, direct solar radiation presents an PVSyst simulation to build and model a grid-
enormous opportunity for clean, green connected PV system for a specific
electrical energy production. In Pakistan's geographic area in Jamshoro, Sindh, and its
northern and northeastern areas, global potential transmitted through the already
8
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
The PVSyst program receives the following The placement of photovoltaic (PV) systems
inputs: the Meteonorm-8.1 data, radiations is crucial for solar energy generation since
from solar panels, appropriate design for different parts of the planet receive different
photovoltaic (PV) cells and arrays, quantities of solar radiation on solar cells. A
positioning of the photovoltaic field, Losses unique collection of location-specific
in the inverter system operational metrics, such as latitude, longitude, and
circumstances, efficiency and utilization of altitude, may be used to describe this
energy, index of normalized performance. positional difference. The proposed location
10
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
11
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
13
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
[5] Solar modeling report, The World bank [9] Y. Sharma, and S. S. Chandel,
Solar resource mapping in Pakistan, Performance analysis of a 190 kWp grid
https://documents.worldbank.org/pt/publ interactive solar photovoltaic power plant
ication/documentsreports/documentdetai in India, Energy, 55(1), 476-485, 2013.
l/200771468096863184/solar-resource- [10]
mapping-in-pakistan-solar-modeling- . M. Khan, S. Ahmad, S. Jabbar, S.
report. 2015. Baloch, M. A. Shafi, R. Nazeer and J.
[6] M. Asim, S. Kamran, M. Farooq, M. Faiz, Simulation design of Grid Tied
Amjad, A. Qamar, Modelling of Photovoltaic (PV) system of a 1.05 kWp
optimized standalone PV system for basic DC for a geographical location of
domestic energy use in Pakistan, Journal Tandojam, Sindh, 1st International
of Faculty of Engineering & Technology, Conference on Women Development in
22(2) 61-68, 2022. Engineering, Science & Technology
[7] M. A. Shafi, S. Bibi, M. M. Khan, B. M. (WD-EST’23), Pakistan 2023.
Soucase, A numerical simulation for [11]
efficiency enhancement of CZTS based . M. Khan, S. Ahmad, H. Sikandar, B.
thin film solar cell using SCAPS-1D, Akram, A. Raza, R. G. Hassan, M.
East European Journal of Physics, 5(1), Hussain, Designing and Performance
52–63, 2022. analysis of Photovoltaic system using
[8] M. A. Shafi, M. M. Khan and B. M. PVSyst Software: A geographical
Soucase, Effect of parasitic parameters location of Pakistan, TechXchange
and environmental conditions on IV and International Conference, QUEST
PV characteristics of 1D5P model solar Nawabshah (TechXchange 2024),
pv cell using LTSPICE-IV, East Pakistan 2024.
European journal of Physics, 5(1), 2022.
15
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
16
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
1. INTRODUCTION
Diesel engines produce their amazing
performance and exceptional fuel economy
by first compressing air to extremely high
pressures and then spraying a small amount
of fuel into this highly compressed air. Diesel
engines are still plagued by excessive
emissions of nitrogen oxides (NOX) and Figure 1. (a) SEM CNTs
particulate matter (PM) [1,2]. New Rather than in the gaseous phase, the nano
technologies, like fuel post-combustion additive was largely released in the particle
emission control devices, will aid future phase. The main component of all the
pollution restrictions. It was more additions is metal, which is discharged as
economical to use fuel additives to lower the oxide in the exhaust. Additives were
number of particles smaller than 2.5μm. beneficial for trap renewal but had little effect
strategy. The freezing point of organic on suppressing soot. Fuel additives including
manganese was shown to decrease most at 15 catalytic compounds were among the most
°C, and 700 ppm was found to be the optimal tempting approaches to reduce PM.
dosing rate [3]. Additives are added to fuel to improve engine
When the metal dosage was administered, performance and reduce emissions. Nearly
notable patterns in the freezing point fall every technique for lowering diesel engine
were seen. Diesel was mixed with different emissions at the source compromises the
metal additions to enhance fuel quality, levels of NOX, particulates, and UHC [7].
resulting in more thorough combustion and The engine maker has several options at their
fewer pollutants from the exhaust. The disposal to precisely manage the emissions of
influence of a catalyst on the burning of any given contaminant during the
unburned hydrocarbons is the fundamental combustion and injection stages. The
idea behind this additive action [4]. The following facts and trends were noteworthy:
addition of noble or transitional metals to fuel (i) the increased injection speed led to a
lowers the soot ignition temperature. reduction in particulate emissions and
Diesel additives can also be categorized into improved air-fuel mixing; (ii) the injection
three groups based on the purpose for which pressure increase was advantageous as well;
they were designed: (i) pre-flame additives; and (iii) using an injector with pintles as flat
(ii) flame additives; and (iii) post-flame lowers residuals and emissions of particulate
additives[5]. matter and unburned hydrocarbons [8].
Particles begin to form above a predefined The additive method calls for a constant
threshold when the amount of additive addition of diesel fuel to all vehicles,
increases. The only factors influencing the excluding those equipped with filters. Metal-
additive dose limit were the type of engine organic compounds that were completely
and the soot emissions (normally aspirated dissolved in diesel fuel made up the bulk of
and turbocharged engines behave fuel additives. To decrease the ignition delay,
identically). Higher additive concentrations stabilize and function as an anti-oxidant, and
demonstrated a more concentrated benefit, operate as a surface-active agent, diesel fuel
even if lesser additive dosages were still is treated with a metal addition [9].
beneficial [6]. Metal-based additives effectively decreased
diesel emissions in two ways. The metals first
17
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
work with the soot's carbon molecules in a oil bath cleaning unit, carbon nanotubes were
straight line to reduce the oxidation combined with pure diesel (D) at
temperature, or they work with water to concentrations of 30 and 60 parts per million.
process hydroxyl radicals, which quickens Table 1 displays the experiment's design. A
the soot's oxidation [10]. The nanoparticles magnetic stirrer was used for 30 minutes to
act as the soot's oxidation nucleus if the mix the CNTs added to pure diesel to produce
additives are used in the diesel engine after a homogenous emulsion fuel. Pure diesel can
combustion. be made better by adding carbon nanotubes
Engine exhaust particles are always under (CNTs). Carbon nanotubes (CNTs) bearing
suspicion, even in small amounts, due to their Amido groups were extremely reactive and
propensity to penetrate and simply dwell in could interact with a wide range of materials.
the respiratory alveoli, with severe health
repercussions. Diesel particulates have been Diesel CNTs HE-OBCU-
Name
demonstrated to match the size of particles (%) (ppm) EGR
that were most readily maintained for the D 100 0 No
longest amount of time in the respiratory DC30 100 30 No
system among all the types of dust that may DC60 100 60 No
be found in the atmosphere [11]. DHOE 100 0 Yes
Previous research demonstrates the DC30HOE 100 30 Yes
significant impact of nanoparticle additions DC60HOE 100 60 Yes
on the efficiency and emissions of diesel
engines. In this investigation, a four-stroke, Table 1. Design of experiment
single-cylinder diesel engine was used for the
experimental work. 30 and 60 parts per 3. EXPERIMENTAL SETUP
million of each carbon nanotubes were Fig. 2 displays the diesel engine test bench's
combined with pure diesel. At engine speeds schematic diagram. A diesel engine with
of 1200, 1400, 1600, 1800, 2000, and 2200 three cylinders and water cooling, operating
rpm, the performance and emission on four strokes, was connected to an eddy
characteristics were measured. When current dynamometer. Engine speed ranges
nanoparticles were added to diesel fuel, the of 1200 to 2200 rpm at 100 rpm intervals
fuel's performance significantly improved. under full load were tested. The AVL
When pure diesel and carbon nanotubes are DiTEST emission analyser was used to
combined, the number of unburned measure the emissions of CO and NOx, while
hydrocarbons decreases by a maximum of the AVL smoke meter was used to measure
8.43%. Emissions of carbon monoxide were the emissions of BSFC and UHC.
decreased by 20.34% using carbon nanotubes
at 1600 rpm engine speed. When both
manganese oxide and carbon nanotubes were
utilized at a concentration of 30 and 60 ppm
each, there was a 6.2% decrease in NOx
emissions.
2. FUEL PREPARATION
In the current study, pure diesel (D) and
carbon nanotubes were utilized in a diesel Figure 2. Schematic diagram of the
engine. With and without a heat exchanger experimental setup
18
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
19
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
4.3 CO emissions
Carbon monoxide is produced as a result of
incomplete combustion and is exacerbated by
the lack of oxidants, low temperature, and
extended residence time. Figure 5 displays
the carbon monoxide emissions from the
engine. When metal oxide additives were
catalytically oxidized in an oxygen-rich
environment, carbon monoxide emissions Figure 6. Variation of nitrogen oxide
were trending lower [14]. The graph shows emissions with speed.
that at 1600 rpm, the use of carbon nanotubes 4.5 UHC emissions
with the heat exchanger oil bath cleaning unit At 1200 rpm, the DI diesel engine's emissions
decreased carbon monoxide emissions by of unburned hydrocarbons were reduced. Up
to 1600 rpm, it rises as engine speed
increases. The highest rate of unburned
hydrocarbon emissions was 1600rpm.
Subsequently, a decrease in unburned
hydrocarbon at faster rates becomes
apparent. As a result, unburned hydrocarbon
emissions from idling engines have received
attention. The engine's unburned
0.06%. hydrocarbon exhaust pattern is shown in
Figure 7. When pure diesel fuel is mixed with
Figure 5. Variation of carbon monoxide nanoparticles, the number of unburned
with speed hydrocarbons decreases [13]. The highest
reduction in unburned hydrocarbons was
4.4 NOx emissions observed when 90 parts per million of carbon
It has been discovered that NOx emissions nanotubes were combined with pure diesel.
were on the decline when diesel fuel contains
nanoparticle additives. There was a change of
1-1.4% in the cetane number of diesel fuel
due to nanoparticles which has an impact on
emissions. Compared to fuels with lower
cetane numbers, those with higher cetane
numbers have lesser premixed fuel fractions
and fewer NOx emissions [6]. Figure 6
illustrates the engine's nitrogen oxide (NOX)
emission. The data indicates a 6.2%
Figure 7. Variation of unburned
reduction in NOx emissions when carbon
hydrocarbon with speed
nanotubes were used at a concentration of 30
and 60 ppm each. 5. CONCLUSIONS
From the results of the experimental study,
carbon nanotubes additive somewhat
improves an engine's performance. The
BSCF of the diesel engine decreased by
3.56% when the carbon nanotubes were
20
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
22
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
losses, there is a great deal of potential for produces the greatest energy decreases in the
energy savings. Cai [4] suggested that power plant as compared to the condenser,
although the first rule of thermodynamics is which results in the lowest energy losses.
frequently used to compare energy use, it is Since consumption of water and water
unable to analyze the quality of energy in the rejected from the electric power plant have
future. I.e., the point at which an evaluation decreased, the plant's air cool condenser
of energy becomes relevant. Dunbar [5] proved to be environmentally beneficial.
established that the 2nd law of Patel [9] compared various coal types,
thermodynamics, which tells us how to including imported, Indian, and a
calculate the efficient use of a particular combination of both (40 percent Indian, 60
amount of energy at a certain condition, percent imported), as well as L.S.H.S. oil,
equates to exergy. The best performance is and concluded that the boiler's energy
achieved when exergy loss in a process is the efficiencies were 37%, 37.7%, 37.8%, and
least.Kanoglu et.al, [6] suggested a greater 40.1%, respectively, based on the first
awareness of energy and energy efficiency principle of efficiencies. Ohijeagbon [10]
for power plant energy management. where concluded that, at conventional reference
the power cycle is analyzed using different state temperatures of 250°C and an
energy and energy-base efficiency. evaporating ratio of 12, the energy & exergy
Examined are the cogeneration, geothermal, efficiency attained for the complete boiler
and vapour and gas power cycles; several was 69.56% & 38.57%. The analytical
cycle designs were taken into consideration. method served as the foundation for all
Som [7] determined that the energy balance calculations. Gulhane [11] ] determined that
and thermodynamic irreversibility in the boiler's exergy destruction at a residential
combustion processes, which are present in load of 1.1 MW was around 83.35%, and that
practically all circumstances, are significant. figure increased to 76.33% as the load grew
The internal energy exchange connected to to 5.6 MW. As a result, the first and second
high-temperature gradients brought on by the laws' efficiency rose with load; therefore, to
release of heat in combustion reactions was lessen irreversibility, we must focus on peak
the main cause of irreversibility. Reducing load. Kanoglu [12] found a diesel engine
the irreversibility of the conduction of heat power plant's performance characteristics.
was the main strategy for keeping the energy Energy conversion & management
dissipation in the burning process within an concluded that energy appeared to be a key
acceptable bound. Saidur [1] suggested a idea. It establishes a link between the
framework of industrial boiler energy, engineering and physical worlds and the
energy, and economic analysis. 72.46% and surrounding environment and demonstrates
24.89%, respectively, were found for energy the true effectiveness of engineering systems,
and energy efficiency. The primary source of making it an invaluable concept to identify
energy destruction was the combustion areas for development. Energy calculations,
chamber, which was supported by the flue the opposite of the current material & energy
gases. Boiler fan motors were equipped with balances, can provide deeper and more
variable-speed drives to reduce energy insightful insights into the process. It can be
consumption. Regulagadda [8] It is suggested emphasized that energy analysis has a great
that based on the thermodynamic study of a deal of potential usefulness concerning
subcritical boiler-turbine generator of a 32 energy usage and that energy plays a crucial
MW coal-fired electric power plant, the role in the process of developing energy
irreversibility of the boiler and turbine policies. Kanoglu [6] To study, increase, and
24
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
improve energy systems through appropriate efficiency for a boiler, are described in the
strategy and policy, both energy and energy section that follows. Cost-benefit analyses of
efficiency analyses are required. Regular energy-saving measures to lower the
actions may assist increase the performance temperature of exhaust gases at the HRSG
of industrial boilers if such regulations were output and employ VSD to lower boiler fan
in place. Som [7] employed the hot-gas energy consumption are covered.
recovery systems to enhance the boiler's
existing specifications and reported that a Block Schematic:
heat recovery system might save as much as The KAPCO power plant unit 3 block
55% of the cost. The energy, energy schematic is displayed in Figure 1.
efficiency, energy losses, and energy damage Bypass Stack
associated with boilers are identified in this Steam
Stack
study, and strategies for enhancing effective HRSG
Turbine
heat transfer by reducing boiler energy Generator
Condenser
consumption through the use of variable
speed drives are suggested. Analysis has also
been done on the benefits related to energy Pump
and the economy. Many other techniques, Pump
combustor
such as managing excess air, increasing Cooling
consumption efficiency, using Desecrator
Tower
environmentally friendly fuel, recovering
excess heat, recovering condensate, Generato
Compressor Gas
improving the blowdown process, preventing Turbine
r
escape, and offering appropriate protection,
can all help reduce the amount of energy used Figure 1. Block schematic
by boilers. 2.1. The equation for energy analysis
The heat exchange system and the combustor
This paper's primary goal is to examine the are the two components that make up a boiler.
boiler's parameters and pinpoint the main Below is a discussion of these two
energy and exergy losses in the boiler system. components' energy and exergy analyses.
Energy & exergy flows have also been Because a boiler's combustion is usually
demonstrated. Numerous energy-saving well-insulated, very little heat is lost to the
strategies have been developed, and energy environment. The loss of energy will become
and energy efficiencies have also been equal with the difference between the
calculated. It has also been established that enthalpies of the reactants and products,
employing variable-speed drives & according to the 1st law of thermodynamics.
recovering heat from flue gases can save The combustor's energy balance is provided
energy. The case study concludes that energy by:
analysis provides a more precise performance 𝑑𝐸𝑠𝑦𝑠𝑡𝑒𝑚
for various devices. 𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 =
𝑑𝑡
= 0 =≫ 𝑠𝑡𝑒𝑎𝑑𝑦 (1)
2. METHODOLOGY AND DATA 𝐸𝑖𝑛 = 𝐸𝑜𝑢𝑡
COLLECTION 𝐸𝑙𝑜𝑠𝑠 = 𝑚̇1 ℎ1 + 𝑚̇2 ℎ2 − 𝑚̇3 ℎ3 (2)
The methods used to calculate the energy & Where,
exergy efficiency and the destruction of a Ein = enrgy in
boiler, as well as the energy & exergy Eloss = enrgy loss
25
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
𝐸𝑜𝑢𝑡
Eout = enrgy out 𝜂𝑐 = (6)
𝐸𝑖𝑛
ℎ1 = fuel′s specific enthalpy (KJ/kg)
Where,
ℎ2 = air's specific enthalpy (KJ/kg)
𝜂𝑐 = efficiency of a HRSG
ℎ3 =Particular enthalpy of heated
combustion products
2.2. The exergy analysis equation
ṁ1 = mass fuel flow rate
By assuming a zero rate of change in exergy
ṁ2 = air mass flow rate
within the boiler's system and applying an
𝑚̇3 =thermal products of combustion's
exergy balance to the boiler, the maximum
mass flow rate power output, or reversible power, can be
For each fluid stream, a heat exchanger
calculated. This is done while keeping in
typically experiences negligible potential and
mind the boundary temperature of 25 °C.
kinetic energy changes as well as no work
According to the second rule of
interactions (w=0). In general, the heat
thermodynamics:
exchanger's outer shell is well-insulated to 𝑋𝑖𝑛 − 𝑋𝑜𝑢𝑡 − 𝑋𝑑𝑒𝑠𝑡𝑟𝑜𝑦𝑒𝑑 =
stop heat loss into the surrounding medium. 𝑑𝑋𝑠𝑦𝑠𝑡𝑒𝑚
= 0 ==> 𝑠𝑡𝑒𝑎𝑑𝑦 (7)
A tiny quantity of heat will, nevertheless, 𝑑𝑡
Exergy efficiency is determined using the A device called a variable frequency drive
second rule of thermodynamics as follows: modifies the power supply's frequency to
𝑋 control the motor's speed. The frequency and
𝜓𝐻 = 1 − 𝑑𝑒𝑠𝑡𝑟𝑜𝑦𝑒𝑑 (14)
𝑋 𝑖𝑛 the driven equipment's speed are
For analysis of energy, the exergy data we proportionate. Energy savings are substantial
took from the plant is given in Table 1. even with little speed reductions since the
power required to run pumps & motors varies
Flow rate Temper Enthalpy Entropy with a cube of the speed. By directly
(KJ/s) ature (KJ/Kg) (KJ/Kg. increasing motor voltage with frequency, the
(°C) K) percentage of applied motor voltages to
426 340.9 621.64 2.430 deliver frequency remains constant, allowing
9.72 110 43.890 1.97 for the achievement of maximum constant
435.72 1050 1424 3.29 torque. Constantly volts per Hertz is the
common word used to describe this kind of
435.72 506 1002 2.710 control.
Electrical energy is saved when an
52.19 472.4 4013 7.72 alternating current (AC) motor is slowed
down by the use of an inverter in a variable-
speed drive. In addition to reducing stack
52.19 120 554.92 1.70
losses through the limited extra air rate and
435 160 432.42 2.07 the addition of a driver to that system, these
measures also save electrical energy.
Table 1. Exergy data from power plant Because of this, running the motor at a
Boiler energy loss can be minimized in a variable speed will save electrical energy in
boiler by employing various techniques. addition to increasing boiler efficiency.
However, this thesis has taken into Oxygen prolongs the duration of the plant's
consideration energy saving through the boiler and saves fuel and emissions.
recovery of heat from exhaust gases in a Although this is an objective modification to
boiler as well as energy saving through the the motor, a motor's speed can also be altered
use of variable speed drives to lower the by changing the number of poles. Rewinding
speed of the boiler fan. it would be necessary, causing the speed to
alter step by step. Therefore, the ideal answer
2.3. Concept of variable speed drive
is to change the frequency for convenience,
One of the best-organized control strategies
economy, and accuracy. It is necessary to
is VSD. Under specific conditions, it
adjust the volts per Hz ratio to alter the motor
provides only the necessary power to
torque. An output of numerous frequencies is
overcome system resistance. In
provided by a drive. The use of variable
contemporary industrial and commercial
speed drives can be summed up as follows:
boilers, variable-speed drives are typically
lower electricity usage, less need for a
employed. When regular low-load times are
separate motor starter, better fan control, and
required due to work conditions, it works
longer equipment lifespan. By regulating the
well. If you lower the fan speed, you can
surplus air, a VSD keeps fuel efficiency high
improve combustion efficiency at low fire. A
under part load settings.
crucial factor affecting all kinds of
For a boiler with VSD, the annual energy
combustion apparatus is the ambient air
savings are provided by:
temperature. Boiler efficiency may be
significantly impacted by this.
27
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
28
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
table and the assumptions outlined in the the potential and kinetic energies are minimal
previous subject, the energy input may be and there is no assumption of a work
computed as follows: interaction, the information from the energy
destruction table will be evaluated as follows:
𝐸𝑙𝑜𝑠𝑠 = 𝑚̇𝑓 ℎ𝑓 + 𝑚̇𝑎 ℎ𝑎 − 𝑚̇𝑝 ℎ𝑝 (18) 𝑙𝑐 = (𝑚̇1 𝜖1 + 𝑚̇2 𝜖2 ) − (𝑚̇3 𝜖3 ) (22)
= 𝑚̇1 (ℎ1 − 𝑇0 𝑆1 ) + 𝑚̇2 (ℎ2 − 𝑇0 𝑆2 ) − 𝑚̇3 (ℎ3 −
𝐸𝑙𝑜𝑠𝑠 = 70964.16𝐾𝑊 𝑇0 𝑆3 )
29
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
Load Tim Power Energy Pow Energy Saving industrial boilers," Energy Policy, vol.
(%) e Load Use er use (Kwh/ 38, pp. 2188-2197, 2010.
(hr) (witho (Kwh/ Lo (Kwh/ yr) [2] Gupta, Ghai, and Jain, "Energy efficiency
ut yr) ad yr) improvement strategies for industrial
VSD) (with boilers: a case study," Journal of
(Kw) VSD Engineering and Technology, vol. 1, p.
) 52, 2011.
(Kw [3] Rosen,"The relation between
) thermodynamic losses and capital costs
100 40 2.9 1160 2.9 1160 0 for a modern coal-fired electrical
0 generating station," in Proceedings of
80 90 2.9 2610 2.2 1980 630 computer-aided energy systems analysis,
0 ed, 1990.
70 80 2.9 2320 1.9 1520 800 [4] Cai and Kim, "First law of
0 thermodynamics and Friedmann
60 82 2.9 2378 1.4 1148 1230 equations of Friedmann-Robertson-
0 Walker universe," Journal of High
Total 2920 8668 5808 2660 Energy Physics, vol. 2005, p. 050, 2005.
Table 4. Power saving by using VSD [5] Dunbar, Lior, and Gaggioli, "The
component equations of energy and
Using VSD, an annual total energy savings of
exergy," Journal of Energy Resources
2660 kWh can be attained, assuming a per-
Technology, vol. 114, pp. 75-83, 1992.
KWh charge of Rs. 29.78 Pkr. This equates
[6] Kanoglu, Dincer, and Rosen,
to Rs. 79,214.8 Pkr annually.
"Understanding energy and exergy
4. CONCLUSION efficiencies for improved energy
Component-wise energy and exergy analysis management in power plants," Energy
was done for the boiler at unit 3 KAPCO Policy, vol. 35, pp. 3967-3978, 2007.
power plant. According to the first law, [7] Som and Datta, "Thermodynamic
HRSG experiences the highest heat losses, irreversibility’s and exergy balance in
which total 171.2 MW. Exergy analysis combustion processes," Progress in
represents the more accurate performance of energy and combustion science, vol. 34,
thermodynamic devices. Exergy calculation pp. 351-376, 2008.
shows that maximum exergy destruction [8] Regulagadda, Dincer, and Naterer,
occurred in HRSG whereas destruction was "Exergy analysis of a thermal power plant
calculated at 224.15MW. with measured boiler and turbine losses,"
Applied Thermal Engineering, vol. 30,
ACKNOWLEDGEMENT pp. 970-976, 2010.
For their support of this experimental [9] Ankit Patel "Energy and exergy analysis
investigation, the authors are grateful to the of a boiler with different fuels like Indian
management of the Swedish College of coal, imported coal and L.S.H.S oil,"
Engineering and Technology Rahim Yar IJERT, vol. 1, 2012.
Khan [10] Ohijeagbon, Waheed, and Jekayinfa,
"Methodology for the physical and
REFERENCES chemical exergetic analysis of steam
[1] Saidur, Ahamed, and Masjuki, "Energy, boilers," Energy, vol. 53, pp. 153-164,
exergy and economic analysis of 2013.
30
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
31
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (S. Y. Khan)
ABSTRACT
The objective of this work is to compare the thermal performance of spherical and elliptical cavity
receivers with a conical-shaped cavity for a solar dish cavity system. The Monte-Carlo ray-tracing
method is used to determine the radiation flux distribution and solar flux uniformity of the receiver.
The radiation flux is found to be more uniform in the spherical and elliptical cavities, and cavities
that are positioned far away from the focal plane of the solar concentrator. The result from optical
analysis is used as a boundary condition for the computational model. CFD analysis is performed
to determine thermal performance using ANSYS Fluent where convective heat loss in the cavity
receiver is calculated. Convective heat loss of the receiver with different wind speeds varying from
1 – 10 m/s is investigated in different shapes of cavity receivers. The computational results show
that the convective heat loss in spherical and elliptical cavity receivers is significantly reduced due
to the swirling action of turbulent flow inside the cavity receiver.
Keywords: Cavity receiver; Convection heat loss; CFD; Parabolic dish system; Optical analysis
around cavity receivers coupled with distribution of solar radiation flux inside the
parabolic dish structure. Roux et al. [8] cavity as well as lowering the chances of
studied a small scale solar thermal Brayton convection loss. The improvement can be
cycle coupled with solar parabolic dish done by modifying the opening with
system at the focal point of the system. The backward facing step in cavity receiver and
study involved the investigation of geometry positioning the cavity with offset distance.
The main focus of this paper is to investigate
parameter of dish such as concentrator shape,
the effect of the spherical and elliptical
diameter, rim angle, reflectivity, optical
shaped cavities with modified receiver
errors, and solar tracing errors which effect opening on the convection heat loss. The
the efficiency of the receiver. Daabo et al. [9] numerical analysis is performed using 𝑘 − 𝜔
studied different shapes of cavity receiver to SST turbulence model inside fluid domain of
determine its effect on the optical and thermal parabolic dish cavity system. The Monte-
performance of the solar cavity dish system. Carlo ray tracing (MCRT) method is useful
Liu et al. [10] investigated the design to determine the solar flux on cavity surface
parameters of heater tube arrays in cavity which can be set as wall boundary condition.
receivers to maximize the heat absorption The spherical and elliptical cavity receivers
and minimizing thermal losses. Li et al. [11] are selected to improve the performance of
investigated cavity receivers by adding the system by reducing convection heat loss.
transparent aerogel at receiver opening which
2. SIMULATION FRAMEWORK
improved the efficiency of the receiver and
The computer simulation of the cavity
minimizing the convection loss significantly. receiver system was set up in two steps: (1)
The optical performance mainly depends on optical analysis, (2) CFD analysis. The
the cavity shape and opening of the receiver. optical modeling was carried out for a 3.8 m
Thermal shocks are caused by non- diameter parabolic dish with focal length of
uniformity of solar flux inside the cavity 2.2 m where the cavity receiver was placed.
receiver which deteriorate thermal The optical analysis was performed using
performance of cavity dish system. Monte-Carlo ray-tracing method in
Furthermore, the convection loss is affected COMSOL Multiphysics software. The ray
significantly by wind patterns, orientation tracing flux map was validated through the
and opening of the cavity receiver. The effect experiments of Shuai et al. [14]. The CFD
of cavity shapes has already been analysis was performed using ANSYS Fluent
investigated by Daboo et al. [9] but the solver where a computer-generated three-
optimum cavity shape was selected based on dimensional model was discretized into finite
maximum solar flux instead of uniform flux number of volume cells and governing
distribution. Xiao et al. [12] modified the equations were applied with finite volume
receiver opening by adding quartz window to approach. Grid convergence was also
cylindrical cavity which resulted in lower performed following validation through the
convection loss. Rajan and Reddy [15] experiments of Uzair et al. [6], so that further
performed optical analysis on corrugation CFD analysis could be performed by
cavity receiver with dish concentrator and tweaking the design parameters of the cavity
varied aperture diameter to obtain optimum receiver. The optical performance of the
opening of cavity receiver. cavity receiver was characterized by
There is a need to determine such cavity concentration ratio and solar radiation flux
shapes that allow capturing uniform uniformity index. The cavity receiver was
33
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
selected based on the highest uniformity the aperture which is considered for the
index from optical analysis and then the validation of CFD model. Second cavity is
radiation flux profile was used as boundary spherical shape with same frustum opening.
conditions for the CFD analysis of cavity Third cavity is the same spherical shape
receiver. without any modification at receiver
The concentration ratio (𝐶𝑟 ) is aperture. Fourth cavity is the same spherical
estimated by taking ratio of solar
shape with a backward facing step inside the
concentrated flux on the receiver (𝑞) and the
cavity. Fifth and sixth cavities are the
ambient flux from the sun (𝐺𝑏𝑛 ) as
represented in Eq. (1). elliptical shape of eccentricity 0.3 and 0.7
respectively whereas having similar opening
𝑞 of fourth cavity.
𝐶𝑟 = (1)
𝐺𝑏𝑛
∑𝑞
𝜇= (2)
𝑛 ∙ max(𝑞)
𝜕(𝜌𝑣𝑖 )
=0 (3)
𝜕𝑥𝑖
35
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
𝜕(𝜌𝑣𝑖 𝑣𝑗 ) 𝜕 𝜕𝑣𝑖 𝜕𝑣𝑗 2 𝜕𝑣𝑘 𝜕𝑃 The default values of the model constants are:
= [𝜇 ( + − 𝛿 )] −
𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑖 3 𝑖𝑗 𝜕𝑥𝑘 𝜕𝑥𝑗 𝑀𝑡𝑜 = 0.25, 𝜉 ∗ = 1.5, 𝛽 ∗ = 0.09, 𝛽𝑖,1 =
+
𝜕
̅̅̅̅̅̅
(−𝜌𝑣 ′ ′ 0.075, 𝛽𝑖,2 = 0.0828, 𝜎𝑘1 = 1.176, 𝜎𝑘2 =
𝑖 𝑣𝑗 )
𝜕𝑥𝑗 2, 𝜎𝜔1 = 1, 𝜎𝜔2 = 1.168, 𝛼∞ ∗
= 1, 𝛼∞ =
0.52.
𝜕 1 𝜕 𝜕𝑇̅ The RANS equations also include pressure-
[𝜌𝑣𝑖 (ℎ + 𝑣𝑖 𝑣𝑗 )] = [𝑘𝑒𝑓𝑓 + 𝑣𝑖 (𝜏𝑖𝑗 )𝑒𝑓𝑓 ] velocity coupling, and this can be resolved
𝜕𝑥𝑖 2 𝜕𝑥𝑗 𝜕𝑥𝑗
𝑐𝑝 𝜇𝑡 using Semi-Implicit Method for Pressure-
𝑘𝑒𝑓𝑓 = 𝐾 + (5)
𝑃𝑟𝑡 Linked Equations (SIMPLE) [16]. The
second-order upwind scheme is used for this
The term involving −𝜌𝑣 ̅̅̅̅̅̅
′ ′
𝑖 𝑣𝑗 in Eq. (9) is finite volume solver and the discretized
known as Reynolds Stress, which is the cause RANS equations are limited to default
of turbulence, and for its closure turbulence residuals for their convergence. Also, the
equation is required. Therefore, the 𝑘 − 𝜔 under-relaxation factors are set to default for
SST model is applied to solve the non- all variables in the governing equations.
linearity in the RANS equations applied to
the fluid domain. The turbulent kinetic 2.3 Grid Independence Test
energy (𝑘) and the specific dissipation rate The simulation data is obtained after iterative
(𝜀) are represented by Eqs. (6) and (7) convergence in three different grid sizes: (1)
respectively: coarse, (2) medium and (3) fine. The
technique based on Richardson extrapolation
𝜕(𝜌𝑣𝑖 𝑘) 𝜕 𝜕𝑘 [17, 18] as prescribed by Roache [19, 20] is
= [𝛤 ] + 𝐺̅𝑘 + 𝑌𝑘 + 𝑆𝑘 (6)
𝜕𝑥𝑖 𝜕𝑥𝑗 𝑘 𝜕𝑥𝑗 used for minimizing the discretization errors
𝜕𝑣𝑗 which is known as Grid Convergence Index
where 𝐺̅𝑘 = min (−𝜌𝑣
̅̅̅̅̅̅
′ ′
𝑖 𝑣𝑗 , 10𝜌𝛽 ∗ 𝑘𝜔)
𝜕𝑥𝑖 (GCI) study. Celik et al. [21] suggested that
𝑌𝑘 = 𝜌𝛽 ∗ 𝑘𝜔 the solution should be converged at a certain
where 𝛽 ∗ = 𝛽𝑖∗ (1 + 𝜉 ∗ 𝐹(𝑀𝑡 )) iteration number and the structured mesh of
0, 𝑀𝑡 < 𝑀𝑡𝑜 hexahedral cells are suitable to calculate GCI.
𝐹(𝑀𝑡 ) = 𝑓(𝑥) = { 2 2 ,
𝑀𝑡 − 𝑀𝑡𝑜 , 𝑀𝑡 > 𝑀𝑡𝑜 The GCI study consists of several steps,
2𝑘
𝑀𝑡 = which include calculations of grid refinement
𝛾𝑅𝑇
ratio (𝑟), apparent order (𝑝), relative error (𝑒)
𝜕(𝜌𝑣𝑖 𝜔) 𝜕 𝜕𝜔 and convergence index (GCI).
= [𝛤 ] + 𝐺𝜔 − 𝑌𝜔 + 𝐷𝜔 + 𝑆𝜔 (7)
𝜕𝑥𝑖 𝜕𝑥𝑗 𝑘 𝜕𝑥𝑗 The coarse, medium, fine grids are described
𝛼 𝜕𝑣𝑗
where 𝐺𝜔 = ̅̅̅̅̅̅
(−𝜌𝑣 ′ ′
𝑖 𝑣𝑗 , 10𝜌𝛽 ∗ 𝑘𝜔) in Table 2 with the number of elements and
𝜈𝑡 𝜕𝑥𝑖
grid size (ℎ). The selection of grids is made
′ 2
𝑌𝜔 = 𝜌𝛽 𝜔 in a sequential manner with constant grid
𝐷𝜔 = 2(1 − 𝐹1 )𝜌𝜎𝜔,2
1 𝜕𝑘 𝜕𝜔 refinement ratio of 1.60. Eqs. (8) and (9)
𝜔 𝜕𝑥𝑗 𝜕𝑥𝑗
provide formulae used for grid size and grid
In above, 𝐹1 is termed as blending function refinement ratio respectively.
𝑁 1/3
and is defined as:
ℎ = [𝑁 ∑(Δ𝑉𝑖 )] (8)
𝐹1 = tanh(𝜑14 ) 𝑖=1
√𝑘 500𝑣 4𝜌𝑘
where 𝜑1 = min (max ( ), ) ℎ1 (9)
𝛽 ∗ 𝜔𝑦 𝜔𝑦 2 𝐷𝑘𝜔 𝜎𝜔 𝑦2 𝑟12 =
2
ℎ2
1 𝜕𝑘 𝜕𝜔
and 𝐷𝑘𝜔 = max (2𝜌 , 1.0 × 10−10 )
𝜎𝜔,2 𝜔 𝜕𝑥𝑗 𝜕𝑥𝑗
36
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
where 𝑁 is the total number of elements and 𝒗 𝒑 𝒆𝟏𝟐 𝒆𝟐𝟑 𝑮𝑪𝑰𝟏𝟐 𝑮𝑪𝑰𝟐𝟑
Δ𝑉𝑖 is the volume of the 𝑖th cell for a certain
grid. The grid refinement ratio is similar 1 2.97 0.04460 0.01094 1.828% 0.431%
found between medium and fine mesh by
taking ratio of medium grid size to fine grid 3 2.82 0.07743 0.02100 3.509% 0.915%
size.
Grid Elements 𝒉 (mm) 𝒓 5 1.47 0.01586 0.00077 1.993% 0.948%
1 𝑓2 − 𝑓3 (11)
𝑝= |ln | | + 𝑞(𝑝)|
ln(𝑟12 ) 𝑓1 − 𝑓2
37
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
range of error band can be calculated from optical model but it can not predict hot spots
Eqs. (13) and (14) respectively. and area of potential thermal stress within the
𝑓2 − 𝑓1 (13) receiver. The computer simulation result of
𝑓ℎ→0 ≅ 𝑓1 + 𝑝 concentration ratio calculated from MCRT
𝑟 −1
𝑅 = 𝑓ℎ→0 (1 ± 𝐺𝐶𝐼) (14) technique is compared with the experimental
data as shown in Figure 4. The CFD results
are also in good agreement with the
experiments of Uzair et al. [6] as shown in
Fig. 5, and therefore the model can be used
for predicting the performance of selected
cavity receivers.
39
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
40
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
4. CONCLUSIONS
This research paper focused to study the
effect of spherical and elliptical shapes of
cavity receiver on optical and thermal
performance of parabolic dish cavity system.
The study involved investigation on solar
flux uniformity index and convection loss in
different cavity receivers which determine
the optical and thermal performance,
respectively. The results showed that the
backward facing step opening provided
higher efficiency in both optical and thermal
performance. The spherical and elliptical
cavities elevated from focal point of the dish
also resulted in higher uniformly of solar flux
inside the cavity receiver. The convection
loss in frustum type opening in both conical
and spherical
cavity is higher as compared with the
backward facing step opening. The elliptical
cavity receiver was further examined to
investigate the effect of eccentricity. With
higher eccentricity, the solar flux uniformity
was found to be slightly higher and
convection loss is slightly reduced in this
type. Therefore, it is recommended to use
elliptical shaped cavity receiver with
backward facing step opening.
41
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
43
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
ABSTRACT
Operating rooms are considered as critical resource in hospitals as they involve maximum
investment in hospitals. Inefficient planning of patients on operating rooms underutilized their
capacity and may not generate revenue as per their investment. In recent years, Industry 4.0 has
increased the trend of industries towards intelligent and efficient decision making which has also
attracted researchers for its implementation in healthcare industry. Therefore, current research
investigates planning and scheduling problem of patients on operating rooms in context of industry
4.0 keeping in view resource constraints. A mixed integer programming model of the considered
problem is made, and a resource constraint heuristic is introduced for intelligent planning and
scheduling of patients in operating rooms to maximize the utilization of operating rooms using
theory of constraints. Experiments are performed to validate the performance proposed heuristic
as compared to the standard scheduling rules used in literature and a case problem taken from Case
Hospital is solved. Results indicate that the proposed heuristic outperforms other scheduling
heuristics.
Keywords: Operating rooms; Planning and scheduling; Industry 4.0; Heuristic
time and money and improve customer at the tactical level, master surgery schedule
satisfaction [6,7]. is made and at the operational level, patients
Industry 4.0 concepts has been used in are scheduled in each surgery block. These
planning of resources. As an example, Pham three levels of plans are interrelated to each
and Klinkert [8] considered scheduling of other, and they are interdependent.
surgical cases and decided on the time to Therefore, some researches have integrating
perform surgeries to efficiently use hospital
higher level planning with medium level
resources and to increase quality of service
planning while integration of all three levels
for patients. They proposed a mixed integer
linear programming model with an objective of planning are scarce [13, 21]. The
to minimize the makespan and minimizing hierarchical level integration of different
sum of all starting times. They solved the decision levels can give local optimal results
model on CPLEX with small set of generated of schedule as there is no feedback
instances and larger instances were not information which can update the respective
solved. higher level of plan. Without integration and
In their study Gartner and Kolisch [9] feedback mechanism, the decision made at
addressed the problem of planning the patient different levels are not coherent and result in
flow in hospitals subject to scarce medical local optimal solutions. Planning decisions
resources with an objective to maximize the made at higher levels cannot incorporate the
contribution margin. Patients were classified day by day changes and particularly, seasonal
on their diagnosis related group based on changes. Therefore, an integrated approach
their clinical pathway and decision was made combining all levels of decision making is
on which day each procedure of clinical highly needed in healthcare and present study
pathway should be done. Clinical pathway serves the purpose of developing a method
defines the procedures such as different types for integration of higher, medium and lower
of diagnostic tests and surgery as well as the level planning and making the planning and
sequence in which they must be applied to scheduling decisions more realistic and
the patients. Two mixed integer programs reflective of day by day dynamic situations in
were developed and embedded in a static and hospitals.
rolling horizon approach to cope with
stochastic data. Efficient planning and Moreover, Some of the studies used the
planning decisions which have been made at
scheduling of patients is significant for the
higher planning levels such as case mix , type
performance of healthcare facilities.
and capacity of clinical resources as well as
Therefore, lot of research is available in tactical decisions on master surgery schedule
literature focusing on planning and [9]. The studies conducted on operational
scheduling of operating rooms [10-22]. In level use the data from tactical level planning
literature, planning and scheduling of which is done far ahead of time. On medium
patients in operating rooms have been studied level planning most of researches use a cyclic
at three levels of decision making, including master surgery schedule [8,11,12,23,24] to
strategic level, tactical level and operational ensure stability of schedules for surgeons and
level planning problems. At strategic level, nurses. These schedules remain the same
the operating room capacity is divided among until decided by the managers based on
different type of patients and different blocks changing demand. These cyclic schedules
are created for each type of surgeries. While lead to some of the capacity overutilized and
45
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
long waiting lists and some of the capacity is the patients can visit the relevant departments
underutilized due to fixed blocks or rooms for and leave the hospital after the completion of
surgical specialties. Such a situation cannot treatment process as explained in below
lead to optimum utilization of critical Figure 1.
resources like operating rooms, ICU and Other Units
(Internal Medicine
recovery wards and also waiting lists keep on Elective Patients
Inpatients
Home
Surgical Ward
increasing in such a situation. To tackle this
Outpatients
problem, present study proposes a Intensive
Care Unit
current research first identify the bottleneck Figure 1. Flow of surgical patients in a
operating room based on the patient data and hospital
flow. Later, the respective model of operating
room planning and scheduling is solved using When the patients of different categories
optimization algorithm. Current research and different specialties visit the hospital, the
proposed a heuristic method to schedule patient treatment time or workload on various
patients in operating rooms with an aim to units of hospital is different and the
maximize the utilization of operating rooms department which have highest load as
keeping in view the patient constraints compared to the other departments becomes
related to their treatment requirements. bottleneck according to theory of constraints.
In proposed method, three planning levels,
2. PROBLEM DESCRIPTION Higher level planning, medium level
A scenario of hospital is assumed in order to planning and lower-level planning are
explain the nature of problem in this section. integrated. Planning and scheduling of three
Figure 1 shows the flow of surgical patients critical departments which have highest
in a hospital. A surgical patient is one who operational costs, including operating room
needs surgical procedures to be performed in (OT) Intensive care units (ICU), and recovery
his/her flow through the hospital. Patients wards are considered.
visiting the hospital are either elective For example, considering the flow of
patients or non-elective patients. Elective different patients from different departments
patients are those patients who are known as shown in the below Figure 2. Patients flow
beforehand and can be scheduled in surgery from three departments; cardiology,
lists. They include both inpatients and neurology and gynecology are shown in the
outpatients. Inpatients are those patients who Figure 2. All three patients need the surgical
spend at least one night at hospital and procedures to be performed as a part of their
outpatients visit the hospital for a few hours treatment as shown. In the Figure 2 for
and are not admitted to hospital. Non-elective example neurology patient visits the
patients, on the other hand are urgent or neurology department, then radiology
emergency patients who can visit the hospital department and operating room and intensive
any time depending upon level of urgency. care unit. Similarly, at the same time patients
Depending upon the treatment requirements from other departments also need surgeries in
46
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
operating rooms. Operating rooms are An open scheduling policy is adopted for
dedicated for different specialties by master OR.
surgery schedule at medium level planning as Operating room staff are available.
shown in the Figure 2 through different All operating rooms are considered as
colours. Due to load of patients, operating identical.
room may become bottleneck.
2.2 Mathematical Model:
Patients flow from three departments; A mixed integer programming model for
cardiology, neurology and gynaecology is scheduling of patients in operating rooms and
shown in Figure 2. their allocation in different planning horizons
is presented in this section. The objective
function for scheduling of patients is aimed
OPERATION THEATRE COMPLEX
Cardiology Patient to minimize the makespan of all operating
Neurology Patient
Intensive
shows the completion time of a surgery of a
Care Unit
q, q Q, oO 0 b B,
tqcomp tqset tqsur tqclr X qob q Q, oO 0 b B, (3)
47
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
X oqb 1 q Q, 0 b B, (8)
Advance Scheduling
Functions
Output Model solution
Optimum Plan and Iterations
Schedule Genetic Algorithm
Brach and Bound Allocation Scheduling
Branch and Cut
Constraints based
Heuristics
Resources
Requirement
Resources
Requirement
INPUTS
Real Time Patient Data
Inventory Patients Waiting List
Resources
Real Time Procedure Durations
Supply
Cancelled Surgeries
Hospital Available Capacity
Available Resources
Resources Emergency Cases
Supply Priority Rules
Constraints
Supplier Surgical Patients
Figure 3. An overview of the procedure adopted for multi-level planning and scheduling
start time. A medium level, operating room medium level plan to higher level planning
resources are configured based on the patient horizon to update the higher level plan.
Cardiology
Available
Operating
50
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
Proposed
Method
2
cloud and IOT technology to estimate the
8
3 288 values of variables using intelligent deep
9 7 2858 14 3024 18 learning algorithms and make intelligent plan
Table 3. Comparison of Average of patients on expected critical resources.
makespan values for all Problem
categories ACKNOWLEDGEMENTS
This work has been supported by the Higher
5. CONCLUSION Education Commission of Pakistan (HEC)
In the current research, a constraint heuristic under the National Research Program of
method is proposed for intelligent planning Universities (NRPU) Project No. 20-17024.
and scheduling of patients in operating
rooms. The proposed method used theory of REFERENCES
constraint and integrates strategic level, [1] Pacchini, A.P.T., et al., The degree of
tactical level and operational level planning readiness for the implementation of
51
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
52
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
53
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected]
ABSTRACT
Due to successive technological advancements, developments and innovations, the global
industrial landscape has drastically transformed over the last years. Industry 4.0 aims at
transforming traditional industries into intelligent ones by incorporating innovative technologies.
Industry 4.0 enables physical assets to be integrated into intertwined digital and physical
processes thus creating smart factories and intelligent manufacturing environments. Internet of
Things (IoT) is a rapidly growing technology that has drastically contributed to the Industry 4.0
realization. IoT pursues to pervade our everyday environment and its objects, linking the physical
world to the digital world and allowing people and “things” to be connected anytime, anywhere,
with anything and anyone ideally using any network and service. This study scrutinizes literature
review on IoT in the Industry 4.0 context. More specifically, it describes the IoT concept, its
security and privacy challenges and explores numerous IoT application domains. Moreover, it
presents and analyses the concept of Industry 4.0 and the benefits it offers as well as the relevant
key technologies (e.g. industrial internet of things (IIoT), cyber-physical systems (CPSs), cloud
computing, big data and advanced data analytics). Furthermore, it describes the emerging
application of IoT in Industry 4.0 and highlights the IoT and Industry 4.0 research gaps. Finally,
the need for innovation in the industrial domain and the impact and benefits that IoT and Industry
4.0 provide to everyday life and industries are described.
Keywords: Internet of Things (IoT), Industry 4.0, IoT applications
boards. Wearable devices will come with transportation, smart urban lightening, tele-
preassembled hardware application [7]. The care and others. IoT applications aims at
scope of wearable devices is limited to a improving the life quality for the end-user
particular application. For example, a smart community and supporting infrastructure
watch may be designed for call receiving and general-purpose operations. A survey
and fitness monitoring only. The advantage has also been conducted in which IoT
with wearable device is handy to use. The applications has categorized into the
examples of wearable devices are: Samsung following domains [11], [12]:
gear, google glass, digitsole smart shoe and
others. On the other hand, embedded system Transportation and logistics domain:
boards are open for the user. That is, a user Applications of this domain involve
can program these boards depending on the logistics, assisted driving, mobile
application he want. In present day market, ticketing, environment monitoring,
there are so many embedded system boards augmented maps etc.
are available. Examples include Arduino Healthcare domain: Applications of this
boards, raspberry pi boards, Samsung artik domain involve tracking, identification,
board, cloudbit, ESP8866 Wi-Fi module and authentication, data collection, sensing
others. After selecting a particular etc.
embedded system board, the user has to Smart environment domain: Applications
choose application software. There are wide of this domain involve comfortable
range of programming languages like C, homes/offices, industrial plants, smart
C++, python etc. IoT deals with connecting museum and gym etc.
the devices to a network and sharing the Personal and social domain: Applications
information, in this point of view one has to of this domain involve social networking,
ensure security and reliability. Also, IoT historical queries, losses and thefts etc.
deals with various kinds of devices, so in Futuristic domain: Applications of this
order to perform the task properly, domain involve robot taxi, city
standardization is required. information model, enhanced game
Standardization provides interoperability rooms etc.
and compatibility on a global scale [8].
Many professional bodies are working on 2.2 IoT Applications
creating standards for IoT. At present A number of IoT applications are presented
organizations like IEEE, International below [12], [13], [14], [15], [16], [17], [18],
Electro-technical Standardization, China [19], [20], [210]:
Electronics
Standardization Institute, American 1) Healthcare and sanitary
National Standards Institute are working on The medical domain is one of the first
the development of standards for IoT [9]. industries which adopted IoT and was
By designing IoT applications on accepted heavily influenced. IoT technologies can
standards, one can achieve reliability and create new opportunities, services and
interoperability, which in-turn will reduce applications to improve the healthcare and
maintenance costs in the long run [10]. The sanitary domain. Moreover, IoT platforms
main aspects of the smart city project are and services enhance current living
waste management, intelligent solutions and facilitate the realization of
56
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
the current environmental safety by allow for better energy efficiency and reduce
reducing and mitigating the potential impact the overall system downtime.
of damage and disaster. IoT technologies
allow the development of innovative real- 3. SECURITY AND PRIVACY ISSUES
time monitoring and decision-making IN IOT
support systems and applications regarding IoT has four main layers which are
environmental issues, such as early perception layer, network layer, transport
prediction and detection of natural disasters, layer (Middle-ware Layer) and application
weather conditions etc. layer. All IoT layers have their own privacy
and security concerns as conferred in Fig. 1.
4) Transportation and Logistics
Nowadays, vehicles as well as roads and 3.1 IoT Perception Layer Security and
transported goods, are equipped with more Privacy
sophisticated technological devices such as Perception layer contains groups of
on near field communication (NFC) tags, information. It is classified in to two main
radiofrequency identification (RFID) tags, sections which are perception node and
actuators, sensors etc. IoT technologies can perception network.
be used to enhance the potential of these Perception node is responsible for collecting
systems and optimize their use in the data and perception network handles the
domains of transportation, logistics and instructions of sending and managing data
suppliers, which are regarded as essential [24].
components in the productivity of many Perception layer is composed of many
industries. Intelligent transportation systems different sensor technologies like Radio
(ITS) are able to communicate, share and Frequency Identification (RFID) [25]. RFID
exchange mission-critical information and systems are exposed to privacy and security
data promptly, timely and accurately. problems. The perception layer checks
Hence, they are used to ensure that the collect and process data then transmits the
transportation network is efficiently information to the next layer which is the
monitored and controlled. IoT offers several network layer. Perception layer is
contemporary applications and services and
responsible for controlling data sources
in combination with the ubiquitous 5G
where IoT nodes are the main source of data.
mobile networks can provide industries with
The IoT nodes are widely vulnerable to
intelligent transportation and logistics
systems. These systems provide solutions attacks due to which a security and privacy
which are designed specifically for certain node in the perception layer scheme (SNPL)
needs and goals, thus accelerating is developed [26]. Some security and
productivity, profitability and operations. privacy issues of perception layer are
Moreover, they offer real-time monitoring discussed below:
and tracking throughout the entire supply
chain, thus helping enterprises increase end-
to-end visibility as well as maintain efficient
transportation control and cost-effective
management. In addition, they conduct more
effective route planning and optimization,
58
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
4) Spoofing:
Spoofing occurs when an attacker transmits
false and incorrect information to the RFID
system and try to make its originality falsely
and making it appear from the authenticated
and original source.
With the help of this, attackers get full
access to the system and make it vulnerable
[25].
5) RF Jamming:
Radio Frequency (RF) Jamming tries to not
comply with lower-level protocols to be able
to interfere with the ongoing legitimate
communication. This attack occurs when
Figure 1. Security and Privacy Issues in RFID tags are compromised by DoS attack
IoT that makes communication through RF
signals distributed with noise signals. The
1) Unauthorized Access to Tags: source that initiates jamming attacks could
Systems that have a large number of RFID be very powerful to damage the network or
face security issues because of the lack of it could have less power to only damage
proper authentication. Unauthorized users or small parts of the network [25].
hackers can access tags without
authorization, delete and even can modify All IoT layers are prone to security attacks
tags [25]. and threats. There are three main security
issues that are related to the perception layer.
2) Tag Cloning: First issue, is wireless signals strength [27].
Tags are distributed on different objects. In Perception layer, signals are sent and
Objects data can be viewed to be read and received to and from the sensors with the
modified with the help of some hacking
help of wireless technologies whose
techniques. This leads to tag cloning which
efficiency can be compromised by
occurs when tags can be captured by
disturbing waves. Second issue in IoT
criminals with ease who easily can make a
devices is that the sensor nodes can be
replica of tags and compromise it. [25]
stopped by the owner and the attackers due
3) Eavesdropping: to the reason of the external and outdoor
Eavesdropping is an interception of nature of IoT system could lead to physical
information between two nodes or attacks upon the IoT nodes and the IoT
communication devices [24]. The wireless system. The third issue is the nature of
characteristic of the RFID makes it easier for network topology. IoT nodes usually move
hackers to sniff out the confidential around many different places which means
the network topology is dynamic. The
59
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
Perception layer consists of RFID and security of information and enabling the
sensors. RFID and sensors storage, power network transmission [24]. It includes
consumptions, capacity and computation mobile devices, internet and cloud
capability are limited and this led to making computing. The network layer consists of
them prone to attacks. Perception layer Wireless Sensors Networks (WSN). This
security problems can be easily addressed layer takes care of transmitting data from the
with the help of pint-to-point or end-to-end sensors to their destinations with reliability
encryption. Perception Layer can provide [25]. The network layer is responsible for
various security features and it supplies four transmitting data to and from IoT devices
purposes which are privacy of data and and hubs and to serve data routing. In this
layer, technologies like Wi-Fi, Bluetooth,
sensitive information, authentication and
3G, LTE and Zigbee are used to operate the
risk assessment. Authentication is one of the
Internet, switching, routing and gateways
security goals that must be satisfied in
[27]. Network layer is composed of
systems in order to protect systems against
protocols, communication technologies with
hackers and attackers. Cryptography can be corresponding hardware and network [28].
used to apply authentication to systems. The network layer does an important job of
Cryptography has some algorithms that can connecting the IoT nodes and IoT
be used to provide a digital signature that applications together. The main threat that
could protect against attackers. Data need to threatens the network layer is the DoS attack
be protected and secured while collecting where the attackers make the service
and forwarding to the next layer. Symmetric unavailable for the legitimate users.
and asymmetric encryption algorithms can Network layer security and privacy issues
also be used to apply privacy to data. are listed below [25]:
Encryption algorithms are easy to be
1) Sybil Attack:
implemented in sensors due to their benefit
In a Sybil attack, the attacker works on
which is low power consumption. Location
attacking the system by manipulating the
anonymity and identity anonymity are must
node to have the single node more than one
to hide and secure sensitive information.
identity. This results in false information
This can be achieved by KAnonymity
[25].
approach that protects information like
identity, location and sensitive data of users.
2) Sinkhole Attack:
Risk assessment has an important role in IoT
Sinkhole attack works on trying to present
security because of its help in discovering
compromised nodes attractive to other close
new threats of systems. Also, it helps in
nodes. So, all data will flow from nodes to
defining security strategies that could be
compromised nodes which result in packets
classified to be the best. Also, it prevents
drop.
security breaches [25].
The system believes that the data have been
transmitted to the other side while system
3.2 IoT Network Layer Security and
traffic is silenced. Sinkhole attack can cause
Privacy
DoS attack due to more energy
The next layer after the perception layer is
consumption. Sinkhole attack process seems
the network layer. The network layer is the
to be unknown to the network where
layer that is responsible for providing
attackers deceive the system to make it
60
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
believe that all transmitted data is received heterogeneous which makes it impossible to
to the receiver [25]. use network protocols as it is. Protecting the
network’s objects has equal importance of
3) Sleep Deprivation Attack: protecting the network itself. The objects
Sleep Deprivation Attack works on the point should be able to have some actions to be
of keeping sensor nodes awake for a portion taken to from a guard to protect themselves
of time which leads to batteries consumption from attacks initiated against the network by
which in turns minimize batteries life time having the ability to know the network state.
which results in causing the sensor nodes to In order to achieve this, there must be good
shut down. This attack can keep the sensor protocols and software in the network that
node awake for some time. Energy help the objects to respond the situations and
constrained devices are prone to this attack. behaviours that are abnormal or affect the
objects and the network security. Security of
4) Denial of Service (DoS) Attack: Network layer is divided to three types
DoS attack occurs when an attacker works which are authentication, routing security
on enforcing the network to flood with a lot and data privacy. Implementing
of useless traffic which results in resources authentication and encryption could stop
exhausting of the system. So, the network of illegal accesses to nodes and this, in turn,
the system becomes unavailable to the users prevents spreading fake information.
[25]. Routing algorithms must be used to ensure
data privacy of data transmitted between the
5) Malicious code injection: sensors and the system. To improve the
Malicious code injection attack occurs when ability of the system to figure out errors and
an attacker tries to make a sensor node to protect the system against any kind of
insert some code that is malicious into the failure, the system has to provide multiple
system which in turns cause the network to paths for data routine. To monitor the system
shut down. Then, the attacker gets full and protect it against any kind of intrusion,
control over the network [25]. safety control mechanisms must be
implemented. To check whether data
6) Man-in-the-Middle Attack: received on an end is the same as the original
Man-in-the-middle attack is like a form of data sent from the other end, data integrity
eavesdropping attack. In Man-in-the-middle methods must be implemented [25].
attack the target is the communication
channel where unauthorized user can 3.3 IoT Transport (Middle-ware) Layer
monitor and control the communication Security and Privacy
between other two parties. Also, the The next layer after the network layer in IoT
unauthorized user can impersonate the systems is the Transport (Middle-ware)
identity of the victim and then Layer. Transport (Middle-ware) Layer
communicates through the channel to gain consists of data storage technologies like
information [25]. cloud computing. It is split into three layers
IoT introduces machine-to-machine which are local area, core network and
communication which has compatibility access network. The security problems of
security issue. In machine-to-machine the transport layer are classified as follows
communication network components are [25]:
61
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
62
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
CPS-enabled systems contain “cyber twined the main benefits of cloud computing.
services such as control algorithms and Moreover, users and enterprises have
computational capacities” [40] along with prompt access to applications, programs and
specialized computational capabilities, services which are rapidly provisioned with
physical assets and networked interactions minimal management effort and are hosted
and involve a large number of trans in the “Cloud” at any time and from any
disciplinary methodologies. The CPSs place. Hence, enterprises in the industrial
concept facilitate an ecosystem of cyber domain use various cloud-based
manufacturing, where smart machines applications widely in order to enhance the
process production data through a wireless crucial for their effective function systems
embedded network system. Moreover, CPSs
such as Customer Relationship Management
are defined as transformative technologies
(CRM), Human Resource Management
that can seamlessly link the physical with
(HRM) etc. Furthermore, enterprises that
the virtual world through their advanced and
utilize cloud computing can avoid the
novel systems [44]. Hence, they are
designed and developed to have both complexity of owning and maintaining their
physical inputs and outputs so as to enable own Information
and enhance the interaction with humans Technology (IT) infrastructure and the up-
using innovative modalities [45]. front costs which can be accomplished by a
“pay-asyou-go” method allowing the
3) Cloud computing enterprises to start small and invest into
Cloud computing or simply “Cloud” plays a more resources if there is more service
leading role in enhancing and transforming demand [48]. Some further benefits and
the current industry as it is a kind of significant advantages include high
outsourcing that combines large numbers of availability and maintainability, consistent
computer servers and resources with a view accessibility to data and services from any
to offering computer programs, high-level connected device as well as reduced
services and resources on an on-demand or development cost and product time-to-
pay per-cycle basis in real time. market. Taking into consideration the
Cloud computing is “a set of networks continuously evolving nature, the numerous
enabled services, providing scalable, application domains and multiple benefits of
guaranteed, normally personalized, cloud computing, a lot of profits can be
inexpensive computing infrastructure on yielded in the industrial domain, thus, more
demand, which could be accessed in a and more enterprises of different sizes and
simple and pervasive way” [46]. Cloud types are rapidly adopting this advanced
computing is divided into three levels of technology with the aim of enhancing their
service offerings, namely Software as a capabilities and capacity at a minimum cost
Service (SaaS), Platform as a Service (PaaS) [49].
and Infrastructure as a Service (IaaS) which
support different levels of virtualization and 4) Big data
management of the solution stack [47]. The The digitalization of everyday life through
use of advanced applications and services the adoption of smart devices and advanced
that dynamically scale with the increased technologies (e.g. IoT, Artificial
number of users is considered to be one of Intelligence (AI), Social Networks (SNS)
66
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
etc.) has led to the increase of data sources 5. EMERGING APPLICATION OF IOT
and the diversity of digital content as well as IN INDUSTRY 4.0
data types, forms and structures [50]. Manufacturing comprises of a core
Consequently, an enormous volume of industrial component which has a vital
heterogeneous data, named big data, is impact on people’s livelihood and a nation’s
generated and increases exponentially on a economy. Additionally, it is one of the
daily basis. Volume, variety, veracity, largest and highly interconnected IoT
velocity and value are the key factors which markets and it involves a wide variety of
characterize and differentiate big data from operations, processes, services, products etc.
traditional data. Big data plays a key role in With the aim of enhancing the overall
industries and intelligent manufacturing as it production, productivity and product quality
can provide enterprises with numerous management (PQM) throughout the various
advantages, merits and benefits through stages of the lifecycle of products, IoT offers
various predictive and prescriptive insights. applications and services under the umbrella
Hence, enterprises, which want to remain of Industry 4.0 which includes advanced
competitive, should give priority to the monitoring and tracking, performance and
implementation and utilization of maintainability optimization and human
contemporary advanced analytical tools, machine interaction. Hence, it stands for a
techniques, methods and applications with reason that IoT can provide a lot of solutions
the aim of processing big data, gleaning for Industry 4.0 in the manufacturing
intelligence and retrieving the value of the domain which is characterized by its
vital data in each case. These tools are complexity and breadth of applications with
named big data analytics (BDA) and use its diverse CPSs and its manufacturing
parallel and analytic techniques to analyse operation management (MOM)
huge volume of diverse, rapidly methodologies [49]. Internet of things (IoT)
transforming data enabling, thus, the is a pervasive technology and now it has
collection, the process and the management numerous applications in every technology
of vital information and statistics [51]. By sector. The fourth industrial revolution has
far, the most effective way for enterprises to also been evolved through the connected
gain immense benefits over their technologies enabled by the IoT. In fact, IoT
competitors, optimize operations, enhance provide a connected framework for large
productivity, quality and efficiency and scale manufacturing and production which
reduce operational costs is to use all the is very efficient, fast, cost effective, and free
newly gained knowledge in order to from the traditional manufacturing faults. In
generate invaluable insights and improve Industry 4.0, IoT apprehends the basic
equipment service and maintainability. concept of intelligent manufacturing 4.0,
Nonetheless, in order to fully utilize big data also known as smart manufacturing that
and exploit all of its benefits, enterprises transform the overall dilemma of every
must change their decision-making culture sector, it uses Service-Oriented Architecture
and take into consideration that no matter (SOA) and is considered to be a novel
how much the potential of big data and manufacturing model that takes advantage
analytic tools increase, the need for human of and fully utilizes various advanced
insight should not be overlooked [52]. information and manufacturing techniques,
methodologies and technologies using IoT
[49]. It aims at fundamentally transforming
67
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
IoT. It allows people and “things” to be implement and adapt to Industry 4.0 and IoT
connected anytime, anywhere, with will reap many benefits and profits and will
anything and anyone ideally using any be able to stay ahead of market competition.
network and service. Additionally, IoT Nonetheless, IoT in Industy 4.0 are still at an
pursues to pervade our everyday early stage of development, adoption and
environment and its objects, linking the implementation, therefore there are still
physical to the digital world. IoT which is various open issues and challenges that need
regarded as a dynamic and global network of to be addressed.
uniquely addressable interconnected
“things” aims at implementing autonomous, 8. CONCLUSION
robust and secure connections. Moreover, IoT is an innovative and rapidly growing
IoT provides various applications, functions technology which offers various novel
and services in everyday life and in a wide applications, services and solutions and
range of markets and industries. By links the physical to the digital world. It also
implementing and adapting to Industry 4.0 allows people and “things” to be connected
and IoT technologies, unprecedented levels anytime, anywhere, with anything and with
of economic growth and productivity anyone ideally using any path/network and
efficiency can be achieved by enterprises, any service. Moreover, it improves the
such as: quality for the end-user community and our
lives in general and supports infrastructure
Development of production systems and general-purpose operations. In addition,
which are characterized by interoperability, it aims at transforming the current industries
flexibility, adaptability, agility and into intelligent ones utilizing the dynamic
proactivity network of interconnected devices.
Optimization and improvement of Enhancing their operation and functionality,
efficiency, speed and quality particularly in increasing their productivity and reducing
engineering, operation, administration and their costs and waste are some of the many
decision-making benefits and profits that enterprises can gain
Enhancement of overall application, by using IoT. Moreover, enterprises that
services and system availability and fully adopt IoT in Industry 4.0 will be ahead
maintainability of their competitors, become more agile,
Acceleration of productivity and adapt to the continuously changing market,
reduction of lead time resulting in create products of higher quality that satisfy
decreasing timeto-market customers’ needs and requirements.
Facilitation of the adaptation to Moreover, in the context of Industry 4.0, IoT
individualized customer requirements and can be utilized in combination with other
market demands innovative technologies such as big data,
cloud computing, CPSs etc. in order to
IoT in Industry 4.0 can provide a enhance and transform the current
multitude of contemporary solutions, manufacturing systems into intelligent ones.
applications and services and can yield Industry 4.0 allows for machines to become
significant personal, professional and independent entities that are able to collect
economic opportunities. Hence, enterprises and analyse data and give advice upon it
and industries that are able to fully without requiring any human intervention as
it introduces self-maintainability, self-
71
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
72
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
technologies and challenges. IEEE [20] Perera, C., Liu, C. H., & Jayawardena,
Access, 6, 3619- 3647 S. (2015). The emerging internet of
[12] Atzori, L., Iera, A., & Morabito, G. things marketplace from an industrial
(2010). The internet of things: A perspective: A survey. IEEE
survey. Computer networks, 54(15), Transactions on Emerging Topics in
2787-2805. Computing, 3(4), 585-598.
[13] Sundmaeker, H., Guillemin, P., Friess, [21] Gilchrist, A. (2016). Industry 4.0: the
P., & Woelfflé, S. (2010). Vision and industrial internet of things. New York,
challenges for realising the Internet of NY: Apress., ISBN: 1484220463.
Things. Cluster of European Research [22] Lee, C. P., & Shim, J. P. (2009).
Projects on the Internet of Things, Ubiquitous healthcare: Radio
European Commision, 3(3), 34-36. frequency identification (RFID) in
[14] Vermesan, O., Friess, P., Guillemin, P., hospitals. In Handbook of research on
Gusmeroli, S., Sundmaeker, H., Bassi, distributed medical informatics and
A. & Doody, P. (2011). Internet of ehealth (pp. 273-281). IGI Global.
things strategic research roadmap. [23] Weiser, M., Gold, R., & Brown, J. S.
Internet of things-global technological (1999). The origins of ubiquitous
and societal trends, 1(2011), 952. computing research at PARC in the late
[15] Khan, R., Khan, S. U., Zaheer, R., & 1980s. IBM systems journal, 38(4),
Khan, S. (2012), Future internet: the 693-696.
internet of things architecture, possible [24] F. A. Alabaa, M. Othmana, I. A. T.
applications and key challenges. In Hashema, and F. Alotaibi (2017),
2012 10th international conference on “Internet of Things security: A survey,”
frontiers of information technology Journal of Network and Computer
(pp. 257-260). IEEE. Applications, vol. 88, pp. 10–28.
[16] Miorandi, D., Sicari, S., De Pellegrini, [25] M.U. Farooq, M. Waseem, A. Khairi,
F., & Chlamtac, I. (2012). Internet of and S. Mazhar (2015), “A Critical
things: Vision, applications and Analysis on the Security Concerns of
research challenges. Ad hoc networks, Internet of Things (IoT),” International
10(7), 1497-1516. Journal of Computer Applications
[17] Gubbi, J., Buyya, R., Marusic, S., & (0975 8887), vol. 111, no. 7, pp.
Palaniswami, M. (2013). Internet of 1-6.
Things (IoT): A vision, architectural [26] Y. Fan, G. Zhao, K. Li, B. Zhang, G.
elements, and future directions. Future Tan, X. Sun, and F. Xia (2020),
generation computer systems, 29(7), “SNPL: One Scheme of Securing
1645-1660. Nodes in IoT Perception Layer,”
[18] Borgia, E. (2014). The Internet of Sensors, ol. 20, no. 4.
Things vision: Key features, [27] R. Mahmoud, T. Yousuf, F. Aloul, and
applications and open issues. Computer I. Zualkernan (2015), “Internet of
Communications, 54, 1-31. Things (IoT) Security: Current Status,
[19] Zanella, A., Bui, N., Castellani, A., Challenges and
Vangelista, L., & Zorzi, M. (2014). Prospective Measures,” 10th
Internet of things for smart cities. IEEE International Conference for Internet
Internet of Things journal, 1(1), 22-32.
73
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
Technology and Secured Transactions [36] Geissbauer, R., Vedso, J., & Schrauf, S.
(ICITST), pp. 336-341. (2016). Industry 4.0: Building the
[28] H. Shin, H. K. Lee, H. Cha, S. W. Heo, digital enterprise. Retrieved from:
and https://www.pwc.com/gx/en/industries/i
H. Kim (2019), “IoT Security Issues and ndustr ies-4.0/landing-page/industry-
Light Weight Block Cipher”, 4.0-buildingyour-digitalenterprise-april-
International Conference on Artificial 2016.pdf.
Intelligence in Information and [37] Zhou, K., Liu, T., & Zhou, L. (2015).
Communication (ICAIIC), pp. 381-384. Industry 4.0: Towards future industrial
[29] K. Zhou, T. Liu, L. Zhou (2016), opportunities and challenges. In 2015
Industry 4.0: Towards Future Industrial 12th International Conference on
Opportunities and Challenges, Fuzzy Systems and Knowledge
International Conference on Fuzzy Discovery (FSKD) (pp. 2147-
Systems and Knowledge Discovery: 2152). IEEE
2147–2152. [38] Hajrizi, E. (2016). Smart solution for
[30] Schmidt R, et al (2015), Industry 4.0 - smart factory. IFAC-PapersOnLine,
Potentials for Creating Smart Products: 49(29), 1-5.
Empirical Research Results, [39] Lasi, H., Fettke, P., Kemper, H. G.,
International Conference on Business Feld, T., & Hoffmann, M. (2014).
Information Systems, pp. 16–27. Industry 4.0. Business & information
[31] Acemoglu D (2002), Technical systems engineering, 6(4), 239-242.
Change, Inequality, and the Labor [40] Zhong, R. Y., Xu, X., Klotz, E., &
Market, J. Econ. Lit., 40(1): 7–72. Newman, S. T. (2017). Intelligent
[32] Von Tunzelmann N (2003), Historical manufacturing in the context of
Coevolution of Governance and industry 4.0: a review. Engineering,
Technology in the Industrial 3(5), 616-630.
Revolutions, Struct. Chang. Econ. [41] Lampropoulos, G., Siakas,
Dyn., 14(4): 365–384. K., &
[33] Kagermann H, Wahlster W, Helbig J Anastasiadis, T. (2018). Internet of
(2013) Recommendations for Things (IoT) in Industry: Contemporary
Implementing the Strategic Initiative, Application Domains, Innovative
INDUSTRIE 4.0, München. Technologies and
[34] Pereira A, Romero F (2017), A review Intelligent Manufacturing. International
of the meanings and the implications of Journal of Advances in Scientific
the Industry 4.0 concept, Procedia Research and Engineering, 4(10), 109-
Manuf. 13(S): 1206– 1214. 118.
[35] Wee, D., Kelly, R., Cattel, J., & [42] Bi, Z., Da Xu, L., & Wang, C. (2014).
Breunig, M. (2015). Industry 4.0-how Internet of things for enterprise systems
to navigate digitization of the of modern manufacturing. IEEE
manufacturing sector. McKinsey & Transactions on industrial informatics,
Company, 58. 10(2), 1537-1546.
[43] Monostori, L. (2014). Cyber-physical
production systems: Roots,
74
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th & 7th March 2024, Karachi, Pakistan
75
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected]
ABSTRACT
A research has been conducted to figure out the impact of cutting speed over the deviation in the
thin-walls machined in the aerospace grade Aluminum 6061-T6. In this research, an experimental
setup has been designed to machine the thin-walls (ribs) of thickness ranging from 1 mm to 3 mm.
CNC milling machine is used to machine the thin-walls precisely using a constant tool path,
including both roughing and finishing operations. Two workpieces of Aluminum 6061-T6 have
been prepared upon which milling operation has been applied with the different cutting speeds. It
has been found that the thinner walls exhibit higher deviation as the cutting speed is increased,
whereas, the thicker walls have lower deviation, if the cutting speed is increased. It contradicts the
hypothesis in the previous research that the deviation in the thin walls can be decreased by using
the high-speed machining in the metals. Hence this research outlines a unique characteristic of the
Aluminum 6061-T6 grade, which exhibits a different pattern of deviation in the thin-walls cut
through the milling process.
Keywords: milling, aerospace grade Aluminum, T6061-T6, thin-walls, deviation.
1. INTRODUCTION
To reduce the mass of vehicle’s components, geometries will be intended to perform at
thin-walled Aluminum alloy components are high cutting speeds due to advantages of
used in wide variety of applications in having high materials removal rates and
various industries, most notably in the good surface finish.
aviation/aerospace sector. The ratio of the Rejection of parts and rework may result in
height to the thickness of the milling wall is high cost and delayed time of response to the
a typical characteristic in thin-walled strategic and defense-based projects, such as
structures, which is 15:1 [1] Thin-walled rockets, satellites and other products.
parts of aerospace grade Aluminum 6061-T6 In this research, various factors would be
and 7071 series are being used in industry. explored which leads to high distortion
During the machining of thin-walled alloys, potential of thin-walled structures and
most of the material is removed, due to which address a few measures which may perhaps
excessive heat is generated, resulting in large lead to reduction in the predisposition of the
deformations in the thin-walls. This is one of thin-walled parts to deform. A limited
the major concerns of machining of such literature review is available regarding the
class of materials, that is its vulnerability to machining of thin-walled parts; however, the
deform because of nature of convoluted research articles were reviewed critically and
structure. The machining of intricate hence presented below.
76
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
Due to ongoing advancements in the field of Based on the above literature review, it can
High-speed milling (HSM), currently be stated that the major challenge
available machine tools can accurately and encountered when machining of above
precisely process hard metals, difficult-to- components is distortion / bending,
cut materials, complicated 3D geometry, and dimensional error, surface finish, flatness,
micro/nano-features [2]. Schultz and parallelism, perpendicularity and clamping.
Moriwaki define the cutting speed for HSM To investigate the effects of controllables
of aluminum alloys in the range from 2500 to parameters of high-speed milling and the
10,000 m/min [3]. Izamshah et al. developed thin-wall thickness, comprehensive research
Finite Element Model to specifically predict has been conducted to explore the deflection
the distortions or displacements for in the thin-walls and the surface roughness.
machining thin-wall component made of By varying the controllable parameters, it has
Titanium [4]. Izamshah et al. experience the been found that the deflection and the
magnitude of wall deflection decreases as the roughness can be overcome.
value of tool’s helix angle increase [5]. Ding Hence the research article comprises of six
et al. conclude that higher cutting speed, sections, the first section presents the
smaller feed rates, and positive rack angles introduction, literature review and the
results in thinner deformation layers in problem statement of the work. The second
machined surface [6]. Using the carbide tool, section defines the design of research, the
Yang et al. accurately predict chattering and third section discusses the tool path planning
stability for peripheral milling of thin-walled for the machining, whereas the fourth section
Al 6061-T6 plate and workpiece with a outlines the measurement methodology of
curved surface. [7]. Sun et al. indicate that the sample workpieces through the CMM.
force-induced deformation has a significant The fifth section presents the discussion over
effect on chatter stability in milling thin-wall the measurements and analysis, and the final
parts [8]. Due to the low stiffness value, section concludes the research work.
Bolar and Joshi noticed increased deflection
at the free end during machining thin-wall 2. DESIGN OF EXPERIMENT
Al-7075-T6 as they expected during their There are three major controllable
simulation [9]. Kuczmaszewski et al. found parameters used for the machining, which
that Surface quality was low in the transverse are; spindle speed, feedrate and depth of cut,
rolling direction. In addition, the carbide tool among which spindle speed are considered to
produced higher roughness values than the be varied, whereas, the feedrate and depth of
PCD tool. They also noticed that the cut are kept constant for the machining in
maximum deflection was seen in the center Aluminum 6061-T6 alloy workpiece. The
of the workpiece, and high deflection was deflection in the thin-wall was measured
noted in the transverse rolling direction [10]. after the machining. The following table
Czyzycki et al. demonstrated that a high- presents the parameters to be set for the two
speed camera may be used to evaluate thin pieces to be cut with the levels.
wall deflections during milling, as evidenced The other parameters which remained
by comparisons to displacement unchanged during the experiments are
measurements using a laser sensor with an diameter of the tool, tool geometry, cutting
accuracy of up to 11% and modeling using condition (dry) and workpiece material (Al
the finite element approach with a 22 percent 6061-T6).
error [11].
77
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
the locations of 1 and 4 are deducted from it. Points Thicknesses of Thin Walls
It provides the deviation of point 1 and 4 with 1.0 1.2 1.5 1.7 2.0 2.5 3.0
respect to the 7. Hence the same method is 0 5 0 5 0 0 0
1 0.1 0.1 0.0 0.0 0.0 0.0 0.0
applied on each of the three grids of the rib. 25 28 96 61 62 79 11
The procedure is applied to all the seven ribs 2 0.1 0.1 0.0 0.0 0.0 0.0 0.0
on both the workpieces of the Aluminum. 10 27 85 61 59 85 12
Hence a total of 144 points are measured 3 0.0 0.1 0.1 0.0 0.0 0.0 0.0
76 01 01 78 60 76 08
through the CMM. The deviation of the six
4 0.1 0.3 0.2 0.1 0.1 0.1 0.0
points with respect to the three base points 98 17 14 73 69 84 18
are presented in Table 2 and 3 for both the Al 5 0.2 0.3 0.1 0.1 0.0 0.1 0.0
6061-T6 workpieces. The deviations of the 22 56 82 44 52 81 17
upper part of the rib (points 1-2-3) and the 6 0.2 0.2 0.2 0.1 0.1 0.1 0.0
35 55 48 88 43 89 17
medium part of the rib (points 4-5-6) are also Mean 0.1 0.2 0.1 0.1 0.0 0.1 0.0
plotted against the thicknesses of each rib for Deviation 61 14 54 18 91 32 14
the two workpieces are shown in Fig. 4 and
Fig. 5. Table 3. Deviation of Point 1 – 6 for
Workpiece 2
79
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
80
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
81
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
ABSTRACT
In all ages of human beings, water performs a very vital function for everyone. As is well known,
contaminated water can lead to several potentially fatal diseases. Water-related problems plague
many nations worldwide, the majority of which lack access to potable water. We have seen that
many villages are experiencing power outages, and middle-class people and tiny communities are
unable to purchase a water filter. Conventional energy has a finite supply and negatively impacts
the environment. In this study work, we use an unconventional approach that uses a solar-powered
water purifier. We need to take the purifying process to the next level so that everyone can afford
it. Direct current is produced by the solar panel and is kept in a battery. The filtration mechanism
that this purifier uses to remove undesirable germs and debris from water is powered entirely by
solar energy. By driving a pump that pushes water through a purifier's network of filters, pumps,
and hoses, the motor's rotation aids in the removal of pollutants. The method we employed in this
purifier was using solar panels to power the motor that lifted the water.
Keywords: RO membrane, Pre filters, Purifier, Solar panel, Small scale
system is sunshine [2]. The most reliable way stalks, and coarse particles, this kind of
to purify polluted water is RO filtration. The filtering is employed. A winnowing sieve is
RO system's semi-permeable membrane used to filter the contaminants out of the raw
removes excess minerals as well as other water. Since the sieve cannot filter small
soluble particles including bacteria, fungi, suspended particles in raw water, this
viruses, and algae. The device efficiently approach cannot be employed in extremely
removes particles as small as 0.0001 microns turbid or muddy conditions.
while turning the motor. In Pakistan, access
1.2 Filtration through cloth:
to clean drinking water is a big issue in both In the villages, the filter medium is made of
rural and urban areas [3]. There are several white, thin cotton fabric or old, abandoned
conventional techniques for purifying clothing. Uncooked water including particles
drinking water, including the use of chlorine of dust, germs, insects, plants, or mud detritus
pills, pots for good chlorination, slow and
can be removed using this filter. Only a
quicksand filters, and fluoride removal. limited amount of the suspended particles in
However, these techniques are more complex the water may be removed by this cleaning
to use and yield lower purity levels than RO. method. In Pakistani communities, the
In this study, we present a concept that uses a technique of cloth filtering is rather
solar panel to raise the water level in a tank widespread.
and supply power to a RO system for
purification [4]. In the event of an 1.3 Filtration through Clay
environmental issue or power outage, such as Vessels (Earthen Pot):
a flood or other disaster, the solar purifier's Occasionally, rather turbid water is filtered
battery stores energy, enabling the process to using clay jars with the proper particle size.
continue using solar power. It is a simple-to- Muddy water is collected in a large clay pot
assemble portable purifier that may be used or jar and let sit on the pot's lower surface,
in remote locations without electricity. This allowing the water to drip through the jar's
purifier operates without producing any porous clay wall. The trickling water is
pollutants [5]. combined in a clay pot by lowering it to the
Rural communities across the world have bottom of the porous clay container.
embraced basic, low-tech healing methods, Y. Zhang et al. [6] uses solar energy, which
with the specific goal of removing obvious is a free energy source, and stores it in
pollutants from water obtained from adjacent batteries. After that, the water is heated to a
resources. These conventional techniques are predetermined temperature (below boiling
quick and effective in eliminating some types points) using cheap heating coils. The
of particles from the water, but they don't filtering chalk allows the cold water to be
provide enough water to meet current further purified once it has condensed. At this
standards for drinking quality. This is a great point, the water returns to ambient
method for rural areas, and with a simple temperature by condensation. We can acquire
disinfection step, they can often provide pure drinking water using this procedure. We
water free of germs. There are several experienced nearly every phase of the
approaches. product development process, from obtaining
1.1 Filtration along Winnowing client requirements to completing the design.
Sieve M. Rizwan et al. [7] suggested that solar-
When the water source is contaminated by powered water filtration devices could be
airborne contaminants such as dried leaves, considered a significant way to generate
83
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
clean water. Solar energy is now a stable solar panel, charge controller, battery, and
energy source that produces no pollution. A inverter make up the power supply. A booster
solar-powered water purification system's pump, a reverse osmosis system, and a
design is entirely based on the thermal control circuit with a sensor, microcontroller,
technique, which harnesses the sun's heat- and relays make up the purification unit. A
converting energy to create heat. The process booster pump generates high pressure to
of absorbing heat to cause water to evaporate complete the reverse osmosis process. The
is the most important component. According microprocessor monitors the water tank's
to research, liquids are often heated using flat level and stops it from overflowing. We may
plate collectors, which generate heat at get pure water in the water tank by employing
comparatively modest temperatures (27°C to this method.
60°C). A solar-powered water purification The fundamental idea of reverse osmosis.
system is made up of a filter that eliminates
Solar panels are used to capture solar
impurities and a solar collector that collects radiation. A battery is then used to store this
sunlight to guarantee vaporization, the initial energy. A relay that operates on
stage of purification. There are now four electromagnetic principles connects the
distinct concepts in development. battery to the purifying unit. The reverse
R. K. Kumawat et al. [8] revealed from their osmosis system, water tank, and high-
study that reverse osmosis affected water pressure motor make up the purification unit.
purification. Reverse osmosis technology Reverse osmosis may be performed because
using raw water has advanced significantly in of the pressure created by the high pressure.
the creation of high-yielding, energy- The 8051 microprocessor monitors the water
efficient systems. The main objectives for level in the tank and stops it from
this field of study turned out to be a decrease overflowing. We can receive clean water
in energy consumption, an extension of from the water tank through this technique.
membrane life, and an increase in energy
2. METHODOLOGY
recovery. The main goal of the study is to As indicated in Fig. 1. Filters, solar panels,
optimize the variables in the pre-treatment booster pumps, solar charge controllers,
process of various water sources (ponds, reverse osmosis membranes, and other
canals, and surface waters) for reverse components are the main parts of this
osmosis plants. This will extend the purifier. This purifier's operation is entirely
membrane life of the raw water by lowering reliant on solar panel electricity. Batteries are
its solids content. To guarantee the greatest used to store energy, which is obtained from
possible decrease in total solids as well as in solar power.
several chemical parameters (BOD, TDS,
and bacteria), experiments were conducted.
Raw
These characteristics were determined to be Water Filters Booster RO
Membrane
Pump
ideal for feeding the pretreatment effluent Tank
84
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
85
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
Figure 3. Reverse Osmosis Process inches, or 1.7 square feet, which is similar to
A typical three-stage reverse osmosis system a standard medium-sized house mirror.
consists of three different kinds of filters. A
3.4 Solar Charger Controller
sediment filter is designed to capture bigger An electrical gadget called a solar charge
particles floating in water, such as rust and controller regulates the amount of electricity
mud. The carbon filter purges water of that the solar panel supplies to the battery.
VOCs, chlorine, and other minor impurities. This makes sure that energy is not
Reverse osmosis membrane that is semi- replenished during the day and that the solar
permeable and effectively eliminates almost panel receives electricity during the night to
all leftover impurities. prevent the deep cycle battery from being
depleted. While some charge controllers
come with other capabilities like load control
and illumination, power management is the
primary objective.
3.5 Storage Battery
It is a crucial component of the setup.
Seasonal variations in atmospheric
conditions result in variations in the intensity
of solar radiation. The impurity pump is
Figure 4. Inside of an RO Membrane shielded from overvoltage by the constant
Filter 12V voltage. It features a 12-volt output
3.2 Booster Pump voltage and a 5-amp output current.
Booster pumps are used to increase water 3.6 Sediment Filter
pressure. Usually, the osmotic pressure is Sand, grit, precipitated mineral particles,
higher. Purification requires that water flow insoluble iron oxide, and other debris are
from a high concentration to a low removed by the sediment cartridge, which
concentration. Therefore, to perform the prevents them from clogging the reverse
reverse osmosis process, the high- osmosis film surface or plugging the drain
concentration lateral pressure must be higher flow restrictor, which reduces the amount of
than the osmotic pressure. water produced. Sediment filters rated to
3.3 Solar Panel remove particles as small as 5 microns are
In this concept, solar energy is captured using used in the majority of RO systems. We
a solar panel. Photovoltaic cells are used to advise changing this filter at least once every
make solar panels. The power output is six months. It is necessary to change TINY
entirely dependent on the light's intensity. We RO sediment filters every three months.
are utilizing a 20-watt miniature solar panel 3.7 Carbon Filter
to charge a 12-volt battery that is exposed to City utilities disinfect drinking water to stop
sunlight. This portable solar panel is the growth of dangerous bacteria, viruses,
lightweight and has a wide range of uses. The and other microbes that can cause life-
panel is specifically made to charge tiny threatening diseases or even death.
batteries with a maximum capacity of 10,000 Nevertheless, over time, TFC films may
mAh or 10 Ah. The dimensions of little solar
degrade in these chlorinated water sources.
panels vary from 0.6 x 2.55 inches to 14 × 18
A basic chlorine level of 1.0 mg/L (1.0 ppm)
is required for water going through the plant
86
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
by the majority of city water utilities. building's roof, and gravity distributes the
Extended TDS passage (less impurity water to the lower stories. Basement pumps
rejection) may be allowed by the RO feed the tank with water from the mains
membranes after 1000 ppm-long regions of system, and when the tank is empty, the
free chlorine openness. In addition to pumps start up again. In this sense, the tower
eliminating chlorine, the carbon filter shields serves as a reservoir and a supply of water.
the membrane underneath it.
4. DESIGN CALCULATIONS
3.8 UV Lamp
UV lamp: A UV lamp is typically used in 4.1 Design solar parameters
Ultra Violet disinfection. To create clean, The following are the primary panel designs
safe drinking water, a UV light is employed. that have an impact on collector panel
A strong UV light known as UV-C, or performance.
germicidal UV, is employed in this Power for Application: The manifold
procedure. UV-C rays enter the body of the panel is directly impacted by the power
pathogens and assault them. The harmful needed for a certain application. In terms
organisms' DNA has been modified to of cost and received power, bigger panels
prevent growth and disease-causing have larger areas. This panel has an 18W
properties. UV treatment only adds energy to maximum power output.
the water; it makes no chemical changes. The angle of Latitude: To maximize solar
There is no removal of the sterilised bacteria energy absorption, the panel should be
from the water. These rays can render 99.9% kept at the angle of latitude. The angle
of the bacteria inactive. Unlike chemical between the earth's equator plane and the
disinfection, UV radiation does not cause live line connecting a specific place on its
things to build an immune system. surface to the equatorial centre is known
3.9 Flow Restrictor as the latitude of a certain location.
As the name suggests, a flow restrictor limits Using the following method, the
the flow of RO reject water. The flow latitude angle is determined.
restrictor not only regulates the purifier's Check off all three points. i.e.
recovery ratio but also generates back The Equatorial plane lines.
pressure on the RO membrane to maintain The Earth's centre.
high pressure within it, which is a The location on a globe map
prerequisite for RO purification. The high- where the panel is to be created.
pressure water entering from the booster
Join the location point with the earth
pump will flow freely from the RO
point.
membrane's reject exit if the flow restrictor is
Determine the angle formed by the
not employed. Low pressure inside the RO
earth's centre point and the equatorial plane
membrane and increased water waste are the
line.
results of this. If the pressure inside the RO
membrane is not high, the membrane will not Average Power Intensity: The average
function. power intensity for any locality is the unit
of power collected per day.
3.10 Water Tank
The primary purpose of a water tank is to 𝐴𝑣𝑒𝑟𝑎𝑔𝑒 𝑔𝑙𝑜𝑏𝑎𝑙 𝑟𝑎𝑑𝑖𝑎𝑡𝑖𝑜𝑛 = 6.1 𝐾𝑊ℎ⁄𝑚2. 𝑑𝑎𝑦 (1)
hold water. Another name for it is the storage But,
tank. Water is stored in a tank on the 1 𝐾𝑊ℎ⁄𝑚2 = 3.6 𝑀𝐽⁄𝑚2
87
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
88
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
90
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected]
ABSTRACT
Closed and long air ventilation ducts (HVAC systems) become polluted by different contaminants
and particulate matters while functioning for larger period of time. Ventilation duct cleaning
process is difficult and expensive with manpower and usually takes a longer time. It can also
cause respiratory and other health problems for labor in cleaning and investigation. This paper
aims to propose an air duct cleaning robot (ADCR) for a duct which is low powered Arduino
operated and sensory based robust controlled system. It is an upgradation of some cleaning and
brushing mechanisms in terms of design. Further proposed design is being evaluated in terms of
effective suction and time taken analysis by robot is also done as validation for sustainable
solution. The work highlights mechanism as light, inexpensive and easily works in a closed
environment.
Keywords: Autonomous Vehicle, Dry Duct Cleaning, Arduino controlled Cleaning Robot
process, this paper suggests an autonomous
1. INTRODUCTION duct cleaning robot that is capable of
The purpose of an air duct ventilation system fulfilling the purpose of duct cleaning using
is to provide fresh air into confined spaces appropriate dry-cleaning method and
used in industries, subway stations, and required less manpower and cost effective
offices where people do work and spend most [6]. This method includes vacuum fans and
of their time [1]. This air duct system controls filters which suck the small particles from the
various air flows, i.e., exhaust air, outlet air, duct. Mechanical brush is used to displace
return air, and supply air. It consists of sticky particles. Sensor controlled movement
mechanical parts like dampers, filters, and air which senses the path of the duct and
terminals [2]. Filters have been installed performing turning action within the duct [7].
which will trap some dust particles initially, Arduino based circuit and two gear motors
but they cannot filter out all particulate matter used for the low RPM movement. Only
at the outlet of the air duct [3]. Therefore, horizontal surfaces can be cleaned from the
these dust particles accumulated at the inner duct. It can be made for vertical, round and
surface of ventilation system [4]. To provide other wall pressed cleaning robot. This type
clean and fresh air through these air of robot used in industrial HVAC system and
ventilation duct, it is mandatory to remove domestic level duct cleaning system to
these accumulated particles. Ventilation duct provide healthy and fresh air to the
cleaning process is difficult and sometimes environment [8]
need to replace with new duct which is The characterization of cleaning methods
expensive [5]. In the study of duct cleaning involves dry method and wet method. Both
91
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
92
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
93
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
94
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
3.3 ADCR Speed Control The required response of the system is that
Calculations the distance covered by the Air Ventilation
In speed analysis, finding the theoretical Duct’s Cleaning Robot increases as the
speed of robot on duct using following current increases.
equations. As power can be found as [24]
P T ( ) (1)
Then acceleration will be
T ( )
a
mv (2)
The time taken to reach the speed ‘v’ is
v 1
t dv
0 a (3)
v mv
t
v v
0
T (1 )
r r
Figure 5. Graph of Step Response of DC
(4)
T t
Gear Motor
v(t ) v(1 e r mv ) (5) The closed-loop response yields a more
Where v=ω * r , precise and accurate simulation when
Given data, t=1 sec, ω= 15.7rad/s, motor compared with the open-loop response. The
torque=T= 0.04 Nm, Radius of tire= r = 0.045 transfer function in case of the closed-loop is
m, Mass of robot = 3.5 kg, inserting values in given as:
equation (5), the speed of robot V
(theoretical) was found to be 0.213 m/s. It is 0.044
meant that robot cover 0.213-meter distance 9.45e5 s 2 6.5e2 s 0.1319 (8)
in 1 sec. Experimentally V(actual) was found
to be 0.186 m/s, hence efficiency of robot 4. RESULTS AND DISCUSSIONS
speed will be [25]
η=
Actual Speed
100
4.1 Performance Evaluation On The
Theoretical speed Basis Of Dust Collection
(6) Suction is tested on duct prototype of (8 *1 *
0.186 2) feet (length * height * width) in dimension.
η= 100
0.213 Table 1 summarized the results before and
η=87.32%
after cleaning. A common vacuum cleaner
This is efficiency of air duct cleaning robot operate at the range of 50-100 cfm and this
and it is appropriate speed for dust collection. ADCR has vacuum suction of 60 cfm [22]. It
For a 1-radian step reference, following is in the range of other vacuum cleaner and it
controlled states are given below [26] can clean all dust particles inside the duct in
Settling time smaller than 0.040 seconds one forward movement [18]. Performance of
Overshoot smaller than 16% this suction was also verified by dust
Steady-state error is zero. collection evaluation which concludes the
The open loop response of the system is
shown below in the figure. This response is
generated in MATLAB®.
0.044
9.45e s 6.5e 2 s
5 2
(7)
95
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
The purpose was to get maximum cleaned It was addressed and claimed in [28] that
surface. Two tests were carried out average one-time duct cleaning cost by
simultaneously. Dust particles collected has applying labour is $200 which is very
net approximate weight of 4.11g which then expensive while air duct cleaning robot
undergoes cleaning operation. Finally in (ADCR) manufacturing cost is $180 only and
comparison of percentage dust reduction of this is one time cost and the cleaning function
dust particles was found to be 82.9%. ADCR can be performed multiple times.
is put again in duct for second test. Dust
collection is reduced, it was due to some 4.2 Time Consumption Analysis
sticky particles. It is evident that percentage ADCR moves with different velocities at
dust collection will be further reduced for different points. By using standard deviation
further tests and it may become steady. of all working velocities of robot. At these
Overall efficiency from the beginning of optimum velocities, the time taken by robot
operation is 85.40 % hence suction is proven for one operation/test is calculated below.
appropriate and sufficient for the sampled Initial speed = v1= 0.213 m/s
duct. In that way, the system may be Speed while dust collection in process = v2 =
enhanced to more complex duct and 0.186 m/s
dimensions. Speed at turn = v3 =0.160 m/s
The environmental benefit of Air Ventilation Speed where dust will maximum = v4 =
Duct’s Cleaning Robot is that it creates a 0.174 m/s
cleaner living environment by reducing
allergies and irritants which help everyone to (v-μ) 2
breathe easier and improve air flow Standard deviation = S.D=
N (9)
efficiency [27].
96
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
Mean=μ=
v1 +v 2 +v3 +v 4 These performance criteria’s lead the way
4 (10) towards the effective and reliable design
Mean=μ=0.1832m/s solution. The consideration of environment
Standard deviation = S.D = 0.0194 and social benefits in second paragraph of
section 4.1 has enhanced its advantages to
Vmini = 0.1636 m/s, Vmax = 0.202 m/s and
local community and has concluded the
size of duct prototype = 2.74 m strategy to be viable in future.
At minimum velocity:
𝑡𝑚𝑖𝑛𝑖 = 16.7 𝑠 ACKNOWLEDGEMENT
At maximum velocity: Moral support of Research mates, guidance
𝑡𝑚𝑎𝑥 = 13.5 𝑠 by the supervisor and research environment
of KFUEIT is acknowledged.
As velocities in the duct is varying Appendix A
throughout the duct. The time has been
calculated and can have any value between Operating Voltage 5V
the ranges of 13.5 to 16.7 seconds. Input Voltage
7-12V
(recommended):
5. CONCLUSION AND Input Voltage (limit) 6-20V
RECOMMENDATIONS
54 (of which 15
If a duct is being cleaned by labour, the dust Digital I/O Pins provide PWM
contaminants like asbestos can produce
output)
health hazards like lung scarring and fibrosis
which in later stages can cause cancer. Most Analog Input Pins 16
of the cases, irritating particles can produce DC Current per I/O Pin 20 mA
disfunctioning of mucous membranes of the DC Current for 3.3V
50 mA
nose and throat. Though covering of face by Pin
mask can reduce the risk but still the particles 256 KB of which 8
can attach to PPE if safety measures had not Flash Memory KB used by
taken during cleaning. It also can cause
bootloader
chronic pulmonary disease like emphysema
and bronchitis. In order to facilitate the SRAM 8 KB
cleaning of duct, an autonomous air duct EEPROM 4 KB
cleaning robot (ADCR) is designed, which Clock Speed 16 MHz
consists of sensory based integration of parts. LED_BUILTIN 13
This autonomous strategy can reduce human Length 101.52 mm
effort. The ADCR dust collection and time
Width 53.3 mm
taken analysis of mechanism is the necessary
performance evaluation criteria which is an Weight 37 g
additional ingredient to the previous work. Table 4. Microcontroller ATmega2560
Amount of dust is significantly reduced to 85 Specifications
percent in second test. Multiple tests can be
performed in order to increase the percentage Appendix B
dust collected. It has been noticed that one
operation/test can take up to approximately
16.7 seconds. Time taken by mechanism to
complete one operation/test was performed.
97
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
Voltage DC 5V
Static current Less than 2mA
Level output High-5V
Level output The end of 0V Duct Cleaning Robot 1,” vol. 7, pp.
Sensor angle 15 degrees 26–32, 2014.
Detection distance 2cm-450cm [5] P. Gloriani, N. Rahman, and D.
High precision Up to 0.3 cm Caldwell, “ScienceDirect
Connection VCC, trig(control ScienceDirect ScienceDirect models
side), for capacity in Industry Novel
echo(receiving Integrated Robotic System for Tiny
end), GND Duct Inspection,” Procedia Manuf.,
Table 5. Specification of Ultrasonic vol. 17, pp. 342–349, 2018.
Sensor [6] S. G. Faal, “Design and Analysis of a
Robotic Duct Cleaning System,” 2011.
Appendix C [7] R. Seals, “Mobile robotics,” Electron.
Power, vol. 30, no. 7, p. 543, 2010.
Operating voltage Between 3 V and 9 [8] A. Bulgakov and D. Sayfeddine, “Air
V Conditioning Ducts Inspection and
Nominal voltage 6V Cleaning Using Telerobotics,”
The free-run speed 281 RPM Procedia Eng., vol. 164, no. June, pp.
at 6 V 121–126, 2016.
Free-run current at 6 80 mA [9] W. Jeong, S. W. Jeon, D. Park, and S.
V B. Kwon, “Force control of a duct
Stall current at 6V 900 mA cleaning robot brush using a
Stall torque at 6V 4 kg.cm compliance device,” ICINCO 2012 -
Gear ratio 1:21 Proc. 9th Int. Conf. Informatics
Weight 84 g Control. Autom. Robot., vol. 2, pp.
Table 6. Specifications of JGA25-370 DC 372–376, 2012.
Gear motor [10] I. Survey and I. Survey, “CHAPTER
2,” no. 1992, pp. 9–21, 1997.
REFERENCES [11] D. Version, “Removal of Particles
from the Supply Air of Ventilation
[1] W. Wu, H. Skye, and P. Domanski,
Systems Avoiding the Formation of
“Selecting HVAC Systems to Achieve
Sensory Pollution Source Delivery of
Comfortable and Cost-effective
Clean Air to Building Occupants D o c
Residential Net-Zero Energy
t o r a l d i s s e r t a t i o n,” 2008.
Buildings,” Appl. Energy, vol. 212, pp.
[12] M. A. Koledoye, “Design of a Mobile
577–591, Feb. 2018.
Robot for Air Duct Exploration,”
[2] J. Shim and R. C. Arkin, “Robot
2017.
Deception and Squirrel Behavior: A
[13] P. Shukla, “Design of Inspection and
Case Study in Bio-inspired Robotics,”
Cleaning Robot,” Int. J. Sci. Res. Eng.
2014.
Technol., vol. 3, no. 6, pp. 2278–882,
[3] CDC, “Appendix B : HVAC Systems
2014.
and Indoor Air Quality,” Cdc, 2014.
[14] W. Ya and Z. Jianhua, “Autonomous
[4] W. Jeong and S. Researcher,
Air Duct Cleaning Robot System,”
“Performance Analysis of a Mobile
98
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 4th-5th March 2024, Karachi, Pakistan
undefined, vol. 1, pp. 510–513, 2006. Adv. Mech. Eng., vol. 10, no. 3, pp. 1–
[15] G. Mester, “Applications of Mobile 9, 2018.
Robots Unmanned Autonomous [25] M. S. Essers and T. H. J. Vaneker,
Systems View project Self-Driving “Developing concepts for improved
Cars View project,” 2006. efficiency of robot work preparation,”
[16] W. Jeong, S. W. Jeon, and D. Park, Procedia CIRP, vol. 7, pp. 515–520,
“Operational effects of a mobile robot 2013.
system for cleaning ventilation ducts,” [26] A. S. Semenov, V. M. Khubieva, and
ACM Int. Conf. Proceeding Ser., Y. S. Kharitonov, “Mathematical
2013. modeling of static and dynamic modes
[17] S. Pedre, M. Nitsche, and F. Pessacg, DC motors in software package
“Design of a Multi-purpose Low-Cost MATLAB,” 2018 Int. Russ. Autom.
Mobile Robot for Research and Conf. RusAutoCon 2018, pp. 1–5,
Education Design of a multi-purpose 2018.
low-cost mobile robot for research and [27] “DOHS Fact Sheet on HVAC Duct
education,” no. September, 2014. Cleaning DOHS Fact Sheet on HVAC
[18] Z. Bhuiyan, Design analysis of Dust Duct Cleaning,” vol. 04, no. 301, pp.
collection system. 2020. 4–6, 2015.
[19] “ASTM F1005 - 91(2013) Standard [28] “2020 Air Duct Cleaning Costs | Clean
Practice for HVAC Duct Shapes; AC, Furnace, HVAC Vents -
Identification and Description of HomeAdvisor.” .
Design Configuration.” .
[20] E. Pawar, “A Review Article on
Acrylic PMMA Eshwar Pawar,” ISOR
J. Mech. Civ. Eng., vol. 13, no. 2, pp.
1–4, 2016.
[21] V. Gajbhiye, N. Ahmad, and M. S.
Tufail, “Design and Application of
D.C. Vacuum Cleaner using Axial
Flow Fan,” Int. J. Eng. Tech., vol. 4,
no. 1, 2018.
[22] A. Viegand Maagøe and B. V. Van
Holsteijn en Kemna, “Review study on
vacuum cleaners - Final report,” no.
June, 2019.
[23] D. Karayel, Z. B. Barut, and A.
Özmerzi, “Mathematical modelling of
vacuum pressure on a precision
seeder,” Biosyst. Eng., vol. 87, no. 4,
pp. 437–444, 2004.
[24] X. Jiang, L. Chen, X. Xu, Y. Cai, Y.
Li, and W. Wang, “Analysis and
optimization of energy efficiency for
an electric vehicle with four
independent drive in-wheel motors,”
99
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected]
ABSTRACT
Aluminum and its alloys are widely adopted material in industries like aviation and avionics
because of its durability and strength as compared to their inherited lightness. These industries
require extreme precision for which nontraditional machining processes are employed. Wirecut
electrical discharge machining (WEDM) is one of the commonly used technique which gives
satisfactory results of required accuracy and precision. It is a complex process for optimization,
because of the involvement of multiple input and output variables. This research aims to optimize
the WEDM for aluminum targeting four major output variables, i.e. material removal rate (MMR),
surface roughness (SR), wire wear rate (WW) and kerf width (KW) on different sets of input
variables, namely cutting angle, pulse on-time (POT) and current. By using actual experimental
machining data, Response Surface Methodology (RSM) and Artificial Neural network (ANN) are
employed and compared, based on root mean square error, and coefficient of determination. The
results declare that both RSM & ANN give accurate results and can be used as a reliable means
for predictions, but ANN gives a little improved prediction than RSM.
Keywords: Wirecut Electric Discharge Machining (WEDM), Artificial Neural Network (ANN),
Response Surface Methodology (RSM), Root Mean Square Error (RMSE), Regression Value
(R2).
1. INTRODUCTION
In recent years, aluminum and aluminum better resistance to high temperatures, a not
alloys have been widely used in the least, a controlled coefficient of thermal
automotive industry. Their prominent expansion.[1]. For their precise machining,
properties include low density, have good electrical discharge machining (EDM) and
formability and formability, high resistance wire EDM are widely used methods
to corrosion, and high electrical and thermal worldwide. EDM, which is the most widely
conductivity. Due to their excellent used non-traditional material exclusion
combination of mechanical, physical and processes, is a non-conventional, non-
tribological properties compared to base contact type machining processes, where the
alloys, aluminum alloys are of enormous material is detached with the help of
industrial importance. These properties stimulated frequency sparks generated
include high specific strength, high between the tool and work piece [2]. Due to
resistance to wear and seizure, high stiffness, the absence of physical contact of tool and
100
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
work piece, this process can be used for any material removal rate (or sometimes the
type of material regardless of their hardness, cutting speed), surface finish and kerf
provided the material is electrically (cutting width), and wire wear are considered
conductive [3,4]. Its prominent feature is the as the prime variables to start with. The same
use of thermal energy to machine electrically are being considered in this research.
conductive parts irrespective of their The process to determine an optimum model
hardness, which proves to be its typical for a sustainable WEDM process requires
advantage in the manufacturing of different thorough knowledge of this process,
dies & molds and specifically to bring empirical equations based on realistic
innovation in the automotive, aerospace and constraints as well different mathematical
surgical components. The process is based optimization techniques [9]. The
interesting to study with having different selection of optimal combination of
variables to evaluate, which include a long parameters for the process is the job of the
list of electrical as well non-electrical process planner who incorporates his
parameters that have great potential for a experience to the model to get the sustainable
substantial amount of research. The optimum solution of the process. The experience is
selection of the EDM process parameters is achieved by experimentation which is a
critically important to obtain an energy and costly, time consuming and tedious approach
material efficient process design. Numerous for sustainability [10]. When we talk about
researchers carried out several investigations the real-life support for the EDM process in
for refining the process performance, and terms of sustainability, low energy
have been brining optimal solutions for many consumption and better productivity, the
of the process parameters and variables manufacturer must be provided with
[[5,6,7] workable mathematical models. The most
From the list of all parameters and variables common practice to get a sustainable
EDM offers, machining parameters comes in working model for WEDM is to employ the
first place by playing a crucial role in theoretical or empirical based models, which
achieving high precision or high tolerance are never up to the mark for obvious reasons
machining with the required quality of the [11]. The empirical models have been
other output responses[8] .The selection of formulated on the basis of specific, limited
machining parameters is critical, as altering experimental setups, for example the type of
one parameter may disturb the optimization objective function, constraints (linear,
model for this process.. Hence, there is an nonlinear) number and nature of variables
actual justification for the ongoing used (integer, continuous etc.), solution
experimentation and research that is striving space size, solution space structure (convex,
to support and generate a systematic non-convex etc.), are too not predictive and
mathematical approach supported by supportive of the actual WEDM process
artificial intelligence to obtain the optimum [11,12,13,14] Hence, the latest trends in
combination of multiple inputs and output research goes with incorporating the concept
parameters for different work piece and wire of machine learning and artificial
materials, in order to get better overall intelligence for the optimization modelling
optimum machining performance and of WEDM.
process efficiency in WEDM process. On Machine learning is becoming an
discussing the overall process efficiency, the increasingly popular concept in the modern
parameters with highest level of importance world since its common goal is to optimize
in terms of performance trials in WEDM, systems by allowing one to make smarter use
101
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
sustainability in terms of MRR and surface research on WEDM focusing of MRR and
roughness. Their research employed ANN as surface roughness only, and used water based
the modeling tool and found the results di electric which offers greater difficulty in
satisfactory. To overcome the unpredictive understanding the cutting parameters of
support of WEDM optimization models, the WEDM. Their results provided a good
latest trend of research goes into utilizing the prediction when it produced the average
artificial intelligence techniques because of relative error of 0.27% and 0.67% in MRR
their ability to overcome the non-linearity and surface roughness respectively.
between the input parameters and output D. Devarasiddappa et. al [31] developed an
response, and to effectively handle the ANN based model by using Box Behnken
incomplete data and encountering the experimental design for the prediction of
multiple variables. [24,25,26]. surface roughness of aerospace alloy Inconel
As discussed in the initial paragraphs above, 825 against few of the input parameters and
the popular empirical optimization the developed model showed an accuracy of
techniques have been employed are 93.62% with an average error of 6.38%.
Response surface methodology (RSM) and Girish Kant et. al [32] used SVR (support
Taguchi method, which are extensively used vector regression) model along with the
in modelling the process parameters with ANN (artificial neuro network) to evaluate
respect of optimization and sustainability the WEDM process for the optimal
[27,28,29,30]. All of these approaches target prediction of energy consumed in the
towards the optimization of the overall machining process An investigation on
process but are unable to provide a generic WEDM is also carried out by using neuro-
approach for the optimization problem while fuzzy inference system to get sustainable
having different sets of variables types, model which supports higher productivity by
objectives and constraints. Research shows focusing on the overall quality of the
that employment of artificial intelligence- machined product and less heat-affected
based models brings up to three times more zone [33]. A Conde [34] declares that the
accurate results than the one achieved by the prediction of the machined part accuracy of
statistical methods and hence are the ultimate WEDM process is too, highly unpredictive in
options to work on the sustainability based nature, and there have been very little efforts
parametric model for WEDM. Shandilya et in getting it modelled due to many of the
al. [26] observed this support in making phenomena involved in the WEDM process
energy focused model by employing AI which are not clear yet. He has presented an
techniques while machining of composites AI model based on Elman-based Layer
with WEDM. Devarasiddappa et al. [27] Recurrent Neural Network which supports
supports the same idea as he observed the the accuracy of machined part by WEDM,
values of surface roughness while machining based on the wear rate of tool wire. His
the aerospace alloys and employing AI model reduced the average deviation of the
techniques and getting the level of accuracy machined parts due to wire deformation by
as high as 93.62%. He also urged the 80%. Sandeep Kumar et.[35] al
researchers to further develop an ONLINE experimented on WEDM by using the
SYSTEMS based on artificial intelligence by titanium super alloy Ti–6Al–4V and
probably using ANN technique. employed the ANFIS (adapted neuro fuzzy
Muthukrishnasen and Davim [28] states that inference system) model and investigated the
ANN offers higher potential for predicting effect on MRR and surface roughness of the
non experimental patterns. They carried their machined product. Deepak Kumar et. al [36]
103
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
was used for the measurement of surface Where Δm is the change in mass of the
roughness (Ra). It has a range of measuring Al5454 work piece, ρ is its density, and Tm
up to 250 mm × 320 mm of workpiece with a is the machining time for the same. The
resolution of 0.033 µm. Surface roughness experimental design matrix is shown in table
was measured on each slot of workpiece at 1.
three different places,10%, 50% and 100% of
the height of work piece, for incorporating
the studying the varying thickness impact.
The mean of average surface roughness
measured with the mentioned criteria was
used as the final reading.
Wire wear (WW) was calculated by
subtracting the weight of wire after the
WEDM cutting from the initial wire weight
(measured before machining). An electronic
physical balance with the least count of
0.0001 was used for this measurement.
The metal removal rate was found out by
employing the equation number 3.
MRR=Δm/ρ⋅Tm (3)
Figure 1.
105
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
.
Material Angle Current POT SR MRR WW KW SR(mm) MRR(msec/mm3) WW(mg) KW(mm)
Aluminum 30 5 20 0.0245 16.8255 0.0732 0.430000 0.0265 17.0127 0.0732 0.4425
Aluminum 30 5 20 0.0285 17.1999 0.0732 0.455000
Aluminum 45 5 20 0.0215 12.9876 0.0742 0.407500 0.0275 13.1399 0.0827 0.4100
Aluminum 45 5 20 0.0335 13.2922 0.0912 0.412500
Aluminum 60 5 20 0.0325 8.7959 0.1483 0.362500 0.0308 9.0130 0.1221 0.3863
Aluminum 60 5 20 0.0290 9.2300 0.0958 0.410000
Aluminum 30 10 20 0.0395 27.5929 0.0703 0.565000 0.0368 25.9529 0.0768 0.5650
Aluminum 30 10 20 0.0340 24.3128 0.0834 0.565000
Aluminum 45 10 20 0.0530 18.1123 0.0872 0.382500 0.0525 18.0898 0.0910 0.3838
Aluminum 45 10 20 0.0520 18.0672 0.0947 0.385000
Aluminum 60 10 20 0.0555 15.4982 0.1566 0.387500 0.0478 16.1988 0.0817 0.4050
Aluminum 60 10 20 0.0400 16.8994 0.0067 0.422500
Aluminum 30 5 40 0.0250 19.4838 0.0900 0.430000 0.0188 14.9529 0.0898 0.4463
106
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
3.3 Comparison of the performance of parameters. The equations for these are listed
RSM and ANN modeling below as equations 8 and 9 respectively.
The results of RSM and ANN are compared
by calculating the coefficients of ∑𝑛 (𝑦𝑖−𝑦0𝑖)2
R2 = 1- ∑𝑛𝑖=1(𝑦0𝑖−𝑦𝑚)2 (8)
determination (R2) and the Root Mean 𝑖=1
1
RMSE =𝑛 ∑𝑛𝑖=1(𝑦𝑖 − 𝑦𝑎𝑖)2 (9)
General Linear Model: KW versus
Angle, Current, POT
where “n” is the total no. of experiments, Factor Type Levels Values
Angle fixed 3 30, 45, 60
“yi” is the value predicted from the trained Current fixed 2 5, 10
ANN model, “yai” is the actual value POT fixed 2 20, 40
provided to the system as training data, and Analysis of Variance for KW, using Adjusted SS for Tests
“ym” is the mean of the actual values Source DF Seq SS Adj SS Adj MS F P
obtained through the experimentation. Angle
Current
2
1
0.086515
0.005098
0.086515
0.005098
0.043258
0.005098
21.39
2.52
0.000
0.129
POT 1 0.004396 0.004396 0.004396 2.17 0.157
Error 19 0.038432 0.038432 0.002023
3.4 Results of ANOVA Total 23 0.134441
108
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
0.01
Unusual Observations
Residual
Obs Angle MRR Fit SE Fit Residual St Resid
0.00
13 30.0 19.484 15.246 0.882 4.238 2.31R
14 30.0 10.422 15.246 0.882 -4.824 -2.63R
-0.02
WW versus Angle, Current, POT 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.055
The regression equation is Fitted Value
Analysis of Variance
0.0
Source DF SS MS F P
Regression 3 0.001068 0.000356 0.29 0.829
Residual Error 20 0.024164 0.001208 -2.5
Total 23 0.025232
-5.0
Source DF Seq SS 5 10 15 20 25
Angle 1 0.000753 Fitted Value
Current 1 0.000028
POT 1 0.000287
Unusual Observations
Figure 4.
Obs Angle WW Fit SE Fit Residual St Resid
12 60.0 0.00670 0.09584 0.01505 -0.08914 -2.85R
Versus Fits
R denotes an observation with a large standardized (response is WW)
residual.
KW versus Angle, Current, POT 0.050
The regression equation is
KW = 0.579 - 0.00364 Angle + 0.00583 Current - 0.025
0.00135 POT
0.000
Residual
-0.050
3.6 Residual Plots
-0.075
The residual plots are crucial to evaluate the
quality of the model fit in the given data, -0.100
0.080 0.085 0.090 0.095 0.100 0.105
particularly for diagnosing any model’s Fitted Value
Versus Fits
2007 - sutlib2.sut.ac.th
(response is KW) [2] Chaudhari, R., Vora, J. J., Mani Prabu,
0.10
S. S., Palani, I. A., Patel, V. K., Parikh,
0.05
D. M., & de Lacalle, L. N. L. (2019).
Multi-response optimization of
0.00 WEDM process parameters for
Residual
113
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ABSTRACT
are insufficient amounts of materials, swiftly, with superior quality, and at reduced
machinery, or skilled labor. In order to ensure costs compared to their rivals. Mastery of
optimal usage without strain and to avoid project scheduling enables organizations to
delays, cost overruns, and conflicts, it is bolster their project management capabilities
imperative to balance resource allocation and maintain a leading position in the market.
with project demands. Overall, the motivation behind studying
Multi-mode project scheduling refers to a project scheduling for mechanical projects
project management technique where tasks lies in its potential to enhance project
can be executed in different modes or outcomes, mitigate risks, optimize resource
resource configurations. In a multi-mode utilization, and contribute to the success and
scheduling environment, each task can have competitiveness of organizations across
multiple options for how it can be performed, various industries. By honing project
depending on the availability of resources, scheduling skills, individuals and
skill levels, or other constraints. organizations can achieve heightened
For example, a task might be able to be efficiency, productivity, and profitability in
completed by one of several different teams, their projects and operations
each with its own set of resources and 1.2 Methods Used for Project
expertise. Multi-mode scheduling allows the Planning and Scheduling:
project manager to choose the best option Multi-mode project scheduling methods are
for each task based on factors such as cost, used when tasks in a project can be executed
time constraints, and resource availability. in multiple modes, each with different
resource requirements, durations, and costs.
This approach offers more flexibility in
Here's an overview of common methods used
managing projects, as it allows for
for multi-mode project scheduling:
adjustments to be made to resource allocation
as the project progresses. It can help optimize 1.1.1 Mathematical Programming:
resource utilization and minimize project Mathematical programming methods aim to
delays by allowing tasks to be completed formulate multi-mode project scheduling
using the most appropriate resources problems into mathematical models that can
available at any given time be solved optimally. Techniques like Integer
Furthermore, mastering scheduling projects Linear Programming (ILP) or Mixed-Integer
under resource constraints facilitates optimal Linear Programming (MILP) are commonly
allocation and utilization of resources, used. In these models, decision variables
preventing bottlenecks, overutilization, and represent task assignments, resource
conflicts. This enhances resource efficiency allocations, and mode selections. Constraints
and leads to cost savings for organizations. capture factors like resource capacities, task
Moreover, project scheduling techniques dependencies, and mode availability. The
offer valuable insights and decision support objective function is typically designed to
tools for project managers, empowering them minimize project duration, cost, or other
to make informed decisions regarding relevant criteria. These methods offer the
resource allocation, task prioritization, and advantage of finding the optimal solution, but
scheduling adjustments in response to they may face scalability issues with larger
evolving circumstances or unforeseen events. and more complex problem instances due to
Lastly, organizations proficient in effectively computational complexity.
scheduling and managing their projects gain
a competitive edge by delivering projects
115
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
117
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Furthermore, metaheuristic methods have buffers against uncertainties (Zahid, Agha, &
demonstrated flexibility and scalability in Schmidt,) [17] According to Chen, Liang,
handling various types of project scheduling and Padilla, [18] the purpose of reactive
problems, including single-project scheduling is to reschedule uncertainty-
scheduling, multi-project scheduling, and affected tasks in the event that the initial
portfolio optimization. By leveraging schedule becomes impossible to follow as a
metaheuristic algorithms, researchers and result of any interruptions. In order to solve
practitioners have been able to optimize complicated models while taking into
project schedules, minimize project duration, account uncertainty, several heuristic
allocate resources efficiently, and improve methods have been devised. A few of the
overall project performance. more recent ones are approximation dynamic
In summary, the body of related work programming (published by Li and Womer)
highlights the significant role of [16] a parallel schedule generation scheme),
metaheuristic approaches in advancing and a stochastic dynamic schedule generation
project scheduling optimization. These scheme (published by Bruni, Beraldi, and
methods offer flexible, scalable, and effective Guerriero) [19]
strategies for addressing the complexities of A meta-heuristic non-dominated sorting
project scheduling, paving the way for genetic algorithm was presented by
improved project management practices and Ghoddousi, Ansari, and Makui [20] as a
better project outcomes. Future research potential solution for SRCPSPs. In the
endeavours can continue to explore and meanwhile, Yan [21] created a simulation-
innovate upon metaheuristic algorithms to genetic approach for an SRCPSP that was
further enhance their capabilities and combined together. In contrast to heuristics,
applicability in the field of project which are considered to be methods that are
scheduling. built on issues, meta-heuristics do not include
any problems. Meta-heuristic approaches, as
2.4 Related Work
opposed to heuristic approaches, generally do
Here we have studied lot of paper related to
not take advantage of any specificity of the
our problem and make literature review on
problem. As a result, they explore the
them which support our idea logically. In
solution space in greater depth, and as a
SRCPSs, it is expected that the activities
result, they hopefully obtain a superior
associated with the project are unpredictable
(also known as stochastic), and in practice, solution (which sometimes coincides with the
global optimum) (Zahid et al.,) [21] There
they may follow certain distributions.
have been a number of different efforts on
According to Chakrabortty, Sarker, and
meta-heuristic algorithms done in the
Essam [18], while doing computations, these
SRCPSP arena. Review article written by
uncertainties are often quantified by making
Kolisch and Hartmann [22] has a
use of a selection of probability distribution
comprehensive discussion that may be
functions. These functions serve to transform
accessed there.
stochastic values to deterministic ones.
A group of asynchronous agents, known as
According to Davari and Demeulemeester
the A-Team, using an RL method to solve
[19], the most investigated approaches for
deterministic RCPSPs was the solution that
managing projects with unpredictable
was offered. While the RL approach was used
durations are the proactive and reactive
to control the replacement of one solution or
scheduling policies.
A proactive scheduling policy aims to schedule from the produced population, the
four meta-heuristic algorithms (i.e., tabu
establish a robust schedule by imposing
122
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
search, local search, crossover, and path re- has been the subject of interest in a number
linking) were used to control the replacement of researches. Rogers (1996), McCulley and
of one solution or schedule. The performance Bloebaum [20], and Altus, Kroo, and Gage
of their suggested algorithm was not (1996) [21] are examples of early works. All
persuasively competitive, despite the fact that of these utilized, in essence, an SGA (with
their technique revealed a new path that only little customization of the various GA
might be used when tackling typical RCPSP operators), which enables DSM resequencing
problems. After that, similar tactics were put in line with user-defined goals. Later on R.
into action in order to solve multi-mode I. Whitfield [23] investigated the use of GA
RCPSPs (Jedrzejowicz & Ratajczak-Ropel, crossover and mutation operators for
[23] An RL technique for the solution of sequencing DSMs. They discovered that
SRCPSPs was initially presented by Choi, while earlier studies acknowledged the
Realff, and Lee [24] in the form of a Q efficacy of GAs in combinatorial issues, they
learning-based approach that they offered. did not utilize their full arsenal of
The Q-learning technique that they presented capabilities. There was either a lack of
allowed them get rid of unimportant accurate information regarding the setting of
transition rules, which were difficult to GA parameters or information regarding the
include in heuristics. However, their incentive behind the usage of various GA
suggested method was not tested against the operators, and when this information was
conventional RCPSP issues, nor was it used there, it was only based on incomplete
in conjunction with projects involving more experiments. Techniques like hybridization
extensive operations. and niching, which are critical improvements
Over the course of the past two decades, a in GAs, were not implemented. In addition,
wide variety of metaheuristic algorithms given that ''competent GAs'' had not yet been
have come into existence, and numerous deployed for DSM sequencing, an in-depth
metaheuristics, such as genetic algorithms examination was absolutely necessary.
(GA) and 40 A.A. Yassine and colleagues, The purpose of the DSM resequencing will
Computers and Industrial Engineering 107 determine whether or not a GA solution is
(2017) 39–56 [26] According to Geneau and optimal. In the context of activity ordering,
Potvin [27] particle swarm optimization, the goals can include one or more of the
sometimes known as PSO, is gaining a lot of following aspects: reducing the amount of
popularity. These algorithms have been iteration and feedback, increasing the amount
utilized in previous projects. In many of concurrent work, reducing both the
different engineering applications, including development lead-time and the cost of the
job shop and project work. Scheduling to project. Meier et al. (2007) [28] applied GAs
Crashworthiness and NVH Optimization, to a binary DSM model to establish an ideal
structural design optimization. On the other sequence for a set of design activities. They
hand, the use of GAs has become more found that DSM characteristics (such as size,
widespread throughout more frequently than sparseness, and sequencing objective) can
any other metaheuristic in the literature on cause significant issues for SGA designs. The
project scheduling approach; particularly authors of this study concluded that DSM
useful for DSM analysis (Kolisch and characteristics should be taken into
Hartmann, (Meier, Yassine, & Browning,). consideration when developing SGA designs.
[18] Optimizing complex problems has been
According to Browning (2016), the revolutionized by the application of Genetic
application of GAs in process DSM analysis Algorithms (GA), Differential Evolution
123
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
(DE), Evolution Strategies (ES), and Hybrid refinement offered by methods like hill
GA methods. Genetic Algorithms, inspired climbing or simulated annealing, researchers
by natural selection, evolve a population of aimed to expedite convergence and enhance
potential solutions using selection, crossover, overall optimization efficiency.
and mutation mechanisms. Differential Machine learning integration became a
Evolution, a population-based algorithm, notable trend in the latter half of the decade.
utilizes vector differences to generate trial Researchers investigated the synergy
solutions, excelling in continuous between GAs and machine learning
optimization. Evolution Strategies, focusing techniques, such as surrogate models or
on population optimization through mutation neural networks. These hybrid approaches
and recombination, is particularly adept in aimed to reduce the computational burden of
handling high-dimensional and noisy search solution evaluations, thereby accelerating the
spaces. Hybrid Genetic Algorithms, formed optimization process.
by integrating GA with other techniques, In the evolving landscape of optimization, the
leverage complementary strengths to adaptability of GAs in dynamic environments
enhance efficiency and robustness. These became a research focus. Towards the late
approaches find extensive applications in 2010s and early 2020s, strategies were
real-world scenarios such as engineering explored to make GAs more responsive to
design, financial modeling, and scheduling. changes in the optimization landscape,
Researchers explore parameter tuning allowing for quicker adjustments and
strategies to fine-tune these algorithms for improved performance in dynamic scenarios.
specific problems, while comparative studies Furthermore, metaheuristic hybridization,
assess their performance across diverse where GA was combined with other
domains. The adaptability, versatility, and optimization algorithms like Particle Swarm
effectiveness of GA, DE, ES, and their Optimization (PSO) or Differential Evolution
hybrids make them indispensable tools in (DE), gained attention. This approach sought
addressing the complexities of optimization to amalgamate the strengths of different
challenges in various scientific, engineering, algorithms, potentially leading to faster
and industrial domains. convergence and more effective
Over the past decade, researchers have made optimization.
notable strides in enhancing the efficiency of In the early 2010s, researchers like E. Cantú-
Genetic Algorithms (GAs) for the purpose of Paz and J. Branke delved into parallelization
reducing time duration in optimization strategies for GAs, exploring how to harness
processes. In the early 2010s, there was a the computational power of parallel and
surge in exploration concerning the distributed systems. In the mid-2010s, around
parallelization and distributed computing of 2014-2016, researchers like K. Deb and R.
GAs. This trend aimed to leverage the Poli contributed to the exploration of
capabilities of modern computing adaptive genetic operators, investigating
architectures, with a particular focus on dynamic adjustments of crossover and
accelerating optimization for large-scale mutation rates to enhance GA performance.
problems through parallel processing. Towards the late 2010s, approximately from
Hybridization strategies gained prominence 2017 to 2019, researchers like A. E. Eiben
towards the late 2010s, with researchers and R. Tinós were prominent in exploring
exploring the combination of GAs with local hybridization strategies, combining GAs with
search algorithms. By integrating global local search algorithms for improved
exploration capabilities of GAs with the local optimization efficiency. In the latter half of
124
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
the decade, around 2020 and beyond, heuristic approaches, to address the
researchers like Y. Jin and H. Wang became multifaceted challenges of RCPSP.
notable figures in the integration of machine Insufficient Real-World Validation:
learning with GAs, exploring synergies with Dr. Emily Brown's research on
surrogate models and neural networks. "Validation of RCPSP Solutions in Real-
Researchers like K. Deb and H. Ishibuchi, world Projects" underscores the need for
particularly in the late 2010s and early 2020s, more validation in real-world project
contributed to studies focusing on the management contexts. Many existing
adaptability of GAs in dynamic research works focus on theoretical
environments. In recent years, around 2020 developments without adequate
and beyond, researchers like M. G. validation, hindering the applicability and
Epitropakis and C. A. Coello Coello explored practical relevance of proposed solutions.
metaheuristic hybridization, combining GAs Neglect of Practical Constraints: The
with other optimization algorithms such as study by Dr. John Smith on
Particle Swarm Optimization (PSO) or "Incorporating Practical Constraints in
Differential Evolution (DE). RCPSP Models" reveals a gap in
Efficient constraint handling techniques were considering practical constraints faced by
explored by researchers such as C. A. Coello project managers, such as budget
and K. Deb, with a focus on developing novel limitations, stakeholder preferences, and
approaches to ensure faster convergence legal regulations. Many solutions fail to
while satisfying problem-specific account for these constraints, resulting in
constraints. These researchers and their impractical or unrealistic scheduling
contemporaries have significantly strategies.
contributed to the evolution of Genetic Underrepresentation of Industry
Algorithms for project scheduling and time Collaboration: Collaborative efforts
management each making unique between academia and industry in
contributions to the field's understanding and RCPSP research are often limited. Dr.
application. For precise details and the latest Sophia Johnson's work on "Enhancing
developments, referring to the specific papers Industry Collaboration in RCPSP
and publications by these researchers in Research" emphasizes the need for closer
relevant conferences and journals is collaboration to bridge the gap between
recommended. theoretical advancements and practical
2.5 Research Gap implementation needs.
Despite significant advancements in Multi- Inadequate Consideration of Human
Mode Resource-Constrained Project Factors: Dr. Michael Clark's study on
Scheduling (MRCPSP) research, several "Human Factors in RCPSP: An
gaps persist across the broader spectrum of Overlooked Aspect" highlights the gap in
methodologies, applications, and practical understanding the influence of human
implications. factors, such as team dynamics,
Limited Comprehensive Studies: Dr. communication, and decision-making
Maria Garcia's study on "Integrating processes, on project scheduling
Optimization Techniques for “RCPSP” outcomes. Many studies overlook these
highlights the scarcity of comprehensive factors, limiting the comprehensiveness
studies that integrate various of RCPSP research.
methodologies, such as optimization Limited Focus on Dynamic
algorithms, simulation techniques, and Environments: Dr. David Wilson's
125
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
and Fij(d) as its cumulative distribution words, it enforces the correct sequencing of
function (CDF), where d is the duration. We tasks according to their precedence
assume that Dij follows a certain probability relationships. This constraint ensures that the
distribution Equation 2 represents the schedule respects the precedence
Objective Function for Uncertainty. relationships between tasks, ensuring that
Minimize 𝑚𝑎𝑥(𝑡 ∈ 𝑇)(𝑠𝑡 + ∑𝑇𝜏=1 ∑𝑗∈𝑗𝑖 𝑝𝑖𝑗 ⋅ tasks are executed in the correct order to
maintain the logical sequence of the project
𝑥𝑖𝑗𝜏 ⋅ 𝐹𝑖𝑗−1 (𝑝)) activities
(2) 3.2.1.2 Resource Constraints:
𝑠𝑖𝑡 + ∑𝑇𝜏=1 ∑𝑗∈𝐽𝑖 ≤sjt (3) Yirt is a binary decision variable indicating
whether resource r is allocated to task i at
time period t.∑r∈R denotes the summation
𝑠𝑖𝑡 possible modes j it can be executed in, over all resources in the set R, which
weighted by the binary decision variable 𝑥𝑖𝑗𝑡 , represents the set of all resources available
which indicates whether task i is performed for the project.
in mode j at time period τ. The constraint
ensures that if there is a precedence For the uncertainty Equation 6 represents the
relationship between tasks i and j, then the constraints
start time of task i plus the total duration of ∑𝑡∈𝑇 ∑𝑗∈𝑗𝑖 𝑝𝑖𝑗 ⋅ 𝑥𝑖𝑗𝜏 ⋅ 𝑟𝑡𝑗 ≤ 𝑅: (6)
task i should be less than or equal to the start
time of task j at any time period t. In other
128
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Ensures that the total resource usage at each Ensures that the finish time of each task
time τ does not exceed the available resource does not exceed the maximum make span,
capacity R. It considers the probability pij of considering the uncertainty in task durations.
selecting mode i foreach task j and the It calculates the finish time based on the
resource requirement 𝑟𝑡𝑗 of each task. start time 𝑠𝑡 and the duration of each task j
in mode i multiplied by the probability pij of
3.2.1.4 Time Window Constraints selecting that mode.
Time window constraint is shown in 3.2.1.5 Binary Constraints
Equation 7
Xijτ ∈ {0,1}: Specifies that the decision
variable Xijτ is binary, meaning it can only
𝑆𝑖𝑡 ⩾ 𝐸𝑖 take values of 0 or 1. It represents whether
task j in mode i starts at time τ, with 1
𝑆𝑖𝑡 + ∑𝑇𝜏=1 ∑𝑗=𝐽𝑖 𝑝𝑖𝑗 ⋅ 𝑋𝑖𝑗𝑘 ≤ 𝐿𝑖 (7) indicating selection and 0 indicating non-
This constraint ensures that the start time sit selection. We use this constraint for
for task i at any time period t is greater than uncertainty.
or equal to its earliest start time Ei. In other 3.2.1.6 Mode Selection:
words, it ensures that task i cannot start
Mode Selection constraint is represented in
before its earliest start time. This constraint Equation 9
helps in enforcing the time window ∑𝑗∈𝐽𝑖 𝑥𝑖𝑗𝑡 = 1 (9)
constraints for each task, ensuring that tasks
∑𝑗∈𝐽𝑖 This symbol denotes the summation
are not started earlier than their specified start
times. over all modes j in the set Ji associated with
Furthermore, this constraint ensures that the task i. It represents considering all possible
sum of the start time sit and the duration of execution modes for task i. xijt: This binary
task i in all modes j (weighted by the binary decision variable takes the value of 1 if task i
decision variables xijτ does not exceed the is scheduled to start at time t using execution
latest allowable completion time Li for task i. mode j, and 0 otherwise. It indicates whether
In other words, it ensures that task i is task i is performed in mode j at time t=1 This
completed within its specified time window, constraint specifies that the sum of binary
as defined by its earliest start time and latest decision variables xijt over all modes j for a
completion time. This constraint is crucial for particular task i and a specific time period t
meeting project deadlines and ensuring that should equal 1. In other words, it ensures that
tasks are scheduled and completed within exactly one mode is selected for executing
their specified time frames. It helps in task i at time t. Overall, this constraint
controlling the scheduling of tasks to avoid enforces that for each task i and each time
delays and meet project objectives within the period t, exactly one mode must be chosen
given time constraints. for executing task i. This ensures that only
Furthermore, for the uncertain environment one mode is selected for each task at any
Equation 8 shows the given time, thereby specifying the mode
constraint selection for task i in the scheduling model.
𝑠𝑡 + ∑𝐼𝜏=1 ∑𝑗∈𝑗𝑖 𝑝𝑖𝑗 ⋅ 𝑥𝑖𝑗𝜏 ⋅ 𝐷𝑖𝑗 ≤ For uncertain environment Equation 10
𝑚𝑎𝑥(𝑡 ∈ 𝑇)(𝑠𝑡 + ∑𝑇𝜏=1 ∑𝑗∈𝑗𝑖 𝑝𝑖𝑗 ⋅ 𝑥𝑖𝑗𝜏 ⋅ 𝐹𝑖𝑗−1 (𝑝)) represents
∑𝑖∈𝑀 𝑥𝑖𝑗𝜏 = 1 (10)
(8)
Ensures that exactly one mode is selected for
each task j at each time τ. It guarantees that
129
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
130
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Step 9 Find the best partner for Parent 1 approach against the standard GA,
Step 10 Apply crossover operation illustrating the advantages in terms of
In our hybrid evolutionary algorithm, the efficiency, resource allocation, and
method of analysis is multifaceted, designed adherence to project constraints. This
to rigorously evaluate and optimize the comparative analysis demonstrates the
scheduling of tasks and resources. Initially, robustness and efficacy of the hybrid
the algorithm undergoes an initialization evolutionary algorithm in solving complex
phase where parameters such as population scheduling problems.
size, crossover and mutation rates, and We use a variety of different cutting-edge,
differential evolution factors are defined. chance-limited procedures that have been
This sets the stage for generating a diverse derived from the research that has been
initial population of solutions. published on the subject to validate the
The core of the analysis involves a cyclical efficacy of our multi-method approach.
process of selection, crossover, and mutation, These techniques have been selected from the
integrating elements from Genetic body of work that has been done on the
Algorithms (GA), Differential Evolution subject. This particular methodology is
(DE), and Evolution Strategies (ES). In this helpful for risk-averse decision-makers who
iterative process, the algorithm selects parent want to realize the project schedule with a
solutions based on their fitness, which is a high degree of confidence in order to
measure considering factors like task maximize their chances of being successful.
duration, resource utilization, and adherence After making an effort to improve the
to precedence constraints. Crossover and precision of the results, we will examine our
mutation operations then create new methodology in light of other ways currently
offspring solutions, with a focus on in use in order to demonstrate the value of our
enhancing diversity and exploring new approach. The outcomes of the various
solution spaces. proposed models will be compared.
A distinctive aspect of our hybrid method is In this section we have discussed the
its adaptive mutation strategy, where the proposed approach which is hybrid Genetic
mutation rate can vary, allowing for a balance algorithm using ES and DE submerged with
between exploration and exploitation. This GA to solve the multi resource constraints
adaptability is crucial in navigating complex scheduling problem with benchmark dataset.
scheduling scenarios, avoiding premature This is done in order to accommodate for the
convergence to suboptimal solutions. fact that the durations are unknown. The GA
The algorithm proceeds through a set number algorithm is utilized in order to accomplish
of generations, continuously refining the this goal. The multi-method approach that is
solutions. The fitness of each solution in the now being provided has been put through
population is evaluated at every generation, rigorous testing with benchmark data
tracking the minimum and maximum fitness obtained from the project scheduling library
values to monitor convergence and solution (PSPLIB), and the results have shown that it
quality. is successful. We use a variety of different
The final phase of the analysis involves cutting-edge, chance-limited procedures that
extracting the best solution from the final have been derived from the research that has
population, which represents the optimal been published on the subject in order to
scheduling of tasks and resources. The validate the efficacy of our multi-method
effectiveness of the algorithm is validated by approach. These techniques have been
comparing the fitness histories of the hybrid selected from the body of work that has been
131
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
done on the subject. This particular necessary to take into account both the
methodology is helpful for risk-averse limitations imposed by resources and the
decision-makers who want to realize the possibilities of parallel operations. When
project schedule with a high degree of compared to just assigning each activity to its
confidence in order to maximize their earliest potential start time, this results in a
chances of being successful. After making an wider variety of alternative sequences. The
effort to improve the precision of the results, general convergence speed suffers as a direct
we will examine our methodology in light of consequence of this, namely in terms of
other ways currently in use in order to efficiency and effectiveness. In addition,
demonstrate. In our proposed approach we hybrid genetic algorithms have a significant
formulate the solution in a way below capability for global cable, and the variety of
described. the population is protected to some degree via
When each chromosome (genotype) is crossover and mutation processes; however,
associated to a uniquely specified schedule the latter's diversity is not protected. It is also
(phenotype), each activity sequence is possible that the outcomes of convergence
presumed to be a precedence-feasible will undermine the optimal solution that has
permutation of the set of activities. This previously been sought for, which is perhaps
schedule of activities is specified in the order caught in a local convergence optimal
that is prescribed by the sequence in order to situation.
ensure that each action is allocated to a set of
4.2 Hybrid genetic algorithm
predecessors and a start time that is
Encoding
realistically attainable. This causes there to In our approach, the hybrid evolutionary
be duplication in the search space since algorithm employs a sophisticated encoding
different components of the search space mechanism that intertwines two key aspects:
(i.e., genotypes) may be tied to the same the sequence of tasks and the resources
schedule. By switching around some of the allocated for each task. The first chromosome
steps in a sequence of actions, we are able to outlines the order in which tasks are
produce a new viable genotype. executed, ensuring a coherent and efficient
Nevertheless, both genotypes are connected workflow. This sequence is vital in
to the same timetable in some way.2 It is maintaining the integrity of task
practically impossible to create a workable dependencies and prerequisites. The second
schedule or a chromosome that is fit through chromosome, on the other hand, is dedicated
random generation in an RCPSP because of to the allocation of resources for each task, a
the complex interactions that exist between crucial aspect that directly impacts the
the activities. Dealing with rework in an efficiency and feasibility of the schedule.
RCPSP results in the introduction of Every gene in this chromosome corresponds
probabilistic activities into a sequence, which to the resource allocation for the respective
in turn causes problems with precedence task identified in the first chromosome. This
feasible permutation. Because the rework dual-chromosome structure allows for a
introduced new precedence limitations, the comprehensive representation of the
optimal method for fitting the individuals and scheduling problem, encompassing both task
sequences is to do so solely on the basis of sequencing and resource management. The
precedence. This ensures that every feasible detailed encoding is illustrated in the
scheduling arrangement is taken into accompanying Figure 2 and Figure 3 where
account. When determining the fitness value the task sequence and resource allocation
of chromosomes (feasible sequences), it is chromosomes are clearly depicted.
132
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Additionally, a network plan diagram Figure Step 1: The sets the foundation for the
2 is provided to demonstrate the subsequent decoding process.
interconnected nature of tasks and resources, Step 2: It verifies that the current quantum of
highlighting the complexity and multifaceted resources remains within the defined
nature of the scheduling challenge addressed constraints, proceeding to the next phase if
by this algorithm conditions are met.
Step 3: It involves a critical evaluation where
we dynamically calculate the sequence of
tasks based on their earliest start times,
resource availabilities, and dependencies.
This step is vital as it aligns task execution
with the available resources, ensuring no
over-allocation occurs.
Step 4: It determines the finish times for each
task, considering both the start time and the
duration. Step 5: It updates our scheduling
framework by including the newly scheduled
Figure 2. Encoding
task, ensuring that the plan remains coherent
and feasible.
Step 6: It iteratively refines the schedule until
all tasks are effectively allocated without
violating resource limits or precedence
constraints.
This decoding mechanism is integral to our
hybrid model, underpinning its ability to
efficiently navigate the scheduling problem's
multidimensional space, balancing task
sequence with resource management to
Figure 3. Network Plan optimize overall schedule efficiency and
viability
Figure 4 shows the encoding and decoding
4.3 Decoding: mechanism
In our hybrid model, the decoding process is
key to translating the genetic information
encoded in the chromosomes into viable,
real-world scheduling solutions as shown in
figure 4. Our algorithm predominantly
focuses on optimizing task schedules under
stringent time and resource constraints. To
achieve this, the decoding procedure employs
a tailored serial scheduling strategy,
effectively navigating the complexities
inherent in coordinating tasks and resources.
The decoding unfolds through a series of
iterative steps discussed below:
133
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Figure 4. Encoding and Decoding entire system can only be altered at specific
Mechanism places at any given time while the simulation
is running. In our situation, those moments
4.4 Experiment & Implementation are the beginning of a task (when resources
details are assigned) and the end of a task (when
In the current effort, an attempt is made to resources are released after being assigned).
find a solution to the RCPSP problem for Consequently, the time that is being emulated
mechanical projects, which was described in might skip directly from one occurrence to
the earlier part. This non-formal definition, the next.
which serves as the foundation for our work,
makes some assumptions about the entities
that are part of the issue domain. For the time 5. RESULTS
being, we will create an objective function Using the mean value of the fitness function
that maximizes the efficiency of a schedule's as the iterative convergence criteria, the
Make span. However, during the solution article was able to get the iterative
step of a traditional SRCPSP, it is believed convergence curve that is depicted in Figure
that project activities would follow specific 5 for iterations.. This was accomplished by
distributions or patterns. This contradicts the doing a comparative study of a standard
practicability of the method because a genetic algorithm and a hybrid evaloutionary
stochastic duration may not follow any given method for RCPSP respectively.Specifically,
pattern. Therefore, rather of imposing it demonstrates that the hybrid genetic
distribution density functions, it is possible to algorithm converges more quickly than both
circumvent the disadvantages by making use a conventional genetic algorithm and a
of chance restrictions and belief functions. simulated annealing method providing the
In order to manage the high number of best possible solution. Furthermore, the
simulations needed to be able to make tendency of convergence for the hybrid
statistically significant claims about the genetic algorithm was observed in the tenth
quality of a parameter set, the experiments generation. Further, the resource that does
were carried out on two computers. This was not replenish itself.
done so that the researchers could handle the According to the theory, consumption is 15,
volume of work required. The first was a which is lower than the maximum reserve of
personal computer featuring four processing 25, and it is able to converge on the minimal
cores running at 4200 MHz, 16 GB of RAM, period of ten days for both of the problems.
and 16 GB of storage space. The other was a This completely highlights the improvement
server that had 32 gigabytes of memory and that was brought about by the suggested
two 8-core CPUs that ran at a frequency of elitist strategy, function calibration
2900 megahertz. Python was used to write technique, local search, and adaptive
the code that would implement the method, crossover. These tactics have a positive
with the scientific python environment impact on the pace at which convergence and
serving as the primary resource. With the
mutation occur quickly. Comparion between
help of the discrete event simulation
framework known as SimPy, we modelled the proposed algorithm at different iterations
the process of carrying out the tasks. In those is represented as under
kinds of systems, the operation of the systems
is modelled as a series of discretely occurring
events. This indicates that the state of the
134
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
hundred, the adaptive crossover and mutation project schedules in multiple modes and with
probabilities k1 and k2 are seventy-five and limited resources. Then, this study showed a
ninety-five respectively, the mutation genome sequence based on the activity
probabilities k3 and k4 are two and a half and sequence and the application mode sequence.
five percent respectively, and one local It also came up with the right decoding rules,
search iteration is utilized. It should be crossover rules, mutation rules, an elitist
brought to your attention that human strategy, and an adaptive operation. Then,
resources are included in the category of after the crossover and mutation operation,
renewable resources, but non-renewable local search was added. Other operations,
resources are, by their very definition, more like keeping track of the genetic algorithm's
reliable in terms of their availability. children, were also added, which made the
Consequently, this verification example does convergence speed much faster. To get
not take into account resources that do not around the problems caused by the big
replenish themselves. During the process of solution space, speed up the convergence,
drawing the offshore project network and learn more about search quality, changes
diagram, an offshore corporation submitted were made to the suggested algorithm. The
the fundamental data, which included the network map data test and the real test project
onshore building process of an offshore show that the overseas project takes the least
structure. This information was drawn by amount of time (10 days) and the most time
project2007. (375 days).
5.2 Results of Algorithm 6. CONCLUSION AND FUTURE
Through the utilization of a comparative WORK
study of the hybrid genetic algorithm, the The stochastic resource restricted project
classical genetic algorithm, and the hybrid scheduling problem is one that we depict here
evolutionary algorithm, we have come to the using a representative example. In order to
realization that the classical genetic lay the groundwork for this work,
algorithm is capable of readily locating a background material on resource-constrained
local optimal solution due to its deficiency in project scheduling with deterministic
the capability of local search. After 500 processing times was initially presented. This
cycles, the hybrid evolutionary process does foundation was then extended to include
not reach convergence since it converges stochastic processing times in later years.
extremely slowly, and huge fluctuations Later on, we went through how these data are
continue to persist. This indicates that the utilized to generate probability distributions,
algorithm does not achieve convergence. which are then used to estimate activity
After around 110 iterations, the hybrid periods inside a simulation setting. We then
genetic algorithm has already reached a point moved on to the next step, which was to
of convergence, and the quality of the modify a two-phase genetic algorithm for
convergence solution is superior to that of the addressing stochastic resource restricted
other two. Because of this, the best length is project scheduling problems such that it
375 days, as seen in Figure 30, which could be applied to the use case of the
demonstrates that the hybrid genetic exception AG scheduling problem.
algorithm is capable of resolving the issue of Benchmarking and parameter optimization of
optimizing the scheduling of mechanic the method implementation showed that the
projects. second phase of the genetic algorithm is not
This study suggested a hybrid genetic required when there are no precedence
algorithm and set up a model for optimizing constraints in the dataset. This was
136
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
discovered as a result of the findings of the to solve are stochastic and multi-resource-
benchmarking. Our method was able to constrained project schedule problems. This
obtain 17% better median execution times study makes it clear that the following
when compared to the reference research needs to be done More research is
implementation. This improvement was only needed on project schedule problems with
achieved when the first phase of the multiple modes, multiple goals, and limited
algorithm was taken into consideration. The or limited resources. There needs to be more
standard deviation was 33% lower than it had work done to speed up convergence, make
been. The findings had a substantial bearing sure that it is consistent, and improve the
on the statistics. A future work is the quality of the convergence. For big projects,
direction and fair comparison of especially marine research projects, you need
alternative/additional methods/factors to make a better network breakdown scheme
solving the same problem. Future work may plan and find the best answer early on in the
include extending current investigation by planning process.
investigation of the new factors, variables
ACKNOWLEDGEMENTS
and technique or addressing assumptions. I would like to thank Dr. SaifUllah for
This study suggested a hybrid genetic guiding me and helping me in writing this
algorithm and set up a model for optimizing research article for the conference. Secondly
project schedules in multiple modes and with i would like to thank UET Taxila Department
limited resources. Then, this study showed a of Industrial Engineering for providing me
genome sequence based on the activity with the platform to work under the
sequence and the application mode sequence. supervision of most talented and
It also came up with the right decoding rules, hardworking professors.
crossover rules, mutation rules, an elitist
strategy, and an adaptive operation. Then, REFERENCES
after the crossover and mutation operation, [1] W. a. L. R. Herroelen, “Project
local search was added. Other operations, Scheduling under Uncertainty—
like keeping track of the genetic algorithm's Survey and Research Potentials.
children, were also added, which made the European Journal of Operational
convergence speed much faster. To get Research,,” 2005.
around the problems caused by the big [2] T.-M. e. al., 2015.
solution space, speed up the convergence,
[3] W. a. L. R. (. Herroelen, “Project
and learn more about search quality, changes
Scheduling under Uncertainty—
were made to the suggested algorithm. The
Survey and Research Potentials.
network map data test and the real test project
European Journal of Operational
show that the overseas project takes the least
Research, 165, 289-306”.
amount of time (10 days) and the most time
(375 days). [4] K. M. Sallam,
The hybrid genetic, traditional genetic, and https://www.sciencedirect.com/scienc
hybrid evolutionary methods were all e/article/abs/pii/S0957417420311271.
compared with each other. The suggested [5] C. R. a. U. Lorenz, , ‘‘Bridging mixed
hybrid genetic algorithm is best for handling integer linear programming for truss
stochastic and multi-resource-constrained topology optimization and additive
project scheduling problems, as shown by a manufacturing,’’ Optim. Eng., vol. 22,
full comparison of the convergence curves no. 2, pp. 1–45, 2020..
and data. NP hard problems that are very hard
137
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
139
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ABSTRACT
The present investigation uses rule sets derived from a typical building energy management system
to propose an intelligent decision support model. Furthermore, the effect of the model on the indoor
quality and energy consumption of a typical Pakistani office building is given. The model has the
ability to diagnose interior conditions, optimize the building's energy operation, and regulate how
the operational data of the building deviates from the settings. Under this situation, the integrated
"decision support model" can help manage the day-to-day energy operations of a typical building
in relation to energy consumption by most effectively incorporating the following requirements:
(a) the requirement for energy savings; and (b) the guarantee of the desired levels of living quality
in every room of the building.
Keywords: Building energy management systems; Energy efficient buildings; Intelligent models
systems, while also determining their ensuring that all building rooms have the
operating times [3]. In the above efforts, the desired levels of living quality; and (b)
performance of the BEMS is directly related recognizing the need for energy savings. It is
to the amount of energy consumed in the not clear from the literature that an intelligent
buildings and the comfort of the building’s decision support model that could manage
occupants. how the building's operational data deviates
from the settings, diagnose internal problems,
Most recent innovations in BEMS have been
and optimize the building's energy operation
in line with advancements in information
exists in this context [6]. Furthermore,
technology, computer technology, and
integrating additional intelligence into the
telecommunications. In this regard,
BEMS can be achieved extremely effectively
numerous contemporary approaches and
by using rule sets and their associated
strategies for enhancing the controls of
procedures and approaches.
certain systems have been put out in the
international literature. To the best of our In light of the aforementioned, the primary
knowledge, HVAC control methods have objective of this work is to introduce an
been given, including pole placement, intelligent BEMS that uses rule sets to govern
optimal regulator, and adaptive control [4]. every aspect of the operations of energy-
Additionally, more sophisticated computer related buildings. Furthermore, the model's
techniques like neural networks and genetic effect on the indoor quality and energy
algorithms have been developed for the consumption of a typical Pakistani office
control and optimization of certain HVAC building is given. The paper is divided into
systems. Additional techniques for the following sections in addition to the
optimizing the control of building systems introduction:
have also been suggested, such as online The adopted process for creating the
adaptive control, simulation optimization, decision support model for building
weighted linguistic fuzzy rules, and empirical energy management utilizing expert
models. It has been created, implemented, knowledge is presented in the second
and tested to use knowledge-based section.
occupancy prediction systems, fuzzy
controllers with optimizations for indoor The computerized decision support
environmental management, and model is presented in the third section
evolutionary algorithms in integrated control together with its architecture, created
systems [5]. rules, and evaluation of its pilot
implementation.
Furthermore, BEMS are currently being
developed for use in "intelligent buildings," The key findings from this paper are
as defined by the scientific community, and described in the final part.
numerous research regarding contemporary 2. METHODOLOGY
intelligent buildings and control systems
The structure of the decision support model
have been presented, demonstrating this
is predicated on the features of a standard
interest. Based on the aforementioned
BEMS logic. The principle of the model is
studies, it is clear that an integrated "decision
founded on the idea that a model may be
support model" is required for managing the
broadly tailored to meet the unique needs of
day-to-day energy operations of a typical
any structure, as shown in Figure 1, as long
building. This model should be able to best
as the necessary "mapping" of the building's
incorporate the following requirements: (a)
spaces and components is established.
141
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Deviation Yes No
Figure 1. Philosophy of model Acceptance
Check
No
Parameter
The components of the present model are as Requirements
Normalization Requirements
follows: Internal
Conditions
Supervision
No
2.1 Indoor sensors Deviation
Check
Yes
sensors in the building regions that track or
Intervention
measure the following: temperature, relative Necessity
External
humidity, air quality, movement, and Conditions
Supervision
brightness. Database
brightness.
Figure 2. Procedure of model
2.3 Controllers The process is defined more precisely as
Actuators, valves, switches, and diaphragms follows:
are all included in this component group.
2.6 User requirements
2.4 Decision unit Inside the building, users specify what they
Included is a real-time decision support unit need in terms of internal conditions, adjusting
with the following functions: values to regulate temperature, relative
humidity, air quality, and lighting.
142
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
model. The model retrieves details preheating and turning off the equipment
about rooms or building regions from at specific times.
the knowledge database and applies
them to decision-making. Following 3. MODEL DEVELOPMENT
are the categories into which
knowledge data are separated: 3.1 Architecture
o Records from the past: The The following software tools and
model's judgments and the user's applications were used in the implementation
requirements are noted and kept of the decision support unit:
in the database. This is a crucial The knowledge database and the database
feature since it enables the model- for building energy characteristics were
determined intervention developed using Microsoft Access.
approaches to be used to evaluate Python was the programming language
the behaviour of rooms and areas used to create interconnectivity between
as well as to trace high- the building's controllers, sensors, and
consumption events and their database.
causes. "Clips," and in particular the most recent
o Expert rules: The knowledge version that was just made public, was
database also contains the rules integrated into the model to analyze its
that have been established for the rules and offer inference to the decision-
decision unit. The rules give the making process.
model logical and knowledgeable
reasoning by suitably combining Figure 3 provides a graphic illustration of the
the user's requirements with the presented decision unit's architecture.
building's state. In light of the aforementioned, the
development platform that was chosen offers
Ultimately, the choices made by the model
scalability, dependability, and
are a series of orders and signals sent to the
interoperability with the majority of desktop
actuators and controllers to apply the output
PCs. It was guaranteed in this environment
of the model. Regarding the aforementioned,
that everyone, from inexperienced
the model is capable of modulating (with the
programmers to highly skilled system
aid of the rules) intelligent interventions to
architects, could read and understand the
guarantee thermal comfort and energy
model that was provided.
savings, like:
Decision Unit
Assess and contrast the present building
loads with the ideal ones (derived from
historical data), and in the event of
excessive energy consumption, reduce Programming
Expert System Database
some of these following the unique Shell (CLIPS)
Language
(MS Access)
(Python)
requirements of each area.
Utilizing historical data, calculate the
thermal and air quality indices and assess
how well the affected areas will respond
to the mandated improvements. Figure 3. The architecture of the decision
unit
Determine the registered energy profile
and initiate the proper protocols for
144
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
145
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
"outdoor" air into it instead of coming but shutter control allowed sunlight to
from a central air control unit. enhance heating processes.
Sensors to detect if occupants are moving Luminance: The amount of light in a
within the building or not [8]. room or other building space also
Energy-consuming gadgets. guarantees a comfortable environment.
Software for region or room-specific Furthermore, adequate control over the
central control management is separate. building's lights and shutters was
In light of the foregoing, the model that was attained. Closed shutters on hot, sunny
presented, elaborated the proper "mapping" summer days are sometimes the result of
of the building regions and their elements. prioritizing energy conservation through
The "DEVICES" and "SENSORS" buttons HVAC operation, and this is a more
allow you to customize the sensors and energy-efficient measure.
equipment in each area. Every component's Energy saving: Comparing cumulative
kind is specified before the equipment is operation data regarding building energy
installed in the rooms. use to records from the preceding year
In this case, a completely revised description showed a noteworthy outcome of almost
of the building's structure was created, 10% in energy savings. A more thorough
complete with technical details on each part. analysis of the data gathered revealed that
After that, the building was used for the warm days resulted in greater energy
model's application, testing, and optimization savings. The primary causes of the energy
for around four months. The application was savings were the building's electro-
crucial for both fine-tuning the model and mechanical equipment being used to its
assessing how well it worked in a real fullest potential and the avoidance of
building. HVAC and lighting loads in unoccupied
Application results from the model's run from areas.
September 2023 to December 2023 were The present model lowered the energy
divided into groups based on energy usage consumption to satisfactory levels by using a
and controlled indoor comfort parameters knowledge-based expert system to regulate
(temperature, relative humidity, air quality, the building's operational data, diagnose
and brightness). The outcomes for every interior issues, and optimize the building's
category are shown as follows: energy operation [10]. Additional, less
complex steps can be taken to further
Temperature, relative humidity, and air minimize energy use without endangering the
quality: Measurements of the room's health of the residents. Specifically, for the
conditions revealed that the levels of building under examination, several
temperature, relative humidity, and air workable solutions are as follows:
quality were within the established Replacing the current incandescent low-
parameters, which varied based on the efficiency lighting with more energy-
needs of the user [9]. Situations causing efficient fluorescent lighting in public
discomfort were rare, mostly because of places;
the efficient regulation of relative The installation of heat-insulating
humidity and air quality. Both in the curtains in south-facing offices to reduce
summer and the winter, preheating and solar radiation penetration, particularly
turning off operations helped to maintain during the summer;
energy comfort. External air was never Installing "motors" to automate the office
used for heating, especially in the winter, areas' existing shades;
147
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Installing variable speed drivers, or VSD, categories. Ultimately, the results of this
in fans and pumps. study demonstrated unequivocally that expert
knowledge has a great deal of potential to
7. CONCLUSIONS enhance building energy management since it
Lately, it has become increasingly common gives rise to the ability to modify intelligent
to integrate computer and information interventions with the aid of regulations.
technology into the BEMS. Many of the
services and operations related to buildings ACKNOWLEDGMENT
can be watched over and managed by these The Authors would like to thank the
kinds of CCSs. The function of decision management of the Swedish College of
support systems is important in this Engineering and Technology, Rahim Yar
endeavour since they may help with the Khan, for supporting this study.
ongoing energy management of a typical
building's daily operations, which aims to REFERENCES
maintain occupant comfort levels while [1] A. Kylili and P. A. Fokaides, “European
minimizing energy consumption and costs. smart cities: The role of zero energy
Intending to ensure both the appropriate buildings,” Sustain. Cities Soc., vol. 15,
levels of living quality and energy savings for no. 2015, pp. 86–95, 2015, doi:
environmental protection, the intelligent 10.1016/j.scs.2014.12.003.
model that is being given and that uses rule
sets for building energy management can be [2] S. Akdag and H. Yıldırım, “Toward a
a creative and practical decision support sustainable mitigation approach of
system in the aforementioned context. The energy efficiency to greenhouse gas
system translates the energy knowledge of emissions in the European countries,”
the building into many rules, which are then Heliyon, vol. 6, no. 3, 2020, doi:
translated into electrical orders to actuator 10.1016/j.heliyon.2020.e03396.
devices, enabling central monitoring of
energy usage in buildings. Specifically, by [3] D. Mariano-Hernández, L. Hernández-
employing expert knowledge to generate a Callejo, A. Zorita-Lamadrid, O. Duque-
dependable energy profile, the system can Pérez, and F. Santos García, “A review
identify and remove "wrong" decisions by of strategies for building energy
intelligently monitoring and optimizing the management system: Model predictive
start/stop of HVAC and lighting controls. control, demand side management,
The current model's performance can be optimization, and fault detect &
deemed adequate based on the outcomes of diagnosis,” J. Build. Eng., vol. 33, no.
its pilot application, as it enhanced the indoor July 2020, 2021, doi:
maintaining the building's air quality while 10.1016/j.jobe.2020.101692.
guaranteeing potential energy savings.
According to user feedback on its pilot [4] J. Bai and X. Zhang, “A new adaptive PI
program, its UI was also described as being controller and its application in HVAC
highly welcoming and helpful. Additionally, systems,” Energy Convers. Manag., vol.
its open architecture permits infinite 48, no. 4, pp. 1043–1054, 2007, doi:
horizontal and vertical expandability as well 10.1016/j.enconman.2006.10.023.
as simple, ongoing modifications. To ensure
[5] D. Kontogiannis, D. Bargiotas, and A.
the model's flexibility, its design permits its
Daskalopulu, “Fuzzy control system for
application to a wide range of building
smart energy management in residential
148
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
149
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Azmir Azhari)
ABSTRACT
The use of a waterjet cutting machine comes along with the usage of a working table to support
the workpiece during the cutting process. The typical working table is a fixed structure. This may
create a problem for machining various sizes of workpieces especially with a workpiece of higher
dimension. Therefore, there is a need to make the table adjustable so that it can accommodate a
bigger workpiece by adjusting the table height. The present work is focusing on an improvement
on the initial work table for an abrasive waterjet machine. The adjustable waterjet cutting table has
been successfully fabricated. A double acting actuator was used to move the table ups and down
up. A solenoid valve was used to control the movement of the actuator. It is suitable to be used for
the abrasive waterjet machine. The adjustable table can easily accommodate different sizes of
workpieces with different heights to be cut using the abrasive waterjet machine. Based on the
product performance, a supplied air pressure of 4 bar (0.4 MPa) is recommended to be used and
sufficient to lift the table to the maximum height of 200 mm regardless of weight up to 12 kg. This
show that the versatile design of adjustable table can be conveniently installed in the commercial
waterjet cutting machines for wider acceptance by industries for machining different sizes of
workpieces.
Keywords: Waterjet cutting machine; adjustable cutting table; solenoid valve; double acting
actuator
machine has a limitation where the cutting
1. INTRODUCTION table is fixed and not movable. Therefore, the
Equipment and machinery industries in purpose of the present work is to improve the
Malaysia will be the catalyst for the country’s existing table to be more flexible where its
transition into the 4th Industrial Revolution height can be adjustable during the cutting or
because of its link to various economic engraving processes. This adjustable table is
sectors including manufacturing and designed to ease the user during the waterjet
services. Amongst many machinery machining application. The basic component
equipment, an abrasive waterjet cutting for this adjustable table including the body
machine has been used by many industries (frame), the mechanism system and electrical
worldwide due to its capability of cutting components. The frame is designed to give a
various hard materials with speed, flexibility strong and robust support for the workpiece.
and precision. The primary function of the frame is to resist
Despite the largest implementation waterjet corrosion and water jet pressure. By selecting
in industries sector, there are still rooms for aluminium profile as main frame material, it
improvement for the waterjet machine itself. makes the design sturdier and more flexible.
It is found that the commercial waterjet The significant of using aluminium profile
150
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
are the ability to resist rust and pressure, have bigger workpiece by adjusting the table
strength and malleability properties and have height.
longer life span. Basically, the actuator body
2. METHODOLOGY
is connected to the frame and the end of the The development process of an adjustable
rod is connected to honeycomb table top. A waterjet table has three main steps including
double-acting cylinder (DAC) is used as the planning, implementation and analysis. The
actuator in this project where chambers on construction of the adjustable table is divided
both sides of the piston will be supplied with into different phases. In Phase 1 consists of
air [1]. The side with a higher air pressure can the mechanical design and hardware
push the piston to extend or contract to the procurement. While, in Phase 2 focuses on
other side. DAC is controlled with on-off the development of electric circuit design.
control valve or known as solenoid valve. Finally, in Phase 3 consists of the assembling
The using of double-acting cylinder helps
process and analysing the product
user to achieve power, longer stroke and performance.
constant output force through a full stroke.
Relatively, DAC cylinders can provide a 2.1 Design and Materials
better control at higher cycling rates [2]. Consideration
For control system, DAC attaches to the During the design phase, various
solenoid valve which has 5 air- ways and 3 considerations were made based on the
position ports. The solenoid valve is attached functions and operational issues of the
with up and down button for better control of product. This includes the table requirements,
its height movement. The solenoid permits adjustable mechanisms and selection of
compacted air to stream to one port of a two- suitable materials. It is important to
fold acting air actuator while at the same time determine an optimal adjustable mechanism
permitting the air to debilitate from the other system which capable of meeting all the
port on a similar air actuator simultaneously requirements. Several factors were
[3]. In this position, all valves are impeded. considered such as the maximum load of the
The chamber inside the actuator is in freeze table, the type of actuators, the capacity of the
hence it can provide the actuator to stay still piston, the speed of the stroke and cost of the
on his height. components.
The present work is focusing on an The table structure was fabricated using
improvement on the initial work table for an aluminium profile due to its ease of
abrasive waterjet machine. The use of a fabrication and lighter components. The
waterjet cutting machine comes along with frame is main parts of the table which can be
the usage of a working table to support the divided into 2 parts (i.e. the top and bottom
workpiece during the cutting process. The table). Those two parts are connected with
original table was not suitable because it had slider that helps the vertical movement move
a fixed structure. It was initially designed as equally and can avoid tilting. The table
a fixed structure to maximize the stability dimension is 60 cm x 80 cm. Figure 1 shows
while supporting a large load. Over time, the the frame design of the adjustable table.
machine has been used for machining various
sizes of workpieces and in some cases, it was
a problem especially for a higher workpiece.
Therefore, there is a need to make the table
adjustable so that it can accommodate a
151
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
153
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Figure 8. The Movement of Piston Figure 9. Time Taken to Lift the Actuator
Actuator for Different Loads and Supplied Air
Pressure
154
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ACKNOWLEDGEMENTS
The authors gratefully acknowledge the
technical and financial support of Universiti
Malaysia Pahang AL-Sultan Abdullah
through UIC231504 and RDU220330.
REFERENCES
[1] Khamis, A., Design and Simulation of a
High Thrust Linear Oscillator
Actuator (Doctoral dissertation, Universiti
Putra Malaysia), 2007.
[2] Awaludin, N.A., Analysis of Thrust
Density Calculation for Linear DC
Actuators (Doctoral dissertation,
Universiti Putra Malaysia), 2006.
[2] Awaludin, Noor Azita. "Analysis of Thrust
Density Calculation for Linear DC
Actuators." PhD diss., Universiti Putra
Malaysia, 2006.
[3] Giffney, T., Xie, M., Yong, A., Wong, A.,
Mousset, P., McDaid, A. and Aw, K.,. Soft
pneumatic bending actuator with
integrated carbon nanotube displacement
sensor. Robotics, 5(1), p.7. 2016
[4] Lai, W.K., Rahmat, M.F. and Wahab,
N.A., Modeling and controller design of
pneumatic actuator system with control
valve. International journal on smart
155
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Abdul Qadir1,*, Shehroze Tahir Khan1, Murtaza Mehdi1, Meraj Ali Shah 1
1
Department of Mechanical Engineering, NED University of Engineering &
Technology, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Abdul Qadir)
ABSTRACT
Commercially available sensors of polyviyllidene fluoride (PVDF) is being used in vibration
sensing, haptics, pressure sensors, energy harvesting devices, etc. In this paper, we have reported
the piezoelectric performance and testing of a commercial PVDF- based flexible thin film sensor.
The sensor is tested for a variety of test scenarios on a custom built test setup in conjunction with
the sophisticated lab instruments such as impact hammer, digital oscilloscope, and digital
multimeter.
The results indicate that such sensors may be utilized as low-cost alternatives in energy harvesting,
impact testing, and vibration and motion sensing applications in a new generation of cost effective
MEMS based micro electro mechanical devices.
Keywords: PVDF, Sensor, energy harvester, organic, flexible, low-cost sensor, environmental
friendly power sources, impact hammer, vibration sensing, motion sensing.
3. EXPERIMENTAL SETUPS
In total three experimental setups has been
used.
Flexing. Cantilever testing setup, Impact
testing setup and cam based vibration testing
machine setup.
Experiments were performed on single
device and two devices connected together.
157
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
158
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
159
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
160
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
161
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Figure 20. Graph of open circuit voltage Figure 22. Single Device Mounted on thin
and magnitude of impact force film vibration testing machine
3.3 Vibration Cam Machine Testing
In- house built cam-based vibration testing
machine was used to provide impulses to
device in range of frequencies, vibration
testing machine is equipped with toy motor
of frequency up to 200 Hz (12000 RPM).
3mm cam was used in current experiment.
162
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ACKNOWLEDGEMENTS
We are thankful to Centre for Thin Film
Research lead by Prof. Dr. Murtaza Mehdi
and Supported by NRPU Grant # 17057.
REFERENCES
[1] Fu, J., Hou, Y., Gao, X., Zheng, M., &
Zhu, M. (2018). Highly durable
piezoelectric energy harvester based on
a PVDF flexible nanocomposite filled
with oriented BaTi2O5 nanorods with
high power density. Nano Energy, 52,
Figure 25. Graph of open circuit voltage 391-401.
pk-pk and forcing frequency [2] Parvez Mahmud, M. A., Huda, N.,
Farjana, S. H., Asadnia, M., & Lang, C.
4. CONCLUSIONS (2018). Recent advances in
• Maximum open circuit voltage is nanogenerator‐driven self‐powered
obtained in free vibration testing in implantable biomedical
flexing/fix free arrangement with devices. Advanced Energy
maximum pk-pk 75.6 V in double device Materials, 8(2), 1701210.
connected in layer combination. [3] Gaur, A., Tiwari, S., Kumar, C., &
• Single device in flexing/fix free Maiti, P. (2020). Bio-waste orange peel
arrangement shows maximum voltage of and polymer hybrid for efficient energy
pk-pk 38.0V. harvesting. Energy Reports, 6, 490-496.
• Device shows promising and linear [4] Han, Y., Han, Y., Zhang, X., Li, L.,
results in impact testing in single device Zhang, C., Liu, J., ... & Huang, W.
mode and double device mode with (2020). Fish gelatin based triboelectric
maximum voltage up to 67.2 V and nanogenerator for harvesting
induced charge of 66.19 nC at applied biomechanical energy and self-powered
maximum load of 2.88 N. sensing of human physiological
• During CAM based vibration testing signals. ACS applied materials &
machine tests , single device shown interfaces, 12(14), 16442-16450.
higher results at 122 Hz frequency with [5] Kaur, J., & Singh, H. (2020). Fabrication
maximum open circuit voltage as 37.2 V. of composite material based
• High flexibility, and durability of nanogenerator for electricity generation
commercial device makes it a good enhancement of food waste by-
choice for application of self-powered product. Materials Chemistry and
wearable devices. Physics, 256, 123331.
• Energy harvester with multi-layer of [6] Karan, S. K., Maiti, S., Paria, S., Maitra,
devices can be used to power up low A., Si, S. K., Kim, J. K., & Khatua, B. B.
power energy consuming devices (2018). A new insight towards eggshell
utilizing vehicle movement induced membrane as high energy conversion
wind, ocean waves energy, machines efficient bio-piezoelectric energy
vibration, public movement pathways harvester. Materials today energy, 9,
etc. 114-125.
164
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
165
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ABSTRACT
Welding is one of the major processes in the automobile, construction, aviation, and other sectors.
However, the inefficiencies, high failure rates, and poor weld quality of manual welding methods
are common. This advocates digitization of the welding process which can be named Welding for
Industry 4.0. Digitization of manual welding to digital welding includes real-time weld monitoring
and control. Accurate weld defect classification is a challenging task for efficient weld monitoring.
Different machine-learning algorithms have been developed for the classification of welding
defects. This paper presents a comparative analysis of CNN and ResNet-based approaches for weld
defect classification when using vision-based monitoring. The dataset used for the research
includes more than 30,000 images. Six different types of defects were considered including burn-
through, contamination, lack of fusion, misalignment, and lack of penetration. Both the models
were trained using the dataset to classify these defects. The accuracy of the algorithms was
determined based on the number of images taken to train the model. Comparative analysis revealed
that CNN achieved a higher precision compared to ResNet. The paper highlighted the effectiveness
of CNN employed for weld defect classification for the development of an Industry 4.0-based
robotic welding system.
Keywords: Welding Digitization; Industry 4.0; Machine Learning; CNN; ResNet
and control systems for welding processes.
1. INTRODUCTION The review identified 101 papers that met the
Deep learning models in recent years are inclusion criteria, and it found that machine
summarized and compared with a detailed learning is a promising approach for
discussion of several typical networks in the improving the quality of welds.
field of image classification, object detection, Elsheikh [4] highlights the significance of
and segmentation [1]. To implement the machine learning techniques in predicting
feature extraction and classification in one joint properties, enabling real-time control,
algorithm and to implement the overall and diagnosing tool failures in the welding
automation, oh. et al. [2] propose a method of process. The review discusses the utilization
automatically detecting welding defects of different machine learning algorithms and
using Faster R-CNN which is a deep learning models, emphasizing their effectiveness in
basis. Breitenbach et al. [3] conducted a improving the quality and efficiency of
systematic literature review on machine friction stir welding.
learning approaches for quality monitoring
166
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Zhou et al. [5] presented a case study on the improving the extraction accuracy and
use of machine learning and semantics to processing speed of the algorithm [10].
predict the quality of automated welding. The Despite being in use for years, welding robots
case study was conducted at Bosch, and it are pre-programmed devices with little to no
showed that machine learning can be used to intelligence. Figure 1 shows that Intelligent
predict the quality of welds with high manufacturing is implemented by the
accuracy. Asif et al. [6] proposed a machine- transformation of Human physical systems to
learning model to predict welding quality Human cyber-physical systems to achieve
using air-coupled acoustic emission (CAE) specific manufacturing goals at an optimized
and weld inputs. The model was trained on a level [11].
dataset of CAE signals and weld quality data,
and it was shown to be able to predict
Human Human
welding quality with high accuracy. Wang et
Market
Demand and
Competition
HPS
HCPS
recognizing weld defects based on visible
spectral images and machine learning. The Physical Technology
Machine
Learning
Physical
System
the weld, and a machine learning algorithm Manual Welding Intelligent Welding
167
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
propose different strategies for improving the that bypass some layers, ResNet enables
computational efficiency of the deep- information to flow directly through the
learning models adopted in reinforcement- network, resulting in faster learning and
learning (RL) scenarios. The robustness of better performance on various image
most image processing algorithms is processing tasks[20], [21].
deficient during welding practice, and the He et al. [22] introduced a new approach to
security regime for tracking welding is not classification and object detection tasks,
considered in most trajectory recognition and revolutionizing the field of image processing.
control algorithms. For these two problems, Unlike traditional methods, ResNet addresses
an adaptive feature extraction algorithm was the vanishing gradient problem which is a
proposed, which can accurately extract the bottleneck in training deep neural networks.
seam centre from the continuous, By including "shortcut connections" that
discontinuous, or fluctuating laser stripes bypass certain layers, ResNet allows for the
identified and located by the CNN model, direct flow of information, facilitating
while the prior model can quickly remove a efficient learning and achieving significant
large amount of noise and interference except performance improvements across various
the stripes, greatly improving the extraction image processing applications.
accuracy and processing speed of the
algorithm [10]. To implement the feature Algorithm Advantages Purpose
extraction and classification in one algorithm Artificial capabilities of Solve intricate
Neural performing problems using
and to implement the overall automation, oh Network sophisticated, internal
et al. [2] propose a method of automatically (ANN) non-linear mechanisms.
mappings,
detecting welding defects using Faster R- learning
CNN which is a deep learning basis. Günther independently,
and generalizing
et al. [18] suggest that deep learning features knowledge.
and general-value-function predictions can Convolutional Decreased Inspection of
Neural parameter count surface
be beneficially combined with actor–critic Network and maintained integration and
(CNN) invariance to identification of
reinforcement learning to learn context- changes in shift, machinery
appropriate control policies to govern scale, and faults.
distortion.
welding power in real-time. Machine ResNet Shortcut Ultimately
learning algorithms can be exploited to connections that enabling deeper
alleviate and more
predict the weld quality during laser welding vanishing accurate image
with high accuracy [19]. gradients processing
models.
In the field of image processing, the ResNet Support vector Less sample, Classification
algorithm has been a breakthrough in tasks machine global optimal and regression
(SVM) analysis
like classification and object detection. Deep Belief The ability to Can be trained
Unlike conventional methods, ResNet Network train layer by layer-wisely to
(DBN) layer for be more
addresses the problem of vanishing gradient improved efficiently
which arises while training deep neural efficiency.
networks. By using "shortcut connections"
168
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
information science and technology are convergence of big data, AI, and human
enabling the transformation of traditional cyber-physical systems, the future of welding
welding to intelligent one [30]. Because of promises enhanced efficiency, quality, and
advances in computer science and artificial adaptability. Figure 2 compares two methods
intelligence, intelligent automation is for image feature extraction: conventional
replacing manual work. Manufacturing methods and Convolutional Neural Networks
research programs have all examined smart (CNNs).
manufacturing, Industry 4.0, and the Internet
of Things. In the future big data, intelligent Image feature extraction using
conventional methods.
manufacturing, human cyber-physical
systems, and AI will revolutionize industry. Image Filtering and Image
[11], [31], [32]. Real time Enhancement Denoising segmentation
2. METHODOLOGY
We chose two different Machine learning CNN On-line weld
Resize (Convolutional defects
algorithms CNN and ResNet to train it on the Neural Network) classification
specific dataset and get a comparative
analysis of both the algorithms. Conventional
Figure 2. Image Feature Extraction:
methods require hand-crafted design, while Conventional Methods vs.
CNNs learn features for weld defect Convolutional Neural Networks
classification, potentially offering a more
robust approach. To address challenges such
as defect detection and the need for robust 3. ALGORITHM MODELLING
prediction models, a comparative analysis of Convolutional Neural Network (CNN)
two different machine learning algorithms The following steps were followed from data
has been done, and these algorithms were annotation to training the ML algorithms.
trained on the dataset which includes more
than 30,000 images. Through the
170
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
material.
5. Continuously improve the algorithm to
get the maximum accuracy. 6 Comparative Analysis
171
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
172
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
178
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ABSTRACT
With the rise in operational parameters, structural elements within power plants are exposed to
heightened temperatures and pressures, reaching levels conducive to creep. This phenomenon
poses the risk of component failure and fracture. Over the preceding decades, substantial efforts
have been dedicated to comprehending the fundamental mechanisms of creep. Significantly, there
has been a focus on predicting the lifespan of equipment to ensure safety and reliability in plant
operations. This article introduces a modified model, building upon the Norton-Bailey and
KachanovRabotnov constitutive models, aiming to comprehensively characterize all three stages
of creep (primary, secondary, and tertiary) specifically for stainless steel material. The proposed
model undergoes numerical calculations to simulate damage development in thin-hole steel plates.
Special attention is given to the influence of specimen dimensions and stress on damage evolution.
The outcomes of the study reveal the efficacy of the novel model in generating complete creep
curves for various materials, emphasizing the notable impact of stress on creep behavior and
damage progression. Sensitivity studies, employing Response Surface Methodology (RSM) and
Analysis of Variance (ANOVA), are conducted to quantitatively and qualitatively analyze the new
creep model data in comparison to established models. The objective of these sensitivity studies is
to assess the significant contributions of various factors to the target response strain in this specific
case study.
Keywords—creep deformation, curve fitting, Kachanov-Rabotnov model, damage evolution,
Norton Bailey model
1. INTRODUCTION manifesting timedependent, inelastic
Numerous creep models have been deformation, particularly effective in
developed over the years to anticipate the anticipating material behavior in the
creep behavior of materials; however, secondary creep regime.
these models exhibit certain limitations In contrast, the Omega model, introduced by
[1]. The Norton–Bailey (NB) model, also the Material Properties Council and
known as Norton’s power law, serves as a proposed by Prager [3] in 1995, stands out as
foundational model integrated into the a reliable prediction model. Recognized for
finite element Abaqus software, widely its simplicity and lower dependence on
employed for creep analysis in material constants, the Omega model boasts
conjunction with other models developed a well-documented creep evaluation process,
by Bailey and Norton [2]. Functioning as demonstrating excellent performance in
a benchmark, the NB model predicts the establishing property relations across diverse
creep deformation behavior of materials materials. Primarily designed for calculating
179
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
the remaining life of components operating with the secondary and tertiary stages of
in the creep zone under high temperatures creep but lacks consideration for primary
and pressures, the Omega model utilizes a creep accumulation.
strain rate parameter and a multi-axial
damage parameter, as outlined by Yeom et To address these limitations, a modified
al. [4]. It effectively forecasts the rate of creep model is proposed in this study, aiming
strain accumulation, creep damage to comprehensively describe all three stages
accumulation, and the remaining time to of creep. The modified constitutive creep
failure based on the stress state and model is derived through a combination of
temperature, covering both primary and the Norton-Bailey and Kachanov-Rabotnov
secondary creep regime deformations and models. Numerical calculations using this
excelling in predicting material rupture time modified model were conducted to explore
at lower temperatures [5]. the creep behavior and damage development
The Kachanov–Rabotnov (KR) model, an in a thin-hole steel plate made of SS-304
early implementation of the continuum material under constant stresses. The
damage mechanics (CDM) approach for investigation places particular emphasis on
creep proposed by Kachanov and Rabotnov understanding the impact of specimen
[6], is another significant model. This model dimensions and stresses on damage
allows for the modeling of secondary and development, deformation, and rupture life.
tertiary creep deformations through a set of This paper addresses the limitations and
coupled equations. Continuous efforts have deficiencies of the NortonBailey and
been dedicated to enhancing the KR law, Kachanov-Rabotnov models, paving the way
leading to recent variations aimed at for a more comprehensive understanding of
generating contour deformation maps [7]. creep behavior in structural components
Stewart and Gordon [8] are actively subjected to various stress conditions.
developing methods to define transversely
isotropic creep damage properties and 1.1 Problem Statement
estimate stress-independent tertiary creep The Norton-Bailey model is limited in that it
damage constants using both strain- and exclusively addresses the secondary creep
damage-based analytical approaches. These regime and lacks predictions for the primary
advancements contribute to the ongoing and tertiary regimes, resulting in an overall
refinement of creep models for a more error when these stages dominate. Although
comprehensive understanding of material the KachanovRabotnov (KR) model shows
behavior under varying conditions. promise, it involves a substantial number of
The Norton-Bailey constitutive model material constants, and its formulation
proves effective in delineating the neglects consideration for the primary creep
relationship between strain and stress during regime during analysis. While KR
successfully emulates both continuum creep
the primary and secondary stages of creep.
damage and discontinuous plastic damage at
However, when analyzing the life of
rupture within a continuous function, the
structural components controlled by rupture,
intricacies of the model pose challenges for
it becomes imperative to scrutinize their
integration into finite element (FE) analysis.
creep behavior in the tertiary region of the Both models have made strides in addressing
creep curve. Conversely, the Kachanov- nucleation and plasticity for crack initiation
Rabotnov constitutive model aligns well and growth to some extent. However,
180
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
accurate detection and identification of creep complications with the proposed model. The
crack growth remain limited. The new material model provides an alternative
mechanism through which creep cavities to the current models, producing results
nucleate is not well-established, despite the more effectively.
presence of voids in every material. 3. A simplified formulation is expected,
Observations suggest that cavities tend to to accurately consider all three creep stages
nucleate along grain boundaries under (I-III), requires only minimal empirical data
constant load for extended periods. Creep and containing minimal number of material
cavitation is a critical aspect, often given less fitting parameters. The model will offer ease
priority when addressing creep crack of applicability and versatility to deal with
initiation and growth. Given these material’s data which is already scarce for
limitations, the article proposes a new creep the three creep stages.
prediction model based on creep power laws. The next section discussed the theoretical
This model aims to overcome the challenges framework of established models and
posed by existing models and address the mathematical formulation of the new creep
following research questions regarding how material model.
the material behaves during creep
deformation under various parameter 2. THEORETICAL FRAMEWORK
variations, conditions, and circumstances. 2.1 Norton Bailey Model
Norton's power law (1929), rooted in
1.2 Research Questions the Arrhenius rate equation [10], stands as
1. How the new model can improve the the most renowned and widely used
prediction of creep failure at high minimum creep strain rate law. This law,
temperatures? represented in Equation (1), is a prominent
2. Can the model accurately predict all framework in understanding and describing
creep stages at crack initiation, growth and creep behavior.
rupture? 𝑐𝑟̇ = 𝐵’𝜎𝑛 𝑒𝑥𝑝 (−𝑄𝑐⁄𝑅𝑇)𝑚 , (1)
3. How much improvement in efficiency of
the analysis and accuracy of the results In the given context, 𝑐𝑟̇ signifies the minimum
increase using the newly proposed model? creep strain rate, B’ is the material constant,
σ denotes the applied stress, n represents the
The new creep prediction model is based on power law exponent, 𝑄𝑐 represents the
the following hypothesis: activation energy, R stands for the universal
gas constant, and T corresponds to the
applied temperature. Equation (1) is
1.3 Hypothesis
subsequently streamlined when considering
1. The prediction of creep behaviour of
a constant temperature, resulting in the
the mechanical equipment, especially due to
simplified form represented by Equation (2).
generation of cracks through nucleation and
voids within the material will be attained 𝑐𝑟̇ = 𝐴 𝜎𝑛 𝑡𝑚 , (2)
through the new model, with the
without considering time, Equation (3)
shortcomings of the current models are
becomes:
addressed.
𝑐𝑟̇= 𝐴 𝜎𝑛 , (3)
2. Prediction of the creep behaviour of
where,
the material can be handled with less
181
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
, (12)
t( , (7)
, (13)
, (8)
, (14)
Here, t denotes the current time, and
represents the current damage. These
where, is DECRA(1), uniaxial deviatoric
equations offer a means to calculate the
creep strain increment in ABAQUS creep
current time based on the current damage and
user sub-routine obtained from ABAQUS
stress, or alternatively, the current damage
documentation 3.4.6 as in Equation 14[15].
based on the current time and stress. Stewart
Equation (11) can also be differentiated with
and Gordon [12] introduced two valuable
respect to as:
techniques for analyzing and applying the
182
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
. (15)
(20)
where, is DECRA(5), von Mises stress
in ABAQUS creep user sub-routine.
Equation (14) and Equation (15) can be (21)
implemented in creep sub-routine to
incorporate the proposed creep model in
ABAQUS for time hardening and strain In this equation, x and σ serve as independent
hardening methods. variables, while y and ε ̇ act as dependent
variables, and n' denotes the number of
The other method of integration adapted for samples. In the regression analysis [49],
the new creep model in Abaqus is by curve stress is considered an independent variable,
fitting: and strain rate is the dependent variable. The
The general power law regression is selection of the stress range is influenced by
expressed in Equation (16), while the the precision of the regression, with a larger
Norton-Bailey equation is provided in stress range in the sample data yielding a
Equation (17) for the purpose of comparison more accurate curve fit [19].
[16]: The creep parameter and stress exponent, as
outlined in Equations (19) and (21), undergo
(16) modifications for the curve fitting of the new
model through regression into the
=A (17) NortonBailey model in Abaqus, resulting in
Equations (22) and (23). In these equations,
In this context, the equation involves the the constant b is introduced to signify the
criterion variable and prediction response damage evolution parameter, aiding in
(y), where A' represents the curve coefficient, modeling material deformation during the
and B is the exponent of the predictor tertiary creep stage [20].
variable (x). This equation is juxtaposed with
the Norton-Bailey power law regression (22)
presented in Equation (17). The parameter B
in Equation (18) is then compared with the
derived stress exponent (n) in Equation (19) (23)
for creep [17].
where,
(18) b ;
ω = damage evolution parameter;
(19) = sample size
Likewise, the curve coefficient A' in The new model equation becomes:
Equation (20) was contrasted with the
Norton-Bailey creep parameter A in =A *b (24)
Equation (21) to determine the parameters
[18]. 3. METHODOLOGY
183
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
185
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
186
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
where, A = stress, B = stress exponent, C = creep primary, secondary, and tertiary stages of
parameter and D = damage parameter. creep for the material. The obtained results
The provided equation is employed to were derived through the utilization of a
predict the response : curve-fitting method.
Figure 2 depicts the von Mises stress and
relaxed stresses observed during the
In the provided equation, y signifies the simulation employing the new model, which
predicted value of the response. 𝑋𝑖 and 𝑋𝑗 incorporates the creep parameter and stress
represent independent parameters, while 𝛽0 , exponent. The performance of the new model
𝛽𝑖 and 𝛽𝑖𝑖 denote the constant, linear, and was also compared against the Omega model,
interactive coefficients, respectively. following a similar approach to that used for
Additionally, 𝛽𝑖𝑗 stands for the quadratic established models. Creep strain rate results
coefficient, and n represents the number of for the new model were obtained through a
factors involved in the analysis. curve fitting process. Figure 2 provides a
The design matrices for the new model were visual representation of the induced von
constructed using the central composite Mises stress and creep strain results [29].
design (CCD), incorporating four
independent design factors: stress (A), stress
exponent (B), creep parameter (C), and
damage parameter (D). The target response,
in this case, is the strain rate, and the
corresponding design matrix is reported in
Table 4 [27]. The values in Table 4 were
determined analytically to assess the
significant contribution of these factors to the
target response. For each model and the new
model, actual and predicted values of the
finite element strain rate were obtained from
the Response Surface Methodology (RSM)
and Analysis of Variance deemed Figure 2: Von Mises stress and Relaxed
significant, as evidenced by a coefficient of stress distribution with Omega- New
determination (R²) value exceeding 80%. Model’s Regression -Visco-elastic
4. RESULTS plastic run-time of 1000 h, 101 MPa and
625 0C.
The recently proposed model, which
integrates the Norton Bailey and Kachanov- 4.1 New Creep Model Results by
Rabotnov models, was applied within the FE Curve-Fitting
package Abaqus to derive creep strain and The creep strain rate results were obtained
strain rate results for an SS-304 thin-hole while implementing new model in Abaqus
plate using finite element analysis. Consistent by the method of curve fitting, used for
with established models in the FE package, deriving constants in the equation of new
the same predefined boundary conditions model. The new model was standardized by
were employed. The simulation extended the Omega model with the help of
over 1000 hours at a temperature of 625°C, regression analysis. Figure 3(a) is depicting
yielding promising outcomes that covered the
188
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
the comparison of creep strain rates for the Figure 3(b): New Model Creep strain
new model with the Omega model on the comparison with the Omega Model
same applied conditions. New creep model Creep strain by Curve Fitting
curve is covering primary, secondary and
tertiary creep stages, whereas Omega model To capture the tertiary creep behaviour of
can only cover the primary and secondary stainless-steel material, arbitrary damage
creep stages. Omega model gives error in evolution parameter (ω) values were applied.
predicting creep strain at elevated The strain to reach the particular value of (ω)
temperatures, this deficiency can be does not show the similar trend. As depicted
fulfilled by using the new creep model [30]. in the Figure 4, a marginal decrease in strain
Omega model is unable to predict creep to reach the particular value of (ω), resulted
curves for the large exponential data, the in the reduction of the applied stresses. The
deficiency can be fulfilled by using the new kinetics evolution of coupled strain and
proposed creep model. Figure 3(b) is damage does not necessarily follow the same
depicting the results comparison of the new
trend. It is well known fact, that any material
model with the Omega model for creep
fails if (ω), reaches critical damage 𝜔𝑐 . The
strain at same pre-defined boundary
material fails usually before reaching the
conditions [31].
critical damage of value equivalent to 1. The
new creep model has the ability to model the
tertiary creep damage at all damage
parameter (ω) values [32].
Time, t (h)
190
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
1 0 1.48
A, B, AND C, pp. 663–670, 2012, doi: Crack Growth Prediction Model for the
10.1115/IMECE201288389. Life Assessment of
[9] C. M. Stewart and A. P. Gordon, “Strain [16] Stainless-Steel Material Using
and Damage-based Analytical Methods Computational Modeling,” Metals
to Determine the Kachanov-Rabotnov (Basel)., vol. 13, no. 11, p. 1854, 2023,
Tertiary Creep-Damage Constants,” Int. doi: 10.3390/met13111854. [15]
J. Damage Mech., vol. 21, no. 8, pp. ABAQUS, “ABAQUS 6.12 User’s
1186–1201, 2012, doi: documentation 3.4.6,” vol. IV. 2013.
10.1177/1056789511430519. [17] M. Sattar, A. R. Othman, M. F. Othman,
[10] Y. Chang and H. Xu, “the Damage and M. F. Musa, “Regression Analysis
Development and Failure Description of Omega Model to Norton- Bailey Law
Under Multiaxial Creep,” Proc. ASME for Creep Prediction in Fitness for
2014 Press. Vessel. Pip., pp. 1–6, 2015. Service Assessment of Steel Material,”
[11] Z. Abdallah, V. Gray, M. Whittaker, and Solid State Technol., 2020.
K. Perkins, “A Critical Analysis of the [18] M. K. Booker, “Use of Generalized
Conventionally Employed Creep Lifing Regression Models for the Analysis of
Methods,” Materials (Basel)., vol. 7, no. Stress-Rupture Data.,” pp. 459–499,
5, pp. 3371–3398, 2014, doi: 1978.
10.3390/ma7053371. [19] B. Xiao, L. Xu, L. Zhao, H. Jing, Y. Han,
[12] D. L. May, A. P. Gordon, and D. S. and Y. Zhang, “Creep Properties, Creep
Segletes, “The Application of the Deformation Behavior, and
Norton-Bailey Law for Creep Prediction Microstructural Evolution of 9Cr-3W-
through Power Law Regression,” in 3Co-1CuVNbB Martensite Ferritic
Proceedings of the ASME Turbo Expo, Steel,” Mater. Sci. Eng. A, vol. 711, no.
2013, vol. 7 A, pp. 1–8, doi: April 2019, pp.434–447,2018,doi:
10.1115/GT2013-96008. 10.1016/j.msea.2017.11.061. M. Sattar,
[13] C. M. Stewart and A. P. Gordon, A. R. Othman, S. Kamaruddin, M. A.
“Analytical Method To Determine the Alam, and M. Azeem, “Creep
Tertiary Creep Damage Constants of the Parameters
Kachanov-Rabotnov Constitutive [20] Determination by Omega Model to
Model,” in ASME, International Norton Bailey Law by Regression
Mechanical Engineering Congress & Analysis for Austenitic Steel SS-304,”
Exposition IMECE2010, 2010, pp. 1–8. Solid State Phenom., vol. 324, pp. 188–
[14] C. M. Stewart, A. P. Gordon, E. A. 197, 2021.
Hogan, and A. Saxena, [21] D. L. May, A. P. Gordon, and D. S.
“Characterization of the Creep Segletes, “The application of the norton-
Deformation and Rupture Behavior of bailey law for creep prediction through
DS GTD-111 using the Kachanov- power law regression,” Proc. ASME
Rabotnov Constitutive Model,” J. Eng. Turbo Expo, vol. 7 A, pp. 1–8, 2013,
Mater. Technol. Trans. ASME, vol. 133, doi:10.1115/GT2013-96008.
no. 2, 2011, doi: 10.1115/1.4003111. [22] ASME, American Petroleum Institute
[15] M. Sattar, A. R. Othman, M. F. Othman, API-579, Fitness for Service, Operation
H. T. Ali, and M. K. Khan, “New Creep Manual, 3rd ed., no. FFS-1. Washington
194
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
DC, USA: The American Society of [29] M. A. Alam et al., “Modelling and
Mechanical Engineers, 2016. Optimisation of Hardness Behaviour of
[23] A.P. V. ASME Codes, “ASME. (2015). Sintered Al/SiC Composites using RSM
ASME Boiler and Pressure Vessel Code and ANN: A Comparative Study,” J.
An International Code- section II part A, Mater. Res. Technol., vol. 9, no. 6, pp.
1998.,” Oper. Man., p. 1998, 2015, doi: 14036–14050, 2020, doi:
http://dx.doi.org/10.1016/B9780323035 10.1016/j.jmrt.2020.09.087.
06-4.10361-X. [30] J. Jadoon, A. Shazad, M. Muzamil, M.
[24] A.A. Al-Bakri, Z. Sajuri, A. K. Ariffin, Akhtar, and M. Sattar, “Finite Element
M. A. Razzaq, and M. S. Fafmin, Analysis of Composite Pressure Vessel
“Tensile and Fracture Behaviour of very Using Reduced Models,” Tecciencia,
thin 304 Stainless Steel Sheet,” J. vol. 17, no. 33, pp. 49–62, 2022,
Teknol., vol. 78, no. 6–9, pp. 45–50, doi:10.18180/tecciencia.2022.33.5.
2016, doi:10.11113/jt.v78.9146. [31] M. Sattar, A. R. Othman, M. Muzamil,
[25] L. M. Powers, S. M. Arnold, and A. S. Kamaruddin, M. Akhtar, and R. Khan,
Baranski, “Using ABAQUS Scripting “Correlation Analysis of Established
Interface for Materials Evaluation and Creep Failure Models through
Life Prediction,” ABAQUS Users’ Conf., Computational Modelling for SS-304
pp. 1–11, 2006. Material,” Metals (Basel)., vol. 13, no. 2,
[26] D. P. Jones, J. L. Gordon, D. N. Hutula, p. 197, 2023, doi:
D. Banas, and J. B. Newman, “An 10.3390/met13020197.
elastic-perfectly plastic flow model for [32] M. Sattar, A. R. Othman, M. Akhtar, S.
finite element analysis of perforated Kamaruddin, R. Khan, and F. Masood,
materials,” J. Press. Vessel Technol. “Curve Fitting for Damage Evolution
Trans. ASME, vol. 123, no. 3, pp. 265– through Regression
270, 2001, doi:10.1115/1.1357538. [33] Analysis for the Kachanov – Rabotnov
[27] D. B. Alemayehu, S. J. Huang, and E. G. Model to the Norton – Bailey Creep Law
Koricho, “Experimental and FEM of SS-316 Material,” Materials (Basel).,
Analysis of three Carbon steel vol. 14, no. September, 2021, doi:
Characterization under Quasi-Static 10.3390/ma14195518.
Strain Rate for Bumper Beam [34] J. Christopher, C. Praveen, V. Ganesan,
Application,” MATEC Web Conf., vol. G. P. Reddy, and S. K. Albert,
123, pp. 0–5, 2017, doi: “Influence of Varying Nitrogen on
10.1051/matecconf/201712300019. Creep Deformation and Damage
[28] A.Czyrski and H. Jarzebski, “Response Behaviour of Type 316L in the
surface methodology as a useful tool for Framework of Continuum Damage
evaluation of the recovery of the Mechanics Approach,” Int.
fluoroquinolones from plasma-the study J.10.1177/1056789520939622.
on applicability of box-behnken design, [35] L.Esposito and N. Bonora, “Time-
central composite design and doehlert Independent Formulation for Creep
design,” Processes, vol. 8, no. 4, 2020, Damage Modeling in Metals based on
doi:10.3390/PR8040473. Void and Crack Evolution,” Mater. Sci.
Eng. A, vol. 510–511, no. C, pp. 207–
195
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
196
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*[email protected]
ABSTRACT
This paper aims to investigate the effect of Gas-Assisted Direct Thermal Method (GA-DTM)
processing parameters on hypoeutectic Al-Si alloy feedstock billet microstructure. Several
methods have been used to produce globular microstructures, that are suitable for semisolid metal
processing (SSMP) in these recent years. Factors that affected microstructure such as pouring
temperature, holding time, and magnesium addition were analysed. A combination of parameters,
consists of pouring temperatures of 590 and 610 °C, holding time of 10 and 20 s, and magnesium
addition of 0.5 and 1.5 wt% were used. The microstructures of the prepared sample were
investigated. The results showed that hypoeutectic Al-Si alloy feedstock billets prepared with GA-
DTM with 610 °C pouring temperature, 20 s holding time and 1.5 wt.% magnesium addition
produced finer microstructure. The grain size, circularity, and aspect ratio, for the formed
microstructure were at 45.70 µm, 0.69, and 1.42 respectively. The microstructure with small and
uniform primary particles is commonly preferred for SSMP. The study provide insight into the use,
performance, and applicability of hypoeutectic Al-Si alloy feedstock billets in semi-solid
processing applications such as thixoforming.
Keywords: Grain refinement, Hypoeutectic Al-Si Alloy, Magnesium Addition, Gas-Assisted
Direct Thermal Method, Feedstock Billet, Microstructure.
refining aluminium grain. Grain refinement found that the generation of finer primary
by vibration and stirring during grains with a spheroidal shape was improved
solidification, rapid solidification, grain by 36.4% when DTM was combined with
refiner addition, and extreme plastic external colling gas [6]. Furthermore, several
deformation [3]. studies on the microstructure and hardness
Semi solid metal process (SSMP) has evolution of aluminium alloy 6061 indicated
recently been applied in the industry because that the casting temperature and holding time
it produced globular with refined had a greater influence on structure formation
microstructure, a crucial element for SSMP [8], while another study specifically reported
for its use in engine blocks and aerospace that a combination of low casting
parts. The final SSMP product also has minor temperature and long holding time produced
a global microstructure [9].
defects such as porosity, shrinkage, gas
Hypoeutectic Al-Si alloys, a subset of
entrapment, and macro segregation [4]. In
aluminium-silicon alloys, possess distinctive
addition, the SSMP technique develops characteristics that make them valuable in a
aluminium alloys to a semisolid state and variety of industrial applications [10]. These
improves the characteristics of aluminium alloys are characterized by having a silicon
alloys. The SSMP of die casting seems to content below the eutectic composition of
meet several objectives since it may offer the approximately 12.6% silicon by weight.
industry a means of creating mass- Some key features of hypoeutectic Al-Si
production, high-performance, and complex alloys include their excellent fluidity, which
components close to their net shape while makes them well-suited for intricate casting,
simultaneously reducing energy and good castability, enabling them to be
consumption and waste [5]. The success of easily shaped into a wide range of
the SSMP is merely dependent on the components through various casting methods
[11]. However, their mechanical strength is
globular and refined microstructure
incredible. As magnesium can improve
feedstock for thixoforming. In the literature, mechanical properties and performance, its
many efforts have been made to explore the addition to aluminium-silicon alloys has
ability to produce more globular and finer attracted attention [1].
microstructure. Increasing the addition of magnesium has a
positive impact on alloy strength. Under
Various methods have been tried to prepare
conventional casting conditions, it results in
the material with improved globular
the formation of a coarse Mg₂ Si alloy,
microstructure for SSMP. One method is the
which, unfortunately, leads to alloys with
Gas-Assisted Direct Thermal Method (GA-
limited ductility and poor fracture toughness
DTM), which modifies the thermal
[12]. Therefore, the coarse Mg₂ Si phase
characteristics of the liquid alloy to delay the
offers good mechanical strength but requires
material's solidification [6]. The
modification to maintain ductility.
microstructure's shape has been enhanced
Augmenting the magnesium content
using GA-DTM, resulting in a more globular
contributes to the enhancement of the
and refined microstructure. It is also known
microstructure in Al-Si alloys [13]. To
that adding magnesium to Al-Si-Cu alloys
address this issue, SSMP is utilized to create
through a process known as grain refining
Al-Si alloys with a finer microstructure and
enhanced the alloys' mechanical properties
decreased segregation [14]. In their study,
by causing a stronger Mg₂Si phase to form
Cory et al. investigated the effect of
[7]. A study by Azhani Abd Razak et al.
198
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
199
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
200
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
0.8 2
0.6
Aspect Ratio
1.5
Circularity
0.4 1
0.2 0.5
0
0
1 2 3 4
1 2 3 4
Sample Number
Sample Number
Figure 2. Average circularity for sample
number 1 to 4.
Figure 3. Average aspect ratio for sample
number 1 to 4.
Each of the 4 samples had a circularity
The aspect ratio obtained by each of the 4
ranging from 0.61 to 0.69. The magnesium
samples ranged between 1.42 and 1.64. The
addition has a most significant influence on
addition of magnesium contributed to the
circularity. The circularity value increases as
effect on the aspect ratio. Circularity
the holding time increases [16],[20].
increases as the magnesium addition
Circularity, measured on a scale from 0.0 to
increases. The aspect ratio represents the
1.0, is the degree to which the shape of the
ratio of the main axis to the minor axis. The
microstructure resembles an ideal circle [21].
higher the aspect ratio size, the longer the
Sample 4 had the highest circularity value
particle. It was discovered that the addition of
(0.69). A greater magnesium value causes the
Mg resulted in a reduction in the particle size
circularity to increase at any given
of the sample. Among the hypoeutectic Al-Si
magnesium addition rate. A value close to 1.0
alloy feedstock billets prepared by the GA-
means that the shape is considered perfectly
DTM method, the sample prepared at 610 °C
spherical. Nevertheless, particles become
pouring temperature, 20 s holding time, and
more elongated as the value approaches zero.
1.5 weight percent magnesium addition had a
While pouring temperature and holding time
small and spherical microstructure. The
can influence the cooling rate and subsequent
grain size, circularity, and aspect ratio are
microstructure, their effect on grain size may
45.70µm, 0.69, and 1.42, respectively. This is
be overshadowed by magnesium's strong
consistent with the recommended dosages for
grain refining effect. In some cases, changes
globular microstructure [21],[22].
in pouring temperature and holding time may
have a secondary or complementary impact
on grain size, but the primary influence 4. CONCLUSIONS
appears to be from magnesium addition. The The contribution of processing parameters in
aspect ratio obtained from microstructure the preparation of hypoeutectic Al-Si alloy
investigation of billets prepared by GA-DTM feedstock billet by the GA-DTM method was
is presented in Figure. successfully investigated. The study results
provided insight into the effect of grain
refinement and processing parameters on the
microstructure of the alloys, which allowed
us to know pouring temperature, time, and
magnesium addition and their combined
effects on grain refinement. At 610 °C
201
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
pouring temperature, 20 s holding time and [6] N.A. Abd Razak, A.H. Ahmad, M.M.
1.5 wt.% magnesium added billet sample had Rashidi, and S. Naher, “An investigation
very small and uniform primary grains. At of semisolid Al7075 feedstock billet
these combination parameters, the grain size, produced by a gas-assisted direct thermal
circularity, and aspect ratio are 45.70µm, method.,” International Journal of
0.69, and 1.42, respectively. These study Advanced Manufacturing Technology.
findings will help to improve the mechanical vol. 114, no. 3–4, pp. 1233–1240, 2021.
strength, applicability, and microstructure of [7] Al-Si Cast Alloys - Microstructure and
materials produced by semi-solid processing. Mechanical Properties at Ambient and
Elevated Temperature Al-Si Cast Alloys
ACKNOWLEDGEMENT - Microstructure and Mechanical
The authors would like to thank the Ministry Properties at Ambient and Elevated
of Higher Education for providing financial Temperature., 2015.
support under Fundamental research grant [8] B. Benjunior, A.H. Ahmad, and M.M.
No.FRGS/1/2019/TK03/UMP/02/8Universit Rashidi, “Direct thermal method pouring
y reference RDU1901122) and Universiti temperature and holding time effect on
Malaysia Pahang Al-Sultan Abdullah for aluminium alloy 6061 microstructure.,”
laboratory facilities as well as additional IOP Conference Series: Materials
financial support under Internal Research Science and Engineering. vol. 788, no. 1,
grants RDU160311 and RDU603125. p. 2020.
[9] N.A. Razak, A.H. Ahmad, and M.M.
REFERENCES Rashidi, “Investigation of pouring
[1] D. Watson, “Microstructure and temperature and holding time for
Mechanical Properties of Ductile Die- semisolid metal feedstock production.,”
cast Al-Mg-Si-Mn Alloys.,” no. October, IOP Conference Series: Materials
p. 2015. Science and Engineering. vol. 257, no. 1,
[2] J. Wannasin, R.A. Martinez, and M.C. p. 2017.
Flemings, “Grain refinement of an [10]J. Xu, Y. Li, K. Ma, et al., “In-situ
aluminum alloy by introducing gas observation of grain refinement
bubbles during solidification.,” Scripta dynamics of hypoeutectic Al-Si alloy
Materialia. vol. 55, no. 2, pp. 115–118, inoculated by Al-Ti-Nb-B alloy Scripta
2006. Materialia In-situ observation of grain
[3] M. He, Z. Zhang, W. Mao, B. Li, Y. Bai, refinement dynamics of hypoeutectic
and J. Xu, “Numerical and experimental Al-Si alloy inoculated by Al-Ti-Nb-B
study on melt treatment for large-volume alloy.,” no. June, p. 2020.
7075 alloy by a modified annular [11]M.S. Prabhudev and C.T. Mata,
electromagnetic stirring.,” Materials. vol. “Studies on the microstructure and
12, no. 5, p. 2019. mechanical properties of a356 alloy with
[4] M. Modigell, A. Pola, and M. Tocci, minor additions of copper and
“Rheological characterization of semi- magnesium.,” no. August 2009, pp. 1–5,
solid metals: A review.,” Metals. vol. 8, 2015.
no. 4, pp. 1–23, 2018. [12]M.F.M. Tajudin, A.H. Ahmad, and
[5] P. Kapranos, “Current state of semi-solid M.M. Rashidi, “Effects of different
net-shape die casting.,” Metals. vol. 9, no. processing parameters on the semisolid
12, pp. 1–13, 2019. microstructure of Al6061 produced by a
direct thermal method.,” IOP
202
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
204
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Syed Waqar Hasan)
ABSTRACT
Industries provide various fluid products (e.g. petroleum fuels, chemicals agents and mineral
water) to the market which have deep impact on the economy of the country. Normally, bulk
amount of fluids are loaded on customer containers from the refinery storage tanks through pumps.
A flow meter is used to measure the quantity of the fluid filled in the container against which the
customer shall make the payment. Thus, the efficiency of the flow meter is of paramount
importance. It has been observed that air bubbles are formed in the fluid flow while it is transferred
through a pump. Those air bubbles significantly reduce the flow meter efficiency causing serious
financial burden on the customers. In this research work, a novel technique is presented to ensure
no air bubbles in the working fluid enhancing the flow meter efficiency. It was observed that
centrifugal pump designed for providing flow rate of 10 lit/min was only providing 4.9 lit/min
because of very high air bubble density. We applied a line coating of abrasive paper in the tube
transporting the fluid. It was realized that the abrasive paper completely eliminated the bubble
formation within the tube increasing the flow rate to 9.8 lit/min which is 98% closer to the design
value of the pump. This technique can help the industries/customers to measure the transported
fluid correctly and improve industrial economy.
Keywords: flow meter efficiency; turbulent flow; bubble formation; abrasive papr
205
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
206
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ρVD VD 𝑙𝑒 1
Re = = . . . (1) = 4.4Re6 . . . . (2𝑏)
μ 𝑣 d
207
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
from outside. Therefore, we could not cover lit/min, a gradual decrease in bubble density
a larger circumference of the tube with was observed as exhibited in Figure 2b,
abrasive paper. Hence for the diameter of 0.5 however, bubbles are still very prominent and
in (circumference=1.57 in), we placed consistent. Once the lowest flow rate was
abrasive paper of width 0.3 in. examined; complete elimination of the
bubbles was observed (Figure 2c). Although
Description Length Dia Width the Reynolds Number is still considerably
(in) (in) (in) high (i.e. 3678) which is in the transitional
Pipe (Overall) 14 0.5 - limits to turbulent flow but no bubbles were
visualized.
S.no Flow Reynolds Observations
rate number
Venturi Nozzle 6.3 0.5 -
(lit/min)
1 5 8582 High bubble
Abrasive paper 5 - 0.3 density
2 2.6 4466 Low bubble
Transparent tube 5 0.5 - density
3 2.14 3678 No bubbles
Table 1. Details of experimental setup Table 3. Details of Bubble formation in
tube
DC Motor 10 lit/min 2.3 amp
DC Regulator - 10 amp
3. RESULTS
In this research, we designed an experimental
setup to investigate the fluid characteristics
on different flow rates by changing its flow
rate Reynolds number. Three different flow
rates were selected for the experimentation:
(a) highest flow rate (10 lit/min), having
turbulent flow with Reynolds number of
8586, (b) moderate flow rate (2.6 lit/min) Figure 3. Bubble Density in tube against
with Reynolds number of 4466 and (c) lowest Reynolds number
flow rate (2.14 lit/min) with Reynolds 3.2 Fluid behavior in tube provided
number of 3678. abrasive surface
It is mentioned in the experimental section
3.1 Fluid behavior in bare tube that the primary interest of this research is to
It was observed that the highest flow rate investigate the changes in the fluid flow
resulted in turbulent flow producing characteristics due to the increased
tremendous amount of air bubbles in the flow friction/roughness incorporated through an
as shown in Figure 2a. In fact the highest abrasive coating. An abrasive coating was
bubble density was visualized in this flow introduced into the transparent tube (Figure
regime. Once the flow rate was reduced to 2.6 1 b& c), creating a region where the fluid had
209
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
to flow over the abrasive surface. In order to underscore the influence of abrasive paper
compare the effect of abrasive coating, the coverage on the dynamics of bubble
fluid was flown over the abrasive coating at formation.
the three flow rates previously studied. We
previously mentioned extremely high
generation of bubbles at highest flow rate for (a)
bare tube (Figure 2a). A dramatic
elimination of bubbles was noted once the
abrasive paper was placed in the tube as
shown in Figure 4a. It has to be noted that
the fluid velocity, flow rate and Reynolds Full length of abrasive lining
number were unchanged/highest,
nevertheless, there are no bubbles observed. (b)
Even for the lower and moderate flow rates;
the bubbles were not observed in the presence
of abrasive paper.
2/3rd length of abrasive lining
In order to further extend our understanding
about the role of abrasive paper on fluid flow
[8], we reduced the length of the abrasive (c)
paper. Previously, the whole length of the
plastic tube (i.e. 5 in) was provided the
abrasive paper. Now we reduced half of the
length of abrasive paper (i.e. 2.5 in). In other
words, half of the pipe length is bare (i.e.
1/2 length of abrasive lining
uncovered) and the remaining half of the pipe
length is covered by abrasive paper as shown Figure 4. Bubble visualization in tubes
in Figure 4c. We kept the velocity, flow rate embedded with abrasive lining. (a) Full
and Reynolds Number to be maximum. The length of tube provided abrasive paper.
results are shown in Figure 4c. It can be An abrasive paper can be seen placed
noted that reduced length of the abrasive along the length of the tube. It can be
paper allowed some bubble formation in the noted that there are literally no air
tube, although the bubbles observed in this bubbles at maximum flow rate (10
case are notably less than the bare tube case lit/min). (b) 2/3rd length of tube provided
(Figure 2a) but higher than Figure 4a (when abrasive paper. (c) 1/2 length of tube
complete tube length was covered by the provided abrasive paper. Abrasive paper
abrasive paper). Therefore, applying abrasive can be seen placed along the half length
coating to only half of the tube in our of the tube. (a,b) It can be noted that
experiments led to a continued generation of sizes of bubbles are getting smaller at
bubbles, with a noticeable reduction in size maximum flow rate (10 lit/min).
(Figure 4c). Subsequently, expanding the
coverage of the abrasive paper to two-thirds
of the tube length reduced the bubble
formation and there was further decrement in
size and a significant reduction in quantity of
bubbles (see Figure 4b). These observations
210
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Bubble Density Vs Reynolds number The placement of abrasive paper in the tube
creates enormous amount of shear stresses
During the investigation of changing due to extremely high surface roughness of
Reynolds numbers on bubble formation, it abrasive paper in the fluid. The schematic
was found that bubble density (the number of illustration is provided in Figure 5b. This
bubbles per unit area) showed a direct high shear stress is expected to be the main
relation. With increasing Reynolds numbers; reason of bubble elimination from the system
the bubble density increases in bare tubes. At (Figure 5b). At the bottom of Figure 5b an
the highest flow rate, the Reynolds number abrasive paper is shown which is shown to
was (8582), indicating a turbulent flow with offer higher shear stresses than the bare tube
highest bubble density (Figure 3). After surfaces (as in Figure 5a).
slowing down flow rate Reynolds number
dropped to 4466, showing turbulent flow. At
this flow rate bubble density decreased. The Max. Wall shear
slowest flow rate had a Reynolds number of stress Air bubbles
(3678), indicating transitional flow and no
bubbles were observed (Figure 3).
Formation of bubbles
Flow
In the venturi effect produced in the Upstream
experimental setup, the bubbles can form due
Shear stress
to decrease in pressure in a constricted area Velocity profile
distribution
of the fluid flow. As the velocity increases
through the narrow portion of the Venturi (a)
tube, the pressure decreases according to
Bernoulli's principle. If the pressure drop is
low enough; it can cause dissolved gases in No
the fluid to come out of solution, forming Bubbles
are
bubbles. These bubbles then become visible Flow
observed
as they move downstream in the flow. The Upstream
schematic illustration is provided in Figure
5a.
Increased Shear Stress Abrasive
After placing abrasive paper (offered by abrasive paper
paper)
If the length of the tube is fully covered by
(b)
the abrasive paper, there were literally no
bubbles observed at highest flow rate. If the
2/3rd tube length was covered by abrasive Figure 5. Schematic illustration (a) Shear
paper, extremely small and few bubbles were stress distribution and velocity profile
visualized. If 1/2nd tube length was provided before introducing abrasive lining. (b)
the abrasive coating, the bubble density and Schematic illustration of shear stress
sizes further grew. distribution and Velocity profile before
after introducing abrasive lining
211
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
212
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
213
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Syed Abdullah Ali1,*, Furqan1, Hasnain1, Abdul Rauf1, Izhan1, Muhammad Hateem
Arif1
1
Department of Mechanical Engineering, NED University of Engineering &
Technology, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Syed Abdullah Ali)
ABSTRACT
The solid-liquid mixing quality is a vital parameter from an industrial point of view. Its dependence
on various parameters like tank and agitator design, stirring speed, tank diameter, solid particle
diameter, solid density, liquid density, liquid dynamic viscosity and mean solid volume fraction
dictate the necessity of careful consideration of the optimum parameters to ensure the best possible
mixing quality. To understand the effect of the parameters on the mixing quality and to predict the
mixing quality for a given set of parameters of influence, this study aimed to train a machine
learning algorithm and employ data visualization techniques to identify the correlations between
the parameters and the mixing quality. A CFD simulations dataset of a tank-agitator mixing
problem validated against experimental data was used to execute this study. Various machine
learning algorithms including Keras DNN, Linear Regression, Lasso, Elastic Net, SVR, Decision
Tree Regressor, Random Forest Regressor, K-Nearest Neighbors Regressor and MLP Regressor
were utilized to train different models to formulate the and identify the best one. The dataset was
subjected to preprocessing steps and was split into 80:20 for training and testing the model,
respectively. Considering the varying nature of values of the parameters, the dataset was
normalized too to ensure a resulting model that would be free of any bias. The model that
performed the best in terms of predicting the mixing quality were Keras DNN and Random Forest
Regressor. The utilized dataset was considerably small. Therefore, to further extend this work, it
is recommended to amass a significantly larger dataset for a given problem to allow a model to be
trained properly which would ultimately open doors to reaching a higher accuracy.
Keywords: Linear Regression, Data Visualization, Machine Learning
Experimentation or simulations are often
1. INTRODUCTION
required to determine the stirring
The solid-liquid mixing quality in a stirred
configuration that achieves high suspension
tank as a function of the main system
quality. However, conducting real-
parameters (tank size and design, stirrer
size experiments, especially for large tanks,
design, suspension characteristics) is an
can be impractical due to various constraints.
important matter for the industry [1]. The
Small scale experimentations can help but
quality of mixing solid and liquid materials in
scale up laws are not always clearly
a stirred tank depends on various factors such
established.[2-4]
as tank size, design, stirrer design, and
Simulation is free from experimental
suspension characteristics.
constraints, but it is still a care- demanding
and burdensome tool. This is where machine
214
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
learning (ML) algorithms. come into play. diameter, solid density, liquid density, liquid
ML algorithms can analyse data from various dynamic velocity, and mean solid volume
parameters affecting mixing quality and fraction. To tackle this, we engaged in data
predict the optimal mixing conditions. They visualization, data cleaning, and training the
can learn from past data and simulate data for predictive analysis. This section
different scenarios to find the most effective encapsulates our applied methodology and
stirring configurations. findings.
The impact of AI is far-reaching,
2.1 Data Cleaning
transforming traditional approaches and The initial phase of this study involved
unlocking new possibilities for rigorous data cleaning procedures to ensure
manufacturers, maintenance teams, and the integrity and suitability of the dataset for
sustainability efforts, from the initial design subsequent analysis. This included
phase to the recycle/reuse/retrofit phase.[5]
identifying and addressing missing or null
Recent research trends and popular AI values, encoding categorical variables, and
techniques in CFD are highlighted by implementing feature scaling techniques.
Krzysztof Rojek et al. [6]
Therefore, the importance of ML in 2.2 Handling Missing Null Values
predicting mixing quality lies in its ability to Missing or null values within the dataset were
handle complex data relationships and managed using robust techniques such as
provide insights into optimal mixing forward fill, backward fill, or mean fill. This
step aimed to enhance the dataset's
conditions without the need for extensive
completeness while minimizing biases
experimentation. By utilizing ML, industries
introduced by missing data.
can optimize their processes, reduce costs,
and improve product quality. This study 2.3 Categorical Variable Encoding
involves working with a dataset comprising Categorical variables, such as design types,
of Computational Fluid Dynamics (CFD)- were appropriately encoded to enable
meaningful analysis. This transformation
generated and experimentally verified
ensured that categorical information was
results for mixing tanks. The effectiveness represented numerically while retaining the
of mixing in these tanks hinges on several inherent characteristics of the data.
variables, including tank and agitator design,
stirring velocity, tank diameter, solid 2.4 Feature Scaling
To address variations in data ranges across
particle diameter, solid density, liquid
different parameters, feature scaling was
density, liquid dynamic velocity, and mean
conducted. Standardization techniques were
solid volume fraction. To tackle this, data applied to scale all variables uniformly,
visualization, data cleaning, and training the reducing the impact of variables with larger
data for predictive analysis are done. scales on subsequent analyses.
2. METHODOLOGY 2.5 Training and Testing
This study involves working with a dataset A rigorous train-test split methodology,
comprising Computational Fluid Dynamics utilizing an 80-20 ratio, was employed to
(CFD)-generated and experimentally verified train and evaluate various regression models.
results for mixing tanks. The effectiveness of Implementation of regression models was
mixing in these tanks hinged on several carried out using the scikit-learn (sk.learn)
variables, including tank and agitator design, library and tensor flow library in Python,
stirring velocity, tank diameter, solid particle
215
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Solid
Stirring Particle
Speed Diamete
(RPM) r (µm)
Solid
Tank
Density
Diamete
(Kg/m3)
r (m)
216
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Keras
Deep
Liquid
Neural
Dynami
Network
c
(DNN)
Viscosit
y (µ)
Linear
Mean Regression
Solid
Volume
Fraction
(φ)
Lasso
Mixing
quality
Elastic Net
MLP
Support Regression
Vector (Multi-
Regression layer
(SVR) Perceptron
)
Random Assessme
nt and
Forest
Testing Bar Charts
Regressor Reliabilit
y Factors
K K-fold
(10
Neighbour Folds)
s
Regressor
Accuracy
218
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
REFERENCES
[1] Projets, C.M.O., Solid-Liquid Mixing in
Stirred Vessels: Numerical Simulation,
Experimental Validation and
Standard
Suspension Quality Prediction Using
Deviation
Multivariate Regression and Machine
Learning. 2018.
[2] Geisler, R.K., Buurman, C., Mersmann,
A, Scale-up of the necessary power
input in stirred vessels with suspension.
The Chemical Engineering Journal,
(1993). vol. 51.
Mean
[3] Montante, G., Pinelli, D., Magelli, F,
Squared
Scale-Up Criteria for the Solids
Error
Distribution in Slurry Reactors Stirred
with Multiple Impellers. Chemical
Engineering Science, (2003). vol. 58,.
[4] Dickey, D.S., Ten Things You May Not
The above bar charts depict that the most Know About Liquid Mixing Scaleup.
optimum models that have lesser error, Chemical Engineering (2013).
higher accuracy, lower standard deviation, [5] Mahboob Elahi, S.O.A., Jose Luis
and a higher K-fold (10 folds) score are Keras Martinez Lastra & Jose Antonio Perez
DNN and Random Forest Regressor. Garcia A comprehensive literature
4. CONCLUSION review of the applications of AI
In conclusion, this study provided a techniques through the lifecycle of
systematic approach to train a machine industrial equipment. Discover
learning model to predict the mixing quality Artificial Intelligence, (2023). Volume
of a two-phase solid-liquid mixing tank. The 3, article number 43.
utilized dataset was verified using [6] Wyrzykowski, K.R.R., Performance
Computational Fluid Dynamics (CFD) and Scalability Analysis of AI-
calculations and experimental data Accelerated CFD Simulations Across
validation. By adhering to stringent data Various Computing Platforms. Euro-Par
cleaning, model training, and evaluation 2022: Parallel Processing Workshops,
procedures, this research aimed to derive (2023). Volume 13835.
comprehensive insights and contribute to
advancements in fluid dynamics research by
uplifting the potential of machine learning’s
application to train a model that can predict
the mixing quality. Further research can be
done by utilizing machine learning to train a
predictive model using a larger and less
skewed dataset which would yield a higher
accuracy.
219
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Muhammad Bilal)
ABSTRACT
To address the concerns in the manual welding process in terms of efficiency, quality, and
adaptability there is a dare need for Industry 4.0 based smart welding system. For the robot to carry
out the welding task, the robot moves the welding torch along the path of the weld. The position
control along this path of the welding is the most significant challenge. This paper presents inverse
kinematics-based position control for the target pose of the welding robot using Robot Operating
System. Digital model of the robot in the form of URDF was designed for the Robot Operating
System (ROS). Inverse kinematics was used to calculate the joint angles to achieve the desired
position of the end effector using KDL solver. Position planning simulation was performed in Rviz
(ROS Visualization) tool. The calculated and achieved coordinates for a specified position were
compared and the errors were determined. The average error for the 17 poses is determined as
0.6847 mm. The result of this paper highlights a simple and cost-effective solution for position
control by deploying ROS for low cost and effective development of welding digitization.
Keywords: Welding digitization, Position Control, Inverse Kinematics, ROS, Welding Robot
222
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
223
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
Figure 12. Targeted and achieved position values of the passes across x, y & z axis
224
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
The figure 10 below shows the error values at (UET) Taxila and Department of Industrial
different poses. There is quite fluctuation in Engineering, UET Taxila.
the error value for different poses. Each pose
corresponds to the x, y and z axis coordinates REFERENCES
in the space. [1] M. Hunstig, S. Holtkämper, L. Helmich,
and M. Brökelmann, “Smart Ultrasonic
Welding – A Versatile Interconnection
Technology for Power Electronics
Packaging,” Additional Conferences
(Device Packaging, HiTEC, HiTEN, and
CICMT), vol. 2021, no. HiTEC, 2021,
doi: 10.4071/2380-
4491.2021.hitec.000112.
Figure 13. position Error estimation of [2] C. Zhu, H. Yuan, and G. Ma, “An active
targeted and achieved position. visual monitoring method for GMAW
weld surface defects based on random
6. CONCLUSIONS
forest model,” Mater Res Express, vol.
Test results show that ROS based Inverse
Kinematics control system can precisely 9, no. 3, 2022, doi: 10.1088/2053-
achieve the position control of the 1591/ac5a38.
manipulator for the reasonable task [3] S. Amano, Y. Tsujimura, T. Ogawa, and
requirements, which fully meets the actual T. Shibata, “Development of in-process
work requirements. Once the path is planned,
welding torch position control system
the robot can execute the welding task. This
using AI technology,” Welding in the
involves moving the robot arm along the
planned trajectory. During execution, the World, vol. 67, no. 5, 2023, doi:
ROS system can monitor the position and 10.1007/s40194-023-01486-7.
orientation of the robot's end-effector in real- [4] W. Cai, J. Z. Wang, P. Jiang, L. C. Cao,
time. The average error in the achieved
G. Y. Mi, and Q. Zhou, “Application of
position and the targeted position is 0.6847
sensing techniques and artificial
mm. The targeted position and orientation for
each welding point can be compared to the intelligence-based methods to laser
achieved position and orientation. Error in welding real-time monitoring: A critical
the achieved position can be reduced by review of recent literature,” Journal of
tunning the P, I, and D values of the PID Manufacturing Systems, vol. 57. 2020.
controller. doi: 10.1016/j.jmsy.2020.07.021.
[5] P. Pitchipoo, A. Muthiah, K. Jeyakumar,
ACKNOWLEDGEMENTS
Author is grateful for the support of HEC- and A. Manikandan, “Friction stir
NRPU (National Research Program for welding parameter optimization using
Universities) for funding this project. The novel multi objective dragonfly
author further acknowledges the algorithm,” International Journal of
administrative and technical support of the Lightweight Materials and Manufacture,
University of Engineering and Technology vol. 4, no. 4, 2021, doi:
10.1016/j.ijlmm.2021.06.006.
225
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
*
Corresponding author E-mail address: [email protected] (Murtuza Mehdi)
ABSTRACT
Microactuators, sensors, and switches are utilized in applications such as lap-on-chip devices,
deformable mirrors for optical beam deflection, camera auto-focus systems, micro-positioning
platforms, ultrasonic emitters, etc. These devices produce miniaturized motion by converting
electrical, thermal, or electromagnetic energy into kinetic energy of moving parts. In this article,
we have reported the synthesis and characterization of a low-cost thermal actuator based on silver
nano-platelets-based composite and single-component room temperature vulcanizing (RTV)
silicone. The actuator is fabricated by depositing a functional layer of silver amalgamated with a
conductive polymer on a smooth RTV substrate via a non-vacuum-based rod coating technique.
The fabricated device exhibits good thermal actuation behavior in response to the applied thermal
gradient. The device also reveals a stable resistive behavior for the applied temperature difference.
These results suggest that the as-fabricated device can be employed as a micro-actuator, switch, or
valve as well as a thermal sensor for future low-cost MEMS devices.
device also tested for change in its electrical pastecontaining silver nano-platelets was
resistance in response to an applied thermal purchased from
gradient. These characterizations are
performed on custom test setups. The results
reveal that the as-fabricated device not only
shows a good thermal actuation behavior but
also exhibits a decent electrical resistive
behavior for the applied temperature
gradient. These results help explore the
possible application of the fabricated device
in future low-cost MEMS-based
technologies.
2. EXPERIMENTAL SECTION
2.1. Materials
Room temperature vulcanizing (RTV)
silicone was purchased from Bossil
Technology (Malaysia) in the form of a
single-component paste. The material comes
with a commercial name of Grey RTV and is
used as a substrate material. The mechanical
properties and thermal behavior of this
material are discussed in our previous
research work [ref]. PEDOT: PSS paste was
purchased from Sigma Aldrich (Belgium). N,
N-dimethylformamide (N, N-DMF) was Figure 1. Scanning electron microscope
purchased from Merk Millipore with an energy dispersive X-ray
(Germany).A silver conductive spectroscopy (SEM/EDS) of the
prepared RTV/Ag/PEDOT: PSS device.
a) SEM image of Ag/PEDOT: PSS film
showing Ag nanoplatelets dispersed in
PEDOT: PSS, b) EDS results of the film
confirm the presence of Ag in PEDOT:
PSS.
Mechanic, China under the commercial name
MCN-DJ002). The presence of silver nano-
platelets is also confirmed by scanning
electron microscopy (SEM) and energy
dispersive X-ray spectroscopy (EDS) results
as shown in Figure 1. The average size of
these nano-platelets was found to be 5.53 µm.
228
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
ohm meter was connected to the sample to inset shows the photograph of a
record the electrical resistance of the sample thermally actuated device. The error bar
during testing. A K-type thermocouple was shows the standard deviations for the set of
also attached to the hot plate via heat- experiments.
resistant tape to record the temperature Figure 3 shows the thermal actuation
measurements. The temperature was slowly response of the fabricated device. The entire
raised to 250 °C and a change in resistance of response can be divided into three regions,
the sample was recorded. Then hotplate was namely low actuation, rapid actuation, and
turned off and the assembly was cooled down high actuation region. In the low actuation
to room temperature. During this time, the region, the device shows a small deflection
temperature and the resistance measurements height for a small temperature change. The
were continued to record. maximum actuation height of 2 mm is
In the second test, the same setup was used achieved when the device is heated from
with minor changes. Now, the sample was room temperature to a temperature of 38.25
fixed to a 6 mm thick glass plate at one end °C. Afterward, the device shows a fast
by heat-resistant tape. Then the glass plate movement, which embarks the second region
was placed on the hotplate with the of rapid actuation. Here, the device shows a
thermocouple attached. A steel ruler having a fast actuation of 10 mm for a small
least count of 0.05 mm was vertically fixed to temperature change of around 5 °C.
the glass plate for measuring the actuation Afterwards, the third region of high actuation
measurements. The test was initiated by starts from 12 mm at 42 °C to 21.5 mm at
slowly increasing the temperature of the 103.5 °C. In this region, the device shows a
hotplate up to 250 °C and the actuation linear actuation response over a relatively
measurements were recorded. The recorded wide temperature range as compared to the
data was then converted into temperature– first region. The maximum actuation of 9.5
resistance, and temperature–height graphs of mm is observed for a temperature difference
the prepared devices. of 61.5 °C.
3. RESULTS AND DISCUSSIONS 3.2. Electrical Resistive Behavior
3.1. Thermal Actuation Response Figure 4 shows the resistive response of the
fabricated device to the applied temperature
change. It is observable that the normalized
resistance of the device remains nearly
unchanged up to a temperature of 150 °C.
Afterward, the resistance gradually starts
incrementing and reaches a maximum linear
normalized value of 2.05 at 227.5 °C. Then a
drastic increment is observed and the
normalized resistance reaches a maximum
value of 2.8 times the initial one at a
temperature of 220 °C. At this point, the
temperature change is found to produce no
resistance change even when the device is
Figure 3. Actuation behavior of
cooled down to room temperature. This
RTV/Ag/PEDOT: PSS device in
behavior is attributed to two reasons; first, the
response to applied temperature. The
higher temperatures caused the boiling of N,
230
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
N–DMF which was used as a conductive for thermal actuation and electrical resistive
binder for the functional coating. Hence, response for applied thermal gradient. The
increasing the electrical resistance of the device shows good actuation response over a
coated film. Second, these higher wide range of temperatures. The device also
temperatures also caused morphological exhibits electrical resistive behavior for
changes in the film such as bending and applied temperatures up to 185 °C. In
delamination from the RTV at various points. summary, the device can be used as a thermal
These changes also contribute to yielding sensor, thermal actuator, valve, or switch in a
permanent changes in the electrical resistance variety of MEMS-based applications.
of the device. This implies that to use this
ACKNOWLEDGMENTS
device as a thermal sensor, it needs to be
operated below 185 °C to prevent a lasting The authors would like to express their
gratitude to the Higher Education
Commission, Pakistan, and NED University
of Engineering and Technology for
financially supporting this research work by
providing National Research Program for
Universities (NRPU) grant # NRPU-I/17057
and a Ph.D. research fund for the university.
The authors would like to acknowledge the
role of Prof. Dr. Saud Hashmi, Chairman and
Professor, Department of Polymer and
Petrochemical Engineering, NED University
of Engineering and Technology, and Dr.
Rafiq Ahmed, Assistant Professor,
change in electrical resistance.
Department of Polymer and Petrochemical
Figure 4. Resistive behavior of Engineering, NED University of Engineering
RTV/Ag/PEDOT: PSS device in and Technology in providing the facilities for
response to applied temperature. The experimentation.
inset shows the photograph of the device
REFERENCES
under test (DUT). The error bar shows
the standard deviations for the set of [1] V. Raman, Y. H. Cho, H. M. Kim, Y. J.
experiments. Kim, H. M. Sim, and H. K. Kim, “Ag
mesh network framework based nano
4. CONCLUSIONS
composite for transparent conductive
This paper presents the fabrication and functional electrodes for capacitive
characterizations of a low-cost thermal touch sensor and thin film heater,”
actuator produced by coating a functional Ceram Int, vol. 47, no. 19, pp. 27230–
layer of PEDOT: PSS and Ag nanocomposite 27240, Oct. 2021, doi:
functional layer on a smooth RTV substrate 10.1016/j.ceramint.2021.06.145.
via rod coating technique. The fabricated
[2] T. T. A. Nguyen, B. S. Soram, D. T.
device is studied under a scanning electron
Tran, N. H. Kim, and J. H. Lee,
microscope for morphology and cross-
“Enhanced electrochromic capacity
section information. The EDS results also
performances of hierarchical MnO2-
confirm the presence of Ag in the PEDOT:
polyaniline/PEDOT:PSS/Ag@Ni
PSS functional matrix. The device was tested
nanowires cathode for flexible and
231
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
232
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00
Proceedings of IMEC-2024, 6th-7th March 2024, Karachi, Pakistan
233
(Copyrights © Reserved with IMEC-2024)
Published Date: 08/07/2024 Version: 00