Abstract: In this paper, the effect of the intake temperature (IT) and intake pressure (IP) on ignition timing and pollutants emission of Homogeneous Charge Compression Ignition (HCCI) engine is investigated. Numerical computations are performed using the CHEMKIN computer code. The numerical temperature obtained using different boundary conditions is compared to published data and a good agreement is assigned. Results show that the HCCI combustion engine is significantly improved by increasing the IT. With a value of IT lower than 390 K, combustion cannot occur. However, with an IT greater than 420 K, the cylinder pressure decreases. An optimum crank rotation angle is achieved by using IT of 420 K. So, we can conclude that the variation of the IT and IP influence notably the emission concentration.
Abstract: In this research, a numerical simulation of an Electrohydrodynamic (EHD) actuator’s effects on the flow around a square cylinder by using a finite volume method has been investigated. This is one of the newest ways for controlling the fluid flows. Two plate electrodes are flush-mounted on the surface of the cylinder and one wire electrode is placed on the line with zero angle of attack relative to the stagnation point and excited with DC power supply. The discharge produces an electric force and changes the local momentum behaviors in the fluid layers. For this purpose, after selecting proper domain and boundary conditions, the electric field relating to the problem has been analyzed and then the results in the form of electrical body force have been entered in the governing equations of fluid field (Navier-Stokes equations). The effect of ionic wind resulted from the Electrohydrodynamic actuator, on the velocity, pressure and the wake behind cylinder has been considered. According to the results, it is observed that the fluid flow accelerates in the nearest wall of the frontal half of the cylinder and the pressure difference between frontal and hinder cylinder is increased.
Abstract: The aim of this work is to modelize the occlusion of a
person with temporomandibular disorders as an evolutionary equation
and approach its solution by the construction and characterizing
of discrete variational splines. To formulate the problem, certain
boundary conditions have been considered. After showing the
existence and the uniqueness of the solution of such a problem, a
convergence result of a discrete variational evolutionary spline is
shown. A stress analysis of the occlusion of a human jaw with
temporomandibular disorders by finite elements is carried out in
FreeFem++ in order to prove the validity of the presented method.
Abstract: Cultural heritage conservation is a challenge for contemporary society. In recent decades, significant resources have been allocated for the conservation and restoration of architectural heritage. Historical buildings were restored, protected and reinforced with the intent to limit the risks of degradation or loss, due to phenomena of structural damage and to external factors such as differential settlements, earthquake effects, etc. The wide diffusion of historic masonry constructions in Italy, Europe and the Mediterranean area requires reliable tools for the evaluation of their structural safety. In this paper is presented a free modal analysis performed on a minor historical architecture located in the village of Bagno Grande, near the city of L’Aquila in Italy. The location is characterized by a complex urban context, seriously damaged by the earthquake of 2009. The aim of this work is to check the structural behavior of a masonry building characterized by several boundary conditions imposed by adjacent buildings and infrastructural facilities.
Abstract: Most of the building materials are considered porous, and composed of solid matrix and pores. In the pores, the moisture can be existed in two phases: liquid and vapor. Thus, the mass balance equation is comprised of various moisture driving potentials that translate the movement of the different existing phases occupying pores and the hygroscopic behavior of a porous construction material. This study suggests to resolve a hygrothermal mathematical model of heat and mass transfers in different porous building materials by a numerical investigation. Thereby, the evolution of temperature and moisture content fields has been processed. So, numerous series of hygrothermal calculation on several cases of wall are exposed. Firstly, a case of monolayer wall of massive wood has been treated. In this part, we have compared the numerical solution of the model on one and two dimensions and the effect of dimensional space has been evaluated. In the second case, three building materials (concrete, wood fiberboard and wooden insulation) are tested separately with the same boundary conditions and their hygrothermal behavior are compared. The evaluation of the exchange of heat and air at the interface between the wall and the interior ambiance is carried.
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 article presents a numerical method to find the
heat flux in an inhomogeneous inverse heat conduction problem with
linear boundary conditions and an extra specification at the terminal.
The method is based upon applying the satisfier function along with
the Ritz-Galerkin technique to reduce the approximate solution of the
inverse problem to the solution of a system of algebraic equations.
The instability of the problem is resolved by taking advantage of
the Landweber’s iterations as an admissible regularization strategy.
In computations, we find the stable and low-cost results which
demonstrate the efficiency of the technique.
Abstract: This paper develops a meshless approach, called Element Free Galerkin (EFG) method, which is based on the weak form Moving Least Squares (MLS) of the partial differential governing equations and employs the interpolation to construct the meshless shape functions. The variation weak form is used in the EFG where the trial and test functions are approximated bye the MLS approximation. Since the shape functions constructed by this discretization have the weight function property based on the randomly distributed points, the essential boundary conditions can be implemented easily. The local weak form of the partial differential governing equations is obtained by the weighted residual method within the simple local quadrature domain. The spline function with high continuity is used as the weight function. The presently developed EFG method is a truly meshless method, as it does not require the mesh, either for the construction of the shape functions, or for the integration of the local weak form. Several numerical examples of two-dimensional static structural analysis are presented to illustrate the performance of the present EFG method. They show that the EFG method is highly efficient for the implementation and highly accurate for the computation. The present method is used to analyze the static deflection of beams and plate hole
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: The aim of presented research was to improve numerical predictions of air parameters distribution in the actual natatorium by the selection of calculation formula of mass flux of moisture emitted from the pool. Selected correlation should ensure the best compliance of numerical results with the measurements' results of these parameters in the facility. The numerical model of the natatorium was developed, for which boundary conditions were prepared on the basis of measurements' results carried out in the actual facility. Numerical calculations were carried out with the use of ANSYS CFX software, with six formulas being implemented, which in various ways made the moisture emission dependent on water surface temperature and air parameters in the natatorium. The results of calculations with the use of these formulas were compared for air parameters' distributions: Specific humidity, velocity and temperature in the facility. For the selection of the best formula, numerical results of these parameters in occupied zone were validated by comparison with the measurements' results carried out at selected points of this zone.
Abstract: The objective of this research work is to investigate
for one dimensional transient radiative transfer equations with
conduction using finite volume method. Within the infrastructure of
finite-volume, we obtain the conservative discretization of the terms
in order to preserve the overall conservative property of finitevolume
schemes. Coupling of conductive and radiative equation
resulting in fluxes is governed by the magnitude of emissivity,
extinction coefficient, and temperature of the medium as well as
geometry of the problem.
The problem under consideration has been solved, for a slab
dominating radiation coupled with transient conduction based on
finite volume method. The boundary conditions are also chosen so as
to give a good model of the discretized form of radiation transfer
equation. The important feature of the present method is flexibility in
specifying the control angles in the FVM, while keeping the
simplicity in the solution procedure.
Effects of various model parameters are examined on the
distributions of temperature, radiative and conductive heat fluxes and
incident radiation energy etc. The finite volume method is considered
to effectively evaluate the propagation of radiation intensity through
a participating medium.
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.
Abstract: Result from the constant dwindle in natural resources,
the alternative way to reduce the costs in our daily life would be urgent
to be found in the near future. As the ancient technique based on the
theory of solar chimney since roman times, the double-skin façade are
simply composed of two large glass panels in purpose of daylighting
and also natural ventilation in the daytime. Double-skin façade is
generally installed on the exterior side of buildings as function as the
window, so there is always a huge amount of passive solar energy the
façade would receive to induce the airflow every sunny day. Therefore,
this article imposes a domestic double-skin window for residential
usage and attempts to improve the volume flow rate inside the cavity
between the panels by the frame geometry design, the installation of
outlet guide plate and the solar energy collection system. Note that the
numerical analyses are applied to investigate the characteristics of flow
field, and the boundary conditions in the simulation are totally based
on the practical experiment of the original prototype. Then we
redesign the prototype from the knowledge of the numerical results
and fluid dynamic theory, and later the experiments of modified
prototype will be conducted to verify the simulation results. The
velocities at the inlet of each case are increase by 5%, 45% and 15%
from the experimental data, and also the numerical simulation results
reported 20% improvement in volume flow rate both for the frame
geometry design and installation of outlet guide plate.
Abstract: Mineral product, waste concrete (fine aggregates),
waste in the optical field, industry, and construction employ separators
to separate solids and classify them according to their size. Various
sorting machines are used in the industrial field such as those operating
under electrical properties, centrifugal force, wind power, vibration,
and magnetic force. Study on separators has been carried out to
contribute to the environmental industry. In this study, we perform
CFD analysis for understanding the basic mechanism of the separation
of waste concrete (fine aggregate) particles from air with a machine
built with a rotor with blades. In CFD, we first performed
two-dimensional particle tracking for various particle sizes for the
model with 1 degree, 1.5 degree, and 2 degree angle between each
blade to verify the boundary conditions and the method of rotating
domain method to be used in 3D. Then we developed 3D numerical
model with ANSYS CFX to calculate the air flow and track the
particles. We judged the capability of particle separation for given size
by counting the number of particles escaping from the domain toward
the exit among 10 particles issued at the inlet. We confirm that
particles experience stagnant behavior near the exit of the rotating
blades where the centrifugal force acting on the particles is in balance
with the air drag force. It was also found that the minimum particle
size that can be separated by the machine with the rotor is determined
by its capability to stay at the outlet of the rotor channels.
Abstract: The modelling of physical phenomena, such as the
earth’s free oscillations, the vibration of strings, the interaction of
atomic particles, or the steady state flow in a bar give rise to Sturm-
Liouville (SL) eigenvalue problems. The boundary applications of
some systems like the convection-diffusion equation, electromagnetic
and heat transfer problems requires the combination of Dirichlet and
Neumann boundary conditions. Hence, the incorporation of Robin
boundary condition in the analyses of Sturm-Liouville problem. This
paper deals with the computation of the eigenvalues and
eigenfunction of generalized Sturm-Liouville problems with Robin
boundary condition using the finite element method. Numerical
solution of classical Sturm–Liouville problem is presented. The
results show an agreement with the exact solution. High results
precision is achieved with higher number of elements.
Abstract: In this study, the three-dimensional cavitating
turbulent flow in a complete Francis turbine is simulated using
mixture model for cavity/liquid two-phase flows. Numerical analysis
is carried out using ANSYS CFX software release 12, and standard k-ε
turbulence model is adopted for this analysis. The computational
fluid domain consist of spiral casing, stay vanes, guide vanes, runner
and draft tube. The computational domain is discretized with a threedimensional
mesh system of unstructured tetrahedron mesh. The
finite volume method (FVM) is used to solve the governing equations
of the mixture model. Results of cavitation on the runner’s blades
under three different boundary conditions are presented and
discussed. From the numerical results it has been found that the
numerical method was successfully applied to simulate the cavitating
two-phase turbulent flow through a Francis turbine, and also
cavitation is clearly predicted in the form of water vapor formation
inside the turbine. By comparison the numerical prediction results
with a real runner; it’s shown that the region of higher volume
fraction obtained by simulation is consistent with the region of runner
cavitation damage.
Abstract: Exact solution of an unsteady MHD flow of elasticoviscous
fluid through a porous media in a tube of elliptic cross
section under the influence of magnetic field and constant pressure
gradient has been obtained in this paper. Initially, the flow is
generated by a constant pressure gradient. After attaining the steady
state, the pressure gradient is suddenly withdrawn and the resulting
fluid motion in a tube of elliptical cross section by taking into
account of the porosity factor and magnetic parameter of the
bounding surface is investigated. The problem is solved in two-stages
the first stage is a steady motion in tube under the influence of a
constant pressure gradient, the second stage concern with an unsteady
motion. The problem is solved employing separation of variables
technique. The results are expressed in terms of a non-dimensional
porosity parameter, magnetic parameter and elastico-viscosity
parameter, which depends on the Non-Newtonian coefficient. The
flow parameters are found to be identical with that of Newtonian case
as elastic-viscosity parameter, magnetic parameter tends to zero, and
porosity tends to infinity. The numerical results were simulated in
MATLAB software to analyze the effect of Elastico-viscous
parameter, porosity parameter, and magnetic parameter on velocity
profile. Boundary conditions were satisfied. It is seen that the effect
of elastico-viscosity parameter, porosity parameter and magnetic
parameter of the bounding surface has significant effect on the
velocity parameter.
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.