Abstract: In order to study the aerodynamic performance of a
semi-flexible membrane wing, Fluid-Structure Interaction simulations
have been performed. The fluid problem has been modeled using
two different approaches which are the vortex panel method and the
numerical solution of the Navier-Stokes equations. Nonlinear analysis
of the structural problem is performed using the Finite Element
Method. Comparison between the two fluid solvers has been made.
Aerodynamic performance of the wing is discussed regarding its
lift and drag coefficients and they are compared with those of the
equivalent rigid wing.
Abstract: We have developed a new computer program in
Fortran 90, in order to obtain numerical solutions of a system
of Relativistic Magnetohydrodynamics partial differential equations
with predetermined gravitation (GRMHD), capable of simulating
the formation of relativistic jets from the accretion disk of matter
up to his ejection. Initially we carried out a study on numerical
methods of unidimensional Finite Volume, namely Lax-Friedrichs,
Lax-Wendroff, Nessyahu-Tadmor method and Godunov methods
dependent on Riemann problems, applied to equations Euler in
order to verify their main features and make comparisons among
those methods. It was then implemented the method of Finite
Volume Centered of Nessyahu-Tadmor, a numerical schemes that
has a formulation free and without dimensional separation of
Riemann problem solvers, even in two or more spatial dimensions,
at this point, already applied in equations GRMHD. Finally, the
Nessyahu-Tadmor method was possible to obtain stable numerical
solutions - without spurious oscillations or excessive dissipation -
from the magnetized accretion disk process in rotation with respect
to a central black hole (BH) Schwarzschild and immersed in a
magnetosphere, for the ejection of matter in the form of jet over a
distance of fourteen times the radius of the BH, a record in terms
of astrophysical simulation of this kind. Also in our simulations,
we managed to get substructures jets. A great advantage obtained
was that, with the our code, we got simulate GRMHD equations in
a simple personal computer.
Abstract: In this paper, numerical solution of system of
Fredholm and Volterra integral equations by means of the Spline
collocation method is considered. This approximation reduces the
system of integral equations to an explicit system of algebraic
equations. The solution is collocated by cubic B-spline and the
integrand is approximated by the Newton-Cotes formula. The error
analysis of proposed numerical method is studied theoretically. The
results are compared with the results obtained by other methods to
illustrate the accuracy and the implementation of our method.
Abstract: In this paper, we have proposed a numerical method
for solving fuzzy Fredholm integral equation of the second kind. In
this method a combination of orthonormal Bernstein and Block-Pulse
functions are used. In most cases, the proposed method leads to
the exact solution. The advantages of this method are shown by an
example and calculate the error analysis.
Abstract: The generalized wave equation models various
problems in sciences and engineering. In this paper, a new three-time
level implicit approach based on cubic trigonometric B-spline for the
approximate solution of wave equation is developed. The usual finite
difference approach is used to discretize the time derivative while
cubic trigonometric B-spline is applied as an interpolating function in
the space dimension. Von Neumann stability analysis is used to
analyze the proposed method. Two problems are discussed to exhibit
the feasibility and capability of the method. The absolute errors and
maximum error are computed to assess the performance of the
proposed method. The results were found to be in good agreement
with known solutions and with existing schemes in literature.
Abstract: In this paper, a backward semi-Lagrangian scheme
combined with the second-order backward difference formula
is designed to calculate the numerical solutions of nonlinear
advection-diffusion equations. The primary aims of this paper are
to remove any iteration process and to get an efficient algorithm
with the convergence order of accuracy 2 in time. In order to achieve
these objects, we use the second-order central finite difference and the
B-spline approximations of degree 2 and 3 in order to approximate
the diffusion term and the spatial discretization, respectively. For the
temporal discretization, the second order backward difference formula
is applied. To calculate the numerical solution of the starting point
of the characteristic curves, we use the error correction methodology
developed by the authors recently. The proposed algorithm turns out
to be completely iteration free, which resolves the main weakness
of the conventional backward semi-Lagrangian method. Also, the
adaptability of the proposed method is indicated by numerical
simulations for Burgers’ equations. Throughout these numerical
simulations, it is shown that the numerical results is in good
agreement with the analytic solution and the present scheme offer
better accuracy in comparison with other existing numerical schemes.
Abstract: The present study is an analysis of the forced convection heat transfer in porous channel with an oriented jet at the inlet with uniform velocity and temperature distributions. The upper wall is insulated when the bottom one is kept at constant temperature higher than that of the fluid at the entrance. The dynamic field is analysed by the Brinkman-Forchheimer extended Darcy model and the thermal field is traduced by the energy one equation model. The numerical solution of the governing equations is obtained by using the finite volume method. The results mainly concern the effect of Reynolds number, jet angle and thermal conductivity ratio on the flow structure and local and average Nusselt numbers evolutions.
Abstract: In this paper, the numerical solution of optimal control problem (OCP) for systems governed by Volterra integro-differential
(VID) equation is considered. The method is developed by means
of the Legendre wavelet approximation and collocation method. The
properties of Legendre wavelet together with Gaussian integration
method are utilized to reduce the problem to the solution of nonlinear
programming one. Some numerical examples are given to confirm the
accuracy and ease of implementation of the method.
Abstract: In this paper, Semi-orthogonal B-spline scaling
functions and wavelets and their dual functions are presented
to approximate the solutions of integro-differential equations.The
B-spline scaling functions and wavelets, their properties and the
operational matrices of derivative for this function are presented to
reduce the solution of integro-differential equations to the solution of
algebraic equations. Here we compute B-spline scaling functions of
degree 4 and their dual, then we will show that by using them we have
better approximation results for the solution of integro-differential
equations in comparison with less degrees of scaling functions
Abstract: Owing to the lack of synchronization between the solar energy availability and the heat demands in a specific application, the energy storing sub-system is necessary to maintain the continuity of thermal process. The present work is dealing with an active solar heating storing system in which an air solar collector is connected to storing unit where this energy is distributed and provided to the heated space in a controlled manner. The solar collector is a box type absorber where the air flows between a number of vanes attached between the collector absorber and the bottom plate. This design can improve the efficiency due to increasing the heat transfer area exposed to the flowing air, as well as the heat conduction through the metal vanes from the top absorbing surface. The storing unit is a packed bed type where the air is coming from the air collector and circulated through the bed in order to add/remove the energy through the charging / discharging processes, respectively. The major advantage of the packed bed storage is its high degree of thermal stratification. Numerical solution of the packed bed energy storage is considered through dividing the bed into a number of equal segments for the bed particles and solved the energy equation for each segment depending on the neighbor ones. The studied design and performance parameters in the developed simulation model including, particle size, void fraction, etc. The final results showed that the collector efficiency was fluctuated between 55%-61% in winter season (January) under the climatic conditions of Misurata in Libya. Maximum temperature of 52ºC is attained at the top of the bed while the lower one is 25ºC at the end of the charging process of hot air into the bed. This distribution can satisfy the required load for the most house heating in Libya.
Abstract: Flows developed between two parallel disks have
many engineering applications. Two types of non-swirling flows can
be generated in such a domain. One is purely source flow in disc type
domain (outward flow). Other is purely sink flow in disc type domain
(inward flow). This situation often appears in some turbo machinery
components such as air bearings, heat exchanger, radial diffuser,
vortex gyroscope, disc valves, and viscosity meters. The main goal of
this paper is to show the mesh convergence, because mesh
convergence saves time, and economical to run and increase the
efficiency of modeling for both sink and source flow. Then flow field
is resolved using a very fine mesh near-wall, using enhanced wall
treatment. After that we are going to compare this flow using
standard k-epsilon, RNG k-epsilon turbulence models. Lastly
compare some experimental data with numerical solution for sink
flow. The good agreement of numerical solution with the
experimental works validates the current modeling.
Abstract: One of the most important parameters in petroleum reservoirs is the pressure distribution along the reservoir, as the pressure varies with the time and location. A popular method to determine the pressure distribution in a reservoir in the unsteady state regime of flow is applying Darcy’s equation and solving this equation numerically. The numerical simulation of reservoirs is based on these numerical solutions of different partial differential equations (PDEs) representing the multiphase flow of fluids. Pressure profile has obtained in a one dimensional system solving Darcy’s equation explicitly. Changes of pressure profile in three situations are investigated in this work. These situations include section length changes, step time changes and time approach to infinity. The effects of these changes in pressure profile are shown and discussed in the paper.
Abstract: Athermal elastic stress analysis of steel fiber reinforced aluminum laminated composite plate is investigated. Four sides of the composite plate are clamped and subjected to a uniform temperature load. The analysis is performed both analytically and numerically. Laminated composite is manufactured via hot pressing method. The investigation of the effects of the orientation angle is provided. Different orientation angles are used such as [0°/90°]s, [30°/-30°]s, [45°/-45°]s, and [60/-60]s. The analytical solution is obtained via classical laminated composite theory and the numerical solution is obtained by applying finite element method via ANSYS.
Abstract: We present in this paper a useful strategy to solve stochastic partial differential equations (SPDEs) involving stochastic coefficients. Using the Wick-product of higher order and the Wiener-Itˆo chaos expansion, the SPDEs is reformulated as a large system of deterministic partial differential equations. To reduce the computational complexity of this system, we shall use a decomposition-coordination method. To obtain the chaos coefficients in the corresponding deterministic equations, we use a least square formulation. Once this approximation is performed, the statistics of the numerical solution can be easily evaluated.
Abstract: Induction heating computer simulation is a powerful tool for process design and optimization, induction coil design, equipment selection, as well as education and business presentations. The authors share their vast experience in the practical use of computer simulation for different induction heating and heat treating processes. In this paper treated with mathematical modeling and numerical simulation of induction heating furnaces with axisymmetric geometries for the numerical solution, we propose finite element methods combined with boundary (FEM) for the electromagnetic model using COMSOL® Multiphysics Software. Some numerical results for an industrial furnace are shown with high frequency.
Abstract: Herpes zoster is a disease that manifests as a dermatological condition. The characteristic of this disease is an irritating skin rash with blisters. This is often limited to one side of body. From the data of Herpes zoster cases in Thailand, we found that age structure effects to the transmission of this disease. In this study, we construct the age structural model of Herpes zoster in Thailand. The local stability analysis of this model is given. The numerical solutions are shown to confirm the analytical results.
Abstract: A novel method based on Genetic Algorithm to solve the boundary value problems (BVPs) of the Falkner–Skan equation over a semi-infinite interval has been presented. In our approach, we use the free boundary formulation to truncate the semi-infinite interval into a finite one. Then we use the shooting method based on Genetic Algorithm to transform the BVP into initial value problems (IVPs). Genetic Algorithm is used to calculate shooting angle. The initial value problems arisen during shooting are computed by Runge-Kutta Fehlberg method. The numerical solutions obtained by the present method are in agreement with those obtained by previous authors.
Abstract: The objective of the present paper is to theoretically investigate the steady-state performance characteristics of journal bearing of finite width, operating with micropolar lubricant in a turbulent regime. In this analysis, the turbulent shear stress coefficients are used based on the Constantinescu’s turbulent model suggested by Taylor and Dowson with the assumption of parallel and inertia-less flow. The numerical solution of the modified Reynolds equation has yielded the distribution of film pressure which determines the static performance characteristics in terms of load capacity, attitude angle, end flow rate and frictional parameter at various values of eccentricity ratio, non-dimensional characteristics length, coupling number and Reynolds number.
Abstract: By means of Sidi-Israeli’s quadrature rules, mechanical quadrature methods (MQMs) for solving the first kind boundary integral equations (BIEs) of steady state Stokes problem are presented. The convergence of numerical solutions by MQMs is proved based on Anselone’s collective compact and asymptotical compact theory, and the asymptotic expansions with the odd powers of the errors are provided, which implies that the accuracy of the approximations by MQMs possesses high accuracy order O (h3). Finally, the numerical examples show the efficiency of our methods.
Abstract: In this paper the Kuramoto-Sivashinsky equation is solved numerically by collocation method. The solution is approximated as a linear combination of septic B-spline functions. Applying the Von-Neumann stability analysis technique, we show that the method is unconditionally stable. The method is applied on some test examples, and the numerical results have been compared with the exact solutions. The global relative error and L∞ in the solutions show the efficiency of the method computationally.