International Conference

on
Mathematical Modelling and Scientific Computing

(ICMMSC 2018)

July 19–21, 2018

Organized by:

Indian Institute of Technology Indore

Supported by:
Science and Engineering Research Board (SERB), New Delhi
National Board for Higher Mathematics (NBHM), Mumbai
Council of Scientific and Industrial Research (CSIR), New Delhi

Sponsored by:
Godrej & Boyce Mfg. Co. Ltd., Mumbai Br Plant-4A, Mumbai
Cambridge University Press India, New Delhi
State Bank of India, Simrol Branch, Indore
Unicomp InfoSolutions Pvt. Ltd., Indore
 International Conference on
Mathematical Modelling and Scientific Computing
(ICMMSC 2018)

Director’s Message

IIT Indore is one of the youngest IITs, which began

its undergraduate programme in July 2009, PhD pro-
gramme in July 2010 and various masters programmes
in July 2013. Today, we have over 1500 students of
which, about 100 have been awarded Ph.D. degree and
more than 400 are pursuing PhD research. IIT In-
dore has placed a heavy emphasis on its research pro-
gramme and is recognized by the label of a research-
cum-teaching institute. In this context, I am glad to
state that IIT Indore is ranked 14th overall amongst all
engineering universities and institutions in the coun-
try, and a very impressive 5th in the teaching and resource category in the NIRF
rankings announced for 2018. I am delighted that the International Conference on
Mathematical Modelling and Scientific Computing will be organized by and held at
IIT Indore during July 19 − 21, 2018. The ICMMSC 2018 aims to bring together
experts from academic, scientific and industrial communities to address new chal-
lenges, present their latest research findings, ideas, development and perspectives
of the future directions in the broad areas of applied mathematics, mathematical
modelling and scientific computing. The conference has attracted more than 150
researchers from IITs, IIITs, NITs, Indian and USA universities, research organi-
zations and industries. I extend my warm welcome to the delegates attending the
conference and wish that the discussions and interactions during the conference will
help them to enhance their knowledge and provide directions for their future re-
search collaborations. I extend my warm wishes and hope that all participants have
a very enjoyable stay at IIT Indore.

I congratulate the organizing committee for conducting the ICMMSC 2018 at

IIT Indore and wish them success.
Director, IIT Indore
Professor Pradeep Mathur
 International Conference on
Mathematical Modelling and Scientific Computing
(ICMMSC 2018)

Foreward
It is a great pleasure to warmly welcome you all to the International Conference on Mathe-
matical Modelling and Scientific Computing (ICMMSC 2018). We are delighted to have you all
here at IIT Indore to share your experience and to discuss the recent progress in the theme of the
conference. It will provide early career researchers with a valuable forum to showcase and obtain
feedback on their work from leading researchers in the field and gain insights into the open chal-
lenges. They will be immensely benefited by having discussions and collaborations with globally
renowned peers. The technical programme of the conference broadly reflects the common research
interests of mathematical modelling, differential equations, applied linear algebra, numerical analy-
sis, computational fluid dynamics, solid mechanics, wavelets, wave Propagation, inverse problems,
scientific computing & simulatio, flow heat and mass transfer in porous media, nano-fluids for
heat transfer, etc. An abstract view of all the papers selected for oral and poster presentation are
presented in this program booklet. The conference enjoys the special support of the Director of
the institute without which it might not have been possible to arrange the conference. We take
the opportunity to thank Prof. Abhinav Kranti, Dean, R & D for his continuous support in orga-
nizing the event. We thank Dr. Sk Safique Ahmad, Head, Discipline of Mathematics, IIT Indore
for his support. We are also thankful to Prof. Neelesh Kumar Jain, Dean of Academic Affairs,
Prof. Subhendu Rakshit, Dean of Faculty Affairs, Dr. Kapil Ahuja, Dean of International Affairs,
Dr. Rajesh Kumar, Dean of Planning and Dr. Swadesh Kumar Sahoo, Dean of Administration,
IIT Indore for their support. We must appreciate the untiring efforts made by the organizers and
secretarial assistants in bringing up this programme booklet in time. Organizing an International
Conference requires contribution from all quarters. We are indebted to the members of the Scien-
tific Advisory Board Members for their prompt attention, support and guidance in this endeavor.
Timely cooperations of Office of Registrar, various other administration sections and purchase
section are appreciated for the practical arrangements of the ICMMSC 2018.
Our special thanks are due to Science and Engineering Research Board (SERB), New Delhi,
National Board for Higher Mathematics (NBHM) of Department of Atomic Energy (DAE), Mum-
bai, Council of Scientific and Industrial Research (CSIR), New Delhi, Godrej & Boyce Mfg. Co.
Ltd., Mumbai Br Plant-4A, Mumbai, Cambridge University Press India, New Delhi, and Unicomp
InfoSolutions Pvt. Ltd., Indore for their financial support to organize the symposium. We also
planning to published all accepted papers for the possible publication in Springer-Nature book
proceedings after the conference and through standard review process. We are glad to have a
dedicated organizing committee who are always there to ensure that the particulars of the events
are executed with precision and quality. We once again extend our gratitude to all the delegates
who have spared their precious time to contribute to the scientific merits of the Conference. We
believe that your stay and interactions with us at IIT Indore will strengthen our confidence. Have
fun and pleasure in the beautiful campus of the Institute, the food, and in your visit to Indore,
the cleanest city of India. Thank you for your interest and active participation in ICMMSC 2018.

Dr. Santanu Manna

Dr. Sk. Safique Ahmad
Dr. Shanmugam Dhinakaran
Dr. Niraj Kumar Shukla
Local Coordinators, ICMMSC 2018
 International Conference on
Mathematical Modelling and Scientific Computing
(ICMMSC 2018)

International Advisory Committee

Prof. Ravi P. Agarwal (Texas A&M University-Kingsville, U.S.A.)
Prof. S. Chakraborti (University of Alabama, U.S.A)
Prof. Biswa N. Datta (Northern Illinois University, U.S.A.)
Prof. Herbert Huppert (University of Cambridge, U.K.)
Prof. Anuar Ishak (Universiti Kebangsaan Malaysia, Malaysia)
Prof. Julius Kaplunov (Keele University, U.K.)
Prof. Volker Mehrmann (Technische Universitt Berlin, Germany)
Prof. Atma R. Sahu (University System of Maryland, U.S.A.)
Prof. Kumar Singh (Miami University, U.S.A.)
Prof. Ivan Slapnicar (University of Split, Croatia)
Prof. Shan Zhao (University of Alabama, U.S.A)

National Advisory Committee

Prof. Rafikul Alam (IIT Guwahati )
Prof. T. Amarnath (University of Hyderabad)
Prof. Shreemayee Bora (IIT Guwahati )
Prof. Rama Bhargawa (IIT Roorkee)
Prof. A. Chattopadhyay (IIT Dhanbad )
Prof. Joydev Chattopadhyay (ISI Kolkata)
Prof. B. V. Rathish Kumar (IIT Kanpur )
Prof. J. C. Misra (IIT KGP/IIEST Shibpur )
Prof. A. K. Nandakumaran (IISc Bangalore)
Prof. Souymendu Raha (IISc Bangalore )
Prof. Mythily Ramaswamy (TIFR Bangalore)
Prof. Govindan Rangarajan, (IISc Bangalore)
Prof. S. Sundar, (IIT Madras)
Prof. R. K. Upadhyay, (IIT Dhanbad)

Organizing Committee
Dr. Santanu Manna, (IIT Indore)
Dr. Sk. Safique Ahmad, (IIT Indore)
Dr. Shanmugam Dhinakaran, (IIT Indore)
Dr. Niraj Kumar Shukla, (IIT Indore)
Contents:
1. Invited Speakers
Sr. No. Speaker Title Page No.

1.1 Biswa Nath Datta Computational and Optimization Methods for … 9

1.2 Volker Mehrmann Energy based Modelling, simulation, and… 10

1.3 Uday Banerjee Generalized Finite Element Method: Successes… 11

1.4 Ranjit K. Upadhyay Modeling Chaotic systems & spatiotemporal… 12

1.5 Santosh Kapuria Mixed-field Multiterm Extended Kantorovich… 13

1.6 Rafikul Alam Fiedler and Generalized Fiedler Pencils… 14

1.7 Govindan Rangarajan A Mathematical Model for Storage and Recall… 15

1.8 S. Gupta and S. Sundar Practicality of using Hierarchical Matrices… 16

1.9 Mythily Ramaswamy Control of some fluid models 16

1.10 S. Raha and S. Dixit Data Assimilation by Directly Projecting… 17

1.11 Rama Bhargav Meshfree Methods: Limitations and Applications 18

1.12 J. C. Misra Mathematical Modelling of Electro-osmotic… 18

1.13 J. Prakash Heat and mass transfer effects on MHD nanofluid… 19

1.14 Debasis Chakraborty Numerical methods in Aeropropulsive characterization.. 20

1.15 Mohammad Sajid Chaos in Mathematical Modelling 21

1.16 Kumar Singh Active Control with the Method of Receptances: Recent… 22

1.17 Mani Mehra Extension of wavelets to topologically complicated… 23

1.18 Muslim Malik Controllability of Fractional Differential Equations… 23

1.19 G. K. Vishwakarma Bio-marker Prediction Using Bayesian State-Space… 24

1.20 P. Sharma and N. Gillis Computing nearest stable matrices 24

2. Contributed Talks:
Sr. No. Paper Id Title Page No.

2.1 ICMMSC -127 Water boundary layer flow over an exponentially… 27

2.2 ICMMSC -82 Mathematical Modelling Of Solar Cell 27
2.3 ICMMSC -42 Reliability assessment of Transformer incorporating … 27
2.4 ICMMSC -229 Mathematical Modeling of Moored Ship Motion… 28
2.5 ICMMSC -162 Stability analysis of a film flow down an incline… 28
2.6 ICMMSC -160 Mathematical Modeling and Aalysis of Single… 29
2.7 ICMMSC -212 Simplified multivariate GARCH (SGARCH) model… 29
2.8 ICMMSC -14 A Green’s function approach to analyze the dispersion… 30
2.9 ICMMSC -180 Sustainable and Smart Manufacturing through… 30
2.10 ICMMSC -159 Development of a dimensionally reduced model … 31
2.11 ICMMSC -84 Unsteady stagnation point nanofluid flow with MHD… 31
2.12 ICMMSC -108 Simulation of Livestock Population in India Using… 32
2.13 ICMMSC -34 Stability Analysis For An Seiqr Epidemic Mode… 32
2.14 ICMMSC -61 Mathematical modelling of UASB reactor: A detailed review 32
2.15 ICMMSC -152 Modelling of Quantum Key Distribution Protocols in… 33
2.16 ICMMSC -16 The role of adaptation in an aquatic predator-prey system… 33
2.17 ICMMSC -93 Analysis of Time-Delayed Non-Linear Equations Using … 33
2.18 ICMMSC -121 Controllability result for a nonlocal nonlinear fractional… 34
2.19 ICMMSC -232 Effect of the heterogeneity, initial stress and viscosity… 34
2.20 ICMMSC -172 Study of Memory Effect in an Inventory Model With… 34
2.21 ICMMSC -75 Switching stability of a delayed e-epidemic model in wireless… 35
2.22 ICMMSC -204 Mathematical Modeling and Analysis of Green Sand Casting… 35
2.23 ICMMSC -153 Modelling of Avalanches and Streamers in Resistive Plate.. 36
2.24 ICMMSC -214 Analysis and Comparative Study of Various Performance… 36
2.25 ICMMSC -178 On the Dynamics of Transverse Domain Walls in … 37
2.26 ICMMSC -177 Numerical study of dragonfly wing section in Gliding Mode 37
2.27 ICMMSC -12 Ball convergence theorems for Chebyshev's method … 38
2.28 ICMMSC -233 Optimal multi-echelon integrated supply chain selecting best … 38
2.29 ICMMSC -190 PSI (?) invariant features for face recognition 39
2.30 ICMMSC -215 Position analysis of five bar mechanism under force 39
2.31 ICMMSC -95 Family of Perfect fluid Balls of Class One in General Relativity 39
2.32 ICMMSC -115 Application of Integral transform in the solution of Fractional… 40
Sr. No. Paper Id Title Page No.

2.33 ICMMSC -21 Investigation of Transient stability in IEEE 9- bus system… 40

2.34 ICMMSC -158 Certain Subclasses of Bi-starlike Functions Associated… 40
2.35 ICMMSC -139 The fruit-basket proof for the theorem on… 41
2.36 ICMMSC -224 Aggregation of LARK Vectors for Facial Image Classification 41
2.37 ICMMSC -101 Personnel Selection using Soft Set Theory 42
2.38 ICMMSC -67 Numerical Solutions of Periodic Burgers’ Equation using… 42
2.39 ICMMSC -81 Solution of Space-Time Fractional Order Partial… 43
2.40 ICMMSC -154 Vibrational Spectra of Ozone (O3) Using Lie Algebraic… 43
2.41 ICMMSC -193 Performance Analysis of a Modified Twisted Savonius… 43
2.42 ICMMSC -201 Statistical Analysis of Cancer Data 44
2.43 ICMMSC -194 Digital Mammogram Analysis using Bendlets 44
2.44 ICMMSC -83 The MHD based numerical analysis of unsteady natural … 45
2.45 ICMMSC -7 CFD simulation of a trickle bed with foaming liquid 45
2.46 ICMMSC -46 Numerical study on the efficiency of magnetophoresis… 46
2.47 ICMMSC -239 Mhd Flow Over A Stretching / Shrinking Sheet With… 46
2.48 ICMMSC -97 Numerical Study on Fluid Dynamics of Collapsible Channels 46
2.49 ICMMSC -238 Analysis of the flow of two immiscible electrically… 47
2.50 ICMMSC -146 Comparative study of boundary conditions in LBM… 47
2.51 ICMMSC -136 Finite Element Analysis of MHD Blood Flow… 48
2.52 ICMMSC -174 Computational Analysis of Buoyant Convection… 48
2.53 ICMMSC -116 A solution to enhance the efficiency of CFD analysis… 49
2.54 ICMMSC -118 Mixed Convective Flow Past an Inclined Square Cylinder… 49
2.55 ICMMSC -208 Case Wise Study Of Heterogeneity On The Traversal… 50
2.56 ICMMSC -24 Exact traveling wave solutions and bifurcation analysis… 51
2.57 ICMMSC -65 Impact of Irregularity on the Propagation of SH wave… 51
2.58 ICMMSC -167 Dynamics of Polarised Shear waves propagating … 52
2.59 ICMMSC -164 Wave interaction with a floating circular porous elastic plate 52
2.60 ICMMSC -109 Investigation of the Diabetic Level in Patients Through… 53
2.61 ICMMSC -192 CT scan based modeling and finite element analysis of… 53
2.62 ICMMSC -22 A Coupled Field Finite element analysis of a Dielectric… 54
2.63 ICMMSC -222 Design and Analysis of Chassis of Hydraulic Test Trolley 54
2.64 ICMMSC -205 Structural analysis of hydraulic aircraft jack using solid… 55
2.65 ICMMSC -209 Kinematic Analysis of Five Bar Mechanism under Tolerance 55
2.66 ICMMSC -98 A nonlinear mathematical study of Tropical Cyclone… 56
2.67 ICMMSC -198 A Finite Difference Approach to Study the Three-Dimensional … 56
2.68 ICMMSC -210 Facial Analysis using Jacobians and Gradient Boosting 57
2.69 ICMMSC -202 Experimental investigation on optimization of concentration… 57
2.70 ICMMSC -76 Unsteady Natural Convection in a Fluid-Saturated Porous… 58
2.71 ICMMSC -234 Haar Wavelet: History and Its Applications 58
2.72 ICMMSC -45 Convective flow in an anisotropic porous wavy channel 58
2.73 ICMMSC -211 Theoretical study of dynamics of the micro-nanomotor 59
2.74 ICMMSC -179 Heat Transfer Analysis of CNT-Nanofluid between Two Rotating… 59
Sr. No. Paper Id Title Page No.

2.75 ICMMSC -199 Numerical investigation of natural convection heat … 60

2.76 ICMMSC -25 Thermal radiation and viscous-Ohmic dissipation… 60
2.77 ICMMSC -156 Analysis of Cooling Effect by Spray Impingement… 61

2.78 ICMMSC -157 Experimental and Numerical Validation of Hybrid… 61

2.79 ICMMSC -47 Analysis and Comparative Study of Various Performance… 62

2.80 ICMMSC -173 CIT Communication Architecture for Improved Security… 62

2.81 ICMMSC -186 Modeling of Viscoelastic Fluid Flow Past a Non-Linearly… 63

2.82 ICMMSC -195 Aggregation of LARK Vectors For Facial Image Classification 63

2.83 ICMMSC -197 Robust and Efficient Face Recognition Using Minimum… 64

2.84 ICMMSC -217 Facial Image Classification using Rotation, Illumination, Scale… 64

2.85 ICMMSC -230 A Copy Detection Technique based on Discrete Cosine… 65

2.86 ICMMSC -124 Improved detection of anomalies via integrated multi-scale… 65

2.87 ICMMSC -36 Numerical Study on the Pennes Bioheat Transfer… 66

2.88 ICMMSC -189 Numerical simulation of melting under the influence 66

2.89 ICMMSC -143 Disciplined Investigations to Fix Number of Neuron in Hidden 67

2.90 ICMMSC -155 Reliable Analysis of Riemann Problem in Magnetogasdynamics 67

2.91 ICMMSC -171 Hand Talk System for Deaf And Dumb Person 67

2.92 ICMMSC -105 Simulation and Prediction of Cereals Rating Using Neura… 68

2.93 ICMMSC -176 Comparative Study of Artificial Neural Network… 68

2.94 ICMMSC -143 Disciplined Investigations to Fix Number of Neuron in Hidden… 69

2.95 ICMMSC -125 Comparison of Regression And Neural Network Model For… 69

2.96 ICMMSC -10 Probabilistic intuitionistic fuzzy set based intuitionistic… 70

2.97 ICMMSC -213 Torque Minimization in Surgical Bone Drilling by Particle… 70

3. Poster:
Sr. No. Paper Id Title Page No.

3.1 ICMMSC -48 Discrete Time Minimum Variance Control of Satellite… 71

3.2 ICMMSC -71 Effect of Environmental Pollutants on Rain due to… 71

3.3 ICMMSC -74 Scrutinization of Discrepancy in Nozzle Throat Diameter… 72

3.4 ICMMSC -123 Image Segmentation of Breast Cancer Histopathalogy… 72

3.5 ICMMSC -135 Study of source influx on interdependent calcium… 73

3.6 ICMMSC -147 Spatiotemporal Dynamics of Immunogenic Tumors 73

3.7 ICMMSC -151 A Survey of Time Series Data Prediction Using Multi Agent… 74

3.8 ICMMSC -166 Stability aspects of the first order Hayes equation with… 74

3.9 ICMMSC -184 Predator-prey model with prey group defense and … 76

3.10 ICMMSC -200 Estimation of Annual Rice Production and Its Methane… 76

3.11 ICMMSC -240 A kinetic Simulation of the expanding plasma in axially… 76

3.12 ICMMSC -216 Design of Novel Efficient circuit switching for Quantum-dot… 77

3.13 ICMMSC -219 Role of subducting slab age and dip on the development… 77

3.14 ICMMSC -225 Simulation of I-V Characteristic of Single Electron … 78

3.15 ICMMSC -27 Reliability and sensitivity estimation of shuffle exchange… 78

3.16 ICMMSC -242 On the convection phenomenon from an impermeable… 78

3.17 ICMMSC -243 Investigating conformational flexibility of N-glycans 79

3.18 ICMMSC -244 Investigating Effectiveness of Currently Used Protease… 79

Invited Speakers
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Computational and Optimization Methods for Quadratic Inverse

Eigenvalue Problems Arising in Mechanical Vibration and Structural
Dynamics
a
Biswa Nath Datta
a
Northern Illinois University DeKalb, Illinois 60115 USA
Email: [email protected]
About the speaker
Biswa Nath Datta is Professor of Mathematical Sciences, Adjunct Pro-
fessor of Electrical and Mechanical Engineering, and Distinguished Research
Professor at Northern Illinois University. He has authored more than 115
interdisciplinary research papers, the books Numerical Methods for Linear
Control Systems Design and Analysis and Numerical Linear Algebra and Ap-
plications, and several associated software packages. He was elected a Fellow
of IEEE in 2000, inducted an Academician of the Academy of Nonlinear Sci-
ences in 2002, and named a recipient of a Senior Fulbright Specialist award,
an IEEE Distinguished Lecturer award, and several IEEE Plaques of Honor.
Abstract. The Quadratic Eigenvalue Problem is to find eigenvalues and eigenvectors a quadratic
matrix pencil of the form P (λ) = M λ2 + Cλ + K, where the matrices M, C, and K are square
matrices. Unfortunately, The problem has not been widely studied because of the intrinsic diffi-
culties with solving the problem in a numerically effective way. Indeed, the state-of-the-art com-
putational techniques are capable of computing only a few extremal eigenvalues and eigenvectors,
especially if the matrices are large and sparse, which is often the case in practical applications.
The inverse quadratic eigenvalue problem, on the other hand, refers to constructing the matrices
M, C, and K given the complete or partial spectrum and the associated eigenvectors. The in-
verse quadratic eigenvalue problem is equally important and arises in a wide variety of engineering
applications, including mechanical vibrations, aerospace engineering, design of space structures,
structural dynamics, etc. Of special practical importance is to construct the coefficient matrices
from the knowledge of only partial spectrum and the associated eigenvectors. The greatest com-
putational challenge is to solve the partial quadratic inverse eigenvalue problem using the small
number of eigenvalues and eigenvectors which are all that are computable using the state-of-the-
art techniques. Furthermore, computational techniques must be able to take advantage of the
exploitable physical properties, such as the symmetry, positive definiteness, sparsity, etc., which
are computational assets for solution of large and sparse problems. This talk will deal with two
special quadratic inverse eigenvalue problems that arise in mechanical vibration and structural
dynamics. The first one, Quadratic Partial Eigenvalue Assignment Problem(QPEVAP), arises
in controlling dangerous vibrations in mechanical structures. Mathematically, the problem is to
find two control feedback matrices such that a small amount of the eigenvalues of the associated
quadratic eigenvalue problem, which are responsible for dangerous vibrations, are reassigned to
suitably chosen ones while keeping the remaining large number of eigenvalues and eigenvectors
unchanged. Additionally, for robust and economic control design, these feedback matrices must be
found in such a way that they have the norms as small as possible and the condition number of the
modified quadratic inverse problem is minimized. These considerations give rise to two nonlinear
unconstrained optimization problems, known respectively, as Robust Quadratic Partial Eigenvalue
Assignment Problem (RQPEVAP) and Minimum Norm Quadratic Partial Eigenvalue Assignment
Problem (MNQPEVAP). The other one, the Finite Element Model Updating Problem (FEMUP),
arising in the design and analysis of structural dynamics, refers to updating an analytical finite
element model so that a set of measured eigenvalues and eigenvectors

9
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

from a real-life structure are reproduced and the physical and structural properties of the original
model are preserved. A properly updated model can be used in confidence for future designs
and constructions. Another major application of FEMUP is the damage detections in structures.
Solution of FEMUP also give rises to several constrained nonlinear optimization problems. I will
give an overview of the recent developments on computational methods for these difficult nonlinear
optimization problems and discuss directions of future research with some open problems for future
research. The talk is interdisciplinary in nature and will be of interests to computational and
applied mathematicians, and control and vibration engineers and optimization experts.

Energy Based Modeling, Simulation, and Control of Coupled Systems

Volker Mehrmanna
a
Institut fur Mathematik, MA 4-5, TU Berlin, Str. des 17. Juni 136, D-10623 Berlin
Email: [email protected]

About the speaker

He has done his doctorate in 1982 with Ludwig Elsner on control theory
in 1987 at the University of Bielefeld. He was an assistant professor from
1990 to 1992 at the RWTH Aachen. In 1993 he accepted the call to the
Technical University of Chemnitz and since 2000 he is a professor at the
Institute of Mathematics of the Technical University of Berlin. From June
2008 to May 2016 he was spokesman for the DFG Research Center Matheon.
From January 2011 to December 2013 he was President of the Society for
Applied Mathematics and Mechanics (GAMM). In 2011, he was granted an
Advanced Grant by the European Research Council (ERC) on modeling,
simulation and control of multi-physics systems. Mehrmann is a member of
the German Academy of Engineering Sciences (Acatech).
Abstract. Motivated from modeling modern energy transport networks, in particular those from
different physical domains, the modeling framework of port-Hamiltonian systems is discussed. The
classical port-Hamiltonian approach is systematically extended to constrained dynamical systems
(partial-differential-algebraic equations). A new algebraically and geometrically defined system
structure is derived, which has many nice mathematical properties. It is shown that this structure
is invariant under Galerkin projections, changes of basis, and that a dissipation inequality holds
and it is automatically stable and passive. Furthermore, the new representation is very robust to
perturbations in the system structure. Using this structure leads to many advantages for tasks like
model reduction or large scale simulation. There exist, however, many open problems associated
with port-Hamiltonian systems, which will also be discussed.

10
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Generalized Finite Element Method: Successes and Challenges

Uday Banerjeea
a
Syracuse University
Email: [email protected]

About the speaker

He is a Professor of Mathematics at Syracuse University and currently
serving as the Chair of the Department of Mathematics. He received his PhD
in Applied Mathematics from University of Maryland at College Park. He
did B.Sc. from Patna University and M.Sc,. from IIT, Kanpur. His research
interests are in numerical solution of partial differential equations. Currently,
he is working of the Generalized Finite Element Method to approximate non-
smooth solutions of elliptic partial differential equations.
Abstract. The Generalized Finite Element Method (GFEM) is used to ap-
proximate non-smooth solutions of PDEs associated with many problems of
engineering interest, e.g., problems involving voids, inclusions, crack-propagation, interface prob-
lems, etc. The GFEM is an extension of the standard Finite Element Method (FEM), where the
standard piecewise polynomial trial space of the FEM is augmented by a suitably chosen, problem
dependent, enrichment space that often contains non-polynomial functions. This method (i) allows
the use of a simple mesh that does not conform to the features of the underlying problem and (ii)
yields optimal order of convergence. However, the GFEM yielding the optimal order of convergence
could be badly conditioned for certain choices of the enrichment functions; it could be difficult to
solve the underlying linear system.
In this talk, we will present sufficient conditions on the enrichment space that will guarantee
that the conditioning of the GFEM is not worse than that of the standard FEM and the conditioning
is robust with respect to the position of the mesh (relative to the features of the underlying
problem). A stably well-conditioned GFEM is called the Stable Generalized Finite Element Method
(SGFEM). We will present a simple local procedure that changes the GFEM into an SGFEM for
a certain class of problems. However, the choice of enrichment functions is non-trivial as this local
procedure may not work for other problems. We will illuminate these ideas through two specific
applications.

11
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Modeling Chaotic Systems & Spatiotemporal Transmission Dynamics of

Recent Ebola Spread and Outbreak
Ranjit Kumar Upadhyaya
a
Department of Applied Mathematics Indian Institute of Technology (Indian School of Mines),
Dhanbad
Email: ranjit [email protected]

International Conference on Mathematical Modelling and Scientific Computing

July 19-21, 2018, Indian Institute of Technology Indore

About the speaker

After completing his Ph.D from Indian Institute of Technology, Delhi
under the supervision of Prof. S.R.K. Iyengar. Prof. Ranjit Kumar Upad-
hyay had joined the Department of Applied Mathematics, Indian Institute
of Technology (ISM), Dhanbad. He has been presently holding the post of
a Professor since 2010 in the Department of Applied Mathematics.
Abstract. The concept of chaos is one of the major discoveries of 20th
century. It (deterministic chaos) has been studied extensively in various
fields, here I will discuss a few important areas. In spite of abundant past
and current work in this direction, there is still no broadly applicable, con-
venient, generally accepted definition of the term chaos. In this talk, I will discuss about different
definition of (opinion about) chaos including my own opinion which can be applied very generally
to situations that are commonly encountered, including attractors, repellers, periodically and non-
periodically forced systems. I will also discuss some related concepts like robust chaos, Ecological
vs Mathematical chaos, crisis, edge of chaos and wave of chaos and examine the observability of
these phenomena in eco-epidemiological systems. I will present the development of a model epi-
demic and eco-epidemic systems (both spatial and temporal) based on sound ecological principle
and discuss why and how the model was developed. What is special about it and what kind of
questions of ecological interest can be answered with the help of this model system? Different
types of paradoxes in ecology and its resolution. Next, I will discuss about my recent work on
spatiotemporal transmission dynamics of recent Ebola spread and outbreak in West Africa. In this
talk, I will try to explain the modeling and simulation for the virus dynamics to better understand
and characterize the transmission trajectories of the Ebola outbreak. I will present the develop-
ment of a model epidemic and eco-epidemic systems (both spatial and temporal) based on sound
ecological/epidemiological principle and discuss why and how the model was developed. What
is special about it and what kind of questions of epidemiological interest can be answered with
the help of this model system? To understand the Ebola transmission dynamics, we formulate
a compartmental epidemic model with exponentially decaying transmission rates and study the
impact of control measures to basic public health using an SEIR and SEIRHD models. We com-
pare the simulated results with the most recent reported data of Ebola infected cases in the three
most affected countries Guinea, Liberia and Sierra Leone. I will discuss about my recent work on
spatiotemporal transmission dynamics of recent Ebola spread and outbreak in West Africa. The
proposed Ebola epidemic model provides an estimate to the potential number of future cases. Two
dimensional simulation experiments show that infectious population and the number of deaths in
Sierra Leone will increase up to one and a half year without control but it will decline after two
years.

12
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Mixed-field Multiterm Extended Kantorovich Method for Accurate

Prediction of Edge Stresses in Piezothermoelastic Laminates
Santosh Kapuriaa,b
a
CSIR-Structural Engineering Research Centre, Taramani, Chennai 600113, India
b
Department of Applied Mechanics, Indian Institute of Technology Delhi, Hauz Khas, New Delhi
110016
Email: [email protected]

About the speaker

Prof. Santosh Kapuria received his B.E. degree in Civil Engineering
from Jadavpur University, Kolkata, in 1989, and M.E. in Structural En-
gineering from Indian Institute of Science, Bangalore, in 1991, securing
first rank in the Department of Civil Engineering in both the institutes.
He obtained his Ph.D. degree in Applied Mechanics from Indian Institute
Technology, Delhi in 1997.

Abstract. Composite and sandwich laminates are widely used in high

performance light weight structures in several applications such as aircraft
wings, marine propeller blades, wind turbine blades, spacecraft and robotic
arms. Such laminated structures also provide the opportunity to accom-
modate surface-bonded or embedded smart piezoelectric transducer layers to introduce self-sensing
and actuation capabilities. However, the material discontinuity at the interfaces of adjacent dis-
similar layers, and geometric discontinuity at the free edges of such laminated structures give rise
to development of large out-of-plane stresses in the vicinity of free edges (free edge effect) under
thermo-electro-mechanical loadings, which are known to be the prime reason for initiation of de-
lamination damage in these laminates. Hence accurate prediction of free edge stress field, which is
truly three-dimensional (3D) in nature, is essential for designing such structures and to ensure its
integrity. Since the 3D elasticity based exact analytical solution for the free edge problem is not
known, there has been a continuous effort to develop various numerical and approximate analytical
solutions for the problem, starting from the finite difference method based solution of Pipes and
Pagano (1970). The extended Kantorovich method (EKM) is a powerful iterative semi-analytical
method for solving partial differential equations (PDEs), which was proposed by Kerr in 1968
based on a solution process given by Kantorovich and Krylov in 1936. The author group has
recently generalized the method by developing the mixed-field multiterm extended Kantorovich
method (MMEKM) to obtain 3D elasticity solutions of laminated plates. In this method, both
displacement and stress variables are considered as field variables, and are represented as n-term
series of the product of separable functions in the inplane and thickness directions. The iterative
solution process alternatively considers the functions in one direction as known (from initial guess
or previous iteration), and transforms the PDEs into set of ordinary differential equations (ODEs)
for the functions in the other direction, using the Reissner-type mixed variational principle. The
ODEs are solved analytically. The iterations over the two directions are continued until the desired
convergence is achieved, which generally happens very rapidly. The mixed-field approach ensures
satisfaction of all boundary conditions at the free edges and interlaminar continuity conditions ex-
actly and pointwise, and also yields the same order of accuracy for both displacements and stresses.
The MMEKM is employed here to obtain accurate 3D piezothermoelasticity solution of the free
edge stress field in infinitely long laminated panels under extension, bending, twisting,

13
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

electric potential and thermal loadings. The imperfect interlaminar bonding is modeled using the
linear spring-layer model, which is incorporated in the variational framework. The numerical results
exhibit rapid convergence with respect to number of terms and iterations and excellent accuracy
for a variety of laminate configurations under all loading conditions. The solution is general, and
is applicable to symmetric, unsymmetric, cross-ply and angle-ply hybrid laminates having surface
bonded or embedded piezoelectric layers. The solution also successfully captures the singular
nature of free edge interlaminar stresses in elastic and hybrid laminates under various loadings.
The results show significant effect of electromechanical and thermoelectric coupling on the free
edge stresses in hybrid laminates. Finally, we show how the free edge stresses due to extension,
bending, twisting and thermal loads can be controlled by applying appropriate actuation potential.

Fiedler and Generalized Fiedler Pencils for Rational Eigenvalue

Problems
Rafikul Alam
Department of Mathematics Indian Institute of Technology, Guwahati - 781 039, INDIA.
[email protected]

About the speaker

Prof. R. Alam received his M.Sc. in Mathematics from Indian Institute
of Technology, Kharagpur in 1991 and Ph.D. in Mathematics from Indian
Institute of Technology, Bombay in 1996. His research interest are Spec-
tral perturbation theory for matrices and operators, Structured perturbation
theory for matrices, Numerics of linear and nonlinear eigenvalue problems,
Numerics of multiparameter eigenvalue problems.
Abstract. Linearization is a standard method for computing eigenval-
ues, eigenvectors, minimal bases and minimal indices of matrix polynomials.
Frobenius companion pencils are examples of linearizations of matrix poly-
nomials. Fiedler and generalized Fiedler companion pencils have been introduced recently which
provide structured and unstructured linearizations of matrix polynomials. Rational eigenvalue
problems arise in many applications and the computation of poles and zeros of rational matrix
functions is an important task. For computing eigenvalues, eigenvectors, poles, minimal bases
and minimal indices of rational matrix functions, we introduce Fiedler and generalized Fiedler
companion pencils for rational matrix functions and show that these pencils are linearizations of
the rational matrix functions in an appropriate sense. We describe the recovery of minimal bases
and minimal indices of rational matrix functions from those of the Fiedler and generalized Fiedler
pencils.

14
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A Mathematical Model for Storage and Recall of Images in the Human

Brain
Govindan Rangarajana
a
Professor, Department of Mathematics and Chairman, Division of Interdisciplinary Research
Indian Institute of Science Bangalore 560 012 India
Email: [email protected]
URL: http://math.iisc.ernet.in/∼rangaraj/

About the speaker

Prof. Govindan Rangarajan obtained his PhD from University of Mary-
land, College Park, USA. He then worked at the Lawrence Berkeley Lab,
University of California, Berkeley before returning to India in 1992. He
is currently a Professor of Mathematics at the Indian Institute of Science,
Bangalore. He is the Chairman of the Division of Interdisciplinary Re-
search. This Division comprises 10 Departments and Centres of IISc. He
is also in-charge of the Office of Development and Alumni Affairs of IISc.
This Office’s mandate includes raising CSR funds from corporates. His
research interests include nonlinear dynamics and chaos, time series anal-
ysis and brain machine interface. He is a J. C. Bose National Fellow. He is also a Fellow of the
Indian Academy of Sciences and the National Academy of Sciences, India. He was knighted by the
Government of France: Knight (Chevalier) of the Order of Palms. He was also a Homi Bhabha
Fellow.
Abstract. We propose a mathematical model for storage and recall of images in the human brain
using coupled maps. We start by theoretically investigating targeted synchronization in coupled
map systems wherein only a desired (partial) subset of the maps is made to synchronize. A simple
method is introduced to specify coupling coefficients such that targeted synchronization is ensured.
The principle of this method is extended to storage/recall of images using coupled Rulkov maps
that are often used to model neuronal networks in the brain. The process of adjusting coupling
coefficients between Rulkov maps to store (display) a desired image mimics the process of adjusting
synaptic strengths between neurons to store memories. Thus our method uses both synchronization
and synaptic weight modification. The stored image can be recalled by providing an initial random
pattern to the dynamical system. The storage and recall of the standard image of Lena is explicitly
demonstrated.

15
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Practicality of using Hierarchical Matrices on a Transient Fluid Solver

Based on Finite Pointset Method (FPM)
Sahil Guptaa and S. Sundara
a
Department of Mathematics, IIT Madras.
Email: [email protected]
About the speaker.
Dr. S. Sundar is currently Professor of Mathematics at IIT Madras
and also the Head of the Department of Mathematics, IIT Madras. His
research interests include PDE Modeling and Numerical Simulation. He is
the DAAD Research Ambassador and Distinguished Alumni of TU Kaiser-
slautern, Germany.
Abstract. An industrial scale computational fluid dynamics simulation
can take a very long time depending on the size and complexity of the
problem. In this work, we use an approach based on H-matrices, which
approximate matrices in a data sparse way and the approximate arithmetic
for H-matrices is almost optimal. One of the most time consuming process
in a FPM solver is solution to the Poisson problem, which needs to be computed at every time step.
Currently, the robust methods for solving Poisson problem are Krylov subspace based iterative
methods. To speed up convergence of the solution, H-arithmetic is used to find a suitable pre-
conditioner for the Krylov space method. Here we analyze the practicality of using H-matrix
techniques to reduce the computation time of a 2D lid driven cavity solver based on FPM.

Control of Some Fluid Models

Mythily Ramaswamya
a
Centre for Applicable Mathematics, Tata Institute of Fundamental Research, Bangalore
Email: [email protected]

About the speaker.

Prof. Mythily Ramaswamy is a Professor in the School of Mathematics
at T.I.F.R, Bangalore and had earlier completed her term as the Dean of
the center. She was awarded PhD from University of Paris 6 in 1990. Her
research interest lies in Analysis of differential equations and applications
to control problems. She has been visiting faculty at many international
institutions. She was Sofia Kowalewsky visiting Professor at University of
Kaiserslautern, Germany, CNRS Poste Rouge in University of Dauphine
and Toulouse. She was awarded Fulbright-Nehru Academic and Profes-
sional Excellence Fellowship, for 2 semesters to visit Virginia Tech., USA.
Recipient of the Karnataka State Kalpana Chawla award for women scien-
tists, 2004. She is a Fellow of the Indian Academy of Sciences, Bangalore, the National Academy
of Sciences, Allahabad and serves on the Editorial Board of the Journal of Ramanujan Mathe-
matical Society, Proceedings of Indian Academy of Sciences-Mathematical Sciences and Boundary
Value Problems. Currently she is also a guest professor at IIT Gandhinagar and serves as a Board
member of IIT Gandhinagar and NIT Calicut Board of Governors and is a member of the IISERs
standing committee and SERB Project Advisory Committee.

16
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Abstract. After introducing controllability (reaching a desired state in finite time) and stabiliz-
ability (reaching a steady state as time tends to infinity), for ODE systems, I will discuss these
issues for some PDE systems, arising in fluid models.

Data Assimilation by Directly Projecting the Nominal Dynamics on the

Data Manifold
Soumyendu Rahaa and Saurabh Dixita
a
Department of Computational and Data Sciences Indian Institute of Science, Bangalore 560012,
India
Email: [email protected]
URL: http://cds.iisc.ac.in/faculty/raha/

About the speaker

Dr. S. Raha is currently working in the Department of Computational
and Data Sciences Indian Institute of Science, Bangalore, India. His re-
search interests include Computational Study of Stochastic Differential-
Algebraic and stiff Stochastic Differential Equations, Math. Software
for Control and Direct Optimization of High Index Differential-Algebraic
Equations, Data Assimilation and Inverse Problems of Stochastic and Con-
strained Dynamics etc.
Abstract. We describe a differential-algebraic equation (DAE) approach
to data assimilation by augmenting the dynamical system model with an
algebraic variable and by treating the observed data as a constraint manifold. The variation of this
additional parameter is analysed with respect to the state variables and the time variable to study
its dynamical and noise characteristics. The augmented dynamics is then used to filter and forecast
the state variables. As an application of this approach, the heat equation for a rod is modelled
on a graph with the algebraic variable. Together with the measurements, the system is posed as
a semi-explicit index-2 DAE. The algebraic variable captures the uncertainty in the measurement
data with respect to temperature and time. This is interpreted as the unmodeled dynamics due
to the inhomogeneities in the thermal diffusivity of the rod. Further statistical analysis is done to
characterize the noisy behaviour of the algebraic variable. The nominal dynamics is then modified
with the uncertainty quantified with the algebraic variable which is now expressed as a random
variable that is used to provide sample paths for the prediction or forecast simulations.

17
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Meshfree Methods: Limitations and Applications

Rama Bhargava
a
Department of Mathematics, I.I.T. Roorkee, Uttarakhand -247667, INDIA.
[email protected]

About the speaker

Dr. Rama Bhargav earned her MSc (Mathematics) University of Roor-
kee, in the year 1975. She completed her Ph.D in Mathematics, University
of Roorkee, in the yaer 1978.Her research interest are Computational Fluid
Dynamics, Finite Element,
Abstract. Mathematical modeling provides a tool for conceiving and re-
solving problems in the real world, ranging from everyday numeracy level
to sophisticated equipments or research. They have great potential for
multidisciplinary work including designers, manufacturers (CAD, CAM)
and academic endeavours. This helps in crossing the boundaries between
a student of mathematics to the world outside the class room. These models are usually either
governed by optimization approach or by a type of differential equations. For solving these com-
plicated systems, one has to deal with the numerical approach, may be grid based or grid Tree
methods. Various grid based methods including Finite difference, Finite element, Boundary ele-
ment etc. sometimes feel themselves limited where the domain is undefined or dealing with phase
change or highly irregular domain. In the past two decades, meshfree methods has emerged into a
new class of computational methods which successfully deal with such problems. The present lec-
ture will highlight an overview of the development, and type and computational approach of these
methods. Given the proper treatment, how these methods can be made more efficient, eliminating
the existing bottleneck of computational expenses in these methods will also be discussed.

Mathematical Modeling of Electro-osmotic Flows in Physiological

Systems
J.C.Misraa
a
Indian Institute of Engineering Science and Technology, Shibpur Howrah-711103, West Bengal
(Formerly, Professor and Head, Department of Mathematics, I.I.T. Kharagpur)
Email: [email protected]

About the speaker.

Professor J. C. Mishra, is currently working in Indian Institute of En-
gineering Science and Technology, Shibpur, Howrah -711103, India. He
was also formerly, Pro Vice-Chancellor, SOA University, Bhubaneswar for-
mer Professor and Head, Department of Mathematics, IIT Kharagpur Ex-
President, Mathematical Sciences Section, Indian Science Congress.
Abstract. Electro-osmotic flows of physiological fluids in micro-vessels
will be discussed during the lecture. After explaining the basic concepts
and the underlying assumptions, an attempt will be made to develop a
mathematical model that would depict the electro-osmotic flow behavior

18
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

of a non-Newtonian fluid in a microchannel. The model will be analyzed in the sequel. Of particular
concern in the investigation will be to examine the effect of heat transfer on electro-osmotic flow
in microvessels under the influence of an external magnetic field. The model and its theoretical
analysis will be based on the consideration of Cattaneo-Christov heat flux model, because of
inherent weaknesses in the Fourier’s law of heat conduction. Considering that many industrial
fluids and most physiological fluids exhibit non-Newtonian behavior, the fluid that flows electro-
osmotically in the microvessel will be considered non-Newtonian. The theoretical analysis will be
carried out by the use of appropriate analytical and numerical methods. The effects of Hartmann
number, surface zeta potential and Joule heating on electroosmotic flow velocity and temperature
will be discussed. The impact of various material parameters on skin friction, local entropy, as
well as on Nusselt number will also be discussed. Computational results will be presented for the
electroosmotic flow of blood in the micro-circulatory system.

Heat and Mass Transfer Effects on MHD nanofluid Flow Over a Semi
Infinite Flat Plate Embedded in a Porous Medium & Two-phase MHD
Flow and Heat Transfer in a Horizontal Channel
J. Prakasha
a
Department of Mathematics, Faculty of Science, University of Botswana, Private Bag 0022,
Gaborone, BOTSWANA.
Email: [email protected]

About the speaker.

Prof. J. Prakash earned his M.Sc.(Mathematics) from Agra Univer-
sity, INDIA in 1965 and M.Sc.(Statistics) from Meerut University, INDIA
in 1971. He has done his DPhil (Applied Mathematics), from Garhwal
University, INDIA, in 1984. His research interest are Bio-Fluid Dynamics,
Fluid Mechanics, Computational Fluid Dynamics (CFD), Mathematical
Modelling, Mathematical Methods.
Abstract. There are two problems. In the first, we have investigated the
effects of radiation absorption, chemical reaction and diffusion thermo on
free convective heat and mass transfer flow of a Nano fluid past a semi-
infinite vertical flat plate embedded in a porous medium. The temperature and concentration at the
surface are assumed to be oscillatory type. We have considered four types of cubic Nano particles
namely, silver, aluminum, copper and titanium oxide which are uniform in size with water as a base
fluid. In the second, we have investigated steady, mixed convective, laminar flow of incompressible,
electrically conducting and heat absorbing two immiscible viscous fluids in a horizontal channel
where the fluid in the region (region I) is saturated with porous medium and the region (region II)
is occupied by a clear viscous fluid. A uniform magnetic field is applied in the transverse direction,
the fluids rise in the channel driven by uniform pressure, the heat transfer is influenced by thermal
radiation. The equations are modeled using the fully developed flow conditions. An exact solution
is obtained for the velocity and temperature distributions. The graphical results are presented, and
the physical aspects are discussed in detail to interpret the effect of various significant parameters
of the problem. The effect of the skin-friction, rate of heat transfer coefficients at the channel walls
is tabulated.

19
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Numerical Methods in Aeropropulsive Characterization of Flight

Vehicles
Debasis Chakrabortya
a
Defence Research and Development Laboratory Hyderabad - 500058, INDIA
Email: debasis [email protected]

About the speaker.

Dr, Debasis Chakraborty earned his M.Sc. degree in Applied Mathe-
matics, from Javadpur University, Kolkata, M. Tech in Computational Me-
chanics from IIT, Kharagpur and Ph.D in Aerospace Engineering from IISc,
Bangalore. His research interest includes CFD, Aerodynamics, Propulsion,
Turbulence and Combustion modeling, Transient and unsteady flows, Nu-
merical methods.

Abstract. Over the last few decades, CFD has developed into a rich and diverse subject and
emerged as a major component of applied and basic fluid dynamic research along with theoretical
and experimental studies. Simultaneous development of new computers, numerical algorithms,
physical and chemical models of flow physics are responsible for the greater impact of CFD in both
basic and applied scientific/engineering problems. Presently, CFD is an integral part of the design
process of air-frame and engines for all major aerospace companies in the world. While experimen-
tal testing will always remain an integral part of the design, CFD is decreasing the dependence on
the more expensive, time consuming experimental testing or rather using experimental work more
effectively and economically. Three dimensional Euler, Reynolds Averaged Navier Stokes (RANS)
and Large Eddy Simulation (LES) solvers are developed in-house using state of art numerical tech-
niques and physical models. Systematic validations were carried out through comparisons against
reliable experimental results to assess their predictive capabilities and range of applications. The
codes are regularly upgraded with new numerical schemes, turbulence models and Magneto-Hydro-
Dynamics model etc. The grid free CFD solver is integrated with 6-DOF trajectory equations to
study the store separation problems. Coupled flow solvers are developed for fluid-structure in-
teraction, hypersonic reentry and high speed aerothermal problems. The indigenously developed
CFD codes are routinely used in DRDL to predict various aerodynamic parameters pertaining
to DRDO missile systems. Aerodynamic characterization of missiles in complete M –δ–φ flight
regime, control surface deflection studies, heat shield separation of hypersonic air breathing mis-
siles, store separation from the aircraft, study of plume impingement on jet deflector, low speed
(incompressible limit) characterization of vehicles etc., are some of the notable applications of the
codes.
Core competence is developed in DRDL to simulate turbulent reacting and nonreacting flows
for missile propulsion system design. Important User Defined Functions (UDFs) are developed to
apply commercial CFD solvers in design. Open source CFD software has been customized to solve
many complex propulsion problems. Systematic validations were carried out through comparisons
against reliable experimental results before applying these indigenous, opensource and commercial
CFD tools in the design exercises. Extensive nonreacting and reacting simulations were carried out
for the development of flightworthy scramjet propulsion system for hypersonic airbreathing cruise
vehicle and CFD simulations guided the experimental testing. Excellent match is obtained between
experiment and pre-test prediction for various performance and flow parameters. Accurate

20
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

estimation of heat flux obtained from high fidelity CFD simulations are used in thermostructural
design of the combustor. Aerodynamic and propulsion parameters obtained from end-to-end sim-
ulation (comprising of nonreacting flow in external surfaces and reacting flow in the combustor)
provide vital input for mission design. CFD based jet vane correlations are adopted in the flight
computer for tactical missiles and forms an integral part of missile control and guidance. The prob-
lem of high temperature in the base cavity caused due to interaction of free stream and exhaust
plume at high altitudes could be analysed only through CFD methods. Performance prediction
of installed air intake of ramjet missiles, Solid Rocket Motor (SRM) port flow field analysis, com-
bustion instability prediction of SRM, Plume-canister / plume-launcher interaction etc. are some
of the other notable applications of CFD methods in propulsion system design. Understanding of
complex flow phenomena has reduced the developmental cost and time of the system significantly.
Development, validation and application of CFD codes for aeropropulsive design of missiles will be
highlighted in the presentation.
Chaos in Mathematical Modelling
Mohammad Sajida
a
College of Engineering, Qassim University, Saudi Arabia
Email: [email protected]

About the speaker.

Dr. Mohammad Sajid was born on September 22, 1974 at Shahja-
hanpur, Uttar Pradesh, India. He is presently an Associate Professor
of Mathematics at Qassim University, Saudi Arabia. He qualified NET
(UGC) in 1995. He received his PhD degree in Mathematics at Indian
Institute of Technology Kanpur, India in 2005. He did also Postdoc under
NBHM fellowship. His current research interests are in the fields Real and
Complex Dynamics associated to fixed points, chaos and Fractals. He has
finished two research projects associated to chaos and fractals in dynami-
cal systems under Deanship of Scientific Research, Qassim University. He
has published about 40 research papers in international/national refereed
journals/conferences. He has delivered 44 talks and presentations in national and international
conferences, workshops and training programs. He gave 10 lectures series on chaos and fractals at
Qassim University. He is referee and editorial board member of several journals of mathematical
sciences and engineering. He has written many reviews for MathSciNet (AMS) and Zentralblatt
Math (zbMath) as a reviewer. Inspite of above, he is actively fully participated about 43 in
national/international conferences or workshops/training including at ICTP, Italy. He has also
chaired various sessions in some conferences. In past and present, he has membership of 15 pro-
fessional societies.
Abstract. The aim of this presentation is to exhibit chaos in mathematical modelling. Various
researches are resolved by using chaos theory in the last few decades. It can easily visualize by
using fractals. It is now applicable in almost all disciplines of science, engineering and technology.
Often, chaotic behavior may get in almost all nonlinear systems. Chaos exists everywhere in the
world since it seems that it is nonlinear world. It is observed that chaotic behavior came into
existence in several mathematical models, like, logistic map, delay model, Henon model, Rossler
model, Lorenz model and so on. Some useful ideas of mathematical models associated to chaos with
their applications are demonstrated here. It may attract mathematicians and scientists towards
chaos in mathematical modelling from different fields of their researches.

21
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Active Control with the Method of Receptances: Recent Progresses and

its Application in Active Aeroelastic Control
Kumar Singha
a
Associate Professor and Graduate Program Director 056 L GAR (Garland Hall) Dept. of
Mechanical and Manufacturing Engineering Miami University, Oxford, OH 45056
Email:[email protected]
URL: http://www.users.muohio.edu/singhkv/

About the speaker.

Dr. Kumar Vikram Singh is an Associate Professor and Graduate Pro-
gram Director in the Department of Mechanical and Manufacturing Engi-
neering. He has received his B.E. (1997) from Birla Institute in Technology,
Mesra, India in Mechanical Engineering and his Ph.D. in Mechanical En-
gineering from Louisiana State University (LSU) in May 2003. He joined
Miami in 2006 after serving as a post-doctoral research associate for three
years at LSU. After his Bachelors degree he was Graduate Engineer at Hin-
dustan Aluminum (HINDALCO) at Renukoot India. He has also held the
rank of visiting assistant professor at LSU (summer 2007), summer faculty
fellow at AFRL-WPAFB (2009- 2011, 2013-2014, 2016), visiting senior fel-
low at University of Liverpool, UK (2009) and visiting researcher at Aerospace Systems Directorate
(AFRL/RQVC), Wright Patterson Air Force Base (WPAFB), Dayton (2014). His research focuses
on inverse eigenvalue problems related to multidisciplinary areas of vibration, aeroelastic control,
structural health monitoring, rotordynamics, smart materials and composite structures. In partic-
ular his interests are in developing technologies for active control of aerospace structures, rotors and
viscoelastic materials, strategies for parameter estimation, design and system identification based
on dynamics or vibration data, and non-destructive testing and evaluation techniques for com-
posite structures. He is an author/co-author of more than 75 peer-reviewed international journal
and conference articles. His research is supported by National Science Foundation (NSF)/TUES,
USDA, and AFRL/DAGSI grants. He is an Associate Fellow of the American Institute of Aeronau-
tics and Astronautics (AIAA) and recipient of the Arthur Olson Teaching Excellence Award and
Outstanding Researcher Award in the College of Engineering. He is also a registered Professional
Engineering (P.E.) Licensee (Mechanical Design) in the state of Louisiana.
Abstract. The design of active control for large and complex engineering structures requires
accurate modeling and prediction of their dynamic response and instabilities. The performance
of traditional model based control may be limited due to the errors in model approximation,
size of the problem and/or availability of limited data for realizing active control. To overcome
some of these challenges the method of receptances is developed. This method allows design
and computation of controller gains based on modest size of receptance matrices which can be
extracted from transfer functions associated with available sensors and actuators. The area of
active aeroelastic control deals with developing wing technology for the next generation aircrafts
to achieve increased performance by controlling and manipulating the aeroelastic response by
active means. In these applications, receptance based controller design is found to be promising
as it eliminates the modeling of complex aeroelastic interaction between elastic structure and
surrounding aerodynamics. In this talk, fundamental of the receptance based control is introduced
and recent progress in this area is highlighted. The effectiveness of the controller designed with on-
board sensors (embedded) and actuators (control surfaces) for suppressing the flutter instabilities

22
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

and flutter boundary extension is demonstrated with numerical examples. The performance of the
controller such as its ability to control prescribed modes of interest without influencing the other
modes including the actuator modes is also presented. Ongoing research in this area and potential
for future research are briefly summarized.
Extension of Wavelets to Topologically Complicated Domains
Mani Mehraa
a
Dept. of Mathematics Indian Institute of Technology Delhi Hauz Khas, New Delhi-110016, India
Email: [email protected]
URL: http://web.iitd.ac.in/∼mmehra/

About the speaker.

Mani Mehra earned her M.Sc. degree in Applied Mathematics from In-
dian Institute of Technology Roorkee, India, July 2000 and Ph.D. degree in
Mathematics from the Indian Institute of Technology Kanpur, April 2005.
From 2005 to 2007, she was a post-doctoral fellow at the Department of
Mathematics, Mcmaster University, Canada. She has been with the Depart-
ment of Mathematics, IIT Delhi, since 2008. Her research interest includes
wavelet methods for PDEs, numerical methods, inverse problems.
Abstract. Differential equations on topologically complicated domains is a
relatively new branch in the theory of differential equations. Some of the
examples include differential equations on manifolds or irregularly shaped domains and differential
equation on network-like structure. Differential equations on manifolds arises in the areas of math-
ematical physics, fluid dynamics, image processing, medical imaging etc. Differential equations
on network-like structure also play a fundamental role in many problems in science and engineer-
ing. The aim of this talk is to show how wavelets could be extended to network to solve partial
differential equations on network-like structure using spectral graph wavelet.
Controllability of Fractional Differential Equations with
Non-instantaneous Impulses
Muslim Malika
a
School of Basic Sciences Indian Institute of Technology Mandi Kamand - 175005, H. P. India
Email:[email protected]

About the speaker.

Dr. M. Malik is working as an Assistant Professor of Mathematics in the
Department of Mathematics, Indian Institute of Technology Mandi, India.
He obtained his Ph.D. degree in 2006 from Indian Institute of Technology
Kanpur, India.
Abstract. Initially, we discussed about the fractional calculus and con-
trollability concept. Latter, we consider a control system governed by an
abstract fractional differential equation of order alpha lies between 1 and 2.
We able to prove the exact controllability of the problem by using the cosine
family of operators and Banach fixed point theorem. Finally, we have given
an example to illustrate the application of these results.

23
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Bio-marker Prediction using Bayesian State-Space Modeling in Gene

Expression Data Analysis
Gajendra K. Vishwakarmaa
a
Department of Applied Mathematics, Indian Institute of Technology Dhanbad, Dhanbad-826004,
India.
Email: [email protected]

About the speaker.

Dr. Gajendra K. Vishwakarma is working as an Assistant Professor
of Statistics in the Department of Applied Mathematics, Indian Institute
of Technology Dhanbad, India. He obtained his Ph.D. degree in 2007
from Vikram University, Ujjain, India. After obtaining his PhD degree, he
worked as an Associate Scientist in novel drug discovery and development,
Lupin Research Park, Pune, India and Cognizant Technology Solutions
Mumbai, India. Prior to joining IIT Dhanbad, he worked as visiting scien-
tist cum faculty in Indian Statistical Institute, North-East Centre, Tezpur
(Assam), India.
Abstract. The advancement in computational biology and statistical modeling help to identify
the genes which causes the disease like cancer by comparing its expression levels in diseased and
healthy people. Bayesian state space modeling is a new advancement in statistics which can
estimate an unobserved values of a process using the information from an observed outcome and
its relationship. Using this two ideas together, we are trying to model and estimate the expression
values of genes which are significantly different among two groups. The complicated integration
of posterior densities are carried out using the Markov Chain Monte Carlo (MCMC) simulations.
The study shed light on the use of Bayesian State Space modeling to elucidate the behavior of
Bio-markers.

Computing Nearest Stable Matrices

Punit Sharmaa and Nicolas Gillis b

a
Department of Mathematics, Indian Institute of Technology Delhi, India b Department of
Mathematics and Operational Research, University of Mons, Belgium.

Email: [email protected]
About the speaker.
Dr. Punit Sharma is working as an Assistant Professor of Mathematics
in the Department of Mathematics, Indian Institute of Technology Delhi,
India. He obtained his Ph.D. degree in 2015 from Indian Institute of Tech-
nology Guwahati, India.
Abstract. The stability of a continuous linear time-invariant (LTI) system
ẋ = Ax + Bu, where A ∈ Rn,n , B ∈ Rn,m solely depends on the eigenvalues
of the stable matrix A. Such a system is stable if all eigenvalues of A are
in the closed left half of the complex plane and all eigenvalues on the
imaginary axis are semisimple. It is important to know that when an
unstable LTI system becomes stable, i.e. when it has all eigenvalues in the stability region, or how

24
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

much it has to be perturb to be on this boundary. For control systems this distance to stability is
well-understood. This is the converse problem of the distance to instability, where a stable matrix
A is given and one looks for the smallest perturbation that moves an eigenvalue outside the stability
region. In this talk, I will talk about the distance to stability problem for LTI control systems.
Motivated by the structure of dissipative-Hamiltonian systems, we define the DH matrix: a
matrix A ∈ Rn,n is said to be a DH matrix if A = (J − R)Q for some matrices J, R, Q ∈ Rn,n
such that J is skew-symmetric, R is symmetric positive semidefinite and Q is symmetric positive
definite. We will show that a system is stable if and only if its state matrix is a DH matrix. This
results in an equivalent optimization problem with a simple convex feasible set. We propose new
algorithms to solve this problem. Finally, we show the effectiveness of our method compared to
the other approaches and to several state-of-the-art algorithms. These ideas can be generalized to
get good approximate solutions to some other nearness problems for control systems like, distance
to stability for descriptor systems, distance to positive realness, and minimizing the norm of static
feedback.

25
International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Contributed Talks

26
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Water boundary layer flow over an exponentially permeable stretching

sheet with variable viscosity and Prandtl number
Abhishek Kumar Singha , N. Govindaraja , S. Royb
a
Mathematics Division, School of Advanced Science, VIT university, Chennai Campus, Chennai
- 6000127, India.
b
Department of Mathematics, Indian Institute of Technology Madras, Chennai - 600036, India
Email: [email protected] (Abhishek Kumar Singh ),
[email protected] (N. Govindaraj), [email protected] (S. Roy ).

Abstract. The present work focus on water boundary layer flow over an exponential permeable
stretching sheet in the presence of suction/injection with variable viscosity and prandtl number.
The nonlinear partial differential equations governing flow and thermal fields are presented in
non-dimensional form using suitable non-similar transformation. Finally non dimensional partial
differential equations are solved by implicit finite difference method in combination with the Quasi-
linearization technique. The numerical results for skin-friction and local Nusselt number are shown
graphically to display effects of physical parameters.

Mathematical Modelling Of Solar Cell

Manjusha Tomara , Pradip Jamkara , Parvin Shaikha
a
Indira college of Engineering and Management, Pune
Email: [email protected] (Manjusha Tomar)

Abstract. This paper presents a mathematical modelling for solar cell with diode in parallel and
series arrangement for photovoltaic effect. The analysis of Current due to voltage, temperature
and other parameters when diodes are connected in series and parallel is studied. MATLAB is
used to develop and present the model. Results are in good agreement with the theory.

Reliability assessment of Transformer incorporating different parameters

Neelam Chantolaa and S. B. Singhb
a
Department of Mathematics, Statistics and Computer Science,
b
G.B. Pant University of Agriculture and Technology, Pantnagar, India
Email: [email protected] (Neelam Chantola)

Abstract. Transformer is one of the most common equipment’s in power system which plays
an important role in providing reliable and efficient electricity. This paper studies the mathe-
matical modelling to analyze reliability of transformer using insulation paper and oil along with
their parameters. Insulation paper is affected by breakdown voltage, moisture content, degree of
polymerization, furanic content whereas insulation oil is affected by breakdown voltage, moisture
content, dissipation factor and inter-facial tension. If any one of the parameters of insulation
oil/paper fails then transformer either functions partially or stops working completely. Various

27
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

reliability characteristics such as state transition diagram, reliability and mean time to failure
(MTTF) of the proposed system (transformer) have been evaluated with the help of supplementary
variable technique, Laplace transformer, Gumbel-Hougaard copula. A numerical example is taken
at last to illustrate the findings.

Mathematical Modeling of Moored Ship Motion in Vishakhapatnam

Port under the Resonance Conditions
Prashant Kumara , Gulshana , Kwang Ik Kim b

a
Department of Applied Sciences (Mathematics), National Institute of Technology,Delhi-110040,
India.
b
Department of Mathematics, Pohang University of Science and Technology,Pohang, South
Korea.
Email: [email protected](Prashant Kumar),
[email protected](Gulshan), [email protected](Kwang Ik Kim)

Abstract. An accurate description of harbor geometry, bathymetry and incident wave phe-
nomenon associated with diffraction, reflection and refraction is required to develop a numerical
model for the analysis of moored ship motion in realistic ports. Most of realistic port connected to
the open experienced extreme oscillation due to high amplitude incident waves under the seasonal
weather conditions. A numerical model is developed based on Helmholtz equation, which solved by
using boundary element method (BEM) in the context of the shallow water waves. A moored ship
motion is described with the six different modes, i.e., surge, sway, heave, roll, pitch and yaw motion.
Simulation results are compared with experimental results obtained by Sawaragi and Takagi for
moored ship motion. Further, the current numerical scheme is implemented on realistic Vishakha-
patnam port, India to analyze the moored ship motion under the resonance conditions. The size
of moored ship is considered 250,000 DWT (Dead Weight Tonnage) inside the Vishakhapatnam
port, India. The countermeasure has been discussed to reduce the wave oscillation at various port
stations.

Stability analysis of a film flow down an incline in the presence of a

floating flexible membrane
S. Ghosha , H. Beheraa , M. Sania
SRM Research Institute and Department of Mathematics,
SRM Institute of Science and Technology,
Kattankulathur, Tamil Nadu, India
Email: [email protected] (S. Ghosh), [email protected]
(H.Behera),[email protected] (M. Sani)

Abstract. The present study deals with the effects of floating flexible membrane on the instabil-
ity of a gravity-driven flow down an incline. Linear stability of the flow system is explored using
normal-mode analysis. Free surface gravity-driven flow is unstable at much lower Reynolds num-
bers. Instability of such a flow can be controlled either by changing behavioral of the lower wall or
by altering the surface waves at the free surface which is done here by including a floating flexible
membrane at the top of the liquid layer. Influence of membrane tension is taken into account in

28
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

terms of stress jump at the free surface. The Orr-Sommerfeld system of the flow is solved nu-
merically using spectral collocation method. The results displays a destabilizing role of membrane
tension for a wide range of parameters. The growth rate of the perturbation waves increases with
an increase of membrane tension and the critical Reynolds number becomes smaller. Therefore, it
is possible to enhance the instability of the flow system with help of membrane properties, which
may be useful in Ocean engineering and coating industries.
Keywords Hydrodynamic stability; film flow; floating membrane; Orr-Sommerfeld analysis.

Mathematical Modeling and Aalysis of Single Walled Carbon Nanotube

Based Mass Sensor using Molecular Structural Mechanics Approach
Harsh Sharmaa , Sandesh Trivedib
a
M. Tech. Scholar, Department of Mechanical Engineering, Govt. Engg. Col.,Ajmer
b
Assistant Professor , Department of Mechanical Engineering, Govt. Engg. Col.,Ajmer
Email:a [email protected],b [email protected]

Abstract. The suitability of single walled carbon Nanotube for detection of mass as small as
of the order of femtogram level (10 -18 kg) has been explored in the present paper using finite
element method based three dimensional atomistic modeling approach. The effect of varying mass
on frequency has been analyzed by considering armchair (5,5) carbon Nanotube for bridged con-
figuration. Further this paper analyses the effect of attached mass lending position on Nanotube.
Based on FEM simulation the mass sensitivity for bridged configuration is found to be higher as
compared to cantilever one. The simulation results revels that shorter length of CNT with bridged
configuration is suitable for sensing mass as small as of the order of femtogram level.
Keywords Carbon Nanotube, Finite Element Method, Mass Sensor, Reson- ance Frequency, Chi-
rality.
Simplified multivariate GARCH (SGARCH) model: A survey
Anuja Guptaa
a
Department of Mathematics, Maulana Azad National Institute of Technology (MANIT),Bhopal
(M.P.) 462007, India
Email: [email protected]

Abstract. Time series models play an important role in portfolio optimization. Various linear
and non-linear time series models were presented in the literature for portfolio optimization. Hedg-
ing is also a portfolio optimization technique in which investment was done to reduce the risk of
adverse price movements in an asset. Some multivariate GARCH models used in the literature
for hedging were CCC (constant conditional correlation), DCC (dynamic conditional correlation),
BEKK (Baba-Engle-Kraft and Kroner), EGARCH (Exponential GARCH), Diagonal VECH etc.
SGARCH was also one of them. This paper contains a survey of the simplified multivariate GARCH
(SGARCH) model of conditional heteroskedasticity. The simplified multivariate GARCH model
and its applications in various stocks and commodities markets are presented. This paper re-
views the model, hedging strategies, applications, and identifies likely directions of future research.
Keywords: Hedging, Hedge ratio and hedging effectiveness, SGARCH model

29
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A Greens function approach to analyze the dispersion characteristics of

Love type wave due to an impulsive point source in a piezoelectric
layered structure
Anusree Raya , Abhishek K. Singhb
a
Junior Research Fellow, Department of Applied Mathematics, Indian Institute of Technology
(Indian School of Mines), Dhanbad, Jharkhand-826004
b
Assistant Professor, Department of Applied Mathematics, Indian Institute of Technology
(Indian School of Mines), Dhanbad, Jharkhand-826004

Abstract. An external source of disturbance in a material, even of point size, give rise to waves
propagating away from the concerned region in its interior at a specified time. Such an impulse
may be best described with the aid of Dirac delta function. Greens function is primarily utilized in
solving these problems of elastodynamics. The present study focuses to investigate the propagation
characteristics of Love-type wave influenced by an impulsive point source in a layered structure
comprised of a heterogeneous piezoelectric layer lying over a heterogeneous isotropic half-space.
Greens function technique is adopted in order to obtain the dispersion equation, which is further
reduced to the classical result of Love wave. For sake of computation, numerical data of PZT-5H
ceramics for the heterogeneous piezoelectric layer is considered. Influence of heterogeneity, piezo-
electricity and dielectric constant associated with the heterogeneous piezoelectric layer; and effect
of heterogeneity parameter and corresponding magnification factor concerned with heterogeneity
in the isotropic half-space has been reported through graphical delineation.

Sustainable and Smart Manufacturing through Distributed

Manufacturing System and Robotics Methodology
Nandini Vishwakarmaa , Kanchan Singhb , Anshul Yadavb , Rajnish Singhb
a
Department of Mechanical Engineering, B.B.S College of Engineering Allahabad, India
b
Department of Mechanical Engineering, KNIT, Sultanpur, India
Email: [email protected](Anshul Yadav)

Abstract. Nowadays, manufacturing companies are working on the sustainable and smart man-
ufacturing systems. In this paper sustainable manufacturing initiatives is the creation of energy-
efficient and manufacturing resources-saving system s. Distributed Decentralized networks of
adaptable and flexible mini-factories are not only helpful to reduce emissions through reduction of
transports, but also serve for the growth and development of regional economic cycles. This work
also makes an attempt at exploring the role of robotics in achieving sustainability in manufac-
turing. Its describes the advancement made so far in applying robots to manufacturing- welding
shot blasting, painting, and so forth,within the sustainability framework. The adoption of new
technologies like Robotic Operating Systems (ROS), in the approaches to the design of industrial
robots, could enhance sustainability in the use of robots in manufacturing at a reduced cost. This
paper gives an overview of trends towards DMS as well as reasons and arguments why Distributed
Manufacturing Systems and Robotics are appropriate concepts for more sustainable and smart
manufacturing.

30
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Development of a dimensionally reduced model for a multilayered

multi-elastic structure using the Variational Asymptotic Method (VAM)
Arun Krishnana , Dineshkumar Harursampathb ,V Eswarana
a
Vikram Sarabhai Space Centre, ISRO, Trivandrum,
b
Department of Aerospace Engineering, IISc, Bangalore
a
k [email protected],b [email protected], a v [email protected]

Abstract. Multi-layered structures are increasingly being used in various fields of engineering.
Many of these structures consist of combinations of hyper-elastic as well as linear elastic mate-
rials whose strain capabilities are vastly different (multi-elastic). 3-D Finite element analyses of
these structures are often computationally intensive and error prone, especially when rubber like
materials are encountered. Variational Asymptotic Method (VAM) is a mathematical technique
based on the calculus of variations. By utilizing the concept of minimizing the total potential
energy, VAM is used to develop asymptotically correct and computationally efficient dimensionally
reduced models for 3-D structures for which small parameters associated with the geometry or
material can be identified. The current work applies VAM to a two-layered structure consisting of
a combination of hyper-elastic and linear elastic materials, in which thickness is small compared to
the other two dimensions. Starting from geometrically exact kinematics, the 3-D problem is split
into a 1-D through thickness analysis and a 2-D reference surface analysis by the application of
VAM. Closed form expressions for the distribution of 3-D strains through the thickness dimension
are obtained in a mathematically consistent manner. Computational efficiency is achieved by per-
forming the 2-D analysis in lieu of the 3-D analysis, followed by recovery of the 3-D fields using
the recovery relations. The current work paves the way for development of dimensionally reduced
models for real-life multi-elastic 3-D structures, like elastomeric bearings and flex seals employed
in solid rocket thrust vector control.

Unsteady stagnation point nanofluid flow with MHD, thermal radiation

and viscous dissipation effects over a stretching sheet
G.S.Setha and B.Kumara
a
IIT (ISM), Dhanbad
Email: [email protected]

Abstract. The present paper is framed to study the magnetohydrodynamic stagnation point flow
and heat transfer of a nanofluid influenced by thermal radiation. Nonlinear ordinary differential
equations are obtained by applying suitable similarity transformation on the governing partial
differential equations. These differential equations have been solved by successive linearization
method (SLM). The effects of various active flow parameters, such as thermal radiation parameter,
stagnation parameter etc. on the flow and heat transfer, tabulated and plotted graphically by
accumulating sufficient data using SLM. Apart from this, quadratic multiple regression analysis is
carried out on numeric data of skin friction coefficient, Nusselt number. Finally obtained results
are compared with previously available results in limiting sense.

31
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Simulation of Livestock Population in India Using Exponential And

Logistic Models
C.Vipina , R.Krishikaa , V. Kuzhal Mozhia

J.Vijayarangamb , M.Sujithrab
a
MSc Data Science, Department of Computing,Coimbatore Institute of Technology.
b
Department of Computing Coimbatore Institute of Technology
Email:[email protected]

Abstract. The aim of the study is to simulate the livestock population of India by using exponen-
tial and logistic models respectively and to compare the actual population values and simulated
population values. In the Livestock Population Cattle, Buffaloes, Sheep, Pigs had been taken for
the simulation. The years that had been taken for the simulation is from 1951 to 2012 .The Mean
Square error is been calculated for individual animal population data for both exponential and
logistic population models and by using that the appropriate model is founded.
Stability Analysis For An Seiqr Epidemic Model With Satur ated
Incidence Rate
Deepti Mokatia , Nirmala Guptab , V.H. Badshaha
a
School of Studies in Mathematics Vikram University, Ujjain (M.P.)
b
Govt. Girls P.G. College, Ujjain(M.P.),India
Email: [email protected]

Abstract. Mathematics plays an important role in study of biological systems through mathemat-
ical models. In the present pap er , we extended the work of Nirwani N. et al. [5] by introducing
the transmission rate from the compartment to and converted the model into an SEIQR epidemic
model with saturated incidence rate . Determine the two equilibrium points : disease-free and
endemic of the model and basic reproduction number is obtained . Stability analysis have bee n
discussed of both equilibrium points by Routh-Hurwitz criteria and Lyapunov function criteria .
Also , Numerical simulations are carried out for the model .
Mathematical modelling of UASB reactor: A detailed review
Varsha Parasara and Dilip K Rajakb
a
Assistant Professor, Department of Petrochemical Engineering,
b
SATI(Degree),Vidisha
Email:[email protected]

Abstract. The main objective of this paper is to review the work already done related to the
mathematical modelling of Up flow anaerobic sludge blanket (UASB) in which substrate degra-
dation, biomass growth and death with time, washout of microorganisms and reaction within the
granules have taken place. Microbial growth was estimated by the empirical model from which the
parameters Ks (substrate concentration) with m (maximum specific growth rate) were predicted
by the theory and comparison with the lab scale experimental data from which mass balance for
substrate degradation and microorganism growth were done that allow to study the consequences
for the reactor performance of different types of substrate, different inlet substrate concentration,
different flow rates and different kinds of biomass.

32
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Modelling of Quantum Key Distribution Protocols in Communicating

Quantum Processes language with Verification and Analyses in PRISM
Satya Kuppama
a
University of Massachusetts Amherst
Email: [email protected]

Abstract. Proof of security of cryptography protocols theoretically establishes the strength of

a protocol and the constraints under which it can perform, it does not take into account the
overall design of the protocol. In the past model checking has been successfully applied to classical
cryptography protocols to weed out design flaws which would have otherwise gone unnoticed.
Quantum key distribution protocols differ from their classical counterparts, in their ability to
detect the presence of an eavesdropper while exchanging the key. Although unconditional security
has been proven for both BB84 and B92 key distribution protocols, in this paper we show that
identifying an eavesdroppers presence is constrained on the number of qubits exchanged. We
first model the protocols in Communicating Quantum Processes (CQP) and then explain the
mechanism by which we have translated this into a PRISM model and how we analysed the
protocols capabilities. We mainly focus on the protocols ability to detect an active eavesdropper
and the extent to which an eavesdropper can retrieve the shared key without being detected by
either party. We then conclude by comparing the performance of the protocols.
The role of adaptation in an aquatic predator-prey system with
Beddington-DeAngelis type functional response
Nilesh Kumar Thakura , Archana Ojhaa , S. K. Tiwarib
a
Department of Mathematics, National Institute of Technology, Raipur 492010, India
b
Department of Mathematics, BIT Sindri, Dhanbad 828123, India
E-mail:[email protected]; [email protected];
[email protected]

Abstract. In this paper two interacting species in presence of adaptation (dormancy of the preda-
tors such as resting eggs) has been discussed. The dormant stage is an equipment to survive in
harsh environment. We have discussed the stability analysis of system without diffusion and in
presence of diffusion. Our numerical investigation reveals that above the critical value of interfer-
ence among the zoo plankton the system become stable. Spatiotemporal pattern shows a transient
complex spatiotemporal pattern by increasing the time and space.

Analysis of Time-Delayed Non-Linear Equations Using HF Functions

Shailesh Vermaa , Aniruddha Ghosha and Anindita Gangulya
a
Indian Institute of Technology, Dhanbad, India

Abstract. The paper deals with the analysis of non-linear time delayed differential equations
solved using HF functions. The analysis is first performed on Mackey-Glass Equation, which is
a standard model for quantitative characterization of chaotic dynamics. The procedure is then
performed on a generalized Human respiratory control model, where for different simulation pa-
rameters the analysis of Cheyne-Stokes Breathing is done. Both models are simulated in MATLAB.
The graphs thus generated are used to provide suitable conclusions.

33
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Controllability Results For a Nonlocal Nonlinear Fractional Functional

Differential Equation
Madhukant Sharmaa , Shruti Dubeyb
a
Mahindra Ecole Centrale Hyderabad,b Indian Institute of Technology Madras
Email: [email protected]
Abstract. This paper is consecrated to study the controllability of solutions to non autonomous
fractional functional differential equations in a general Banach space along with a nonlocal con-
dition. The main results are proved by using fixed point techniques, classical semi-group theory
and tools of fractional calculus. As an example, a nonlocal retarded elliptic evolution equation of
fractional order has been given for demonstration.
Effect of the heterogeneity, initial stress and viscosity on the
propagation characteristics of shear wave
Raju Kumhara and Santimoy Kundua
a
Department of Applied Mathematics
Indian Institute of Technology (Indian School of Mines)
Dhanbad, Jharkhand-826004, India
Email: [email protected]
Abstract.This paper irradiates the influence of initial stress, heterogeneity and viscoelasticity on
the propagation behavior of shear wave in a heterogeneous orthotropic layer under initial stress over
a heterogeneous visco-elstic half-space. A variable separable method has been adopted to obtain
the analytical solutions for the layer and half-space separately. A complex frequency equation has
been derived by using suitable boundary conditions. Thereafter, the complex frequency equation
has been separated into dispersion and absorption equations. Effects of heterogeneities, initial
stress and attenuation coefficient on the phase velocity and dissipation function of shear wave have
been shown graphically.

Study of Memory Effect in an Inventory Model With Quadratic Type

Demand Rate and no Shortage
Rituparna Pakhiraa , Uttam Ghosha and Susmita Sarkara
a
Department of Applied Mathematics, University of Calcutta,Kolkata,India
Email: [email protected]
Abstract.In this paper we have developed a memory dependent inventory system using fractional
derivative.The model has been formulated with quadratic type demand rate and no shortage. The
minimized total average costs and the optimal ordering interval have been found for this fractional
order model i.e.memory dependent inventory model using fractional order integration. Those
quantities depend on the fractional order of derivative of the inventory level (i.e. on the memory
parameter). Finally the numerical results are presented graphically to justify our motivation.

34
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Switching stability of a delayed e-epidemic model in wireless sensor

network
Sangeeta Kumaria and Ranjit Kumar Upadhyaya
a
Department of Applied Mathematics
Indian Institute of Technology (Indian School of Mines)
Dhanbad, Jharkhand-826004, India
Email: [email protected]
Abstract. We have formulated an e-epidemic energy efficient susceptible-infected-terminally
infected-recovered (SITR) model to analyse the attacking behaviour of worms in wireless sensor
network using cyrtoid type functional response. In this model, once a sensor node has been
attacked by the worms, the terminally infected node spreads the worms to its neighbouring nodes
using normal communications, which further spread it to their neighbouring nodes and the process
continues. To tackle this issue, we have propose an SITR model by considering the sleep mode
concept of wireless sensor network in which the operational capabilities and power consumption
of the motes decreases. Stability and bifurcation analysis are analysed for the proposed model
system with or without delay. The impact of the control parameters on the system dynamics are
investigated. Numerical study of the delayed system shows double Hopf bifurcation dynamics.
Double Hopf bifurcation causes stability switching i.e., instability to stability to instability, or the
reverse transition of the solution of the proposed system, which shows a fragile security system of
wireless sensor network.

Mathematical Modeling and Analysis of Green Sand Casting Process

Parameters Based On Taguchi Design of Experiments
Sanjay Alonia , S. S. Chaudharib and R. L. Shrivastavc
a
PhD Scholar,b Head & Asso. Professor,c Professor
Department of Mechanical Engineering, Yeshwantrao Chavan College of
Engineering,Maharashtra, India
Email: [email protected]
Abstract. Green sand casting process is one of the most popular sand casting processes for man-
ufacturing cast components using machine moulding in medium and large-scale foundries. Wide
range of automotive castings and castings required in other engineering applications is produced
in various grades of cast irons by a green sand casting process. The quality of castings produced
in green sand casting process, is always a result of a proper combination and settings of a large
number of parameters in the process. The main objective of this work is to develop a mathemat-
ical model of green sand casting process to predict the occurrences of the persisting defects that
are shrinkage porosity, blow holes,sand inclusion, scabbing and low hardness in gray iron cylinder
head components. Regression analysis is used to investigate the relationship between input and
responses.19 potential process parameters selected are taken as input and persisting defects are as
responses. Experiments are conducted as per the plan given by Taguchi L36 OA. The outcome of
the study clearly indicated that mathematical models developed are helpful for foundry industries
to envisage the defects under different combination of parameters.

35
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Modelling of Avalanches and Streamers in Resistive Plate Chambers

Jaydeep Dattaa , Sridhar Tripathya and Nayana Majumdara
a
Homi Bhabha National Institute, SINP, INO, India
Email:[email protected]
Abstract. India-based Neutrino Observatory (INO) is a proposed underground facility to host
a magnetized Iron Calorimeter (ICAL) (doi.org/10.1007/s12043-017-1373-4) for studying atmo-
spheric neutrino properties using Resistive Plate Chambers (RPC). RPCs, gaseous detectors with
glass or bakelite electrodes, deliver position and timing information of propagating muons. Charged
particles passing through the RPC gas gap produce (primary) electron ion pairs due to the ioniza-
tion of gas molecules. These primary electrons, under the influence of applied electric field, gather
energy and create further secondaries. The extent of the growth of secondaries within the RPC
determines avalanche or streamer mode of operation. Due to the high Global Warming Potential
(GWP 1430) of the present gas mixture (R134a+isobutane+SF6), a more eco-friendly one capable
of working in avalanche mode is being sought for. In this work,we have tried to model avalanches
and streamers inside an RPC using a hydrodynamic model. The hydrodynamic equations for the
electrons and ions and the electrostatic equations for the given system have been solved by us-
ing COMSOL Multiphysics (comsol.co.in). The number of electrons and ions and the resulting
modified electric field have been calculated for the necessary time interval. The study fulfills the
purpose of finding out the working voltage region for operating an RPC in pure avalanche mode
for a given gas mixture.

Analysis and Comparative Study of Various Performance Measures of

M/G/1 and M/G/S Queuing Models
Pradeep K. Joshia , Shejal Gupta b
and Dr.K.N.Rajeshwaric
a
Reader, Department of Mathematics, IPS Academy, Indore, India,b Research Scholar, School of
Mathematics, DAVV, Indore, India,c Professor Head, School of Mathematics, DAVV, Indore,
India
Email:[email protected]

Abstract. Queuing theory is a mathematical tool, which deals with the phenomenon of waiting
lines. Keeping customers waiting for too long, affects various parameters of queuing theory. In this
paper, optimization technique has been applied to calculate the performance measures like waiting
time, traffic intensity, waiting cost and proportion of time the server is idle for single server M/G/1
model as well as for multiple servers M/G/S model. Also, we analysis the comparative study of
various performance measures of single server M/G/1 model and multiple servers M/G/S queuing
model. This investigation deals with M/G/1 and M/G/S Models in which arrival follows Poisson
distribution, service time follows general distribution, and Queue discipline is FCFS. The result of
analysis showed that as we increase the number of servers, the waiting time in the queue decreases,
total cost of the system decreases, traffic intensity also decreases whereas idle time increases. The
results are optimize and effective. The graphical representation of various performance measures
and their comparison are shown by MATLAB software.

36
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

On the Dynamics of Transverse Domain Walls in Ferromagnetic

Heterostructures with Nonlinear Dissipative Effects
Sharad Dwivedi
SRM Research Institute & Department of Mathematics, SRM IST, Che nnai, India.
[email protected]
Abstract. This work delineates the investigation of static and dynamical behavior of transverse
domain walls in ferromagnetic heterostructure in which the thin ferromagnetic layer is sandwiched
between a heavy metal and a metal oxide layers. We consider the metallic ferromagnet with high
perpendicular magnetic anisotropy and exhibits structural inversion asymmetry. The presence of
a strong inversion asymmetry in the ferromagnetic material leads to additional spin-orbit coupling
torque terms into the dynamics qualitatively different from the standard spin-transfer torques.
The evolution of magnetization inside the ferromagnetic layer is governed by the one-dimensional
model of generalized Landau-Lifshitz-Gilbert-Slonczewski equation consisting of the standard spin
transfer and spin-orbit coupling torque terms along with the nonlinear dissipation factors viz.
viscous and dry-friction. Under this framework, first, we establish the static magnetization profile
in the two faraway domains in the presence of transverse magnetic field once the electric and
longitudinal magnetic fields are switched off. Then, by means of regular perturbation expansion
met hod, we derive the zero order traveling wave solutions under the influence of small applied
magnetic field and spin-polarized electric currents.

Numerical study of dragonfly wing section in Gliding Mode

Mohd imran Ansaria and Syed Fahad Anwerb
a
Department of mechanical Engineering, Zakir Husain College of Engineering and Technology,
b
Aligarh Muslim University
Email: [email protected]

Abstract. The computational study of insect dragonfly Aeshna cyanea forewing is carried out
for aerodynamic performance and compare it with smooth airfoil and flat plate. For the study of
wing, a 2d section of Aeshna cyanea at mid of the wing, which has pleated surface is considered.
The smooth airfoil is constructed by drawing a smooth curve from the edge of corrugated wing.
A range of angle of attack with 5◦ increment is considered from 0 to 30 and Re=1400, which is
actual gliding Reynolds number of dragonflys. The unsteady form of Navier-Stokes equation for
incompressible Newtonian flow was solved using fully implicit NITA scheme for time advancement.
The convective terms were evaluated using Quadratic Upstream of Convection kinetics (QUICK).
The solution is obtained using fractional step method as pressure and velocity coupling scheme
on the staggered structured grid. It is found that there is the reduction in shear drag for pleated
airfoil. The pleated airfoil and flat plate have same lift up to 25?, thereafter flat plate have highest
lift at 30◦ angle of attack. The best performance in terms of gliding ratio is found for flat plate
airfoil at 5? angle of attack.

37
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Ball convergence theorems for Chebyshev’s method in Banach space

Chandni Kumaria and Pradip Kumar Paridaa
Central University Of Jharkhand, Brambe, Ranchi-835205
Email: [email protected]

Abstract. This article is concerned with the local convergence analysis of Chebyshev’s method
in order to approximate a locally unique zero of an nonlinear operator in Banach space. Under
some new type of majorant condition, a new local convergence analysis for Chebyshev’s method
is presented. Moreover, we will present two important special cases about the convergence results
based on the Kantrovich-type and Smale-type assumptions.

Optimal multi-echelon integrated supply chain selecting best supplier

and distributor using multi-objective genetic algorithm
Isha Taltia and Poonam Mishra a

a
Department of Mathematics, School of technology, Pandit Deendayal Petroleum University,
Raisan, Gandhinagar, 382007, India
Email: [email protected]

Abstract. Over last few decades, mathematics has played a crucial role in developing efficient
algorithms for Face Recognition (FR) used in biometric systems. FR using Machine Learning
(ML) techniques has impacted FR systems tremendously, towards efficient and accurate models
for FR. Existing FR systems used in biometrics use ML techniques to learn patterns in the images
by extracting various features from them and often require pre-processed face image data for the
learning process. In this paper, we have used various pre-processing techniques and compared them
in the deployed FR framework. It was observed that the Steerable Pyramid (SP) filter was the
most efficient pre-processing technique among all techniques used for pre-processing in this work.
Though existing feature extraction methods such as SIF T (Scale-Invariant Feature Transform),
SURF (Speeded-Up Robust Features), ORB (Oriented FAST and Rotated BRIE F) have been
used in the past, they have not been accurate enough in various vision based biometric systems.
Hence, a novel PSI (Pose Scale and Illumination) invariant SURF-RootSIFT technique is proposed
by extending the well kn own SIFT-RootSIFT feature extraction technique which is achieved by
calculating the Bhattacharyya Cofficient between the feature vectors. A framework which uses the
proposed novel feature extraction technique is deployed in this work. This paper demonstrates that
the novel SURF-Root SIFT based framework is proven to perform more accurately and efficiently
than the other techniques, with 99.65%, 99.74% and 97.93?curacy on the Grimace, Faces95 and
Faces96 databases respectively.

38
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

PSI (?) invariant features for face recognition

Ajaykumar S Cholina , A. Vinaya , Aditya D Bhata , Arnav Ajay Deshpandea and K N
Balasubramanya Murthya , S Natarajana
a
PES University
Email: [email protected]

Abstract. The study is to compare the efficiency of the trend models with that of neural networks
models and to identify which model is well suited for the actual and the simulated data set using
the mean square error values. The trend is fitted for the data. The data set is divided into test
and train data set, then the regression is performed and the test data set is also predicted and
the mean square error term also found to identify how much the prediction gets vary from the
actually obtained values. Same procedure is carried in terms of fitting the neural networks model.
The neural network model is fitted using the build in functions in R and based on the results
inferences are made. Also the model that is well suited for the given data set is found. The second
part of the study is to simulate the data. So, understanding the nature of data and based on the
column descriptions of the data, appropriate simulation is done. Then the simulated data set is
normalized for processing. Then the above procedure is followed to the simulated dataset whether
which model gets well suited.

Position analysis of five bar mechanism under force

Akashkumar Singh Rajputa and H.P Jawalea
a
VNIT nagpur
Email: [email protected]

Abstract. This paper presents the position analysis of a planer five bar mechanism, when the
mechanism is taken as symmetrical. Forward kinematics is used to determine the position analysis
and workspace of mechanism. In this paper the links has considered as elastic . This method used
when the mechanism undergoes to some external disturbance such as force, moment and torque.
When the force is applied on the mechanism, the mechanism under goes some deflection because
of considering link to be elastic and then comparing the deflection with rigid link. Here this paper
also includes the significances of mechanism when the link would be considered as elastic link and
rigid link. Due to that property there is the some deviation and this paper represents the behaviour
of mechanism and analysis of position of mechanism under deflection.

Family of Perfect fluid Balls of Class One in General Relativity

Amit Kumar Prasada and Jitendra Kumara
a
Central University of Jharkhand
Email: [email protected]

Abstract. In this paper, we consider a spherical symmetric metric in 5D flat form. We have
derived non conformal accelerating perfect fluid solutions to Einstein field equations by applying
similarity transformations method. The pressure and de nsity are monotonically decreasing with
increasing radius ’r’.The value of metric potentialg( 11) and g( 44) is non-zero positive.The expression
of fluid parameter is monotonically decreasing with increasing radius ’r’.

39
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Application of Integral transform in the solution of Fractional Kinetic

Equations
a
Garima Agarwal
a
Manipal university Jaipur
Email: [email protected]

Abstract. In the present paper we use pathway transform to find the solution of Fractional Kinetic
equations by using τ -Gauss Hypergeometric functions. This transform contains many classes of
transforms including Laplace Transform. Fur ther we solve the Fractional Kinetic equations by
using derivative of τ -Gauss Hypergeometric functions.

Investigation of Transient stability in IEEE 9- bus system using Power

World Simulator
Ravishankar Tiwaria , Mayank Goyala and Prashant Prakasha
a
GLA University
Email: [email protected]

Abstract. With the day-by-day increase in power consumption, the transmission & distribution
networks are operating on the verge of stability. The transient stability studies is one of the
important study that helps to maintain the more reliable power supply with increasing demand in
complex network configuration. This paper helps to understand the problem of transient stability,
its effect on the power system. In this paper IEEE 9- bus system is considered to analyze a list
of contingencies & its affect on the system. The contingency analysis is done using a powerful
simulation tool for different generator bus.

Certain Subclasses of Bi-starlike Functions Associated with q-derivative

Operator
P N Kamblea and M G Shrigana
a
Dr D Y Patil SOET
Email: [email protected]

Abstract. Making use of q-derivative operator, in this paper, we introduce two new subclasses
of the function class Σ of normalized analytic and bi-univalent functions defined in the open disk
U. Furthermore, we find estim ates on the first two Taylor-Maclaurin coefficients |a2 | and |a3 | .
Moreover, we obtain Fekete-Szegö inequalities for the new function classes.

40
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

The fruit-basket proof for the theorem on disconnected graph stating

the maximum number of edges it can have and its application in viral
marketing
Mohd. Hasnain Meera and Shivram Sharmaa
a
Medi-Caps University, Indore
Email:[email protected]

Abstract. This study deals with the proof of a theorem on graph connectivity that tells about the
maximum number of edges that a simple disconnected graph with n vertices’s and k components
can have. The approach that this proof follows is quite logical and lucid that uses some basic
concepts of the combination theory along with the fundamental principle of counting. The proof
begins with an interesting problem on fruits and baskets, analogous to the constraints of th e
theorem under consideration, and then proceeds in a way that leads us to the final result that
is aimed to be prove d. The methodology employed in proving the theorem is such that it not
only proves the theorem but also gives the physical orientation of the graph depicting how many
vertices’s does each component have and how are they connected to each other. Ultimately, the
result that is obtained is used to propose a model for increasing the profit of a firm via viral
marketing.

Aggregation of LARK Vectors for Facial Image Classification

A Vinaya , Vinayaka R Kamatha , Varun M, Nidheesha and S Natarajana , KNB
Murthya
a
PES University
Email: [email protected]

Abstract. Face recognition is prevailing to be a key aspect wherever there is a need for interaction
between humans and machines. This can be achieved by containing a set of sketches for all the
possible individuals and then cross-validating at necessary circumstances. We propose a mechanism
to fulfill this task which is centered on locally adaptive regression kernels. A comparative study
has been presented at encoding stages as well as at the classification stages of the pipeline. The
results are cautiously examined and analyzed to deduce the best mechanism out of the proposed
methodologies. All the ideologies have been tested for multiple iterations on benchmark datasets
like ORL, grimace and faces 95. The vectorized descriptors have been subjected to encoding using
slightly refined methods of feature aggregation and clustering to assist classifiers in imputing the
test subjects to their respective classes. The encoded vectors are classified using Gaussian Naive
Bayes, Stochastic Gradient Descent classifier, linear discriminant analysis and K Near est Neighbor
to accomplish face recognition. An inference on sparse nature of locally adaptive regression kernels
was made from the experimentation. A rigorous study on discrepancies of the performance of
LARK descriptors is reported.

41
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Personnel Selection using Soft Set Theory

Leena Sonia and Pragati Jaina
a
School of Mathematics, DAVV, Indore
Email: [email protected]

Abstract. Soft Set theory was formerly proposed as a general mathematical tool for dealing with
uncertainty; not clearly defined objects. Selection of a suitable person for a leading position in
any organization is very important strategic decision for a Personnel Department. The role of the
Personnel Department is to screen a number of received applications and select an ideal candidate
among those who are applying for the position of a particular post. Different job needs different
ability and criteria to measure. Thus, it needs a suitable and flexible method to evaluate the
performance of each candidate accordingly. The paper presents the notion of Fuzzy Soft Set and
its Application in Human Resource Management.

Key words: Soft Set, Fuzzy Soft Set, Optimization Technique, and HR Management.

Numerical Solutions of Periodic Burgers Equation using Collocation of

Cubic B-spline Finite Elements
Amit Tripathia and R. C. Mittalb
a
Department of Applied Science and Humanities Government Engineering College, Atarra,
Banda Uttar Pradesh-210201
b
Indian Institute of Technology Roorkee, Roorkee
Email: [email protected]

Abstract. In this work, we present a novel method to obtain numerical solutions of periodic
Burgers equation. Our approach is based on collocation of cubic B-spline basis functions over the
finite elements. First, we discretized the given equation in space and then we use collocation of
cubic B-spline basis functions for spatial variable and its derivatives, which produces a system of
first order ordinary differential equations. This system of differential equations has been solved
by SSP-RK54 scheme. The efficiency of the proposed method has been confirmed with numerical
experiments, which shows that the obtained results are acceptable and are in good agreement with
earlier studies. Ease of implementation and very small size of computational work are two major
advantages of this approach. The computational complexity of the method is linear in size of
partition. Moreover, this method provides approximate solutions not only at the grid points but
also at any point in the solution range.

42
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Solution of Space-Time Fractional Order Partial Differential Equations:

Modified Adomian Decomposition Method
Vidya Nilesh Bhadgaonkara and D.B.Dhaigudeb
a
Vidnyan Mahavidyalaya Malkapur, Dist Buldhana
b
Dr. Babasaheb Ambedkar Marathwada University, Aurangabad
Email: [email protected]

Abstract. The purpose of this paper is to suggest a new modied Adomian decomposition method
for finding a solutions of non homogeneous space-time fractional order partial differential equations
by using a fractional Taylor’s series expansion.T his method is more effective than the general
Adomian decomposition method. Keywords: nonlinear nonhomogeneous space-time fractional
order advection equation ; New modified Adomian decomposition method; Mittage-Leffler function;
Fractional Taylor’s series expansion.
Vibrational Spectra of Ozone (O3) Using Lie Algebraic Method
J. Vijayasekhara
a
Department of Mathematics, School of Technology, GITAM (Deemed to be University),
Hyderabad
Email: [email protected]

Abstract. In this paper, vibrational spectra (infrared and Raman) of Silane (SiH4) is calculated
using theoretical method (Vibron model). The determined vibrational energy levels by Vibron
model are compared with experimental results. The obtained results from our method are consis-
tent with experimental results.
Performance Analysis of a Modified Twisted Savonius Hydrokinetic
Turbine
Sourish Singhaa and R. P. Sainia
a
Indian Institute of Technology Roorkee , Roorkee, Uttarakhand, India
Email: [email protected]

Abstract. Extraction of hydrokinetic energy from in stream flow can be an effective, sustainable
and environmental friendly rep lacement of the conventional energy resources. Savonius hydroki-
netic turbine is a vertical axis turbine which has a poor efficiency and its operating range under
various loading conditions is very less. The coefficient of performance value reported under different
investigations on Savonius hydrokinetic turbine for increasing its efficiency are highly fluctuated in
respect of tip speed ratios (TSR). In order to enhance the efficiency of the turbine and to get less
fluctuations in coefficient of power output for smooth operation under higher range of tip speed
ratios, this paper aims to model a modified design of Savoni-us turbine with twisted blades and
analyse its performance on various operating conditions. For CFD analysis commercial unsteady
Reynolds-Averaged Navier-Stokes (URANS) solver in conjunction w ith SST k- ? turbulence
model have been used. Due to complex geometry 3D transient simulations are done to find o ut
the average torque. Coefficient of torque and coefficient of power are analysed and discussed. The
result of this study concluded as very less deviation in power coefficient up to tip speed ratio 1.4
which indicates the increment of operating range and beyond that TSR value it falls down rapidly.
The maximum power coefficient obtained 0.30 corresponding to a TSR value 1.4 for input water
velocity of 2 m/sec.

43
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Statistical Analysis of Cancer Data

C.P. Sri Chidambarama , R. Mohan Kumara and S. Sreenatha
a
Coimbatore Institute of Technology
Email: [email protected]

Abstract. Cancer is a group of diseases involving abnormal cell growth with the potential to
invade or spread to other parts of the body. It is the uncontrolled growth of abnormal cells in the
body. It develops when the bodys normal control mechanism stops working. Old cells do not die
and instead grow out of control, forming new, abnormal cells. These extra cells may form a mass of
tissue, called a tumor. In cancerous tumors or malignant tumors, the cells have lost the ability to
stop growing. In other words, they have gone rogue and will not stop dividing. To better explain,
nearly every cell in the body is able to grow and divide to make new cells, to a certain extent. This
is important for all living organisms. When cells go rogue, however, there is a problem with the
DNA. When mutations, which are changes in th e DNA sequence, occur, they cause the cells to
forget how to stop dividing. After some time, the mass of cells becomes a tumor. The tumor can
either be malignant or it can be benign, which means it is not cancerous. The cancer data analysed
in this project work comprises of state wise distribution of cancer cases, sex wise distribution of
cancer cases in each state and age wise distribution of cancer cases from 2011 to 2016. The aim of
this work is to find whether there is a significant difference between the mean individuals affected
by cancer in all the States, in both the sex group (Male and Female) and also in all the age groups
by Analysis of Variance. If there is a significant difference, we find which state, sex group and
age group is prominently affected by cancer. Finally, the cancer cases are predicted by fitting a
suitable trend equation and the results are presented in a graphical form and interpreted by means
of Statistical analysis.

Digital Mammogram Analysis using Bendlets

N. Syamalaa , T. Krishna Chaitanyaa and B. Chandra Mohana
a
Dept. of ECE, Bapatla Engineering College, Bapatla, Andhra Pradesh
Email: [email protected]

Abstract. Worldwide, breast cancer in woman is the most alarming disease causing deaths in
recent times. Proper screening a nd diagnosis at an early stage, can save the lives of thousands
of woman. Detecting and diagnosis of abnormal tissue is the best way to cure breast cancer.
Analysis through biopsy is a costly and invasive and often not preferred by many patients. An
alternative to this is utilization of digital mammograms. Mammogram screening is non-invasive.
However, it requires domain expertise for accurate analysis. Hence, the classification efficiency
has to be improved. In this paper, an algorithm is proposed for digital mammogram analysis
using Bendlets with improved accuracy. Bendlets are extended versions of shearlets and are able
to capture the curvature of the discontinuities in the images accurately. Benchmark database
MIAS is used for assessing the performance of the algorithm. Several classifiers are explored for
the analysis. Performance of the conventional KNN classifier is superior compared to the other
classifiers used in this work. Simple statistical features are used in the work. Both 10 fold and 20
fold cross validation is done. Precision, sensitivity, specificity and accuracy are considered as the
metrics for the analysis. Performance of the proposed algorithm is compared with the state of the
art existing algorithms.

44
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

The MHD based numerical analysis of unsteady natural convection flo

w and heat transfer in a square cavity with rotating heated cylinder u
sing OpenFOAM
Ranjit J. Singhaa and Trushar B. Gohila
a
Department of Mechanical Engineering, Visvesvaraya National Institute of Technology, Nagpur
Email: [email protected]

Abstract. The present analysis is concentrated on the numerical simulation of two-dimensional

unsteady free convection flow a nd heat transfer arisen by the temperature alteration produced by
the rotating heated cylinder inside the cold square cavity at the fixed Rayleigh number of Ra =
10 7 and Prandtl number of Pr = 0.71. The transport properties of the fluid are expected to be
incompressible, viscous, non-magnetic and electrically conducting in nature. The magnetic field is
applied in the vertical (parallel to the y-axis) direction with intensity in terms of Hartmann number
of Ha = 0 an d 50. The cylinder is rotated in the clockwise (CW) and counter-clockwise (CCW)
direction about its axis with an angular velocity of ? = 100 (rad/s). The magnetohydrodynamics
based buoyancy- driven natural convection flow solver i s developed in the open source CFD tool
OpenFOAM. The direct numerical simulation is executed to explore the influence of the Lorentz
force on the quashing out the undulation generated by the unsteady flow. The fluid motion in the
existence of the magnetic field produces the electric current (j) which further interrelates with the
magnetic field to generate the Lorentz force. The Lorentz force forcefully clashes with the buoyancy
force and hence the undulation in the flow of fluid is reduced. The study of the development of the
Lorentz force and its directional effect on the flow a nd heat transfer is reported. Furthermore, the
relevant transport quantity as Nusselt number also calculated to observe the effect of the existence
of the Lorentz force in the fluid domain. It is observed that the undulation in the flow, heat
transfer, drag and lift coefficient is greatly influenced by the Lorentz force.

CFD simulation of a trickle bed with foaming liquid

Rahul Omara
a
Department of Cement Technology, AKS University, Satna (M.P.)-485001
Email: [email protected]

Abstract. The objective of the present study is to numerically investigate the effect of liquid and
gas velocities of foaming liquid solutions with varying liquid surface tension on pressure drop and
dynamic liquid saturation in a trickle bed. CFD simulations were carried out on the system by
using two-dimensional Eulerian two phase porous media model. Simulation results indicate higher
pressure drop with higher superficial liquid and gas velocities, lower surface tension of the liquid,
the dynamic liquid saturation increase with increase in liquid velocity, decrease with increase in
gas velocity and decrease with decrease in surface tension. These general hydrodynamic behaviors
are quantified with foaming liquids as the liquid phase. Pressure drop is constant along the radial
direction and it increases with increase in bed height. There is linear variation of pressure drop
along the column height which indicates that there is uniform flow of the fluid in the in the column.

45
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Numerical study on the efficiency of magnetophoresis in human

vasculature like conditions
Rekha Sa and Sarbari Bhattacharyaa
a
Department of Physics, Bangalore University, Bengaluru - 560056, In dia.
Email: [email protected]

Abstract. Magnetic Nanoparticles (MNPs) have been proposed as therapeutic agents for treat-
ment of medical conditions like cancer using them either as drug carriers to specific locations in the
body or as agents for localised heating. For these treatments to be effective a substantial quantity
of MNPs have to reach the target site from the site of injection. This may not be the case if
blood flow alone is relied on to transport the MNPs as losses to branching blood vessels leading
to transport in directions away from the target site can occur. We have carried out numerical
simulations on the flow of a dilute concentration of magnetic particles (MPs), where inter particle
interactions can be ignored, in channels that mimic human vasculature like conditions using Open
Field Operation And Manipulation (OpenFOAM). We have treated blood as a Newtonian fluid
and have considered a laminar flow. The size of the channel considered and the direction of low
(towards target site) would mimic a situation in an arteriole. We have considered the situation
where the main channel gives rise to a single daughter channel with a similar cross section at dif-
ferent branching angles to the direction of flow. While we find that a substantial number of MPs
are lost into the daughter channel when the branching angle is less than 45◦ with moderate loss
for branching angles greater than 45◦ , an optimum external magnetic field gradient significantly
improves the fraction of particles in the desired direction with no loss for branching angles greater
than 90◦ .
Mhd Flow Over A Stretching / Shrinking Sheet With Variable Heat
Flux
Kalpna Sharmaa
a
Department of Mathematics, Manipal University Jaipur, Jaipur, Rajasthan
Email: [email protected]

Abstract. In this paper, we have analyzed MHD forced convective boundary layer flow and heat
transfer through a stretching/shrinking sheet with variable heat flux in the presence of heat source.
Using a similarity transformation, the governing equations are transformed into a system of non-
linear ordinary differential equations, which are then solved by homotopy analysis method (HAM).
The effects of various parameters, namely the magnetic field parameter, heat source parameter,
suction parameter, Prandtl number, stretching/shrinking sheet parameter and heat flux parameter
on the velocity and temperature are discussed graphically. Moreover, a comparative study between
the previously published and the present study is made.

Numerical Study on Fluid Dynamics of Collapsible Channels

Vedant Dhruva , Ujwal Mishraa and Ranjith Maniyeria
a
Department of Mechanical Engineering, National Institute of Technology Karnataka, Surathkal,
Mangalore-575025,
Email: [email protected]

46
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Abstract. The fluid flow in collapsible channels or tubes is an interesting problem with several
physiological applications; for ex ample, blood flow in veins, air flow in lungs and wheezing. In this
paper we present a fluid-structure interaction based model for a single phase fluid flow through a
microchannel containing two elastic walls. A two-dimensional model is developed and simulations
have been performed on COMSOL Multiphysics 5.2a. The deforming geometry is analyzed using
moving mesh. The flow field and deformation of the elastic walls for different boundary loads and
inlet flow conditions are presented and discussed.

Analysis of the flow of two immiscible electrically conducting fluid

through parallel plates.
Ruchika Mehtaa
a
Department of Mathematics and Statistics Manipal University Jaipur, Jaipur
Email: [email protected]

Abstract. An unsteady flow of two immiscible fluids through horizontal parallel plates with time-
dependent oscillatory wall transpiration velocity is investigated. Both the fluids are assumed to
be electrically conducting , and region II is filled with porous medium.Separate solutions for each
fluid are obtained. The expressions for velocity and temperature distributions are obtained using
regular perturbation technique, discussed numerically and shown through graphs. The expressions
of skin-friction and Nusselt number are derived, discussed numerically and their numerical values
for various values of physical parameters are shown through graphs.

Comparative study of boundary conditions in LBM for incompressible l

aminar flow
Alankar Agarwala and Akshay Prakashb
a
VDepartment of Mechanical Engineering, Indian Institute of Technology Jodhpur, Rajasthan,
India - 342037, b Indian Institute of Technology Kharagpur, West Bengal, India
Email: [email protected]

Abstract. In this paper, we conduct a comparative study amongst different boundary condi-
tions with two-dimensional single-relaxatation time lattice Boltzmann method (SRT-LBM), for
incompressible laminar flow. Three types of boundary condition are considered for the simulation:
including full-way bounce-back, half-way bounce-back, and modified bounceback for the imple-
mentation of no-slip boundary condition on the wall with, pressure (density) boundary condition
proposed by Zuo and He [1] for inlet/oulet. The benchmark fluid flow problem of steady plane
Poiseuille’s flow with Reynolds number, Re = 75 is chosen . The numerical simulations are vali-
dated with the analytical solution, and grid convergence test are performed to compare accuracy
of different boundary conditions.

47
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Finite Element Analysis of MHD Blood Flow in Stenosed Coronary

Artery with the Suspension of Nanoparticles
Ankita Dubeya and B. Vasua
a
Department of Mathematics, MNNIT Allahabad, Allahabad
Email: [email protected]

Abstract. The numerical study presents for two dimensional mathematical modelling and com-
putational simulation of flow and heat transfer analysis of non-Newtonian nanofluid in stenosed
coronary artery in the presence of magnetic field. We considered the blood flow model based on
second grade fluid flow and heat transfer equation. The blood fluid viscosity is considered as a
function of temperature based on Vogels model. In this study the blood is considered as base fluid
with the suspension of nanoparticles. In order to complete our model, the variability in design and
size of stenosis is considered. The transformed conservation equations are solved numerically by
using finite element method in conjunction of variational approach and FreeFEM++. The results
show that increasing the thermophoresis parameter (Nt) and the Brownian motion parameter (Nb)
caused an increase in temperature values and nanoparticles concentration at the stenosis and a
decreases in nanoparticles concentration at the boundary of the vessel after the stenosis. These
properties changes in the domain by changing the shapes and designs of the stenosis in the domain.

Computational Analysis of Buoyant Convection with Top-Heated and

Bottom-Cooled Inclined Cylindrical Enclosures Using OpenFOAM
A. K. Raja , K. Dileepa and S. Jayaraja
a
National Institute of Technology Calicut
Email: [email protected]

Abstract. Gravity-assisted convection helps in the re-circulation of fluid within the cylindrical
enclosures when heated on their sides by passive solar radiation. This paper primarily aims at de-
veloping a transient, incompressible, laminar, threedimensional generalized numerical solver using
OpenFOAM which specifically deals with the heat transfer aspects of a non-pressurized ETC-Type
SWH system. A scaled model of the SWH system with an unstructured domain was discretized by
using Finite Volume Method (FVM) for the present analysis. The inclination angle of the tubes
and Prandtl number of the fluid was set to be 30◦ and 6.3, respectively. The effect of Re and
dimensionless temperature on non-dimensional variables such as Gr, Ra and Nu with respect to a
fixed aspect ratio (length to diameter) has been analyzed. The flow contours concluded that the
re-circulation is significant adjacent to the top wall while the approaching fluid cools the bottom
wall. The solver was validated with the works of Hayek et al. and Thant et al.

48
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A solution to enhance the efficiency of CFD analysis results for real

time usage.
Daniel Nayagam Ja , Shalini R Nairb and Jessy Roobyb
a
Structural Engineering, Hindustan Institute of Technology and Science, Chennai, b Department
of Civil Engineering, Hindustan Institute of Technology and Science, Chennai.
Email: [email protected]

Abstract. Importance of wind effects on buildings has grown with its height. Though wind itself
is invisible and innocuous, the effect it can have on tall structures can be visible and catastrophic.
Unless proper effort is taken for study in this area, proper mitigation techniques cannot be imple-
mented. Over the years researchers have developed Wind Tunnel as an experimental method to
replicate the response and behavior of building to wind effects. Though wind tunnel test is the
only source that can provide the required data, it requires skilled professionals for modelling, plac-
ing and reading the results. Moreover, the experimental procedure is uneconomical for multiple
tests. Any iterative research method requires more than two iterations. Hence this stands as great
hindrance to any aspiring researcher. With technological development computer-based simulation
software incorporated with computational fluid dynamics was created as an alternative method.
Though this method eliminates several disadvantages of wind tunnel it has not yet become source
of reliable data provider especially in the field of structural engineering. The validation of a method
can be established through a relation of a common point or a constant deviation. This paper aims
to validate CFD based softwares by establishing a relation between the experimental data and the
data obtained through simulation.

Mixed Convective Flow Past an Inclined Square Cylinder at Low

Reynolds Number
Md Reyaz Arifa and Nadeem Hasana
a
Zakir Husain College of Engg. and Tech., Aligarh Muslim University, Aligarh, UP, India
Email: [email protected]
Abstract. The present work aims to numerically investigate the mixed convective flow past an
inclined square cylinder with incidence angle equal to 45o and to determine the critical Richardson
number ( Ric ) at which the vortex shedding (VS) i s suppressed. The flow parameters are Re
= 100, Pr = 0.71, ? = 0o, ? = 45o. Also the aerodynamics parameters such as coefficient of
drag, CD and coefficient of lift, CL are studied with Richardson number varying in the range
of 0.2 ? Ri ? 1.0. The Boussinesq approximation is used to model the governing equations.
The continuity, momentum and energy equations with suitable boundary conditions are solved
using FVM solver commonly known as ANSYS FLUENT. The convective terms in the momentum
and energy equations are discretized using QUICK scheme. To determine the VS suppression
phenomenon, a numerical simulations are performed for a set of Ri close to the bifurcation point.
The simulations are performed in the Ri range 0.2 ? Ri ? 1.0. By observing the instantaneous
streamline patterns an d the coefficient of lift plot with time it is observed that VS suppression
occurs at a critical Richardson number, Ric = 0.78. The variation of aerodynamics parameters
such as coefficient of Drag, CD and coefficient of lift, CL with Ri are also studied in the mixed
convection regimes. It is found that CD increases with the increase in Ri and no lift is developed
for any Ri.

49
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Investigating Behaviour of Multi-Clutch Plate Frictional Materials

using ANSYS
Nishant Dhengrea , Ashish Mograb and M k Pradhanc
a
NMIMS b MPSTME Shirpur, a MANIT Bhopal
Email: [email protected]
Abstract. Clutches plays an essential part in automobile transmission in modern high-speed
vehicle. Clutches are devices which are highly required as it used to control comfort, economy and
sportiness Clutches are mainly use to engage and disengage the driving shaft to a driven shaft,
so that the driven shaft may be stopped or start without effecting the driving shaft. It can be
used in any machine where power is being transmitted and subjected to variable loading like in
automobile. Various studies are being carried out in analysis and increasing the torque carrying
capacity of clutch in automobile. But when we use this clutch in any customize vehicle it dosed not
transmit same efficiency as it was designed for. In this study, an attempt is made to investigate
the behavior of clutch plate under different loading. Currently used material of friction disc used
in analysis are Cork and Kevlar fiber. During engagement and dis engagement of clutch plate
deteriorate and get erode at every cycle. The aim of the research is to study the effect of axial
force offered by the spring for engagement using theories based on uniform pressure and uniform
wear. Clutch plate of CT-100 was designed in solidworks and for its FEA analysis ANSYS software
was used.

Case Wise Study Of Heterogeneity On The Traversal Characterstics Of

Torsional Surface Waves.
Chandni Kumaria and Santimoy Kundua
a
Department of Applied Mathematics, IIT (ISM) Dhanbad
Email: [email protected]

Abstract. Within the framework of a heterogeneous anisotropic layer laid over a pre-stressed half-
space, propagation of torsional surface wave is investigated. We have considered two different cases
of heterogeneity function in the anisotropic layer, namely exponential and trigonometric hyperbolic
functions. We have also obtained two different dispersion equations for both cases of heterogeneities
on the propagation of torsional wave in the said model. The effects of heterogeneities and pre-stress
on the phase velocity of torsional wave have been shown graphically for their different numerical
values.

50
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Exact traveling wave solutions and bifurcation analysis for time

fractional dual power Zakharov-Kuznetsov-Burgers equation
Amiya Dasa and Niladri Ghosha
a
Kazi Nazrul University, Asansol, 713340, WB
Email: [email protected]

Abstract. In this paper, we introduce the time fractional dual power Zakharov-Kuznetsov-Burgers
equation in the sense of modified Riemann-Liouville derivative. We apply the (G’/G)-expansion
method in adherence of fractional complex transformation and obtain various form of exact trav-
eling wave solutions including solitary wave, kink-type wave, breaking wave and periodic wave
solutions of the equation. Next we investigate the dynamical behavior, bifurcations and phase
portrait analysis of the exact traveling wave solutions of the system in presence and absence of
damping effect. Moreover, we demonstrate the exceptional features of the traveling wave solutions
and phase portraits of dynamical system with the help of interesting figures.

Impact of Irregularity on the Propagation of SH wave in a Slightly

Compressible Composite Structure
Sayantan Guhaa , Mriganka Shekhar Chakia and Abhishek Kumar Singha
a
Dept. of Applied Mathematics, IIT (ISM), Dhanbad
Email: [email protected]

Abstract. The present paper enunciates the propagation behavior of SH wave in a composite
structure comprised of a slightly compressible layer overlying a slightly compressible half-space
with rectangular shaped irregularity and parabolic sh aped irregularity at the common interface in
two distinct cases. Dispersion relation has been deduced analytically in closed form. The impact
of various parameters, viz. wave number, rectangular shaped irregularity and parabolic shaped
irregularity associated at the common interface, on the phase velocity of SH wave has been studied
for the composite structure. Numerical computation and graphical illustration has been done in
order to analyse the impact of irregularity parameters (rectangular and parabolic) more profoundly
for slightly compressible varga material which serves as one of the major highlight of the present
study.

51
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Dynamics of Polarised Shear waves propagating in Al/BaTiO3

composite structure
Sonal Nirwala and Sanjeev Anand Sahua
a
Dept. of Applied Mathematics, IIT (ISM), Dhanbad
Email: [email protected]

Abstract. The present article investigates the transference of shear waves in a plate surface
bonded with a Functionally Grade d Piezoelectric Material (FGPM) layer. Phase velocity of
horizontally polarised shear wave (SH wave) is calculated for both, electrically open and short
case. The method of separation of variables is used to solve the problem. Some particular cases
are derived. Significant influence of parameters (like material gradient and layer width) on phase
velocity of shear wave has been highlighted. Numerical example is provided in support of obtained
analytical solution. Findings may be applicable in the structural health monitoring and ultrasonic
inspection techniques.

Wave interaction with a floating circular porous elastic plate

Harekrushna Beheraa and Vinay Kumar Guptaa
a
SRM Institute of Science and Technology, Chennai, Tamil Nadu
Email: [email protected]

Abstract. In this paper, wave interaction with a circular porous elastic floating plate in finite
water depth is investigated under the linearized theory of water waves. The eigenfunction expansion
method–along with the Darcy’s law for wave past porous structure under the assumption that the
flexible structures are having fine pores and are of homogeneous structural rigidity–are exploited
to obtain an analytic solution of the problem. The plate deflection and flow distribution around
the plate are computed, and the results are analyzed to study the effects of various wave and
structural para meters. Moreover, the results obtained by considering porosity are compared with
the results in case of impermeable flexible plate. The study reveals that due to the porosity of
the plate, major part of the wave can be dissipated. Th e findings of this work are expected to
make substantial contributions to the safety/serviceability assessment and to the designs of coastal
structures for creating tranquility zones in order to protect coastal areas and marine facilities.

52
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Investigation of the Diabetic Level in Patients Through Computational

Fluid Dynamic
Pardeep Bishnoia and M.K. Sinhaa
a
NIT Jamshedpur
Email: [email protected]

Abstract. This article presents a new method to measure diabetes in the patient. Conventionally,
the experiments are perform ed by measuring and analyzing the constituents of urine and their
diabetic level which are mainly based on tradition al clinical methodologies. In the current work,
the diabetic level of a patient can be analyzed by using a drop dynamics process. This process
includes a urine sample of the diabetic patient being used and a drop test created from a capillary
nozzle with the help of the FLUENT tool. Through this mechanism, after setting the thread
length of a normal persons urine droplet as the lowest limit, a comparative analysis is drawn and
the variation of thread length in the diabetic patients urine sample is determined. When the
diabetic level increases in the body, there is an increase in the viscous nature of the urine. For the
range of glucose concentration level 182 mg/dl -342 mg/dl, the variation recorded in the thread
length is 14 %, which is a measurable quantity. A variation of thread length is plotted for different
urine samples with detachment time of the drop. For the entire analysis, thread length is the major
factor to compare the diabetic level between a normal person and a diabetic persons urine. The
results of the analysis and actual values of blood sugar glucose in the samples were found to be in
good correlation with each other.

CT scan based modeling and finite element analysis of femur bone for
different loading conditions
Hasin Khana and P M Padolea
a
Department of Mechanical Engineering, VNIT, Nagpur
[email protected]

Abstract. Femur is the longest, heaviest and strongest bone of the human body. It has been
observed that sudden accident and continuous fatigue loading is the main reason for the femur
bone failure. While many studies have analyzed the bone considering static loading, only a few
have simulated the bone using dynamic loading of the gait cycle. In this paper, a 3D model of
the femur was generated from bone CT (Computed Tomography) scan of a 60-year-old female
using MIMICS Materialize software. The femur was assumed to consist of only cortical bone since
its modulus of elasticity is dozen to hundred times that of cancellous bone. Gait cycles were
studied extensively and the femur was loaded dynamically using the time-dependent load data
obtained from Orthoload considering 75 kg as the body weight. FE Analysis was done in ANSYS
Workbench for loading conditions such as level walking, stair climbing, and other human routine
activities. Maximum normal stress and fatigue life corresponding to the above loading conditions
were obtained and compared with results from previous studies.

53
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A Coupled Field Finite element analysis of a Dielectric Elastomer-based

Buckling Actuator
Chinmay Jivania , Atul Kumar Sharmaa and M. M. Joglekara
a
Department of Mechanical and Industrial Engineering Indian Institute of Technology Roorkee,
247667 Uttarakhand, India
Email: [email protected]

Abstract. This paper presents finite element analysis of circular disk type dielectric elastomer
actuator. This type of DEAs are circumferentially fixed so that when actuated, they tend to buckle
out in their thickness direction. The nonlinear, three dimensional finite element model is developed
for DEA subjected to large deformation. The behavior of the DE material is characterized by
employing isotropic hyperelastic neo-Hookean material model. This along with the Maxwell stress
generated due to applied voltage, are used to formulate this electromechanically coupled model.
The circular geometry of the DEA is meshed by using eight noded hexahedron elements. The
model is solved numerically by using Newton-Raphson method for different boundary and pre-
stretch conditions. The numerical results for a freely expandable DEA are validated by using the
analytical lumped parameter model. Then a parametric analysis is done for buckling deformation
of DEA disk with different pre-stretch values. The electromechanical instability parameters are
obtained for different pre-stretch values to define the maximum achievable bulging deformation.

Design and Analysis of Chassis of Hydraulic Test Trolley

Roshan Gawandea and M.S. Kotambkara
a
VNIT Nagpur
Email: [email protected]

Abstract. The Hydraulic Test Trolley (HTT) is one the ground support equipment of the aircraft
and helicopters. It is used to supply hydraulic power to aircraft under test and remote hydraulic
consoles. Each unit consists of a reservoir, heat exchanger, pump/motor groups, filters, manifolds,
and electrical control box mounted on a rigid frame. Controls are manually operated and labelled
for ease of use. Maintenance items such as filters and suction strainers are located for easy access.
The hydraulic circuit control components (relief valve, check valve, filter, and differential pressure
switch) is mounted and interconnected through a remote mount filter manifold system. The present
work focuses mainly on the design and analysis of the chassis of the Hydraulic Test Trolley (HTT).
The material of structure is ASTM A710 steel with 450 MPa of yield strength and 515 MPa of
tensile strength. This paper aims to model, simulate and perform the stress analysis of the ladder
chassis of C-section beam. The scope of this study is to get the information about the maximum
stress and deformation of C-section beams used for chassis design, which will be used for further
design improvement. A 3-D model is created by using a modelling software i.e. CATIA P3 V5R17
and then transported to the ANSYS WORKBENCH for the analysis. The results obtained by
ANSYS are then compared to analytical calculations.

54
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Structural analysis of hydraulic aircraft jack using solid mechanics and

mechanical APDL
Sandeep Rathora and M.S. Kotambkara
a
Dept. of Mechanical Engineering, VNIT, Nagpur
Email: [email protected]

Abstract. The aircraft jack is a portable, self-contained, hydraulically operated jack. These jacks
are used for raising the wing, nose, or tail of aircraft. It is a specific type of jack used to lift a
heavy load of aircraft. When it is used in sufficient numbers and at the required jacking points,
the jack can lift the complete aircraft off the deck. The aim of this paper is to determine the load
and stresses on each structural member of the Jack when the vertical load is acting on the top of
the telescopic cylinder assembly having semi spherical adoptor on its top by using the concept of
solid mechanics and simulation on mechanical APDL. Using mechanical APDL, we got the forces
on each member of the jack and bending moment on all fixed nodes. By using APDL results
and using the concept of solid mechanics we got all the stresses on the structural members of
the jack which is checked by using maximum shear stress theory, whether the jack is safe or not.
The project study is limited to the design of structural member only. It also gives best possible
arrangement of structural members of jack and nature of joints. Using all different load condition,
selection of safest design will be done. Based on the factor of safety obtained cost optimization can
be done. This study is done on existing jack which is one of the products of a group of industries
in the region dealing with the ground support equipment and t heir maintenance. The principle
parameters for the jack are maximum lifting capacity of 15 tons and stroke of 1030 mm.

Kinematic Analysis of Five Bar Mechanism under Tolerance

Aditya Rathorea and H.P Jawalea
a
VNIT Nagpur
Email: [email protected]

Abstract. Geometrical deviations can influence the functioning of manipulator systems in posi-
tion. Geometrical deviations due to manufacturing discrepancies are limited by tolerances which
are defined during the design process.Inappropriate dimensional tolerances can cause output po-
sition to fall below the expected desired position resulting in diminished performance.Since the
position of end-effector produced will be different for each set of dimensional value.This paper
presents the effect of tolerance on positional deviation of end-effector due tolerance in individual
link length..And fin ding that link which is responsible for maximum error for given input rotation
angle.Since it is necessary to predict the position of end-effector , because there exist the proba-
bility that the output position will fall below a desired posit ion.This paper shows the positional
analysis of five bar mechanism by algebraic method by taking tolerance into account.This paper
shows the comparison of workspace of mechanism when taking tolerance in individual links and
when there is no tolerance in individual link.The result obtained shows the variation of error on
varying the input rotation angle. The result obtained through Matlab is validated through CATIA
and analytically.

55
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A nonlinear mathematical study of Tropical Cyclone Eye: A review

Shashi Kanta
a
India Meteorological Department, Bhubaneswar, India 751020
Email: [email protected]

Abstract. Tropical cyclones (TCs) are one of the common and most devastating weather activities
around the globe. The high speed winds associated with the tropical cyclone usually generate
ocean waves and tides resulting into heavy rain, which generally results in the destruction or much
damage. Both human and economy is are the victims of tropical cyclones. This disastrous nature
of cyclones attracts policymakers and governments to observe this phenomenon seriously like other
disastrous events. calamitous properties are the main concern of scientists and meteorologists to
observe them and accordingly improve the early prediction of the moving trajectories. The motive
of this article is to review those mathematical theories associated with the TC eye. This is a
nonlinear mathematical study.

A Finite Difference Approach to Study the Three-Dimensional Rotating

Flow of Carbon Nanotubes Past a Stretching Surface
Rajan Kumara
a
Indian Institute of Technology (Indian School of Mines) Dhanbad
Email: [email protected]

Abstract. This research article explains the effects of volume fraction of CNTs and magnetohy-
drodynamics on the flow over a stretching sheet in rotating frame. We analysed the behavior of base
fluid such as water with both the kind of CNTs i.e multi-wall carbon nanotubes (MWCNTs) and
single-wall carbon nanotubes (SWCNTs). The nanofluid transport mechanism is obtained follow-
ing Xue model. The mathematical model is developed by order of magnitude analysis. Governing
partial differential equations are first converted into a system of ordinary differential equations
making the us e of suitable similarity transformation and then solved by developing finite differ-
ence scheme using Newtons linearization technique, which helps us to deal with the nonlinear terms
in a very smooth way. Numerical simulation for the fluid velocity and fluid temperature along with
skin friction coefficient and Nusselt number are carried out for a range of values of regulatory flow
parameters that characterize the physics of the flow by graphically. An outstanding correlation
is perceived in tabular form between the two numerical approaches of finite difference scheme i.e
second order method and fourth order method for both momentum and energy equations.

56
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Facial Analysis using Jacobians and Gradient Boosting

A Vinaya , Abhijay Guptaa , Aprameya Bharadwaja , Arvind Srinivasana , K N B
Murthya and Vinayak Kamatha , S Natarajana
a
PES University
Email: [email protected]

Abstract. Security and identity have become one of the primary concerns of the people in this
digital world. Person authentication and identification is transforming the way these services are
provided. Earlier it was mainly achieved through passwords and patterns but with significant
advancements in face recognition technologies, it has been exploited in providing authentication in
smart phones and computers. Face Recognition(FR) extends its use in applications like face tagging
in social media, surveillance system at theaters, airports and so on. The proposed mathematical
model employs linear algebra and mathematical simulations for the task of person identification.
Kernel singular value decomposition is used to denoise the facial image which is then passed to
a feature detector and descriptor based on nonlinear diffusion filtering. The obtained descriptors
are quantized into a vector using an encoding model called VLAD which uses k-means++ for
clustering. Further, classification is done using Gradient boosting decision trees. The pipeline
proposed aims at reducing the average computational power and also enhancing the efficiency of
the system. The proposed system has been tested on the three benchmark datasets namely Face
95, Face 96 and Grimace .

Experimental investigation on optimization of concentration of

nanofluid(MgO) for enhancing heat transfer coefficient and comparison
of heat transfer coefficient with distilled water by using smooth
aluminium pipe domain.
Sagar Bhingarea , Shubham Lidhadea and Vinayak Jadhavb
a
Discipline of Mechanical Engineering,VJTI,Mumbai, 400019,india
b
Mechanical Deparment of Rajarambapu Institute of Technology, India
Email: [email protected]

Abstract. The development of high performance thermal systems to compacting the size of heat
exchanging devices is basic need of industries. Heat transfer coefficient is not material property and
can be improved by increasing the thermal properties of fluid like specific heat ,thermal conductivity
,viscosity,density .High thermal conductivity of solids can be used to increase the thermal conduc-
tivity of a fluid by adding small solid particles to that fluid[2]. In this work different concentrations
of nanofluid are prepared[1] and experimental trials on (0.005%,0.01%,0.05%,0.1%,0.5%,1%, 1.5%)
vol./vol. concentrations of nanofluid and distilled water are conducted. Analysis has been car-
ried out on this trials to calculate heat transfer coefficient and obtained results are represented
by graphs. From graphical analysis for turbulent flow maximum heat transfer rate among all the
concentrations of nanofluids is achieved for 0.01% concent ration fig.(2). Nanofluid gives better
performance for enhancing heat transfer coefficient in turbulent zone than laminar zone with re-
spect to water. Maximum percentage increase in heat transfer coefficient for 0.01% concentration
at 7 lpm is 61.88%. with respect to water.

57
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Unsteady Natural Convection in a Fluid-Saturated Porous Enclosure

with a pair of l-Shaped Diathermal Partition
Chordiya Jayesh Subhasha and R. V. Sharmaa
a
NIT Jamshedpur
Email: [email protected]

Abstract. Unsteady natural convection in a fluid-saturated square porous enclosure with a pair
of L- shaped diathermal partition located symmetrically at bottom-left and top-right corner of
enclosure is numerically analysed in this paper. The fundamental objective is devising a way to
reduce the heat transfer rate across a differentially heated porous enclosure. The fact that the
bottom-left and top-right corner of enclosure shows higher temperature gradient is the reason of
choosing L-shaped partition such that these particular regions are obstructed and isolated during
the fluid flow; which in turn causes significant reduction in Nusselt number. The governing equa-
tions are based on Darcys model and numerically dealt with Successive Accelerated Replacement
and explicit scheme. Streamlines and isotherms are observed for partition ratio 0, 0.1 and 0.2,
where partition ratio is the ratio of distance between horizontal or vertical partition to horizon-
tal or vertical center of enclosure to the total length of the enclosure; while Nusselt number is
estimated to manifest the heat transfer rate. Modified Rayleigh number is taken as 100, 500 and
1000.

Haar Wavelet: History and Its Applications

Mahendra Kumar Jenaa and Kshama Sagar Sahua
a
Department of Mathematics, Veer Surendra Sai University of Technology Burla,
Sambalpur-768018, Odisha, India
Email: [email protected], [email protected]

Abstract. In this paper, we have given a brief history of the Haar wavelet. Later the operational
matrix which is obtained from Haar wavelet is used to find the numerical solutions of some differen-
tial equations. The solutions thus obtain from operational matrix method is compared with exact
solution as well as solution from Runge-Kutta method and Modified Euler method is presented.

Convective flow in an anisotropic porous wavy channel

V. K. Narlaa
a
Department of Mathematics, GITAM School of Technology, Hyderabad
Email: [email protected]

Abstract. Heat transfer analysis for fully developed convective flow in a wavy channel filled with
an anisotropic porous material subject to constant heat flux has been investigated. Flow inside the
anisotropic porous bed is assumed to be governe d by the Brinkman model. The main objective of
the work is to study the effect of anisotropy parameter on the hydrodynamic and convective heat
transfer inside the wavy channel.

58
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Theoretical study of dynamics of the micro-nanomotor

Sunny Kumara
a
Department of Chemical Engineering, MaNIT, Bhopal
Email: [email protected]

Abstract. Tiny machines are the next generation micro-nanomotor [1-3] which can perform many
tasks such as in medical application, digital microfluidic application, sensor in chemical industries.
There are so many materials which can be use d for the fabrication of micro-nanomotor and can be
of different shapes, and sizes. These tiny machines can perform dynamics such as rotation, linear
motion, zigzag motion, and splitting. The micro-nano motors dynamics required so me external
forces which are generated by magnetic field, chemical gradient, electric field, light sensitive, surface
tension and ultrasonic. These external forces are required in tiny magnitude such as micro to Pico
newton. Here, we summarize the used mathematical equation for external forces which can help
to calculate the theoretical velocity and drag force. We have also explained dynamics (velocity) of
tiny machines in term of dimensionless number and comp are the dynamics of different size and
shape machines (in terms of body length per second). References: 1. S. Kumar et al., Graphene
based multifunctional superbots, Carbon, 2015, 89, 31- 40. 2. S. Kumar et al., Magnetic field
Induced push-pull motility of liquibots, RSC Advance, 2016, 6, 107049-56. 3. J. Wang et al.,
Synthetic micro/nanomotors in drug delivery, Nanoscale, 2014, 6, 10486-10494.

Heat Transfer Analysis of CNT-Nanofluid between Two Rotating Plates

in the Presence of Viscous Dissipation Effect
Amit Kumara , Ramayan Singha and Rajat Tripathib
a
NIT Jamshedpur, b Rewa Engineering College, Rewa
Email: [email protected]

Abstract. In this research study, an investigation of three-dimensional (3D) CNT based nanofluid
flow through a horizontal rotating channel, is carried out. The Mathematical model of the problem
is developed and is presented in the form of a system of nonlinear partial differential equations.
Similarity transformation technique is used to convert these partial differential equations into or-
dinary differential equations. The transformed equations are then solved numerically by bvp4c
routine of MATLAB. Computations for the nanofluid velocity and nanofluid temperature along
with skin fraction coefficient and Nusselt number are, carried out for relevant flow parameters. A
comparative analysis of single-wall carbon nanotubes (SWCNTs) and multiwall carbon nanotubes
(MWCNTs) on velocity and temperature distribution is carried out for different values of nanopar-
ticle volume fraction. It is observed that increasing values of suction/injection parameter leads to
decrease in the temperature profile.

59
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Numerical investigation of natural convection heat transfer of Nano

fluids in a vertical annulus
Suhail Ahmad Khana , M Altamush Siddiquia and Shahid Husaina
a
Heat Transfer Research Lab, Mech. Engg. Deptt., AMU
Email: [email protected]

Abstract. The objective of present study is to carried out numerical investigation of natural
convection heat transfer of Al2O3 water Nano fluids in partially heated vertical annulus of high as-
pect ratio. Numerical investigation has been carried out by using the computational fluid dynamics
(CFD) solver Ansys Fluent 14.0. The convective terms of the momentum and energy equations
are discretized by the third order MUSCL and second order upwind scheme respectively. The pres-
sure?velocity decoupling is done by using semi implicit method for the pressure-linked equations
(SIMPLE). Effect of Nanofluids on natural convection as function of volume fraction of nanoparti-
cles is investigated and results show enhancement of heat transfer coefficient and deterioration of
Nusselt number with increase of volume fraction of nanoparticles.

Thermal radiation and viscous-Ohmic dissipation influence on Ag-water

nanofluid flow over a flat porous plate with suction/injection
Ashish Mishraa and Manoj Kumara
a
G. B. Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, India, 263145
Email: [email protected]

Abstract. The present article focus on numerical study of the impact of thermal radiation,
viscous-Ohmic dissipation and heat generation/absorption on MHD flow of Ag-water nanofluid
over a flat porous plate in the existence of suction/injection. The fundamental PDEs are con-
verted into dimensionless ODEs employing appropriate transformation and explained numerically
by four-fifth order Runge-Kutta-Fehlberg method through shooting algorithm. Two shape depend-
ing models of nanoparticles with unlike formulae of thermal conductivity and dynamic viscosity
are considered. The velocity and temperature profiles are plotted for assorted values of thermal
radiation parameter, magnetic field, heat generation/absorption parameter, Eckert number, suc-
tion/injection parameter and permeability parameter. Moreover, the impact of these pertinent
factors on skin friction coefficient and local Nusselt number are created in a tabular way. It has
been inspected that as enhancing the value of suction/injection parameter, the heat transfer coef-
ficient boosts while it decreases as an enhancement in the porosity parameter. Also, the upshots
portray that viscous dissipation and thermal radiation have a tendency to accelerate the boundary
layer thickness.

60
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Analysis of Cooling Effect by Spray Impingement of Nanofluids

Bikash Pattanayaka , Abhishek Munda and Jayakumar J Sa
a
Amrita School of Engineering, Amritapuri Campus, Amrita Vishwa Vidyapeetham, Kerala
Email: [email protected]

Abstract. Spray impingement technology acts as an effective cooling media for high heat flux
application. It plays a vital role i n industries like heat treatment of steel plates, nuclear re-
actor, coal gasification etc. From research papers, it has be en concluded that nanofluids have
enhanced thermal properties than base fluids. So in this experiment Al2O3, TiO2, and Al2O3
doped TiO2nanofluids at different volume concentrations were used. The nanoparticles were syn-
thesized by using HEBM (High Energy Ball Milling) technique at 300 rpm with the ball to powder
ratio of 10:1. These nanoparticles were characterized by using XRD, SEM, and TEM and were
found to be in the range of 30nm-50nm. These nano particles were suspended in water in addition
to polyvinyl alcohol of 3% for making it stable. The densities, viscosity,thermal conductivity, the
specific heat of the nanofluid were calculated using different models. The heat transfer study was
carried out on an electrically preheated iron plate of dimensions 100mm x 100mm x 8mm at dif-
ferent temperatures of 200oC, 150oC, 100oC. The cooling rate, the effect of air pressure on cone
angle and cooling rate were analyzed. The uncertainties in the experiment were due to temperature
fluctuations and location of the thermocouple. Hence from the experiment, an optimized nanofluid
for effective cooling was determined.

Experimental and Numerical Validation of Hybrid Double Pipe Heat

Exchanger using Nanofluid
Abhishek Munda , Bikash Pattanayaka and Jayakumar J Sa
a
Amrita School of Engineering, Amritapuri Campus, Amrita Vishwa Vidyapeetham, Kerala
Email: [email protected]

Abstract. Nanofluids are the colloidal suspension (solid-liquid mixture) of nanoparticles in a

base fluid; basically, it is a two-ph ase fluid. In comparison with base fluid, nanofluid has more
effective thermal conductivity and enhanced thermal properties. Many researchers have studied
about hybrid or composite nanofluid and they also have good heat transfer characteristics. In this
experiment, CuO, ZnO, and ZnO doped CuO nanoparticles were prepared by using high energy ball
milling technique. These nanoparticles were characterized by using XRD, SEM, and TEM and were
found to be in the range of 30nm-50nm. Polyvinyl alcohol of 3g stable nanofluid. The stability was
observed for 32 hours which was good to conduct an experiment. The density, viscosity, thermal
conductivity, the specific heat of the nanofluid was calculated using various models. The overall
heat transfer coefficient, logarithmic mean temperature difference, friction factor and effectiveness
of the hybrid double pipe heat exchanger using the nanofluid were calculated using NTU method.
At same input conditions the data obtained using ANSYS (FLUENT) 18.2 were compared with
the experiment. An optimized volume concentration of the nanofluid was found out to be used as
an effective cooling fluid in the hybrid heat exchanger.

61
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Analysis and Comparative Study of Various Performance Measures of

M/G/1 and M/G/S Queuing Models
Shejal Guptaa and Pradeep K. Joshia
a
Shri Vaishnav Institute of Management, Indore
Email: [email protected]

Abstract. Queuing theory is a mathematical tool, which deals with the phenomenon of waiting
lines. Keeping customers waiting for too long, affects various parameters of queuing theory. In
this paper, optimization technique has been applied to calculate the performance measures like
waiting time, traffic intensity, waiting cost and proportion of time the server is idle for single server
M/G/1 model as well as for multiple servers M/G/S model. Also, we analysis the comparative
study of various performance measures of single server M/G/1 model and multiple servers M/G/S
queuing model. This investigation deals with M/G/1 and M/G/S Models in which arrival follows
Poisson distribution and service time follows general distribution and Queue discipline is FCFS.
The result of analysis showed that as we increase the numb er of servers, the waiting time in the
queue decreases, total cost of the system decreases, traffic intensity also decreases whereas idle
time increases. The results are optimize and effective. The graphical representation of various
performance measures and their comparison are shown by MATLAB software.

CIT Communication Architecture for Improved Security and Data

Privacy in IOT Devices.
Dhileepan Thangamanimarana and S. Sanjay Kumara
a
Coimbatore Institute of Technology
[email protected]

Abstract. In this world of IOT devices growing exponentially comes with it challenges involved
it, which includes security and data privacy as they are deployed in fields of health care, home
automation etc. Our aim is to mitigate security issues in authentication of devices and privacy in
data transfer. The primary concern that comes to mind while dealing with these challenges are
limited amounts of storage, memory, processing capability and power as most of them operate with
batteries. These IOT devices operates in public networks which can lead to compromise in security
but building separate network for these devices involves cost for infrastructure. Thus, we propose
a communication architecture f or client server interaction which we call as CIT (Cash -In Transit)
architecture which involves sharing of encryption key to the server through Cellular network and
data transfer through the Internet which will be encrypted with the k ey mentioned above like
splitting the freight and locks key in cash-transits. Using analogue noise generator random number
is generated which is then used as the encryption key and gets shared to the server through Cellular
networ k followed by encrypted data through internet. The authentication of device is by a two-
factor authentication which involves one-time-password sent by server to device through Cellular
network which in turn sends back to server through Internet. This two-factor authentication will
be done for log on to critical applications and services on the serve r. This will improve security
by reducing probability of security compromise over network and data privacy issues.

62
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Modeling of Viscoelastic Fluid Flow Past a Non-Linearly Stretching

Surface with Convective Heat Transfer: OHAM Analysis
CArnab Bhattacharyyaa , Gauri Shanker Setha and Rajan Kumara
a
Indian Institute of Technology (Indian School of Mines) Dhanbad
Email: [email protected]

Abstract. This paper inspects the two-dimensional boundary layer hydromagnetic flow of an
incompressible, electrically conducting and thermally radiating viscoelastic fluid past a non-linearly
stretching surface in the presence of viscous and Joule dissipation. Similarity transformation is
used to convert the governing partial differential equations into a system of ordinary differential
equations. The trans-formed non-linear and coupled equations along with the boundary conditions
are then explored with the help of Optimal Homotopy Analysis Method (OHAM). Plots for velocity,
temperature a nd concentration fields are portrayed and their behavior are deliberated for several
sets of values of governing flow parameters. Numerical values of local skin friction coefficient, local
Nusselt number and local Sherwood number are presented in tabular form. It is concluded that
viscoelasticity reduces the velocity field whereas the convective heat transfer condition enhances
the fluid temperature.

Aggregation of LARK Vectors For Facial Image Classification

A Vinaya , Vinayaka Kamatha , Varun Ma , Nidheesha , S Natarajana and KNB
Murthya
a
PES University
Email: [email protected]

Abstract. Face recognition is prevailing to be a key aspect wherever there is a need for interaction
between humans and machines. This can be achieved by containing a set of sketches for all the
possible individuals and then cross-validating at necessary circumstances. We propose a mechanism
to fulfill this task which is centered around locally adaptive regression kernels. A comparative study
has been presented at encoding stages as well as at the classification stages of the pipeline. The
results are cautiously examined and analyzed to deduce the best mechanism out of the proposed
methodologies. All the ideologies have been tested for multiple iterations on benchmark datasets
like ORL, grimace an d faces95. The vectorized descriptors have been subjected to encoding using
slightly refined methods of feature aggregation and clustering to assist classifiers in imputing the
test subjects to their respective classes. The encoded vectors are classified using Gaussian Naive
Bayes, Stochastic Gradient Descent classifier, Linear discriminant analysis and K Nearest Neighbor
to accomplish face recognition.

63
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Robust and Efficient Face Recognition Using Minimum Barrier

Detection with XNORNet.
A Vinaya , Rahul R Bharadwaja , Gadicherla Sameera , Sagar V Belavadia , K.N.B.
Murthya and S Natarajana
a
PES University
Email: [email protected]

Abstract. Deep convolutional neural networks (CNNs) have been very successful in many com-
puter vision applications and are state-of-the-art. It has been very successful in the area of face
recognition where many models report almost perfect performance on long-standing and well known
benchmarks. A huge challenge being faced by CNNs is the fact that they require a huge amount
of memory and processing power to deliver these impressive results hindering their deployment on
smaller devices like smart-phones whose ubiquity is unparalleled. This paper critically explores a
CNN bas ed architecture in which convolutions are approximated by primarily binary operations
of XNOR and bitcounts making it much faster and more memory efficient than standard CNN
architectures. The goal of this paper is to evaluate the performance of the above architecture
for facial recognition and study the effect of self-normalizing properties induced by the activation
function Scaled Exponential Linear Unit (SELU) on the convergence of the network and compare it
with convergence performance of the network using the well known Rectified Linear Unit (RELU)
activation function. We have used two well known datasets Faces95 and Faces96 preprocessed
using Minimum Barrier Detection (MBD) for saliency detection to benchmark the performance of
the network. Our results indicate that the network using SELU converged faster than RELU for
the Faces95 dataset and both converged at similar rates for Faces96 dataset. The network with
SELU gave better accuracy than RELU for both the datasets.

Facial Image Classification using Rotation, Illumination, Scale and

Expression Invariant Dense Features and ENN
Vinay Aa , Ankura , Mayank Raja , Nikhil Ananda , K N Balasubramanya Murthya ,
Aniket Bharatia and S Natarajana
a
PES University, Bengaluru
Email: [email protected]

Abstract. Face Recognition is immensely proliferating as a research area in the paradigm of

Computer Vision as it provides an extensive choice of applications in surveillance and commercial
domains.This paper throws light upon the comparison of various dense feature descriptors(Dense
SURF, Dense SIFT, Dense ORB) with each other and also with their classical counterparts(SURF,
SIFT, ORB) using a novel technique for recognition. This proposed technique uses Laplacian Of
Gaussian filter for enhancement of the image. It applies various dense and classical feature descrip-
tors on the enhanced image and outputs a feature vector. In order to achieve high performance,
this feature vector is given to Fisher vector since Fisher Vector is a feature patch-aggregation
method. Finally, extended nearest neighbor Classifier is used for classification over the orthodox
k-nearest classifier. Experiments were carried out on three diverse datasets- ORL, Faces94, and
Grimace. On scrutinizing the results, Dense SIFT and Dense ORB were found to be preeminent
as measured by various performance metrics. Dense ORB furnished an accuracy of 98.15% on
Faces94 and Dense SIFT with an accuracy of 98.44 % on Grimace dataset

64
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A Copy Detection Technique based on Discrete Cosine Transformation

(DCT)
Shikha Srivastavaa
a
Department of Mathematics, GLA University, Mathura (UP), India
Email: [email protected]

Abstract. The increase in various kinds of image editing software resulted in the generation of
large number of duplicate copies of original images. In such a scenario there is an urgent need
of hashing based image authentication technique, which can easily identify duplicate copies of
original images while differentiating it from other images. In this paper, we have proposed an
image hashing technique based on Discrete Cosine Transform and Statistical features, which is
robust to number of image processing attacks including Gaussian low pass filtering. The input
image is initially pre-processed to remove any of kind of minor effects. Discrete cosine transform is
then applied to the pre-processed image to produce DCT coefficients from which statistical features
are extracted to form the image hash which is used for image identification. Different experiments
were conducted to compare the proposed technique with the state-of-the-art techniques. Receiver
Operating Characteristics (ROC) curve and Average mean value difference shows that the proposed
technique out performs many of the state-of-the-art techniques in terms of various performance
parameters. The result also indicates that proposed technique is robust to various kinds of image
manipulation operation and performs particularly well against Gaussian low-pass filtering.

Improved detection of anomalies via integrated multi-scale principal

component analysis: Application to a bench-mark quadruple tank
system.
K Ramakrishna Kinia and Muddu Madakayarub
a
Department of Instrumentation and Control engineering, Manipal Institute of Technology,
Manipal Academy of Higher Education, Manipal
b
Department of Chemical engineering, Manipal Institute of Technology, Manipal Academy of
Higher Education, Manipal
Email: [email protected]

Abstract. In process monitoring domain, continuous monitoring of industrial processes for diag-
nosing any anomalies is crucial since they are the deciding factors for the betterment of produc-
tion quality and safety. To tackle these anomalies or faults, Principal Component Analysis (PCA)
method has been successfully applied vowing to its capability to handle complex multiple variables.
However, presence of measurement noise in the data tends to degrade the quality of a developed
fault detection strategy. The measurement noise is distributed in time-frequency domain which
requires efficient feature extraction from time as well as frequency scale. Wavelet-based multi-scale
representation of data has shown to be a powerful data analysis tool since it can extract useful data
simultaneously from time as well as frequency domain. In this paper, integrated multi-scale PCA
(IMSPCA) using wavelets is used as modeling framework while T 2 and Q statistics are used as fault
detection indicators. The efficiency of the developed strategy is demonstrated on a bench-mark
quadruple tank system. The results from the simulation indicates that the integrated multi-scale
PCA is found to be better in detecting incipient faults in comparison to the conventional PCA
strategy.

65
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Numerical Study on the Pennes Bioheat Transfer Equation using finite

difference method
Akula Sri Charana , Ranjith Maniyeria
a
Department of Mechanical Engineering, National Institute of Technology Karnataka, Surathkal,
Srinivasnagar, 575025
Email: [email protected], [email protected]

Abstract. Many of the cancer treatment methods involve application of energy in some form and
cause ablation of tumor cells. Precise temperature reading need to be obtained to avoid damage
to normal tissues. Numerical methods provide a better tool to estimate temperature when such
readings are difficult to be obtained experimentally. The paper aims at analyzing the temperature
distribution in a skin with a tumor in it, when subjected to constant and sinusoidal heat flux. It
was found that Pennes Model Bioheat Transfer overestimates temperature as compared to that
estimated by Thermal Wave Model of Bioheat Transfer.

Numerical simulation of melting under the influence of magnetic field at

various Rayleigh numbers
Sumer Bharat Dirbudea and Vivek Kumar Mauryaa
a
Department of Mechanical Engineering, National Institute of Technology Delhi, Delhi-40
Email: [email protected]

Abstract. Melting (and solidification) occurs in metal castings of turbine blades, environmental
engineering, PCM-based thermal storage devices, etc. To improve the heat transfer efficiency,
understanding of the important flow processes is necessary. In this work, 2D unsteady numerical
simulations, using ANSYS-Fluent with the enthalpy-porosity formulation, are performed to study
an effect of a uniform transverse magnetic field on melting inside a cavity filled initially with a solid
gallium at various Rayleigh numbers. The horizontal walls of the cavity are considered insulated
and vertical walls are, respectively, considered hot and cold. The magnetic field is characterized
by the Hartmann number (Ha) and the results are shown for the Ha = 0, 30 and 50. It is observed
that the role of natural convection during the melting is in deducing the temperature distribution
and solid-liquid interface. The increased magnetic field found to have a suppressing effect on the
dominance of natural convection at all Rayleigh numbers. Keywords: Magneto-hydrodynamics,
melting and solidification, CFD simulations, natural convection, Hartmann number.

66
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Disciplined Investigations to Fix Number of Neuron in Hidden Layer of

Multilayer Perceptron Used to Forecast Length of Brick for Manually
Driven Rectangular Brick Making Machine
Arati R. Lendea , Pawan A. Chandaka and J. P. Modakb
a
PES’S Modern College of Engineering, Shivajinagar, Pune
b
D.Sc. and Dean (R and D), Priyadarshini College of Engineering, Nagpur
Email: [email protected]

Abstract. The manually driven brick making machine is developed by Dr. J. P. Modak and their
associates. Lot of modifications and investigation in this arena has been distinguished in last three
decades through numerous research papers. The aim and objectives of these investigations was to
develop a precise and accurate model to predict the length of brick in respect to several in-put
parameters of the machine system. The present research paper puts light on the experiments
carried out to fix the number of neuron in hidden layer of a feed forward neural network. The
performance of the ANN model is evaluated for every alteration in size of neuron of the network
through performance plot, regression plot and results are discussed in detail. The article finally
discovers the method to fix the optimal value of neuron size in hidden layer for manually driven
rectangular brick making machine and concludes the functioning characteristics of the ANN model.

Reliable Analysis of Riemann Problem in Magnetogasdynamics

Mithilesh Singha , Nidhi Handab and Shakuntla Sharmab
a
Applied Sciences Departments, REC, Sonbhadra.231206, b Department of Mathematics, Kanya
Gurukula Campus, Haridwar, Uttrakhand 249404
Email: [email protected]

Abstract. Abstract In this paper, we studied Riemann problem in an isentropic fluid flow in
the presence of magnetic field. Here we studied shock, rarefaction and stationary waves in the
magnetic field. We used Newton Raphson method to determine maximum and minimum values of
density function for different stationary curves.

Hand Talk System for Deaf And Dumb Person

Vikash kumara , Sanjeev Kumar Raghuwanshb and Praful Singha
a
Gla university mathura, b IIT Dhanbad
Email: [email protected]

Abstract. This project is a prototype of data glove which has ability to convert the movement of
figure into visual dis play and audible sound of predefined language. I have seen many people who
find it difficult to communicate with ot her due to language problem or due to problem in vocal
chord. In a diverse country like INDIA where number of mot her tongues spoken is 1652.People
form one region of country when traveled to other region find it difficult to convey their message,
which make it tedious to survive them in other region. So this data glove can convert the speaker
language into the listener language which make it possible to easy communication between different

67
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

language speaker. There is approximately 70 million people in world who are deaf and dumb.
These people use sign language to communicate with each other .But the person with no disability
in vocal chord does not learn sign language, this means that not everyone can understand sign
language. So it will be a tedious task for deaf or dumb person to communicate with other person
having no knowledge of sign language.This data glove can be used to replace code language which
is used during war. A particular command can be converted into a gesture and the commanding
officer have to make particular gesture using his hand and the command is displayed on the
command receiving personnel. Data glove consists of flex sensor . In this prototype I have used 3
flex sensor which produce 8 output result.
Simulation and Prediction of Cereals Rating Using Neural Networks
Model and Multiple Regressions
Monisha Ga , Vishnu Karthiklua and Keerthivasan Sa
a
MSc Data Science, Department of Computing, Coimbatore Institute of Technology
Email: [email protected]

Abstract. The objective of the present study is to compare the efficiency and consistency of the
regression models with that of neural networks models and in addition to identify which model
is well suited for the actual and the simulated dataset. The mean square error term is found to
identify how much the prediction gets varied from the actually obtained values. Same procedure is
carried in terms of fitting the neural networks model. The neural network model is fitted using the
build-in packages and functions in R and based on the results, inferences are made. Also the model
that is we ll suited for the actual dataset is found. The second part of the study is to simulate
the data. So, based on the column descriptions of the data, appropriate simulation is done. Then
the simulated dataset is normalised for processing. Then the above procedure is followed to the
simulated dataset to identify which model is well suited.
Comparative Study of Artificial Neural Network (ANN) and
Multilinear Regression Model for Predicting Thermal Characteristics of
Double Pa ss Solar Air Heater (DPSAH)
A. K. Raja , G. Kunala and S. Jayaraja
a
National Institute of Technology Calicut
Email: [email protected]

Abstract. This paper represents modelling and comparison of Artificial Neural Network (ANN)
and regression model for predicting the thermal characteristics of a counter-flow double pass so-
lar air heater (DPSAH) system. ANN and regression models were developed based on the ex-
periments conducted on air heater located at the Solar Energy Centre of NIT Calicut, India
(11.15o N, 75.49o E) during November/December, 2017 under the clear sky conditions. The exper-
iments were performed during 09:00 to 17:00 hrs on a regular basis. Based on the comparison
of both the models employed, it was concluded that the ANN model with Feed Forward Back
Propagation (FFBP) training algorithm with its LM (Lavenberg-Marquardt) learning variant and
LS (LogisticSigmoidal) transfer function having 10 hidden neurons shows closer results than mul-
tilinear regression model with the experimental findings. The values of coefficient of determination
(R2 ), correlation coefficient (R), Root Mean Squared Error (RMSE) and Mean Absolute Error
(MAE) for both the models affirms the findings.

68
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Disciplined Investigations to Fix Number of Neuron in Hidden Layer of

Multilayer Perceptron Used to Forecast Length of Brick for Manually D
riven Rectangular Brick Making Machine
Arati R. Lendea , Pawan A. Chandaka and J. P. Modakb
a
PES’S Modern College of Engineering, Shivajinagar, Pune, b Priyadarshini College of
Engineering, Nagpur
Email: [email protected]

Abstract. The manually driven brick making machine is developed by Dr. J. P. Modak and their
associates. Lot of modifications and investigation in this arena has been distinguished in last three
decades through numerous research papers. The aim and objectives of these investigations was to
develop a precise and accurate model to predict the length of brick in respect to several in-put
parameters of the machine system. The present research paper puts light on the experiments
carried out to fix the number of neuron in hidden layer of a feed forward neural network. The
performance of the ANN model is evaluated for every alteration in size of neuron of the network
through performance plot, regression plot and results are discussed in detail. The article finally
discovers the method to fix the optimal value of neuron size in hidden layer for manually driven
rectangular brick making machine and concludes the functioning characteristics of the ANN model.

Comparison of Regression And Neural Network Model For Indian

Carbon-di-oxide emission rate, Internet Uers And Literacy Rate
Karthick Kumar Ga , Sai Krishna Na and Gokul Ranjith Sa
a
Department of computing, Coimbatore Institute of Technology
Email:[email protected]

Abstract. The study is to compare the efficiency of the trend models with that of neural networks
models and to identify which model is well suited for the actual and the simulated data set using
the mean square error values. The trend is fitted for the data. The data set is divided into test
and train data set, then the regression is performed and the test data set is also predicted and
the mean square error term also found to identify how much the prediction gets vary from the
actually obtained values. Same procedure is carried in terms of fitting the neural networks model.
The neural network model is fitted using the build in functions in R and based on the results
inferences are made. Also the model that is well suited for the given data set is found. The second
part of the study is to simulate the data. So, understanding the nature of data and based on the
column descriptions of the data, appropriate simulation is done. Then the simulated data set is
normalized for processing. Then the above procedure is followed to the simulated dataset whether
which model gets well suited.

69
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Probabilistic intuitionistic fuzzy set based intuitionistic fuzzy time

series forecasting method
Krishna Kumar Guptaa and Sanjay Kumara
a
G. B. Pant University of Agriculture and Technology, Pantnagar, 263145, Uttarakhnd, India
Email: [email protected]

Abstract. Intuitionistic fuzzy set (IFS) can handle non-stochastic non-determinism that arises
due to single membership function for both membership and non-membership grade. Probabilistic
intuitionistic fuzzy sets (PIFS) can handle non-determinacy that arises due to both stochastic and
non-stochastic reasons. In the present study we propose PIFS based fuzzy time series forecasting
model to handle both types of non-determinism along with non-stochastic uncertainty in time
series forecasting. The proposed method describes issue of non-determinism that arises due to both
randomness and linguistic representation of time series data. An aggregation operator to aggregate
the PIFS into intuitionistic fuzzy set (IFS) is also proposed in this study. The developed method
has been simulated using the financial time series data of TAIEX to confirm its outperformance
using root mean square error (RMSE).

Torque Minimization in Surgical Bone Drilling by Particle Swarm

Optimization
Ponnusamy Pandithevana , Varatharajan Prasannavenkadesana and Natarajan
Vinayaga Muruga Pandya
a
Department of Mechanical Engineering, IIITDM Kancheepuram
Email: [email protected]

Abstract. Machining operations such as grinding, milling, drilling, and sawing are used in or-
thopaedic surgery. Out of these, bone drilling is an essential process. Thrust force and torque
generated during bone drilling play a vital role in the success of the process being carried out. The
enormous quantity of force and torque will lead to osteonecrosis through the increase of temper-
ature around the drilling region. As a result, an irreversible state of thermal failure will occur if
the temperature increases above 47◦ C. Thus it affects the quality of the drilling and it will cause
adverse effects such as bone-implant failure, cracks in bone and nerve damage. So, it is very much
essential to have an optimal level of torque and thrust force. In this work, bone drilling was car-
ried out by considering the spindle rotational speed and feed rate at three levels. Optimal spindle
rotational speed (1269.8 rpm) and feed rate (51.73 mm/min) were identified through the particle
swarm optimization to obtain minimum torque during bone drilling.

70
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Discrete Time Minimum Variance Control of Satellite System

D Y Dubea , H G Patela
a
NIT Surat
[email protected]

Abstract. This paper is concerned with the types of stochastic disturbances affecting the potential
of the aerial system. The satellite system for continuous and discrete time domain is discussed. A
phase lead compensator completes the orientation successfully. Astroms single-input single-output
(SISO) model is implemented with using the minimum variance control strategy. The separation
principle then provides the optimal control law which curtails the cost function to a value as small
as possible. The satellite system is positioned for one quarter revolution with the co-ordination of
generalized minimum variance controller (GMVC) and standard generalized dual controller (GDC)
based on certainty equivalence assumption. The revolutions in radians are tracked as output of the
system for the input specified in degrees to the system. The controller proved useful in reducing the
overshoot and atmospheric disturbances which allows a stable motion even for larger time delays.

Effect of Environmental Pollutants on Rain due to Stakeholders

Nita H. Shaha , Moksha H. Satiaa , and Foram A. Thakkara
a
Department of Mathematics, Gujarat University, Navarangpura, Ahmedabad, Gujarat, India
[email protected]

Abstract. Human density, plant density and animal density give their contribution to change
global climate which is termed as global warming. This includes melting of ice, rising temperature
and sea level etc. The reason behind these events is mainly dependent on the level of environmental
pollutants. Increase in environmental pollutants also diminishes rain in the world. To compare
the impact of these three stakeholders on environmental pollutants and rain, we have developed
a dynamical system. The model is formulated using the system of non-linear differential equa-
tions. In the proposed paper, threshold in terms of environmental pollutants is computed for each
stakeholder. Stability is worked out. The model is validated through numerical simulation.

71
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Scrutinization of Discrepancy in Nozzle Throat Diameter for Supersonic

Aerospace Vehicles
Rahul Shuklaa , N Sankea
a
Osmania University
[email protected]

Abstract. A nozzle is a machine representing the regulation of drift, velocity, route, bulk, ap-
pearance and the thrust. All these exhaust from the same instrument i.e. nozzle. Nozzles come
in a variety of shapes and sizes depending on the mission of the airspace vehicle, this is very im-
portant in understanding the mission of the rocket. Supersonic nozzles have many applications in
the aerospace industry, including high-speed military and combat jets, rocket nozzles and missiles.
The traditional Convergent Divergent (CD) nozzle is a relatively simple devise used commonly
among most jets and rocket nozzles. Depending on the requirement of the thrust and its appli-
cations, the geometrical configuration of the nozzle varies. The shape of the nozzle is key to the
expansion process, and plays a vital part towards designing the flow to minimize the thrust lost.
It is crucial during engine performance the thrust generated by expanding the exhaust gas retains
its maximum potential. When gas is expanded through a CD nozzle supersonically, the flow un-
dergoes many forms of unique phenomena, including flow separation, unsteadiness, flow mixing,
turbulence, Shock Induced Boundary Layer (SIBL) separation and Mach Shock Diamonds. Some
of these phenomena lead to energy loss, thereby reducing the overall thrust generated by the noz-
zle. The thrust loss due to shock waves and boundary layer separations generated internally in
the nozzle region remains poorly understood, hence failed to reach maximum potential of an en-
gine. To reach its functional potential, the design of an exhaust nozzle serves two main purposes.
Firstly, to control engine exit pressure to reach to ideal design conditions, which could be achieved
by varying the Nozzle throat diameter. Secondly, to convert the potential energy of the expanding
gas into kinetic energy efficiently.
Image Segmentation of Breast Cancer Histopathalogy Images using
PSO based Clustering Technique
Vandana Katea and Pragya Shuklaa
a
DAVV
[email protected]

Abstract. Image segmentation has key influence in numerous medical imaging uses. An image
segmentation method that is based on the Particle Swarm Optimizer (PSO) is developed in this
paper for Breast Cancer Histopathalogy images of different magnification level(40X, 100X, 200X,
400X), thus simplifying image representation and making it meaningful and easier for future anal-
ysis. As lower the magnification level, the bigger is the field of view and as it can provide greater
detail, thus more time and care must be taken to use such images. Thus finding a segmentation
method that works equally well for all zoom levels of images is a big challenge. To illustrate the
better performance of the proposed method and its applicability on breast cancer images, results of
applications and tests are augmented which shows PSO image clustering approach using intra clus-
ter distance as an optimization function, performs better than state of-the-art techniques namely,
K-means and Genetic Algorithm(GA). The algorithms when given specified number of clusters,
finds the centroids, thus grouping similar image primitives.The influence of different values of PSO
control parameters on performance is also illustrated.

72
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Study of source influx on interdependent calcium and inositol- 1; 4; 5

triphosphate dynamics in a cardiac myocytecell
Nisha Singha and Neeru Adlakhaa
a
Sardar Vallabhbhai National Institute of Technology, Surat
[email protected]

Abstract. Intracellular calcium governs the most versatile and universal signaling mechanism
in living systems which includes contraction of the cardiac tissues, information processing in the
brain, release of digestive enzymes by the liver etc. Various investigators have been made on study
of calcium signaling in cardiac myocyte to understand its mechanisms. But most of existing inves-
tigators have mainly focused on study of calcium signaling in cardiac myocyte cell without paying
attention on interdependence of calcium signaling and inositol- 1; 4; 5 triphosphate signaling. In
this paper, we propose a mathematical model to understand the impact of source influx of both
calcium and inositol-1; 4; 5 triphosphate on these coupled signaling processes. This study suggests
that source influx plays a vital role in these signaling processes. Also, calcium and inositol-1; 4; 5
triphosphate shows a beautiful coordination with each other, which explains the role of inositol-1;
4; 5 triphosphate in calcium signaling in cardiac myocyte cell. Such studies will provide the better
understanding of various factors involved in calcium signaling in cardiac myocytes, which as a
result will be of great use to biomedical scientists for making protocols for various heart diseases
and their cure.

Spatiotemporal Dynamics of Immunogenic Tumors

Teekam Singha , Sandip Banerjeea
a
Indian Institute of Technology Roorkee
[email protected]

Abstract. A mathematical model, consisting of a system of two coupled reaction-diusion partial

differential equations describing the interaction between solid tumor and immune system (termed
as effector cells), is proposed here. The main focus is on the analysis of the spatiotemporal dynamics
of tumor cells and immune cells. The resulting system is analyzed and numerical simulations are
presented. Different types of spatial patterns with respect to different initial conditions, and time
are observed. Their analysis and mechanism of spatiotemporal pattern formation in immunogenic
tumor are studied. Spatiotemporal perturbation around non-spatial steady state beyond the linear
regime are obtained based on the analysis of higher order perturbation terms.

73
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

A Survey of Time Series Data Prediction Using Multi Agent

Architecture for Unification of Association Rules
Mohammed Ali Shaika
a
Chhatrapati Shahu Ji Maharaj University
[email protected]

Abstract. Incredible amount streams of data are regularly produced by dynamic environments in
a time series database that comprises of many arrays of values obtained in a precise period of time
and these are normally measured at various time stamps towards occurrence of events to mine all
the transactional data sets that consists of transactional behavior and association rules identify
occurrence of regular data set as a group of different item sets or their subsequences. Apriori
and FP Growth algorithms are used to mine all the yielded frequent patterns of item sets and
major expensive phase of FP Growth algorithm is candidate set generation and database mining
of association rules unification process [1].
Stability aspects of the first order Hayes equation with scalable
hysteresis
Santhosh Kumar Balijaa , Saurabh Biswasa , and Anindya Chatterjeea
a
IIT Kanpur
[email protected]

Abstract. Hysteresis is a strongly nonlinear rate-independent phenomenon wherein the slope of

the output curve changes discontinuously with the sign of the rate of change of input. Stability
analysis of hysteretic systems is challenging. Recently, a scalable hysteresis model has been reported
in [1]. Scalability in the hysteretic damping allows us to study interesting stability aspects of
otherwise-linear systems.
We take up for study the Hayess delay differential equation [2], with hysteretic damping of [1],
as below:

ż(t) = − pz(t) − qz(t − 1) − γθ(t)z(t), (1)

κ
θ̇(t) = {θa + β sgn(z(t)ż(t)) − θ(t)}|ż(t)|, (2)
|z(t)| + 

where p and q are real valued parameters; γ > 0 is a parameter governing the strength of hysteretic
dissipation; θ is an internal variable; κ, θa , and β are positive real parameters; and 0 <  << 1 is
a small regularizing parameter to smooth the singularity at z = 0. The internal variable θ in Eq.
(2) is the non-scalable part, while the hysteretic force f = θz is scalable.
Stability of the system of Eqs. (1) and (2) is examined on the parameter plane p − q. Following
[3], we use Galerkin projection to convert Eq. (1) to a system of 15 ODEs. We then use numerically
obtained Lyapunov-like exponents λ [4] for stability analysis. For a chosen point in the p − q
plane,λ > 0 implies the system is unstable. In computations, we assign instability to points where
exceeds 10−3 .
Some stability boundaries in the p − q plane contain periodic solutions. We compute periodic
solutions numerically by continuation and also approximately using harmonic balance and related
approximations.

74
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

There is an extended region in the p − q plane for which nonzero equilibria exist. To study the
stability of nonzero equilibrium solutions in those parameter regions, we let z = zc and θ = θc be
equilibrium states and consider small variations z(t) = zc +η(t) and θ(t) = θc (t)+ϕ(t). Substituting
these in Eqs. (1) and (2) gives

η̇(t) =q{η(t) − η(t − 1)} − γzc φ(t), (3)

κ
θ̇(t) = {θa + β sgn(zc η̇(t)) − θc − Φ(t)}|η̇(t)|. (4)
|zc + η(t)| + 

We consider zc > 0. Then |zc + η(t)| ≈ zc . The regularizing parameter  is no longer needed. Since
the system is scalable, we can take zc = 1. We obtain from Eqs. (3) and (4)

η̇(t) =q{η(t) − η(t − 1)} − γzφ(t), (5)

φ̇(t) =κ(θa − θc )|η̇(t)| + κβ η̇(t). (6)

Equations (5) and (6) is a system of pseudolinear delay differential equations. Stability regions
of the original (Eqs. (1) and (2)) and the pseudolinearized system (Eqs. (5) and (6)) overlap, but
are not identical. The reason can be explained in terms of differences in sets of initial conditions
used for computing Lyapunov-like exponents. On the stability boundary of the pseudolinear sys-
tem, there are solutions with linear drift in one state and periodicity in all other states.

References

[1] Biswas, S., Jana, P. and Chatterjee, A., (2016), ’Hysteretic damping in an elastic body with
frictional microcracks’, International Journal of Mechanical Sciences, 108, pp. 61-71.
[2] Hayes, N. D., (1950), Roots of the transcendental equation associated with a certain difference-
differential equation, Journal of the London Mathematical Society, 25, pp. 226-232.
[3] Wahi, P. and Chatterjee, A., (2005), Galerkin projections for delay differential equations, Jour-
nal of Dynamic Systems, Measurement, and Control, 127 (1), pp. 80-87.
[4] Marathe, A. and Chatterjee, A., (2006), Asymmetric Mathieu equations, Proceedings of the
Royal Society A, 462(2070), pp. 1643-1659.

75
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Predator-prey model with prey group defense and non-linear predator

harvesting
Rajat Kaushika and Sandip Banerjeea
a
Indian Institute of Technology Roorkee
[email protected]

Abstract. This paper is concerned with a predator-prey system with a prey group defense and
non-linear harvesting of the predator incorporating deterrence hypothesis for predators. Inclusion
of predator deterrence rate makes the modelling approach more practicable and exhibits significant
impact on the net predation. Taking all possible interactions into account, model equations are
formulated. In brief qualitative analysis, existence of interior equilibrium and stabilities of all equi-
librium points of the system are discussed to investigate the dynamical behavior of the ecosystem.
Hopf, transcritical and saddle-node bifurcations are illustrated for various parameters. Numerical
simulations are ecological justified and supportive of theoretical results.

Estimation of Annual Rice Production and Its Methane Emission in

Paddy Fields
C.S.Rahula , G.Prabhakarana and B.Hariharana
a
Coimbatore Institute of Technology, Department of Computing
[email protected]

Abstract. The project aim is to estimate the annual rice production of India for the successive
years and also to estimate the a mount of methane emission from the paddy fields. The outcome
of the project contributes as a quite useful analytic model for solving the environmental issues
regarding methane gas emission. India is one of the leading producers of rice crop, the basic
food crop. Methane is one of the important greenhouse gases accounting for 15% of the total
enhanced global warming. Scientifically, it is found that there arise methane emission in the paddy
fields during the process of cultivation. This is due to the waterlogged condition in rice fields
which creates an anoxic environment, which is conducive for methane production by anaerobic
methanogenic bacteria. Therefore Methane emission from rice fields is a matter of concern. The
amount of methane emission along the paddy fields in India is also estimated. The historical data is
collected and the simulation is processed using Monte Carlo Simulation method. These estimations
are calculated and visualized using R language.

A kinetic Simulation of the expanding plasma in axially nonuniform

magnetic field for space propulsion thruster.
Sneha Guptaa and Devendra Sharmaa
a
Institute for plasma research
[email protected]

Abstract. A fully parallized high performance computing module HTSM-01 is developed to

simulate expanding magnetized plasma flow in the cavity of helicon source based space propulsion
thruster. This grid-free thrusting mechanism is preferred for deep space propulsion application of

76
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

advance space-crafts for its durability and capacity of refuelling from the heterogenous planetory
atmosphere. The 2X-3V Lagrangian phase space simulation procedure achieves Poisson-Boltzmann
equilibrium in an expanding plasma to simulate the plasma acceleration through a magnetic nozzel
and estimate thrust with a range of fuel gas species and plasma parameters. The phase space dis-
tribution functions are computed for individual plasma species allowing to compute the boundary
fluxes and thrust produced by outflowing plasma sourced by a helicon device.

Design of Novel Efficient circuit switching for Quantum-dot Cellular

Automata
Mahaboob Subani Shaika and B Chandra Mohana
a
Bapatla Engineering College Islampet
[email protected]

Abstract. Circuit switched network is basic component to transfer the input signal within a
distributed communication network where several users exist . In traditional communication circuit
switching is paramount important for faster communication. this paper presents and evaluates a
modular design methodology to build high performance circuit switched network by using MG and
RMG and inverters. It is proved that proposed layout design is less complex and power efficient
with nearly accurate results and less latency than previous designs. Our proposed circuit switched
network is 40 percent less Complex in terms of cell count and nearly 70 percent less complex in
layout area then existing.

Role of subducting slab age and dip on the development of fore-arc

topography
Ritabrata Dasguptaa , Changyeol Leeb and Nibir Mandala
a
Department of Geological Sciences, Jadavpur University, Kolkata, India, b Faculty of Earth
Systems and Environmental Sciences, Chonnam Nati onal University, Gwangju, Republic of
Korea
[email protected]

Abstract. Lithospheric subduction zones are fascinating features in the field of geodynamics.
Dynamic topographic development associated with the subduction zones are diverse depending
upon some factors, among them slab dip and age are pivotal. Variation in slab dip results variation
on the down going motion of the subducting slab. However, slab dehydration reduces the viscosity
of the wedge so the extent of viscous reduction in the wedge works an inverse function of the
slab dip. On the other hand older slab posses, larger negative buoyancy compared to younger
slabs. To check the effect of these two factors on the development of subduction zone topography
we developed a computational fluid dynamic model of an oceanic subduction zone. In this finite
element model, we imply an Arbitrary Lagrangian Eulerian scheme to track the deformed mesh
configuration, which enables better calculation of the topographic development above the plates.
Results show that shallow dipping slabs (15?) develop deeper fore-arc depression compared to
high dipping (60?) slabs. The longer hydrated wedge for shallow dipping slabs facilitates the basin
deepening being more weaker rheologically. However, older slabs being more negatively buoyant
tends to roll back more compared to younger slabs. This roll backing phenomenon opens the fore
arc basin more hence forms a broader basin. We compared our experimental findings with the
Cretaceous Gyeongsang basin of Korean Peninsula. It resulted from the subduction of Izanagi

77
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

plate. Plate tectonic reconstructions compared with our results explains that constant younging
of the subducting Izanagi plate cease the opening of Gyeongsang basin.
Simulation of I-V Characteristic of Single Electron Transistor (SET)
With Scilab’s Model Based Design Codes
S. Kumara and H. Jeevakhana
a
National Institute of Technical Teachers Training and Research, Bhopal

Abstract. Scilab’s X-cos model-based simulation blocks have been used to simulate the equation
of drain current of single electron transistor (SET). The drain current single electron transistor
(SET) can be calculated by using the tunneling rate of an electron through the tunneling junction.
Effect of temperature on drain current has also been analyzed.
Reliability and sensitivity estimation of shuffle exchange network
Soni Bishta and S.B Singha
a
Department of Mathematics, Statistics and Computer Science, G.B. Pant University of
Agriculture and Technology, Pantnagar, India

Abstract. Multistage interconnection network is the network which provides fast and efficient
communication at the reasonable cost. A shuffle- exchange network is the special class of multistage
interconnection networks in which signal is transmitted through source nodes to the sink nodes with
some specific connection. It is the challenge to increase the reliability of multistage interconnection
networks in an efficient way using different methods. The present study explores the idea about
the reliability of the shuffle-exchange network in the three aspects viz. terminal, broadcast and
network reliability with the help of universal generating function in which all the components are
non-identically distributed. With the help of structure-function, here we also propose to investigate
the sensitivity with the help of Birnbaum measure so that the importance of different components
of the considered network can be examined.
On the convection phenomenon from an impermeable bluff body
sheathed by a porous layer
Rajendra S. Rajpoota and S. Dhinakarana
The Centre for Fluid Dynamics, Department of Mechanical Engineering, Indian Institute of
Technology Indore
a
[email protected], b [email protected]

Abstract. Hydrodynamic and thermal performance of porous media have interested researchers
for the past several decades as they facilitate exceptional flow control and heat transfer augmenta-
tion. The present study aims to investigate the effect of porous media on heat transfer enrichment
from a heated impermeable bluff body sheathed by a porous layer. Forced convective heat trans-
fer is analysed for different values of Darcy number (10−6 − 10−2 ), Reynolds number ( 40), and
non-dimensional porous layer thickness (0.1 − 0.3), employing D2Q9 lattice model of the lattice
Boltzmann method. The effects of considered parameters on coefficient of drag and lift, and local
and overall Nusselt numbers of the cylinder are analysed. Numerical results indicate the addition
of a porous layer can change the flow pattern, suppresses unsteadiness of the flow and enhances
the heat transfer rate significantly.

78
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

Investigating conformational flexibility of N-glycans via all-atom

molecular dynamics simulations
Rajarshi Roy and Parimal Kara
a
Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore,
Khandwa Road, Simrol, Indore - 453552, Madhya Pradesh, India
a
[email protected]

Abstract. Understanding three dimensional structures and conformational preferences of glycans

at a molecular level is fundamental to cognizance the critical events of biological processes. Multi-
ple conformations of a single glycan chain pose a more significant question to the specificity and the
binding mechanism to their specific receptor proteins. However, the flexible nature of the glycans
makes their characterization by experimental techniqes a challenging task. Higher flexibility in the
glycosidic linkage between the rigid monosaccharides has created a restriction to reveal structural
information using X-ray crystallography. All atom molecular dynamics (MD) simulation gives
atomic details of glycan conformations. Here, we simulate four different hybrid N-linked glycans of
virus glycoproteins, with varying length of mannose and complex chain. Multiple distinct confor-
mations are observed in all different sequences which are separated by a high energy barrier. We
also found different puckering states in some distinct monosaccharides throughout our simulation
timescales. The specific pattern of inter-residue and water-mediated hydrogen bonds are observed
which further govern the global conformation of the glycans. Identifying a relationship between
glycan chain length and branching, glycosidic linkages, puckering, and water-mediated interaction
facilitate engineering of biologically significant receptor, especially antibodies to deal with different
viruses and throw light on conformational landscapes of the whole glycome.

Investigating Effectiveness of Currently Used Protease Inhibitors

against HIV-1 Subtype D by Using Multiscale Simulations
MD Fulbabu Sk, Nisha Amarnath Jonniya and Parimal Kara
a
Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore,
Khandwa Road, Simrol, Indore - 453552, Madhya Pradesh, India
a
[email protected]

Abstract. Acquired immune deficiency syndrome (AIDS) is caused by the human immunodefi-
ciency virus type 1 (HIV-1) and type 2 (HIV-2). However, HIV-1 is more prevalent world wide
than HIV-2. The diversity in HIV-1 has given rise to multiple subtypes and recombinant strains.
A majority of research focussing on antiretroviral agents and inhibitor resistance mechanism has
been performed for subtype B viruses, although non-subtype B strains are responsible for 90%
of global infections in developing countries. One of such subtypes, D is spreading rapidly and
infecting an enormous population in North Africa and Middle East. HIV-1 protease plays a cru-
cial role in viral replication and maturation, which makes it one of the most attractive targets for
anti-retroviral therapy. However, drugs have been developed targeting subtype B. In this work, the
popular molecular mechanics PoissonBoltzmann surface area (MM-PBSA) method has been used
to investigate the effectiveness of the HIV-1/subtype B PR inhibitors darunavir (DRV), tipranavir

79
 International Conference on Mathematical Modelling and Scientific Computing
July 19-21, 2018, Indian Institute of Technology Indore

(TPV), atazanavir (ATV) and lopinavir (LPV) against subtype D PR. Our simulations suggest
that TPV is more effective compared to other inhibitors. An increase in favorable intermolecular
van der Waals and electrostatic interactions and a decrease in unfavorable polar solvation free
energy compared to other cases result in better binding affinity. In all cases, binding is favoured
by the intermolecular van der Waals and electrostatics interactions as well as nonpolar solvation
free energy while the configurational entropy and desolvation of polar groups always disfavor the
complex formation. Furthermore, we have noticed that the favourable intermolecular electrostatic
interaction is overcomplensated by the unfavorable polar solvation free energy. This suggests that
the binding is mainly driven by the van der Waals interactions. We have also performed dynamics
cross-correlations (DCC) analysis suggesting that the PR-TPV complexation tends to make the flap
regions tightly closed. We expect that our study will serve as a standardized platform optimizing
the design of future more potent drugs targeting protease of various subtypes of HIV-1.

80