Abstract: This article presents an approach with regards to the Functional Testing of Space System (SS) that could be a space vehicle (spacecraft-S/C) and/or its equipment and components – S/C subsystems. This test should finalize the Space Qualification Tests (SQT) campaign. It could be considered as a generic test and used for a wide class of SS that, from the point of view of System Dynamics and Control Theory, may be described by the ordinary differential equations. The suggested methodology is based on using semi-natural experiment laboratory stand that does not require complicated, precise and expensive technological control-verification equipment. However, it allows for testing totally assembled system during Assembling, Integration and Testing (AIT) activities at the final phase of SQT, involving system hardware (HW) and software (SW). The test physically activates system input (sensors) and output (actuators) and requires recording their outputs in real time. The data are then inserted in a laboratory computer, where it is post-experiment processed by the MATLAB/Simulink Identification Toolbox. It allows for estimating the system dynamics in the form of estimation of its differential equation coefficients through the verification experimental test and comparing them with expected mathematical model, prematurely verified by mathematical simulation during the design process. Mathematical simulation results presented in the article show that this approach could be applicable and helpful in SQT practice. Further semi-natural experiments should specify detail requirements for the test laboratory equipment and test-procedures.
Abstract: SARS-CoV-2 virus is currently one of the most
infectious pathogens for humans. It started in China at the end of
2019 and now it is spread in all over the world. The origin and
diffusion of the SARS-CoV-2 epidemic, is analysed based on the
discussion of viral phylogeny theory. With the aim of understanding
the spread of infection in the affected countries, it is crucial to
modelize the spread of the virus and simulate its activity. In this
paper, the prediction of coronavirus outbreak is done by using SIR
model without vital dynamics, applying different numerical technique
solving ordinary differential equations (ODEs). We find out that ABM
and MRT methods perform better than other techniques and that the
activity of the virus will decrease in April but it never cease (for
some time the activity will remain low) and the next cycle will start
in the middle July 2020 for Norway and Denmark, and October 2020
for Sweden, and September for Finland.
Abstract: Parkinson's disease (PD) is a heterogeneous disorder with common age of onset, symptoms, and progression levels. In this paper we will solve analytically the PD model as a non-linear delay differential equation using the steps method. The step method transforms a system of delay differential equations (DDEs) into systems of ordinary differential equations (ODEs). On some numerical examples, the analytical solution will be difficult. So we will approximate the analytical solution using Picard method and Taylor method to ODEs.
Abstract: The problem of nonlinear oscillations of a two-layer liquid completely filling a limited volume is considered. Using two basic asymmetric harmonics excited in two mutually perpendicular planes, ordinary differential equations of nonlinear oscillations of the interface of a two-layer liquid are investigated. In this paper, hydrodynamic coefficients of linear and nonlinear problems in integral relations were determined. As a result, the instability regions of forced oscillations of a two-layered liquid in a cylindrical tank occurring in the plane of action of the disturbing force are constructed, as well as the dynamic instability regions of the parametric resonance for different ratios of densities of the upper and lower liquids depending on the amplitudes of liquids from the excitations frequencies. Steady-state regimes of fluid motion were found in the regions of dynamic instability of the initial oscillation form. The Bubnov-Galerkin method is used to construct instability regions for approximate solution of nonlinear differential equations.
Abstract: In this paper, an enhancement of the heat transfer using non-Newtonian nanofluids by magnetohydrodynamic (MHD) mixed convection along stretching sheets embedded in an isotropic porous medium is investigated. Case of the Maxwell nanofluids is studied using the two phase mathematical model of nanofluids and the Darcy model is applied for the porous medium. Important effects are taken into account, namely, non-linear thermal radiation, convective boundary conditions, electromagnetic force and presence of the heat source/sink. Suitable similarity transformations are used to convert the governing equations to a system of ordinary differential equations then it is solved numerically using a fourth order Runge-Kutta method with shooting technique. The main results of the study revealed that the velocity profiles are decreasing functions of the Darcy number, the Deborah number and the magnetic field parameter. Also, the increase in the non-linear radiation parameters causes an enhancement in the local Nusselt number.
Abstract: Three structure-dependent integration methods have been developed for solving equations of motion, which are second-order ordinary differential equations, for structural dynamics and earthquake engineering applications. Although they generally have the same numerical properties, such as explicit formulation, unconditional stability and second-order accuracy, a different performance is found in solving the free vibration response to either linear elastic or nonlinear systems with high frequency modes. The root cause of this different performance in the free vibration responses is analytically explored herein. As a result, it is verified that a weak instability is responsible for the different performance of the integration methods. In general, a weak instability will result in an inaccurate solution or even numerical instability in the free vibration responses of high frequency modes. As a result, a weak instability must be prohibited for time integration methods.
Abstract: In this work, a five step continuous method for the solution of third order ordinary differential equations was developed in block form using collocation and interpolation techniques of the shifted Legendre polynomial basis function. The method was found to be zero-stable, consistent and convergent. The application of the method in solving third order initial value problem of ordinary differential equations revealed that the method compared favorably with existing methods.
Abstract: Due to many applications and problems in the fields of plasma physics, geophysics, and other many topics, the interaction between the strain field and the magnetic field has to be considered. Adomian introduced the decomposition method for solving linear and nonlinear functional equations. This method leads to accurate, computable, approximately convergent solutions of linear and nonlinear partial and ordinary differential equations even the equations with variable coefficients. This paper is dealing with a mathematical model of generalized thermoelasticity of a half-space conducting medium. A magnetic field with constant intensity acts normal to the bounding plane has been assumed. Adomian’s decomposition method has been used to solve the model when the bounding plane is taken to be traction free and thermally loaded by harmonic heating. The numerical results for the temperature increment, the stress, the strain, the displacement, the induced magnetic, and the electric fields have been represented in figures. The magnetic field, the relaxation time, and the angular thermal load have significant effects on all the studied fields.
Abstract: The water-based bioconvection of a nanofluid
containing motile gyrotactic micro-organisms over nonlinear
inclined stretching sheet has been investigated. The governing
nonlinear boundary layer equations of the model are reduced to a
system of ordinary differential equations via Oberbeck-Boussinesq
approximation and similarity transformations. Further, the modified
set of equations with associated boundary conditions are solved using
Finite Element Method. The impact of various pertinent parameters
on the velocity, temperature, nanoparticles concentration, density of
motile micro-organisms profiles are obtained and analyzed in details.
The results show that with the increase in angle of inclination δ,
velocity decreases while temperature, nanoparticles concentration,
a density of motile micro-organisms increases. Additionally, the
skin friction coefficient, Nusselt number, Sherwood number, density
number are computed for various thermophysical parameters. It
is noticed that increasing Brownian motion and thermophoresis
parameter leads to an increase in temperature of fluid which results
in a reduction in Nusselt number. On the contrary, Sherwood number
rises with an increase in Brownian motion and thermophoresis
parameter. The findings have been validated by comparing the
results of special cases with existing studies.
Abstract: In this study, we have investigated the strict stability
of fuzzy differential systems and we compare the classical notion of
strict stability criteria of ordinary differential equations and the notion
of strict stability of fuzzy differential systems. In addition that, we
present definitions of stability and strict stability of fuzzy differential
equations and also we have some theorems and comparison results.
Strict Stability is a different stability definition and this stability
type can give us an information about the rate of decay of the
solutions. Lyapunov’s second method is a standard technique used
in the study of the qualitative behavior of fuzzy differential systems
along with a comparison result that allows the prediction of behavior
of a fuzzy differential system when the behavior of the null solution
of a fuzzy comparison system is known. This method is a usefull
for investigating strict stability of fuzzy systems. First of all, we
present definitions and necessary background material. Secondly, we
discuss and compare the differences between the classical notion
of stability and the recent notion of strict stability. And then, we
have a comparison result in which the stability properties of the null
solution of the comparison system imply the corresponding stability
properties of the fuzzy differential system. Consequently, we give
the strict stability results and a comparison theorem. We have used
Lyapunov second method and we have proved a comparison result
with scalar differential equations.
Abstract: Evaporative coolers has a minimum potential to reach the wet-bulb temperature of intake air which is not enough to handle a large cooling load; therefore, it is not a feasible option to overcome cooling requirement of a building. The invention of Maisotsenko (M) cycle has led evaporative cooling technology to reach the sub-wet-bulb temperature of the intake air; therefore, it brings an innovation in evaporative cooling techniques. In this work, we developed a mathematical model of the Maisotsenko based air cooler by applying energy and mass balance laws on different air channels. The governing ordinary differential equations are discretized and simulated on MATLAB. The temperature and the humidity plots are shown in the simulation results. A parametric study is conducted by varying working air inlet conditions (temperature and humidity), inlet air velocity, geometric parameters and water temperature. The influence of these aforementioned parameters on the cooling effectiveness of the HMX is reported. Results have shown that the effectiveness of the M-Cycle is increased by increasing the ambient temperature and decreasing absolute humidity. An air velocity of 0.5 m/sec and a channel height of 6-8mm is recommended.
Abstract: A numerical technique in a boundary-fitted curvilinear grid model is developed to simulate the extent of inland inundation along the coastal belts of Peninsular Malaysia and Southern Thailand due to 2004 Indian ocean tsunami. Tsunami propagation and run-up are also studied in this paper. The vertically integrated shallow water equations are solved by using the method of lines (MOL). For this purpose the boundary-fitted grids are generated along the coastal and island boundaries and the other open boundaries of the model domain. A transformation is used to the governing equations so that the transformed physical domain is converted into a rectangular one. The MOL technique is applied to the transformed shallow water equations and the boundary conditions so that the equations are converted into ordinary differential equations initial value problem. Finally the 4th order Runge-Kutta method is used to solve these ordinary differential equations. The moving boundary technique is applied instead of fixed sea side wall or fixed coastal boundary to ensure the movement of the coastal boundary. The extent of intrusion of water and associated tsunami propagation are simulated for the 2004 Indian Ocean tsunami along the west coast of Peninsular Malaysia and southern Thailand. The simulated results are compared with the results obtained from a finite difference model and the data available in the USGS website. All simulations show better approximation than earlier research and also show excellent agreement with the observed data.
Abstract: Numerical Methods is a course that can be conducted using workshops and group discussion. This study has been implemented on undergraduate students of level two at the Faculty of Engineering, International Islamic University Malaysia. The Numerical Method course has been delivered to two Sections 1 and 2 with 44 and 22 students in each section, respectively. Systematic steps have been followed to apply the student centered learning approach in teaching Numerical Method course. Initially, the instructor has chosen the topic which was Euler’s Method to solve Ordinary Differential Equations (ODE) to be learned. The students were then divided into groups with five members in each group. Initial instructions have been given to the group members to prepare their subtopics before meeting members from other groups to discuss the subtopics in an expert group inside the classroom. For the time assigned for the classroom discussion, the setting of the classroom was rearranged to accommodate the student centered learning approach. Teacher strength was by monitoring the process of learning inside and outside the class. The students have been assessed during the migrating to the expert groups, recording of a video explanation outside the classroom and during the final examination. Euler’s Method to solve the ODE was set as part of Question 3(b) in the final exam. It is observed that none of the students from both sections obtained a zero grade in Q3(b), compared to Q3(a) and Q3(c). Also, for Section 1(44 students), 29 students obtained the full mark of 7/7, while only 10 obtained 7/7 for Q3(a) and no students obtained 6/6 for Q3(c). Finally, we can recommend that the Numerical Method course be moved toward more student-centered Learning classrooms where the students will be engaged in group discussion rather than having a teacher one man show.
Abstract: Volterra integro-differential equations appear in many models for real life phenomena. Since analytical solutions for this type of differential equations are hard and at times impossible to attain, engineers and scientists resort to numerical solutions that can be made as accurately as possible. Conventionally, numerical methods for ordinary differential equations are adapted to solve Volterra integro-differential equations. In this paper, numerical solution for solving Volterra integro-differential equation using extended trapezoidal method is described. Formulae for the integral and differential parts of the equation are presented. Numerical results show that the extended method is suitable for solving first order Volterra integro-differential equations.
Abstract: In this study, dynamic responses of composite plates on elastic foundations subjected to impact and moving loads are investigated. The first order shear deformation (FSDT) theory is used for moderately thick plates. Pasternak-type (two-parameter) elastic foundation is assumed. Elastic foundation effects are integrated into the governing equations. It is assumed that plate is first hit by a mass as an impact type loading then the mass continues to move on the composite plate as a distributed moving loading, which resembles the aircraft landing on airport pavements. Impact and moving loadings are modeled by a mass-spring-damper system with a wheel. The wheel is assumed to be continuously in contact with the plate after impact. The governing partial differential equations of motion for displacements are converted into the ordinary differential equations in the time domain by using Galerkin’s method. Then, these sets of equations are solved by using the Runge-Kutta method. Several parameters such as vertical and horizontal velocities of the aircraft, volume fractions of the steel rebar in the reinforced concrete layer, and the different touchdown locations of the aircraft tire on the runway are considered in the numerical simulation. The results are compared with those of the ABAQUS, which is a commercial finite element code.
Abstract: The triple diffusive boundary layer flow of nanofluid under the action of constant magnetic field over a non-linear stretching sheet has been investigated numerically. The model includes the effect of Brownian motion, thermophoresis, and cross-diffusion; slip mechanisms which are primarily responsible for the enhancement of the convective features of nanofluid. The governing partial differential equations are transformed into a system of ordinary differential equations (by using group theory transformations) and solved numerically by using variational finite element method. The effects of various controlling parameters, such as the magnetic influence number, thermophoresis parameter, Brownian motion parameter, modified Dufour parameter, and Dufour solutal Lewis number, on the fluid flow as well as on heat and mass transfer coefficients (both of solute and nanofluid) are presented graphically and discussed quantitatively. The present study has industrial applications in aerodynamic extrusion of plastic sheets, coating and suspensions, melt spinning, hot rolling, wire drawing, glass-fibre production, and manufacture of polymer and rubber sheets, where the quality of the desired product depends on the stretching rate as well as external field including magnetic effects.
Abstract: This article addresses the boundary layer flow and heat transfer of Casson fluid over a nonlinearly permeable stretching surface with chemical reaction in the presence of variable magnetic field. The effect of thermal radiation is considered to control the rate of heat transfer at the surface. Using similarity transformations, the governing partial differential equations of this problem are reduced into a set of non-linear ordinary differential equations which are solved by finite difference method. It is observed that the velocity at fixed point decreases with increasing the nonlinear stretching parameter but the temperature increases with nonlinear stretching parameter.
Abstract: MHD chemically reacting viscous fluid flow towards
a vertical surface with slip and convective boundary conditions has
been conducted. The temperature and the chemical species
concentration of the surface and the velocity of the external flow are
assumed to vary linearly with the distance from the vertical surface.
The governing differential equations are modeled and transformed
into systems of ordinary differential equations, which are then solved
numerically by a shooting method. The effects of various parameters
on the heat and mass transfer characteristics are discussed. Graphical
results are presented for the velocity, temperature, and concentration
profiles whilst the skin-friction coefficient and the rate of heat and
mass transfers near the surface are presented in tables and discussed.
The results revealed that increasing the strength of the magnetic field
increases the skin-friction coefficient and the rate of heat and mass
transfers toward the surface. The velocity profiles are increased
towards the surface due to the presence of the Lorenz force, which
attracts the fluid particles near the surface. The rate of chemical
reaction is seen to decrease the concentration boundary layer near the
surface due to the destructive chemical reaction occurring near the
surface.
Abstract: In this study, one dimensional phase change problem
(a Stefan problem) is considered and a numerical solution of this
problem is discussed. First, we use similarity transformation to
convert the governing equations into ordinary differential equations
with its boundary conditions. The solutions of ordinary differential
equation with the associated boundary conditions and interface
condition (Stefan condition) are obtained by using a numerical
approach based on operational matrix of differentiation of shifted
second kind Chebyshev wavelets. The obtained results are compared
with existing exact solution which is sufficiently accurate.
Abstract: In this study, out-of-plane free vibrations of a circular
rods is investigated theoretically. The governing equations for
naturally twisted and curved spatial rods are obtained using
Timoshenko beam theory and rewritten for circular rods. Effects of
the axial and shear deformations are considered in the formulations.
Ordinary differential equations in scalar form are solved analytically
by using transfer matrix method. The circular rods of the mass matrix
are obtained by using straight rod of consistent mass matrix. Free
vibrations frequencies obtained by solving eigenvalue problem. A
computer program coded in MATHEMATICA language is prepared.
Circular beams are analyzed through various examples for free
vibrations analysis. Results are compared with ANSYS results based
on finite element method and available in the literature.