Exploring efficient: numerical methods for differential equations

Numerical analysis is a way to do higher mathematical problems on a computer, a technique widely used by scientists and engineers to solve their problems. A major advantage of numerical analysis is that a numerical answer can be obtained even when a problem has no “analytical” solution. Results from numerical analysis are an approximation, which can be made as accurate as desired. The analysis of errors in numerical methods is a critically important part of the study of numerical analysis. Hence, we will see in this research that computation of the error is a must as it is a way to measure the efficiency of the numerical methods developed. Numerical methods require highly tedious and repetitive computations that can only be done using the computer. Hence in this research, it is shown that computer programs must be written for the implementation of numerical methods. In the early part of related research the computer language used was Fortran. Subsequently more and more computer programs used the C programming language. Additionally, now computations can also be carried out using softwares like MATLAB, MATHEMATICA and MAPLE. Many physical problems that arise from ordinary differential equations (ODEs) have magnitudes of eigenvalues which vary greatly, and such systems are commonly known as stiff systems. Stiff systems usually consist of a transient solution, that is, a solution which varies rapidly at the beginning of the integration. This phase is referred to as the transient phase and during this phase, accuracy rather than stability restricts the stepsize of the numerical methods used. Thus the generally the structure of the solutions suggests application of specific methods for non-stiff equations in the transient phase and specific methods for stiff equations during the steady-state phase in a manner whereby computational costs can be reduced. Consequently, in this research we developed embedded Runge-Kutta methods for solving stiff differential equations so that variable stepsize codes can be used in its implementation. We have also included intervalwise partitioning, whereby the system is considered as non-stiff first, and solved using the method with simple iterations, and once stiffness is detected, the system is solved using the same method, but with Newton iterations. By using variable stepsize code and intervalwise partitioning, we have been able to reduce the computational costs. With the aim of increasing the computational efficiency of the Runge-Kutta methods, we have also developed methods of higher order with less number of stages or function evaluations. The method used is an extension of the classical Runge-Kutta method and the approximation at the current point is based on the information at the current internal stage as well as the previous internal stage. This is the idea underlying the construction of Improved Runge-Kutta methods, so that the resulting method will give better accuracy. Usually higher order ordinary differential equations are solved by converting them into a system of first order ODEs and using numerical methods suitable for first order ODEs. However it is more efficient, in terms of accuracy, number of function evaluations as well as computational time, if the higher order ODEs can be solved directly (without being converted to a system of first order ODEs), using numerical methods. In this research we developed numerical methods, particularly Runge-Kutta type methods, which can directly solve special third order and fourth order ODEs. Special second order ODE is an ODE which does not depend on the first derivative. The solution from this type of ODE often exhibits a pronounced oscillatory character. It is well known that it is difficult to obtain accurate numerical results if the ODEs are oscillatory in nature. In order to address this problem a lot of research has been focused on developing methods which have high algebraic order, reduced phase-lag or dispersion and reduced dissipation. Phaselag is the angle between the true and approximate solution, while dissipation is the difference between the approximate solution and the standard cyclic solution. If a method has high algebraic order, high order of dispersion and dissipation, then the numerical solutions obtained will be very accurate. Hence in this research we have developed numerical methods, specifically hybrid methods which have all the above mentioned properties. If the solutions are oscillatory in nature, it means that the solutions will have components which are trigonometric functions, that is, sine and cosine functions. In order to get accurate numerical solutions we thus phase-fitted the methods using trigonometric functions. In this research, it is proven that trigonometrically-fitting the hybrid methods and applying them to solve oscillatory delay differential equations result in better numerical results. These are the highlights of my research journey, though a lot of work has also been done in developing numerical methods which are multistep in nature, for solving higher order ODEs, as well as implementation of methods developed for solving fuzzy differential equations and partial differential equations, which are not covered here

Analytical and Numerical Methods for Differential Equations and Applications

The book is a printed version of the Special issue Analytical and Numerical Methods for Differential Equations and Applications, published in Frontiers in Applied Mathematics and Statistic

The Effect of Malaysia General Election on Financial Network: An Evidence from Shariah-Compliant Stocks on Bursa Malaysia

Instead of focusing the volatility of the market, the market participants should consider on how the general election affects the correlation between the stocks during 14th general election Malaysia. The 14th general election of Malaysia was held on 9th May 2018. This event has a great impact towards the stocks listed on Bursa Malaysia. Thus, this study investigates the effect of 14th general election Malaysia towards the correlation between stock in Bursa Malaysia specifically the shariah-compliant stock. In addition, this paper examines the changes in terms of network topology for the duration, sixth months before and after the general election. The minimum spanning tree was used to visualize the correlation between the stocks. Also, the centrality measure, namely degree, closeness and betweenness were computed to identify if any changes of stocks that plays a crucial role in the network for the duration of before and after 14th general election Malaysia

Mathematical Modeling with Differential Equations in Physics, Chemistry, Biology, and Economics

This volume was conceived as a Special Issue of the MDPI journal Mathematics to illustrate and show relevant applications of differential equations in different fields, coherently with the latest trends in applied mathematics research. All the articles that were submitted for publication are valuable, interesting, and original. The readers will certainly appreciate the heterogeneity of the 10 papers included in this book and will discover how helpful all the kinds of differential equations are in a wide range of disciplines. We are confident that this book will be inspirational for young scholars as well

Numerical treatment for mathematical model of farming awareness in crop pest management

The most important factor for increasing crop production is pest and pathogen resistance, which has a major impact on global food security. Pest management also emphasizes the need for farming awareness. A high crop yield is ultimately achieved by protecting crops from pests and raising public awareness of the devastation caused by pests. In this research, we aim to investigate the intricate impacts of nonlinear delayed systems for managing crop pest management (CPM) supervised by Ordinary Differential Equations (ODEs). Our focus will be on highlighting the intricate and often unpredictable relationships that occur over time among crops, pests, strategies for rehabilitation, and environmental factors. The nonlinear delayed CPM model incorporated the four compartments: crop biomass density [B(t)], susceptible pest density [S(t)], infected pest density [I(t)], and population awareness level [A(t)]. The approximate solutions for the four compartments B(t), S(t), I(t), and A(t) are determined by the implementation of sundry scenarios generated with the variation in crop biomass growth rate, rate of pest attacks, pest natural death rate, disease associated death rate and memory loss of aware people, by means of exploiting the strength of the Adams (ADS) and explicit Runge-Kutta (ERK) numerical solvers. Comparative analysis of the designed approach is carried out for the dynamic impacts of the nonlinear delayed CPM model in terms of numerical outcomes and simulations based on sundry scenarios

MATLAB

This excellent book represents the final part of three-volumes regarding MATLAB-based applications in almost every branch of science. The book consists of 19 excellent, insightful articles and the readers will find the results very useful to their work. In particular, the book consists of three parts, the first one is devoted to mathematical methods in the applied sciences by using MATLAB, the second is devoted to MATLAB applications of general interest and the third one discusses MATLAB for educational purposes. This collection of high quality articles, refers to a large range of professional fields and can be used for science as well as for various educational purposes

Simulation and fabrication of micro magnetometer using flip-chip bonding technique

Magnetic field detection has been widely accepted in many applications such as military systems, outer space exploration and even in medical diagnosis and treatment. Low magnetic field detection is particularly important in tracking of magnetic markers in digestive tracks or blood vessels. The presence of magnetic fields’ strength and direction can be detected by a device known as magnetometer. A magnetometer that is durable, room temperature operation and having non-movable components is chooses for this project. Traditional magnetometer tends to be bulky that hinders its inclusion into micro-scaled environment. This concern has brought the magnetometer into the trend of device miniaturization. Miniaturized magnetometer is usually fabricated using conventional microfabrication method particularly surface micromachining in which micro structures are built level by level starting from the surface of substrates upwards until completion of final structure. Although the miniaturization of magnetometer has been widely researched and studied, the process however is not. Thus, the process governing the fabrication technique is studied in this paper. Conventional method of fabrication is known as surface micromachining. Besides time consuming, this method requires many consecutive steps in fabrication process and careful alignment of patterns on every layer which increase the complexity. Hence, studies are done to improve time consuming and reliability of the microfabrication process. The objective of this research includes designing micro scale magnetometer and complete device fabrication processes. A micro-scale search coil magnetometer of 15 windings with 600μm thickness of wire and 300μm distance between each wire has been designed. Keywords: Magnetometer, microfabrication, miniaturization, micro-scale

LHD vibrations analysis and numerical modeling during operations

Load-haul-dump vehicles (LHDs) are extensively used as primary loaders in mining operations. LHDs have proven to be vigorous, extremely productive and reliable in mining applications. They have a wide range of tramming capacities that have enabled them to become an essential component in the hard rock mining industry. Increased mining economic challenges and global competition means the mining industry has to maximize productivity by cutting down operating and capital costs. Also, improvements in safety standards have led to the demand for safer and efficient machines. LHD operators are at a high risk of whole-body vibrations (WBVs) exposure leading to musculoskeletal disorders (MSDs) over long exposure periods, and elevated lower back and neck injuries. Thus, there is a health and safety concern among LHD operators. Despite manufacturer’s emphasis on ergonomics, there is lack of adequate fundamental vibration models of large mining equipment accessible to the public. This research focused on developing valid analytical and numerical models for determining the vibration propagation in LHDs. Also, this research pioneered the development and analysis of comprehensive dynamic virtual models of LHDs with detailed vibration analysis of the operator-seat interface. The introduced LHD virtual prototype has a total of 24-DOF and captures the complex vibration mechanics of the LHD, with emphasis on vibrations reaching the operator seat-interface in the three dimensions (3D), x, y and z-directions. The RMS accelerations recorded at the operator-seat interface are 0.62 m/s² in the x-direction, 0.51 m/s² in the y-direction, and 1.01 m/s² in the z-direction which exceed the ISO-2631 comfort level --Abstract, page iii

A Multi-Agent Architecture for the Design of Hierarchical Interval Type-2 Beta Fuzzy System

This paper presents a new methodology for building and evolving hierarchical fuzzy systems. For the system design, a tree-based encoding method is adopted to hierarchically link low dimensional fuzzy systems. Such tree structural representation has by nature a flexible design offering more adjustable and modifiable structures. The proposed hierarchical structure employs a type-2 beta fuzzy system to cope with the faced uncertainties, and the resulting system is called the Hierarchical Interval Type-2 Beta Fuzzy System (HT2BFS). For the system optimization, two main tasks of structure learning and parameter tuning are applied. The structure learning phase aims to evolve and learn the structures of a population of HT2BFS in a multiobjective context taking into account the optimization of both the accuracy and the interpretability metrics. The parameter tuning phase is applied to refine and adjust the parameters of the system. To accomplish these two tasks in the most optimal and faster way, we further employ a multi-agent architecture to provide both a distributed and a cooperative management of the optimization tasks. Agents are divided into two different types based on their functions: a structure agent and a parameter agent. The main function of the structure agent is to perform a multi-objective evolutionary structure learning step by means of the Multi-Objective Immune Programming algorithm (MOIP). The parameter agents have the function of managing different hierarchical structures simultaneously to refine their parameters by means of the Hybrid Harmony Search algorithm (HHS). In this architecture, agents use cooperation and communication concepts to create high-performance HT2BFSs. The performance of the proposed system is evaluated by several comparisons with various state of art approaches on noise-free and noisy time series prediction data sets and regression problems. The results clearly demonstrate a great improvement in the accuracy rate, the convergence speed and the number of used rules as compared with other existing approaches