Abstract: We develop a method based on polynomial quintic
spline for numerical solution of fourth-order non-homogeneous
parabolic partial differential equation with variable coefficient. By
using polynomial quintic spline in off-step points in space and
finite difference in time directions, we obtained two three level
implicit methods. Stability analysis of the presented method has been
carried out. We solve four test problems numerically to validate the
derived method. Numerical comparison with other methods shows
the superiority of presented scheme.
Abstract: We present a robust nonlinear parabolic partial
differential equation (PDE)-based denoising scheme in this article.
Our approach is based on a second-order anisotropic diffusion model
that is described first. Then, a consistent and explicit numerical
approximation algorithm is constructed for this continuous model by
using the finite-difference method. Finally, our restoration
experiments and method comparison, which prove the effectiveness
of this proposed technique, are discussed in this paper.
Abstract: By using fixed point theorems for a class of
generalized concave and convex operators, the positive solution of
nonlinear fractional differential equation with integral boundary
conditions is studied, where n ≥ 3 is an integer, μ is a parameter
and 0 ≤ μ < α. Its existence and uniqueness is proved, and an
iterative scheme is constructed to approximate it. Finally, two
examples are given to illustrate our results.
Abstract: In population dynamics the study of both, the
abundance and the spatial distribution of the populations in a
given habitat, is a fundamental issue a From ecological point of
view, the determination of the factors influencing such changes
involves important problems. In this paper a mathematical model to
describe the temporal dynamic and the spatiotemporal dynamic of the
interaction of three populations (pollinators, plants and herbivores) is
presented. The study we present is carried out by stages: 1. The
temporal dynamics and 2. The spatio-temporal dynamics. In turn,
each of these stages is developed by considering three cases which
correspond to the dynamics of each type of interaction. For instance,
for stage 1, we consider three ODE nonlinear systems describing
the pollinator-plant, plant-herbivore and plant-pollinator-herbivore,
interactions, respectively. In each of these systems different types of
dynamical behaviors are reported. Namely, transcritical and pitchfork
bifurcations, existence of a limit cycle, existence of a heteroclinic
orbit, etc. For the spatiotemporal dynamics of the two mathematical
models a novel factor are introduced. This consists in considering
that both, the pollinators and the herbivores, move towards those
places of the habitat where the plant population density is high.
In mathematical terms, this means that the diffusive part of the
pollinators and herbivores equations depend on the plant population
density. The analysis of this part is presented by considering pairs of
populations, i. e., the pollinator-plant and plant-herbivore interactions
and at the end the two mathematical model is presented, these models
consist of two coupled nonlinear partial differential equations of
reaction-diffusion type. These are defined on a rectangular domain
with the homogeneous Neumann boundary conditions. We focused
in the role played by the density dependent diffusion term into
the coexistence of the populations. For both, the temporal and
spatio-temporal dynamics, a several of numerical simulations are
included.
Abstract: This paper revisits the free vibration problem of delaminated composite beams. It is shown that during the vibration of composite beams the delaminated parts are subjected to the parametric excitation. This can lead to the dynamic buckling during the motion of the structure. The equation of motion includes time-dependent stiffness and so it leads to a system of Mathieu-Hill differential equations. The free vibration analysis of beams is carried out in the usual way by using beam finite elements. The dynamic buckling problem is investigated locally, and the critical buckling forces are determined by the modified harmonic balance method by using an imposed time function of the motion. The stability diagrams are created, and the numerical predictions are compared to experimental results. The most important findings are the critical amplitudes at which delamination buckling takes place, the stability diagrams representing the instability of the system, and the realistic mode shape prediction in contrast with the unrealistic results of models available in the literature.
Abstract: A designing of a structure requires its realization on rough or sloping ground. Besides the problem of the stability of the landslide, the behavior of the foundations that are bearing the structure is influenced by the destabilizing effect of the ground’s slope. This article focuses on the analysis of the slope stability exposed to loading by introducing the different factors influencing the slope’s behavior on the one hand, and on the influence of this slope on the foundation’s behavior on the other hand. This study is about the elastoplastic modelization using FLAC 2D. This software is based on the finite difference method, which is one of the older methods of numeric resolution of differential equations system with initial and boundary conditions. It was developed for the geotechnical simulation calculation. The aim of this simulation is to demonstrate the notable effect of shear modulus « G », cohesion « C », inclination angle (edge) « β », and distance between the foundation and the head of the slope on the stability of the slope as well as the stability of the foundation. In our simulation, the slope is constituted by homogenous ground. The foundation is considered as rigid/hard; therefore, the loading is made by the application of the vertical strengths on the nodes which represent the contact between the foundation and the ground.
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: The problem under research is that of unpredictable modes occurring in two-stage centrifugal hydraulic pump as a result of hydraulic processes caused by vibrations of structural components. Numerical, analytical and experimental approaches are considered. A hypothesis was developed that the problem of unpredictable pressure decrease at the second stage of centrifugal pumps is caused by cavitation effects occurring upon vibration. The problem has been studied experimentally and theoretically as of today. The theoretical study was conducted numerically and analytically. Hydroelastic processes in dynamic “liquid – deformed structure” system were numerically modelled and analysed. Using ANSYS CFX program engineering analysis complex and computing capacity of a supercomputer the cavitation parameters were established to depend on vibration parameters. An influence domain of amplitudes and vibration frequencies on concentration of cavitation bubbles was formulated. The obtained numerical solution was verified using CFM program package developed in PNRPU. The package is based on a differential equation system in hyperbolic and elliptic partial derivatives. The system is solved by using one of finite-difference method options – the particle-in-cell method. The method defines the problem solution algorithm. The obtained numerical solution was verified analytically by model problem calculations with the use of known analytical solutions of in-pipe piston movement and cantilever rod end face impact. An infrastructure consisting of an experimental fast hydro-dynamic processes research installation and a supercomputer connected by a high-speed network, was created to verify the obtained numerical solutions. Physical experiments included measurement, record, processing and analysis of data for fast processes research by using National Instrument signals measurement system and Lab View software. The model chamber end face oscillated during physical experiments and, thus, loaded the hydraulic volume. The loading frequency varied from 0 to 5 kHz. The length of the operating chamber varied from 0.4 to 1.0 m. Additional loads weighed from 2 to 10 kg. The liquid column varied from 0.4 to 1 m high. Liquid pressure history was registered. The experiment showed dependence of forced system oscillation amplitude on loading frequency at various values: operating chamber geometrical dimensions, liquid column height and structure weight. Maximum pressure oscillation (in the basic variant) amplitudes were discovered at loading frequencies of approximately 1,5 kHz. These results match the analytical and numerical solutions in ANSYS and CFM.
Abstract: In this paper, the influence of van der Waals, as well as electrostatic forces on the structural behavior of MEMS and NEMS actuators, has been investigated using of a Euler-Bernoulli beam continuous model. In the proposed nonlinear model, the electrostatic fringing-fields and the mid-plane stretching (geometric nonlinearity) effects have been considered. The nonlinear integro-differential equation governing the static structural behavior of the actuator has been derived. An original Galerkin-based reduced-order model has been developed to avoid problems arising from the nonlinearities in the differential equation. The obtained reduced-order model equations have been solved numerically using the Newton-Raphson method. The basic design parameters such as the pull-in parameters (voltage and deflection at pull-in), as well as the detachment length due to the van der Waals force of some investigated micro- and nano-actuators have been calculated. The obtained numerical results have been compared with some other existing methods (finite-elements method and finite-difference method) and the comparison showed good agreement among all assumed numerical techniques.
Abstract: In this work we present a family of new convergent
type methods splitting high order no negative steps feature that
allows your application to irreversible problems. Performing affine
combinations consist of results obtained with Trotter Lie integrators
of different steps. Some examples where applied symplectic
compared with methods, in particular a pair of differential equations
semilinear. The number of basic integrations required is comparable
with integrators symplectic, but this technique allows the ability
to do the math in parallel thus reducing the times of which
exemplify exhibiting some implementations with simple schemes for
its modularity and scalability process.
Abstract: In this paper, we present an analytical method for
analysis of nano-scale spherical shell subjected to thermo-mechanical
shocks based on nonlocal elasticity theory. Thermo-mechanical
properties of nano shpere is assumed to be temperature dependent.
Governing partial differential equation of motion is solved
analytically by using Laplace transform for time domain and power
series for spacial domain. The results in Laplace domain is
transferred to time domain by employing the fast inverse Laplace
transform (FLIT) method. Accuracy of present approach is assessed
by comparing the the numerical results with the results of published
work in literature. Furtheremore, the effects of non-local parameter
and wall thickness on the dynamic characteristics of the nano-sphere
are studied.
Abstract: Torrefaction of biomass pellets is considered as a
useful pretreatment technology in order to convert them into a high
quality solid biofuel that is more suitable for pyrolysis, gasification,
combustion, and co-firing applications. In the course of torrefaction,
the temperature varies across the pellet, and therefore chemical
reactions proceed unevenly within the pellet. However, the
uniformity of the thermal distribution along the pellet is generally
assumed. The torrefaction process of a single cylindrical pellet is
modeled here, accounting for heat transfer coupled with chemical
kinetics. The drying sub-model was also introduced. The nonstationary
process of wood pellet decomposition is described by the
system of non-linear partial differential equations over the
temperature and mass. The model captures well the main features of
the experimental data.
Abstract: By using a fixed point theorem of a sum operator, the
existence and uniqueness of positive solution for a class of
boundary value problem of nonlinear fractional differential equation
is studied. An iterative scheme is constructed to approximate it.
Finally, an example is given to illustrate the main result.
Abstract: Lateral torsional buckling is a global buckling mode
which should be considered in design of slender structural members
under flexure about their strong axis. It is possible to compute the
load which causes lateral torsional buckling of a beam by finite
element analysis, however, closed form equations are needed in
engineering practice for calculation ease which can be obtained by
using energy method. In lateral torsional buckling applications of
energy method, a proper function for the critical lateral torsional
buckling mode should be chosen which can be thought as the
variation of twisting angle along the buckled beam. Accuracy of the
results depends on how close is the chosen function to the exact
mode. Since critical lateral torsional buckling mode of the cantilever
I-beams varies due to material properties, section properties and
loading case, the hardest step is to determine a proper mode function
in application of energy method. This paper presents an approximate function for critical lateral
torsional buckling mode of doubly symmetric cantilever I-beams.
Coefficient matrices are calculated for concentrated load at free end,
uniformly distributed load and constant moment along the beam
cases. Critical lateral torsional buckling modes obtained by presented
function and exact solutions are compared. It is found that the modes
obtained by presented function coincide with differential equation
solutions for considered loading cases.
Abstract: In the present study we have investigated axial
buckling characteristics of nanocomposite beams reinforced by
single-walled carbon nanotubes (SWCNTs). Various types of beam
theories including Euler-Bernoulli beam theory, Timoshenko beam
theory and Reddy beam theory were used to analyze the buckling
behavior of carbon nanotube-reinforced composite beams.
Generalized differential quadrature (GDQ) method was utilized to
discretize the governing differential equations along with four
commonly used boundary conditions. The material properties of the
nanocomposite beams were obtained using molecular dynamic (MD)
simulation corresponding to both short-(10,10) SWCNT and long-
(10,10) SWCNT composites which were embedded by amorphous
polyethylene matrix. Then the results obtained directly from MD
simulations were matched with those calculated by the mixture rule
to extract appropriate values of carbon nanotube efficiency
parameters accounting for the scale-dependent material properties.
The selected numerical results were presented to indicate the
influences of nanotube volume fractions and end supports on the
critical axial buckling loads of nanocomposite beams relevant to
long- and short-nanotube composites.
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: This paper deals with nonlinear vibration analysis
using finite element method for frame structures consisting of elastic
and viscoelastic damping layers supported by multiple nonlinear
concentrated springs with hysteresis damping. The frame is supported
by four nonlinear concentrated springs near the four corners. The
restoring forces of the springs have cubic non-linearity and linear
component of the nonlinear springs has complex quantity to represent
linear hysteresis damping. The damping layer of the frame structures
has complex modulus of elasticity. Further, the discretized equations in
physical coordinate are transformed into the nonlinear ordinary
coupled differential equations using normal coordinate corresponding
to linear natural modes. Comparing shares of strain energy of the
elastic frame, the damping layer and the springs, we evaluate the
influences of the damping couplings on the linear and nonlinear impact
responses. We also investigate influences of damping changed by
stiffness of the elastic frame on the nonlinear coupling in the damped
impact responses.
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: Bringing forth a survey on recent higher order spline
techniques for solving boundary value problems in ordinary
differential equations. Here we have discussed the summary of the
articles since 2000 till date based on higher order splines like Septic,
Octic, Nonic, Tenth, Eleventh, Twelfth and Thirteenth Degree
splines. Comparisons of methods with own critical comments as
remarks have been included.
Abstract: This paper deals with the study of reflection and
transmission characteristics of acoustic waves at the interface of a
semiconductor half-space and elastic solid. The amplitude ratios
(reflection and transmission coefficients) of reflected and transmitted
waves to that of incident wave varying with the incident angles have
been examined for the case of quasi-longitudinal wave. The special
cases of normal and grazing incidence have also been derived with
the help of Gauss elimination method. The mathematical model
consisting of governing partial differential equations of motion and
charge carriers’ diffusion of n-type semiconductors and elastic solid
has been solved both analytically and numerically in the study. The
numerical computations of reflection and transmission coefficients
has been carried out by using MATLAB programming software for
silicon (Si) semiconductor and copper elastic solid. The computer
simulated results have been plotted graphically for Si
semiconductors. The study may be useful in semiconductors,
geology, and seismology in addition to surface acoustic wave (SAW)
devices.