Search published articles


Showing 69 results for Lar

J. Ahadiyan, B. Mardasi,
Volume 22, Issue 1 (6-2018)
Abstract

Discharge of contaminants in the acceptor ambient has negative environmental impacts. Extremely shallow acceptor ambient conditions will have a significant impact on the diffusion of the contaminants flow. To achieve the effect of the hydraulic, geometric and environmental conditions of the contaminant flow in the acceptor ambient, an experimental model of surface draining was applied. The model consisted of a flume with 3.2 meters length, 0.6 meters width and 0.9 height. Accordingly, by considering three simple single drainers with the diameters of 1.2, 1.6 and 2.2 centimeters, a dense flow with 20, 40 and 60 g/L concentrations occurred in an acceptor ambient with H/d = 2. The results showed that equilibrium length was extremely decreased by increasing the concentrations; this was such that with increasing the concentration three times from 20 to 60 g/L, this length was decreased from 25 to 10. However,  based on  the results, it was found that the effect of  the changes of  the viscosity parameter over the relative length was significant, as its effect was high in all concentrations to 100; then it tended  to become a constant value.. In this regard, the effect of changes in the surface tension parameters over the relative length was investigated; the results showed that this parameter was always effective in the extremely shallow acceptor ambient at all stages from the beginning to the end of the  progress. It is important to note that this had a constant slope in all concentrations and surface tension had an effect on flow diffusion with a certain trend in all concentrations. In fact, surface tension in all concentrations reached to 5 after the equilibrium length and jet energy dissipation area.

M. Dorosti, M. Shafai Bajestan, F. Amirsalari Maymani,
Volume 22, Issue 2 (9-2018)
Abstract

Scour around structures in the river is the most destructivr factor of these structures. Therefore, different methods have been studied to reduce it. A creation slot is one of these methods in the suitable position that can be modified to control the flow pattern scour. In this study, the effects of using two parallel chord rectangular slots  in the single spur dike  with the opening of the effective spur dike area of 10% and a ratio of a/b = 4 (a and b are the length and width of the slot) on two occasions were investigated, one near the bed (near the slot spur dike) and the other near the water surface (away from the nose of spur dike) on the pattern of erosion and sediment around the nose triangular spur dike in clear water conditions. In all experiments, flow depth and angle of the triangular spur dike installation angle were constant. The experiment was done in different hydraulic conditions (Froude number 0.287, 0.304, 0.322). Eventually the results of spur dike without slot (control experiments) were compared. The results showed the slot near the bed toward the slot near water surface and control experiments had a better performance in reducing the maximum depth of scour and deposition of sediments washed to the outer bank. The reduction of the maximum depth of scour at the slot near the bed in Froude numbers 0.287, 0.302 and 0.322, respectively, was  compared to the control experiments , which were 23%, 13% and 24%, respectively, and then compared to the slot near the water surface, which was 60%, 40% and 32% respectively.

S. Ghasemi, A. Hosseinpur, Sh. Kiani,
Volume 22, Issue 3 (11-2018)
Abstract

The rate of metal transfer from the solid phase to solution is an important factor governing their concentration in the soil solution and its availability. In this research, the release rate of Zn in contaminated soils from Isfahan was studied using solutions citric acid, oxalic acid and malic acid 0/01 M during the period of 2 - 504 hours and its relationship with soil characteristics was investigated. The results showed that low molecular weight organic acids could release Zn in the contaminated soils. The Zn released by acids was in the order citric acid> oxalic acid> malic acid. Variation range of Zn released a solution of citric acid, oxalic acid and malic acid, which was 38/9 -21173, 25/2 - 26761 and 25/5 – 20650
mg/kg of soil. Zn released in citric acid solution was higher than that of the two acids. Based on the determination coefficient and standard error estimates done by the kinetic equations, the release of Zn of the contaminated soils and three acid solution was described by the first order equation, elovich, parabolic diffusion and power function, with a high coefficient of determination and a low standard error. Correlation results showed that for the Zn release with the index of Zn usability, there was a significant correlation at 5% level. Showed Multivariate regression model showed that Calcium carbonate, pH and EC affected characteristics of Zn desorption.

E. Yabbarepour, M. Shafai Bajestan, S. M. Kashefipour,
Volume 22, Issue 4 (3-2019)
Abstract

Channels and surface water are ways for the transfer of pollutants to the environment and human. When any pollutant is spilled into the channel, the pollutant concentrations are decreased after the travel. Reducing the distance is an engineering expedient. To reduce the distance, mixing in water should be increased. Thus, the main goal of the present study was to investigate the effect of the triangular vane on transverse mixing used for control bank erosion. Experiments were carried out in an 80cm width flume. The vane, which was triangular, was made of Plexi-glass with a 30% width of the flume length, 15cm height and 50cm far from the tracer injection. Salt solution was used as the soluble tracer. The experiment was carried out to investigate the effect of the triangular vane for two conditions: with and without vane. The transverse mixing (ez) and complete mixing length were estimated for the two conditions of with and without vane. The results showed with installing the triangular vane, the transverse mixing was increased up to 2.5 and the length of mixing was decreased by 60%, as compared with the tests of no vane.

M. Majedi Asl, M. Fuladipanah,
Volume 22, Issue 4 (3-2019)
Abstract

A labyrinth weir is a nonlinear weir folded in the plan-view which increases the crest length and the flow rate for a given channel width and an upstream flow depth. Nowadays, a labyrinth weir is an attractive alternative for those weirs that have a problem in passing the probable maximum flood. The three-dimensional flow pattern and unlimited geometric parameters provide a major challenge to the designers of these weirs. The present study aimed at determining discharge coefficients of sharp-crested triangular labyrinth weirs using the support vector machine (SVM). The results were compared with the experimental data. For this purpose, 123 laboratory test data including  geometric and hydraulic parameters such as vertex angle (θ), magnification ratio (L/B), head water ratio (h/w), Froude number (Fr), Weber Number (We) and Reynolds number (Re) were used. The results showed that the SVM-based model produced the most accurate results when only three geometric parameters, e.g. (h/w, θ, L/B), were introduced as the input parameters (R2 = 0.974, Root mean square error [RMSE] = 0.0118, mean absolute error [MAE] =0.0112 and mean normal error [MNE] =0.017 for the test stage). Also, for these weirs, polynomials linear and nonlinear regression equations were presented. Finally, the discharge coefficient of sharp-crested triangular labyrinth weirs based on the Rehbock equation was evaluated and compared with the SVM using nonlinear and linear regression methods.

S. Chavoshi,
Volume 22, Issue 4 (3-2019)
Abstract

Regional flood frequency studies are initialized by the delineation of the homogeneous catchments. This study was based on "Region of Influence" concept, aiming to find the similar catchments in the south of Caspian Sea. The methodology utilized the Particle Swarm Optimization Algorithm, PSO, to optimize the fuzzy system over a dataset of catchment properties. The main catchment variables in relation to flood were determined by the principle component analysis method and employed as the inputs in the fuzzy system. Catchments grouping was performed over these fuzzy input variables by the iterative process. The optimum similar groups were obtained by PSO, and the heterogeneous L-moment index was used as the termination criterion for the optimization process. A total of 61 hydrometric stations located in the study area were selected and their relevant catchments' physical, climatic and hydrologic properties in relation to flood were studied. Principle Component Analysis by Variomax Rotation Factor over the catchments datasets tended to four out of 16 physical variables, including area, mean elevation, Gravelious Factor and Form Factor, as the main parameters in terms of homogeneity with 84 percent of accumulative variance. These variables, as well as mean annual rainfall, were used as the input data to define the fuzzy system. PSO algorithm was then employed to optimize the developed fuzzy system. The developed algorithm tended to yield the best result in the 9th iteration with 26 and 22 for the minimum average and the optimum values of cost function, respectively. The topology of the resulting algorithm included inertia weight, local and acceleration rates, the number of generations and population size, with the values of 0.7298, 1.4962, 1.4962, 10 and 5, respectively. This study tended to a total of 61 regions of influence, proportional to the relevant 61 sites. According to the geographical location of the catchments in the region, it could be concluded that the geographical proximity doesn't necessarily involve homogeneity. The obtained results indicated the efficient potential of PSO-FES in the delineation of the homogenous catchments in the study area.

R. Monjezi, M. Heidarnejad, A. R. Masjedi, M. H. Pourmohammadi, A. Kamanbedast,
Volume 23, Issue 2 (9-2019)
Abstract

Nonlinear weirs are regarded as important hydraulic structures for water level adjustment and flow control in channels, rivers and dam reservoirs. One example of non-linear weirs is shaped as curved-zigzag. The crest axis of these weirs is non-linear. At a given width, the crest length is greater than that of the conventional linear weirs. Thus, they achieve a higher flow rate for an identical hydraulic load. This research experimentally focused on the discharge coefficient and flow rate of curved weirs with three different curve radii in two triangular linear and zigzag shapes. The discharge coefficients of these weirs were comparatively explored in terms of the hydraulic performance as a function of the total hydraulic load to weir crest height ratio (hd/P) and curvature angle (θ) (or curve radius). The results indicated that for the same hydraulic load, the increase of θ (the decrease in curve radius) led to a lower discharge coefficient; this was first because of the increased topical rise of water level, and then the more indirect path with a greater curvature through which the flow had to transport. Both factors could negatively affect the water discharge coefficient. In practice, the runoff coefficient at a weir with a curve radius of R/w=1.25 was approximately 8.5% greater than that of a weir with a curve radius of R/w=0.75 under a hydraulic load of 0.2.

M. J. Asadi, S. Shabanlou, M. Najarchi, M. M. Najafizadeh,
Volume 23, Issue 3 (12-2019)
Abstract

In this study, the discharge coefficient of the circular side orifices was predicted using a new hybrid method. Combinations made in this study were divided into two sections: 1) the combination of two algorithms including Particle Swarm Optimization (PSO) and Genetic Algorithm (GA) and providing the PSOGA algorithm 2) using the PSOGA algorithm in order to optimize the Adaptive Neuro Fuzzy Inference Systems (ANFIS) network and providing the ANFIS-PSOGA method. Next, by identifying the parameters affecting on the discharge coefficient of the circular side orifices, 11 different combinations were provided. Then, the sensitivity analysis conducted by ANFIS showed that the Froude number and the ratio of the flow depth to the orifice diameter (Ym/D) were identified as the most effective parameters in modeling the discharge coefficient. Also, the best combination including the Froude number (Fr), the ratio of the main channel width to the side orifice diameter (B/D), the ratio of the orifice crest height to its diameter (W/D) and the ratio of the flow depth to the orifice diameter (Ym/D) for estimating the discharge coefficient was introduced. For this model, the values of Mean Absolute Percentage Error (MAPE), Root Mean Square Error (RMSE) and correlation coefficient (R) were obtained 0.021, 0.020 and 0.871, respectively. Additionally, the performance of the ANFIS-PSOGA method was compared with the ANFIS-PSO and ANFIS methods. The results showed that the ANFIS-PSOGA method for predicting the discharge coefficient was the superior model

Gh. Safarinejadi, M. Heidarnejad, A. Bordbar, M. H. Pourmohammadi, A. Kamanbedast,
Volume 23, Issue 4 (12-2019)
Abstract

The use of free launch jets in flip bucket structures with associated submerged ponds, in the appropriate geological, and topographic and hydraulic conditions, could have significant economic and safety benefits. In this research, the downstream scour phenomenon of a flip bucket jet was investigated in free conditions, as well as in the presence of a trapezoidal and triangular slot in the coastal manner with different layout intervals at flow rate of 9, with a total of 45 experiments. The results of this study showed that the presence of the slot had significant effects on the depth and range of scour, so that the fit bucket jet with alternate triangular slots reduced the scour by about 12.7%, as compared to the no slot mode. Moreover, the maximum scour depth occurred in the bucket mode with the alternate trapezoidal slots in the more favorable interval than the rest of the models. Then, the results were compared with several empirical formulas and Veronese A relation showed closer results to the actual values.

S. Janatrostami, H. Mahmoudpour,
Volume 23, Issue 4 (2-2020)
Abstract

One of the main sources of greenhouse gas (GHG) emissions is the use of energy for groundwater pumping. Reducing energy consumption is very important to achieve the environmental sustainability and decrease the climate change impacts. In this paper, the amount of greenhouse gas emissions from groundwater pumping in the Guilan’s aquifers was investigated. Firstly, groundwater depletion and the types of pumps for water pumping were examined in the current condition; then, the values of consumed energy, GHG emissions and climate change indicator of the current condition were estimated. The primary investigations showed that 55 percent of wells in the studied region had a diesel engine pump, while 51.3 percent of the required energy for groundwater pumping was supplied by electric pump. Calculated total GHG emissions and the value of climate change indicator in the current condition were equal to 8.98 and 7.59 Milion kg CO2 eq , respectively. In order to achieve environmental sustainability and energy security, scenarios of replacing electric pumps and applying solar energy were examined. The results of the scenarios showed that diesel fuel wells had no significant effect on the reduction of greenhouse gas emissions, but the use of solar energy reduced them. GHG emissions, in comparison with the base scenarios, were decreased by 44.4% in June, July and August, respectively, by applying the scenario of using the solar pump in agricultural section. Therefore, it is preferred to apply policies in future planning to use renewable energies such as solar energy instead of diesel and electricity energy.

A. Alizadeh, B. Yaghoubi, S. Shabanlou,
Volume 24, Issue 2 (7-2020)
Abstract

In this study, the discharge coefficient of sharp-crested weirs located on circular channels was modeled using the ANFIS and ANFIS-Firefly (ANFIS-FA) algorithm. Also, the Monte Carlo simulations (MCs) were used to enhance the compatibilities of the soft computing models. However, the k-fold cross validation method (k=5) was used to validate the numerical models. According to the input parameters, four models of ANFIS and ANFIS-FA were introduced. Analyzing the numerical results showed that the superior model simulated the discharge coefficient as a function of the Froude number (Fr) and the ratio of flow depth over weir crest to the weir crest height) h/P(. The values of the mean absolute relative error (MARE), root mean square error (RMSE) and correlation coefficient (R) for the superior model were calculated 0.001, 0.002 and 0.999, respectively. However, the maximum error value for this study was less than 2%. 

S. Khalilian, M. Sarai Tabrizi, H. Babazadeh, A. Saremi,
Volume 24, Issue 4 (2-2021)
Abstract

In the present study, the SWAT hydrological model was developed for the upstream of the Zayandehrood dam to evaluate the inflow to this dam. Accordingly, after entering the meteorological and hydrometric information of the region, the runoff simulation was performed. Due to the high volume of entrances to the Zayandehrood Dam, Shahrokh Castle hydrometric stations were selected as the base station for calibration and validation during the statistical period of 1990-2015. After hydrological simulation and accuracy of results, climate prediction was performed using the fifth model of the climate change for the RCP scenarios. According to the forecast, by using climate change models, the temperature could be assumed to increase in all models and the highest rate of increase would occur under the RCP 8.5 climate scenario. After evaluating climate change in different diffusion scenarios, the runoff of the basin was simulated in the SWAT model. The simulation results of runoff in the catchment area showed that although the amount of rainfall was increased in the region, increasing the temperature had a greater effect, reducing the amount of runoff in the basin. Based on the results of climate change, hydrological simulation was performed using the SWAT model. The results showed that the effect of diffusion scenarios in the region was different, causing an increase in temperature and precipitation. The highest increase was observed in the RCP8.5 scenario, which was consistent with the nature of this emission scenario, with the highest emission of greenhouse gases and carbon dioxide. Then, the evaluation of the hydrological model was done; the results showed that although the amount of rainfall in the region had been increased, the increase in temperature of this basin had a greater effect and efficiency in reducing the amount of runoff.

R Mousavi Zadeh Mojarad, S. H. Tabatabaei, N. Nourmahnad,
Volume 25, Issue 2 (9-2021)
Abstract

The contact angle is a numerical index of differentiation between hydrophilic and hydrophobic soils. The objectives of this research are: 1) assessing different methods such as capillary rise, the molarity of ethanol droplet, repellency index, and sessile drop, and 2) Determining the most efficient method in a typical soil with sandy loam texture. In this study, hydrophobic soil was hydrophobized artificially using stearic acid and according to the water drop penetration time classification method. Calculated contact angles of hydrophilic soil with capillary rise method, the molarity of ethanol droplet method, repellency index (two methods of calculation), and sessile drop method were 89.9, 75.41, (57.81), 56.28, and 58.91, respectively. Using the contact angle measuring device, the contact angle of five hydrophobic levels were 58.91, 104.92, 120.48, 129.96, and 173.07, respectively. According to the precession of the device where the operator is capable to control data and processes and there is no limitation in usage, therefore, the sessile drop method is the most suitable method to measure contact angle. The contact angle of the late method and water drop penetration data are positively correlated (R2 = 0.975).

H. Kazemizadeh, M. Saneie, H. Haji Kandi,
Volume 25, Issue 2 (9-2021)
Abstract

To prevent demolishing bridge piles due to developing the scour hole under the foundation of these piles some solution has been proposed in the literature. One of the important approaches could be installing different geometric of roughness at the downstream and upstream piles sections. This causes the downward flows which are performing the main role in developing scour holes to be marginally decreased. The present study explores the effect of geometric roughness and also, continuity and un-continuity of roughness length on maximum scour holes around bridge pile. Results indicate that due to increasing the length of roughness the developed scour holes were formed by less scour hole depths. Furthermore, continuity of roughness increases the scour hole depths; however, un-continuity causes the height of scour holes to be developed by fewer values. Also, the comparison shows that the length of installed roughness in maximum value is decreasing the scour hole depth constitute 34 percent. Based on the non-linear regression technique an equation has been proposed to predict the maximum scour hole due to different conditions. Comparison between experimental and proposed values shows that the accuracy of the proposed equation has an acceptable error which has been calculated less than 11 percent.

Prof. J. Abedi-Koupai, S. Rahimi, S. Eslamian,
Volume 25, Issue 3 (12-2021)
Abstract

Changing the date of the first fall frost and the last spring frost is an important phenomenon in agriculture that can be one of the consequences of global warming. Using general circulation models (GCMs) is a way to study future climate. In this study, observations of temperature and precipitation were weighted by using Mean Observed Temperature-Precipitation (MOTP) method. This method considers the ability of each model in simulating the difference between the mean simulated temperature and mean precipitation in each month in the baseline period and the corresponding observed values. The model that had more weight, selected as the optimum model because it is expected that the model will be valid for the future. But, these models are not indicative of stationary climate change due to their low spatial resolution. Therefore, in this research, the outputs of GCM models are based on the three emission scenarios A2 and B1 and A1B, downscaled by LARS-WG for Isfahan station. The data were analyzed by SPSS software at a 95% confidence level (P<0.05). The results indicated that in the Isfahan in the future period 2020-2049 based on the three scenarios, as compared with baseline period 1971-2000, the first fall frost will occur later and the last spring frost will occur earlier. The first fall frost will occur later for 2 days (based on the A1B emission scenario) to 5 days (based on the A2 emission scenario) and the last spring frost will occur earlier for 2 days (based on the and B1 emission scenario) to 4 days (based on the A2 emission scenario). Finally, the best distribution functions for the first fall frost and the last spring frost for the baseline period and under climate change were selected and compared using the EasyFit software.


F.z. Asadi, R. Fazloula, A. Emadi,
Volume 25, Issue 3 (12-2021)
Abstract

Investigating and understanding river change issues is one of the important factors in sediment hydraulic sciences and river engineering. These studies can be done with the help of physical, mathematical models, or both, but due to financial and time constraints, mathematical models are more general and often used. In this study, the GSTARS model was used to investigate erosion and sedimentation and select the most appropriate function in 12.5 km in length from the Talar river in Mazandaran Province. Simulation using the 55 sections taken in 2006, the daily flow data of the hydrometric station of the Shirgah, located at the beginning of the rich and characteristics of the river sediment, was done. The calibration and validation of the model with cross sections taken in 2012 showed that Yang's sediment transport equation has the highest correlation with reality and can be used to predict river change. The amount of sediment depleted from the case study using the Yang equation is estimated at 8590 tons per year. Also, the study of longitudinal profiles of the river with different sediment transfer functions showed that the study reach at the end range has an erosion trend and is not capable of sand and gravel mining.

R. Sadeghi Talarposhti, R K. Ebrahimi, A. Horfar,
Volume 25, Issue 4 (3-2022)
Abstract

Protection of rivers’ water quality as the most accessible source of the water supply has always been considered. In this paper, self-purification and the pollution decay coefficient values of Talar River, IRAN were studied based on field measurement of DO, BOD, pH, EC, Nitrate, Phosphate, and Temperature, in four seasons of the year 2018, in tandem with the river simulation and its calibration using QUAL2Kw model and the Streeter-Phelps method. In addition to the modeling and analysis results, the measured laboratory data values of the river water samples are also presented. Based on the results, the DO variations were ranged from 5.15 in summer to 7.47 mg/l in spring and BOD variations ranged from 1.88 in fall to 7.9 mg/l in summer. Also, according to the Streeter-Phelps method the decay coefficient values varied from 1.57 (1/day) in spring to 9.63 (1/day) in fall. The values of the Talar River decay coefficient also varied from 2 in fall to 7.7 (1/day) in summer involving the QUAL2Kw model.

M. Sabouri, A.r. Emadi, R. Fazloula,
Volume 26, Issue 2 (9-2022)
Abstract

A compound sharp-crested weir is often used to measure a wide range of flows with appropriate accuracy in open channels. In this study, experiments were performed to investigate the hydraulic flow through a compound weir of circular-rectangular with changes in hydraulic and geometric parameters in free and submerged flow conditions. The characteristics of the weirs include rectangular spans width of 39 cm, a circular radius of 5, 7.5, and 12.5 cm, and heights of 10 and 15 cm. The results showed that by increasing the radius and height of the Weir, upstream water depth increases around 28.4%. At a constant h/p, the discharge coefficient increases with the increasing radius of the circular arc. Also, in the submerged conditions, the discharge coefficient is less (around 40%) than in the free flow condition, which is due to the resistance of the depth of the created stream against the passage of the flow.

J. Abedi Koupai, A.r. Vahabi,
Volume 27, Issue 2 (9-2023)
Abstract

Awareness of water resources status is essential for the proper management of resources and planning for the future due to the occurrence of climate change in most parts of the world and its impact on different parts of the water cycle. Hence, many studies have been carried out in different regions to analyze the effects of climate change on the hydrological process in the coming periods. The present study examined the effects of climate change on surface runoff using the Atmosphere-Ocean General Circulation Model (AOGCM) in Khomeini Shahr City. The maximum and minimum temperatures and precipitation of the upcoming period (2020-2049) were simulated using a weighted average of three models for each of the minimum and maximum temperatures and precipitation parameters based on the scenario A2 and B1 (pessimistic and optimistic states, respectively) of the AOGCM-AR4 models. The LARS-WG model was also used to measure the downscaling. The HEC-HMS was used to predict runoff. The effects of climate change in the coming period (2020-2049) compared with the observation period (1971-2000), in the A2 scenario, the minimum and maximum temperatures would increase by 1.1 and 1.6 Degrees Celsius, respectively, and the precipitation would decrease 17.8 percent. In the B1 scenario, the minimum and maximum temperatures would increase by 1.1 and 1.4 degrees Celsius, respectively, and the precipitation would decrease by 13 percent. The results of runoff were different in the six scenarios in the way the most runoff reduction is related to the scenario of fixed land use and scenario A2 (22.2% reduction), and the most increase is related to the scenario of 45% urban growth and scenario B1 (5.8% increase). So, according to increase urban texture in the future and consequently enhance the volume of runoff, this volume of runoff can be used to feed groundwater, irrigate gardens, and green space in the city.

R. Samadi, Y. Dinpashoh, A. Fakheri-Fard,
Volume 27, Issue 3 (12-2023)
Abstract

A hydrological parameter affecting the management of water resource systems is changes in the amount and occurrence time of extreme precipitation (OTEP). In this research, the seasonality of precipitation in the Lake Urmia (LU) basin was analyzed using the daily extreme precipitation data of 30 rain gauges in the statistical period of 1991-2018. The uniformity of OTEP was tested by Rayleigh and Kuiper’s tests at 0.1, 0.05, and 0.01 levels. The slope of the trend line for OTEP was estimated using the modified Sen slope estimator. The uniformity of OTEP was rejected at each level. The results revealed two strong seasons: late winter and early spring (S1) and autumn (S2) for OTEP. The results showed a general median seasonality index of 0.3, which changed to 0.82 and 0.9 for S1 and S2, respectively, after dividing the whole year into two seasons. The seasonal strength of S1 was similar in both the western and eastern parts of LU, but the west of the lake was stronger than the eastern part in S2. In S1, negative and positive trends in the OTEP were observed on average in 40% and 60% of the stations, respectively, with corresponding values of 77% and 27% for S2, respectively.


Page 3 from 4     

© 2024 CC BY-NC 4.0 | JWSS - Isfahan University of Technology

Designed & Developed by : Yektaweb