Abstract: In this paper, an advanced Nonlinear Exponential
Model (NEM), able to simulate the uniaxial dynamic behavior of
seismic isolators having a continuously decreasing tangent stiffness
with increasing displacement in the relatively large displacements
range and a hardening or softening behavior at large displacements, is
presented. The mathematical model is validated by comparing the
experimental force-displacement hysteresis loops obtained during
cyclic tests, conducted on a helical wire rope isolator and a recycled
rubber-fiber reinforced bearing, with those predicted analytically.
Good agreement between the experimental and simulated results
shows that the proposed model can be an effective numerical tool to
predict the force-displacement relationship of seismic isolation
devices within the large displacements range. Compared to the
widely used Bouc-Wen model, unable to simulate the response of
seismic isolators at large displacements, the proposed one allows to
avoid the numerical solution of a first order nonlinear ordinary
differential equation for each time step of a nonlinear time history
analysis, thus reducing the computation effort. Furthermore, the
proposed model can simulate the smooth transition of the hysteresis
loops from small to large displacements by adopting only one set of
five parameters determined from the experimental hysteresis loops
having the largest amplitude.
Abstract: In this paper, a one-dimensional (1d) Parallel Elasto-
Plastic Model (PEPM), able to simulate the uniaxial dynamic
behavior of seismic isolators having a continuously decreasing
tangent stiffness with increasing displacement, is presented. The
parallel modeling concept is applied to discretize the continuously
decreasing tangent stiffness function, thus allowing to simulate the
dynamic behavior of seismic isolation bearings by putting linear
elastic and nonlinear elastic-perfectly plastic elements in parallel. The
mathematical model has been validated by comparing the
experimental force-displacement hysteresis loops, obtained testing a
helical wire rope isolator and a recycled rubber-fiber reinforced
bearing, with those predicted numerically. Good agreement between
the simulated and experimental results shows that the proposed
model can be an effective numerical tool to predict the forcedisplacement
relationship of seismic isolators within relatively large
displacements. Compared to the widely used Bouc-Wen model, the
proposed one allows to avoid the numerical solution of a first order
ordinary nonlinear differential equation for each time step of a
nonlinear time history analysis, thus reducing the computation effort,
and requires the evaluation of only three model parameters from
experimental tests, namely the initial tangent stiffness, the asymptotic
tangent stiffness, and a parameter defining the transition from the
initial to the asymptotic tangent stiffness.
Abstract: The solution of the nonlinear dynamic equilibrium equations of base-isolated structures adopting a conventional monolithic solution approach, i.e. an implicit single-step time integration method employed with an iteration procedure, and the use of existing nonlinear analytical models, such as differential equation models, to simulate the dynamic behavior of seismic isolators can require a significant computational effort. In order to reduce numerical computations, a partitioned solution method and a one dimensional nonlinear analytical model are presented in this paper. A partitioned solution approach can be easily applied to base-isolated structures in which the base isolation system is much more flexible than the superstructure. Thus, in this work, the explicit conditionally stable central difference method is used to evaluate the base isolation system nonlinear response and the implicit unconditionally stable Newmark’s constant average acceleration method is adopted to predict the superstructure linear response with the benefit in avoiding iterations in each time step of a nonlinear dynamic analysis. The proposed mathematical model is able to simulate the dynamic behavior of seismic isolators without requiring the solution of a nonlinear differential equation, as in the case of widely used differential equation model. The proposed mixed explicit-implicit time integration method and nonlinear exponential model are adopted to analyze a three dimensional seismically isolated structure with a lead rubber bearing system subjected to earthquake excitation. The numerical results show the good accuracy and the significant computational efficiency of the proposed solution approach and analytical model compared to the conventional solution method and mathematical model adopted in this work. Furthermore, the low stiffness value of the base isolation system with lead rubber bearings allows to have a critical time step considerably larger than the imposed ground acceleration time step, thus avoiding stability problems in the proposed mixed method.
Abstract: Many types of controllers were applied on the continuous stirred tank reactor (CSTR) unit to control the temperature. In this research paper, Proportional-Integral-Derivative (PID) controller are compared with Fuzzy Logic controller for temperature control of CSTR. The control system for temperature non-isothermal of a CSTR will produce a stable response curve to its set point temperature. A mathematical model of a CSTR using the most general operating condition was developed through a set of differential equations into S-function using MATLAB. The reactor model and S-function are developed using m.file. After developing the S-function of CSTR model, User-Defined functions are used to link to SIMULINK file. Results that are obtained from simulation and temperature control were better when using Fuzzy logic control compared to PID control.
Abstract: The heat and mass transfer characteristics of
contaminants in groundwater subjected to a biodegradation reaction
is analyzed by taking into account the thermal diffusion (Soret)
effects. This phenomenon is modulated mathematically by a
system of partial differential equations which govern the motion
of fluid (groundwater) and solid (contaminants) particles. The
numerical results are presented graphically for different values of
the parameters entering into the problem on the velocity profiles of
fluid, contaminants, temperature and concentration profile.
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: We examine two-dimensional oil displacement by water in a petroleum reservoir. The pore fluid is immiscible, and the porous media is homogenous and isotropic in the horizontal direction. Buckley-Leverett theory and a combination of Laplacian and Darcy’s law are used to study the fluid flow through porous media, and the Laplacian that defines the dispersion and diffusion of fluid in the sand using heavy oil is discussed. The reservoir is homogenous in the horizontal direction, as expressed by the partial differential equation. Two main factors which are observed are the water saturation and pressure distribution in the reservoir, and they are evaluated for predicting oil recovery in two dimensions by a physical and mathematical simulation model. We review the numerical simulation that solves difficult partial differential reservoir equations. Based on the numerical simulations, the saturation and pressure equations are calculated by the iterative alternating direction implicit method and the iterative alternating direction explicit method, respectively, according to the finite difference assumption. However, to understand the displacement of oil by water and the amount of water dispersion in the reservoir better, an interpolated contour line of the water distribution of the five-spot pattern, that provides an approximate solution which agrees well with the experimental results, is also presented. Finally, a computer program is developed to calculate the equation for pressure and water saturation and to draw the pressure contour line and water distribution contour line for the reservoir.
Abstract: The effects of thermal radiation, Soret and Dufour
parameters on mixed convection and nanofluid flow over a stretching
sheet in the presence of a magnetic field are investigated. The flow is
subject to temperature dependent viscosity and a chemical reaction
parameter. It is assumed that the nanoparticle volume fraction at the
wall may be actively controlled. The physical problem is modelled
using systems of nonlinear differential equations which have been
solved numerically using a spectral relaxation method. In addition
to the discussion on heat and mass transfer processes, the velocity,
nanoparticles volume fraction profiles as well as the skin friction
coefficient are determined for different important physical parameters.
A comparison of current findings with previously published results
for some special cases of the problem shows an excellent agreement.
Abstract: This article considers the problem of evaluating
infinite-time (or finite-time) ruin probability under a given compound
Poisson surplus process by approximating the claim size distribution
by a finite mixture exponential, say Hyperexponential, distribution. It
restates the infinite-time (or finite-time) ruin probability as a solvable
ordinary differential equation (or a partial differential equation).
Application of our findings has been given through a simulation study.
Abstract: The magnetohydrodynamic (MHD) Falkner-Skan
equations appear in study of laminar boundary layers flow over
a wedge in presence of a transverse magnetic field. The partial
differential equations of boundary layer problems in presence of
a transverse magnetic field are reduced to MHD Falkner-Skan
equation by similarity solution methods. This is a nonlinear ordinary
differential equation. In this paper, we solve this equation via
spectral collocation method based on Bessel functions of the first
kind. In this approach, we reduce the solution of the nonlinear
MHD Falkner-Skan equation to a solution of a nonlinear algebraic
equations system. Then, the resulting system is solved by Newton
method. We discuss obtained solution by studying the behavior
of boundary layer flow in terms of skin friction, velocity, various
amounts of magnetic field and angle of wedge. Finally, the results
are compared with other methods mentioned in literature. We can
conclude that the presented method has better accuracy than others.
Abstract: We analyze stochastic integrals associated with a mutation process. To be specific, we describe the cell population process and derive the differential equations for the joint generating functions for the number of mutants and their integrals in generating functions and their applications. We obtain first-order moments of the processes of the two-way mutation process in first-order moment structure of X (t) and Y (t) and the second-order moments of a one-way mutation process. In this paper, we obtain the limiting behaviour of the integrals in limiting distributions of X (t) and Y (t).
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 work, functionally graded materials (FGMs), subjected to loading, which varies with time has been studied. The material properties of FGM are changing through the thickness of material as power law distribution. The conical shells have been chosen for this study so in the first step capability equations for FGM have been obtained. With Galerkin method, these equations have been replaced with time dependant differential equations with variable coefficient. These equations have solved for different initial conditions with variation methods. Important parameters in loading conditions are semi-vertex angle, external pressure and material properties. Results validation has been done by comparison between with those in previous studies of other researchers.
Abstract: This research aims at obtaining the equations of pulse propagation in nonlinear plasmonic waveguides created with As2S3 chalcogenide materials. Via utilizing Helmholtz equation and first-order perturbation theory, two components of electric field are determined within frequency domain. Afterwards, the equations are formulated in time domain. The obtained equations include two coupled differential equations that considers nonlinear dispersion.
Abstract: The paper describes the availability analysis of milling system of a rice milling plant using probabilistic approach. The subsystems under study are special purpose machines. The availability analysis of the system is carried out to determine the effect of failure and repair rates of each subsystem on overall performance (i.e. steady state availability) of system concerned. Further, on the basis of effect of repair rates on the system availability, maintenance repair priorities have been suggested. The problem is formulated using Markov Birth-Death process taking exponential distribution for probable failures and repair rates. The first order differential equations associated with transition diagram are developed by using mnemonic rule. These equations are solved using normalizing conditions and recursive method to drive out the steady state availability expression of the system. The findings of the paper are presented and discussed with the plant personnel to adopt a suitable maintenance policy to increase the productivity of the rice milling plant.
Abstract: The diffusion-reaction equations are important Partial Differential Equations in mathematical biology, material science, physics, and so on. However, finding efficient numerical methods for diffusion-reaction systems on curved surfaces is still an important and difficult problem. The purpose of this paper is to present a convergent geometric method for solving the reaction-diffusion equations on closed surfaces by an O(r)-LTL configuration method. The O(r)-LTL configuration method combining the local tangential lifting technique and configuration equations is an effective method to estimate differential quantities on curved surfaces. Since estimating the Laplace-Beltrami operator is an important task for solving the reaction-diffusion equations on surfaces, we use the local tangential lifting method and a generalized finite difference method to approximate the Laplace-Beltrami operators and we solve this reaction-diffusion system on closed surfaces. Our method is not only conceptually simple, but also easy to implement.
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: In this paper, a spectral decomposition method is developed for the direct integration of stiff and nonstiff homogeneous linear (ODE) systems with linear, constant, or zero right hand sides (RHSs). The method does not require iteration but obtains solutions at any random points of t, by direct evaluation, in the interval of integration. All the numerical solutions obtained for the class of systems coincide with the exact theoretical solutions. In particular, solutions of homogeneous linear systems, i.e. with zero RHS, conform to the exact analytical solutions of the systems in terms of t.
Abstract: A nonlinear model of the mathematical fluid dynamics which describes the motion of an incompressible viscous rotating fluid in a homogeneous gravitational field is considered. The model is a generalization of the known Navier-Stokes system with the addition of the Coriolis parameter and the equations for changeable density. An explicit algorithm for the solution is constructed, and the proof of the existence and uniqueness theorems for the strong solution of the nonlinear problem is given. For the linear case, the localization and the structure of the spectrum of inner waves are also investigated.
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.