Abstract: This work is the modeling and simulation of fluid flow (liquid) through porous media. This type of flow occurs in many situations of interest in applied sciences and engineering, fluid (oil) consists of several individual substances in pure, single-phase flow is incompressible and isothermal. The porous medium is isotropic, homogeneous optionally, with the rectangular format and the flow is two-dimensional. Modeling of hydrodynamic phenomena incorporates Darcy's law and the equation of mass conservation. Correlations are used to model the density and viscosity of the fluid. A finite volume code is used in the discretization of differential equations. The nonlinearity is treated by Newton's method with relaxation coefficient. The results of the simulation of the pressure and the mobility of liquid flowing through porous media are presented, analyzed, and illustrated.
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: In this work, we present an efficient approach for
solving variable-order time-fractional partial differential equations,
which are based on Legendre and Laguerre polynomials. First, we
introduced the pseudo-operational matrices of integer and variable
fractional order of integration by use of some properties of
Riemann-Liouville fractional integral. Then, applied together with
collocation method and Legendre-Laguerre functions for solving
variable-order time-fractional partial differential equations. Also, an
estimation of the error is presented. At last, we investigate numerical
examples which arise in physics to demonstrate the accuracy of the
present method. In comparison results obtained by the present method
with the exact solution and the other methods reveals that the method
is very effective.
Abstract: In this paper, the DC-link voltage control of DC-DC boost converter–inverter system is proposed. The mathematical model is developed from four different sub-circuits that depended on the switch positions. The developed differential equations are combined to develop the dynamic model. Transfer function is generated from the switched function model. Fluctuation of DC-link voltage causes connected loads malfunction. For this problem, a kind of traditional controller, the PI controller is applied to achieve constant DC-link voltage. The PI controller gains are obtained based on transfer function step response. The simulation work has been studied by using MATLAB/Simulink software and hardware prototype is implemented with a low-cost microcontroller Arduino Nano. Experimental results are collected by using ArduinoIO library package. Closed-loop DC-link voltage control system is tested with various line and load disturbances. It is found that the experimental results give equal responses with the simulation results.
Abstract: In this paper, by constructing a special set and utilizing
fixed point theory, we study the existence and multiplicity of the
positive solutions for systems of nonlinear third-order differential
equations with p-laplacian, which improve and generalize the result
of related paper.
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: The paper discusses the subinterval-based numerical
method for fractional derivative computations. It is now referred
to by its acronym – SubIval. The basis of the method is briefly
recalled. The ability of the method to be applied in time stepping
solvers is discussed. The possibility of implementing a time step size
adaptive solver is also mentioned. The solver is tested on a transient
circuit example. In order to display the accuracy of the solver –
the results have been compared with those obtained by means of a
semi-analytical method called gcdAlpha. The time step size adaptive
solver applying SubIval has been proven to be very accurate as
the results are very close to the referential solution. The solver is
currently able to solve FDE (fractional differential equations) with
various derivative orders for each equation and any type of source
time functions.
Abstract: The impact of the storm leads to accidents even in the case of vessels that meet the computed safety criteria for stability. That is why, in order to clarify the causes of the accident and shipwreck, it is necessary to study the dynamics of the ship under the complex sudden impact of external forces. The task is to determine the movement and landing of the ship in the complex and sudden impact of external forces, i.e. when the ship's load changes over a relatively short period of time. For the solution, a technique was used to study the ship's dynamics, which is based on the compilation of a system of differential equations of motion. A coordinate system was adopted for the equation of motion of the hull and the determination of external forces. As a numerical method of integration, the 4th order Runge-Kutta method was chosen. The results of the calculation show that dynamic deviations were lower for high-altitude vessels. The study of the movement of the hull under a difficult situation is performed: receiving of cargo, impact of a flurry of wind and subsequent displacement of the cargo. The risk of overturning and flooding was assessed.
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: A mathematical model for knowledge acquisition in
teaching and learning is proposed. In this study we adopt the
mathematical model that is normally used for disease modelling
into teaching and learning. We derive mathematical conditions which
facilitate knowledge acquisition. This study compares the effects
of dropping out of the course at early stages with later stages of
learning. The study also investigates effect of individual interaction
and learning from other sources to facilitate learning. The study fits
actual data to a general mathematical model using Matlab ODE45
and lsqnonlin to obtain a unique mathematical model that can be
used to predict knowledge acquisition. The data used in this study
was obtained from the tutorial test results for mathematics 2 students
from the Central University of Technology, Free State, South Africa
in the department of Mathematical and Physical Sciences. The study
confirms already known results that increasing dropout rates and
forgetting taught concepts reduce the population of knowledgeable
students. Increasing teaching contacts and access to other learning
materials facilitate knowledge acquisition. The effect of increasing
dropout rates is more enhanced in the later stages of learning
than earlier stages. The study opens up a new direction in further
investigations in teaching and learning using differential equations.
Abstract: Minimizing the weight in flexible structures means
reducing material and costs as well. However, these structures could
become prone to vibrations. Attenuating these vibrations has become
a pivotal engineering problem that shifted the focus of many research
endeavors. One technique to do that is to design and implement
an active control system. This system is mainly composed of a
vibrating structure, a sensor to perceive the vibrations, an actuator
to counteract the influence of disturbances, and finally a controller to
generate the appropriate control signals. In this work, two different
techniques are explored to create two different mathematical models
of an active control system. The first model is a finite element model
with a reduced number of nodes and it is called a super-element.
The second model is in the form of state-space representation, i.e.
a set of partial differential equations. The damping coefficients are
calculated and incorporated into both models. The effectiveness of
these models is demonstrated when the system is excited by its first
natural frequency and an active control strategy is developed and
implemented to attenuate the resulting vibrations. Results from both
modeling techniques are presented and compared.
Abstract: This paper focuses on the study of two dimensional magnetohydrodynamic (MHD) steady incompressible viscous Williamson nanofluid with exponential internal heat generation containing gyrotactic microorganism over a stretching sheet. The governing equations and auxiliary conditions are reduced to a set of non-linear coupled differential equations with the appropriate boundary conditions using similarity transformation. The transformed equations are solved numerically through spectral relaxation method. The influences of various parameters such as Williamson parameter γ, power constant λ, Prandtl number Pr, magnetic field parameter M, Peclet number Pe, Lewis number Le, Bioconvection Lewis number Lb, Brownian motion parameter Nb, thermophoresis parameter Nt, and bioconvection constant σ are studied to obtain the momentum, heat, mass and microorganism distributions. Moment, heat, mass and gyrotactic microorganism profiles are explored through graphs and tables. We computed the heat transfer rate, mass flux rate and the density number of the motile microorganism near the surface. Our numerical results are in better agreement in comparison with existing calculations. The Residual error of our obtained solutions is determined in order to see the convergence rate against iteration. Faster convergence is achieved when internal heat generation is absent. The effect of magnetic parameter M decreases the momentum boundary layer thickness but increases the thermal boundary layer thickness. It is apparent that bioconvection Lewis number and bioconvection parameter has a pronounced effect on microorganism boundary. Increasing brownian motion parameter and Lewis number decreases the thermal boundary layer. Furthermore, magnetic field parameter and thermophoresis parameter has an induced effect on concentration profiles.
Abstract: In this paper, a theoretical study on the forced vibration of one degree of freedom system equipped with inerter, working under load-type or displacement-type excitation, is presented. Differential equations of movement are solved under cosinusoidal excitation, and explicit relations for the magnitude, resonant magnitude, phase angle, resonant frequency, and critical frequency are obtained. Influence of the inertance and damping on these dynamic characteristics is clarified. From the obtained results, one concludes that the inerter increases the magnitude of vibration and the phase angle of the damped mechanical system. Moreover, the magnitude ratio and difference of phase angles are not depending on the actual type of excitation. Consequently, such kind of similitude allows for the comparison of various theoretical and experimental results, which can be broadly found in the literature.
Abstract: In this paper, the problem proposed by Von Karman is treated in the attendance of additional flow field effects when the liquid is spaced above the rotating rigid disk. To be more specific, a purely viscous fluid flow yield by rotating rigid disk with Navier’s condition is considered in both magnetohydrodynamic and hydrodynamic frames. The rotating flow regime is manifested with heat source/sink and chemically reactive species. Moreover, the features of thermophoresis and Brownian motion are reported by considering nanofluid model. The flow field formulation is obtained mathematically in terms of high order differential equations. The reduced system of equations is solved numerically through self-coded computational algorithm. The pertinent outcomes are discussed systematically and provided through graphical and tabular practices. A simultaneous way of study makes this attempt attractive in this sense that the article contains dual framework and validation of results with existing work confirms the execution of self-coded algorithm for fluid flow regime over a rotating rigid disk.
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: Dependencies between diverse factors involved in probabilistic seismic loss evaluation are recognized to be an imperative issue in acquiring accurate loss estimates. Dependencies among component damage costs could be taken into account considering two partial distinct states of independent or perfectly-dependent for component damage states; however, in our best knowledge, there is no available procedure to take account of loss dependencies in story level. This paper attempts to present a method called "modal cost superposition method" for decoupling story damage costs subjected to earthquake ground motions dealt with closed form differential equations between damage cost and engineering demand parameters which should be solved in complex system considering all stories' cost equations by the means of the introduced "substituted matrixes of mass and stiffness". Costs are treated as probabilistic variables with definite statistic factors of median and standard deviation amounts and a presumed probability distribution. To supplement the proposed procedure and also to display straightforwardness of its application, one benchmark study has been conducted. Acceptable compatibility has been proven for the estimated damage costs evaluated by the new proposed modal and also frequently used stochastic approaches for entire building; however, in story level, insufficiency of employing modification factor for incorporating occurrence probability dependencies between stories has been revealed due to discrepant amounts of dependency between damage costs of different stories. Also, more dependency contribution in occurrence probability of loss could be concluded regarding more compatibility of loss results in higher stories than the lower ones, whereas reduction in incorporation portion of cost modes provides acceptable level of accuracy and gets away from time consuming calculations including some limited number of cost modes in high mode situation.
Abstract: Detailed numerical calculations are illustrated in our investigation for unsteady natural convection heat and mass transfer of non-Newtonian Casson fluid along a vertical wavy surface. The surface of the plate is kept at a constant temperature and uniform concentration. To transform the complex wavy surface to a flat plate, a simple coordinate transformation is employed. The resulting partial differential equations are solved using the fully implicit finite difference method with SUR procedure. Flow and heat transfer characteristics are investigated for a wide range of values of the Casson parameter, the dimensionless time parameter, the buoyancy ratio and the amplitude-wavelength parameter. It is found that, the variations of the Casson parameter have significant effects on the fluid motion, heat and mass transfer. Also, the maximum and minimum values of the local Nusselt and Sherwood numbers increase by increase either the Casson parameter or the buoyancy ratio.
Abstract: The movement of points feet of the anthropomorphous robot in space occurs along some stable trajectory of a known form. A large number of modifications to the methods of control of biped robots indicate the fundamental complexity of the problem of stability of the program trajectory and, consequently, the stability of the control for the deviation for this trajectory. Existing gait generators use piecewise interpolation of program trajectories. This leads to jumps in the acceleration at the boundaries of sites. Another interpolation can be realized using differential equations with fractional derivatives. In work, the approach to synthesis of generators of program trajectories is considered. The resulting system of nonlinear differential equations describes a smooth trajectory of movement having rectilinear sites. The method is based on the theory of an asymptotic stability of invariant sets. The stability of such systems in the area of localization of oscillatory processes is investigated. The boundary of the area is a bounded closed surface. In the corresponding subspaces of the oscillatory circuits, the resulting stable limit cycles are curves having rectilinear sites. The solution of the problem is carried out by means of synthesis of a set of the continuous smooth controls with feedback. The necessary geometry of closed trajectories of movement is obtained due to the introduction of high-order nonlinearities in the control of stabilization systems. The offered method was used for the generation of trajectories of movement of point’s feet of the anthropomorphous robot. The synthesis of the robot's program movement was carried out by means of the inverse method.
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: The convective and radiative heat transfer performance and entropy generation on forced convection through a direct absorption solar collector (DASC) is investigated numerically. Four different fluids, including Cu-water nanofluid, Al2O3-waternanofluid, TiO2-waternanofluid, and pure water are used as the working fluid. Entropy production has been taken into account in addition to the collector efficiency and heat transfer enhancement. Penalty finite element method with Galerkin’s weighted residual technique is used to solve the governing non-linear partial differential equations. Numerical simulations are performed for the variation of mass flow rate. The outcomes are presented in the form of isotherms, average output temperature, the average Nusselt number, collector efficiency, average entropy generation, and Bejan number. The results present that the rate of heat transfer and collector efficiency enhance significantly for raising the values of m up to a certain range.