Abstract: One of the major challenges faced in solving initial and boundary problems is how to find approximate solutions with minimal deviation from the exact solution without so much rigor and complications. The Taylor series method provides a simple way of obtaining an infinite series which converges to the exact solution for initial value problems and this method of solution is somewhat limited for a two point boundary problem since the infinite series has to be truncated to include the boundary conditions. In this paper, the Ying Buzu Shu algorithm is used to solve a two point boundary nonlinear diffusion problem for the fourth and sixth order solution and compare their relative error and rate of convergence to the exact solution.
Abstract: In the present study, the free vibration of magnetostrictive nano-plate (MsNP) resting on the Pasternak foundation is investigated. Firstly, the modified couple stress (MCS) and nonlocal elasticity theories are compared together and taken into account to consider the small scale effects; in this paper not only two theories are analyzed but also it improves the MCS theory is more accurate than nonlocal elasticity theory in such problems. A feedback control system is utilized to investigate the effects of a magnetic field. First-order shear deformation theory (FSDT), Hamilton’s principle and energy method are utilized in order to drive the equations of motion and these equations are solved by differential quadrature method (DQM) for simply supported boundary conditions. The MsNP undergoes in-plane forces in x and y directions. In this regard, the dimensionless frequency is plotted to study the effects of small scale parameter, magnetic field, aspect ratio, thickness ratio and compression and tension loads. Results indicate that these parameters play a key role on the natural frequency. According to the above results, MsNP can be used in the communications equipment, smart control vibration of nanostructure especially in sensor and actuators such as wireless linear micro motor and smart nano valves in injectors.
Abstract: This paper presents an iteration method for the numerical solutions of a one-dimensional problem of generalized thermoelasticity with one relaxation time under given initial and boundary conditions. The thermoelastic material with variable properties as a power functional graded has been considered. Adomian’s decomposition techniques have been applied to the governing equations. The numerical results have been calculated by using the iterations method with a certain algorithm. The numerical results have been represented in figures, and the figures affirm that Adomian’s decomposition method is a successful method for modeling thermoelastic problems. Moreover, the empirical parameter of the functional graded, and the lattice design parameter have significant effects on the temperature increment, the strain, the stress, the displacement.
Abstract: In this paper, we present a fast and accurate numerical scheme for the solution of a Laplace equation with Dirichlet boundary conditions. The non-standard finite difference scheme (NSFD) is applied to construct the numerical solutions of a Laplace equation with two different Dirichlet boundary conditions. The solutions obtained using NSFD are compared with the solutions obtained using the standard finite difference scheme (SFD). The NSFD scheme is demonstrated to be reliable and efficient.
Abstract: A formulation of postbuckling analysis of end supported rods under self-weight has been presented by the variational method. The variational formulation involving the strain energy due to bending and the potential energy of the self-weight, are expressed in terms of the intrinsic coordinates. The variational formulation is accomplished by introducing the Lagrange multiplier technique to impose the boundary conditions. The finite element method is used to derive a system of nonlinear equations resulting from the stationary of the total potential energy and then Newton-Raphson iterative procedure is applied to solve this system of equations. The numerical results demonstrate the postbluckled configurations of end supported rods under self-weight. This finite element method based on variational formulation expressed in term of intrinsic coordinate is highly recommended for postbuckling analysis of end-supported rods under self-weight.
Abstract: This article presents a study of the effect of
a contra-rotating component on the efficiency of centrifugal
compressors. A contra-rotating centrifugal compressor (CRCC) is
constructed using two independent rotors, rotating in the opposite
direction and replacing the single rotor of a conventional centrifugal
compressor (REF). To respect the geometrical parameters of the REF
one, two rotors of the CRCC are designed, based on a single rotor
geometry, using the hub and shroud length ratio parameter of the
meridional contour. Firstly, the first rotor is designed by choosing
a value of length ratio. Then, the second rotor is calculated to be
adapted to the fluid flow of the first rotor according aerodynamics
principles. In this study, four values of length ratios 0.3, 0.4, 0.5, and
0.6 are used to create four configurations CF1, CF2, CF3, and CF4
respectively. For comparison purpose, the circumferential velocity at
the outlet of the REF and the CRCC are preserved, which means that
the single rotor of the REF and the second rotor of the CRCC rotate
with the same speed of 16000rpm. The speed of the first rotor in this
case is chosen to be equal to the speed of the second rotor. The CFD
simulation is conducted to compare the performance of the CRCC
and the REF with the same boundary conditions. The results show
that the configuration with a higher length ratio gives higher pressure
rise. However, its efficiency is lower. An investigation over the entire
operating range shows that the CF1 is the best configuration in this
case. In addition, the CRCC can improve the pressure rise as well
as the efficiency by changing the speed of each rotor independently.
The results of changing the first rotor speed show with a 130% speed
increase, the pressure ratio rises of 8.7% while the efficiency remains
stable at the flow rate of the design operating point.
Abstract: This paper presents the Finite Element Method (FEM) for analyzing the failure pattern of rectangular slab with various edge conditions. Non-Linear static analysis is carried out using ANSYS 15 Software. Using SOLID65 solid elements, the compressive crushing of concrete is facilitated using plasticity algorithm, while the concrete cracking in tension zone is accommodated by the nonlinear material model. Smeared reinforcement is used and introduced as a percentage of steel embedded in concrete slab. The behavior of the analyzed concrete slab has been observed in terms of the crack pattern and displacement for various loading and boundary conditions. The finite element results are also compared with the experimental data. One of the other objectives of the present study is to show how similar the crack path found by ANSYS program to those observed for the yield line analysis. The smeared reinforcement method is found to be more practical especially for the layered elements like concrete slabs. The value of this method is that it does not require explicit modeling of the rebar, and thus a much coarser mesh can be defined.
Abstract: Knee collateral ligaments play a significant role in restraining excessive frontal motion (varus/valgus rotations). In this investigation, a multiscale frame was developed based on structural hierarchies of the collateral ligaments starting from the bottom (tropocollagen molecule) to up where the fibred reinforced structure established. Experimental data of failure tensile test were considered as the principal driver of the developed model. This model was calibrated statistically using Bayesian calibration due to the high number of unknown parameters. Then the model is scaled up to fit the real structure of the collateral ligaments and simulated under realistic boundary conditions. Predications have been successful in describing the observed transient response of the collateral ligaments during tensile test under pre- and post-damage loading conditions. Collateral ligaments maximum stresses and strengths were observed near to the femoral insertions, a results that is in good agreement with experimental investigations. Also for the first time, damage initiation and propagation were documented with this model as a function of the cross-link density between tropocollagen molecules.
Abstract: In this paper, longitudinal vibration of a micro-beam in micro-scale fluid media has been investigated. The proposed mathematical model for this study is made up of a micro-beam and a micro-plate at its free end. An AC voltage is applied to the pair of piezoelectric layers on the upper and lower surfaces of the micro-beam in order to actuate it longitudinally. The whole structure is bounded between two fixed plates on its upper and lower surfaces. The micro-gap between the structure and the fixed plates is filled with fluid. Fluids behave differently in micro-scale than macro, so the fluid field in the gap has been modeled based on micro-polar theory. The coupled governing equations of motion of the micro-beam and the micro-scale fluid field have been derived. Due to having non-homogenous boundary conditions, derived equations have been transformed to an enhanced form with homogenous boundary conditions. Using Galerkin-based reduced order model, the enhanced equations have been discretized over the beam and fluid domains and solve simultaneously in order to obtain force response of the micro-beam. Effects of micro-polar parameters of the fluid as characteristic length scale, coupling parameter and surface parameter on the response of the micro-beam have been studied.
Abstract: In the current study, two-dimensional unsteady heat conduction in a functionally graded cylinder is studied analytically. The temperature distribution is in radial and longitudinal directions. Heat conduction coefficients are considered a power function of radius both in radial and longitudinal directions. The proposed solution can exactly satisfy the boundary conditions. Analytical unsteady temperature distribution for different parameters of functionally graded cylinder is investigated. The achieved exact solution is useful for thermal stress analysis of functionally graded cylinders. Regarding the analytical approach, this solution can be used to understand the concepts of heat conduction in functionally graded materials.
Abstract: Rockfills are widely used in civil engineering, such as dams, railways, and airport foundations in mountain areas. A significant long-term post-construction settlement may affect the serviceability or even the safety of rockfill infrastructures. The creep behavior of rockfills is influenced by a number of factors, such as particle size, strength and shape, water condition and stress level. However, the effect of particle shape on rockfill creep still remains poorly understood, which deserves a careful investigation. Particle-based discrete element method (DEM) was used to simulate the creep behavior of rockfills under different boundary conditions. Both angular and rounded particles were considered in this numerical study, in order to investigate the influence of particle shape. The preliminary results showed that angular particles experience more breakages and larger creep strains under one-dimensional compression than rounded particles. On the contrary, larger creep strains were observed in he rounded specimens in the direct shear test. The mechanism responsible for this difference is that the possibility of the existence of key particle in rounded particles is higher than that in angular particles. The above simulations demonstrate that the influence of particle shape on the creep behavior of rockfills can be simulated by DEM properly. The method of DEM simulation may facilitate our understanding of deformation properties of rockfill materials.
Abstract: The main purpose of this study is to show differences between the numerical solution of the flow through the artificial heart valve using Newtonian or non-Newtonian fluid. The simulation was carried out by a commercial computational fluid dynamics (CFD) package based on finite-volume method. An aortic bileaflet heart valve (Sorin Bicarbon) was used as a pattern for model of real heart valve replacement. Computed tomography (CT) was used to gain the accurate parameters of the valve. Data from CT were transferred in the commercial 3D designer, where the model for CFD was made. Carreau rheology model was applied as non-Newtonian fluid. Physiological data of cardiac cycle were used as boundary conditions. Outputs were taken the leaflets excursion from opening to closure and the fluid dynamics through the valve. This study also includes experimental measurement of pressure fields in ambience of valve for verification numerical outputs. Results put in evidence a favorable comparison between the computational solutions of flow through the mechanical heart valve using Newtonian and non-Newtonian fluid.
Abstract: Numerical computation of wave propagation in a large
domain usually requires significant computational effort. Hence, the
considered domain must be truncated to a smaller domain of interest.
In addition, special boundary conditions, which absorb the outward
travelling waves, need to be implemented in order to describe the
system domains correctly. In this work, the linear one dimensional
wave equation is approximated by utilizing the Fourier Galerkin
approach. Furthermore, the artificial boundaries are realized with
absorbing boundary conditions. Within this work, a systematic work
flow for setting up the wave problem, including the absorbing
boundary conditions, is proposed. As a result, a convenient modal
system description with an effective absorbing boundary formulation
is established. Moreover, the truncated model shows high accuracy
compared to the global domain.
Abstract: The biaxial buckling behavior of single-layered graphene sheets (SLGSs) is studied in the present work. To consider the size-effects in the analysis, Eringen’s nonlocal elasticity equations are incorporated into classical plate theory (CLPT). A Generalized Differential Quadrature Method (GDQM) approach is utilized and numerical solutions for the critical buckling loads are obtained. Then, molecular dynamics (MD) simulations are performed for a series of zigzag SLGSs with different side-lengths and with various boundary conditions, the results of which are matched with those obtained by the nonlocal plate model to numerical the appropriate values of nonlocal parameter relevant to each type of boundary conditions.
Abstract: In this paper, we consider a geometric inverse source
problem for the heat equation with Dirichlet and Neumann boundary
data. We will reconstruct the exact form of the unknown source
term from additional boundary conditions. Our motivation is to
detect the location, the size and the shape of source support.
We present a one-shot algorithm based on the Kohn-Vogelius
formulation and the topological gradient method. The geometric
inverse source problem is formulated as a topology optimization
one. A topological sensitivity analysis is derived from a source
function. Then, we present a non-iterative numerical method for the
geometric reconstruction of the source term with unknown support
using a level curve of the topological gradient. Finally, we give
several examples to show the viability of our presented method.
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: 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: Tsunami and inundation modelling due to far field tsunami propagation in a limited area is a very challenging numerical task because it involves many aspects such as the formation of various types of waves and the irregularities of coastal boundaries. To compute the effect of far field tsunami and extent of inland inundation due to far field tsunami along the coastal belts of west coast of Malaysia and Southern Thailand, a formulated boundary condition and a moving boundary condition are simultaneously used. In this study, a boundary fitted curvilinear grid system is used in order to incorporate the coastal and island boundaries accurately as the boundaries of the model domain are curvilinear in nature and the bending is high. The tsunami response of the event 26 December 2004 along the west open boundary of the model domain is computed to simulate the effect of far field tsunami. Based on the data of the tsunami source at the west open boundary of the model domain, a boundary condition is formulated and applied to simulate the tsunami response along the coastal and island boundaries. During the simulation process, a moving boundary condition is initiated instead of fixed vertical seaside wall. The extent of inland inundation and tsunami propagation pattern are computed. Some comparisons are carried out to test the validation of the simultaneous use of the two boundary conditions. All simulations show excellent agreement with the data of observation.
Abstract: The adjoint method has been used as a successful tool to
obtain sensitivity gradients in aerodynamic design and optimisation
for many years. This work presents an alternative approach to the
continuous adjoint formulation that enables one to compute gradients
of a given measure of merit with respect to control parameters other
than those pertaining to geometry. The procedure is then applied to
the steady 2–D compressible Euler and incompressible Navier–Stokes
flow equations. Finally, the results are compared with sensitivities
obtained by finite differences and theoretical values for validation.