This paper introduces a non-conventional approach with multi-dimensional random sampling to solve a cocaine abuse model with statistical probability. The mean Latin hypercube finite difference (MLHFD) method is proposed for the first time via hybrid integration of the classical numerical finite difference (FD) formula with Latin hypercube sampling (LHS) technique to create a random distribution for the model parameters which are dependent on time t . The LHS technique gives advantage to MLHFD method to produce fast variation of the parameters’ values via number of multidimensional simulations (100, 1000 and 5000). The generated Latin hypercube sample which is random or non-deterministic in nature is further integrated with the FD method to complete one cycle of LHS-FD simulation iteration. This process is repeated until n final iterations of LHS-FD are obtained. The means of these n final solutions (MLHFD solutions) are tabulated, graphed and analyzed. The numerical simulation results of MLHFD for the SEIR model are presented side-by-side with deterministic solutions obtained from the classical FD scheme and homotopy analysis method with Pade approximation (HAM-Pade). The present MLHFD results are also compared with the previous non-deterministic statistical estimations from 1995 to 2015. Good agreement between the two is perceived with small errors. MLHFD method can be used to predict future behavior, range and prediction interval for the epidemic model solutions. The expected profiles of the cocaine abuse subpopulations are projected until the year 2045. Both the statistical estimations and the deterministic results of FD and HAM-Pade are found to be within the MLHFD prediction intervals for all the years and for all the subpopulations considered.
There have been many advances in the solar chimney power plant since 1930 and the first pilot work was built in Spain (Manzanares) that produced 50 KW. The solar chimney power plant is considered of a clean power generation that needs to be investigated to enhance the performance by studying the effect of changing the area of passage of air to enhance the velocity towards the chimney to maximize design velocity. In this experimental and numerical study, the reduction area of solar collector was investigated. The reduction area that mean changing the height of glass cover from the absorbing plate (h1=3.8cm, h2=2.6cm and h3=1.28cm). The numerical study was performed using ANSYS Fluent software package (version 14.0) to solve go
... Show MoreIn this study, experimental and numerical applied of heat distribution due to pulsed Nd: YAG laser surface melting. Experimental side was consists of laser parameters are, pulse duration1.3
This study was chosen because of the entry of our regions into the seismic zone recently, where Diyala governorate was hit by the Halabja earthquake in 2017 by 7.3Mw. Therefore, the impact of earthquakes will be studied on the AL-Mafraq bridge foundations piles located in Iraq- east of Baghdad in Diyala Governorate and the extent of its resistance to the Halabjah, EL-Centro, and Kobe earthquakes with acceleration 0.1g, 0.34g, and 0.58g respectively. After modeling and performing the analysis by using Midas Gts-Nx software, the settlement (mm) results at nine nodes (four nodes for the pile cap and five nodes for the piles) were obtained for each of Halabjah, EL-Centro, and Kobe earthquakes to know the resistance of the br
... Show MoreCurrent numerical research was devoted to investigating the effect of castellated steel beams without and with strengthening. The composite concrete asymmetrical double hot rolled steel channels bolted back to back to obtain a built-up I-shape form are used in this study. The top half part of the steel is smaller than the bottom half part, and the two parts were connected by bolting and welding. The ABAQUS/2019 program employed the same length and conditions of loading for four models: The first model is the reference without castellated and strengthening; the second model was castellated without strengthened; the third model was castellated and strengthened with reactive powder concrete encased in the
... Show MoreThis paper is dealing with non-polynomial spline functions "generalized spline" to find the approximate solution of linear Volterra integro-differential equations of the second kind and extension of this work to solve system of linear Volterra integro-differential equations. The performance of generalized spline functions are illustrated in test examples
In this work, a flat-plate solar air heater (FSAH) and a tubular solar air heater (TSAH) were designed and tested numerically. The work investigates the effect of increasing the contact area between the flowing air and the absorber surface of each heater and predicts the expected results before the fabrication of the experimental rig. Three-dimensional two models were designed and simulated by the ANSYS-FLUENT 16 Program. The solar irradiation and ambient air temperature were measured experimentally on December 1st 2022, at the weather conditions of Baghdad City- Iraq, at three air mass flow rates, 0.012 kg/s, 0.032 kg/s, and 0.052 kg/s. The numerical results showed the advantage in the thermal performance of
... Show MoreThe determination of aerodynamic coefficients by shell designers is a critical step in the development of any projectile design. Of particular interest is the determination of the aerodynamic coefficients at transonic speeds. It is in this speed regime that the critical aerodynamic behavior occurs and a rapid change in the aerodynamic coefficients is observed. Two-dimensional, transonic, flow field computations over projectiles have been made using Euler equations which were used for solution with no special treatment required. In this work a solution algorithm is based on finite difference MacCormack’s technique for solving mixed subsonic-supersonic flow problem. Details of the asymmetrically located shock waves on the projectiles hav
... Show MoreNumerical simulations are carried out to assess the quality of the circular and square apodize apertures in observing extrasolar planets. The logarithmic scale of the normalized point spread function of these apertures showed sharp decline in the radial frequency components reaching to 10-36 and 10-34 respectively and demonstrating promising results. This decline is associated with an increase in the full width of the point spread function. A trade off must be done between this full width and the radial frequency components to overcome the problem of imaging extrasolar planets.
Gurney flap (GF) is well-known as one of the most attractive plain flaps because of the simple configuration and effectiveness in improving the lift of the airfoil. Many studies were conducted, but the effects of GF on the various airfoil types need to be further investigated. This study aimed to clarify the effect of GF in the case of the supercritical airfoil RAE2822. This research includes a steady, two-dimensional computational investigation carried out on the supercritical airfoil type RAE-2822 to analyze Gurney flap (GF) effects on the aerodynamic characteristics of this type of airfoil utilizing the Spalart-Allmaras turbulence model within the commercial software Fluent. The airfoil with the Gurney flap was analyz
... Show MoreThe problem of water scarcity is becoming common in many parts of the world, to overcome part of this problem proper management of water and an efficient irrigation system are needed. Irrigation with a buried vertical ceramic pipe is known as a very effective in the management of irrigation water. The two- dimensional transient flow of water from a buried vertical ceramic pipe through homogenous porous media is simulated numerically using the HYDRUS/2D software. Different values of pipe lengths and hydraulic conductivity were selected. In addition, different values of initial volumetric soil water content were assumed in this simulation as initial conditions. Different value
... Show More