• صفحه اصلی
  • بررسی آزمایشگاهی تاثیر زاویه واگرائی و همگرائی بر خصوصیات پرش هیدرولیکی در تبدیل‌های واگرا-همگرا
کد مقاله : WEJ-2402-2409 (R3) بازدید : 58 PDF XML

نوع مقاله: پژوهشی

بررسی آزمایشگاهی تاثیر زاویه واگرائی و همگرائی بر خصوصیات پرش هیدرولیکی در تبدیل‌های واگرا-همگرا

reza Mohammadpour 1 ( هیئت علمی دانشگاه آزاد اسلامی واحد استهبان )
محمد حسین جعفری ابنوی 2 ( دانشگاه آزاد اسلامی - واحد استهبان )
محمد کریم بیرامی 3 ( دانشکده مهندسی عمران- دانشگاه اصفهان )

تاریخ ارسال : 1402/12/10 تاریخ تایید : 1403/05/27

کلید واژه: پرش هیدرولیکی, حوضچه واگرا, تبدیلهای واگرا-همگرا, عمق ثانویه پرش,

چکیده مقاله :

پرش هیدرولیکی بعنوان یکی از عوامل مهم جهت استهلاک انرژی جریان در پایین دست سازه‌های هیدرولیکی نظیر دریچه‌ها و سرریز‌ها مورد استفاده قرار می‌گیرند. نوآوری این تحقیق بررسی پرش هیدرولیکی در حوضچه‌های واگرا-همگرا با در نظر گرفتن تاثیر زاویه واگرائی می‌باشد. کلیه آزمایش‌ها در کانالی به طول 11، عمق 7/0 و عرض 48/0 متر با استفاده از دو مدل حوضچه که بترتیب دارای زاویه واگرایی و همگرایی 8 و 6 درجه و برعکس انجام شد. نتایج نشان میدهد که در حداکثر دبی جریان، نسبت عمق جریان y_2/y_1 و طول جریان L_j/y_1 در مدل نخست به ترتیب به میزان 6/48 و 9/122 درصد و در مدل بعد به ترتیب به میزان 5/35 و 8/95 درصد نسبت به حوضچه کلاسیک کاهش داشته است. از طرفی نسبت افت انرژی در هر دو مدل‌ به ترتیب به میزان 8/26 و 7/23 درصد نسبت به حوضچه‌های کلاسیک افزایش یافته است. مقایسه نتایج بین دو حوضچه واگرا-همگرا نشان داد که افزایش نسبت زاویه واگرایی به همگرایی (θ1/θ2) موجب افزایش عملکرد خصوصیات حوضچه‌ها خواهد شد.

چکیده انگلیسی :

Hydraulic jumps are used as one of the key factors for energy dissipation in the downstream of hydraulic structures such as gates and spillways. This research focuses on studying hydraulic jumps in divergent-convergent basins while considering the influence of the divergence angle.. All experiments were conducted in a channel with a length of 11, depth of 0.7 and width of 0.48 meters using two basins models, which have a divergence and convergence angle of 8 and 6 degrees, respectively, and vice versa. The results show that at the maximum flow rate, the ratio of flow depth y_2/y_1 and flow length L_j/y_1 in the first model is 48.6 and 122.9% respectively and in the second model it is 35.5% and 95.8% respectively has decreased compared to the classic basin. On the other hand, the ratio of energy loss in both models has increased by 26.8 and 23.7 percent, respectively, compared to classic basins. The comparison of the results between the two divergent-convergent basins, showed that increasing the ratio of the angle of divergence to convergence θ1/θ2 will increase the performance of the basins.

منابع و مأخذ:
متن کامل:


*بررسی آزمایشگاهی تاثیر زاویه واگرائی و همگرائی بر خصوصیات پرش هیدرولیکی در تبدیل‌های واگرا-همگرا

 

Exploring the Influence of Divergence and Convergence Angles on Hydraulic Jump Characteristics in Divergent-Convergent Transformations: A Laboratory Study

 

Abstract

Hydraulic jumps are used as one of the key factors for energy dissipation in the downstream of hydraulic structures such as gates and spillways. This research focuses on studying hydraulic jumps in divergent-convergent basins while considering the influence of the divergence angle.. All experiments were conducted in a channel with a length of 11, depth of 0.7 and width of 0.48 meters using two basins models, which have a divergence and convergence angle of 8 and 6 degrees, respectively, and vice versa. The results show that at the maximum flow rate, the ratio of flow depth  and flow length  in the first model is 48.6 and 122.9% respectively and in the second model it is 35.5% and 95.8% respectively has decreased compared to the classic basin. On the other hand, the ratio of energy loss in both models has increased by 26.8 and 23.7 percent, respectively, compared to classic basins. The comparison of the results between the two divergent-convergent basins, showed that increasing the ratio of the angle of divergence to convergence  will increase the performance of the basins.

 

Keywords:Divergent basin, hydraulic jump, divergent-convergent transformations, Second depth.

 

Introduction

Downstream of hydraulic structures, such as valves, overflows, and spillways, the flow velocity is high, and the flow energy is destructive. Therefore, stilling basins are used to control the energy of the flow associated with the hydraulic jump. During a hydraulic jump, the flow transitions from supercritical to subcritical, resulting in the dissipation of flow energy. The design of basins directly influences energy dissipation and the length of the hydraulic jump. Consequently, various geometric models of stilling basins are employed to optimize these factors. One of the most common types of hydraulic jumps is the classic jump that is formed in a rectangular channel with a horizontal bed. The classic hydraulic jump was investigated by Blanger (3) who presented the equations for calculating the secondary depth of the hydraulic jump according to the equation of continuity and the momentum.

Hassanpour et al. (8) studied the characteristics of the hydraulic jump in a gradually expanding rectangular stilling basin. They showed that the sequent depth ratio and relative length of the jump decrease with decreasing divergence ratio. Daneshfraz et al. (5) investigated S-type hydraulic jump in stilling basins with sudden divergence and rough bottom. The results showed that the secondary depth and jump length in the rough bed are reduced by 20 and 16%, respectively, and the secondary depth is reduced by 58.5%. In this type of channels, the energy loss for smooth and rough bed increases by 23.5 and 28.7 percent, respectively, compared to the classical jump.

According to the studies, it can be concluded that most of the researches on hydraulic jump are related to divergent basins and some of them are related to convergent basins. The innovation of this research is to investigate the influence of divergence and convergence angles in divergent-convergent transformations (in supercritical flow it reduces divergence and in subcritical flow convergence reduces the depth) to investigate the characteristics of hydraulic jump in stilling basins. Experiments were performed on two divergent-convergent basins with different angles and their results were compared with each other. At the end, the rate of decrease and increase of hydraulic jump characteristics in this type of basins was checked with classic basins.

 

Methodology

All experiments were conducted in a channel with a length of 11.0 m, depth of 0.7 m, and width of 0.48 m. In this study, to investigate the performance of stilling basin, fourteen experiments with different discharges were conducted on both divergent-convergent basins. As shown in Figure (1), the length of the expansion and contraction channel was 65-85 cm and 85-65 cm, respectively. Moreover, the channel width was 30 cm and 48 cm at the first of expansion and contraction, respectively.

 

img0 

Fig. 1-a) Experimental stilling basin b) divergent-convergent basin model-1; c) divergent-convergent basin model-2

 

Results and Discussion

As shown in Table (1), in maximum discharge(Q=66.1 lit/s ), the maximum ratio of flow depth () and jump length () in the divergent-convergent basin (model 1) are 5.32 and 18.30, while these values in the divergent-convergent channel (model 2) are 5.83 and 20.83, respectively. On the other hand, the maximum value of  in the model 1 and model 2 are 0.77 and 0.73, respectively. The results show that with increasing the Froud number the flow depth ratio () and jump length () increase in both basins, but both of these values in the divergent-convergent basin (model 1) are always less than those in the divergent-convergent basin (m0del 2). For example, for the Fr1 = 6.96, the flow depth ratio in the model 1 is 5.31 and increases to 5.88 in the divergent-convergent basin model 2. Furthermore, for Fr1 = 6.96, the  is 18.30 and 21.18 in the divergent-convergent (model 1) and model 2, respectively.

To evaluate the performance of the divergent-convergent (model 1) and divergent-convergent (model 2) basins, the results were compared to the classic basin. The maximum  and  in the divergent-convergent model 1 basin decrease by 48.6% and 122.9%, respectively, compared to the classical basin, while for the divergent-convergent basin model 2 the reduction is 35.5% and It is 95.7%. Furthermore, the energy loss ratio  in model 1 and model 2 has increased by 26.89% and 23.79%, respectively. The results indicate that the performance of divergent-convergent basins (model 1) is much better than divergent-convergent (model 2) and classic basins.

 

 

 

 

 

Table 1-The hydraulic jump characteristics

 

Q

(l/s)

 

 

Model 1

 

 Model 2

 

Fr1

(ΔE/E1)

(y2/y1)

(Lj/y1)

Fr1

(ΔE/E1)

(y2/y1)

(Lj/y1)

66.1

 

6.90

0.77

5.32

18.30

 

6.68

0.73

5.83

20.83

58.0

 

6.85

0.77

5.31

18.48

 

6.05

0.71

5.32

18.09

56.0

 

7.07

0.78

5.31

19.32

 

6.24

0.72

5.44

19.11

54.5

 

7.25

0.79

5.38

17.50

 

6.50

0.73

5.58

19.53

53.5

 

7.39

0.79

5.38

16.67

 

6.61

0.75

5.48

19.05

49.5

 

7.55

0.80

5.62

17.81

 

6.58

0.74

5.50

18.75

45.0

 

8.36

0.82

6.15

18.75

 

6.46

0.74

5.39

18.95

41.0

 

8.19

0.81

6.03

19.67

 

6.96

0.75

5.88

21.18

40.0

 

8.61

0.81

6.20

18.97

 

7.43

0.76

6.41

21.88

38.0

 

8.49

0.81

6.07

19.43

 

7.40

0.76

6.45

21.61

32.0

 

8.12

0.80

6.15

21.15

 

6.89

0.74

6.03

20.69

29.5

 

9.00

0.82

6.73

21.74

 

7.94

0.78

6.80

24.00

25.5

 

10.3

0.84

7.7

26.32

 

8.31

0.78

7.55

25.91

23.5

 

10.3

0.84

7.8

27.78

 

8.84

0.79

8.00

28.50

 

Conclusions

Hydraulic jumps serve as a fundamental element for energy dissipation downstream of hydraulic structures such as spillways. In this research, the effect of divergence and convergence angle on the characteristics of hydraulic jump in divergent-convergent stilling basins has been investigated and the results were compared with each other. Experiments were conducted in a channel with a length of 11 meters, a width of 0.48 meters and a depth of 0.7 meters. In order to investigate the effect of divergence and convergence angle on hydraulic jump characteristics, from two types of basins with divergence and convergence angle   and  respectively (model 1 and ) and   and   respectively (Model 2 and ) has been used. The results show that at the maximum flow rate, the flow depth ratio () and the flow length ratio () in model 1 are equal to 5.32 and 18.30 respectively and in model 2 equal to 5.83 and It is 20/83. On the other hand, the ratio of energy loss in models 1 and 2 is equal to 0.77 and 0.73, respectively. Comparing the results obtained with classic basins shows that the ratio of flow depth () and flow length () in model 1 is 48.6% and 122.9% respectively, and in model 2 it is 35.5 and 95.8 percent have decreased. Also, the energy loss ratio  in models 1 and 2 has increased by 26.8 and 23.7 percent, respectively, compared to classic basins. Considering that the ratio of flow depth and jump length in the divergent-convergent ponds of models 1 and 2 compared to classical ponds is significantly reduced and the energy loss is increased, it can be concluded that this type of basins reduces the dimensions and also increases the performance of the stilling basin. On the other hand, increasing the ratio of divergence to convergence angle  will increase the performance of diverging-converging basins.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

چكيده

پرش هیدرولیکی بعنوان یکی از عوامل مهم جهت استهلاک انرژی جریان در پایین دست سازه‌های هیدرولیکی نظیر دریچه‌ها و سرریز‌ها مورد استفاده قرار می‌گیرند. نوآوری این تحقیق بررسی پرش هیدرولیکی در حوضچه‌های واگرا-همگرا با در نظر گرفتن تاثیر زاویه واگرائی می‌باشد. کلیه آزمایش‌ها در کانالی به طول 11، عمق 7/0 و عرض 48/0 متر با استفاده از دو مدل حوضچه که بترتیب دارای زاویه واگرایی و همگرایی 8 و 6 درجه و برعکس انجام شد. نتایج نشان می‌دهد که در حداکثر دبی جریان، نسبت عمق جریان  و طول جریان  در مدل نخست به ترتیب به میزان 6/48 و 9/122 درصد و در مدل بعد به ترتیب به میزان 5/35 و 8/95 درصد نسبت به حوضچه کلاسیک کاهش داشته‌ است. از طرفی نسبت افت انرژی در هر دو مدل‌ به ترتیب به میزان 8/26 و 7/23 درصد نسبت به حوضچه‌های کلاسیک افزایش یافته است. مقایسه نتایج بین دو حوضچه واگرا-همگرا نشان داد که افزایش نسبت زاویه واگرایی به همگرایی () موجب افزایش عملکرد خصوصیات حوضچه‌ها خواهد شد.

 

واژه‌هاي كليدي: حوضچه واگرا، پرش هیدرولیکی، تبدیل‌های واگرا-همگرا، عمق ثانویه پرش.

 

 

مقدمه

در انتهای سازه‌های هیدرولیکی مانند دریچه‌ها، سرریز‌ها، تندآب‌ها و غیره سرعت جریان زیاد استبالا بوده  و انرژی جریان مخرب هست می باشد. بنابراین از همین رو برای کنترل انرژی جریان که عموما با پرش هیدرولیکی همراه است از حوضچه‌های آرامش استفاده می‌شود. در پرش هیدرولیکی جریان از حالت فوق بحرانی به زیربحرانی تبدیل و انرژی جریان مستهلک می‌گردد. طراحی حوضچه‌ها تاثیر مستقیم بر روی افت انرژی و طول پرش هیدرولیکی دارد، از این رو حوضچه‌های آرامش با مدل‌های هندسی متفاوت مورد استفاده می‌شود. یکی از معمول ترین نوع پرش‌های هیدرولیکی پرش نوع کلاسیک است که در کانال مستطیلی با بستر افقی شکل می‌گیرد. پرش هیدرولیک کلاسیک توسط بلانگر (34) مورد بررسی قرار گرفت که با توجه به معادله پیوستگی و اندازه حرکت رابطه‌ای برای محاسبه عمق ثانویه پرش هیدرولیکی ارائه کرد. در چند دهه اخیر، پرش هیدرولیکی در کانال‌های مستطیلی با مقاطع واگرای تدریجی و ناگهانی از جمله مواردی بوده که مورد علاقه محققین واقع شده است. نیسی و همکاران (1014) مطالعاتی بر روی پرش در حوضچه‌های آرامش واگرای ناگهانی با بستر زبر و نسبت بازشدگی 33/0 داشتند. نتایج تحقیق آنها نشان داد که  و به  این نتیجه رسیدند که در چنین مقاطعی و در محدوده عدد فرود بین 2 تا 10 عمق ثانویه پرش نسبت به پرش کلاسیک 53 درصد کاهش و افت انرژی 25 درصد افزایش داشته است. ورکی و همکاران (1216) مطالعاتی پیرامون پرش هیدرولیکی در کانال‌های واگرا تدریجی با شیب معکوس بصورت تحلیلی و آزمایشگاهی انجام دادند. آنها روابطی را برای خصوصیات پرش بر اساس پارامترهای شیب کف کانال، زاویه واگرایی و عدد فرود اولیه ارائه کردند. نتایج این تحقیق نشان داد که با افزایش توام درجه واگرایی و شیب کف، عمق ثانویه و طول پرش کاهش و افت انرژی افزایش می‌یابد. دانشفراز و همکاران (54) خصوصیات پرش هیدرولیکی در همگرایی و واگرایی را بصورت آزمایشگاهی بررسی کردند. آنها به این نتیجه رسیدند که در کانالهای همگرا مقدار عمق ثانویه و طول پرش به ترتیب 8/1 و 5/2 سانتیمتر افزایش و افت انرژی حدود 74/8 درصد کاهش داشته است. حسن‌پور و همکاران (812) پرش هیدرولیکی در حوضچه آرامش واگرا با کف زبر را مورد بررسی قرار دادند. این محققین رابطه‌هایی برای اعماق مزدوج، طول و افت انرژی پرش ارائه کردند، . این روابط دارای متغیرهایی از جمله عدد فرود اولیه، نسبت واگرائی و زبری کف کانال می باشد. نتایج حاصله نشان می­دهد که با افزایش عدد فرود اولیه عمق ثانویه افزایش می­یابد. همچنین با افزایش زبری، عمق ثانویه و طول پرش کاهش می­یابد. اشکو و همکاران (710) تاثیر بلوک‌های زاویه دار را بر خصوصیات پرش در کانال واگرای تدریجی با شیب معکوس را بررسی کردند. نتایج تحقیقات آن‌ها نشان داد که بلوک‌ها، نسبت اعماق متوالی و طول پرش را به‌طور قابل ملاحظه‌ای کاهش می‌دهد. مقایسه دو چیدمان بلوک‌ها نشان داد که استفاده از بلوک‌های بافل در چیدمان همگرا خصوصیات پرش هیدرولیکی را بهبود می‌بخشد و همچنین موثر ترین زاویه همگرایی 30 درجه است که نسبت طول و عمق پرش را بترتیب 35 و 16 درجه کاهش می‌دهد. دانشفراز و همکاران (56) پرش هیدرولیکی که پنجه پرش کاملا در مقطع واگرا است را در حوضچه‌های آرامش با واگرایی ناگهانی و کف زبر مورد بررسی قرار دادند. نتایج نشان داد که عمق ثانویه و طول پرش در بستر زبر نسبت به بستر صاف بترتیب 20 و 16 درصد کاهش یافته و همچنین عمق ثانویه به میزان 5/58 درصد کاهش می یابد. در این نوع کانالها افت انرژی برای بستر صاف و زبر بترتیب 5/23 و 7/28 درصد نسبت به پرش کلاسیک افزایش می‌یابد. ترکمن زاد و همکاران (1115) خصوصیات پرش هیدرولیکی را در حوضچه آرامش نامتقارن با واگرایی ناگهانی و بستر زبر مورد بررسی قرار دادند. نتایج این تحقیق نشان داد. که عمق ثانویه و طول پرش در اینگونه حوضچه‌ها در مقایسه با حوضچه‌های متقارن با بستر افقی کاهش می‌یابد همچنین زبری بستر نقش تعیین کننده‌ای در محل وقوع پرش دارد. از طرفی نتایج در مورد افت انرژی پرش بدین گونه بود که در حوضچه‌های آرامش واگرای ناگهانی با بستر زبر نسبت به پرش هیدرولیک کلاسیک کارآمدتر بوده است. دزفولی و همکاران (2) تاثیر زبری بر خصوصیات پرش هیدرولیکی در حوضچه‌های واگرای ناگهانی رابررسی کردند. نتایج این مطالعه نشان داد که زبری‌ها افزایش 31 درصدی طول پرش حاصل از واگرایی ناگهانی نسبت به پرش کلاسیک،موجب کاهش طول پرش به را به میزان 12درصد نسبت کاهش داده استمی گردد و همچنین و از نظر طول حوضچه مشابه با حوضچه USBR-II می‌باشد. همچنین از طرفی زبری‌ها باعث کاهش طول غلتاب به میزان 20 درصد نسبت به پرش در مقطع ناگهانی ولی بستر صاف خواهد شد. دانشفراز و نجیبی (67) استهلاک انرژی جریان فوق بحرانی در تنگ شدگی ناگهانی به همراه زبری دیواره را بررسی کردند. نتایج آزمایشگاهی نشان داد که با افزایش عدد فرود بالادست، افت انرژی نسبی افزایش می‌یابد. آنها نشان دادند که با افزایش زبری، میزان افت انرژی نسبی به‌طور چشمگیری افزایش می یابد. حاجی عزیزی و همکاران (1) ساختار تلاطم پرش هيدروليکي مستغرق در کانال واگراي تدريجي با بستر زبر را بررسی کردند. نتایج تحقیق آنها نشان داد که با افزايش ارتفاع اجزاي زبر، نرخ رشد لاية مرزي افزايش يافته که اين خود باعث افزايش تنش برشي و درنهايت افزايش افت انرژي می‌شود. قادری و همکاران (11) خصوصیات پرش هیدرولیکی آزاد و مستغرق را در زبری‌های مختلف بررسی کردند. نتایج نشان داد که مدل عددی به خوبی قادر به شبیه سازی ویژگی‌های پرش آزاد و غوطه ور است. زبری مثلثی نسبت به سایر اشکال بیشترین تاثیر را بر طول پرش دارد و بیشترین تنش برشی و افت انرژی در هر دو پرش در زبری مثلثی رخ می‌دهد. دانشفراز و همکاران (8) پرش هیدرولیکی در کانال واگرای ناگهانی زبر را بررسی کردند. نتایج نشان داد که تنش برشی در بستر ناهموار یک کانال واگرای ناگهانی بیش از 12 برابر بیشتر از یک کانال منشوری صاف است. همچنین عمق ثانویه و طول پرش S در بستر ناهموار نسبت به بستر صاف به ترتیب 22% و 9-13% کاهش یافت. زینی وند و همکاران (3) مقایسه آزمایشگاهی تاثیر زبری مصنوعی و شیب معکوس بستر در کنترل پرش

هیدرولیکی را بررسی کردند. نتایج این تحقیق نشان داد که نوع زبری در افزایش میزان افت انرژی پرش هیدرولیکی موثر می باشد. تاثیر همزمان اعمال زبری بستر و ایجاد شیب معکوس می تواند میزان افت انرژی در پدیده پرش هیدرولیکی را افزایش دهد. حمیدی‌ نژاد و همکاران (913) خصوصیات پرش هیدرولیکی در حوضچه‌های واگرای ناگهانی با پله منفی را بررسی کردند. محدوده عدد فرود انتخابی آنها بین 5/4-9 بود. نتایج نشان داد که پرش های S نامتقارن بوده و انبساط ناگهانی باعث افزایش طول پرش و اتلاف انرژی  و همچنین کاهش عمق ثانویه می شود، در حالی که عمق متوالی کاهش می یابد.. آنها نشان دادند که با استفاده از پله منفی، پرش هیدرولیک به شکل متقارن تغییر می یابدیافت. همچنین، طول پرش هیدرولیک در مقایسه با کانالها با انبساط ناگهانی و بدون پله، به طور قابل توجهی کاهش داشته اما و همچنین عمق متوالیثانویه پرش افزایش و افت انرژی کاهش یافته است. دانشفراز و همکاران (9) خصوصیات پرش هیدرولیکی و افت انرژی در پایین دست سرریزهای پلکانی با پله‌های ناهموار را بررسی کردند. نتایج نشان داد که سرریزهای پلکانی با پله‌های ناهموار نسبت به پله‌های صاف افت انرژی را در حدود 97/20 درصد افزایش می‌دهند.

با توجه به مطالعات صورت گرفته می‌توان نتیجه گرفت که اکثر تحقیقات انجام شده بر روی پرش هیدرولیکی مربوط به حوضچه‌های واگرا و تعدادی از آنها مربوط به حوضچه‌های همگرا بوده است. نوآوری این تحقیق بررسی تاثیر زوایای واگرائی و همگرائی در تبدیل‌های واگرا-همگرا (در جریان زیربحرانی و فوق بحرانی تبدیل واگرا و همگرا به ترتیب عمق جریان را کاهش می دهد (34) جهت بررسی خصوصیات پرش هیدرولیکی در حوضچه‌های آرامش به منظور کاهش هزینه و ایجاد حوضچه اقتصادی است. آزمایش‌ها بر روی دو حوضچه واگرا-همگرا با زاویه‌های مختلف انجام شده و نتایج آنها با یکدیگر مقایسه شده است. در انتها میزان کاهش و افزایش خصوصیات پرش هیدرولیکی در این نوع حوضچه‌ها با حوضچه‌های کلاسیک بررسی شد.

 

مواد و روش‌ها

آزمایش‌ها در فلوم مستطیلی واقع در آزمایشگاه هیدرولیک دانشگاه آزاد اسلامی-واحد استهبان انجام گرفت. کانال آزمایشگاهی دارای طول 11 متر، عرض 48/0 متر و ارتفاع 7/0 متر می‌باشد. مقطع طولی کانال در شکل 1-الف نشان داده شده است. برای تشکیل پرش هیدرولیکی از یک سرریز اوجی شکل استاندارد با ارتفاع  متر و عرض  متر استفاده شد و تبدیل‌های مورد نظر بلافاصله بعد از سرریز نصب شدند. برای کنترل موقعیت و تشکیل پرش هیدرولیکی از یک دریچه در انتهای کانال استفاده شد. لازم به توضیح است که با تنظیم عمق جریان شرایط آزمایش‌ها به گونه‌ای تنظیم شد که شروع پرش هیدرولیکی در انتهای سرریز و ابتدای حوضچه آرامش تشکیل ‌شود. با شروع آزمایش‌ها در ابتدا اعماق جریان برداشت شد برای برداشت اعماق جریان از عمق سنج (Point Gage) با دقت 1/0 میلی‌متر استفاده شد به این صورت که در سه نقطه از عرض کانال عمق جریان برداشت شده و عمق متوسط آن‌ها به‌عنوان عمق نهایی در نظر گرفته شد. طول پرش هیدرولیکی که عبارت است از فاصله شروع پرش تا نقطه‌ای روی سطح آب، بلافاصله پس از آخرین موج غلطان  اندازه گیری شد(34). طول پرش هیدرولیکی توسط یک خط کش فلزی که در دیواره کانال نصب شده بود با دقت مناسب اندازه گیری شد. پلان حوضچه آرامش با تبدیل واگرا-همگرا و همچنین سرریز اوجی در شکل 1-ب نشان داده شده است. حوضچه­های آرامش واگرا-همگرا به گونه­ای طراحی شد که واگرائی در ابتدا و همگرائی پس از آن قرار گیرد. 

 

img0 

شکل 1-الف و ب مقطع طولی و پلان حوضچه آرامش با تبدیل واگرا-همگرا در انتهای سرریز اوجی در کانال آزمایشگاهی

 

در این تحقیق، برای بررسی عملکرد حوضچه آرامش، آزمایش‌ها بر روی دو نوع حوضچه واگرا- همگرا انجام شد. مشخصات هر دو حوضچه در جدول 1 و شکل 2 نشان داده شده است. شکل2-الف تبدیل واگرا-همگرا را در کانال آزمایشگاهی نشان می‌دهد. در شکل2-ب  جزئیات تبدیل واگرا-همگرا (مدل 1) نشان داده شده است. در این تبدیل جریان ابتدا وارد یک تبدیل واگرا و سپس همگرا خواهد شد. طول‌ تبدیل‌های واگرا و همگرا بترتیب  و  و عرض آنها به ترتیب  و  می‌باشد. لازم به ذکر است که عرض کانال پس از تبدیل همگرا مجددا به cm30 خواهد رسید. در این مدل زاویه واگرایی و همگرایی بترتیب  و  و نسبت زاویه واگرایی به همگرایی برابر با  می­باشد. در شکل 2-ج نیز تبدیل واگرا-همگرا (مدل 2) نشان داده شده است. این مدل برعکس مدل قبل بوده بدین صورت که طول تبدیل‌های واگرا و همگرا به ترتیب  و  و مقدار زاویه واگرایی و همگرایی بترتیب  و   درجه و نسبت زاویه واگرایی به همگرایی برابر با   می­باشد . با توجه به محدودیت‌های کانال آزمایشگاهی و همچنین مطالعات پیشینیان که معمولا با زوایای بیش از چهار درجه انجام گرفته (1216)، در این تحقیق از دو زاویه شش و هشت درجه استفاده شده است.

 

جدول 1-خصوصیات حوضچه‌های آزمایشگاهی

 

نوع حوضچه آرامش

عرض ابتدایی b1(m)

عرض ثانویه b2(m)

طول واگرایی

L1 (m)

طول همگرایی L2 (m)

زاویه واگرایی θ1 (Deg.)

زاویه همگرایی θ2 (Deg.)

 

واگرا-همگرا-مدل 1

3/0

48/0

65/0

85/0

8

6

واگرا-همگرا-مدل 2

3/0

48/0

85/0

65/0

6

8

 

img0 

شکل 2-حوضچه آرامش: الف) حوضچه آرامش در کانال آزمایشگاهی ب) تبدیل واگرا-همگرا با زاویه واگرائی 8 درجه  ج) تبدیل واگرا-همگرا با زاویه واگرائی 6 درجه 

 

تحلیل ابعادی

با توجه به اینکه در این تحقیق حوضچه‌های واگرا-همگرا مورد بررسی و ارزیابی قرار گرفته است، لازم است در ابتدا تحلیل ابعادی به منظور يافتن پارامترهای مؤثر بر پرش هیدرولیکی در اینگونه حوضچه‌ها انجام گیرد. به طورکلی مشخصات پرش هیدرولیکی در این نوع حوضچه‌ها تابعی از پارامترهای زير است.

 

(1)                                  

 

در این رابطه  جرم حجمی سیال،  گرانروی دينامیکی،  عمق اولیه پرش هیدرولیکی،  عمق ثانویه پرش هیدرولیکی،  شتاب ثقل،   سرعت قبل از پرش هیدرولیکی،  طول غلطانی پرش،  افت انرژی پرش،  انرژی در مقطع ابتدایی،   زاویه واگرایی،   زاویه همگرایی،  و   به ترتیب عرض کانال در ابتدا و انتهای تبدیل واگرا می‌باشد. با توجه به اینکه عمق آب روی سرریز بیشتر از 5 سانتیمتر بوده اثر کشش سطحی بروی مشخصات هیدرولیکی پرش تاثیری ندارد. با استفاده از نظریه باکینگهام پارامتر‌های بی بعد بصورت زیر انتخاب شد.

 

(2)

 

که در این رابطه  و  بترتیب عدد فرود و رینولدز قبل از پرش می‌باشد. با توجه به اینکه در این تحقیق مقادیر عدد رینولدز بین  91300 تا 256900 بوده بنابراین جریان کاملا آشفته می‌باشد و می‌توان از تاثیر عدد رینولدز چشم پوشی کرد. با توجه به اینکه نسبت  در کلیه آزمایش‌ها ثابت می‌باشد بنابراین، این نسبت را می‌توان از رابطه فوق حذف کرد. در نتیجه می‌توان رابطه نهایی بی بعد را برای نسبت عمق مزدوج، طول و افت انرژی پرش بصورت زیر خلاصه کرد:

 

(3)                                                              

(4)                                               

(5)                               

 

محدوده پارامتر‌های برداشت شده در آزمایشگاه در جدول 2 نشان داده شده است. پارامترهای عمق اولیه پرش هیدرولیکی  ، عمق ثانویه پرش هیدرولیکی و طول پرش هیدرولیکی  برداشت شد. لازم به توضیح است که با جمع آوری این مقادیر میزان سرعت جریان در ابتدا و انتها با استفاده از رابطه پیوستگی بدست آمده و سپس عدد فرود و پارامترها هیدرولیکی پرش محاسبه شد.

 

نتایج و بحث

جدول 3 نتایج بدست آمده برای هر دو حوضچه‌ واگرا-همگرا مدل 1 و مدل 2 را نشان می‌دهد. نتایج نشان می دهد که با کاهش دبی، مقادیر عمق جریان و همچنین طول پرش هیدرولیکی در هر دو حوضچه کاهش می‌یابد.

 

جدول 2- محدوده پارامتر‌های آزمایشگاهی

V1

(m/s)

y2

(cm)

y1

(cm)

Fr1

Q

(l/s)

Variables

3.6-4.7

14.2-30

1.8-5

6.05-10.5

23.5-66.1

Range of parameters

 

جدول 3-مقادیر برداشت شده برای حوضچه‌های واگرا- همگرا

 

NO

 

 

Q

(Lit/s)

 

Model 1

 

 Model 2

y1

(m)

y2

(m)

Lj

(m)

y1

(m)

y2

(m)

Lj

(m)

1

66.1

 

0.047

0.25

0.86

 

0.048

0.28

1

2

58.0

 

0.0433

0.23

0.8

 

0.047

0.25

0.85

3

56.0

 

0.0414

0.22

0.8

 

0.045

0.245

0.86

4

54.5

 

0.04

0.215

0.7

 

0.043

0.24

0.84

5

53.5

 

0.039

0.21

0.65

 

0.042

0.23

0.8

6

49.5

 

0.0365

0.205

0.65

 

0.04

0.22

0.75

7

45.0

 

0.032

0.197

0.6

 

0.038

0.205

0.72

8

41.0

 

0.0305

0.184

0.6

 

0.034

0.2

0.72

9

40.0

 

0.029

0.18

0.55

 

0.032

0.205

0.7

10

38.0

 

0.0283

0.172

0.55

 

0.031

0.2

0.67

11

32.0

 

0.026

0.16

0.55

 

0.029

0.175

0.6

12

29.5

 

0.023

0.155

0.5

 

0.025

0.17

0.6

13

25.5

 

0.02

0.148

0.5

 

0.022

0.166

0.57

14

23.5

 

0.018

0.142

0.5

 

0.02

0.16

0.57

 

جدول 4 مقادیر پارمترهای جریان از جمله نسبت عمق جریان ، نسبت افت انرژی و همچنین نسبت طول پرش هیدرولیکی   را نشان می‌دهد. همانگونه که در این جدول نشان داده شده است در حداکثر دبی(lit/s  )، نسبت عمق جریان و طول پرش در حوضچه واگرا-همگرا مدل 1 برابر 32/5 و 30/18 و مدل 2 بترتیب برابر با 83/5 و 83/20 بوده است. که نشان دهنده کاهش این پارامترها در مدل 1 نسبت به مدل 2 می‌باشد. از طرفی مقدار  برای حداکثر دبی، در حوضچه‌ واگرا-همگرا مدل 1 برابر 77/0 و در مدل 2 برابر 73/0بوده است که نشان دهنده افزایش افت در مدل 1 می‌باشد. در شکل 3 تغییرات نسبت عمق جریان بر حسب عدد فرود نشان داده شده است. نتایج نشان می‌دهد که با افزایش عدد فرود نسبت عمق در هر دو مدل افزایش می‌یابد. همچنین در تمامی مقادیر عدد فرود، نسبت عمق در حوضچه‌ واگرا-همگرا مدل 1 کمتر از حوضچه واگرا-همگرا مدل 2 می‌باشد. بطور مثال برای عدد فرود  این مقدار در مدل 1 برابر 31/5 و در مدل 2 برابر 88/5 می باشد.نتایج حاصله نشان دهنده کاهش نسبت عمق در حوضچه واگرا-همگرا مدل 1 می‌باشد.

 

جدول 4- پارامترهای مربوط به جریان در حوضچه‌های واگرا- همگرا

 

Q

(l/s)

 

 

Model 1

 

 Model 2

 

Fr1

(ΔE/E1)

(y2/y1)

(Lj/y1)

Fr1

(ΔE/E1)

(y2/y1)

(Lj/y1)

66.1

 

6.90

0.77

5.32

18.30

 

6.68

0.73

5.83

20.83

58.0

 

6.85

0.77

5.31

18.48

 

6.05

0.71

5.32

18.09

56.0

 

7.07

0.78

5.31

19.32

 

6.24

0.72

5.44

19.11

54.5

 

7.25

0.79

5.38

17.50

 

6.50

0.73

5.58

19.53

53.5

 

7.39

0.79

5.38

16.67

 

6.61

0.75

5.48

19.05

49.5

 

7.55

0.80

5.62

17.81

 

6.58

0.74

5.50

18.75

45.0

 

8.36

0.82

6.15

18.75

 

6.46

0.74

5.39

18.95

41.0

 

8.19

0.81

6.03

19.67

 

6.96

0.75

5.88

21.18

40.0

 

8.61

0.81

6.20

18.97

 

7.43

0.76

6.41

21.88

38.0

 

8.49

0.81

6.07

19.43

 

7.40

0.76

6.45

21.61

32.0

 

8.12

0.80

6.15

21.15

 

6.89

0.74

6.03

20.69

29.5

 

9.00

0.82

6.73

21.74

 

7.94

0.78

6.80

24.00

25.5

 

10.3

0.84

7.7

26.32

 

8.31

0.78

7.55

25.91

23.5

 

10.3

0.84

7.8

27.78

 

8.84

0.79

8.00

28.50

 

img0 

شکل 3) تغییرات نسبت عمق جریان بر حسب عدد فرود

در شکل 4 تغییرات نسبت طول پرش  بر حسب عدد فرود نشان داده شده است. نتایج نشان می‌دهد که در هر دو حوضچه با افزایش عدد فرود میزان طول پرش افزایش می‌یابد اما همواره طول پرش در حوضچه‌های واگرا-همگرا مدل 1 از حوضچه‌های واگرا- همگرا مدل 2 کمتر می‌باشد. بطور مثال برای عدد فرود  این مقدار در مدل 1 برابر 45/18 و در مدل 2 برابر 18/21 می‌باشد.  همانگونه که درجدول 4 نشان داده شده است، در حداکثر دبی، مقدار نسبی طول پرش  در مدل‌ 1 برابر با 30/18 و در مدل 2 برابر 83/20 می‌باشد. نتایج نشان می‌دهد که طول پرش در حوضچه واگرا- همگرا مدل 1 نسبت به حوضچه‌ واگرا-همگرا مدل 2 کاهش یافته است. با توجه به اینکه میزان افت انرژی یکی از پارامترهای مهم در طراحی حوضچه‌های آرامش می‌باشد لذا در شکل 5 تغییرات افت انرژی نسبت به عدد فرود نشان داده شده است. نتایج نشان می دهد که  با افزایش عدد فرود نسبت افت انرژی در هر دو حوضچه افزایش می‌یابد اما این مقدار در حوضچه‌های واگرا-همگرا مدل 1 نسبت به حوضچه‌های واگرا-همگرا مدل 2 بیشتر است. نتایج ارائه شده در جدول 4 نشان می‌دهد که در حداکثر دبی، میزان افت انرژی  در حوضچه واگرا-همگرا مدل 1 برابر 77/0 و در حوضچه مدل 2 برابر  73/0 بوده که نشان دهنده افزایش افت در مدل 1 می‌باشد. 

 

img0 

شکل 4) تغییرات نسبت طول پرش بر حسب عدد فرود

 

img0 

شکل 5) تغییرات نسبت افت انرژی بر حسب عدد فرود

شکل 6  مقایسه‌ بین نتایج این تحقیق و محققان پیشین را نشان می‌دهد. با توجه به شکل برای مقادیر مختلف عدد فرود، نسبت عمق جریان  در مدل 1 و 2 از تمامی مدل‌ها کمتر می‌باشد. بعبارت دیگر برای یک عدد فرود ثابت بیشترین مقدار  مربوط به پرش هیدرولیک کلاسیک و کمترین آن بترتیب متعلق به حوضچه آرامش واگرا-همگرا مدل 1 و مدل 2 می­باشد می‌توان نتیجه گرفت که حوضچه‌های واگرا-همگرا بهترین تبدیل برای پرش هیدرولیکی با حداقل عمق ثانویه می­باشد.

 

 

img0 

شکل 7- مقایسه نسبت عمق پرش در تحقیق حاضر و تحقیقات دیگران

 

در انتهای این تحقیق مقایسه‌ میزان در صد کاهش یا افزایش پارامترهای پرش نسبت به حوضچه کلاسیک (مولفه‌های پرش کلاسیک با استفاده از رابطه بلانگر و معادله انرژی محاسبه شده است) در جدول 5 ارائه شده است. در این جدول مقادیر مثبت و منفی به ترتیب نشان دهنده میزان افزایش و کاهش می‌باشد. نتایج نشان می‌دهد که در حداکثر دبی، مقادیر  و در حوضچه واگرا-همگرا مدل 1 به ترتیب به میزان 6/48 و 9/122 و در مدل 2 به ترتیب به میزان 5/35 و 8/95 درصد نسبت به حوضچه کلاسیک کاهش داشته‌اند. از طرفی مقدار افت انرژی در مدل‌ 1 و 2 به ترتیب  8/26 و 7/23 درصد نسبت به حوضچه‌های کلاسیک افزایش یافته است.

 

جدول 5- مقایسه عمق ثانویه، طول پرش و افت انرژی در تحقیق حاضر با پرش هیدرولیک کلاسیک

Q

(Lit/s)

Variation between divergence-convergence and classic basins

(model 1) (%)

 

Variation between divergence-convergence and classic basins

(model 2) (%)

 

y2/y1

Lj/y1

ΔE/E1

 

y2/y1

Lj/y1

ΔE/E1

66.1

-48.61

-122.9

26.89

 

-35.51

-95.79

23.79

58.0

-37.88

-100.4

30.54

 

-37.69

-104.7

24.96

56.0

-41.71

-98.55

30.35

 

-38.31

-100.7

24.55

54.5

-40.76

-120.5

30.89

 

-35.55

-97.58

25.86

53.5

-40.14

-130.8

31.68

 

-37.79

-101.9

27.22

49.5

-25.84

-98.63

34.89

 

-28.51

-88.65

30.26

45.0

-15.53

-90.18

36.26

 

-31.84

-88.20

29.49

41.0

-19.84

-85.14

35.18

 

-22.91

-71.99

30.05

40.0

-31.27

-123.3

30.78

 

-27.19

-93.64

25.23

38.0

-35.72

-121.3

30.10

 

-27.86

-99.00

24.35

32.0

-22.68

-81.97

32.64

 

-25.11

-86.05

26.79

29.5

-13.52

-80.06

33.69

 

-12.50

-63.10

29.69

25.5

-17.13

-82.85

28..88

 

-17.13

-76.43

23.78

23.5

-31.21

-103.4

23.25

 

-29.38

-98.26

17.80

 

برای محاسبه عمق ثانویه پرش با استفاده از روش تحلیلی می‌توان رابطه زیر را نوشت:

 

(6)                                                                                                                                                

 (7)                                                                                                                                                    

 

با توجه به معادله اندازه حرکت رابطه زیر را می‌توان معرفی نمود.

(8)                                                                                          

(9)                                              

(10)                                                                                 

 

که  نیروی مومنتوم قبل از حوضچه،  نیروی وارده از دیوار واگرا،  نیروی وارده از دیوار همگرا است که می‌توان بصورت زیر نوشت.

 

(11)                                                                                                                                 

(12)                                                                                  

(13)                                                                                

(14)

(15)

که α و β بترتیب زوایای واگرایی و همگرایی می‌باشند و  عمق جریان در محل تلاقی دیواره واگرا و همگرا و بخشی از  عمق ثانویه پرش می‌باشد که بصورت سعی و خطا محاسبه می‌شود.

 

نتیجه­گیری

برای تنظیم و مهار پرش هیدرولیکی در پایین‌دست سازه­های هیدرولیکی نظیر سرریز‌ها، تندآب‌ها، حوضچه‌های آرامش با اشکال مختلف مورد استفاده قرار می­گیرد. در این تحقیق تاثیر زاویه واگرائی و همگرائی بر خصوصیات پرش هیدرولیکی در حوضچه‌های آرامش واگرا-همگرا بررسی و نتایج آن با یکدیگر مقایسه شده است. آزمایش‌ها در کانالی به طول 11 متر، عرض 48/0 متر و عمق 7/0 متر انجام گرفت. جهت بررسی عملکرد زاویه واگرایی و همگرایی بر خصوصیات پرش هیدرولیکی، از دو نوع حوضچه به ترتیب با زاویه واگرایی و همگرایی  و  (مدل 1 و ) و  و  درجه (مدل 2 و ) استفاده شده است. نتایج نشان می‌دهد که در حداکثر دبی، نسبت عمق جریان  و نسبت طول جریان  در مدل 1 به ترتیب برابر 32/5 و 30/18 و در مدل 2 به ترتیب برابر با 83/5 و 83/20 می‌باشد. از طرفی نسبت افت انرژی در مدلهای 1 و 2 به ترتیب برابر با 77/0 و 73/0 می باشد.مقایسه نتایج بدست آمده با حوضچه های کلاسیک نشان می‌دهد که نسبت عمق جریان  و طول جریان  در مدل 1 به ترتیب به میزان 6/48 و 9/122 درصد، و در مدل 2 به ترتیب به میزان 5/35 و 8/95 درصد کاهش داشته‌ است. همچنین نسبت افت انرژی در مدل‌های‌ 1 و 2 به ترتیب 8/26 و 7/23 درصد نسبت به حوضچه‌های کلاسیک افزایش یافته است. با توجه به به اینکه نسبت عمق جریان و طول پرش در حوضچه‌های واگرا- همگرا مدل 1 و 2 نسبت به حوضچه‌های کلاسیک به طور چشمگیری کاهش و افت انرژی افزایش یافته است میتوان نتیجه گرفت که این نوع حوضچه ها موجب کاهش ابعاد و همچنین افزایش عملکرد حوضچه‌های آرامش می‌گردد. از طرفی افزایش نسبت زاویه واگرایی به همگرایی () موجب افزایش عملکرد حوضچه‌های واگرا-همگرا خواهد شد.

 

 

 تقدیر و تشکر

اين تحقیق در مرکز تحقیقات آب دانشگاه آزاد اسلامی واحد استهبان انجام شده است. نويسندگان این مقاله  از مسئولین مربوطه  جهت  مساعدت و استفاده تجهیزات آزمایشگاهی تشکر و قدرداني مي­نمايند.

 

منابع:

1.         عزیزی ح, ساجده, زاده ف, داود, عباسپور, اکرم, et al. بررسی آزمایشگاهی ساختار تلاطم پرش هیدرولیکی مستغرق در کانال واگرای تدریجی با بستر زبر. پژوهش آب ایران. 2022;16(3).

2.         دزفولی ن, بدست ک, عباس ا, مسجدی, حیدرنژاد, محمد, پورمحمدی. تعیین اثر زبری بر مشخصات پرش هیدرولیکی در حوضچه‌های آرامش واگرا. مجله علوم و فنون دریایی. 2022;21(1):81-91.

3. وند ز, سجادی, نژاد ح, مصیب, پری ا, امین س, بجستان ش. مقایسه آزمایشگاهی تاثیر زبری مصنوعی و شیب معکوس بستر در کنترل پرش هیدرولیکی. فصلنامه علمی مهندسی منابع آب. 2023;16(57):19-30.

34.        Chow. V.T. Open channel hydraulics. McGraw-Hill, New York, 1959.

45.        Daneshfaraz R, Sadeghı H, Joudı AR, Abraham J. Experimental investigation of hydraulic jump characteristics in contractions and expansions. Sigma Journal of Engineering and Natural Sciences. 2017;35(1):87-98.

56.        Daneshfaraz R, MajediAsl M, Mirzaee R, Ghaderi A. The S-jump’s characteristics in the rough sudden expanding stilling basin. AUT Journal of Civil Engineering. 2020;4(3):349-56.

67.        Daneshfaraz R, Najibi A. Experimental investigation of supercritical flow energy dissipation in sudden contraction with wall roughness. Journal of Hydraulics. 2021;16(4):79-92.

8.        Daneshfaraz R, MajediAsl M, Mirzaee R, Tayfur G. Hydraulic jump in a rough sudden symmetric expansion channel. AUT Journal of Civil Engineering. 2021;5(2):245-56.

9.        Daneshfaraz R, Sadeghi H, Ghaderi A, Abraham JP. Characteristics of hydraulic jump and energy dissipation in the downstream of stepped spillways with rough steps. Flow Measurement and Instrumentation. 2024;96:102506.

710.        Eshkou Z, Dehghani A, Ahmadi A. Forced hydraulic jump in a diverging stilling basin using angled baffle blocks. Journal of Irrigation and Drainage Engineering. 2018;144(8):06018004.

11.        Ghaderi A, Dasineh M, Aristodemo F, Ghahramanzadeh A. Characteristics of free and submerged hydraulic jumps over different macroroughnesses. Journal of Hydroinformatics. 2020;22(6):1554-72.

812.        Hassanpour N, Hosseinzadeh Dalir A, Farsadizadeh D, Gualtieri C. An experimental study of hydraulic jump in a gradually expanding rectangular stilling basin with roughened bed. Water. 2017;9(12):945.

913.        Hamidinejad AE, Heidarpour M, Ghadampour Z. Hydraulic Jump Control Using Stilling Basin with Abruptly Expanding and Negative Step. Iranian Journal of Science and Technology, Transactions of Civil Engineering. 2023;47(6):3885-94.

1014.        Neisi K, Shafai Bejestan M, Ghomshi M, Kashefipoor S. Investigation of Hydraulic Jump Characteristics at Roughened Bed of Sudden Expansion Stilling Basin. Irrigation Sciences and Engineering. 2014;37(2):83-93.

1115.        Torkamanzad N, Hosseinzadeh Dalir A, Salmasi F, Abbaspour A. Hydraulic jump below abrupt asymmetric expanding stilling basin on rough bed. Water. 2019;11(9):1756.

1216.        Varaki ME, Kasi A, Farhoudi J, Sen D. Hydraulic jump in a diverging channel with an adverse slope. Iranian Journal of Science and Technology Transactions of Civil Engineering. 2014;38(C1):111.

 

 

 


مقالات مرتبط

حقوق این وب‌سایت متعلق به سامانه مدیریت نشریات دانشگاه آزاد اسلامی است.
حق نشر © 1403-1400

صفحه اصلی| عضویت/ ورود| درباره سند| تماس با ما|
[English] [العربية] [en] [ar]