Numerical study of the dam-break waves and Favre waves down sloped wet rigid-bed at laboratory scale

Numerical study of the dam-break waves and Favre waves down sloped wet rigid-bed at laboratory scale

WenjunLiua  BoWangb  YakunGuoc

a State Key Laboratory of Hydraulics and Mountain River Engineering, College of Water Resource and Hydropower, Sichuan University, Chengdu, 610065, China
State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, College Of Water Resource and Hydropower, Chengdu, 610065, China
faculty of Engineering & Informatics, University of Bradford, BD7 1DP, UK

Abstract

The bed slope and the tailwater depth are two important ones among the factors that affect the propagation of the dam-break flood and Favre waves. Most previous studies have only focused on the macroscopic characteristics of the dam-break flows or Favre waves under the condition of horizontal bed, rather than the internal movement characteristics in sloped channel. The present study applies two numerical models, namely, large eddy simulation (LES) and shallow water equations (SWEs) models embedded in the CFD software package FLOW-3D to analyze the internal movement characteristics of the dam-break flows and Favre waves, such as water level, the velocity distribution, the fluid particles acceleration and the bed shear stress, under the different bed slopes and water depth ratios. The results under the conditions considered in this study show that there is a flow state transition in the flow evolution for the steep bed slope even in water depth ratio α = 0.1 (α is the ratio of the tailwater depth to the reservoir water depth). The flow state transition shows that the wavefront changes from a breaking state to undular. Such flow transition is not observed for the horizontal slope and mild bed slope. The existence of the Favre waves leads to a significant increase of the vertical velocity and the vertical acceleration. In this situation, the SWEs model has poor prediction. Analysis reveals that the variation of the maximum bed shear stress is affected by both the bed slope and tailwater depth. Under the same bed slope (e.g., S0 = 0.02), the maximum bed shear stress position develops downstream of the dam when α = 0.1, while it develops towards the end of the reservoir when α = 0.7. For the same water depth ratio (e.g., α = 0.7), the maximum bed shear stress position always locates within the reservoir at S0 = 0.02, while it appears in the downstream of the dam for S0 = 0 and 0.003 after the flow evolves for a while. The comparison between the numerical simulation and experimental measurements shows that the LES model can predict the internal movement characteristics with satisfactory accuracy. This study improves the understanding of the effect of both the bed slope and the tailwater depth on the internal movement characteristics of the dam-break flows and Favre waves, which also provides a valuable reference for determining the flood embankment height and designing the channel bed anti-scouring facility.

댐붕괴 홍수와 파브르 파도의 전파에 영향을 미치는 요인 중 하상경사와 후미수심은 두 가지 중요한 요소이다. 대부분의 선행 연구들은 경사 수로에서의 내부 이동 특성보다는 수평층 조건에서 댐파괴류나 Favre파동의 거시적 특성에만 초점을 맞추었다.

본 연구에서는 CFD 소프트웨어 패키지 FLOW-3D에 내장된 LES(Large Eddy Simulation) 및 SWE(Shallow Water Equation) 모델의 두 가지 수치 모델을 적용하여 댐-파괴 흐름 및 Favre 파도의 내부 이동 특성을 분석합니다.

수위, 속도 분포, 유체 입자 가속도 및 층 전단 응력, 다양한 층 경사 및 수심 비율로. 본 연구에서 고려한 조건하의 결과는 수심비 α = 0.1(α는 저수지 수심에 대한 tailwater 깊이의 비율)에서도 급경사면에 대한 유동상태 전이가 있음을 보여준다. 유동 상태 전이는 파면이 파단 상태에서 비정형으로 변하는 것을 보여줍니다.

수평 경사와 완만한 바닥 경사에서는 이러한 흐름 전이가 관찰되지 않습니다. Favre 파의 존재는 수직 속도와 수직 가속도의 상당한 증가로 이어집니다. 이 상황에서 SWE 모델은 예측이 좋지 않습니다.

분석에 따르면 최대 바닥 전단 응력의 변화는 바닥 경사와 꼬리 수심 모두에 영향을 받습니다. 동일한 바닥 경사(예: S0 = 0.02)에서 최대 바닥 전단 응력 위치는 α = 0.1일 때 댐의 하류에서 발생하고 α = 0.7일 때 저수지의 끝쪽으로 발생합니다.

동일한 수심비(예: α = 0.7)에 대해 최대 바닥 전단 응력 위치는 항상 S0 = 0.02에서 저수지 내에 위치하는 반면, S0 = 0 및 0.003에 대해 흐름이 진화한 후 댐 하류에 나타납니다. 수치적 시뮬레이션과 실험적 측정을 비교한 결과 LES 모델이 내부 움직임 특성을 만족스러운 정확도로 예측할 수 있음을 알 수 있습니다.

본 연구는 댐 파절류 및 Favre파의 내부 이동 특성에 대한 하상 경사 및 후미 수심의 영향에 대한 이해를 향상 시키며, 이는 또한 제방 높이를 결정하고 수로 저반위 설계를 위한 귀중한 참고자료를 제공한다.

Keywords

Figure Numerical study of the dam-break waves and Favre waves down sloped wet rigid-bed at laboratory scale
Figure Numerical study of the dam-break waves and Favre waves down sloped wet rigid-bed at laboratory scale

Dam-break flow, Bed slope, Wet bed, Velocity profile, Bed shear stress, Large eddy simulation
댐파괴유동, 하상경사, 습상, 유속분포, 하상전단응력, 대와류 시뮬레이션

Fig. 3. Breakwaters model in Flow-3D with meshing geometry and boundary (a) circular slots (b) square slots.

Study of Unconventional Alternatives to Vertical Breakwater

수직 방파제에 대한 비전통적 대안 연구

Karim Badr Hussein and Mohamed Ibrahim
Lecturer of Irrigation and Hydraulics, Faculty of Engineering, Al-Azhar University
Corresponding author E-mail: badrkarim713@yahoo.com

Abstract

방파제의 주요 목적은 항만 내부의 안정을 유지하여 선박의 안전과 운영의 용이성을 달성하는데 도움이 되기 때문에 강한 파도와 폭풍으로부터 항만, 해변 또는 해변 시설을 보호하는 것입니다.

이 연구는 수직 방파제에 대한 비전통적인 대안을 연구하는 것을 목표로 합니다. 이 연구에서는 유체역학적 성능의 연구 및 평가를 위해 구현된 수직파 장벽의 두 가지 다른 모델을 선택했습니다.

첫 번째 모델은 원형 슬롯이 있는 수직 벽이고 두 번째 모델은 사각형 슬롯이 있는 수직 벽입니다. 두 모델을 비교한 결과 정사각형 슬롯은 원형 슬롯보다 파동의 전송을 5~20% 감소시키는 것으로 나타났습니다.

두 개의 원형 홈이 있는 벽을 사용하면 단일 벽에 비해 파동 전송이 최대 30% 감소하고 파동 에너지 분산이 최대 40% 증가합니다. 상대 길이(h/L)가 증가함에 따라 수평파력이 증가합니다.

다공성 = 0.25에서 상대파력(F/Fo)은 다공성 = 0.50에서보다 10~30% 더 컸습니다. 개구부에서 파동 속도가 높고 파동 에너지 소산 계수도 높습니다. 파동 진폭이 클수록 파동 에너지 소산 계수가 커집니다.

Key words: Coastal, Breakwater, FLOW-3D, Numerical Models, Energy Dissipation, Vertical Wall.

Introduction

모든 국가에서 해안 지역은 가장 중요하고 중요한 지역 중 하나입니다. 연안지역과 항만은 대외무역 촉진, 연안관광 개발 및 활성화 등 다양한 분야에 기여하고 있어 경제적 파급효과가 매우 크며, 일자리 창출은 물론 도시근린 정착 및 안정에도 기여한다. 젊은이들에게 강력한 수익을 제공하는 가능성과 어항을 건설하여 어획량을 늘리는 것입니다. [1].

그러나 해안선 부근의 파도, 바람, 조수, 조류 등의 자연 현상은 해변과 해안 지역의 안정성에 영향을 미칩니다. 따라서 연안 보전 서비스는 연안 환경의 균형을 유지하고 보존하는 데 중요한 역할을 합니다. 거센 파도로부터 항구와 해변 시설을 보호하는 방파제 방파제. 방파제는 선박이 안전하게 정박할 수 있는 조용한 지역을 제공하고 건설 및 석유 및 광물 발견 동안 임시 보호를 제공합니다.

파도는 방파제에 부딪힐 때 많은 에너지를 잃습니다. 방파제는 눈에 보이거나 떠 있거나 수중일 수 있으며 다양한 크기, 재료 및 출력 표준이 있습니다[11]. 전통적인 장벽 또는 눈에 보이는 격벽은 매우 효율적이지만 해변의 미적 비전을 가립니다. 많은 건축 자재가 필요하고 건설 비용이 증가합니다[9].

이에 반해 부유방벽은 자재가 필요없고 공사비가 저렴하지만 그 효과는 제한적입니다. 결과적으로 수중 파티션은 이러한 종류의 단점을 방지하기 때문에 더 나은 옵션 중 하나로 간주됩니다.

수중 방벽은 가장 중요한 해변 방어 시설 중 하나이며, 수중 방벽의 장점 중 하나는 투명 방벽에 비해 건설 비용이 비교적 저렴하고 물이 앞에서 뒤로 흐를 수 있다는 것입니다[3].

멤브레인 아래에서 물이 재생됩니다. 또한 바다의 미적 이미지를 왜곡하지 않고 조망을 방해하지 않아 인근 해변에 미치는 영향도 미미하다[18]. 반면에 잠긴 방파제는 건설 후 가라앉으면서 파도 에너지를 분산시키고 해안선을 방어하는 효과를 잃습니다. 장벽의 품질은 높은 수위의 영향도 받습니다.

결과적으로 해안 보호의 가장 중요한 측면 중 하나는 수중 방파제의 효율성을 향상시키는 것입니다. 수직 방파제 이러한 유형의 방파제는 바다를 향한 수직면이 있는 설비입니다[10]. 이러한 장벽은 파도 에너지의 일부가 해안이나 보호할 수역에 도달하는 것을 방지하여 파도를 진정시키는 역할을 합니다[16].

수직 방파제는 블록, 케이슨, 시트 파일 또는 셀룰러로 구성될 수 있습니다. 이 연구는 정사각형 및 원형 구멍이 있는 천공된 수직 방파제의 유체역학적 성능에 대한 연구를 제시하는 것을 목적으로 합니다.

이 논문은 또한 제안된 모델의 유체역학적 효율뿐만 아니라 이 분야의 유사한 연구와 비교되었습니다. 이것은 다음 헤드라인으로 이 백서에 나와 있습니다.

 Materials and methods.
 Results and discussion.
 Conclusions and recommendations.

Fig. 1. The open channel
Fig. 1. The open channel
Fig. 2. Breakwaters model (a) perforated wall with circular slots and (b) perforated wall with square slots.
Fig. 2. Breakwaters model (a) perforated wall with circular slots and (b) perforated wall with square slots.
Fig. 3. Breakwaters model in Flow-3D with meshing geometry and boundary (a) circular slots (b) square slots.
Fig. 3. Breakwaters model in Flow-3D with meshing geometry and boundary (a) circular slots (b) square slots.
Fig. 4. Details and dimensions of proposed breakwater
Fig. 4. Details and dimensions of proposed breakwater
Fig 5 .Wave profiles using (Flow-3D) at wave period (T) = 1.2 sec for perforated walls with circular slots at behind model (Ht).
Fig 5 .Wave profiles using (Flow-3D) at wave period (T) = 1.2 sec for perforated walls with circular slots at behind model (Ht).
Fig. 11. Velocity distribution through slots at (a) quarter wave period, (b) half wave period and (c) three quarters wave period.
Fig. 11. Velocity distribution through slots at (a) quarter wave period, (b) half wave period and (c) three quarters wave period.
Fig. 13. Velocity vectors at front, between and behind barriers.
Fig. 13. Velocity vectors at front, between and behind barriers.

Conclusion & Recommendations

얻어진 결과에 대한 이전 분석을 바탕으로 도달한 결론은 다음과 같습니다.
 결과와 연구에 따르면 FLOW-3D는 수직으로 구멍이 뚫린 벽이 있는 선형 파동과 파동의 관계를 설명하는 강력한 능력을 가지고 있습니다. 또한 실험실 데이터 및 반분석 결과의 가장 중요한 측면을 복제할 수 있습니다. FLOW-3D에 의해 생성된 수치적 결과는 훌륭합니다.
 사각슬롯은 원형슬롯에 비해 파동의 투과율이 5:20% 감소합니다.
 한 쌍의 원형 슬롯 벽을 사용하면 단일 벽에 비해 파동 투과율이 최대 30% 감소하고 파동 에너지 분산이 최대 40% 증가합니다.
 수평파력은 상대길이(h/L)가 증가할수록 증가한다. 다공성 = 0.25에서 상대파력(F/Fo)은 다공성 = 0.50에서보다 10~30% 더 높았다.
 파도가 원 모양으로 움직이고 큰 원이 위쪽에 있었다가 점차 아래쪽으로 내려갑니다.  개구부에서 파동 속도가 높았고 파동 에너지 소산 계수도 높았습니다. 파동 진폭이 높을수록 파동 에너지 소산 계수가 높아집니다.

REFERENCES

[1] Bahaa Elsharnouby and Mohamed, E. (2012). “Study of environment
friendly porous suspended breakwater for the Egyptian Northwestern
Coast” J. of Ocean Engineering, Vol. 48, 47-58.
[2] Huang Z. (2007) “Wave interaction with one or two rows of closely
spaced rectangular cylinders” J. Ocean Eng Vol. 34,1584–1591.
[3] Huang, C. J.; Chang, H. H.; and Hwung, H. H., 2003. “Structural
permeability effects on the interaction of a solitary wave and a
submerged breakwater,” Coastal Engineering. Vol. 49, pp. 1-24.
[4] Hsu, H-H. & Wu, Y-C., 1999. “Numerical solution for the second-order
wave interaction with porous structures.” International Journal for
Numerical Methods in Fluids, Vol. 29 Issue 3, pp. 265-288.
[5] Isaacson, M., Baldwim, J., Premasiro, S. and Yang, G., (1999) “Wave
interaction with double slotted barriers.” J. Applied Ocean Research,
Vol. 21, No. 2, pp. 81-91.
[6] Isaacson, M., Premasiro, S. and Yang, G. (1998) “Wave Interaction with
Vertical Slotted Barrier” J. Waterway, Port, Coastal and Ocean Eng.,
ASCE, Vol. 124, No. 3.
[7] Ji, C.H. and Suh, K.D. (2010) “Wave interactions with multiple-row
curtainwall-pile breakwaters” J.Coastal Engineering vol. 57 issue 5, p.
500-512.
[8] Koraim, A. S., Iskander, M.M. and Elsayed, W. R. (2014)
“Hydrodynamic performance of double rows of piles suspending
horizontal c shaped bars” J. Coastal Engineering, Vol. 84, P. 81-96.
[9] Koraim, A. S., Iskander, M.M. and Elsayed, W. R. (2013)
“Hydrodynamic performance of double rows of piles suspending
horizontal c shaped bars” J. Coastal Engineering, Vol. 84, P. 81-96.
[10] Koraim, A. S. and Salem, T. N. (2012) “The hydrodynamic
characteristics of a single suspended row of half pipes under regular
waves” J. Ocean Engineering, Vol. 50, P. 1-9.
[11] Laju, K., Sundar, V. & Sundaravadivelu, R., 2011. “Hydrodynamic
characteristics of pile supported skirt breakwater models.” Journal of
Ocean Re, 33,12-22.
[12] Lin, P.; and Karunarathna, .S.A., 2007. “Numerical study of solitary
wave interaction with porous breakwaters,” J. of waterway, port,
coastal and ocean engineering. , pp. 352-363.
[13] Moh. Ibrahim (2017) “Linear Wave Interaction with Permeable
Breakwaters” A Thesis Submitted for Partial Fulfillment of Doctor of
Philosophy Degree in Civil Eng., al-Azhar University.
[14] Mansard, E .P. D. & Funke, E. R., 1980. “The measurement of incident
and reflected spectra using a least squares method.” In Proc. 17th
Coastal Eng. Conf., Sydney, Australia, pp 159-174.
[15] Nadji Chioukh et al (2017) “Reflection and Transmission of Regular
Waves from/Through Single and Double Perforated Thin Walls”
China Ocean Eng., 2017, Vol. 31, No. 4, P. 466–475.
[16] Rageh, O., Koraim, A. (2010b). “Hydraulic performance of vertical
walls with horizontal slots used as breakwater”. J.Coastal Engineering,
Vol. 57, 745–746. 12.
[17] Suh KD, Jung HY and Pyun CK (2007) “Wave reflection and
transmission by curtain wall–pile breakwaters using circular piles”. J.
Ocean Eng,Vol. 34(14–15), 2100–2106.
[18] Suh, K. D., Shin, S. & Cox, D. T., 2006. “Hydrodynamic
characteristics of Pile-Supported vertical wall breakwaters.” J. of
Waterways, Port, Coastal and Ocean Engineering, Vol.132, No.2,
pp.83-96.

Figure 1. The push barge model in 1:20 geometrical scale during field experiments.

Experimental Method for the Measurements and Numerical Investigations of Force Generated on the Rotating Cylinder under Water Flow

by Teresa Abramowicz-Gerigk 1,*,Zbigniew Burciu 1,Jacek Jachowski 1,Oskar Kreft 2,Dawid Majewski 3,Barbara Stachurska 3,Wojciech Sulisz 3 andPiotr Szmytkiewicz 3

1Faculty of Navigation, Gdynia Maritime University, 81-225 Gdynia, Poland
2AREX Ltd., 81-212 Gdynia, Poland
3Institute of Hydro-Engineering of Polish Academy of Sciences, 80-328 Gdansk, Poland
*Author to whom correspondence should be addressed.
Academic Editor: Remco J. WiegerinkSensors202121(6), 2216; https://doi.org/10.3390/s21062216
Received: 20 January 2021 / Revised: 9 March 2021 / Accepted: 18 March 2021 / Published: 22 March 2021(This article belongs to the Special Issue Sensing in Flow Analysis)

Abstract

본 논문은 자유 표면 효과를 포함한 균일한 흐름 하에서 회전하는 실린더 (로터)에 발생하는 유체 역학적 힘의 실험 테스트 설정 및 측정 방법을 제시합니다. 실험 테스트 설정은 고급 유량 생성 및 측정 시스템을 갖춘 수로 탱크에 설치된 고유 한 구조였습니다.

테스트 설정은 로터 드라이브가 있는 베어링 장착 플랫폼과 유체 역학적 힘을 측정하는 센서로 구성되었습니다. 낮은 길이 대 직경 비율 실린더는 얕은 흘수 강 바지선의 선수 로터 방향타 모델로 선택되었습니다. 로터 역학은 최대 550rpm의 회전 속도와 최대 0.85m / s의 수류 속도에 대해 테스트되었습니다.

실린더의 낮은 종횡비와 자유 표면 효과는 생성 된 유체 역학적 힘에 영향을 미치는 현상에 상당한 영향을 미쳤습니다. 회전자 길이 대 직경 비율, 회전 속도 대 유속 비율 및 양력에 대한 레이놀즈 수의 영향을 분석했습니다. 실험 결과에 대한 계산 모델의 유효성이 표시됩니다. 결과는 시뮬레이션 및 실험에 대한 결과의 유사한 경향을 보여줍니다.

The paper presents the experimental test setup and measurement method of hydrodynamic force generated on the rotating cylinder (rotor) under uniform flow including the free surface effect. The experimental test setup was a unique construction installed in the flume tank equipped with advanced flow generating and measuring systems.

The test setup consisted of a bearing mounted platform with rotor drive and sensors measuring the hydrodynamic force. The low length to diameter ratio cylinders were selected as models of bow rotor rudders of a shallow draft river barge. The rotor dynamics was tested for the rotational speeds up to 550 rpm and water current velocity up to 0.85 m/s. The low aspect ratio of the cylinder and free surface effect had significant impacts on the phenomena influencing the generated hydrodynamic force. The effects of the rotor length to diameter ratio, rotational velocity to flow velocity ratio, and the Reynolds number on the lift force were analyzed. The validation of the computational model against experimental results is presented. The results show a similar trend of results for the simulation and experiment.

Keywords: rotating cylinderforce sensor with built-in amplifierstrain gauge sensorCFD analysis

Figure 1. The push barge model in 1:20 geometrical scale during field experiments.
Figure 1. The push barge model in 1:20 geometrical scale during field experiments.
Figure 2. Scheme of the measurement area.
Figure 2. Scheme of the measurement area.
Figure 3. The force measuring part of the experimental test setup: (a) side view: 1—bearing-mounted platform, 2—drive system, 3—cylinder, 4—support frame, 5—force sensors, and 6—adjusting screw; (b) top view.
Figure 3. The force measuring part of the experimental test setup: (a) side view: 1—bearing-mounted platform, 2—drive system, 3—cylinder, 4—support frame, 5—force sensors, and 6—adjusting screw; (b) top view.
Figure 4. Location of the rotor, rotor drive, and supporting frame in the wave flume.
Figure 4. Location of the rotor, rotor drive, and supporting frame in the wave flume.
Figure 5. Lift force obtained from the measurements in the wave flume for different flow velocities and cylinder diameters.
Figure 5. Lift force obtained from the measurements in the wave flume for different flow velocities and cylinder diameters.
Figure 6. Variation of the lift coefficient with rotation rate for various free stream velocities and various cylinder diameters—experimental results.
Figure 6. Variation of the lift coefficient with rotation rate for various free stream velocities and various cylinder diameters—experimental results.
Figure 7. Boundary conditions for rotor-generated flow field simulation—computing domain with free surface level.
Figure 7. Boundary conditions for rotor-generated flow field simulation—computing domain with free surface level.
Figure 8. General view and the close-up of the rotor wall sector applied for the rotor simulation.
Figure 8. General view and the close-up of the rotor wall sector applied for the rotor simulation.
Figure 9. Structured mesh used in FLOW-3D and the FAVORTM technique—the original shape of the rotor and the shape of the object after FAVOR discretization technique for 3 mesh densities.
Figure 9. Structured mesh used in FLOW-3D and the FAVORTM technique—the original shape of the rotor and the shape of the object after FAVOR discretization technique for 3 mesh densities.
Figure 10. Parameter y+ for the studied turbulence models and meshes.
Figure 10. Parameter y+ for the studied turbulence models and meshes.
Figure 11. Results of numerical computations in time for the cylinder with D2 diameter at 500 rpm rotational speed and current speed V = 0.82 m/s using LES model in dependence of mesh density: (a) FX and (b) FY
Figure 11. Results of numerical computations in time for the cylinder with D2 diameter at 500 rpm rotational speed and current speed V = 0.82 m/s using LES model in dependence of mesh density: (a) FX and (b) FY
Figure 12. Results of 3D flow simulation for V = 0.40 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 12. Results of 3D flow simulation for V = 0.40 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 13. Results of 3D flow simulation for V = 0.50 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 13. Results of 3D flow simulation for V = 0.50 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 14. Results of 3D flow simulation for V = 0.82 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 14. Results of 3D flow simulation for V = 0.82 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.
Figure 15. Flow chart of validation of the computational model against experimental results.
Figure 15. Flow chart of validation of the computational model against experimental results.
Figure 16. Measured (EXP) and computed (CFD) lift force values.
Figure 16. Measured (EXP) and computed (CFD) lift force values.

결론

결론은 다음과 같습니다.
계산 결과가 일반적으로 실험 데이터와 일치하는 경우 계산 결과는 검증 된 것으로 간주되며 추가 예측에 사용할 수 있습니다. 검증 실험을 통해 메쉬 밀도와 난류 모델을 결정할 수있었습니다.
작은 전류 속도 0.4m / s 및 0.5m / s에서 직경 D3의 로터에 대해 계산 된 양력 값은 회전 속도가 200rpm 이상일 때의 실험 값과 달랐습니다. 그 이유는 실험 중에 관찰 된 강한 진동과 수치 시뮬레이션에서 모델링되지 않은 유동 분리 때문이었습니다.
D2 직경을 가진 로터의 경우 작은 rpm에서 양력의 반대 부호가 관찰되었습니다. 이 현상은 시뮬레이션 중에 관찰되지 않았습니다.
제시된 실험 테스트 설정은 드라이브,지지 구조물 및 측정 장치에 손상을 주지 않고 진동을 포함한 모든 현상을 관찰 할 수 있도록 구성되었습니다. Wang et al. [14]는 동일한 α 값에서 실린더 종횡비가 증가함에 따라 와류 유발 진동이 증가하는 것을 관찰했습니다.
실험의 원활한 진행은 장치 손상 가능성과 함께 약 4의 α에 영향을 미쳤습니다. 본 연구에서는 α = 4.8에서 시작하는 가장 큰 직경의 실린더에서 가장 강한 진동이 관찰되었습니다.
제시된 연구는 로터 생성 흐름의 능동적 제어에 대한 추가 연구의 첫 번째 부분으로 유체 역학적 힘의 신뢰할 수 있는 실험적 예측 방법을 설명했습니다 [22]. , 바람, 파도 [23].
논문의 참신함은 저상 실린더에 대해 회 전자에서 생성 된 유체 역학적 힘을 모델링 할 수있는 가능성에 대한 조사입니다.
이 방법의 주요 장점은 자유 표면 효과 및 유동 유도 회 전자 진동과 관련된 현상을 포함하여 회 전자 생성 유동장 및 유체 역학적 힘을 관찰 할 수 있다는 것입니다. 제안 된 테스트 설정 구성은 유체 역학적 힘의 매개 변수 연구, 스케일 효과 조사 및 낮은 전류 속도와 큰 회전 속도에서 큰 불일치가 확인 된 CFD 시뮬레이션 모델의 검증에 사용될 것입니다.

References

  1. Abramowicz-Gerigk, T.; Burciu, Z.; Jachowski, J. An Innovative Steering System for a River Push Barge Operated in Environmentally Sensitive Areas. Pol. Marit. Res. 201724, 27–34. [Google Scholar] [CrossRef]
  2. Abramowicz-Gerigk, T.; Burciu, Z.; Krata, P.; Jachowski, J. Steering system for a waterborne inland unit. Patent 420664, 2017. [Google Scholar]
  3. Abramowicz-Gerigk, T.; Burciu, Z.; Jachowski, J. Parametric study on the flow field generated by river barge bow steering systems. Sci. J. Marit. Univ. Szczec. 201960, 9–17. [Google Scholar]
  4. Gerigk, M.; Wójtowicz, S. An Integrated Model of Motion, Steering, Positioning and Stabilization of an Unmanned Autonomous Maritime Vehicle. TransnavInt. J. Mar. Navig. Saf. Sea Transp. 20159, 591–596. [Google Scholar] [CrossRef]
  5. Thouault, N.; Breitsamter, C.; Adams, N.A.; Seifert, J.; Badalamenti, C.; Prince, S.A. Numerical Analysis of a Rotating Cylinder with Spanwise Disks. AIAA J. 201250, 271–283. [Google Scholar] [CrossRef]
  6. Badr, H.M.; Coutanceau, M.; Dennis, S.C.R.; Menard, C. Unsteady flow past a rotating circular cylinder at Reynolds numbers 10 3 and 10 4. J. Fluid Mech. 1990220, 459. [Google Scholar] [CrossRef]
  7. Karabelas, S.; Koumroglou, B.; Argyropoulos, C.; Markatos, N. High Reynolds number turbulent flow past a rotating cylinder. Appl. Math. Model. 201236, 379–398. [Google Scholar] [CrossRef]
  8. Chen, W.; Rheem, C.-K. Experimental investigation of rotating cylinders in flow. J. Mar. Sci. Technol. 201924, 111–122. [Google Scholar] [CrossRef]
  9. Zhou, B.; Wang, X.; Guo, W.; Gho, W.M.; Tan, S.K. Experimental study on flow past a circular cylinder with rough surface. Ocean Eng. 2015109, 7–13. [Google Scholar] [CrossRef]
  10. Tokumaru, P.T.; Dimotakis, P.E. The lift of a cylinder executing rotary motions in a uniform flow. J. Fluid Mech. 1993255, 1–10. [Google Scholar] [CrossRef]
  11. Wong, K.W.L.; Zhao, J.; Jacono, D.L.; Thompson, M.C.; Sheridan, J. Experimental investigation of flow-induced vibration of a rotating circular cylinder. J. Fluid Mech. 2017829, 486–511. [Google Scholar] [CrossRef]
  12. Bourguet, R.; Jacono, D.L. Flow-induced vibrations of a rotating cylinder. J. Fluid Mech. 2014740, 342–380. [Google Scholar] [CrossRef]
  13. Carstensen, S.; Mandviwalla, X.; Vita, L.; Schmidt, P. Lift of a Rotating Circular Cylinder in Unsteady Flows. J. Ocean Wind Energy 20141, 41–49. Available online: http://www.isope.org/publications (accessed on 15 January 2021).
  14. Wang, W.; Wang, Y.; Zhao, D.; Pang, Y.; Guo, C.; Wang, Y. Numerical and Experimental Analysis of the Hydrodynamic Performance of a Three-Dimensional Finite-Length Rotating Cylinder. J. Mar. Sci. Appl. 202019, 388–397. [Google Scholar] [CrossRef]
  15. Mobini, K.; Niazi, M. Simulation of unsteady flow around a rotating circular cylinder at various Reynolds numbers. JMEUT 201746, 249–257. Available online: https://www.researchgate.net/publication/323447030_Simulation_of_Unsteady_Flow_Around_a_Rotating_Circular_Cylinder_at_Various_Reynolds_Numbers (accessed on 15 January 2021).
  16. Babarit, A.; Delvoye, S.; Arnal, V.; Davoust, L.; Wackers, J. Wave and Current Generation in Wave Flumes Using Axial-Flow Pumps. In Proceedings of the 36th International Conference on Ocean, Offshore and Artic Engineering (OMAE2017), Trondheim, Norway, 25–30 June 2017; pp. 1–10. [Google Scholar] [CrossRef]
  17. Nortek Manuals. The Comprehensive Manual for Velocimeters. 2018. Available online: https://support.nortekgroup.com/hc/en-us/articles/360029839351-The-Comprehensive-Manual-Velocimeters (accessed on 15 January 2021).
  18. Stachurska, B.; Majewski, D. Propagation of Surface waves under currents—Analysis of measurements in wave flume of IBW PAN. IMiG 20144, 280–290. [Google Scholar]
  19. Lohrmann, A.; Cabrera, R.; Kraus, N. Acoustic-Doppler Velocimeter (ADV) for laboratory use. In Fundamentals and Advancements in Hydraulic Measuremensts and Experimentation; Buffalo: New York, NY, USA, 1994. [Google Scholar]
  20. Stachurska, B.; Majewski, D. Experimental Measurements of Current Velocity in Wave Flume of IBW PAN; Internal Report; Institute of Hydro-Engineering of Polish Academy of Sciences: Gdańsk, Poland, 2013. (In Polish) [Google Scholar]
  21. FLOW-3D. Available online: https://www.flow3d.com/ (accessed on 15 January 2021).
  22. He, J.W.; Glowinski, R.; Metcalfe, R.; Nordlander, A.; Periaux, J. Active control and drag optimization for flow past a circular cylinder: Oscillatory cylinder rotation. J. Comput. Phys. 2000163, 83–117. [Google Scholar] [CrossRef]
  23. Lebkowski, A. Analysis of the Use of Electric Drive Systems for Crew Transfer Vessels Servicing Offshore Wind Farms. Energies 202013, 1466. [Google Scholar] [CrossRef]
Figure 3. Flow velocity on seawall in A2-3 modeling.

Modeling of the Changes in Flow Velocity on Seawalls under Different Conditions Using FLOW-3D Software

Open Journal of Marine Science
Vol.06 No.02(2016), Article ID:65874,6 pages
10.4236/ojms.2016.62026

FLOW-3D 소프트웨어를 사용하여 다양한 조건에서 Seawalls의 흐름 속도 변경 모델링

Maryam Deilami-Tarifi1, Mehdi Behdarvandi-Askar2*, Vahid Chegini3, Sadegh Haghighi-Pour4
1Department of Coastal Engineering, Khorramshahr University of Marine Science and Technology, Khorramshahr, Iran

2Department of Marine Structures, Khorramshahr University of Marine Science and Technology, Khorramshahr, Iran
3Iran National Center for Oceanography and Atmospheric Sciences, Tehran, Iran
4Department of Civil Engineering, Excellence in Education Center of Jihad University of Khuzestan, Ahvaz, Iran
Copyright © 2016 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

ABSTRACT

방파벽은 파도힘의 수준을 감소시키고 다른 구조물로부터 보호하기 위해 건설되는 보호 구조물 중 하나입니다. 이와 관련하여 이러한 구조에 대한 보다 정확한 조사는 다른 관점에서 매우 중요합니다. 이 연구는 다른 레이아웃과 경사면에서 장애물을 고려하여 방파제 크라운의 속도 변화를 조사합니다. FLOW-3D는 모델링을 위한 이 연구에서 사용되었습니다. 모델링의 결과는 장애물의 존재가 방파벽의 크라운의 유량을 줄이는 결정적인 역할을 한다는 것을 보여줍니다. 또한, 예상대로, 상류 방파의 경사계는 벽의 가장 낮은 속도가 D-상태 레이아웃과 45°의 경사에서 발생하므로 이 속도를 줄이는 데 매우 결정적입니다.

Keywords: 플로우 속도, 방파제 크라운, 모델링, Flow Velocity, Seawall Crown, Modeling, FLOW-3D

1. 소개

방파벽은 파도의 속도를 감소시키고 다른 구조물을 보호하기 위해 건설되는 보호 구조물 중 하나입니다. 등대는 일반적으로 방파벽에 의해 보호되는 구조 중 하나입니다. 따라서, 방파성상에 통과하는 물의 부피의 중요성 외에도, 이 구조물에 대한 크라운의 통과-흐름의 속도는 이러한 벽 뒤에 있는 구조물에 추진력과 충동을 만드는 속도 요인의 중요성 때문에 매우 중요하다. 기본적으로 업스트림 경사면에서 장애물을 생성하고 업스트림 경사의 속도는 이 속도의 양을 줄이는 데 매우 효과적일 수 있습니다. 그러나 특정 경사면에서 최적의 장애물 레이아웃에 도달하기 위해 모델링하여 이 문제를 정확하게 조사해야 합니다. 본 연구에서는, FLOW-3D의 3차원 모델이 언급된 문제점을 조사하는 데 사용된다 [1].

2. 연구 역사

여러 연구는 파도가 해양 구조물을 덮어 넘나는 데 초점을 맞추고 있습니다. 이러한 방법은 지속적으로 바다 파도로부터 해안을 보호하기 위해 구조물의 오버 토핑을 정확하게 예측했다. 2002년까지 거의 6,500건의 시험이 실시되었습니다. 일반 파도의 물리적 모델도 미국에서 수행되었습니다 [2] . 무작위 파도의 가장 완벽한 세트는 오웬에 의해 완료되었다 (1980). 오웬은 오버 토핑과 바다 벽의 높이와 오버 토핑의 정도 사이의 관계를 연구하기 위해 물리적 모델 테스트의 번호를 수행 [3] . 그는 오버 토핑의 정도는 파도 높이 및 파도 기간과 같은 환경 조건뿐만 아니라 구조 재료의 기하학 및 유형에 따라 달라지며 있음을 보여주었습니다. 이러한 요인의 조합을 조사해야 합니다. 폰 마이어와 듀발 (1992) 연구의 또 다른 시리즈를 수행 [4] .

3. 재료 및 방법

이 연구에서는 68개의 다양한 형상이 모델링용 소프트웨어에 제공되며 다음 표 1에간단히 소개됩니다. 이 68 개의 다른 기하학에는 4 개의 다른 슬로프, 4 개의 다른 레이아웃 및 4 개의 다른 장애물 높이및 장애물이없는 4 개의 상태및 다른 경사에서만 포함 [5] . 그런 다음, 이러한 서로 다른 형상 및 상태는 FLOW-3D 3차원 모델을 사용하여 동일한 조건에서 평가 및 분석됩니다.

표 1. 변수지정.

4. 숫자 모델

FLOW-3D 소프트웨어는 3차원 유동 필드 분석을 통해 유체 역학 분야에서 강력한 유압 시뮬레이터 응용 프로그램입니다. 모델에서 지배하는 방정식은 다른 유사한 모델과 마찬가지로 Navier-Stokes 방정식과 질량 방정식의 보존[6]입니다.

이 응용 프로그램의 채널을 모델링하려면 일반 조건(모든 시스템의 시뮬레이션 포함), 물리적 조건, 형상 및 모델 해결 네트워크, 출력 및 관련 옵션을 조정해야 합니다. 온도도는 시스템 단위, SI 및 온도에 대해 선택되었습니다.

물리적 인 측면에서, 소프트웨어는 현상을 지배하는 물리학의 원칙에 따라 관련 조건을 선택할 수 있습니다. 이 연구를 지배하는 물리적 조건은 중력과 점도와 난기류입니다. 이 소프트웨어의 난기류는 5 가지 모델에 의해 자극되고이 연구에 사용되는 모델은 재정상화 그룹 (RNG)이었습니다. 난기류의 이 모델에서, K-모델에서 실험적으로 계산된 상수값은 암시적으로 파생된다[7].

그 후 유체를 정의해야 합니다. 이 연구의 선택된 유체는 섭씨 20도물[ 8]이다.

다음 단계는 형상을 정의하고 시뮬레이션에서 중요한 네트워크를 해결하는 것입니다 [9]. FLOW3D를 사용하면 소프트웨어에서 사용할 수 있는 도구로 많은 유체 현상을 묘사할 수 있습니다. 채널 형상을 정의하면 네트워크를 해결해야 합니다. 소프트웨어의 정의된 해결 네트워크는 네트워크 크기, 셀 수 및 X, Y 및 Z 및 경계 조건의 세 가지 좌표에서 해당 치수를 포함한 일반(입방) 해결 네트워크의 형태입니다. 네트워크 셀 치수의 크기가 작을수록 시뮬레이션을 위한 프로그램의 기능과 정밀도가 높을수록[10]이됩니다.

5. 결과

다른 그림에서 관찰할 수 있으므로 다이어그램은 두 가지 유형으로, 먼저 그림 1-4를 포함하는 소프트웨어의 직접 출력과 다른 숫자 5-7을 변경 프로세스의 다이어그램으로 포함합니다. 그러나 그림 1-4에서는 경사면 중 하나에서 출력이 소프트웨어 출력에서 직접 가져온다는 점을 언급해야 합니다.

언급된 수치와 관련하여, 이러한 속도는 장애물없이 상태의 상류 경사면에서 최대인 반면 방파제의 상류 경사면에서 가장 높은 속도 비율이 발생한다는 것을 이해할 수 있다. 흥미로운 점은 가장 낮은 속도는 일반적으로 방파제 크라운에 존재한다는 것입니다.

그림 5-8에서 볼 수 있듯이, 상류 방파제의 모든 다른 경사 상태에서, 가장 높은 유량 속도는 10cm 높이와 가장 낮은 속도의 장애물과 관련이 있으며 50cm 높이의 장애물과 관련이 있다. 그 이유는 장애물과의 충돌로 인해 잠재적 에너지로 변환되는 유동 운동 에너지의 가치가 장애물의 높이를 증가시켜 증가하기 때문입니다. 따라서, 높이가

그림 1. A1 모델링의 방파제의 흐름 속도.

그림 2. A2-1 모델링의 방파제의 흐름 속도.

Figure 3. Flow velocity on seawall in A2-3 modeling.

그림 4. A3-1 모델링의 방파제의 흐름 속도.

그림 5. 방파제 유형 A(61° 경사)의 흐름 속도 의 변화.

그림 6. 방파제 형 B (56 ° 경사)의 흐름 속도의 변화.

그림 7. 방파제 유형 C(51° 경사)의 흐름 속도 의 변화.

그림 8. 방파제 유형 D(45° 경사)의 흐름 속도 변경입니다.

해당 유동 운동 에너지는 각 장애물에 대한 흐름의 충돌에서 잠재적 에너지의 해당 높이로 변환되며, 흐름 속도가 잠시 0이 되고 장애물을 건너면 속도가 증가한다. 장애물의 높이가 낮은 것이든, 순간적인 제로 속도 상태가 줄어들고 흐름은 더 높은 속도와 함께 계속 움직입니다.

6. 결론

Also, as it can be observed, the highest difference of velocity in all the figures is between the obstacles with 10
cm height and the obstacles with 50 cm height. Also, this amount of difference in velocity for difference between the obstacles with 10 cm and 20 cm heights is higher than that of the differences in the obstacles with 20
cm and 30 cm heights which can be related to the special conditions in flow hydraulic in that range of height.

또한, 관찰할 수 있으므로 모든 수치에서 속도의 가장 높은 차이는 높이 가 10cm의 장애물과 높이가 50cm인 장애물 사이에 있습니다. 또한, 10cm와 20cm 높이의 장애물 사이의 차이에 대한 속도차이는 20cm 및 30cm 높이의 장애물의 차이보다 높으며, 이는 그 높이 범위에서 유압의 특별한 조건과 관련이 있을 수 있다.

이 논문 인용

메리암 데일라미-타리피, 메디 베다르반디-아스카르, 바히드 체기니, 사데 그 하그하이-부어(2016) FLOW-3D 소프트웨어를 사용하여 다양한 조건하에서 해벽에 흐르는 속도의 변화를 모델링한다. 해양 과학의 오픈 저널,06,317-322. doi: 10.4236/ojms.2016.62026

참조

  1. 1. Owen, M.W. (1980) Design of Seawalls Allowing for Wave Overtopping. Rep. EX924, Hydraulics Research Wallingford, England.   [Citation Time(s):1]
  2. 2. van der Meer, J.W. and Janssen, J.P.F.M. (1995) Wave Run-Up and Wave Overtopping at Dikes. In: Kobayashi, N. and Demirbilek, Z., Eds., Wave Forces on Inclined and Vertical Wall Structures, ASCE, New York.   [Citation Time(s):1]
  3. 3. CIRIA/CUR (1995) Manual on the Use of Rock in Hydraulic Engineering. CUR/RWS Report 169, A.A. Balkema, Rotterdam.   [Citation Time(s):1]
  4. 4. Pullen, T., Allsop, N.W.H., Bruce, T., Kortenhaus, A., Schuttrumpf, H. and van der Meer, J.W. (2007) EurOtop— Wave Overtopping of Seadefences and Related Structures Assessment Manual.
    http://www.overtopping-manual.com/manual.html?   [Citation Time(s):1]
  5. 5. De Wall, J.P. and Van der Meer, J.W. (1992) Wave Run-Up and Overtopping at Coastal Structures. ASCE, Proceeding of 23rd ICCE, Venice, 1758-1771.   [Citation Time(s):1]
  6. 6. De Gerloni, M., Franco, L. and Passoni, G. (1991) The Safety of Breakwaters against Wave Overtopping. Proceedings of ICE Conference on Breakwaters and Coastal Structures, Thomas Telford, London.   [Citation Time(s):1]
  7. 7. Fenton, J.D. (1988) The Numerical Solution of Steady Water Wave Problems. Computers & Geosciences, 14.
    http://dx.doi.org/10.1016/0098-3004(88)90066-0   [Citation Time(s):1]
  8. 8. Owen, M.W. (1982) Design of Seawalls Allowing for Wave Overtopping. Rep. EX924, Hydraulics Research Wallingford, England.   [Citation Time(s):1]
  9. 9. Allsop, W., Bruce, T., Pearson, J. and Besley, P. (2006) Wave Overtopping at Vertical and Steep Seawall.   [Citation Time(s):1]
  10. 10. TAW (1974) Technical Advisory Committee on Protection against Inundation, Wave Run-Up and Overtopping. Government Publishing Office, The Hague.   [Citation Time(s):1]
Abf - Three-dimensional view of the abbot from short to long to short

Flow-3D 수치 모형을 이용한 파동 감소에 대한 규칙적인 레이아웃으로 식생 고도 변화 효과 연구

세예드 아마드가 헤리 네 자드 1 , Mehdi Behdarvandi Askar  2 , 모하마드 안사리 고이 가르 3, 에산 파르시 4
1 공학, 해안, 항만 및 & amp; 해양 구조물 _ 코람 샤르 해양 과학 기술 대학교
2 코람 샤르 해양 과학 기술 대학교 해양 공학부 해양 구조학과
3 이란 카라 지 테헤란 대학교 농업 및 천연 자원 대학 관개 및 매립 공학과.
4 연구 전문가, Arvand Water and Energy Consulting Engineers Company, Ahvaz, Iran.

Abstract

The development of water waves through submerged and non-submerged vegetation is accompanied by a loss of energy through the resistive force of the vegetation, resulting in a decrease in wave height. Wave damping by vegetation is a function of cover characteristics such as geometry and structure, immersion ratio, density, hardness, and spatial arrangement, as well as wave conditions such as input wave height, duration, and wave direction. In the present study, the effect of geometric arrangement of vegetation with variable height on wave damping has been investigated using the Flow 3D numerical model. For this purpose, a channel with a length of 480 cm and a width of 10.8 cm, which has been previously used by Cox and Wu (2015) to study the effect of plant density with variable height on wave damping, is modeled. The operation of the three arrangements, including long to short arrangement, short to long arrangement, and zigzag arrangement, is examined under four different waves, all of which are linear waves. It should be noted that in this study, wave height is considered as an damping index. The results obtained by measuring the height of the waves at four different points along the channel show that the behavior of the waves in dealing with different arrangements follows a fixed pattern and also changes in the geometry of the vegetation can greatly lead to Increase the damping of the waves. The results show that a change in height arrangement can cause a change in damping of up to 7.1%.

Keywords : Green belt , wave , geometric layout , vegetation

물에 잠긴 초목과 물에 잠기지 않은 초목을 통한 물결의 발달은 초목의 저항력을 통한 에너지 손실을 동반하여 파고가 감소합니다. 식생에 의한 파동감쇠는 기하와 구조, 몰입도, 밀도, 경도, 공간배열 등 커버 특성과 입력파동 높이, 지속시간, 파동방향 등의 파동조건의 함수입니다.

본 연구에서는 Flow 3D 수치 모델을 사용하여 가변 높이 식물이 파동 댐핑에 미치는 기하학적 배치가 조사되었습니다. 이를 위해 Cox와 Wu (2015)가 이전에 파동 댐핑에 대한 가변 높이의 발전소 밀도가 미치는 영향을 연구하기 위해 사용한 길이 480cm, 폭 10.8cm의 채널을 모델링합니다.

장파에서 단파, 단파에서 장파까지, 지그재그 배열을 포함한 세 가지 배열의 작동은 4개의 다른 파장에서 조사됩니다. 모두 선형파입니다.

본 연구에서는 파고가 감쇠 지수로 간주된다는 점에 유의해야 합니다.

채널을 따라 네 곳의 서로 다른 지점에서 파도의 높이를 측정하여 얻은 결과는 다른 배열을 다루는 파도의 동작이 고정된 패턴을 따르며 또한 초목의 기하학적인 변화가 파도의 감쇠를 증가 시키는 것으로 크게 이어질 수 있다는 것을 보여줍니다.

결과는 높이 배열의 변화가 최대 7.1%의 댐핑 변화를 일으킬 수 있음을 보여줍니다.

Figure 1 - Geometry used by Cox and Wu (2015) to study the effect of plant density on wave damping
Figure 1 – Geometry used by Cox and Wu (2015) to study the effect of plant density on wave damping
Figure 2 - Schematic of Erie wave
Figure 2 – Schematic of Erie wave
Abf - Three-dimensional view of the abbot from short to long to short
Abf – Three-dimensional view of the abbot from short to long to short

References

خلیلی نفت­چالی، آ. خزیمه­نژاد، ح. اکبرپور، ا. ورجاوند، پ. 1394. بررسی آزمایشگاهی تأثیر تراکم پوشش گیاهی بر مشخصه‌های جریان غلیظ. نشریه آبیاری و زهکشی ایران. 9 (1): 95-83.
زارعی، م. فتحی­مقدم، م. داوودی، ل. 1395. بررسی اثر پوشش گیاهی ساحلی بر میرایی نیروی مخرب امواج منفرد ناشکنا در سواحل شیبدار. نشریه مهندسی آبیاری و آب ایران. 7 (26): 75-62.
گرمئی، ا. امامی، ح. خراسانی، ر. 1396. اثر تراکم سه نوع پوشش گیاهی بر میزان رواناب و رسوب در حاشیه شهر مشهد. نشریه آبیاری و زهکشی ایران. 11 (1): 20-11.
فضلی، س. نور، ح. 1396. شبیه‌سازی و ارزیابی اثر سناریوهای مختلف درصد پوشش گیاهی بر فرسایش خاک. نشریه آبیاری و زهکشی ایران. 11 (4): 571-562.
قنبری عدیوی، ا. فتحی مقدم، م. 1393. مروری بر تحقیقات استهلاک و میرایی امواج دریا از طریق پوشش گیاهی ساحلی. فصلنامه علوم و فناوری دریا. 18 (70): 62-54.
معتمدی­نژاد، ع. فتحی­مقدم، م. زارعی، م. 1394. بررسی آزمایشگاهی اثر پوشش گیاهی ساحلی بر کاهش نیروی امواج شکنا. دهمین سمینار بین المللی مهندسی رودخانه. دانشگاه شهید چمران اهواز، اهواز، ایران.
میرزاخانی، گ. قنبری عدیوی، ا. فتاحی­نافچی، ر. 1398. میرایی موج توسط پوشش گیاهی صلب در سواحل. دومین همایش ملی مدیریت منابع طبیعی با محوریت آب، سیل و محیط زیست. دانشگاه گنبد کاووس، گنبد کاووس، ایران.
Asano, T. S. Sutsui, T. and Sakai.T. 1988. Wave damping characteristics due to seaweed. Proceedings of the 35th Coastal Engineering Conference in Japan. JSCE. 138-142 (in Japanese).
Asano, T., Deguchi, H. and N. Kobayashi. 1992. Interactions between water waves and vegetation. Proceedings of the 23rd International Conference on Coastal Engineering. ASCE. 2710-2723.
Augustin, L.N., Irish, J.L. and Lynett, P. 2009. Laboratory and numerical studies of wave damping by emergent and nearemergent wetland vegetation. Coastal Engineering. 56(3): 332-340.
Cavallaro L., Re, C.L., Paratore, G., Viviano, A. and Foti, E. 2010. Response of Posidonia oceanic to wave motion in shallowwaters: Preliminary experimental results. Proceedings of the 32nd International Conference on Coastal Engineering. Coastal Engineering Research Council. 1-10.
Cook, H.L. and Campbell, F.B. 1939. Characteristics of some meadow strip vegetation. Agricultural Engineering. 20:345-348.
Cooper, N.J. 2005. Wave dissipation across intertidal surfaces in the Wash Tidal inlet, Eastern England. Journal of Coastal Research. 21(1): 28-40.
Dean, R.G. 1979. Effects of vegetation on shoreline erosional processes. Wetland Function and Values: The State of Our Understanding. 1: 415-426.
Dean, R.G., and Dalrymple, R.A. 1991. Water Wave Mechanics for Engineers and Scientist. World Scientific Publishing.Singapore.
Dubi, A. 1995. Damping of water waves by submerged vegetation: A case study on Laminaria hyperborea. PhD thesis. University of Trondheim, the Norwegian Institute of Technology, Trondheim, Norway.
Fathi Moghadam, M., Drikundi, K.h., Masjidi, A. and M. 2012. Investigation of the Effect of Vegetation Density and Flexibility on Roughness Coefficients in Riverside and Flood Plains, Iranian Water Resources Research Quarterly, Year 8, Issue 2, Fall 91.
Fathi Moghadam, M. and Zaraei, M. 2016. Investigation of the Effect of Coastal Vegetation on the Damping of Destructive Force of Unbreakable Individual Waves on Shabidar Coasts, Journal of Irrigation and Water Engineering, Year 7, No. 26.
Furukawa, K., Wolanski, E. and Mueller, H. 1997. Currents and sediment transport in mangrove forests. Estuar Coast Shelf Sci 44:301–310.
Harada, K. and Imamura, F. 2006. Experimental study on the resistance by mangrove under unsteady flow, Proc. Congress. Asian and Pacific Coastal Engineering Dalia, 984-975.
Jellilund, R., M. Zeid Ali, L. Nouri Hindi and M. 2012. Investigating the advantages and disadvantages of protection and organization of beaches with vegetation against morphological changes, Fifth National Conference and Specialized Environment Exhibition, 90.
Journal 629, Guide to the Design and Implementation of a Coastal Protection Structure.
Kongko, W. 2004. Study on tsunami energy dissipation in mangrove forest, Master Thesis Report, wate University, Japan, 43 pages.
Kutija, V. and Erduran, K. S. 2003. Quasi-three-dimensional numerical model for flow through flexible, rigid, submerged and non-sub merged vegetation. Journal of Hydro informatics. 35(3): 189-202.
Li, R.M. and Shen, H.W. 1973. Effect of tall vegetations on flow and sediment. Journal of the Hydraulics Division, ASCE. 99(5):739-814.
Wu, W.C. and Cox, D, T. 2015. Effects of Vertical Variation in Vegetation Density on Wave Attenuation. Journal of Waterway, Port, Coastal and Ocean Engineering. Volume 142 Issue 2.

Interaction between oblique waves and arc-shaped breakwater

Interaction between oblique waves and arc-shaped breakwater: Wave action on the breakwater and wave transformation behind it

XinyuHanaShengDongaYizhiWangb
aCollege of Engineering, Ocean University of China, Qingdao, 266100, China
bShandong Harbour Engineering Group Co., Ltd., Rizhao, 276826, China

Highlights

Interaction of oblique waves and the arc-shaped breakwater was simulated.

Wave force and pressure distribution along central axis were analysed.

Arc curvature has little effect on the maximum wave force of different sections.

Overtopping-induced Hmax behind breakwater up to 0.7 times of incident wave height.

Abstract

The hydrodynamic interaction between oblique waves and an arc-shaped breakwater and the wave field behind it. A three-dimensional computational fluid dynamic model was used to simulate the interaction between the oblique waves and arc-shaped breakwater. The pressure distribution and wave force in the different sections under different wave directions were measured by experiments to validate the numerical results. The pressure distribution and wave force in the arc-shaped vertical part of the breakwater along the central axis were further analysed using numerical model. The maximum positive and negative forces in each section along the central axis were compared. The results indicated that the arc curvature exerted little effect on the maximum wave force in the different sections. The wave height behind the breakwater was obviously smaller than that at the front. With the decrease in the incident angle, the influence of diffraction on the wave field gradually decreased. Under east–southeast waves, the maximum wave height behind the breakwater caused by overtopping was approximately 0.7 times the incident-wave height. In the spatial distribution of the wave period behind the breakwater, some areas with smaller periods existed, which may be caused by the overtopping flow that broke behind the breakwater.

경사파와 호 모양의 방파제와 그 뒤에 있는 파동 장 사이의 유체 역학적 상호 작용. 3 차원 전산 유체 역학 모델을 사용하여 사선 파와 호 모양의 방파제 사이의 상호 작용을 시뮬레이션했습니다.

서로 다른 파동 방향에서 서로 다른 섹션의 압력 분포와 파력은 수치 결과를 검증하기 위해 실험을 통해 측정 되었습니다. 방파제 중심 축을 따라 호 모양의 수직 부분의 압력 분포와 파력은 수치 모델을 사용하여 추가로 분석되었습니다.

중심 축을 따라 각 섹션에서 최대 양의 힘과 음의 힘을 비교했습니다. 결과는 아크 곡률이 다른 섹션에서 최대 파력에 거의 영향을 미치지 않음을 나타냅니다. 방파제 뒤의 파도 높이는 정면보다 분명히 작았습니다. 입사각이 감소함에 따라 파동 장에 대한 회절의 영향이 점차 감소했습니다.

동-남동 파 하에서 오버 탑으로 인한 방파제 뒤의 최대 파고는 입사 파고의 약 0.7 배였다. 방파제 뒤의 파동주기의 공간적 분포에는 방파제 뒤에서 파열 된 과잉 흐름에 의해 발생할 수 있는 더 작은주기를 가진 일부 지역이 존재했습니다.

Keywords

Arc-shaped breakwater3D numerical modelWave forcePressure distributionWave height and period behind breakwater

Figures -Interaction between oblique waves and arc-shaped breakwater
Figures -Interaction between oblique waves and arc-shaped breakwater
Figures-Interaction between oblique waves and arc-shaped breakwater2
Figures-Interaction between oblique waves and arc-shaped breakwater2

Recommended articles

Figure 1. Sketch map of the port Laozi on Lake Hongze

FLOW-3D software for substantiation the layout of the port water area

항구 수역의 레이아웃을 입증하기 위한 FLOW-3D 소프트웨어

B Pan1

  • and N Belyaev2
    1 Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia
    pan3.b@edu.spbstu.ru

Abstract

방파제 설계시 항만 내 수역 및 접안 시설에서 중대한 파도 발생을 배제 할 필요가 있다. 이 기사는 항구 지역의 물 이동 계산 결과를 제시합니다. 항구 입구의 위치와 주어진 물체의 크기에 대한 가장 안전한 옵션이 확인되었습니다. FLOW – 3D 프로그램을 사용하여 항구 수역의 예비 배치 단계에서 계획되고 건설적인 솔루션을 선택할 가능성이 입증되었습니다.

Introduction

항구를 설계 할 때 계류선의 가장 합리적인 구성과 항구 수역의 레이아웃을 선택하여 항구의 영토를 배치하는 것이 필요합니다. 러시아 연방의 항구 수역 배치는 항구 수역 요소에 대한 사양을 포함하는 해로, 페어웨이 및 기동 구역에 대한 설계 표준의 요구 사항에 따라 수행됩니다 [1, 2].

항구 물은 파도, 퇴적물 축적 및 얼음으로부터 보호되어야 합니다. 항구 계획을 작성할 때, 선박의 기동 및 연안 계류 중뿐만 아니라 선적 및 하역 및 기타 항구 운영 중에 선박 계류에 대한 정상적인 조건이 생성되도록 파도로부터 수역을 보호하는 정도를 제공해야 합니다.

설계 결정은 새로운 포트를 설계하거나 기존 포트를 개발할 때 물리적 또는 수학적 모델링을 기반으로 합니다 [2]. 항구 수역에서 계산 된 물 흐름의 매개 변수는 수문 기상 조사, 장기 현장 관찰 및 실험실 연구의 결과를 기반으로 하도록 권장됩니다.

공학 수문 기상 측량 데이터가 불충분하면 계산 방법을 기반으로 설계 폭풍의 파도 매개 변수를 결정할 수 있습니다. 사용된 계산 방법이 국제 실무에서 동일한 목적으로 채택된 방법 (모델)에 부합하는지 표시하는 것이 좋습니다 [3].

Figure 1. Sketch map of the port Laozi on Lake Hongze
Figure 1. Sketch map of the port Laozi on Lake Hongze
Figure 2. The location of the port entrance on Lake Hongze: a – variant 1; b – variant 2; c – variants 3-5
Figure 2. The location of the port entrance on Lake Hongze: a – variant 1; b – variant 2; c – variants 3-5
Figure 3. Port water area plan
Figure 3. Port water area plan
Figure 4. Modeling of variant 1 with the movement of waves in the port water area
Figure 4. Modeling of variant 1 with the movement of waves in the port water area
Figure 5. Modeling of variant 2: a is prevailing movement of water towards the enclosed water area; b is prevailing reverse movement of water
Figure 5. Modeling of variant 2: a is prevailing movement of water towards the enclosed water area; b is prevailing reverse movement of water
Figure 6. Modeling of variant 3
Figure 6. Modeling of variant 3
Figure 7. Modeling of variant 4
Figure 7. Modeling of variant 4
Figure 8. Modeling of variant 5
Figure 8. Modeling of variant 5
Figure 9. Plan of the port water area with design points
Figure 9. Plan of the port water area with design points
Figure 10. Change in water depth at point A: a – variant 1; b – variant 2
Figure 10. Change in water depth at point A: a – variant 1; b – variant 2
Figure 11. Change in water depth at point A: a – variant 3; b – variant 4; c – variant 5
Figure 11. Change in water depth at point A: a – variant 3; b – variant 4; c – variant 5
Figure 12. Change in water depth at points A (a) and C (b) for variant 3
Figure 12. Change in water depth at points A (a) and C (b) for variant 3
Figure 13. Change in water depth at points A (a) and B (b) for variant 3
Figure 13. Change in water depth at points A (a) and B (b) for variant 3
Figure 14. Scheme of vessel traffic: a – variant 3; b – variant 4
Figure 14. Scheme of vessel traffic: a – variant 3; b – variant 4

References
[1] SP 350.1326000.2018. 2018 Norms for technological design of sea ports (Moscow:
Standartinform) p 226
[2] SP 444.1326000.2019. 2019 Standards for the design of sea channels, fairways and
maneuvering areas (Moscow: Standartinform) p 62
[3] SP 38.13330.2012. 2014 Loads and impacts on Hydraulic structures (from wave, ice and ships)
(Moscow: Ministry of Regional Development of the Russian Federation) p 112
[4] Rijnsdorp D P Smit P B and Zijlema M 2012 Non-hydrostatic modelling of infragravity waves
using SWASH. Proceedings of 33rd Conference on Coastal Engineering. pp 1287–1299
[5] Kantardgi I G Zheleznyak M J 2016 Laboratory and numerical study of waves in the port area.
Magazine of Civil Engineering No 6 pp 49–59 DOI: 10.5862/MCE.66.5
[6] Zheleznyak M J Kantardgi I G Sorokin M S and Polyakov A I 2015 Resonance properties of
seaport water areas Magazine of Civil Engineering № 5(57) pp 3-19 DOI:10.5862/MCE.57.1
[7] Kantarzhi I Zuev N Shunko N 2014 Numerical and physical modelling of the waves inside the
new marina in Gelendjik (Black Sea) Application of physical modelling to port and coastal
protection. Proceedings of 5th international conference Coastlab (Varna) Vol 2 pp 253–262
[8] Makarov K N and Chebotarev A G 2015 Breakwater placement at the root of a seawall
Magazine of Civil Engineering № 3(55) pp 67-78 DOI: 10.5862/MCE.55.8
[9] Belyaev N D Lebedev V V and Alexeeva A V 2017 Investigation of the soil structure changes
under the tsunami waves impact on the marine hydrotechnical structures V 10 № 4 pp 44-52
DOI: 10.7868/S2073667317040049
[10] Lebedev V V Nudner I S and Belyaev N D 2018 The formation of the seabed surface relief near
the gravitational object Magazine of Civil Engineering No 79(3) pp 120–131 DOI:
10.18720/MCE.79.13
[11] Kofoed-Hansen H Sloth P Sørensen O R Fuchs J 2000 Combined numerical and physical
modelling of seiching in exposed new marina Proceedings of 27th international conference of
coastal engineering pp 3600–3614
[12] Smit P Stelling G and Zijlema M 2011 Assessment of nonhydrostatic wave-flow model
SWASH for directionally spread waves propagating through a barred basin Proceedings of
ACOMEN 2011 pp 1–10
[13] Zijlema M Stelling G Smit P 2011 SWASH: An operational public domain code for simulating
wave fields and rapidly varied flows in coastal waters. Coastal Engineering. № 10(58). pp 992–
1012
[14] FLOW-3D® 2008 User’s Manual Version 9.3 Flow Science Inc p 821
[15] Pan Bayan and Belyaev N D 2019 Week of Science SPbPU: Proceedings of an international
scientific conference The best reports. pp 3-7
[16] Girgidov A A 2011 Hybrid simulation in hydrotechnical facilities design and FLOW-3D as a
tool its realization Magazine of Civil Engineering №3 pp 21-27
[17] Girgidov A A 2010 Proceeding of the VNIIG vol 260. pp 12-19

[18] Vasquez J A Walsh B W 2009 CFD simulation of local scour in complex piers under tidal flow,
33rd IAHR Congress: Water Engineering for a Sustainable Environment, © 2009 by
International Association of Hydraulic Engineering & Research (IAHR) ISBN: 978-94-90365-
01-1.
[19] Shan-Hwei Ou Tai-Wen Hsu and Jian-Feng Lin 2010 Experimental and Numerical Studies on
Wave Transformation over Artificial Reefs Proceedings of the International Conference on
Coastal Engineering (Shanghai, China) No 32
[20] Hirt C and Nichols B 1980 Volume of Fluid Method for the Dynamics of Free Boundaries
Journal Comp. Phys 39 p 201.

Figure 7. Formation of incident and reflected waves.

Investigate Impact Force of Dam-Break Flow against Structures by Both 2D and 3D Numerical Simulations

2D 및 3D 수치 시뮬레이션에 의한 댐 붕괴유동의 구조물 충격력 조사

1 Faculty of Water Resources Engineering, Thuyloi University, 175 Tay Son, Dong Da, Ha Noi 116705, Vietnam
2 Hydraulic Construction Institute, 3/95 Chua Boc, Dong Da, Ha Noi 116705, Vietnam
* Author to whom correspondence should be addressed.
Academic Editor: Costanza Aricò
Water 2021, 13(3), 344;

Abstract

본 논문의 목적은 일부 2D 및 3D 수치 모델이 침수 지역에 고립된 건물 또는 건물 배열이 있는 곳에서 홍수 파동을 시뮬레이션하는 능력을 조사하는 것이었습니다.

먼저, 제안된 2D 수치 모델은 구조화된 메시에서 2D 얕은 물 방정식(2D-SWEs)을 해결하기 위한 유한 볼륨 방법(FVM)을 기반으로 했습니다.

FDS (flux-difference splitting)은 정확한 질량 균형을 얻기 위해 사용되었고 Roe 체계는 Riemann 문제를 근사하기 위해 호출되었습니다.

둘째, 상업적으로 이용 가능한 3D CFD 소프트웨어 패키지가 선택되었으며, 여기에는 두 가지 난류 모델이 포함된 Flow 3D 모델이 포함되어 있습니다.

RNG(Renormalized Group) 및 LES(Large-eddy Simulation)를 사용하는 레이놀즈 평균 Navier-Stokes(RAN)입니다. 댐 붕괴 흐름으로 인한 장애물에 대한 충격력의 수치 결과는 3D 솔루션이 2D 솔루션보다 훨씬 낫다는 것을 보여주었습니다.

건물 배열에 작용하는 충격력의 3D 수치 힘 결과를 보유하고 있는 실험 데이터와 비교함으로써, 속도 유도력이 동적 힘에 미치는 영향은 Froude 숫자의 함수와 사고 파동의 수심 함수에 의해 정량화 되었습니다. 또한, 우리는 힘의 강도의 피크 값의 3D 계산 결과에 대한 초기 물 단계와 댐 붕괴 폭의 영향을 조사했습니다.

The aim of this paper was to investigate the ability of some 2D and 3D numerical models to simulate flood waves in the presence of an isolated building or building array in an inundated area. Firstly, the proposed 2D numerical model was based on the finite-volume method (FVM) to solve 2D shallow-water equations (2D-SWEs) on structured mesh. The flux-difference splitting method (FDS) was utilized to obtain an exact mass balance while the Roe scheme was invoked to approximate Riemann problems. Secondly, the 3D commercially available CFD software package was selected, which contained a Flow 3D model with two turbulent models: Reynolds-averaged Navier-Stokes (RANs) with a renormalized group (RNG) and a large-eddy simulation (LES). The numerical results of an impact force on an obstruction due to a dam-break flow showed that a 3D solution was much better than a 2D one. By comparing the 3D numerical force results of an impact force acting on building arrays with the existence experimental data, the influence of velocity-induced force on a dynamic force was quantified by a function of the Froude number and the water depth of the incident wave. Furthermore, we investigated the effect of the initial water stage and dam-break width on the 3D-computed results of the peak value of force intensity.

Keywords: dam-break wave2D numerical modelFlow 3D modelstructuresimpact force

Introduction

홍수 위험 분석에 따른 도시 계획은 최근에 큰 연구 과제였습니다.

건물 또는 건물 그룹에 대한 홍수 파동의 영향에 대한 연구는 하류 지역에 대한 조기 경고 또는 안전 의식 향상에 중요한 역할을 했습니다. 기본적으로 댐 파괴 흐름에 대한 연구는 실험 측정이나 수치 시뮬레이션을 통해 추정 할 수 있습니다 [1,2,3,4,5,6].

컴퓨터 처리 능력의 증가로 인해 불연속 흐름에 대한 수치 연구가 비용 효율적이되었습니다. 지난 10 년 동안 얕은 물 솔버는 정확성과 계산 능력면에서 크게 향상되었습니다. 침수 가능 지역의 수심 및 속도 프로파일과 같은 유체 역학적 매개 변수에 많은주의를 기울였습니다 [1,2,3,4,5,6,7,8].

Migot et al. [9]는 도시 홍수의 실험적 모델링에 관한 많은 기사를 검토했습니다. 그 논문에 언급 된 45 개의 작품 중 단 4 개의 프로젝트 만이 장애물에 가해지는 일정한 또는 비정상적인 흐름의 힘 또는 압력을 측정했습니다.

또한 물리적 및 2D 수치 모델에서 건물 또는 건물 그룹에 돌발 홍수가 미치는 영향에 대한 연구는 거의 없었습니다. 얕은 물 모델은 [10,11]에서 고립된 장애물에 대한 충격의 힘을 예측하는데 사용되었습니다.

한편 Shige-eda [12]는 액체와 건물 배열 간의 상호 작용을 결정하기 위해 물리적 모델과 2D 수치 체계를 선택했습니다. Aureli와 Shige-eda는 수직 속도와 가속도를 무시하기 때문에 댐 파괴 흐름의 힘을 추정하기 위한 2D 얕은 물 방정식 (SWE)의 단점을 보여주었습니다 [10,12].

Migot [9]은 또한 장애물 주변의 시뮬레이션된 홍수 흐름에 대한 2D SWE에 대한 여러 출판물이 있었지만 이 주제에 대한 3D 수치 모델에 대한 연구는 거의 없다고 지적했습니다. 최근 전산 유체 역학 (CFD) 3D 시뮬레이션은 유체 흐름과 관련된 문제를 해결하기위한 광범위한 도구가되었습니다.

댐 파괴 파의 특성은 [13,14,15,16]에 의해 주목되었고 Issakhov [17]는 다양한 종류의 장애물이 압력 분포에 미치는 영향을 조사하기 위해 CFD 방법을 사용했습니다. 그들은 분포가 댐 표면에서 3 배 더 낮다는 것을 밝혔다.

Aureli [10]는 댐 파괴 파가 구조물에 미치는 영향의 정적 힘을 평가하기 위해 실험 테스트와 2D 및 3D 수치 모델을 사용했습니다. Mokarani [18]는 댐 브레이크 흐름 영향의 VOF 시뮬레이션에서 피크 압력 안정성 조건을 연구했습니다.

앞서 언급한 작품에서 구조물이나 구조물 군에 작용하는 힘은 압력에 의한 정 수력 또는 정력이었다. 한편, 급류에서 속도로 인한 힘은 압력 력보다 크거나 같았습니다 [19]. Armanini [20]는 정상 흐름에 대해이 항을 추정하기 위한 분석적 표현 만을 제시했습니다. 우리가 아는 한, 건물 그룹에 작용하는 비정상 흐름의 동적 힘을 생성하기 위해 2D 및 3D 수학적 모델을 모두 사용하는 작업은 없습니다.

따라서 본 연구에서는 제안된 2D 수치 모델과 3D 수학적 모델 모두에 의해 고립 된 장애물 또는 장애물 그룹에 대한 급격한 비정상 흐름의 테스트 사례를 재현했습니다. 수심 및 유속 수문 그래프와 같은 몇 가지 수력 학적 특성이 추정되었으며 측정 된 데이터와 매우 잘 일치했습니다.

특히 댐 브레이크 흐름이 서로 다른 건물에 가하는 동적인 힘도 시뮬레이션했습니다. 속도 유도 힘이 동적 힘에 미치는 영향 수준을 나타내는 매개 변수는 Froude 수와 입사 파동의 수심의 함수인 것으로 밝혀졌습니다. 또한 붕괴된 댐 사이트 폭 (b)과 초기 수위 (h0)는 충격력의 최대 값에 영향을 미치는 변수로 고려되었습니다.

Figure 1. (a) Configuration of experiment test (dimension in meters); (b) Gauges on the vertical front face of building.
Figure 1. (a) Configuration of experiment test (dimension in meters); (b) Gauges on the vertical front face of building.
Figure 2. (a) Distributed pressure profiles at centerline of front face of column; (b) Comparison of load-time histories simulated by different numerical models
Figure 2. (a) Distributed pressure profiles at centerline of front face of column; (b) Comparison of load-time histories simulated by different numerical models
Figure 3. Group of buildings in flooded area.
Figure 3. Group of buildings in flooded area.
Figure 4. Water depth and u-velocity profiles at gauge b.
Figure 4. Water depth and u-velocity profiles at gauge b.
Figure 5. Water hydrographs at gauges a and c.
Figure 5. Water hydrographs at gauges a and c.
Figure 6. Velocity component profiles at gauges a and c.
Figure 6. Velocity component profiles at gauges a and c.
Figure 7. Formation of incident and reflected waves.
Figure 7. Formation of incident and reflected waves.
Figure 8. Snapshots of streamlines of Froude number at different times: 1.0 s, 2.0 s, 5.0 s and 10 s.
Figure 8. Snapshots of streamlines of Froude number at different times: 1.0 s, 2.0 s, 5.0 s and 10 s.
Figure 9. Force in the flow direction exerted on 6 buildings.
Figure 9. Force in the flow direction exerted on 6 buildings.
Figure 10. The linear regression between forces per unit width (F) and q2b/h0.
Figure 10. The linear regression between forces per unit width (F) and q2b/h0.

Conclusions

댐 붕괴 흐름으로 인한 홍수 파도는 높은 속도 또는 큰 깊이가 관련되었을 때 건물에 큰 영향을 미칩니다. 본 논문에서는 2D 및 3D 수치 모델의 건물 및 건물 그룹에 대한 빠른 흐름에 의해 발생하는 유압 특성과 충격 부하를 추정할 수 있는 능력을 조사했습니다. 얕은 물 방정식에 기초한 2D 수학 모델은 FDS 방법으로 해결되었으며, FDS 방법은 최신 버전의 Flow 3D 유체 역학 모델과 함께 사용되었습니다. 연구의 주요 발견은 다음과 같습니다.
(1) 수심 또는 속도 프로파일을 공식화하기 위해 2D 및 3D 수치 솔루션은 모두 매우 유사합니다. 제안된 2D 수치 모델은 정적 힘의 최대 값 뿐만 아니라 수심 및 속도 구성 요소를 포함하는 유압 특성을 예측하는 데 적합합니다. 그러나 LES 및 RAN 난류 모듈이 포함된 3D 유체역학 모델은 2D 얕은 흐름 모델이 1개만 제공하는 동안 두 개의 최고 충격 부하를 잘 포착할 수 있습니다. 일반적으로 3D 결과는 실험 결과와 더 가깝습니다.
(2) 여러 건물에 대한 정적 및 동적 힘은 모두 LES 모듈을 사용하여 Flow 3D에 의해 계산되었습니다. 건물에서 속도에 의한 힘과 압력의 역할은 위치에 따라 다릅니다. 댐 현장 근처에서, 속도 유도 힘은 댐 파괴 파동의 주 방향에서 멀리 떨어져 있거나 두 번째 배열에서 압력 힘이 더 중요합니다. 속도 유도 힘의 영향은 매개 변수 α에 의해 정량화되며, 이는 사고파의 Froude 숫자와 수심 함수로 수행됩니다. q2b/h0과 정적 힘과 동적 힘의 피크 강도 사이의 선형 회귀 관계는 합리적인 R-제곱 양으로 해결됩니다.

추가 연구에서, 제시된 2D 수치 모델의 견고성과 효과는 더 명확하게 드러날 것입니다. 대규모 도메인에 대한 홍수 흐름을 시뮬레이션하는 데 쉽게 적용할 수 있습니다. 또한, α 매개변수의 제안된 방정식(21)은 실제 사례 연구에서 다운스트림 영역의 건물에 대한 속도 유도 힘의 영향을 정확하게 평가하기 위한 매우 의미가 있습니다. 이 매개 변수의 정확도 수준을 높이려면 서로 다른 조건에서 장애물에 작용하는 여러 가지 힘 실험이 구현되어야 합니다.

References

  1. Testa, G.; Zuccala, D.; Alcrudo, F.; Mulet, J.; Frazao, S.S. Flash flood flow experiment in a simplifed urban district. J. Hydraul. Res. 200745, 37–44. [Google Scholar] [CrossRef]
  2. Soares-Frazao, S.; Zech, Y. Dam-break flow through an idealized city. J. Hydraul. Res. 200846, 648–665. [Google Scholar] [CrossRef]
  3. Soares-Frazão, S.; Zech, Y. Experimental study of dam-break flow against an isolated obstacle. J. Hydraul. Res. 200745, 27–36. [Google Scholar] [CrossRef]
  4. Soares-Frazão, S. Experiments of dam-break wave over a triangular bottom sill. J. Hydraul. Res. 200745, 19–26. [Google Scholar] [CrossRef]
  5. di Cristo, C.; Evangelista, S.; Greco, M.; Iervolino, M.; Leopardi, A.; Vacca, A. Dam-break waves over an erodible embankment: Experiments and simulations. J. Hydraul. Res. 201856, 196–210. [Google Scholar] [CrossRef]
  6. Evangelista, S. Experiments and numerical simulations of dike erosion due to a wave impact. Water 20157, 5831–5848. [Google Scholar] [CrossRef]
  7. Li, Y.L.; Yu, C.H. Research on dam break flow induced front wave impacting a vertical wall based on the CLSVOF and level set methods. Ocean Eng. 2019178, 442–462. [Google Scholar] [CrossRef]
  8. Özgen, I.; Zhao, J.; Liang, D.; Hinkelmann, R. Urban flood modeling using shallow water equations with depth-dependent anisotropic porosity. J. Hydrol. 2016541, 1165–1184. [Google Scholar] [CrossRef]
  9. Mignot, E.; Li, X.; Dewals, B. Experimental modelling of urban flooding: A review. J. Hydrol. 2019568, 334–342. [Google Scholar] [CrossRef]
  10. Aureli, F.; Dazzi, A.; Maranzoni, A.; Mignosa, P.; Vacondio, R. Experimental and numerical evaluation of the force due to the impact of a dam break wave on a structure. Adv. Water Resour. 201576, 29–42. [Google Scholar] [CrossRef]
  11. Milanesi, L.; Pilotti, M.; Belleri, A.; Marini, A.; Fuchs, S. Vulnerability to flash floods: A simplified structural model for masonry buldings. Water Resour. Res. 201854, 7177–7197. [Google Scholar] [CrossRef]
  12. Shige-eda, M.; Akiyama, J. Numerical and experimental study on two dimensional flood flows with and without structures. J. Hydraul. Eng. 2003129, 817–821. [Google Scholar] [CrossRef]
  13. Cagatay, H.O.; Kocaman, S. Dam break flows during initial stage using SWE and RANs approaches. J. Hydraul. Res. 201048, 603–611. [Google Scholar] [CrossRef]
  14. Yang, S.; Yang, W.; Qin, S.; Li, Q.; Yang, B. Numerical study on characteristics of dam break wave. Ocean Eng. 2018159, 358–371. [Google Scholar] [CrossRef]
  15. Robb., D.M.; Vasquez., J.A. Numerical simulation of dam break flows using depth averaged hydrodynamic and three dimensional CFD models. In Proceedings of the 22nd Canadian Hydrotechnical Conference, Ottawa, ON, Canada, 28–30 April 2015. [Google Scholar]
  16. Kocaman, S.; Evangelista, S.; Viccione, G.; Guzel, H. Experimental and Numerical analysis of 3D dam break waves in an enclosed domain with a single oriented obstacles. Environ. Sci. Proc. 20202, 35. [Google Scholar] [CrossRef]
  17. Issakhov, A.; Zhandaulet, Y.; Nogaeva, A. Numerical simulation of dambreak flow for various forms of the obstacle by VOF method. Int. J. Multiph. Flow 2018109, 191–206. [Google Scholar] [CrossRef]
  18. Mokarani, C.; Abadie, S. Conditions for peak pressure stability in VOF simulations of dam break flow impact. J. Fluids Struct. 201662, 86–103. [Google Scholar] [CrossRef]
  19. Liu, L.; Sun, J.; Lin, B.; Lu, L. Building performance in dam break flow—an experimental sudy. Urban Water J. 201815, 251–258. [Google Scholar] [CrossRef]
  20. Armanini, A.; Larcher, M.; Odorizzi, M. Dynamic impact of a debris flow front against a vertical wall. In Proceedings of the 5th international conference on debris-flow hazards mitigation: Mechanics, prediction and assessment, Padua, Italy, 14–17 June 2011. [Google Scholar] [CrossRef]
  21. Hubbard, M.E.; Garcia Navarro, P. Flux difference splitting and the balancing of source terms and flux gradients. J. Comput. Phys. 2000165, 89–125. [Google Scholar] [CrossRef]
  22. Roe, P.L. A basis for upwind differencing of the two-dimensional unsteady Euler equations. In Numerical Methods in Fluids Dynamics II; Oxford University Press: Oxford, UK, 1986. [Google Scholar]
  23. Bradford, S.F.; Sander, B. Finite volume model for shallow water flooding of arbitrary topography. J. Hydraul. Eng. (ASCE) 2002128, 289–298. [Google Scholar] [CrossRef]
  24. Brufau, P.; Garica-Navarro, P. Two dimensional dam break flow simulation. Int. J. Numer. Meth. Fluids 200033, 35–57. [Google Scholar] [CrossRef]
  25. Hien, L.T.T. 2D Numerical Modeling of Dam-Break Flows with Application to Case Studies in Vietnam. Ph.D. Thesis, Brescia University, Brescia, Italy, 2014. [Google Scholar]
  26. Hien, L.T.T.; Tomirotti, M. Numerical modeling of dam break flows over complex topography. Case studies in Vietnam. In Proceedings of the 19th IAHR-APD Congress 2014, Hanoi, Vietnam, 21–24 September 2014; ISBN 978-604-82-1383-1. [Google Scholar]
  27. Flow-3D, Version 12.0; User Mannual; Flow Science Inc.: Santa Fe, NM, USA, 2020.
  28. Guney, M.S.; Tayfur, G.; Bombar, G.; Elci, S. Distorted physical model to study sudden partial dam break flow in an urban area. J. Hydraul. Eng. 2014140, 05014006. [Google Scholar] [CrossRef]
  29. Shige-eda, M.; Akiyama, J. Discussion and Closure to “Numerical and experimental study on two dimensional flood flows with and without structures” by Mirei Shige-eda and Juichiro Akiyama. J. Hydraul. Eng. 2005131, 336–337. [Google Scholar] [CrossRef]
  30. Ritter, A. Die Fortpflanzung der Wasserwelle (Generation of the water wave). Z. Ver. Dtsch. Ing. 189236, 947–954. [Google Scholar]
Figure 1. The bathymetry provided with the benchmark problem.

Performance Assessment of NAMI DANCE in Tsunami Evolution and Currents Using a Benchmark Problem

1Civil Engineering Department, Middle East Technical University, Ankara 06800, Turkey
2Ocean Engineering Department, University of Rhode Island, Narragansett, RI 02882, USA
3Civil Engineering Department, Middle East Technical University, Ankara 06800, Turkey
4Department of Applied Mathematics, Nizhny Novgorod State Technical University, Nizhny Novgorod 603950, Russia
*
Author to whom correspondence should be addressed.
Academic Editor: Richard P. Signell
J. Mar. Sci. Eng. 20164(3), 49; https://doi.org/10.3390/jmse4030049
Received: 5 July 2016 / Revised: 2 August 2016 / Accepted: 12 August 2016 / Published: 18 August 2016

Abstract

쓰나미 진화, 전파 및 침수의 수치 모델링은 현상에 관련된 수많은 매개 변수로 인해 복잡합니다. 쓰나미 모션을 해결하는 숫자 코드의 성능과 흐름 및 속도 패턴을 평가하는 것이 중요합니다. NAMI DANCE는 긴 파도 모델링을 위해 개발된 계산 도구입니다.

쓰나미 생성, 전파 및 침수 메커니즘의 수치 모델링 및 효율적인 시각화를 제공하고 쓰나미 매개 변수를 계산합니다. 긴 파도 이론에서, 물 입자의 수직 움직임은 압력 분포에 영향을 미치지 않습니다.

이러한 근사치와 소홀히 하는 수직 가속을 기반으로 질량 보존 및 모멘텀 방정식은 2차원 깊이 평균 방정식으로 줄어듭니다. NAMI DANCE는 유한차 계산 방법을 사용하여 긴 파도 문제에서 선형 및 비선형 형태의 깊이 평균 얕은 수식을 해결합니다.

이 연구에서 NAMI DANCE는 미국 포틀랜드에서 열린 2015 년 국립 쓰나미 위험 완화 프로그램 (NTHMP) 연례 회의에서 논의된 벤치 마크 문제에 적용됩니다.

벤치마크 문제는 하나의 독방 파도가 해양 섬 특징이 있는 삼각형 모양의 선반을 전파하는 일련의 실험을 특징으로 합니다. 이 문제는 섬 부근에서 상세한 무료 표면 고도 및 속도 의 타임 시리즈를 제공합니다. 결과를 비교한 결과, NAMI DANCE는 긴 파도 진화, 전파, 증폭 및 쓰나미 전류를 만족스럽게 예측할 수 있음을 보여주었습니다.

키워드: 수치 모델링;쓰나미 전류;깊이 평균 방정식;벤치마크,numerical modelingtsunami currentsdepth-averaged equationbenchmark

쓰나미는 해저 지진, 수중 산사태, 화산 폭발 또는 큰 운석 파업으로 인한 해저의 갑작스런 움직임에 의해 생성되는 큰 파도입니다. 쓰나미 파도는이 현상의 가장 파괴적인 매개 변수로 받아 들여진다; 그러나 큰 파도 움직임에 의해 트리거되는 전류는 경우에 따라 매우 치명적일 수 있습니다.

분지 공명 및 기하학적 증폭은 폐쇄 된 분지에서 쓰나미 영향의 지역 배율에 대한 두 가지 합리적으로 잘 이해된 메커니즘이며, 일반적으로 항구 또는 항구에서 쓰나미 위험 잠재력을 추정 할 때 조사 되는 메커니즘입니다. 반면에 전류에 대한 이해력과 예측 능력은부족하다[1]. 

이 연구는 수치 도구를 사용하여 쓰나미 진화, 전파 및 증폭뿐만 아니라 쓰나미 전류의 추정에 2 차원 깊이 평균 얕은 물 방정식의 충분성을 조사하는 것을 목표로; 즉 나미 댄스. 1970 년대 이후, 독방 파도는 일반적으로 실험 및 수학 연구에서, 쓰나미를 모델링하는 데 사용되었습니다[2]. 

이러한 점에서 수치 코드는 복잡한 목욕을 통해 단일 독방 파도의 진화와 전파에 초점을 맞춘 벤치마크 문제에 적용됩니다. 이 문제는 선반의 근해에 위치한 섬 특징이 있는 삼각형 모양의 선반을 전파할 때 단일 고독한 파도의 변형을 분석하는 일련의 실험을 설명합니다. 섬 부근에 형성되는 해류도 실험에서 조사된다.

이 연구에 사용된 벤치마크 문제는 미국 포틀랜드에서 개최된 2015 년 국립 쓰나미 위험 완화 프로그램 (NTHMP) 워크샵의 벤치마크 문제 #5.3]. 벤치마크 데이터와 수치 결과를 비교하여 2차원 깊이 평균 얕은 수식은 쓰나미 파도 진화와 해류에 대해 만족스러운 결과를 제공하므로 쓰나미 완화 전략을 결정하는 동안 사용하기에 충분한 도구임을 관찰합니다.

Figure 1. The bathymetry provided with the benchmark problem.
Figure 1. The bathymetry provided with the benchmark problem.
Figure 2. Model parameters: (a) bathymetry of the numerical model, NAMI DANCE; (b) incoming wave.
Figure 2. Model parameters: (a) bathymetry of the numerical model, NAMI DANCE; (b) incoming wave.
Figure 3. Comparison of free surface elevation (FSE) results: (a) X = 7.5 m and Y = 0.0 m at Gage 1; (b) X = 13.0 m and Y = 0.0 m at Gage 2; (c) X = 21.0 m and Y = 0.0 m at Gage 3; (d) X = 7.5 m and Y = 5.0 m at Gage 4; (e) X = 13.0 m and Y = 5.0 m at Gage 5; (f) X = 21.0 m and Y = 5.0 m at Gage 6; (g) X = 25.0 m and Y = 0.0 m at Gage 7; (h) X = 25.0 m and Y = 5.0 m at Gage 8. Black line represents benchmark data, red line represents numerical results.
Figure 3. Comparison of free surface elevation (FSE) results: (a) X = 7.5 m and Y = 0.0 m at Gage 1; (b) X = 13.0 m and Y = 0.0 m at Gage 2; (c) X = 21.0 m and Y = 0.0 m at Gage 3; (d) X = 7.5 m and Y = 5.0 m at Gage 4; (e) X = 13.0 m and Y = 5.0 m at Gage 5; (f) X = 21.0 m and Y = 5.0 m at Gage 6; (g) X = 25.0 m and Y = 0.0 m at Gage 7; (h) X = 25.0 m and Y = 5.0 m at Gage 8. Black line represents benchmark data, red line represents numerical results.
Figure 4. Comparison of results: (a) horizontal velocity in x-direction, U, recorded at X = 13.0 m, Y = 0.0 m and Z = 0.75 m at Gage 2; (b) horizontal velocity in y-direction, V, recorded at X = 13.0 m, Y = 0.0 m and Z = 0.75 m at Gage 2; (c) horizontal velocity in x-direction, U, recorded at X = 21.0 m, Y = −5.0 m and Z = 0.77 m at Gage 9; (d) horizontal velocity in y-direction, V, recorded at X = 21.0 m, Y = −5.0 m and Z = 0.77 m at Gage 9. Black line represents benchmark data, red line represents numerical results.
Figure 4. Comparison of results: (a) horizontal velocity in x-direction, U, recorded at X = 13.0 m, Y = 0.0 m and Z = 0.75 m at Gage 2; (b) horizontal velocity in y-direction, V, recorded at X = 13.0 m, Y = 0.0 m and Z = 0.75 m at Gage 2; (c) horizontal velocity in x-direction, U, recorded at X = 21.0 m, Y = −5.0 m and Z = 0.77 m at Gage 9; (d) horizontal velocity in y-direction, V, recorded at X = 21.0 m, Y = −5.0 m and Z = 0.77 m at Gage 9. Black line represents benchmark data, red line represents numerical results.

References

  1. Lynett, P.J.; Borrero, J.C.; Weiss, R.; Son, S.; Greer, D.; Renteria, W. Observations and modeling of tsunami-induced currents in ports and harbors. EPSL 2012327, 68–74. [Google Scholar]
  2. Madsen, P.A.; Fuhrman, D.R.; Schaffer, H.A. On the solitary wave paradigm for tsunamis. J. Geophys. Res. 2008113. [Google Scholar] [CrossRef]
  3. NTHMP Mapping & Modeling Benchmarking Workshop: Tsunami Currents. Benchmark #5. Available online: http://coastal.usc.edu/currents_workshop/problems/prob5.html (accessed on 2 August 2016).
  4. Onat, Y.; Yalciner, A.C. Initial stage of database development for tsunami warning system along Turkish coasts. Ocean Eng. 201374, 141–154. [Google Scholar] [CrossRef]
  5. Kian, R.; Yalciner, A.C.; Aytore, B.; Zaytsev, A. Wave Amplification and Resonance in Enclosed Basins; A Case Study in Haydarpasa Port of Istanbul. In Proceedings of the 2015 IEEE/OES Eleventh Current, Waves and Turbulence Measurement, St. Petersburg, VA, USA, 2–6 March 2015; Volume 11, pp. 1–7.
  6. Patel, V.M.; Dholakia, M.B.; Singh, A.P. Emergency preparedness in the case of Makran tsunami: A case study on tsunami risk visualization for the western parts of Gujarat, India. Geomat. Nat. Hazards Risk 20167, 826–842. [Google Scholar] [CrossRef]
  7. Yalciner, A.C.; Pelinovsky, E.; Zaytsev, A.; Kurkin, A.; Ozer, C.; Karakus, H.; Ozyurt, G. Modeling and visualization of tsunamis: Mediterranean examples. In Tsunami and Nonlinear Waves, 1st ed.; Kundu, A., Ed.; Springer: Berlin, Germany, 2007; pp. 273–283. [Google Scholar]
  8. Synolakis, C.E.; Bernard, E.N.; Titov, V.; Kanoglu, U.; Gonzalez, F. Validation and verification of tsunami numerical models. PAGEOPH 2008165, 2197–2228. [Google Scholar] [CrossRef]
  9. Yalciner, A.C.; Zaytsev, A.; Kanoglu, U.; Velioglu, D.; Dogan, G.G.; Kian, R.; Sharghivand, N.; Aytore, B. NTHMP Mapping and Modeling Benchmarking Workshop: Tsunami Currents. Available online: http://coastal.usc.edu/currents_workshop/presentations/Yalciner.pdf (accessed on 2 August 2016).
  10. Ozer, C.; Yalciner, A.C. Sensitivity study of hydrodynamic parameters during numerical simulations of tsunami inundation. PAGEOPH 2011168, 2083–2095. [Google Scholar]
  11. Sozdinler, C.O.; Yalciner, A.C.; Zaytsev, A. Investigation of tsunami hydrodynamic parameters in inundation zones with different structural layouts. PAGEOPH 2014172, 931–952. [Google Scholar] [CrossRef]
  12. Sozdinler, C.O.; Yalciner, A.C.; Zaytsev, A.; Suppasri, A.; Imamura, F. Investigation of hydrodynamic parameters and the effects of breakwaters during the 2011 Great East Japan Tsunami in Kamaishi Bay. PAGEOPH 2015172, 3473–3491. [Google Scholar] [CrossRef]
  13. Velioglu, D.; Kian, R.; Yalciner, A.C.; Zaytsev, A. Validation and Performance Comparison of Numerical Codes for Tsunami Inundation. In Proceedings of the 2015 American Geophysical Union Fall Meeting, San Francisco, CA, USA, 14–18 December 2015.
  14. Velioglu, D.; Kian, R.; Yalciner, A.C.; Zaytsev, A. Validation and Comparison of 2D and 3D Codes for Nearshore Motion of Long Waves Using Benchmark Problems. In Proceedings of the 2016 European Geosciences Union, Vienna, Austria, 17–22 April 2016.
  15. Dilmen, D.I.; Kemec, S.; Yalciner, A.C.; Düzgün, S.; Zaytsev, A. Development of a tsunami inundation map in detecting tsunami risk in Gulf of Fethiye, Turkey. PAGEOPH 2015172. [Google Scholar] [CrossRef]
  16. Heidarzadeh, M.; Krastel, S.; Yalciner, A.C. The state-of-the-art numerical tools for modeling landslide tsunamis: A short review. In Submarine Mass Movements and Their Consequences, 6th ed.; Sebastian, K., Jan-Hinrich, B., David, V., Michael, S., Christian, B., Roger, U., Jason, C., Katrin, H., Michael, S., Carl, B.H., Eds.; Springer: Bern, Switzerland, 2013; Volume 37, pp. 483–495. [Google Scholar]
  17. Yalciner, A.C.; Gülkan, P.; Dilmen, D. I.; Aytore, B.; Ayca, A.; Insel, I.; Zaytsev, A. Evaluation of tsunami scenarios for western Peloponnese, Greece. Boll. Geofis. Teor. Appl. 201455, 485–500. [Google Scholar]
  18. Zahibo, N.; Pelinovsky, E.; Kurkin, A.; Kozelkov, A. Estimation of far-field tsunami potential for the Caribbean Coast based on numerical simulation. Sci. Tsunami Hazards 200321, 202–222. [Google Scholar]
  19. Swigler, D.T. Laboratory Study Investigating the Three-dimensıonal Turbulence and Kinematic Properties Associated with a Breaking Solitary Wave. Master’s Thesis, Texas A&M University, College Station, TX, USA, August 2009. [Google Scholar]
  20. National Tsunami Hazard Mitigation Program. Proceedings and Results of the 2011 NTHMP Model Benchmarking Workshop. Available online: http://nws.weather.gov/nthmp/documents/nthmpWorkshopProcMerged.pdf (accessed on 21 July 2016).
圖1. 1 南海孤立內波空間分布圖(Hsu et al., 2000)

Numerical Modeling on Internal Solitary Wave propagation over an obstacle using Flow-3D

Keyword: Internal solitary waves, Numerical, Flow-3D, Computational Fluid Dynamics

연구자 : Yu-Ren Chen
지도교수 : Dr John R C Hsu
June 2012

기술과 수치 알고리즘의 발전으로 파도가 해양이나 항만 구조물에 미치는 영향에 대한 많은 연구가 개발되었으며,보다 정확한 결과를 얻기 위해 고효율 수치 계산 소프트웨어를 사용할 수 있습니다. 현재 내부 파 생성, 전송, 파동의 물리적 메커니즘은 국내외 해양 분야에서 중요한 연구 주제 중 하나입니다.

이 연구는 FLOW-3D 전산 유체 역학 (Computational Fluid Dynamics, CFD) 소프트웨어를 사용하여 상층의 담수와 하층의 담수를 시뮬레이션합니다. 바닷물의 밀도 계층화 유체는 중력 혼합 붕괴 방식을 사용하여 내부 파도를 생성하고 긴 경사와 같은 일반적인 장애물을 통해 파형 진화 및 유동장 분포를 탐구합니다.

짧은 플랫폼 사다리꼴 경사와 이등변 삼각형. 이 기사에서는 또한 소프트웨어 작동 설정과 FLOW-3D를 내부 파 실험에 적용하는 방법을 소개하고, 이전 실험 조건과 결과를 참조하여 내부 파 전송 과정을 시뮬레이션합니다. 시뮬레이션 결과는 실험 데이터를 확인하고 첫 번째 분석을 시뮬레이션합니다.

중력 붕괴 방식의 게이트의 개방 속도가 내부 파의 전송 시간 및 진폭에 미치는 영향; 시뮬레이션 결과는 게이트 개방 속도가 빠르고 내부 파의 진폭이 크고 전송 속도가 빠릅니다. ; 반대로 게이트 개방 속도가 느리면 내부 파의 진폭이 작고 전송 속도가 느리지 만 둘 다 비선형 비례 관계.

이 연구는 또한 다양한 장애물 (긴 기울기, 사다리꼴 기울기가있는 짧은 플랫폼, 이등변 삼각형)을 통한 내부 고독 파의 전송 과정을 시뮬레이션하고 단일 장애물을 통과하는 내부 파도의 파형 진화, 와류 및 유동장 변화를 논의합니다.

연구를 통해 우리가 매우 미세한 그리드를 사용하고 수치 시뮬레이션의 그래픽 출력을 열심히 분석 할 수 있다면 실험실 실험 관찰보다 내부 고독 파의 전송 특성을 더 잘 이해할 수 있다고 믿습니다.

요약

서로 다른 특성을 가진 두 유체의 계면에있는 파동을 계면 파라고합니다. 바다에서는 표층의 기압 변화에 의해 형성된 바람 장이 공기와 바다의 경계 파인 해면에 불어 올 때 변동을 일으킨다. 기체 또는 유체의 밀도 층화가 발생할 때 외부 힘 (예 : 바람, 압력, 파도 및 조류, 중력 등)에 의해 교란되면 내부 파도라고하는 경계면에서 변동이 발생할 수 있으므로 내부 파도가 발생할 수 있습니다. 웨이브는 밀도가 다른 층화 된 유체의 웨이브 현상입니다.

대기의 내부 파도와 같이 일상 생활에서 볼 수있는 내부 파도는 특히 오후 또는 비가 내리기 전에 깊고 얕은 altocumulus 구름 층으로 하늘에 자주 나타납니다. 대기 중의 내부 파의 움직임은 공기의 흐름에 영향을 주어 기류를 상승시키고 공기 중의 수증기가 물방울로 응축되어 구름이되도록합니다.

반대로 기류가 가라 앉으면 수증기가 응결이 쉽지 않습니다. 구름이 있어도 내부의 파도가 응결하기 어렵습니다. 소산되어 루버와 같은 altocumulus 구름을 형성합니다. 안정된 밀도와 층화 상태의 자연 수체는 외부 세계에 의해 교란 될 때 내부 파동 운동을 겪게됩니다.

예를 들어, 밀도가 안정되고 층화가 분명한 호수에서 바람 장은 수면에 파도에서 파생 된 내부 파동을 일으켜 물의 질량이 전달되고 호수 가장자리로 물이 축적되어 수위가 높아집니다. 위치 에너지를 형성하는 축적 영역; 수역이 가라 앉기 시작하면 위치 에너지를 운동 에너지로 변환하고 남미 콜롬비아의 Babine Lake의 내부 파동 거동과 같은 내부 파동 운동을 생성 할 수도 있습니다 (Farmer, 1978). ). 염분, 밀도 또는 온도가 안정된 바다에서는 조수와 지형의 영향으로 수역이 행성의 중력에 따라 움직입니다.

격렬한 기복이있는 지형을 통과 할 때 내부 파동이 발생합니다. ; 중국 해에서 발견되는 남쪽 내부 파도에서와 같이 (Hsu et al., 2000). 파동은 심해에서 얕은 물로 전달되며, 얕아 짐, 깨짐, 혼합, 소용돌이, 굴절, 회절 및 반사가있을 것입니다. 내부 파 전달은 일종의 파동이기 때문에 위에서 언급 한 파동 특성도 갖습니다.

해양 내부 파도는 길이가 수백 미터에서 수십 킬로미터에 이르는 광범위한 파장을 가지고 있으며,주기는 몇 분 정도 빠르며 수십 시간 정도 느리며 진폭은 몇 미터에서 수백 미터. 해양 내부 파도가 움직일 때 층화 위와 아래의 물 흐름 방향이 반대가되어 현재 전단 작용으로 인해 층화 경계면에서 큰 비틀림 힘이 발생합니다.

바다에 기초 말뚝과 같은 구조물이있는 경우 석유 시추 플랫폼의 고정 케이블은 큰 비틀림을 견딜 수 없어 파손될 가능성이 매우 높습니다 (Bole et al. 1994). 빽빽한 클라인 경계 근처에서 항해하는 잠수함이 해양 내부 파도 활동을 만나게되면 내부 파도에 의한 상승 전류로 인해 잠수함이 해저에 수면에 닿거나 충돌하여 잠수함이 손상 될 수 있습니다.

그러나 바다의 내부 파는 바람직하지 않으며 매우 중요한 역할을합니다. 예를 들어, 내부 파가 심해 지역에서 근해 대륙붕으로 전달되면 상하수 체가 교환됩니다. 해저에 영양분을 운반합니다. 선반 가장자리까지 생물학적 성장을 촉진하고 해당 지역의 생태 환경을 조절하며 (Osborne and Bruch et al., 1980; Sandstorm and Elliot et al., 1984) 어업 자원을 풍부하게합니다.

위에서 언급 한 항목 외에도 해저에 대한 케이블 및 파이프 라인, 수중 음파 탐지기, 해양 생물 환경, 군사 활동 등이 해양 내부 파도의 영향에 포함되므로 해양 내부 파도에 대한 연구가 매우 중요합니다.

최근 내부 파를 연구하는 방법에는 분석 이론 도출, 현장 조사 및 관찰, 실험실 실험 분석이 포함됩니다. 그러나 과학 기술의 급속한 발전, 발전과 발전, 컴퓨터의 대중화, 수치 계산 방법의 진화로 해양 공학과 관련된 많은 파동 효과는 일반적으로 수치 시뮬레이션 방법으로 해결됩니다.

또한 수치 연산 방법의 비용이 현장 조사 관측 및 실험실 실험 해석보다 저렴하고 시뮬레이션 결과를 더 빨리 얻을 수 있기 때문에 본 논문에서는 전산 유체 역학 (전산 유체 역학, 참조)의 FLOW-를 선정 하였다. 3D 소프트웨어는 내부 파 생성, 전송, 장애물 통과, 점차 소멸하는 움직임 과정을 시뮬레이션하고, 내부 파의 변화 과정을 분석하고 비교하기 위해 이전 실험실 모델 실험을 참조합니다.

圖1. 1  南海孤立內波空間分布圖(Hsu et al., 2000)
圖1. 1 南海孤立內波空間分布圖(Hsu et al., 2000)
圖1. 2  障礙高度與分層流體厚度關係之示意圖
圖1. 2 障礙高度與分層流體厚度關係之示意圖
圖3. 1 下沉型內孤立波通過梯形障礙的實驗配置圖(鄭明宏,2011)
圖3. 1 下沉型內孤立波通過梯形障礙的實驗配置圖(鄭明宏,2011)
圖3. 3  實驗室下沉型內孤立波經過13°斜坡梯形障礙物的連續組圖(鄭明宏,2011)
圖3. 3 實驗室下沉型內孤立波經過13°斜坡梯形障礙物的連續組圖(鄭明宏,2011)
圖3. 3 (a) 實驗室下沉型內孤立波(鄭明宏,2011;θ=13°,T = t0 = 42 s)
圖3. 3 (a) 實驗室下沉型內孤立波(鄭明宏,2011;θ=13°,T = t0 = 42 s)
圖3. 5 比較實驗室(上圖)內孤立波(圖3. 3 (a))與FLOW-3D模擬(下圖)的傳遞波形(θ=13°,t = 42 s)
圖3. 5 比較實驗室(上圖)內孤立波(圖3. 3 (a))與FLOW-3D模擬(下圖)的傳遞波形(θ=13°,t = 42 s)
圖4. 6閘門開啟速率0.14 m/s之等密度線及流場
圖4. 6閘門開啟速率0.14 m/s之等密度線及流場

圖4. 53 內波在三角形前坡反轉為順時針渦流,後坡面上形成逆時針渦流(t = 63 s)
圖4. 53 內波在三角形前坡反轉為順時針渦流,後坡面上形成逆時針渦流(t = 63 s)

Reference

Apel, J.R., Holbrook, J.R, Tsai, J. and Liu, A.K. (1985). The Sulu Sea internal soliton experiment. J. Phys. Oceanography, 15(12): 1625-1651. Ariyaratnam, J. (1998). Investigation of slope stability under internal wave action. B.Eng. (Hons.) thesis, Dept. of Environmental Eng., University of Western Australia, Australia. Baines, P.G. (1983). Tidal motion in submarine canyons – a laboratory experiment. J. Physical Oceanography, 13: 310-328. Benjamin, T.B. (1966). Internal waves of finite amplitude and permanent form. J. Fluid Mech., 25: 241-270. Bole, J.B., Ebbesmeyer, J.J. and Romea, R.D. (1994). Soliton currents in South China Sea: measurements and theoretical modelling. Proc. 26th Annual Offshore Tech. Conf., Houston, Texas. 367-375. Burnside, W. (1889). On the small wave-motions of a heterogeneous fluid under gravity. Proc. Lond., Math. Soc., (1) xx, 392-397. Chen C.Y., J.R-C. Hsu, H.H. Chen, C.F. Kuo and Cheng M.H (2007). Laboratory observations on internal solitary wave evolution on steep and inverse uniform slopes. Ocean Engineering, 34: 157-170. Cheng M.H., J.R-C. Hsu, C.Y. Chen (2005). Numerical model for internal solitary wave evolution on impermeable variable seabad, Proc.27th Ocean Eng, pp.355-359. Choi, W. and Camassa, R. (1996). Weakly nonlinear internal waves in a two-fluid system. J. Fluid Mech., 313: 83-103. Ebbesmeyer, C.C., and Romea, R.D. (1992). Final design parameters for solitons at selected locations in South China Sea. Final and supplementary reports prepared for Amoco Production Company, 209pp. plus appendices. Ekman, V. M., (1904). “On dead-water, Norwegian North Polar Expedition”, 1893-1896. Scientific Results, 5(15):1-150. Farmer, D.M. (1978). Observation of long nonlinear internal waves in a lake. J. Phys. Oceanography, 8(1): 63-73. Garret, C. and Munk, W. (1972). Space-time scales of internal waves. Geophys. Fluid Dyn., 3: 225-264. Gill, A.E. (1982). Atmosphere-Ocean Dynamics. International Geophysical Series, Vol. 30, San Diego, CA: Academic Press. Harleman, D.R.F. (1961). Stratified flow. Ch. 26 in Handbook of Fluid Dynamics (ed., V. Streeter), NY: McGraw-Hill, (26): 1-21. Helfrich, K.R. (1992). Internal solitary wave breaking and run-up on a uniform slope. J. Fluid Mech., 243: 133-154.

Helfrich, K.R. and Melville, W.K. (1986). On long nonlinear internal waves over slope-shelf topography. J. Fluid Mech., 167: 285-308. Honji, H., Matsunaga, N., Sugihara, Y. and Sakai, K. (1995). Experimental observation of interanl symmetric solitary waves in a two-layer fluid. Fluid Dynamics Research, 15 (2): 89-102. Hsu, M.K., Liu, A.K., and Liu, C. (2000). A study of internal waves in the China Sea and Yellow Sea using SAR. Continental Shelf Research, 20: 389-410. Johns, K. (1999). Interaction of an internal wave with a submerged sill in a two-layer fluid. B.Eng. (Hons.) thesis, Dept. of Environmental Eng., University of Western Australia, Australia Kao, T.W., Pan, F.S. and Renouard, D. (1985). Internal solitions on the pycnocline: generation, propagation, shoaling and breaking over a slope. J. Fluid Mech. 159: 19-53. Koop, C.G. and Butler, G. (1981). An investigation of internal solitary waves in a two-fluid system. J. Fluid Mech., 112: 225-251. Lin, T.W. (2001). A study on internal waves characteristics in north of South China Sea, Master Thesis, Institute of Oceanography, National Taiwan Univ., Taiwan. (In Chinese). Lynett, P., Wu, T.-R. and Liu, P. L.-F. (2002), Modeling wave runup with depth-integrated equations, Coastal Engineering, Vol. 46, pp. 89-107. Ming-Hung Cheng,John R.-C. Hsu, Chen-Yuan Chen and Cheng-Wu Chen (2009). Modelling the propagation of an internal solitary wave across double ridges and a shelf-slope.Environ Fluid Mech,9:321–340. Ming-Hung Cheng and John R.C. Hsu (2011). Effect of frontal slope on waveform evolution of a depression interfacial solitary wave across a trapezoidal obstacle. Ocean Engineering. Matsuno, Y. (1993). A unified theory of nonlinear wave propagation in two-layer fluid systems. J. Phys. Soc. Japan, 62: 1902-1916. Michallet, H. and Barthelemy, E. (1998). Experimental study of interfacial solitary waves. J. Fluid Mech., 366: 159-177. Muller, P. and X. Liu (2000). Scattering of internal waves at finite topography in two dimensions. Part I: Theory and case studies, J. Phys. Oceanogr., 30: 532-549 Nagashima, H. (1971). Reflection and breaking of internal waves on a sloping beach. J. Oceanographical Soc. Japan, 27(1): 1-6. Nansen, F. (1902). The oceanography of the north polar basin. Sci. Results, Norwegian North Polar Expedition 1893-1896, 3: 9. Osborne, A.R. and Burch, T.L. (1980). Internal solitons in the Andaman Sea. Science, 208 (43): 451-460

82 Russell, J.S. (1844). On waves. Report of the 14th Meeting of the British Association for the Advancement of Science, York, 311-390. Sandstrom, H. and Elliot J. A. (1984). Internal tide and solitons on the Scotian Shelf: a nutrient pump at work. Journal of Geophysical Research, 89 (C4): 6415-6428. Stokes G.G. (1847). On the Theory of Oscillatory Waves. Transactions of the Cambridge Philosophical Society, 8: 441–455. Strutt, J. W., Lord Rayleigh. (1883). Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density.Proceedings of the London mathematical society, 8: pp. 170-177. Sveen, J.K., Guo, Y., Davies, P.A. and Grue, J. (2002). On the breaking of internal solitary waves at a ridge. J. Fluid Mech., 469 (25): 161-188. Vlasenko, V., and Hutter, K. (2002). Numerical experiments on the breaking of solitary internal waves over a slope-shelf topography. J. of Physical Oceanography, 32(6), pp.1779-1793. Wessels F. and Hutter K. (1996). Interaction of internal waves with a topographic sill in a two-layered fluid. J. Phys. Oceanogr , 26: 5-20

ANSI/HI 9.8 Pump Intake Design

Hydraulic Jump in a Trench Type Pump Sump

트렌치 형 펌프 배수 조의 유압 점프

이 기사는 Ibis Group의 대표인 Steve Saunders가 기고했습니다.

유압 점프는 개방형 채널 애플리케이션으로 작업하는 사람들에게 친숙한 흐름 현상입니다. Wikipedia는 수력 점프를 “개방형 채널 흐름이 초-임계에서 아임계로 갑자기 변환되는 조건”으로 정의합니다. 점프가 발생하는 위치에서 속도 헤드가 수면 상승으로 거래되는 것을 관찰 할 수 있습니다. 방수로와 같은 흐름 제어 응용 분야에서 수압 점프는 침식을 완화하기 위해 에너지를 소산하는 수단으로 의도적으로 설정됩니다. 또한 레크리에이션 목적으로 사용됩니다. 유압 점프로 생성된 정상 파도는 어떤 바다에서든 수천 마일 떨어진 서핑 공원에서 타는 방법을 서퍼를 훈련시키는데 사용됩니다. 유압 점프의 새로운 응용 분야는 점프의 에너지 전달이 다시 중단되고 정상적인 펌핑 작업 중에 침전된 고형물을 제거하는 자가 세척 트렌치 유형 펌프 섬프(sump)입니다.

트렌치 유형 집수 펌프 시뮬레이션
FLOW-3D는 유압 점프 시뮬레이션에서 신뢰할 수 있는 도구로 입증되었으며 자가 세척 트렌치 유형 펌프 섬프의 설계 및 시연에 사용되었습니다. 트렌치 형 펌프 섬프는 펌프 흡입 라인이 있는 좁은 채널로 구성됩니다. 일반적인 응용 분야는 들어오는 물에서 모래와 자갈을 걸러내는 입구 스크린이 없는 빗물 수집입니다. 아래 회로도에 예가 나와 있습니다.

ANSI/HI 9.8 Pump Intake Design
ANSI/HI 9.8 Pump Intake Design

이 수치는 ANSI / HI 9.8 펌프 흡기 설계 매뉴얼에서 발췌한 것이며 4 개의 펌프가 설치된 섬프의 평면도 및 입면도를 보여줍니다. 유입 암거, 웅덩이 바닥 및 펌프 흡입 바닥을 벗어난 높이의 배열은 이 설계 유형의 자체 청소 기능에 매우 중요합니다. 유입 암거는 최소 작동 웅덩이 수위보다 높은 고도에 있습니다. 또한 유입단의 ​​트렌치 벽은 Ogee 모양입니다. 마지막으로, 트렌치의 맨 끝에 있는 펌프 흡입 벨은 상류 펌프의 절반 높이에 설정됩니다.

Designing for Storm Events

폭풍이 닥친 후 모래와 자갈이 웅덩이 바닥에 쌓입니다. 그들은 점진적인 유압 점프를 통해 다시 매달리고 빠져 나갑니다. 청소 주기 동안 물은 유입 암거를 통해 유입되는 것보다 더 빠른 속도로 트렌치의 맨 끝에 있는 하부 펌프에 의해 배출됩니다.

이 시퀀스 동안 유압 점프는 두 가지 중요한 역할을 수행합니다. 점프 업스트림의 초임계 부분은 섬프 바닥의 모래와 자갈을 휘감아 펌핑이 되도록 다시 일시 중단합니다. 애니메이션의 색상 스케일을 보면 ogee 바닥의 수색 속도가 약 9ft/sec에 가깝다는 것을 알 수 있습니다. 한편, 점프 하류의 계단식 수면 상승은 하단 펌프에 충분한 잠수를 제공하여 섬프가 펌핑 될 때까지 계속 작동합니다.

물이 최소 정상 작동 수준 아래로 떨어지면 유입이 Ogee 모양의 벽 아래로 가속되어 궁극적으로 초임계가됩니다. 섬프의 수위가 바닥에 가까워지면 수압 점프가 형성되고 하단 원단 펌프가 흡입력을 잃을 때까지 섬프를 따라 진행됩니다. 아래 애니메이션에서 이런 일이 일어나는 것을 관찰 할 수 있습니다.

The Magnolia Storm Water Pumping Station

이 자체 세척 섬프 응용 분야에 FLOW-3D를 사용하면 트렌치 형상을 쉽게 조정하여 유압 점프 동작을 최적화 할 수 있습니다. 텍사스 엘파소에있는 Magnolia Storm Water Pumping Station은 FLOW-3D가 설계 및 평가 도구로 사용 된 예입니다. 2016 년에 시운전 된 Magnolia Storm Water Pumping Station은 폭우시 고속도로 10 번의 홍수를 방지하기 위해 건설되었습니다.

Magnolia 스테이션은자가 세척 트렌치 유형 섬프에 3 개의 대형 수직 터빈 펌프로 구성됩니다. 섬프 설계 과정에서 FLOW-3D를 사용하여 몇 가지 기하학적 변형을 평가하여 자체 세척 기능을 통해 펌프 작동 효율성 및 유지 보수 용이성에 이상적인 구성에 도달했습니다.

Evaluation of the Wind Effects on the Iron-Ore Stock Pile

Energy

Energy

전 세계 에너지 부문의 엔지니어는 전산 유체 역학(CFD)을 통해 해결책을 찾기 위해 광범위한 프로세스에서 매일 복잡한 설계 문제에 직면합니다. 특히 자유 표면 흐름과 관련이 높은 이러한 문제의 대부분은 FLOW-3D가 매우 정확한 분석을 제공하여 문제 해결에 적합합니다.

  • 공해에서 컨테이너 내부의 연료 또는화물 슬로싱 / Fuel or cargo sloshing inside containers on the high seas
  • 해양 플랫폼에 대한 파도 효과 / Wave effects on offshore platforms
  • 6 자유도 모션을 받는 분리 장치의 성능 최적화 / Performance optimization for separation devices undergoing 6 DOF motion
  • 파동 에너지 포착 장치 / Design of devices to capture energy from waves

Energy Case Studies

천연자원이 계속 감소함에 따라, 대체 자원과 방법을 탐구하고 가능한 한 효과적으로 현재 공급량을 사용하고 있습니다. 엔지니어는 사고를 예방하고 채굴 및 기타 에너지 수확 기법으로 인한 환경적 영향을 평가하기 위해 FLOW-3D를 사용합니다.

Tailing Breach Simulation – CFD Analysis with FLOW-3D

점성이 높은 유체, 비 뉴턴 흐름, 슬러리 또는 심지어 세분화 된 흐름의 형태를 취할 수있는 많은 채광 응용 프로그램의 잔여 물인 테일링은 악명 높은 시뮬레이션 전제를 제공합니다. FLOW-3D  는 비 뉴턴 유체, 슬러리 및 입상 흐름에 대한 특수 모델을 포함하여 이러한 분석을 수행하는 데 필요한 모든 도구를 제공합니다. FLOW-3D 의 자유 표면 유동 모델링 기능 과 결합되어  이러한 어렵고 환경 적으로 민감한 문제에 대한 탁월한 모델링 솔루션을 제공합니다.

관련 응용 분야에는 바람 강제 분석에 따른 광석 비축 더미 먼지 드리프트가 포함되며, 여기서 FLOW-3D 의 드리프트 플럭스 모델을 통해 엔지니어는 광석 침착 및 유입 패턴과 개선 솔루션의 효과를 연구 할 수 있습니다.

액화와 기계적 방해가 물과 같은 뉴턴 흐름과는 대조적으로 입자 흐름의 매우 독특한 속성 인 결국 저절로 멈추는 위반의 동적 특징의 일부라는 점에 유의하십시오.

오일 및 가스 분리기

FLOW-3D  는 기름과 물과 같은 혼합 불가능한 유체를 모델링 할 수 있으며 개방 된 환경 (주변 공기)과 관련된 구성 요소 간의 뚜렷한 인터페이스를 정확하게 추적 할 수 있습니다. 유체는 전체 도메인에 영향을 미치는 역학으로 인해 자유롭게 혼합 될 수 있습니다. 시간이 지남에 따라 유체는 연속 상과 분산 상 간의 드리프트 관계에 따라 다시 분리됩니다. 중력 분리기의 성능은 CFD 모델링을 통해 향상 될 수 있습니다.

  • 기체 및 액체 흐름의 균일성을 개선하고 파도에 의한 슬로싱으로 인한 오일과 물의 혼합을 방지하기 위해 용기 입구 구성을 개발합니다.
  • 유압 효율 및 분리 성능에 대한 내부 장비의 영향을 결정합니다.
  • 작동 조건 변화의 영향 측정
  • 소규모 현상 (다상 흐름, 방울, 입자, 기포)을 정확하게 모델링

생산 파이프 | Production Pipes

생산에 사용되는 공정 파이프의 청소 과정에서 유체가 위로 흘러도 고밀도 입자가 침전될 수 있습니다. 침전 입자를 포착하도록 장치를 설계 할 수 있습니다. 파이프 중앙에 있는 “버킷”이 그러한 잠재적 장치중 하나 입니다. 흐름 변위로 인해 버킷 외부의 상류 속도는 고밀도 입자에 대한 침전 속도보다 높으며 버킷 내부에 모여 있습니다. 표시된 디자인에서 버킷 주변의 상향 유체 속도는 입자 안정화 속도보다 높습니다. 이로 인해 입자가 버킷과 파이프 벽 사이의 틈새를 통해 빠져 나갈 수 없습니다. 따라서 시뮬레이션된 입자는 버킷을 통과하여 아래에 정착하지 않습니다.

파동 에너지 장치 모델링 | Modeling Wave Energy Devices

포인트 흡수 장치 | Point Absorber Devices

이 시뮬레이션은 상단에 부력이있는 구형 구조가있는 점 흡수 장치를 보여 주며, 들어오는 파도의 볏과 골과 함께 위아래로 이동합니다. FLOW-3D 의 움직이는 물체 모델은 x 또는 y 방향으로의 움직임을 제한하면서 z 방향으로 결합 된 움직임을 허용하는 데 사용됩니다. 진폭 5m, 파장 100m의 스톡 스파를 사용했다. RNG 모델은 파도가 점 흡수 장치와 상호 작용할 때 발생하는 난류를 포착하는 데 사용되었습니다. 예상대로 많은 난류 운동 에너지가 장치 근처에서 생성됩니다. 플롯은 난류로 인해 장치 근처의 복잡한 속도 장의 진화로 인해 질량 중심의 불규칙한 순환 운동을 보여줍니다.

다중 플랩, 하단 경첩 파동 에너지 변환기 | Multi-Flap, Bottom-Hinged Wave Energy Converter

진동하는 플랩은 바다의 파도에서 에너지를 추출하여 기계 에너지로 변환합니다. Arm은 물결에 반응하여 피벗된 조인트에 장착된 진자로 진동합니다. 플랩을 배열로 구성하여 다중 플랩 파동 에너지 변환기를 만들 수 있습니다. 아래 상단에 표시된 CFD 시뮬레이션에서 3 개의 플랩 배열이 시뮬레이션됩니다. 모든 플랩은 바닥에 경첩이 달려 있으며 폭 15m x 높이 10m x 두께 2m입니다. 어레이는 30m 깊이에서 10 초의 주파수로 4m 진폭파에서 작동합니다. 시뮬레이션은 중앙 평면을 따라 복잡한 속도 등 가면을 보여줍니다. 이는 한 플랩이 어레이 내의 다른 플랩에 미치는 영향을 연구하는 데 중요합니다. 3 개의 플랩이 유사한 동적 동작으로 시작하는 동안 플랩의 상호 작용 효과는 곧 동작을 위상에서 벗어납니다. 유사한 플랩 에너지 변환기가 오른쪽 하단에 표시됩니다. 이 시뮬레이션에서 플랩은 가장 낮은 지점에서 물에 완전히 잠 깁니다. 이러한 에너지 변환기를 Surface Piercing 플랩 에너지 변환기라고합니다. 이 두 시뮬레이션 예제는 모두 미네르바 역학 .

진동 수주 | Oscillating Water Column

진동하는 수주는 부분적으로 잠긴 중공 구조입니다. 그것은 물의 기둥 위에 공기 기둥을 둘러싸고 수면 아래의 바다로 열려 있습니다. 파도는 물 기둥을 상승 및 하강시키고, 차례로 공기 기둥을 압축 및 감압합니다. 이 갇힌 공기는 일반적으로 기류의 방향에 관계없이 회전 할 수 있는 터빈을 통해 대기로 흐르게 됩니다. 터빈의 회전은 전기를 생성하는 데 사용됩니다.

아래의 CFD 시뮬레이션은 진동하는 수주를 보여줍니다. FLOW-3D에서 포착한 물리학을 강조하기 위해 중공 구조에서 물기둥이 상승 및 하강하는 부분만 모델링  합니다. 시뮬레이션은 다른 파형 생성 선택을 제외하고 유사한 결과를 전달합니다. 아래의 시뮬레이션은 웨이브 유형 경계 조건을 사용하는 반면 그 아래의 시뮬레이션은  움직이는 물체 모델  을 사용하여 실험실에서 수행한 것처럼 차례로 웨이브를 생성하는 움직이는 플런저를 생성합니다. 각 시뮬레이션에 대해 속이 빈 구조의 압력 플롯이 표시됩니다. 결국 그 압력에 기초하여 터빈이 회전 운동으로 설정되기 때문에 챔버에서 얼마나 많은 압력이 생성되는지 아는 것이 중요합니다.

사례 연구

eadership-in-energy-and-environmental-design

Architects Achieve LEED Certification in Sustainable Buildings

LEED (Leadership in Energy and Environmental Design)는 제 3자가 친환경 건축물 인증을 제공하는 자발적 인증 시스템입니다.

FLOW-3D는 보고타(콜롬비아)의 사무실 건물에서 “IEQ-Credit2–환기 증가”라는 신뢰를 얻는 데 큰 도움을 주었습니다. 이러한 인정을 받기 위해서는 실외 공기가 ASHRAE의 표준 비율인 30%를 초과한다는 것을 증명해야만 합니다. 이 건물에서 실외 공기는 태양 광선에 의해, 가열되는 지붕 위의 2개의 유리 굴뚝에 의해 발생되는 온도 차이에 의해 발생하는 열 부력의 영향으로 제공됩니다. 이것은 바람이 불지 않는 조건에서 이루어져야 합니다.

Comparing HVAC System Designs

최근 프로젝트에서 Tecsult의 HVAC(난방, 냉방 및 환기)시스템 엔지니어는 강력한 에어컨 시스템의 두 가지 다른 구성을 고려해야 했고 노동자들에게 어떤 것이 가장 쾌적함을 제공하는지 보여주기를 원했습니다. FLOW-3D는 대체 설계를 시뮬레이션하고 비교하는 데 사용되었습니다.

이 발전소는 대형(길이 90m, 너비 33m, 높이 26m)건물로 변압기, 전력선, 조명 등 열 발생 장비를 갖추고 있습니다. 에어컨 시스템의 목적은 건물 내 최대 온도를 35ºC로 제한하는 것입니다. 디퓨저가 하부 레벨에 위치하고 천장 근처의 환기구가 있기 때문에 천장 근처에서 최대 공기 온도가 발생하고 바닥 레벨은 반드시 몇도 더 낮습니다.

Modeling velocity of debris types

Debris Transport in a Nuclear Reactor Containment Building

이 기사는 FLOW-3D가 원자력 시설에서 봉쇄 시설의 성능을 모델링하는데 사용된 방법을 설명하며, Alion Science and Technology의 Tim Sande & Joe Tezak이 기고 한 바 있습니다.

가압수형 원자로 원자력 발전소에서 원자로 노심을 통해 순환되는 물은 약 2,080 psi 및 585°F의 압력과 온도로 유지되는 1차 배관 시스템에 밀폐됩니다. 수압이 높기 때문에 배관이 파손되면 격납건물 내에 여러 가지 이물질 유형이 생성될 수 있습니다. 이는 절연재가 장비와 균열 주변 영역의 배관에서 떨어져 나가기 때문에 발생합니다. 생성될 수 있는 다양한 유형의 이물질의 일반적인 예가 나와 있습니다(오른쪽).

Evaluation of the Wind Effects on the Iron-Ore Stock Pile

Evaluation of the Wind Effects on the Iron-Ore Stock Pile

바람이 개방형 골재 저장소에 미치는 영향은 전 세계적으로 환경 문제가 되고 있습니다. 2.7km철골 저장소 부지에서 이런 문제가 관찰되었습니다. 이 시설은 철도 운송차량를 통해 광석을 공급받는데, 이 운송차량은 자동 덤프에 의해 비워집니다. 그런 다음 이 광석은 일련의 컨베이어와 이송 지점을 통과하여 저장 장소중 하나로 운송됩니다. 비산먼지 배출은 풍력이 비축된량에 미치는 영향의 결과로 관찰된 결과입니다.

관련 기술자료


Figure 1. Photorealistic view of an inclined axis TAST (photo A. Stergiopoulou).

그리스 수로의 작은 수력 전위를 활용하는 관형 아르키메데스 스크류 터빈의 CFD 시뮬레이션

CFD Simulations of Tubular Archimedean Screw Turbines Harnessing the Small Hydropotential of Greek Watercourses Alkistis Stergiopoulou1, Vassilios Stergiopoulos21Institut für Wasserwirtschaft, ...
더 보기
Fig. 1. Nysted Offshore Wind Farm

FLOW-3D 모형을 이용한 해상풍력기초 세굴현상 분석

박영진1, 김태원2*1 서일대학교 토목공학과, 2 (주)지티이 Analysis of Scour Phenomenon around Offshore Wind Foundation using Flow-3D Mode Abstract 국내․외에서 다양한 ...
더 보기
Figure 2.1. Test Setup.The test setup consists of a clear plastic scale model tank attached to a rigid aluminum frame by three multi-axis load cells driven by a position-controlled servo hydraulic system.(Data acquisition cabling removed for clarity).

Coupled Simulation of Vehicle Dynamics and Tank Slosh. Phase 1 Report. Testing and Validation of Tank Slosh Analysis

Prepared byGlenn R. WendelSteven T. GreenRussell C. Burkey Abstract: 차량 동력학의 컴퓨터 시뮬레이션은 차량 설계에서 귀중한 도구가 되었다. 그러나 그들은 ...
더 보기
Figure 6. Scour depth (in negative value) at different views around pier

Three-dimensional numerical simulation of local scour around circular bridge pier using Flow-3D software

Flow-3D 소프트웨어를 이용한 원형 교각 주변 지역 scour의 3 차원 수치 시뮬레이션 To cite this article: Halah Kais Jalal and ...
더 보기
Abb. 3 Detail des Rechens am Vorversuch zum Seilrechen – Blick in Fließrichtung

Implementation of an angled trash rack in the 3D-numerical simulation with FLOW-3D

Abstract Sebastian Krzyzagorski · Roman Gabl · Jakob Seibl · Heidi Böttcher · Markus AuflegerOnline publiziert: 17. Februar 2016© Die ...
더 보기
Figure 1: Die configuration for a multi-attribute composite die for high die life and self-lubricating surface

Innovative Die Material and Lubrication Strategies for Clean and Energy Conserving Forging Technologies

청정 및 에너지 절약 단조 기술을 위한 혁신적인 다이 재료 및 윤활 전략 이 최종 기술 보고서에는 수상 번호 DE-FC07-01ID14206에 ...
더 보기
(a) Moving Reference Frame

Study on Swirl and Cross Flow of 3D-Printed Rotational Mixing Vane in 2×3 Subchannel

A thesis/dissertationsubmitted to the Graduate School of UNISTin partial fulfillment of therequirements for the degree ofMaster of ScienceHaneol Park07/09/2019Approved by_________________________AdvisorIn ...
더 보기
Figure 8 Evaluation test of thermal sprayed coatings

Development of Advanced Materials and Manufacturing Technologies for High-efficiency Gas Turbines

고효율 가스 터빈용 신소재 및 제조 기술 개발 Mitsubishi Heavy Industries Technical Review Vol. 52 No. 4 (December 2015) 가스 ...
더 보기
Figure 2. Diagram. Schematic design of a living snow fence. Source: Wyatt et al., 2012b

Design of Living Barriers to Reduce the Impacts of Snowdrifts on Illinois Freeways

눈사태가 일리노이 고속도로에 미치는 영향을 줄이기 위한 생활장벽 설계 John Petrie, et al. (2020) 일리노이 교통 센터 시리즈 번호 20-019, 연구 ...
더 보기
aerospace-sloshing-simulation

Aerospace Sloshing Dynamics

Sloshing Dynamics 우주선의 연료 탱크에서 추진체의 움직임에 대한 지식은 작동 및 성능의 다양한 측면을 이해하는 데 필수적입니다. 추진체 운동은 액체 ...
더 보기
Numerical modelling of a two-degree-of-freedom Wave Energy Converter

Energy Presentations | 에너지 프레젠테이션

Energy Presentations | 에너지 프레젠테이션

 지난 사용자 컨퍼런스에서 에너지 산업을위한 FLOW-3D의 응용에 초점을 맞춘 사용자 프레젠테이션을 다운로드하십시오  .

2019 년

Numerical modelling of a two-degree-of-freedom Wave Energy Converter: Creation, validation and utilization of the model

2 자유도 파동 에너지 변환기의 수치 모델링 : 모델의 생성, 검증 및 활용

Eliseo Marchesi, Politecnico di Milano / Studio Frosio srl
Marco Negri 및 Stefano Malavasi, Politecnico di Milano
Filippo Palo, XC Engineering Srl

Numerical modelling of a two-degree-of-freedom Wave Energy Converter
Numerical modelling of a two-degree-of-freedom Wave Energy Converter

이 연구의 목적은 FLOW-3D를 통한 수치 모델링입니다., Politecnico di Milano에서 실험실 규모 테스트를 거친 특정 Wave Energy Converter (WEC) : 두 개의 자유도 진동 체 시스템 인 EDS (Energy Double System), 급증하는 패들. 두 본체는 서로 연결되어 있으며 각각은 지상에 반응하는 PTO (Power Take-Off)에 연결됩니다. 수치 모델은주기적인 파동으로 사용 가능한 실험 테스트에 대해 검증되었습니다. 첫 번째 시뮬레이션은 실험실 테스트의주기적인 파동을 재현하는 것을 목표로했습니다. 그런 다음 실험 시스템에 해당하는 EDS의 수치 모델을 생성하고 이전에 모델링 한주기 파를 적용했습니다. 수치 방법의 품질이 확인되면 EDS 시스템의 새로운 구성에 대한 시뮬레이션이 수행되었습니다. 첫 번째 시뮬레이션 시리즈에서는이 매개 변수가 실험 모델에서 최적화되지 않았기 때문에이를 최적화하기 위해 패들 PTO의 댐핑이 변경되었습니다. 그 후, EDS 동작은 이전에 시뮬레이션 된 주기적 파동과 에너지 적으로 동일한 임의 파동에서 조사되었습니다.

다운로드

2015 년

Numerical simulation of extreme wave loading on an axisymmetric point absorber wave energy converter in a survival sea state

생존 해상 상태에서 축 대칭 점 흡수 파 에너지 변환기에 대한 극한 파 하중의 수치 시뮬레이션

Peter Arnold, Minerva Dynamics Limited

생존 해상 상태에서 파력 에너지 변환기 (WEC)가 경험할 수있는 힘과 모멘트의 초기 평가는 개념 설계 단계에서 특히 중요합니다. 현재까지 WEC의 생존 가능성을 평가하는 데 사용되는 주요 방법은 모델 규모 탱크 테스트이지만 일반적으로 10m에서 15m 범위의 상당한 파도 높이를 갖는 생존 파도의 크기로 인해 탱크 테스트 프로그램은 소규모를 사용해야합니다. 관련 계측 및 물리적 확장 문제가있는 스케일 프로토 타입 또는 관련 비용이 더 큰 대형 프로토 타입에 의존합니다. 보다 최근에 CFD는 더 작은 계산 비용으로 인해 불규칙한 파동 스펙트럼과는 반대로, 단독 집중 또는“New Waves”를 사용하여 정적 및 부동 구조에 대한 비선형 파동 부하를 평가하는 데 사용되었습니다. 그러나 이러한 초기 연구는 WEC 설계 엔지니어가 결과 부하의 통계적 분포를 필요로하기 때문에 WEC 생존에 필요한 조건 만 제공하지만 충분하지는 않습니다. 이 연구의 목적은 해법을 얻기 위해 스펙트럼 성분 수와 메시 미세 조정 수준 및 파동 탱크 폭이 감소 된 불규칙 파 스펙트럼을 활용하여 CFD 모델에서 모델링 가정에 첨부 된 중요성을 재분배하는 것입니다. 합리적인 시간 척도로. 그 결과 WEC 운동 및 하중은 주요 통계 매개 변수 측면에서 35 스케일 축 대칭 점 흡수기의 생존에 대한 고 충실도 탱크 테스트 결과와 비교됩니다. 이 연구의 목적은 해법을 얻기 위해 스펙트럼 성분 수와 메시 미세 조정 수준 및 파동 탱크 폭이 감소 된 불규칙 파 스펙트럼을 활용하여 CFD 모델에서 모델링 가정에 첨부 된 중요성을 재분배하는 것입니다. 합리적인 시간 척도로. 그 결과 WEC 운동 및 하중은 주요 통계 매개 변수 측면에서 35 스케일 축 대칭 점 흡수기의 생존에 대한 고 충실도 탱크 테스트 결과와 비교됩니다. 이 연구의 목적은 해법을 얻기 위해 스펙트럼 성분 수와 메시 미세 조정 수준 및 파동 탱크 폭이 감소 된 불규칙 파 스펙트럼을 활용하여 CFD 모델에서 모델링 가정에 첨부 된 중요성을 재분배하는 것입니다. 합리적인 시간 척도로. 그 결과 WEC 운동 및 하중은 주요 통계 매개 변수 측면에서 35 스케일 축 대칭 점 흡수기의 생존에 대한 고 충실도 탱크 테스트 결과와 비교됩니다.

다운로드

Wave propagation and reflection at an inclined plane – simulations and experiments

경사면에서의 파동 전파 및 반사 – 시뮬레이션 및 실험

Boris Huber, 비엔나 기술 대학교

20m 길이의 수로에서 물리적 모델 테스트를 수행하여 수로 끝의 경사면에서 파동 전파 및 반사를 관찰했습니다. 웨이브 생성은 익스텐더 휠에 의해 앞뒤로 움직이는 바닥에 장착 된 패들로 이루어졌습니다. 파도의 전파는 수위 측정에 의해 여러 지점에서 기록되었습니다. 실험은 다양한 파동으로 실행 된 다음 FLOW-3D 로 시뮬레이션되었습니다  . 또한 CFD 시뮬레이션에서 적절한 경계 조건을 얻기 위해 돌과 천공 시트로 구성된 파동 흡수 경계를 사용한 실험을 수행하고 다른 경계 조건에서 시뮬레이션을 실행했습니다.

다운로드

2013 년

Flap type wave power device in near shore conditions

해안 근처에서 플랩 형 파력 장치

Ibis Group, Inc의 Stephen Saunders

FLOW-3D  v10.1을 솔루션 코드로 사용하여 단일 이동 플랩 파 전력 캡처 장치의 CFD 분석을 수행했습니다  . 이 작업의 목적은 프로토 타입을 제작하고 배치하기 전에 의도 한 작업 환경에서 플랩이 겪는 힘을 예측하는 것입니다. FLOW-3D  는 이동하는 공기 / 물 인터페이스의 역학을 캡처하는 데 필요한 강력한 VOF 모델 때문에 경쟁 업체보다이 프로젝트에 선택되었습니다. 또한  FLOW-3D움직이는 고체 물체를 표현하는 FAVOR ™ 방법은 움직이는 플랩을 시뮬레이션하는 데 중요합니다. 플랩 형상의 성능을 시뮬레이션하고 다가오는 파도가 플랩 표면에 수직 인 해안 근처 조건에서 평가되었습니다. 테스트 된 모델은 플랩 끝 주변의 흐름 특성을 평가하기 위해 3D로 구성됩니다. 현재까지 두 개의 바다 상태가 테스트되었습니다. 이들은 플랩에 도달하는 즉시 깨지는 깨지지 않는 팽창 및 파도입니다. 예상대로 깨지지 않는 팽창은 중립의 양쪽에서 동일한 플랩 편향과 거의 대칭 인 부드러운 플랩 동작을 유도합니다. 파동 사례를 깨는 결과는 비대칭 동작과 하중으로 훨씬 더 극적입니다.

다운로드

Ocean waves resonance analysis of an oscillating water column energy converter

진동 수주 에너지 변환기의 해양 파도 공명 분석

José Manuel Grases ; 센데 키아

SDK Marine은 진동하는 수주 챔버 내의 물에 잠긴 수력 터빈을 기반으로 파도에서 전기 에너지를 수집하는 새로운 방법을 개발하고 있습니다. FLOW-3D  는 챔버 내부와 외부의 흐름 동작을 이해하는 데 사용되었습니다. 이 프로젝트의 주요 목표는 내부 수위를 측정하고이를 외부 파 여기와 비교하여 장치의 응답을 얻는 것이 었습니다. 또한 장치의 수력을 계산하기 위해 서로 다른 수력 터빈의 거동을 시뮬레이션하기 위해 서로 다른 다공성 멤브레인을 구현했습니다.

다운로드

collapsed-raised-fluid-column-figure-1-1

Steady-State Accelerator for Free-Surface Flows

자유 표면 흐름을 위한 정상 상태 가속기

이 기사에서 Tony Hirt 박사는 다가오는 FLOW-3D  v12.0 릴리스에서 사용할 수있는 새로운 Steady-State Accelerator에 대해 설명합니다  .

일시적인 흐름의 점근 적 상태를 계산하는 것보다 안정된 자유 표면 흐름을 생성하는 더 빠른 방법이 자주 필요합니다. 상황은 압축성 흐름 솔버를 사용하여 비압축성 흐름을 해결하는 것과 유사합니다. 후자의 경우 압축 파는 붕괴하는 데 오랜 시간이 걸리고 결과적으로 비압축성 흐름을 남길 수 있습니다. 이에 따라 자유 표면 흐름에서 유체는 비압축성이지만 표면 파동은 안정된 자유 표면 구성을 생성하는 데 오랜 시간이 걸릴 수 있습니다.

비압축성 흐름의 경우, 압축 파를 심각하게 감쇠시키는 반복적 인 프로세스 (즉, 압력-속도 반복)를 사용합니다. 물리적으로 반복은 압력과 같은 파동이 국부적 인 영역에 영향을 미치는 짧은 거리를 이동하도록 허용하지만 압력 장에 상당한 노이즈를 유발할 수있는 장거리 전파 및 반사를 피할 수있을만큼 빠르게 감쇠됩니다.

이 노트에서 자유 표면 셀에 적용된 간단한 압력 조정은 표면 교란에 대한 감쇠력으로 작용합니다. 이 댐핑은 안정적인 자유 표면 구성에 대한 접근을 가속화합니다.

Steady-State Accelerator Idea

유체 인터페이스 또는 자유 표면은  VOF (Volume-of-Fluid) 기술을 사용하여 FLOW-3D 에서 추적됩니다 . 유체 변수 F의 비율은 유체가 차지하는 영역을 찾습니다. 유체에 고정 된 자유 표면이있는 경우 유체를 정의하는 F 값도 안정된 값을 유지해야합니다. F가 일정하려면 표면에 수직 인 유체 속도가 0이어야합니다. 물론 표면에서의 접선 유체 속도는 0 일 필요는 없습니다. 예를 들어, 위어 위의 흐름에는 일정한 흐름이 있지만 계단에서 나오는 흐름의 위치와 모양은 변하지 않습니다.

자유 표면 흐름에 대한 정상 상태 솔버를 사용하려면 흐름의 비압축성을 유지하면서 정상 표면 속도를 0으로 유도하는 방법을 찾아야합니다.

이를 수행하는 한 가지 방법은 정상 속도를 0으로 유도하는 방식으로 표면 압력을 조정하는 것입니다. 특히 정상 속도에 비례하는 총 표면 압력에 “댐핑”압력 기여를 추가하는 것입니다. 속도는 표면 밖으로 향하고 그렇지 않으면 음수입니다.

정상 속도가 0에 가까워지면 수정 압력도 0이되어야 표면이 고정 위치를 초과하지 않게됩니다. 물론 보정이 너무 크면 오버 슈트가 발생할 수 있습니다. 따라서 안정적인 보정 적용을 위해서는 몇 가지 제한 요소가 있어야합니다.

계수 약어 ssacc 을 나타내며, S는 teady- S 테이트 액세서리 elerator이 새로운 옵션을 활성화하는 프로그램 입력에 추가되었다. ssacc 의 값 은 편리한 상한 인 1.0보다 작거나 같아야합니다. 프로그램 내에서 댐핑 압력에 자동으로 적용되는 여러 제한 기가 불안정 해 지거나 일시적인 현상에 악영향을 미치는 것을 방지합니다.

안정성 및 댐핑 리미터에 대한 이전 문제는 강조되어야합니다. 정상 상태 가속기를 사용하면 자유 표면 흐름의 모든 과도 현상이 더 이상 완전히 사실적인 것으로 볼 수 없습니다. 댐핑 압력은 물리적 인 힘이 아니라 파동 전파와 반사를 줄이는 메커니즘입니다. 댐퍼는 큰 과도 현상의 발생을 방해하지 않도록 고안되었으며 흐름이 안정됨에 따라 안정된 결과를보다 빠르게 얻는 데에만 기여해야합니다. 그러나 사용자는 리미터가 예상하지 못한 초과 댐핑에 대해 주의를 기울여야 합니다. 이는 댐핑 계수 ssacc 의 입력 값을 줄임으로써 제거 할 수 있습니다 .

두 가지 예는 정상 상태 가속기의 댐핑 메커니즘이 어떻게 작동하는지 보여줍니다.

Steady-State Accelerator Examples

Collapse of Raised Fluid Column

첫 번째 예는 길이 100cm, 깊이 5cm의 2 차원 물 웅덩이로 구성됩니다. 물을 담은 탱크의 모든 경계는 대칭 경계입니다. 수영장 중앙에는 폭 10cm, 높이 3cm의 수영장 위에 물 블록이 있습니다. 이 블록은 중력으로 인해 물에 떨어지고 충돌 지점에서 멀리 이동 한 다음 탱크 끝에서 반사되는 파도를 생성합니다. 100 초 후에도 반복되는 반사 때문에 여전히 상당한 파동 작용이 있습니다 (그림 1).

새로운 정상 상태 가속기를 계수 ssacc = 1.0 과 함께 사용하면 모든 파동이 빠르게 감쇠되어 거의 평평한 표면이됩니다. 일부 잔류 흐름은 표면 아래에 남아 있지만 점도의 작용으로 서서히 감쇠됩니다 (그림 2). 이 예에서 추가 된 댐핑은 특히 인상적입니다.

Figure 1. Column collapse without damping. Times of flow plots are 0.0, 10.0, and 100.0s. Bottom figure is the mean kinetic energy vs. time.
Figure 2. Column collapse with damping coefficient ssacc=1.0 at times of 0.0, 10.0 and 100.0s. Bottom figure is the mean kinetic energy vs. time.

 

사각형 격자에서 45 °의 정사각형 채널에서 모세관 상승

수직 채널에서 유체의 모세관 상승은 간단한 분석할 수 있으며 솔루션이 있는 양호한 정상 상태 문제입니다. 중력에 대해 상승 된 유체의 양은 벽의 접착력, 즉 접촉각의 코사인에 표면 장력 곱하기 접촉 선 길이에 의해 결정됩니다. 이 예에서 유체는 물이며 표면 장력은 70 dynes / cm이고 접촉각은 30 °입니다. 채널은 단면이 정사각형이며 가장자리 길이가 0.707cm이고 직사각형 격자에서 45 ° 회전합니다. 문제가 x 및 y 방향으로 대칭을 이루기 때문에 그리드의 사분면 만 모델링됩니다. 그리드의 바닥에는 제로 게이지 압력의 물이 있으며 그리드의 가장자리 길이는 0.0125cm (41x41x80 셀)입니다. 상승시켜야하는 이론적 유체 량은 0.04373cc입니다. 그림 3a는 정상 상태 결과를 보여줍니다. 이는 감쇠 사용 여부와 비슷합니다. 댐핑없이 계산된 유체의 양은 이론 값보다 1.74 % 높습니다. 그림 3b와 같이 댐핑이 있는 경우에는 2.24 %가 너무 높습니다. 가속기를 사용하면 정상 상태는 약 0.15 초에 도달하는 반면 표준 솔버는 0.8 초 후에 만 ​​정상 상태 솔루션을 생성하므로 5 배 이상 더 오래 걸립니다.

Figure 3a. Capillary rise in square channel without damping pressures.
Figure 3b. Histories of fluid volume in the two simulations (blue is with damping).

ssacc가 1.0보다 작으면 댐핑이 적어 수렴에 더 빨리 도달합니다. 1.0을 포함한 모든 ssacc 값은 댐핑되지 않은 ssacc = 0.0 경우와 비교하여 이론과 밀접하게 일치하고 후면 벽에 적은 양의 유체를 나타내는 수렴된 솔루션을 만듭니다.

뒤쪽 벽에있는 작은 유체 조각은 평형 위치를 초과하는 유체의 오버 슈트에서 발생하며, 이는 점성력으로 인해 정착하는 데 오랜 시간이 필요한 소량의 유체를 벽에 남기고 뒤로 떨어집니다. 이 오버 슈트는 ssacc 가 0이 아닐 때 제거됩니다 .

Figure 1. Alaska requires minimum water depth for fish passage to be 2.5 times the height of the caudal fin (D) (Hotchkiss and Frei 2007).

EFECTS OF HYDRAULIC STRUCTURES ON FISH PASSAGE: AN EVALUATION OF 2D VS 3D HYDRAULIC ANALYSIS METHODS

물고기 통로 유압 구조효과 :2D VS 3D 유압 분석 방법의 평가

ABSTRACT

채널 스패닝 유압 구조물은 상류 물고기 이동에 대한 장벽 역할을 할 수 있습니다. 이러한 종단 적 서식지 연결의 중단과 관련된 부정적인 결과는 정확하고 실행 가능한 평가 기술의 필요성을 강조합니다.

3 차원 평가 방법은 인스트림 구조에서 복잡한 흐름을 해결하고 물고기 움직임을 정확하게 예측하는 것으로 나타났습니다. 그러나 3 차원 모델링은 시간과 리소스 요구 사항으로 인해 비실용적 일 수 있습니다.

이 연구는 2 차원 전산 유체 역학 모델과 통계 분석을 사용하여 콜로라도 주 리옹에있는 화이트 워터 공원 구조의 수력 조건을 설명하는 것을 조사합니다. 물고기의 움직임 관찰은 잠재적 인 수영 경로를 나타내는 공간적으로 명시적이고 연속적인 경로를 따라 결과 수력 변수와 쌍을 이룹니다.

로지스틱 회귀 분석은 흐름 깊이와 속도가 어류 통과와 밀접한 관련이 있음을 나타냅니다. 결합 된 깊이 및 속도 변수무지개 송어 (92 %를 정확하게 예측Oncorhynchus mykiss) 및 갈색 송어 (Salmo trutta)는 이 유압 구조에서) 움직임 관찰을 합니다.

이 연구의 결과는 2 차원 분석 방법이 3 차원 분석이 불가능한 경우 유사한 수력 학적 구조가 어류 통과에 미치는 영향을 평가하는 비용 효율적인 접근 방식을 제공할 수 있음을 시사합니다. 또한,이 연구의 결론은 비교적 낮은 수영 성능을 가진 송어와 물고기 모두에 대한 관리 및 설계 결정을 안내하는 데 사용할 수 있습니다.

서문

수력 구조물은 수생 생물의 종 방향 서식지 연결을 의도적으로든 우연히든 효과적으로 차단할 수 있습니다. 의도적 장벽은 일반적으로 침입성 종의 도입 또는 교잡을 방지하기 위해 관리자에 의해 배치됩니다 (Holthe et al. 2005; Fausch et al. 2006). 그러나 구조물을 설계하고 설치할 때 물고기 통행 촉진을 고려하지 않았기 때문에 장벽이 더 자주 생성됩니다. 따라서 인위적 장애로 인해 전 세계 수로가 분열되었습니다 (Williams et al. 2012). 철새 어종의 성공적인 수명주기를 위해서는 종단 서식지 연결이 필수적입니다 (Schlosser and Angermeier 1995). 상류 이동에 대한 지연 또는 종료는 인구에 부정적인 영향을 미치고 생태계 기능을 방해 할 수 있습니다 (Beechie et al. 2010). 

수로를 가로 지르는 수력 구조물은 어류 통행에 미치는 영향을 철저히 평가하지 않고 하천과 강에 계속 배치됩니다 (Cada 1998; Noonan et al. 2012). 그러나 강 조각화와 관련된 문제에 대한 인식이 높아짐에 따라 설계 프로세스 전반에 걸쳐 물고기 통과 문제가 해결되는 방식에서 패러다임 전환이 일어나고 있습니다 (Katopodis and Williams 2012). 비 연어 종은 경제적 가치가 높은 종을 선호하는 경우가 많지만, 칼륨 종의 상류 이동 요구가 점점 더 중요하게 고려되고 있습니다 (Santos et al. 2012; Silva et al 2012) (Katopodis 2005; Roscoe and Hinch 2010). . 천연 자원 관리자는 제안 된 수력 구조물에 대해 의견을 제시하고 허용하도록 자주 요청받으며 (Kondratieff 2015),이 검토 과정에서 엔지니어와 과학자는 설계에 대한 예상 어류 통과 성능에 대한 모델 기반 증거를 제공하도록 요청받을 수 있습니다. 어류 통행과 관련하여 기존의 수력 구조물을 평가하고 우선 순위를 정하는 여러 방법이 현재 사용 가능하지만 (Kemp et al. 2010), 이전에 이 중요한 지점에서 제안된 구조물의 통행 효율성을 평가할 수있는 정확하고 실행 가능한 승인 및 설치 도구가 필요합니다.  

이러한 요구를 해결하는 데 초점을 맞춘 이전 작업은 3D 수력 모델링 기술이 상류 어류 이동을 평가할 목적으로 채널 스패닝 구조의 복잡한 유체 역학을 적절하게 해결할 수 있음을 보여주었습니다 (Stephens 2014).

이러한 새로운 3D 분석 방법은 전체 예측 정확도가 80 % 이상 (Stephens 2014)으로 매우 효과적 일 수 있지만 3D CFD (전산 유체 역학) 모델을 개발하는 데는 시간과 리소스가 많이 사용됩니다.

추가 데이터 수집, 소프트웨어 라이선스, 모델링 전문 지식 등에 대한 필요성은 많은 하천 관리 결정에 3D 분석을 비실용적으로 만들 수 있습니다. 다양한 2D 모델 플랫폼이 홍수 배출을 추정하고 (Horritt and Bates 2002; Merwade et al. 2008) 인스 트림 평가에 광범위하게 사용 되었기 때문에 실무 엔지니어와 과학자는 대부분의 수력 구조물 프로젝트에서 2D 수력 모델링을 수행 할 가능성이 더 높습니다.

물고기 서식지 (Clark et al. 2008; Katopodis 2012). 2D 및 3D 유압 모델의 실제 비교가보고되었지만 (Lane et al. 1999; Shen and Diplas 2008; Kolden 2013), 어류 통과에 대한 2D 및 3D 모델 기반 평가의 효능을 조사한 연구는 현재에서 발견되지 않았습니다.

목표

천연 자원 관리자와 설계 엔지니어가 Stephens (2014)의 매우 효과적인 3D 방법에 더 쉽게 접근 할 수 있도록하기 위해이 연구는 자유롭게 사용할 수있는 산업 표준 2D CFD 모델 인 River2D (Steffler and Blackburn 2002)를 사용하여 타당성을 조사합니다. 수력 구조가 어류 통로에 미치는 영향을 평가합니다.

유사한 접근 방식을 기반으로하고 이전의 수력 학 및 어류 이동 데이터 세트 (Fox 2013, Kolden 2013, Stephens 2014)를 사용하여 이 2 개의 연구는 2D 분석 방법을 사용하여 St. Vrain River의 WWP (화이트 워터 파크) 구조를 평가합니다. Lyons, CO.이 연구의 구체적인 목표는 다음과 같습니다. 

1. WWP 구조에서 복잡한 유압 환경을 설명하는 2D CFD 모델을 개발합니다. 

2.이 2D CFD 모델의 결과를 사용하여 WWP 구조를 통해 잠재적 인 물고기 이동 경로를 따라 연속적이고 공간적으로 명시적인 수력 학적 설명을 생성합니다. 

3. 무지개 송어 (대해 사용 가능한 어류 이동 데이터와 가장 밀접하게 관련된 수리적 변수를 결정Oncorhynchus mykiss) 및 갈색 송어 (Salmo trutta)에합니다. 

4. 이전에 개발 된 3D 접근 방식 (Stephens 2014)의 PIT (Passive Integrated Transponder) 태그 연구의 움직임 데이터를 기반으로 한 예측 평가 능력을이 연구의 2D 접근 방식과 비교합니다. 

5. 어류 통행의 관점에서 수력 구조물에 대한 비용 효율적인 평가를 통해 천연 자원 관리자 및 설계자를 지원하기위한 권장 사항을 제공합니다.

배경

상류 어류 이동에 대한 장벽은 유속 깊이, 유속 또는 유속과 거리의 조합을 포함한 다양한 물리적 조건에 의해 생성 될 수 있습니다 (Coffman 2005; Cahoon et al. 2005). 깊이 장벽은 일반적으로 흐름 깊이가 너무 얕아 통과 시도를 허용하지 않을 때 생성됩니다.

깊이 장벽은 또한 자리 잡은 구조물의 낙하 높이 및 플런지 풀 깊이가 도약 제약으로 인해 통과를 허용하지 않을 때 존재할 수 있습니다. 유속이 구조물을 통과하려는 물고기의 수영 능력을 초과 할 때 속도 장벽이 생성되어 상류 진행을 방해합니다. 수력 구조물에 의해 생성 된 난류는 물고기의 통과에도 역할을 할 수 있습니다. 조건에 따라 난류는 물고기 수영에 긍정적 인 영향과 부정적인 영향을 모두 미칠 수 있습니다 (Liao 2007; Cotel and Webb 2012; Lacey et al. 2012).  

수영 성능 지표는 종종 기존의 수력 학적 구조가 물고기 통행의 장벽으로 작용하는지 여부를 평가하는 데 사용됩니다. 이러한 메트릭 중 가장 일반적인 것은 달리기 속도라고도 하는 버스트 수영 속도와 지구력 곡선입니다 (Castro-Santos et al. 2013).

물고기는 지속, 연장, 파열의 세 가지 수영 모드를 나타냅니다 (Peake et al. 1997). 지속적인 수영은 이론적으로 무한정 유지 될 수 있지만 장시간 및 버스트 수영 속도는 제한된 시간 동안만 유지 될 수 있습니다.

지구력 곡선은 세 가지 수영 모드 (Videler and Wardle 1991)에 걸쳐 연속적으로 수영 속도와 피로 시간 사이의 역 관계를 설명하여 생성됩니다. 버스트 수영 속도는 속도 장벽을 식별 할 때 유용하며 (Haro et al. 2004) 지구력 곡선은 잠재적 인 완전 장벽을 식별하는 데 도움이됩니다 (Castro-Santos et al. 2013). 현재 물고기 수영 성능과 난류 임계 값 또는 분포 사이의 물리적 관계는 잘 알려져 있지 않습니다 (Liao 2007).

그러나 총 운동 에너지 (TKE), 총 수력 변형, 레이놀즈 전단 응력 및 와도와 같은 일부 프록시 변수는 난류가 어류에 미치는 영향을 정량화 할 때 유용한 것으로 나타났습니다 (Nestler et al. 2008; Cotel and Webb 2012; Lacey et al. 2012; Silva et al. 2012). 

장벽은 완전 할 수 있으며, 물고기 통행을 허용하지 않거나 선택적 통행 성공이 생리적 또는 수리적 특성에 따라 결정되는 경우 부분적 일 수 있습니다. 이 연구의 목적을 위해 총 시도 횟수에 대한 성공적인 통과 횟수를 기반으로 한 인구 수준의 통과 효율을 사용하여 유압 구조로 인한 상류 이동 억제 정도를 정량화합니다 (Haro et al. 2004). 다양한 방법 개발되었습니다. 

장벽이 물고기 통로 (켐프와 O’Hanley 2010)에 영향을 미치는 방법을 정량화하기 위해  한 가지 접근 방식은 통계 모델을 사용하여 통과 효율 추정치를 0 ~ 100 %의 연속 척도로 표현할 수 있습니다. 과거에는 규칙 기반 또는 회귀 기법을 사용하여 암거 (Coffman 2005; Burford et al. 2009), 도로 횡단 (Warren and Pardew 1998) 또는 수로 실험 설정 (Haro et al. 2004)을 다양한 성공으로 평가했습니다.

통계적 방법은 다양한 척도에서 수리적 변수에 대한 정보를 결합하여 통과에 큰 영향을 미치는 변수를 식별 할 수 있습니다 (Kemp and O’Hanley 2010). 이러한 모델은 현장 기반 어류 이동 관찰을 사용하여 검증 할 수도 있습니다 (Coffman 2005; Burford 2009).

2014 년에 Stephens는 3D CFD 모델 출력 (Kolden 2013)을 활용하여 수력 구조물에서 물고기 통과를 평가하기위한 연속적이고 공간적으로 명시적인 분석 방법을 만드는 새로운 통계 방법을 개발했습니다. 이 방법은 콜로라도에있는 3 개의 파도 생성, 인공 화이트 워터 파크 (WWP) 구조물에서 수집 한 수력 측정 및 PIT 태그 통과 관찰 (Fox 2013)을 통해 검증되었습니다. 통계 결과에 따르면 Stephens (2014) 방법은 전체 정확도가 80 % 이상인 통과 효율을 예측할 수 있습니다. 

Stephens는 3D CFD 모델의 결과를 사용했지만 다른 연구에서는 2D CFD 모델을 사용하여 물고기와 관련된 규모의 복잡한 흐름을 설명하는 데 초점을 맞추 었습니다 (Lane et al. 1999; Crowder and Diplas 2000; Shen and Diplas 2008). 2D CFD 모델링의 주요 관심사는 물고기 서식지 및 수영 성능에 중요한 중간 규모 기능과 관련된 복잡성을 포착 할 수 있는지 여부였습니다 (Crowder and Diplas 2000).

혼합된 결과는 서식지 평가를 위해 모델링되는 도달 범위의 특성에 따라 2D CFD 모델이 수력 조건에 대한 적절한 설명을 제공하거나 제공하지 않을 수 있음을 보여줍니다 (Clark et al. 2008; Shen and Diplas 2008; Kozarek et al. 2010) . 서식지 또는 지형 모델링에 중점을 두는 경우 깊이 평균 2D 모델과 직접 비교할 때 3D 모델 사용이 선호되었습니다 (Lane et al. 1999; Shen and Diplas 2008). 그러나 수력 구조물에서 상류 어류의 움직임을 평가할 때 2D 및 3D 모델의 성능을 비교 평가 한 연구는 거의 없습니다. 

이 연구에서 CFD 모델의 비교는 2D 소프트웨어 River2D와 3D 소프트웨어 FLOW-3D에 중점을 둡니다 (Flow Science, 2009). 2D 모델과 3D 모델의 가장 큰 차이점은 2D 모델은 각 계산 노드에서 유압 변수의 값을 깊이 평균한다는 것입니다. 이 깊이 평균은 구조물의 물고기 친화성에 큰 영향을 미칠 수있는 중요한 흐름 특징과 경계층 효과를 배제 할 수있는 잠재력을 가지고 있습니다.

예를 들어, 수심 평균 속도 값은 WWP 구조 하류의 수력 조건이 동일한 도달 범위 내의 자연 풀에있는 것과 유사하다고 잘못 제안 할 수 있습니다. 실제로 두 유동장은 어류 개체군에 다르게 영향을 미칠 수있는 고유 한 특성을 가지고 있습니다 (Kolden 2013). River2D는 또한 정수압과 일정한 수평 속도 분포를 가정하는 반면 FLOW-3D는 이러한 가정을 피할 수 있습니다.

대부분의 2D CFD 모델링 프로그램 (Toombes and Chanson 2011)에서 요구하는 정수압 가정은 가파른 경사 (> 10 %)와 급변하는 경사 (Steffler and Blackburn 2002)에서 계산 정확도를 제한합니다. 속도 분포가 일정하다는 가정은 수직 속도 구성 요소가 무시할 수 있음을 의미하며 본질적으로 2D CFD 모델을 사용하여 2 차 흐름 및 강한 순환을 분석하는 기능을 제거합니다 (Steffler and Blackburn 2002; Toombes and Chanson 2011).

이러한 가정과 2D 물리적 표현의 단순화 된 특성을 고려할 때 2D CFD 모델이 물고기 통과 예측 평가를 위해 채널 스패닝 구조의 복잡한 유체 역학을 적절하게 해결할 수 있는지 여부는 불분명합니다.

Figure 1. Alaska requires minimum water depth for fish passage to be 2.5 times the height of the caudal fin (D) (Hotchkiss and Frei 2007).
Figure 1. Alaska requires minimum water depth for fish passage to be 2.5 times the height of the caudal fin (D) (Hotchkiss and Frei 2007).
Figure 2. Depth (m) and velocity magnitude (m/s) River2D contours for 0.42 cms.
Figure 2. Depth (m) and velocity magnitude (m/s) River2D contours for 0.42 cms.

구체적인 내용은 아래 원문을 참고하시기 바랍니다.

Sketch of a subaerial landslide-induced tsunami wave

NUMERICAL SIMULATION OF THREE-DIMENSIONAL TSUNAMI GENERATION BY SUBAERIAL LANDSLIDES

SUBAERIAL LANDSLIDES에 의한 3 차원 쓰나미 생성의 수치 시뮬레이션

A Thesis by GYEONG-BO KIM
Submitted to the Office of Graduate Studies of
Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE

초록

쓰나미는 해저 지진으로 인해 종종 발생하는 해안 지역에 영향을 미치는 가장 치명적인 자연 현상 중 하나입니다. 그럼에도 불구하고 밀폐된 분지, 즉 피요르드, 저수지 및 호수에서, 수중 또는 해저 산사태는 유사한 결과로 파괴적인 쓰나미를 일으킬 수 있습니다. 큰 수역에 충돌하는 수중 또는 해저 산사태가 쓰나미를 발생시킬 수 있지만, 해저 산사태는 대응하는 것보다 훨씬 더 위협적인 쓰나미 발생원입니다.

이 연구에서 우리는 지하 산사태에 의한 쓰나미 발생에 대한 실험실 규모의 실험을 수치 모델과 통합하는 것을 목표로 합니다. 이 작업은 2 개의 3 차원 Navier-Stokes (3D-NS) 모델, FLOW-3D 및 당사가 개발 한 모델 TSUNAMI3D의 수치 검증에 중점을 둡니다.

이 모델은 Georgia Institute of Technology의 Hermann Fritz 박사가 이끄는 쓰나미 연구팀이 수행 한 이전의 대규모 실험실 실험을 기반으로 검증되었습니다. 일련의 실험실 실험에서 세 가지 대규모 산사태 시나리오, 즉 피요르드 유사, 곶 및 원거리 해안선이 선택되었습니다. 이러한 시나리오는 복잡한 파도 장이 지하 산사태에 의해 생성 될 수 있음을 보여주었습니다.

파동 장의 정확한 정의와 진화는 뒤 따르는 쓰나미와 해안 지역에서의 영향을 정확하게 모델링하는 데 중요합니다. 이 연구에서는 수치 결과와 실험실 실험을 비교합니다. 토양 유변학에 대한 방법론과 주요 매개 변수는 모델 검증을 위해 정의됩니다. 모델의 결과는 쓰나미 수치 모델의 검증을 위해 National Tsunami Hazard Mitigation Program (NTHMP), National Oceanic and Atmospheric Administration (NOAA) 지침에 명시된 허용 오차 미만일 것으로 예상됩니다.

이 연구의 궁극적 인 목표는 멕시코만과 카리브해 지역의 침수지도를 구축하는 데 필요한 해저 산사태 쓰나미에 대한 3D 모델의 실제 적용을 위한 더 나은 쓰나미 계산 도구를 얻는 것입니다.

주요 분석 이미지

 Sketch of a subaerial landslide-induced tsunami wave
Figure 1.4: Sketch of a subaerial landslide-induced tsunami wave: (a) cross section
defining parameters in the direction of slide motion; (b) plan view defining coordinate
system to reference and quantify the generated tsunami wave
A typical computational domain with moving and stationary objects
Figure 2.1: A typical computational domain with moving and stationary objects. Courtesy Dr. Juan J. Horrillo, Texas A&M at Galveston.
A typical tsunami computational domain
Figure 2.2: A typical tsunami computational domain: (a) Location of variables in a computational cell. The horizontal (ui,j ) and vertical (vi,j ) velocity components are located at the right cell face and top cell faces, respectively. The pressure pi,j and VOF function Fi,j are located at the cell center; (b) Volume and side cell apertures. Courtesy Dr. Juan J. Horrillo, Texas A&M at Galveston.
Figure A.1: Configurations of boundary conditions for fjord case: FLOW-3D
Figure A.1: Configurations of boundary conditions for fjord case: FLOW-3D

<자료 안내>

원문 다운로드

Water & Environmental 논문 자료보기

Wave Energy Devices

파동 에너지 장치 모델링
최근 몇 년 동안 파력 에너지와 같은 재생 가능 자원을 사용하여 환경 영향이 적은 에너지를 생산하는 신기술 개발에 대한 국제적인 관심이 기하 급수적으로 증가했습니다. 바다 (해류, 파도 등)에서 전기를 유도하는 파동 에너지 장치는 특히 중요하며 FLOW-3D로 정확하게 모델링 할 수 있습니다.

포인트 흡수 장치
점 흡수 장치는 수면의 파도를 사용하여 에너지를 생성하는 많은 파도 장치 중 하나입니다. 포인트 흡수 에너지 장치는 기본적으로 파도에서 에너지를 흡수하고 바닥에 대한 부력 상단의 움직임을 전력으로 변환하는 부동 구조입니다.

이 시뮬레이션은 부력 구형 구조가 위에 있는 포인트 흡수기 장치를 보여주고, 들어오는 파동의 파고와 수조에 따라 위아래로 움직입니다. FLOW-3D의 이동 객체 모델은 x 또는 y 방향으로 이동을 제한하면서 z 방향으로 커플링 모션을 허용하는 데 사용됩니다. 스톡스 유형의 파장은 진폭 5m, 파장은 100m로 사용되었습니다. RNG 모델은 파동이 포인트 업소버 장치와 상호작용할 때 발생하는 난류를 포착하기 위해 사용되었습니다. 예상대로, 많은 난류 운동 에너지가 장치 근처에서 생성됩니다. 그림은 난류로 인해 장치 근처의 복잡한 속도장이 진화하기 때문에 질량 중심의 불규칙한 순환 운동을 보여줍니다.

Multi-Flap, Bottom-Hinged Wave Energy Converter

Oscillating flap은 바다의 파동으로부터 에너지를 추출하여 기계 에너지로 변환합니다. 암은 Water wave에 반응하여 피벗 조인트에 장착된 진자로 진동합니다. 플랩을 배열로 구성하여 멀티플랩파 에너지 변환기를 만들 수 있습니다. 3개의 플랩 배열이 아래 왼쪽에 표시된 CFD 시뮬레이션에서 시뮬레이션됩니다. 모든 플랩은 하단에 힌지로 연결되며 폭 15m x 높이 10m x 두께 2m입니다. 어레이는 깊이 30m에서 주파수가 10초인 4m 진폭 파형으로 작동 중입니다. 시뮬레이션은 한 플랩이 배열 내의 다른 플랩에 미치는 영향을 연구하는 데 중요한 중심 평면을 따라 복잡한 속도 ISO 표면을 보여줍니다. 3개의 플랩이 유사한 동적 모션으로 시작하는 동안, 곧 플랩의 상호 작용 효과가 모션을 위상 밖으로 렌더링합니다. 우측에는 유사한 플랩 에너지 변환기가 표시되어 있습니다. 이 시뮬레이션에서 플랩은 가장 낮은 지점에서 완전히 물에 잠깁니다. 이러한 에너지 변환기를 표면 천공 플랩 에너지 변환기라고 합니다. 이 두 시뮬레이션 예는 모두 미네르바 다이내믹스에 의해 제공되었습니다.

Oscillating Water Column

진동하는 물 기둥은 부분적으로 잠긴 속이 빈 구조입니다. 그것은 물의 기둥 위에 공기 기둥을 둘러싸고 수선 아래의 바다로 열려 있습니다. 파도는 물 기둥을 상승 및 하강시키고, 차례로 공기 기둥을 압축 및 감압합니다. 이 갇힌 공기는 일반적으로 기류의 방향에 관계없이 회전 할 수있는 터빈을 통해 대기로 흐르게됩니다. 터빈의 회전은 전기를 생성하는 데 사용됩니다.

위의 CFD 시뮬레이션은 진동하는 water columns를 보여줍니다. FLOW-3D로 포착된 물리학을 강조하기 위해 물기둥이 중공 구조에서 상승 및 하강하는 부분만 모델링합니다. 시뮬레이션은 파형 생성의 다른 선택을 제외하고 유사한 결과를 전달합니다. 왼쪽의 시뮬레이션은 웨이브 유형 경계 조건을 사용하고 오른쪽의 시뮬레이션은 움직이는 물체 모델을 사용하여 실험실에서 수행한 것처럼 차례로 웨이브를 생성하는 움직이는 플런저를 생성합니다. 각 시뮬레이션에 대해 속이 빈 구조의 압력 플롯이 표시됩니다. 결국 그 압력에 기초하여 터빈이 회전 운동으로 설정되기 때문에 챔버에서 얼마나 많은 압력이 생성되는지 아는 것이 중요합니다.

Sloshing & Slamming

Sloshing & Slamming

LNG 선과 같은 전 세계 항해 선박의 내부 슬로싱 운동에 의해 유발되는 가속 및 하중은 이러한 선박의 안전 설계에 중요한 요소입니다. 선박은 이러한 파도 슬로싱 역학으로 인해 상당한 내부 강제력을 경험할 수 있으며,화물 탱크는 격리 시스템 및 지지 구조물에 추가 하중을 경험할 수 있습니다. FLOW-3D의 비관성 기준 프레임 모델을 사용하면 컨테이너에서 유체 모션을 정확하게 추적하기 위해 복잡한 모션 매개 변수를 규정 할 수 있습니다. FLOW-3D는 연료 탱크의 액체화물과 추진체 움직임을 정확하게 예측합니다.

슬래밍 분석과 “그린 워터”(파동 스프레이 또는 선체 오버 토핑의 결과로 데크를 덮는 물)는 FLOW-3D의 일반적인 응용 분야이며, 이러한 문제에 대한 FLOW-3D의 TruVOF 접근 방식은 엔지니어에게 고유 한 기능을 제공합니다.

이 시뮬레이션은 FLOW-3D로 캡처한 것처럼 파도로 움직이는 LNG 유조선의 복잡하고 완전히 결합된 동작을 보여줍니다. 탱크에서 유체의 슬로싱은 선박의 역학에 영향을 미치며, 여기에는 bilge keels 과 bulbous bow가 포함되어 입사파에서 안정성을 유지하는데 도움이됩니다. 또한 탱크의 유체 슬로싱은 예상대로 선박의 움직임에서 위상 이동이 있습니다. FLOW-3D의 움직이는 물체 모델, 밀도 평가 모델 및 파동 생성 경계 조건은 해양 엔지니어가 유사한 문제를 분석하는 데 도움이 될 수 있습니다.

Coastal applications using FLOW-3D/연안(해변) FLOW-3D해석사례

배수(Backwaters)이론

  • Air entrainment(공기 혼입 모델)
  • Turbulence(난류 모델)
  • Waves(파동 모델)
  • Sediment scour and deposition(세굴 모델)

하구 매커니즘(Estuarine mechanisms)

  • Air entrainment(공기 혼입 모델)
  • Turbulence(난류 모델)
  • Density Evaluation(밀도 유동 모델)
  • Wind

Wave generation

  • Solitary Wave
  • Linear Wave
  • JONSWAP
  • Pierson-Moskowitz
  • Stokes Wave

Wave 생성하는 모델은 크게 위의 5가지 모델이 있습니다. 아래는 위의 5가지 모델에 대한 해석 사례를 보여줍니다. 이를 참고하시면 해석에 도움이 됩니다.

부두에서 파도 부하 추정, 물리적 모델링 및 수치테스트, 새로운 Wave에 대한 고유 2차원 비선형 접근방식 등의 FLOW-3D결과는 실제 실험 데이터와 잘 일치함을 보여줍니다.

  • Eillott, T., and Fullarton, M., “Cyclone wave loads on wharf structure using the new wave approach”, FLOW-3D Americas User conference, 2014

세굴 모델(Sediment scour and deposition)

Sediment scour and deposition
  • Critical Shields number definition(임계 Shields 수) : 0.05
  • Bed Load Transport Rate equation : Meyer-Peter & Muller equation
  • Richardson-Zaki coefficient multipller : 1
  • 다음과 같이 Wave와 세굴(Sediment)를 같이 고려해서 해석하는 것을 추천합니다.
    – 퇴적물 탱크의 파동(Solitary wave)
    – 무연탄 및 모래
    – 움직이는 물체 모델을 사용하여 생성된 파도

해석 결과

Wave Forces on Coastal Bridges

Wave Forces on Coastal Bridges

This article was contributed by Jun Jin, Assistant Professor in the Maritime Systems Engineering Department at Texas A&M University at Galveston.

 

최근 몇 년 동안 허리케인은 멕시코만 연안의 4개 주에서 여러 해안 고속도로 교량에 구조적 피해를 입혔습니다. 허리케인 피해를 방지하기 위해 연안 교량에 가해지는 파도 하중의 크기는 합리적인 정확성을 가지고 결정되어야 합니다. 일반적으로 구조에 대한 파도 하중의 계산은 구조 상호작용의 여러 특성으로 인해 복잡해집니다. 이 본문은 이러한 복잡한 문제들 중 해결할 수 있는 몇 가지를 보여줍니다.

Validation of Stokes Nonlinear Wave Solution in FLOW-3D

FLOW-3D는 해안 교량의 파도 하중을 계산하고 구조물 고도 및 녹색갑판 위의 물 하중에 대한 영향을 조사하기 위해 사용되었습니다.  FLOW-3D에서 StokesWave의 정확도를 평가하기 위해 계산 도메인의 세개 지점을 선택했습니다(그림 1). 시뮬레이션을 하는 동안, 각 입자 소스는 초당 10개의 입자를 방출했습니다. 입자 궤적은 그림 1에 표시되어 있습니다. 한 지점의 입자 속도를 이론적 결과와 비교했습니다(그림 2a 및 2b). FLOW-3D에 의해 계산된 입자의 역학은 Stokes비선형 파동이론과 잘 일치한다고 결론내렸습니다.

Figure 1. Particle trajectories and pressure variation during a wave simulation

Figure 2a. Comparison of computed particle velocities of a point with theoretical results

Figure 2b. Comparison of computed particle velocities of a point with theoretical results

 

Wave Loads vs. Superstructure Elevations

다리 갑판을 다른 높이로 설정하고 2m파형 높이와 6초 주기의 파형으로 시뮬레이션을 수행했습니다. 그림 3은 파형력이 최대치에 있을 때의 파동 단면과 압력 등고선을 보여 줍니다.

Figure 3. Wave profile and pressure contours at different moments (Pressure unit: Pa).

Flow Field under Bridge Deck

시뮬레이션은 다른 파도 높이로 실행되었습니다. 그림 4에서 볼 수 있듯이 교량 바닥판 밑의 유동장은 복잡합니다. 대들보와 상부 구조물에 의해 형성된 개방된 공간내에서, 물의 흐름은 경계를 따라 원형패턴으로 있었습니다. 또한 유동의 복잡성은 웨이브 높이가 증가함에 따라 증가하는 것으로 관찰되었다.

Observations

시뮬레이션 결과에서 다음과 같은 결론을 도출할 수 있습니다.

  1. 상부 구조물의 높이가 증가함에 따라 교량 상부 구조물의 수직 방향 파동력이 감소하는 반면, 수평 방향 파형력은 교량의 상부 구조물의 상부 표면이 웨이브 높이의 아래에 있을 때 큰 영향을 미치지 않습니다.
  2. 교량 상층부의 물 하중은 파고와 함께 증가합니다. 파손이 거의 발생하는 경우 최대 물 하중이 최대 수직 파력의 50 %를 초과했습니다. 그러나 그린 워터 하중과 상향 파력 사이의 위상차로 인해 수직 파력의 최대 값은 물의 부하에 의해 현저히 감소하지는 않았습니다.
  3. 교량 상부 구조물 주위의 유동장은 평행한 대들보에 의해 차단됩니다. 그러므로 웨이브 힘에 대한 순전히 기반한 방법은 웨이브 힘을 계산하는데 사용할 수 없습니다.

Figure 4. Vectors of water particle velocity at the moment of maximum wave force (Velocity unit: m/s)

Learn more about the power and versatility of modeling coastal and maritime applications with FLOW-3D>

 

 

Landslide-Induced Wave Hazard

Landslide-Induced Wave Hazard 

Figure 1. The outskirts of Chungtangh village

인도 Sikkim에 위치한 The Teesta III Hydropower Project는 가파르고 좁은 히말라야 계곡에 위치한 60m의 Concrete Face Rockfill Dam (CFRD)이 포함되어 있습니다. 이 계곡은 지진 활동이 활발하며 가파른 경사면은 산사태를 발생시킬 수 있습니다. 댐 상류 저수지의 산사태로 CFRD를 범람할 수 있다는 우려가 있었습니다. 몇 초 이상 과도하게 지속되면 오버플로우로 인해 CFRD가 잘못될 수 있습니다. 비록 댐이 무너지지 않았지만, 여전히 Chungtangh에 있는 상류쪽 작은 마을은 홍수가 날 것이라는 우려가 있었습니다.

Teesta강 계곡의 가장 가파른 경사면은 댐의 바로 상류에 위치해 있는데, 댐의 산사태가 가장 일어날 가능성이 높은 지역입니다. 이 분석의 목적은 저수지에 대한 산사태를 시뮬레이션하고 그 결과로 발생하는 파도가 댐에 넘치는지 여부를 결정하는 것이었습니다.

Moving Objects Model Used to Simulate Landslide                                      

Tecsult는 저수지의 침전물과 퇴적물을 모델링하는데 성공적이었기 때문에 FLOW-3D를 선택하여 이를 시뮬레이션하였습니다. 저수지의 시뮬레이션은 시작점으로 사용되었습니다. FLOW-3D의 Moving Objects모델은 산사태를 시뮬레이션하는데 사용되었으며 VOF모델은 웨이브 생성을 시뮬레이션하는 데 사용되었습니다.

저수지의 산사태를 추정하기 위해서는 여러가지 방법이 고려되었습니다. 경험적 방법은 흔히 산사태가 발생한 파도를 평가하는데 사용되지만, 이러한 방법은 여러가지 면에서 부족합니다. 이러한 방법은 근접 필드 또는 스플래시 영역에 대한 정보를 제공하지 않습니다. 댐은 슬라이드 면과 매우 가깝기 때문에 스플래시 영역을 아는 것이 중요했습니다. CFRD는 몇 초 이상 overflow를 견딜 수 없었습니다. FLOW-3D는 미끄러운 지형 질량과 물 사이의 완전 결합된 상호 작용을 계산하여 시나리오를 3 차원에서 시뮬레이션하는 방법을 제공합니다.

이 문제를 시뮬레이션하기 위해 간단하고 작은 크기의 자유 낙하 블록으로 구성된 실험과 비교하였습니다. 이 경우는 아래 동영상에 나와 있습니다. 그 결과로 생긴 파도 높이는 그 실험과 잘 맞았습니다.

이 모델의 STL파일은 FLOW-3D로 직접 가져옵니다. 예상 산사태 지역의 크기는 지질 정보와 주변 산사태 관측치를 바탕으로 결정되었습니다. 30,000m³, 100m높이의 산사태가 310만 셀의 메쉬로 시뮬레이션 되었습니다. 높이가 1m인 측면 3m의 균일한 셀을 사용했습니다. 최대 슬라이딩 속도는 진입 지점에서 23m/s에 도달했습니다. 파도는 높이 8m, 속도 10m/s로 댐에 도달하여 몇 초 동안 범람했습니다. 그 결과로 상류 마을에서는 홍수가 나타나지 않았습니다.

Figure 3. Prediction of wave height in the splash zone and near field in a small reservoir, with refraction.

Figure 4. Wave heights plotted against each other

Figure 5. Downstream view of TEEST III dam and water intake CATIA model

Conclusions

이 작업의 주된 관심사는 댐의 범람으로 인해 댐과 Chungtangh 마을이 파괴될 수 있었다는 것입니다. 그러나 시뮬레이션에 따르면 댐은 잠시 동안만 범람했고 파도는 마을에 닿지 않았습니다. Chungtangh마을은 강 위에 충분히 높기 때문에, 그것을 범람시키기 위해서는 상당한 파도의 높이가 필요할 것입니다.

 

Aerial Landslide Generated Wave Simulations

Aerial Landslide Generated Wave Simulations

Aerial Landslide Generated Wave(ALGW)는 수역에 영향을 미치는 빠른 슬라이드의 결과이다. 이것은 암석에 의해 생성된 작은 파도 이거나, 3000만 입방 미터의 암석으로 인한 500m를 초과하는 파도 일 수도 있다.
공학적 관점에서 보면 ALGW는 근접한 해안을 따라 인간이 거주하는 인구/자산이 있는 수역에서 발생할 때 큰 관심을 가진다. 여기서 파동이 발생하면 해안선이 파손되고 홍수가 날수 있으며, 댐붕괴로 인한 사망까지 일으킬 수 있다(Müller-Salzburg, 1987). 결과적으로, ALGW에 의해 야기되는 최대 파도 상승을 예측하는 것은 경제적, 환경적, 안전상의 이유로 매우 중요합니다.
안타깝게도 분석적인 솔루션이 없는 매우 복잡한 문제로, 유체 역학적인 측면에서뿐만 아니라 지질학적인 관점(즉, 크기/기하학적인 슬라이드의 밀도 프로파일)에서도 마찬가지입니다. 이와 같이, 대부분의 현장 별 ALGW 최대 파형 예측은 확장된 물리적 모델을 사용하여 평가되었다. 일부는 전산유체역학(CFD) 소프트웨어를 기반으로 할 수도 있지만 비용이 많이 들며, 특히 풀 스케일 3차원 문제의 경우 정확성에 대한 논쟁의 대상이 되고 있습니다.
그러나 컴퓨터 하드웨어와 CFD소프트웨어가 계속 발전함에 따라 이제 CFD를 사용하여 ALGW를 실제로 시뮬레이션할 수 있게 되었습니다. 이와 같이 본 연구는 고 충실도의 물리적 모델 데이터를 FLOW-3D와 비교하여 ALGW를 CFD시뮬레이션을 검증하기 위한 지속적인 노력으로 진척시키는 것을 목표로 한다.
다음 절에서는 실제 및 수치 모델 설정에 대한 개요를 제공한다. 뿐만 아니라, 생성된 데이터와 간단한 비교를 제공한다.

Experimental Setup
물리적 실험은 Northwest Hydraulic Consultants 노스 밴쿠버, 캐나다 실험실에서 만들어졌고 실험을 거쳤다. 그것은 30° 경사의 서쪽 벽을 가진 0.5미터 폭의 수로, 45°의 경사진 동쪽 벽, 그리고 두개의 북쪽과 남쪽 측면에 수직 벽, 그리고 1.025m의 수평 단면을 가진 0.610m 너비의 수로로 구성되었다. ALGW를 생성하고 평가하기 위해, 45° 경사 노즈를 가진 0.177×0.305×0.305m의 아크릴 박스를 사용한 6초 시험을 사용했다.
이 슬라이드를 놓았을 때, 슬라이드는 (중력에 의해) 0.607m 심층수에 충돌하기 전에 서쪽 경사면에서 0.768m 아래로 이동했다. 그 후, 물을 통해 또 다른 1.05m를 이동하여 정지 블록을 치기 시작했다. 슬라이드 가속 및 변위뿐만 아니라 파고 높이는 6 초 실험 전체에 대해 100Hz의 주파수에서 기록되었다. 이 데이터를 수집하는 데 사용 된 도구는 다음과 같다.

  • 컴퓨터화된 데이터 수집 시스템
  • 슬라이드의 시간에 따라 이동 한 거리를 측정하는 문자열 가변 저항기
  • 슬라이드 가속도를 측정하는 1 차원 가속도계
  • 물의 주요 본체 내에 배치 된 3 개의 1 차원 커패시턴스 웨이브 – 프로브
  • 웨이브 런업을 캡처하기 위해 동쪽 경사면을 따라 사용되는 저항 사다리꼴 웨이브 프로브
  • 타이밍 스위치 캡처 슬라이드 릴리스 시간 사용
  • 흑백 비디오 카메라

테스트가 반복 가능하고 오작동이 발생하지 않았는지 확인하기 위해 테스트를 5 번 반복하고 각 장비에 대해 평균을 구했다.

Numerical Model Setup
물리적 실험의 전산화 된 3 차원 모델을 제작한 STL 파일을 FLOW-3D로 가져왔다. 일단 FLOW-3D에 들어간 3D 모델은 약 1,370 만개의 0.0075m 크기의 정사각형 셀로 이산화되었고, 벽을 둘러싸고있는 6 개의면 각각에 ‘wall’경계가 사용되었다.
슬라이드를 일반적인 이동 물체로 설정하고, 물리 모델로부터 수집 된 데이터(즉, 가속 및 변위 데이터의 후 처리)에 기초하여 속도가 주어졌다. 동서면 경사면의 표면 거칠기는 0.00025m으로 설정되었다. 모델링 된 유체는 293k의 물이었고, 동적 RNG 난류 모델이 기본 설정과 함께 사용되었다(implicit pressure solve; and, explicit viscous stress, free surface pressure, advection, moving object/fluid coupling solvers).
물리적 모델과 마찬가지로 FLOW-3D는 6 초의 시간을 시뮬레이트하지만 실제 모델과 같이 매 0.01 초가 아닌 0.02 초마다 데이터를 저장하였다(데이터 관리 관점에서 선택하였음).

Result

FLOW-3D 실험의 결과는 그림에 나와 있다. 4개의 웨이브 각각에 대해 실험 시간 동안 파고를 보여준다. 이와 같이, 제시된 파도 높이는 단순히 flume을 통해 전파되는 파도의 구현(즉, 2 차원의 경우에서 볼 수있는 것)이 아니라 오히려 여러 파도의 상호 작용으로 인한 파도 높이를 초래한다.

  • 슬라이드 충격시 발생하는 충격파(1차 신호)
  • 슬라이드 뒤의 충격파 충돌(2차 신호)
  • 북쪽, 동쪽, 서쪽 및 남쪽 벽에서의 웨이브 반사(3차 신호)

또한 길이 방향의 FLOW-3D 데이터(중심선에서)를 실제 모델 비디오 위에 겹쳐서 자유 표면의 FLOW-3D 글로벌 예측을 평가했다. 이것은 아래의 동영상에서 볼 수 있다.
그림과 위의 비디오를 보면 FLOW-3D 데이터가 웨이브 프로브 1, 2 및 3의 경우 물리적 데이터를 매우 잘 일치한다는 것을 알 수 있다. 하지만 웨이브 프로브 4에 대해서는 정확도가 떨어진다.
FLOW-3D 시간 데이터와 관련된 오류는 각 웨이브 프로브에 대한 RMSE (root-mean-square-error)를 취하여 평가된다.

Discussion
이 조사에서 실제 모델의 고 충실도 데이터는 ALGW로 최대 파도 상승에 대한 FLOW-3D 예측과 비교되었다. RNG 모형의 기본 설정을 사용하여 FLOW-3D는 주요 수역 내에서 파고를 정확하게 재현 할 수 있었다. 그러나 최대 파동은 약 43%가 넘었다.
최대 웨이브 런업을 줄이기 위해 몇 가지 대안인 FLOW-3D 물리 설정이 사용되었다. 그러나 43 % 이하로 떨어지는 것은 불가능했다. 이러한 대체 시뮬레이션에 대한 주목할만한 관찰은 다음과 같다.

  • first-order momentum advection scheme의 0.01m 메쉬는 최대 파동 상승 오차가 96% 인 반면 동일하게 0.0075m 메쉬의 오차는 130%였다. 그러나 second-order로 변경하면 0.01 m 및 0.0075 m 메시의 경우 각각 55% 및 43%의 오차가 발생한다. 또한 메쉬 셀 크기를 0.005m으로 줄이면 80%의 오차가 발생한다.
  • 이 테스트 케이스에서 가장 중요한 매개 변수는 momentum advection scheme이다. 평균적으로 second-order를 사용하면 first-order대비 오차가 약 50% 감소한다.
  • FLOW-3D의 MP 버전을 사용하여 0.005m의 메쉬 셀 크기를 사용해야 한다. 해석 시 CPU 시간은 33 시간이었다. 비교를 위해 FLOW-3D의 SMP 버전은 0.0075m의 메쉬 셀 크기로 시뮬레이션을 실행하는 데 26시간이 필요했지만 MP 버전은 4.5시간 밖에 걸리지 않았다.

[1] 3.5GHz 8 코어 AMD FX-8320 프로세서에서 약 6초의 시뮬레이션 시간이 대략 26시간 소요되었다.

References
Fritz, H. M., Hager, W. H., & Minor, H.-E. (2004). Near Field Characteristics of Landslide Generated Impulse Waves. Journal of Waterway, Port, Coastal & Ocean Engineering, 130(6), 287–302. doi:10.1061/(ASCE)0733-950X(2004)130:6(287)
Miller, D. J. (1960). Giant Waves in Lituya Bay Alaska (Geological Survey Professional Paper No. 354-C). Washington, D.C.: United States Government Printing Office.
Müller-Salzburg, L. (1987). The Vajont catastrophe— A personal review. Engineering Geology, 24(1–4), 423–444. doi:10.1016/0013-7952(87)90078-0

Wave Forces on Coastal Bridges

Wave Forces on Coastal Bridges

This article was contributed by Jun Jin, Assistant Professor in the Maritime Systems Engineering Department at Texas A&M University at Galveston.

 

최근 몇 년 동안 허리케인은 멕시코만 연안의 4 개 주에서 여러 해안 고속도로 교량에 구조적 피해를 입혔습니다. 허리케인 피해를 방지하기 위해 연안 교량에 가해지는 파도 하중의 크기는 합리적인 정확성을 가지고 결정되어야 합니다. 일반적으로 구조에 대한 파도 하중의 계산은 구조 상호작용의 여러 특성으로 인해 복잡해집니다. 이 본문은 이러한 복잡한 문제들 중 해결할 수 있는 몇 가지를 보여줍니다.

Validation of Stokes Nonlinear Wave Solution in FLOW-3D

FLOW-3D는 해안 교량의 파도 하중을 계산하고 구조물 고도 및 녹색갑판 위의 물 하중에 대한 영향을 조사하기 위해 사용되었습니다.  FLOW-3D에서 StokesWave의 정확도를 평가하기 위해 계산 도메인의 세개 지점을 선택했습니다(그림 1). 시뮬레이션을 하는 동안, 각 입자 소스는 초당 10개의 입자를 방출했습니다. 입자 궤적은 그림 1에 표시되어 있습니다. 한 지점의 입자 속도를 이론적 결과와 비교했습니다(그림 2a 및 2b). FLOW-3D에 의해 계산된 입자의 역학은 Stokes비선형 파동이론과 잘 일치한다고 결론내렸습니다.

Figure 1. Particle trajectories and pressure variation during a wave simulation

Figure 2a. Comparison of computed particle velocities of a point with theoretical results

Figure 2b. Comparison of computed particle velocities of a point with theoretical results

 

Wave Loads vs. Superstructure Elevations

다리 갑판을 다른 높이로 설정하고 2m파형 높이와 6초 주기의 파형으로 시뮬레이션을 수행했습니다. 그림 3은 파형력이 최대치에 있을 때의 파동 단면과 압력 등고선을 보여 줍니다.

Figure 3. Wave profile and pressure contours at different moments (Pressure unit: Pa).

Flow Field under Bridge Deck

시뮬레이션은 다른 파도 높이로 실행되었습니다. 그림 4에서 볼 수 있듯이 교량 바닥판 밑의 유동장은 복잡합니다. 대들보와 상부 구조물에 의해 형성된 개방된 공간내에서, 물의 흐름은 경계를 따라 원형패턴으로 있었습니다. 또한 유동의 복잡성은 웨이브 높이가 증가함에 따라 증가하는 것으로 관찰되었다.

Observations

시뮬레이션 결과에서 다음과 같은 결론을 도출할 수 있습니다.

  1. 상부 구조물의 높이가 증가함에 따라 교량 상부 구조물의 수직 방향 파동력이 감소하는 반면, 수평 방향 파형력은 교량의 상부 구조물의 상부 표면이 웨이브 높이의 아래에 있을 때 큰 영향을 미치지 않습니다.
  2. 교량 상층부의 물 하중은 파고와 함께 증가합니다. 파손이 거의 발생하는 경우 최대 물 하중이 최대 수직 파력의 50 %를 초과했습니다. 그러나 그린 워터 하중과 상향 파력 사이의 위상차로 인해 수직 파력의 최대 값은 물의 부하에 의해 현저히 감소하지는 않았습니다.
  3. 교량 상부 구조물 주위의 유동장은 평행한 대들보에 의해 차단됩니다. 그러므로 웨이브 힘에 대한 순전히 기반한 방법은 웨이브 힘을 계산하는데 사용할 수 없습니다.

Figure 4. Vectors of water particle velocity at the moment of maximum wave force (Velocity unit: m/s)

Learn more about the power and versatility of modeling coastal and maritime applications with FLOW-3D>

 

 

Interaction Between Waves and Breakwaters

Interaction Between Waves and Breakwaters

This article is an adapted version of an article  published in the journal of the Engineering Association for Offshore and Marine in Italy by Fabio Dentale, E. Pugliese Carratelli, S.D. Russo, and Stefano Mascetti. The first three authors are users at the University of Salerno; Mr. Mascetti is an engineer at XC Engineering, Flow Science’s associate for Italy and France.

 

방파제의 설계는 복잡한 자연 시스템 (바다와 해안)과 인공 구조물 (방파제)의 상호 작용에 대한 완전한 이해가 필요합니다. 일반적으로 설계 작업은 광범위한 물리적 모델링을 수반하므로 비용이 많이 들고 시간이 오래 걸릴 수 있습니다. 최근까지 방파제의 복잡한 측면은 상세한 수치 시뮬레이션에 너무 어려웠습니다. 이것은 물이 비정상적인 동작으로 복잡한 경로를 통해 흐르는 콘크리트 또는 암석 블록으로 구성된 방파제의 경우 특히 그렇습니다.

컴퓨팅 기술의 진보로 수치, 물리적 조사 간의 격차가 좁혀졌습니다. 상호 작용하는 개별 블록으로 구성된 견고한 구조를 정확하게 표현할 수 있으므로 블록 사이의 빈 공간 내에 수치적으로 유동 영역을 생성 할 수 있습니다. 이것은 방류수가 균일한 다공성 매질로 근사되는 Classical Darcy 주제에 고려될 수 없는 대류항 및 난류의 영향을 포함한 전체 유체 역학적 거동의 영향을 평가할 수 있게 합니다

Modeling Rubble Mound Breakwaters

The following examples describe cases where rubble mound breakwaters are modelled on the basis of their real geometry, taking into account the hydrodynamic interactions with the wave motion.

잔재물 분쇄기 모델링

다음의 사례는 잔해 분쇄물이 파도 운동과의 유체 역학적 상호 작용을 고려하여 실제 형상을 기반으로 모델링된 경우를 설명합니다.

Figure 1: Artificial blocks

Figure 2a: Submerged Breakwaters

Figures 2b and 2c: Emerged Breakwater – Accropode regular & Accropode irregular

 

방파제의 개략적인 표현을 고려하여 구체 세트로 재현한 것으로 the cube, the modified cube, the antifer, the tetrapod, the accropode, the accropode II, the coreloc, the xbloc,and the xbloc base 등과 같은 일반적으로 사용되는 인공 블록을 고려하기 위해 개발되었습니다. (그림 1).

방파제는 물에 잠기거나 잠긴 경우에는 문헌에 나와 있는 표준 실험식을 사용하여 크기를 결정하고 실제 기하학적 패턴을 따르는 전체 크기, 구조 및 물리적 모델링과 같이 수치적으로 구성했습니다 (그림. 2).

제안된 절차의 품질을 검증하기 위해 침수된 방파제에 대해 세 가지 기하학적 구조를 고려했다. 즉, 부유, 다공성, 고형물과 부유물(그림 2a)이 출현한 방파제의 경우, 두 가지 다른 기하학적 구조를 사용했다(Fig. 2b – 2c).

방파제가 결정되면 기하학적 구조을 FLOW-3D로 가져 와서 유체 역학적 작용을 평가 및 Wave propagation의 연구를 위해 테스트했습니다. 시뮬레이션은 RNG 난류 모델과 coarse격자 안쪽에 중첩된 미세한 격자가 있는 전산메쉬를 사용하여 Navier-Stokes 방정식을 3 차원으로 통합하여 수행되었습니다.

수중 장벽 (계산 영역: 90 × 1.9 × 6.5m)의 경우, 포함된 메쉬 블록은 동일한 크기 (0.30 × 0.27 × 0.30m)의 46,200 개의 요소로 구성되며 중첩된 블록은 2,353,412 개의 요소로 구성하였습니다(0.061 × 0.055 × 0.061m).

방파제에도 동일한 기준이 적용되었습니다. 포함된 격자 블록은 150,000개의 요소(0.50×0.20×0.30m)로 구성되며, 중첩된 블록은 2,025,000개의 요소(0.10× 0.10×0.10m)로 생성되었습니다.

Figures 3a and 3b: Mesh views of submerged breakwater (3a above) & emerged breakwater (3b below)

Figures 4b: Emerged Breakwater – Accropode regular

Figures 4a: Submerged breakwater

결과 중 일부는 다음 이미지에 요약되어 있습니다. 그림 4에서 3 차원 영역의 2 차원 단면을 따른 압력 및 난류 에너지가 나타납니다. 그림 5에는 서로 다른 순간에 잡힌 자유 표면의 3 차원 형상이 나타나있습니다.

유동경로를 따라 개별 솔리드 요소의 윤곽의 유체 역학에 의한 유동 변화는 쉽게 검출 가능합니다. 이것은 자유 표면의 3 차원 재구성에서 가장 잘 드러나며 (그림 5) 방파제에 대한 파동 작용의 효과가보다 자세하게 표현됩니다.

Figures 5a: Submerged breakwater.

Figures 5b: Emerged Breakwater – Accropode regular.

Figures 5c: Emerged Breakwater – Accropode irregular  

Conclusions

잠수함이나 해상 구조물 간의 상호 작용을 정확히 표현하기 위한 Navier-Stoke기반 수치 시뮬레이션을 활용한 방법, 그리고 유체 움직임이 입증되었습니다. 시뮬레이션은 난류 시뮬레이션을 위한 RANS와 자유 표면 계산을 포함하는 첨단 컴퓨터 유체 동적 소프트웨어 시스템(FLOW-3D)을 사용하여 수행되었습니다.

이 결과는 블록 사이의 경로 내에서 유체 운동의 상세한 그림을 제공함으로써 기존의 흐름 방법보다 더 정확한 시뮬레이션을 제공함을 보여줍니다. 블록을 사용하여 기존의 누설 흐름 방법보다 더 정확한 시뮬레이션을 제공합니다. 원칙적으로 모든 관련 부품(필터, 코어 및 토우)에서 구조물이 물에 잠기거나 나타나는 경우 시뮬레이션이 가능하며 제한은 없습니다.

Further studies will be aimed at assessing the stability of individual blocks through the use of the Moving Object model in FLOW-3D.

 

조선/해양 분야

Coastal & Maritime

FLOW-3D 는 선박설계, 슬로싱 동역학, 파도에 미치는 영향 및 환기를 포함하여 해안 및 해양 관련 분야에 사용할 수 있는 이상적인 소프트웨어입니다.

자유 표면 유체 역학, 파동 생성, 움직이는 물체, 계선 및 용접 공정과 관련한 FLOW-3D 의 기능은 해양 및 해양 산업에서 CFD 공정을 모델링하는 데 매우 적합한 도구입니다. 해안 응용 분야의 경우  FLOW-3D  해안 응용 분야의 경우 FLOW-3D  는 해안 구조물에 대한 심한 폭풍 및 쓰나미 파동의 세부 사항을 정확하게 예측하고 돌발 홍수 및 중요 구조물 홍수 및 피해 분석에 사용됩니다. 기능은 다음과 같습니다.

  • 자유 표면 – 파동 유체 역학 및 오버 토핑 : 규칙 및 불규칙파 및 파동 스펙트럼 (Pierson Moskowitz, JONSWAP)
  • Seakeeping – slamming, planing, porpoising 및 선체 선체 변위 : 완전히 결합된 선박 및 수중 차량 유체 역학
  • 선체 – Resistance, stability and dynamics: surging, heaving, pitching and rolling motion (response amplitude operators or RAOs)
  • 슬로싱 – LNG / 밸러스트 탱크
  • 해양 공학 – 파동 에너지 변환기
 

해안 응용 분야의 경우, FLOW-3D 는 강력한 폭풍과 쓰나미 현상에 의한 해안 구조물이 받는 영향에 대한 세부 사항 예측, 돌발 홍수에 의한 중요한 시설물에 대한 정확한 피해 분석 등을 위해 사용됩니다.

Mooring Lines, Springs and Ropes

FLOW-3D (계류선 및 스프링 등)의 특수 물체를 다른 움직이는 물체에 부착하면 엔지니어가 선박 런칭, 부유 장애물 역학, 부표, 파도에너지 변환기 등을 정확하게 포착할 수 있습니다.

Welding

FLOW-3D 용접 모듈이 추가되면서 조선업계의 용접분야에서는 다공성 등 용접 결함을 최소화할 수 있어 선체의 품질을 크게 높이는 동시에 생산 시간을 최적화할 수 있습니다.

Coastal & Maritime Case Studies

FLOW-3D 사용자들은 연약한 해안선 보호, 구조물에 대한 파장 시뮬레이션, 선체 설계 최적화, 선박 내 환기 연구 등 해안 및 해양 애플리케이션에 FLOW-3D를 사용합니다.

우리는 보트가 세계 항해를 하면서 마주칠 것 같은 다양한 조건에서 항해를 할 수 있는지를 볼 수 있었습니다. 그리고 속도뿐만 아니라 연료 효율과 안전도 고려하도록 설계를 수정할 수 있었습니다.
– Pete Bethune, skipper of Earthrace

Lateral wave impact in waterWave resultsEarthrace vessel
Validation of Sloshing Simulations in Narrow Tanks / Aerial Landslide Generated Wave Simulations / Earthrace: Speed, Fuel Efficiency and Safety
Wave impact vertical displacementEmerged breakwater accropodeStokes theory horizontal velocity
Wave Impact on Offshore StructuresInteraction Between Waves and BreakwatersWave Forces on Coastal Bridges

기타

Bibliography

Models


관련 기술자료

Watershed area

Analysis on inundation characteristics by compound external forces in coastal areas

연안 지역의 복합 외력에 의한 침수 특성 분석 Taeuk KangaDongkyun SunbSangho Leec*강 태욱a선 동균b이 상호c*aResearch Professor, Disaster Prevention Research Institute, ...
더 보기
Fig. 3. Breakwaters model in Flow-3D with meshing geometry and boundary (a) circular slots (b) square slots.

Study of Unconventional Alternatives to Vertical Breakwater

수직 방파제에 대한 비전통적 대안 연구 Karim Badr Hussein and Mohamed IbrahimLecturer of Irrigation and Hydraulics, Faculty of Engineering, Al-Azhar ...
더 보기
Figure 3. Flow velocity on seawall in A2-3 modeling.

Modeling of the Changes in Flow Velocity on Seawalls under Different Conditions Using FLOW-3D Software

Open Journal of Marine ScienceVol.06 No.02(2016), Article ID:65874,6 pages10.4236/ojms.2016.62026 FLOW-3D 소프트웨어를 사용하여 다양한 조건에서 Seawalls의 흐름 속도 변경 모델링 Maryam Deilami-Tarifi1, ...
더 보기
Interaction between oblique waves and arc-shaped breakwater

Interaction between oblique waves and arc-shaped breakwater: Wave action on the breakwater and wave transformation behind it

XinyuHanaShengDongaYizhiWangbaCollege of Engineering, Ocean University of China, Qingdao, 266100, ChinabShandong Harbour Engineering Group Co., Ltd., Rizhao, 276826, China Highlights Interaction ...
더 보기
Figure 6. Maximum inundation field in simulations with (a) no barrier on the seawall (red line), (b) a 1 m barrier across the entire sea wall, and (c) a 1.7 m barrier partially installed on the seawall.

Storm surge inundation simulations comparing three-dimensional with two-dimensional models based on Typhoon Maemi over Masan Bay of South Korea

Jae-Seol Shim†, Jinah Kim†, Dong-Chul Kim‡, Kiyoung Heo†, Kideok Do†, Sun-Jung Park ‡† Coastal Disaster Research Center,Korea Institute of Ocean ...
더 보기
Fig. 3. Mesh and depth map for the storm surge model of ADCSWAN model.

ADCSWAN과 FLOW-3D 모델을 이용한 태풍 차바 내습 시 부산 마린시티의 침수범람 재현

최흥배․엄호식†․박종집․강태욱*, *** ㈜지오시스템리서치 선임, ** ㈜지오시스템리서치 책임, **** 부경대학교 박사 Reproduction of Flood Inundation in Marine City, Busan During the ...
더 보기
A photo of HeMOSU-1.

FLOW-3D를 이용한 해상 자켓구조물 주변의 세굴 수치모의 실험

Numerical Simulation Test of Scour around Offshore Jacket Structure using FLOW-3D J Korean Soc Coast Ocean Eng. 2015;27(6):373-381Publication date (electronic) ...
더 보기
Figure 6. Scour depth (in negative value) at different views around pier

Three-dimensional numerical simulation of local scour around circular bridge pier using Flow-3D software

Flow-3D 소프트웨어를 이용한 원형 교각 주변 지역 scour의 3 차원 수치 시뮬레이션 To cite this article: Halah Kais Jalal and ...
더 보기
圖1. 1 南海孤立內波空間分布圖(Hsu et al., 2000)

Numerical Modeling on Internal Solitary Wave propagation over an obstacle using Flow-3D

Keyword: Internal solitary waves, Numerical, Flow-3D, Computational Fluid Dynamics 연구자 : Yu-Ren Chen지도교수 : Dr John R C HsuJune 2012 ...
더 보기
유압 헤드 계산에서는 유선이 평행하다고 가정

FLOW-3D Output variables(출력 변수)

Output variables(출력 변수) FLOW-3D에서 주어진 시뮬레이션의 정확한 출력은 어떤 물리적 모델, 출력 위젯에 정의된 추가 출력 및 특정 구성 요소별 ...
더 보기

파도 / Waves

파도 / Waves

FLOW-3D 는 비정형 파뿐만 아니라 일반 선형 및 비선형파 표면을 시뮬레이션 할 수 있는 기능이 있습니다. 선형파는 작은 진폭 및 급경사를 갖는 사인파 표면 프로파일을 가지며, 비선형파는 선형 파보다 더 큰 진폭 (유한 진폭), 더 뾰족한 볏 및 평탄한 골짜기를 갖는다. 비선형 파는 파동 문자와 그 해를 구하기 위해 사용 된 수학적 방법에 따라 스톡 (stookes), 코니이드 (cnoidal) 파 및 독방 파로 분류 될 수 있습니다.

그림 1. 다른 진행파의 프로파일 비교
도 1 및도 2에 도시 된 바와 같이, 스톡스 파는 심층 및 과도수의 주기적인 파이다. Cnoidal 파는 얕은 물과 중간 물에서 긴주기적인 파이고 Stokes 파보다 더 뾰족한 볏과 평평한 골짜기를 가지고 있습니다. 스톡스와 코니 형 파와 달리 독방 파는 얕은 물과 과도 수에서 존재하는 비 주기적 파이다. 그것은 하나의 산마루와 물마루를 가지며 완전히 방해받지 않은 수면 위입니다. 수학적으로 파장이 무한대가 될 때 그것은 코니 형 파의 제한적인 경우입니다. 심층수, 과도 수 및 파도에 대한 얕은 물의 분류는 표 1에서 찾아 볼 수있다.

그림 2. 다양한 파도의 적용 범위 (Le Méhauté, 1976, Sorensen, 2005 및 USACE, 2008). d : 평균 수심; H : 파고; T : 파주기; g : 중력 가속도

선형 파 이론 (Airy, 1845)이 많은 응용 분야에서 사용되었지만 비선형 파 이론은 파동의 진폭이 작지 않은 경우 선형 파 이론보다 정확도가 크게 향상되었습니다. FLOW-3D 에서 3 개의 비선형 파 이론이 5 차 스톡스 파 이론 (Fenton, 1985), 스톡스 및 코니이드 파에 대한 푸리에 급수 방법 (Fenton, 1999), McCowan의 독방 파 이론 (McCowan, 1891, Munk, 1949). 그 중에서 Fenton의 Fourier 시리즈 방법은 선형 물, 스톡 (Stokes) 및 코니형 (cnoidal) 파를 포함하여 심층수, 과도 수 및 얕은 물에서 모든 종류의 주기적 전파 파들에 유효합니다. 또한 다른 웨이브 이론보다 정확도가 높습니다 (USACE, 2008). 따라서 모든 수심에서 선형 및 비선형 주기파의 모든 유형을 생성하는 것이 권장되는 방법입니다. solitary wave의 경우, FLOW-3D 에 사용 된 McCowan의 이론은 Boussinesq (1871)에 의해 개발 된 다른 널리 사용되는 이론보다 더 높은 주문 정확도를 갖는다.

그림 3. PM과 JOHNSWAP 스펙트럼 (USCE, 2006에서 적응)

Classifications d /\lambda
Deep water 1/2 to ∞
Transitional water 1/20 to 1/2
Shallow water 0 to 1/20

불규칙한 물결은 파도의 물성이 일정하지 않은 자연적인 바다의 상태를 나타냅니다. FLOW-3D에서 불규칙한 파동은 다양한 진폭과 주파수 및 임의의 위상 변이를 갖는 많은 선형 성분 파의 중첩으로 표현됩니다. Pierson-Moskowitz (Pierson and Moskowitz, 1964)와 JONSWAP 파력 에너지 스펙트럼 (Hasselmann, et al., 1973)은 FLOW-3D에서 구성 요소 파를 생성하기 위해 구현된다. 다른 웨이브 에너지 스펙트럼은 사용자 정의 데이터 파일을 가져와서 사용할 수 있습니다.

계산 시간을 절약하기 위해 웨이브는 메시 블록 경계에서뿐만 아니라 초기 조건으로 정의 될 수 있습니다.

아래의 애니메이션은 웨이브 초기화가 있거나없는 웨이브의 모든 유형에 대한 예제를 보여줍니다.
선형 및 비선형 수위 시뮬레이션을 위해 FLOW-3D 의 성공적인 적용이 이루어졌습니다. Bhinder 외의 예를 참조하십시오. al (2009), Chen (2012), Hsu et. al (2012) Thanyamanta et. al (2011) 및 Yilmaz et. 자세한 내용은 알 (2011)을 참조하십시오.






References

Airy, G. B., 1845, Tides and Waves, Encyc. Metrop. Article 102.

Bhinder, M. A., Mingham, C. G., Causon, D. M., Rahmati, M. T., Aggidis, G. A. and Chaplin, R.V., 2009, A Joint Numerical And Experimental Study Of a Surging Point Absorbing Wave Energy Converter (WRASPA), Proceedings of the ASME 28th International Conference on Ocean, Offshore and Arctic Engineering, OMAE2009-79392, Honolulu, Hawaii.

Boussinesq, J., 1871, Theorie de L’intumescence Liquide Appelee Onde Solitaire ou de Translation se Propageant dans un Canal Rectangulaire, Comptes Rendus Acad. Sci. Paris, Vol 72, pp. 755-759.

Chen, C. H., 2012, Study on the Application of FLOW-3D for Wave Energy Dissipation by a Porous Structure, Master’s Thesis: Department of Marine Environment and Engineering, National Sun Yat-sen University.

Fenton, J. D., 1985, A Fifth-Order Stokes Theory for Steady Waves, Journal of Waterway, Port, Coastal and Ocean Engineering, Vol. 111, No. 2.

Fenton, J. D., 1999, Numerical Methods for Nonlinear Waves, Advances in Coastal and Ocean Engineering, Vol. 5, ed. P.L.-F. Liu, pp. 241-324, World Scientific: Singapore, 1999.

Hasselmann, K., Barnet, T. P., Bouws, E., Carlson, H., Cartwright, D. E., Enke, K., Ewing, J. A., Gienapp, H., Hasselmann, D. E., Kruseman, P., Meerburg, A., Muller, P., Olbers, D. J., Richter, K., Sell, W., and Walden, H., 1973, Measurement of Wind-Wave Growth and Swell Decay During the Joint North Sea Wave Project (JONSWAP), German Hydrographic Institute, Amburg.

Hsu, T. W., Lai, J. W. and Lan, Y., J., 2012, Experimental and Numerical Studies on Wave Propagation over Coarse Grained Sloping Beach, Proceedings of the International Conference on Coastal Engineering, No 32 (2010), Shanghai, China.

Kamphuis, J. M., 2000, Introduction to Coastal Engineering and Management, World Scientific, Singapore.

Le Méhauté, B., 1976, An Introduction to Hydrodynamics and Water Waves, Springer-Verlag.

McCowan, J., 1891, On the solitary wave, Philosophical Magazine, Vol. 32, pp. 45-58.

Munk, W. H., 1949, The Solitary Wave Theory and Its Application to Surf Problems, Annals New York Acad. Sci., Vol 51, pp 376-423.

Pierson W. J. and Moskowitz, L., 1964, A proposed spectral form for fully developed wind seas based on the similarity theory of S.A. Kitiagordskii, J. Geophys. Res. 9, pp. 5181-5190.

Thanyamanta, W., Herrington, P. and Molyneux, D., 2011, Wave patterns, wave induced forces and moments for a gravity based structure predicted using CFD, Proceedings of the ASME 2011, 30th International Conference on Ocean, Offshore and Arctic Engineering, OMAE2011, Rotterdam, The Netherlands.

USACE (U.S. Army Corps of Engineers), 2006, Coastal Engineering Manual, EM 1110-2-1100, Washington, DC.

Yilmaz, N., Trapp, G. E., Gagan, S. M. and Emmerich, T., R., 2011 CFD Supported Examination of Buoy Design for Wave Energy Conversion, IGEC-VI-2011-173, pp. 537-541

난류 모델링

본 자료는 국내 사용자들의 편의를 위해 원문 번역을 해서 제공하기 때문에 일부 오역이 있을 수 있어서 원문과 함께 수록합니다. 자료를 이용하실 때 참고하시기 바랍니다.

Turbulence Modeling

The majority of flows in nature are turbulent. This raises the question, is it necessary to represent turbulence in computational models of flow processes? Unfortunately, there is no simple answer to this question, and the modeler must exercise some engineering judgment. The following remarks cover some things to consider when faced with this question.

난류 모델링

자연에서의 흐름은 대부분은 난류입니다. 이것은 유동의 수치해석 모델에서 난류를 표현할 필요가 있는가? 에 대한 의문이 생깁니다.  불행히도이 질문에 대한 답은 모델링을 할 경우 엔지니어가 공학적인 판단을 내려야합니다.  다음에 이 질문에 직면했을 때 고려해야 할  몇 가지를 설명합니다.

Definitions and Orders of Magnitude

The possibility that turbulence may occur is generally measured by the flow Reynolds number:

난류가 발생할 가능성은 일반적으로 흐름의 레이놀즈 수에 의해 측정됩니다.

where ρ is fluid density and μ is the dynamic viscosity of the fluid. The parameters L and U are a characteristic length and speed for the flow. Obviously, the choice of L and U are somewhat arbitrary, and there may not be single values that characterize all the important features of an entire flow field. The important point to remember is that Re is meant to measure the relative importance of fluid inertia to viscous forces. When viscous forces are negligible the Reynolds number is large.

여기서 ρ는 유체 밀도이고 μ는 유체의 동적 점도입니다. 매개 변수 L과 U는 흐름의 특성 길이와 속도입니다. 분명히 L과 U의 선택은 다소 임의적이며, 전체 유동장의 모든 중요한 특징을 특징 짓는 단일 값이 없을 수도 있습니다. 기억해야 할 중요한 점은 Re가 점성력에 대한 유체 관성의 상대적 중요성을 측정한다는 것입니다. 점성력을 무시할 수있는 경우 레이놀즈 수가 큽니다.

A good choice for L and U is usually one that characterizes the region showing the strongest shear flow, that is, where viscous forces would be expected to have the most influence.

L과 U에 대한 좋은 선택은 일반적으로 가장 강한 전단 흐름을 나타내는 영역, 즉 점성 힘이 가장 큰 영향을 미칠 것으로 예상되는 영역을 특징 짓는 것입니다.

Roughly speaking, a Reynolds number well above 1000 is probably turbulent, while a Reynolds number below 100 is not. The actual value of a critical Reynolds number that separates laminar and turbulent flow can vary widely depending on the nature of the surfaces bounding the flow and the magnitude of perturbations in the flow.

대략적으로 말하면, 1000을 훨씬 넘는 레이놀즈 수는 아마도 난류 일 수 있지만 100 미만의 레이놀즈 수는 그렇지 않습니다. 층류와 난류를 분리하는 임계 레이놀즈 수의 실제 값은 유동을 경계하는 표면의 특성과 유동의 섭동의 크기에 따라 크게 달라질 수 있습니다.

In a fully turbulent flow a range of scales exist for fluctuating velocities that are often characterized as collections of different eddy structures. If L is a characteristic macroscopic length scale and l is the diameter of the smallest turbulent eddies, defined as the scale on which viscous effects are dominant, then the ratio of these scales can be shown to be of order L/l≈Re3/4. This relation follows from the assumption that, in steady-state, the smallest eddies must dissipate turbulent energy by converting it into heat.

완전 난류 흐름에서는 다양한 와류 구조의 집합으로 특징 지어지는 변동 속도에 대해 다양한 스케일이 존재합니다. L이 거시적 길이 특성 척도이고, l을 점성 효과가 우세한 척도로 정의되는 가장 작은 난류 소용돌이의 직경인 경우, 이러한 척도의 비율은L/l≈Re3/4 정도인 것으로 표시 될 수 있습니다.  이 관계는 정상 상태에서 가장 작은 소용돌이가 난류 에너지를 열로 변환하여 발산해야한다는 가정에서 비롯됩니다.

Turbulence Models

From the above relation for the range of scales it is easy to see that even for a modest Reynolds number, say Re=104, the range spans three orders of magnitude, L/l=103. In this case, the number of control volumes needed to resolve all the eddies in a three-dimensional computation would be greater than 109. Numbers of this size are well beyond current computational capabilities. For this reason, considerable effort has been devoted to the construction of approximate models for turbulence.

난류 모델

스케일의 범위에 대한 위의 관계를 보면 적당한 레이놀즈 수 (예 : Re = 10 4 )에서도 범위가 세 자릿수인 L/l=103에 걸쳐 있음을 쉽게 알 수 있습니다. 이 경우 3 차원 계산에서 모든 소용돌이를 해결하는데 필요한 제어 볼륨의 수는 109보다 커집니다.이 크기의 수는 현재 계산 능력을 훨씬 뛰어 넘습니다. 이러한 이유로 난류에 대한 대략적인 모델을 구성하는 데 상당한 노력을 기울였습니다.

We cannot describe turbulence modeling in any detail in this short article. Instead, we will simply make some basic observations about the types of models available. Be forewarned, however, that no models exist for general use. Every model must be employed with discretion and its results cautiously treated.

이 짧은 기사에서는 난류 모델링에 대해 구체적으로 설명 할 수 없습니다.  대신 사용 가능한 모델의 유형에 대한 몇 가지 기본적인 설명만 합니다.  그러므로 일반 모델은 존재하지 않는 것을 미리 양해 바랍니다.  어떤 모델도 신중하게 선택하고 결과를 주의 깊게 처리해야 합니다.

The original turbulence modeler was Osborne Reynolds. Anyone interested in this subject should read his groundbreaking work (Phil. Trans. Royal Soc. London, Series A, Vol.186, p.123, 1895). Reynolds’s insights and approach were both fundamental and practical.

난류를 처음으로 모델링 한 인물은 Osborne Reynolds 입니다.  이 건에 관심이있는 분은 Reynolds 의 획기적인 저서 (Phil. Trans. Royal Soc. London, Series A, Vol.186, p.123,1895)를 참조하십시오.  Reynolds 의 통찰력과 접근 방식은 기본이며 동시에 실용적인 것입니다.

The Pseudo-Fluid Approximation

In a fully turbulent flow it is sometimes possible to define an effective turbulent viscosity, μeff, that roughly approximates the turbulent mixing processes contributing to a diffusion of momentum (and other properties). Thinking of a turbulent flow as a pseudo-fluid having increased viscosity leads to the observation that the effective Reynolds number for a turbulent flow is generally less than 100:

의사 유체 근사

완전 난류 흐름에서는 운동량 (및 기타 특성)의 확산에 기여하는 난류 혼합 공정에 대략적으로 근접하는 효과적인 난류 점도 μ eff를 정의 할 수 있습니다. 난류 흐름을 점도가 증가 된 유사 유체로 생각하면 난류 흐름에 대한 유효 레이놀즈 수가 일반적으로 100 미만이라는 관찰이 가능합니다.

This observation is particularly useful because it suggests a simple way to approximate some turbulent flows. In particular, when the details of the turbulence are not important, but the general mixing behavior associated with the turbulence is, it is often possible to use an effective turbulent (eddy) viscosity in place of the molecular viscosity. The effective viscosity can often be expressed as

이 관찰 결과는 몇 가지 난류를 근사하는 간단한 방법을 제시하고 있기 때문에 특히 유용합니다.  특히 난류 대한 자세한 내용은 중요하지 난류와 관련된 일반적인 혼합 거동이 중요한 경우에는 분자 점성 대신 사용 난류 (소용돌이) 점성을 사용할 수있는 경우가 있습니다.  유효 점성은 다음의 식으로 나타낼 수 있습니다.

where α is a number between 0.02 and 0.04. This expression works well for the turbulence associated with plane and cylindrical jets entering a stagnant fluid. The effective Reynolds number associated with this model is Re=1/α, a number between 25 and 50.

α는 0.02에서 0.04 사이의 숫자입니다.  이 수식은 정체 유체에 들어가는 평면 제트 및 원통형 분류 관련 난류에 대하여 효과가 있습니다.  이 모델에 대한 사용 레이놀즈 수는 Re = 1 / α 25에서 50 사이의 숫자입니다.

While this model is often adequate for predicting the gross features of a turbulent flow, it may not be suitable for predicting local details. For example, it would predict a parabolic flow (i.e., laminar) profile in a pipe instead of the measured logarithmic profile.

이 모델은 종종 난류의 전반적인 특징을 예측하는데는 적합하지만, 로컬 세부 사항을 예측하는 데는 적합하지 않을 수 있습니다.  예를 들어, 측정된 대수 프로필 대신 파이프의 포물선 흐름 (층류 등)의 프로파일을 예측합니다.

Local Viscosity Model

The next level of complexity beyond a constant eddy viscosity is to compute an effective viscosity that is a function of local conditions. This is the basis of Prandtl’s mixing-length hypothesis where it is assumed that the viscosity is proportional to the local rate of shear. The proportionality constant has the dimensions of a length squared. The square root of this constant is referred to as the “mixing length.”

This model offers an improvement over a simple constant viscosity. For example, it predicts the logarithmic velocity profile in a pipe. However, it is not used much because it doesn’t account for important transport effects.

국소 점성 모델

일정한 소용돌이 점성보다 복잡한 것은 국소적 조건의 함수인 유효 점성을 계산하는 것입니다.  이것은 점성이 국소적 전단 속도에 비례한다고 가정된다는 프란틀 혼합 길이 가설(Prandtl’s mixing-length hypothesis )의 기초가됩니다.  비례 상수의 차원은 길이의 제곱입니다.  이 상수의 제곱근은 “혼합 장”이라고합니다.

이 모델은 간단한 일정한 점성 개선을 제공합니다.  예를 들어, 파이프의 대수 속도 프로파일을 예측할 수 있습니다.  그러나 중요한 수송 효과를 지원하지 않기 때문에 그다지 많이 사용되지 않습니다.

Turbulence Transport Models

For practical engineering purposes the most successful computational models have two or more transport equations. A minimum of two equations is desirable because it takes two quantities to characterize the length and time scales of turbulent processes. The use of transport equations to describe these variables allows turbulence creation and destruction processes to have localized rates. For instance, a region of strong shear at the corners of a building may generate strong eddies, while little turbulence is generated in the building’s wake region. The strong mixing observed in the wakes of buildings (or automobiles and airplanes) is caused by the advection of upstream generated eddies into the wake. Without transport mechanisms, turbulence would have to instantly adjust to local conditions, implying unrealistically large creation and destruction rates.

난류 수송 모델

실용 공학의 목적인 가장 뛰어난 수치 모델에는 2 개 이상의 수송 방정식이 있습니다.  난류 과정의 길이와 시간의 스케일을 특징으로는 2 개 분량이 필요하므로 최소한 2 개의 방정식이있는 것이 바람직 할 것입니다.  수송 방정식을 사용하여 이러한 변수를 표현하면 난류의 생성 속도와 파괴율을 국소적으로 할 수 있습니다.  예를 들어, 건물의 모서리의 전단력이 강한 영역에서 강력한 소용돌이가 생성 된 건축물의 후류 영역에서 난류는 거의 생성되지 않습니다.  건축물 (또는 자동차 나 비행기)의 후류에서 관찰되는 강력한 혼합은 상류에서 생성된 소용돌이 후류의 이류에 의해 발생합니다.  수송 메커니즘이 없는 경우, 난류는 국소적 조건에 즉시 적응해야하므로 생성 속도와 파괴율이 비현실적인 크기입니다.

Nearly all transport models invoke one or more gradient assumptions in which a correlation between two fluctuating quantities is approximated by an expression proportional to the gradient of one of the terms. This captures the diffusion-like character of turbulent mixing associated with many small eddy structures, but such approximations can lead to errors when there is significant transport by large eddy structures.

거의 모든 수송 모델에서 하나 이상의 경사 가정을 이루어 두 변동하는 양의 상관 관계가 하나의 항 기울기에 비례하는 식으로 근사됩니다.  이를 통해 다수의 작은 소용돌이 구조와 관련된 난류 혼합 확산적인 특징을 파악할 수 있지만, 큰 소용돌이 구조에 의해 상당한 전송이 존재하는 경우, 이러한 근사 오류가 발생할 수 있습니다.

Large Eddy Simulation

Most models of turbulence are designed to approximate a smoothed out or time-averaged effect of turbulence. An exception is the Large Eddy Simulation model (or Subgrid Scale model). The idea behind this model is that computations should be directly capable of modeling all the fluctuating details of a turbulent flow except for those too small to be resolved by the grid. The unresolved eddies are then treated by approximating their effect using a local eddy viscosity. Generally, this eddy viscosity is made proportional to the local grid size and some measure of the local flow velocity, such as the magnitude of the rate of strain.

Large Eddy 시뮬레이션

난류의 대부분의 모델은 매끄럽게 또는 시간 평균된 난류의 효과를 근사하도록 설계되어 있습니다.  예외는 큰 에디 시뮬레이션 모델 (또는 서브 그리드 스케일 모델)입니다.  이 모델의 배경에는 너무 작은 격자에 의해 해결할 수 없는 것을 제외하고는 난류의 모든 변동 내용은 계산에 의해 직접 모델링 할 수 있어야 한다는 생각이 있습니다.  미해결 소용돌이는 로컬 점성을 사용하여 효과를 근사하여 처리됩니다.  일반적으로이 소용돌이 점성은 국소적인 격자 크기 및 어떤 국소적인 흐름의 속도 측정 (변형 속도의 크기 등)에 비례합니다.

대부분의 난류 모델은 난류의 평활화 또는 시간 평균 효과에 근접하도록 설계되었습니다. 예외는 Large Eddy Simulation 모델 (또는 Subgrid Scale 모델)입니다. 이 모델의 이면에있는 아이디어는 계산이 격자에 의해 해결 되기에는 너무 작은 것을 제외하고, 난류 흐름의 모든 변동 세부 사항을 직접 모델링 할 수 있어야 한다는 것입니다. 해결되지 않은 소용돌이는 로컬 소용돌이 점도를 사용하여 효과를 근사화하여 처리됩니다. 일반적으로, 이 와류 점도는 로컬 격자 크기와 변형률의 크기와 같은 로컬 유속 측정치에 비례하여 만들어집니다.

Such an approach might be expected to give good results if the unresolved scales are small enough, for example, in the viscous sub-range. Unfortunately, this is still an uncomfortably small size. When these models are used with a minimum scale size that is above the viscous sub-range, they are then referred to as Coherent Structure Capturing models.

이러한 접근 방식은 미해결 스케일이 충분히 작은 경우, 예를 들어 점성이 작은 영역에 있는 경우에 좋은 결과를 얻을 수 있을 것으로 기대됩니다.  불행히도 아직은 여전히 불편한 작은 크기 입니다.  이러한 모델을 점성 작은 영역보다 높은 최소 스케일 사이즈로 사용하는 경우는 CSC (Coherent Structure Capturing) 모델이라고합니다.

The advantage of these more realistic models is that they provide information not only about the average effects of turbulence but also about the magnitude of fluctuations. But, this advantage is also a disadvantage, because averages must actually be computed over many fluctuations, and some means must be provided to introduce meaningful fluctuations at the start of a computation and at boundaries where flow enters the computational region.

이보다 현실적인 모델의 장점은 난류의 평균 효과에 대한 정보뿐만 아니라 변동의 크기에 대한 정보도 제공 될 것입니다.  그러나 이와같은 장점은 단점도 있습니다.  평균적으로 실제로 다수의 변동에 대해 계산해야 하며, 계산의 시작 및 흐름이 계산 영역에 들어가는 경계에서 상당한 변화를 도입하기위한 수단을 제공 할 필요가 있기 때문입니다.

Turbulence from an Engineering Perspective

We have seen that it is probably not reasonable to attempt to compute all the details of a turbulent flow. Furthermore, from the perspective of most applications, it’s not likely that we would be interested in the local details of individual fluctuations. The question then is how should we deal with turbulence, when should we employ a turbulence model, and how complex should that model be?

공학적 관점에서의 난류

지금까지 난류의 모든 세부 사항을 계산하려고하는 것은 아마도 합리적이지 않다는 것을 확인했습니다.  또한 많은 적용례의 관점에서 개별 변동의 국소적인 세부 사항이 관심의 대상이 될 수는 없을 것입니다.  거기서 생기는 의문은 난류를 어떻게 처리해야 할지 난류 모델을 언제 선택할지 그 모델이 얼마나 복잡할지에 있다는 것입니다.

Experimental observations suggest that many flows become independent of Reynolds number once a certain minimum value is exceeded. If this were not so, wind tunnels, wave tanks, and other experimental tools would not be as useful as they are. One of the principal effects of a Reynolds number change is to relocate flow separation points. In laboratory experiments this fact sometimes requires the use of trip wires or other devices to induce separation at desired locations. A similar treatment may be used in a numerical simulation.

실험적 관찰에 따르면 특정 최소값이 초과되면 많은 흐름이 레이놀즈 수와 무관하게됩니다. 그렇지 않다면 풍동, 파도 탱크 및 기타 실험 도구는 그다지 유용하지 않을 것입니다. 레이놀즈 수 변경의 주요 효과 중 하나는 흐름 분리 지점을 재배치하는 것입니다. 실험실 실험에서이 사실은 때때로 원하는 위치에서 분리를 유도하기 위해 트립 와이어 또는 기타 장치를 사용해야합니다. 유사한 처리가 수치 시뮬레이션에서 사용될 수 있습니다.

Most often a simulation is done to determine the dominant flow patterns that develop in some specified situation. These patterns consist of the mean flow and the largest eddy structures containing the majority of the kinetic energy of the flow. The details of how this energy is removed from the larger eddies and dissipated into heat by the smallest eddies may not be important. In such cases the dissipation mechanisms inherent in numerical methods may alone be sufficient to produce reasonable results. In other cases it is possible to supply additional dissipation with a simple turbulence model such as a constant eddy viscosity or a mixing length assumption.

대부분의 경우 특정 상황에서 발생하는 지배적 인 흐름 패턴을 결정하기 위해 시뮬레이션이 수행됩니다. 이러한 패턴은 평균 흐름과 흐름의 대부분의 운동 에너지를 포함하는 가장 큰 소용돌이 구조로 구성됩니다. 이 에너지가 더 큰 소용돌이에서 제거되고 가장 작은 소용돌이에 의해 열로 소산되는 방법에 대한 세부 사항은 중요하지 않을 수 있습니다. 그러한 경우 수치 적 방법에 내재 된 소산 메커니즘만으로도 합리적인 결과를 얻을 수 있습니다. 다른 경우에는 일정한 소용돌이 점도 또는 혼합 길이 가정과 같은 간단한 난류 모델을 사용하여 추가 소산을 제공 할 수 있습니다.

Turbulence transport equations require more CPU resources and should only be used when there are strong, localized sources of turbulence and when that turbulence is likely to be advected into other important regions of the flow.  When there is reason to seriously question the results of a computation, it is always desirable to seek experimental confirmation.

An excellent introduction to fluid turbulence can be found in the book Elementary Mechanics of Fluids by Hunter Rouse, Dover Publications, Inc., New York (1978).

난류 전송 방정식은 더 많은 CPU 리소스를 필요로하며 강력하고 국부 화 된 난기류 소스가 있고 그 난류가 흐름의 다른 중요한 영역으로 전파 될 가능성이있는 경우에만 사용해야합니다. 계산 결과에 매우 의문이 생길 경우는 실험에 의해 확인하는 것이 좋습니다.

유체 난류에 대한 훌륭한 소개는 Hunter Rouse, Dover Publications, Inc., New York (1978)의 책 Elementary Mechanics of Fluids에서 찾을 수 있습니다.

Breakwater Structures

Breakwater Structures

복잡한 지형과 해안 구조를 통해 큰 진폭의 비선형 파를 모델링 할 수 있는 FLOW-3D의 능력은 간단한 선형 파에서 복잡한 바다 상태에 이르기까지 파도의 강제 유체 역학적 방파제 구조물의 성능을 분석 할 수 있는 이상적인 도구입니다. 수신 파로부터 운동 에너지를 제거하여 구조물 의 강도를 측정하는 것 이외에도 FLOW-3D는 사용자로 하여금 퇴적물 전송의 상세한 분석을 수행 할 수 있게 하며 방파제 요소에 작용하는 힘의 평가도 가능합니다. 생태학적 요구로 인해 통합 연안 구조가 더 중요해 지는 상황에서 FLOW-3D는 엔지니어로 하여금 구조물의 주위에 흐름 패턴의 세부 사항을 평가할 수 있습니다.

Rubble-mound Breakwaters Simulation

Rubble-mound Breakwaters Simulation

테트라포드 방파제는 종종 극심한 파동 하중의 영향으로부터 해안선을 보호하고 긴 해안 표류를 막기 위해 사용됩니다. 엔지니어는 FLOW-3D를 사용하여 보호 레이어, 파동 반사 계수 및 이러한 구조와 관련된 볼륨 overtopping volumes 을 추정할 수 있습니다.