전산유체역학을 이용한 FLNG의 풍하중 추정에 관한 연구

이상의 한국항해항만학회 2022 한국항해항만학회지 Vol.46 No.6

큰 상부 형상을 가지는 FLNG (Liquefied Natural Gas Floating Production Storage Offloading Units, LNG FPSOs) 등의 해양구조물은 안정적인 운동성능 확보 및 계류라인 설계에 있어 정도 높은 풍하중 추정이 필수적이다. 따라서 본 연구의 목적은 FLNG의 풍하중 추정을 위한 수치해석 기법을 개발하는 데 있다. 특히, 본 연구에서 개발한 수치해석 기법은 저자의 이전연구를 FLNG에 맞추어 수정하였다. 풍하중 추정을 위한 수치해석은 15° 간격으로 0-360° 범위에서 균일 풍속 조건과 풍속 프로파일을 적용한 NPD (Norwegian Petroleum Directorate) 조건에서 수행하였다. 먼저, NPD 모델 풍속 프로파일 모델 개발을 위해 Sand-Grain Roughness 변화에 따른 풍속 프로파일을 분석하였다. 개발된 NPD 모델을 이용하여 3가지 풍향 (Head, Quartering & Beam)에 대한 메쉬 수렴성 시험을 수행하였다. 최종적으로 개발된 NPD 모델과 메쉬를 이용하여 균일한 풍속 조건과 NPD 조건에서의 풍하중을 평가하고 비교하였다. 본 연구에서는 RANS (Reynolds-averaged Navier-Stokes) 기반 Solver인 STAR-CCM+ (17.02)를 이용하였다. 결과를 요약하면, 풍속 프로파일을 적용한 NPD 모델에서의 풍하중은 균일 풍속() 조건과 비교하여, Surge와 Yaw 하중이 최대 20.35% 와 34.27% 증가하였다. 특히, 특정 일부 구간에서만 큰 하중의 차이를 보인 Sway (, )와 Roll (, )은 구간별 평균 증가율이 15.60%와 10.89% 수준으로 나타났다. It has been noted that an accurate estimation of wind loads on offshore structures such as an FLNG (Liquefied Natural Gas Floating Production Storage Offloading Units, LNG FPSOs) with a large topside plays an important role in the safety design of hull and mooring system. Therefore, the present study aims to develop a computational model for estimating the wind load acting on an FLNG. In particular, it is the sequel to the previous research by the author. The numerical computation model in the present study was modified based on the previous research. Numerical analysis for estimating wind loads was performed in two conditions for an interval of wind direction (), 15° over the range of 0° to 360°. One condition is uniform wind speed and the other is the NPD model reflecting the wind speed profile. At first, the effect of sand-grain roughness on the speed profile of the NPD model was studied. Based on the developed NPD model, mesh convergence tests were carried out for 3 wind headings, i.e. head, quartering, and beam. Finally, wind loads on 6-degrees of freedom were numerically estimated and compared by two boundary conditions, uniform speed, and the NPD model. In the present study, a commercial RANS-based viscous solver, STAR-CCM+ (ver. 17.02) was adopted. In summary, wind loads in surge and yaw from the wind speed profile boundary condition were increased by 20.35% and 34.27% at most. Particularly, the interval mean of sway (, ) and roll (, ) increased by 15.60% and 10.89% against the uniform wind speed () boundary condition.

풍속 분포곡선이 어선의 풍하중에 미치는 영향에 관한 연구

이상의 해양환경안전학회 2020 해양환경안전학회지 Vol.26 No.7

지난 10년간 복원력 상실에 의한 어선의 해양사고가 지속해서 증가하고 있으며, 갑작스러운 강풍이 주요 원인으로 지적되고 있다. 이러한 강풍에도 견딜 수 있는 어선의 운동·조종성능을 확보하기 위해서는 정밀한 풍하중 예측 기법이 우선되어야 한다. 따라서 본 연구에서는 전산유체역학 기법을 이용한 어선의 풍하중 평가기법을 개발하고자 한다. 특히, 고도 변화에 따라 풍속이 변화하는 계산환경을 모사하여 그 결과를 균일한 속도분포를 가정한 수치해석 결과와 비교 분석하고자 한다. 본 연구에서는 0-180°까지 15° 간격으로 13개의 방향에 대해 풍하중을 계산하였으며, 계산에 사용된 메쉬 모델은 메쉬 의존성 시험을 수행하여 개발하였다. 전산수치해석은 RANS(Reynolds-averaged Navier-Stokes) 기반 상용 해석 Solver인 STAR-CCM+(Ver. 13.06)와 k-ω 난류 모델을 이용하여 정상상태(Steady State) 유동해석을 수행하였다. 수치해석결과를 간략히 살펴보면 Surge, Sway 및 Heave에서 39.5%, 41.6% 및 46.1% 풍하중이 감소하였으며 Roll, Pitch 및 Yaw에서 48.2%, 50.6% 및 36.5% 감소하였다. 결론적으로 본 연구에서는 고도에 따른 풍속 변화 모델을 통해 기존보다 정밀한 수준의 풍하중 추정이 가능한 것을 확인하였으며, 그 결과가 선박의 풍하중 추정 평가기법 발전에 이바지하길 기대한다. Marine accidents involving fishing boats, caused by a loss of stability, have been increasing over the last decade. One of the main reasons for these accidents is a sudden wind attacks. In this regard, the wind loads acting on the ship hull need to be estimated accurately for safety assessments of the motion and maneuverability of the ship. Therefore, this study aims to develop a computational model for the inlet boundary condition and to numerically estimate the wind load acting on a fishing boat. In particular, wind loads acting on a fishing boat at the wind speed profile boundary condition were compared with the numerical results obtained under uniform wind speed. The wind loads were estimated at intervals of 15° over the range of 0° to 180°, and i.e., a total of 13 cases. Furthermore, a numerical mesh model was developed based on the results of the mesh dependency test. The numerical analysis was performed using the RANS-based commercial solver STAR-CCM+ (ver. 13.06) with the k-ω turbulent model in the steady state. The wind loads for surge, sway, and heave motions were reduced by 39.5%, 41.6%, and 46.1% and roll, pitch, and yaw motions were 48.2%, 50.6%, and 36.5%, respectively, as compared with the values under uniform wind speed. It was confirmed that the developed inlet boundary condition describing the wind speed gradient with respect to height features higher accuracy than the boundary condition of uniform wind speed. The insights obtained in this study can be useful for the development of a numerical computation method for ships.

Field monitoring of boundary layer wind characteristics in urban area

Q.S. Li,Lunhai Zhi,Fei Hu 한국풍공학회 2009 Wind and Structures, An International Journal (WAS Vol.12 No.6

This paper presents statistical analysis results of wind speed and atmospheric turbulence data measured from more than 30 anemometers installed at 15 different height levels on 325 m high Beijing Meteorological Tower and is primarily intended to provide useful information on boundary layer wind characteristics for wind-resistant design of tall buildings and high-rise structures. Profiles of mean wind speed are presented based on the field measurements and are compared with empirical models’ predictions. Relevant parameters of atmospheric boundary layer at urban terrain are determined from the measured wind speed profiles. Furthermore, wind velocity data in longitudinal, lateral and vertical directions, which were recorded from an ultrasonic anemometer during windstorms, are analyzed and discussed. Atmospheric turbulence information such as turbulence intensity, gust factor, turbulence integral length scale and power spectral densities of the three-dimensional fluctuating wind velocity are presented and used to evaluate the adequacy of existing theoretical and empirical models. The objective of this study is to investigate the profiles of mean wind speed and atmospheric turbulence characteristics over a typical urban area.

Field monitoring of boundary layer wind characteristics in urban area

Li, Q.S.,Zhi, Lunhai,Hu, Fei Techno-Press 2009 Wind and Structures, An International Journal (WAS Vol.12 No.6

This paper presents statistical analysis results of wind speed and atmospheric turbulence data measured from more than 30 anemometers installed at 15 different height levels on 325 m high Beijing Meteorological Tower and is primarily intended to provide useful information on boundary layer wind characteristics for wind-resistant design of tall buildings and high-rise structures. Profiles of mean wind speed are presented based on the field measurements and are compared with empirical models' predictions. Relevant parameters of atmospheric boundary layer at urban terrain are determined from the measured wind speed profiles. Furthermore, wind velocity data in longitudinal, lateral and vertical directions, which were recorded from an ultrasonic anemometer during windstorms, are analyzed and discussed. Atmospheric turbulence information such as turbulence intensity, gust factor, turbulence integral length scale and power spectral densities of the three-dimensional fluctuating wind velocity are presented and used to evaluate the adequacy of existing theoretical and empirical models. The objective of this study is to investigate the profiles of mean wind speed and atmospheric turbulence characteristics over a typical urban area.

HeMOSU-1 풍속자료를 이용한 연직 분포함수의 매개변수 추정 및 분석

고동휘(Dong-Hui Ko),조홍연(Hong-Yeon Cho),이욱재(Uk-Jae Lee) 한국해안해양공학회 2021 한국해안해양공학회 논문집 Vol.33 No.3

다양한 목표 고도에서의 풍속 추정은 해상풍력 구조물 설계 및 풍파 추정 등의 분야에서 매우 중요한 요소이다. 그러나 풍속 관측 자료가 특정 고도에 한정되어 있기 때문에 다른 고도에서의 풍속 추정은 일반적으로 사용되는 연직 분포함수와 평균적인 매개변수를 이용하여 추정한다. 본 연구에서는 HeMOSU-1 관측타워의 다양한 고도에서 측정한 풍속 자료를 이용하여 Power 함수, 대수함수의 매개변수를 추정하고 그 변동 양상을 분석하였다. 매개변수 추정 결과, Power 함수의 지수 매개변수는 일반적으로 제안되는 0.14(= 1/7) 보다 작은 평균 0.10 정도로 추정되었으며, 변동 범위도 0.0~0.3 정도로 파악되었다. 대수분포함수의 경우, 매개변수는 마찰속도와 조도 길이로 그 범위가 풍속에 따라 차이를 보이고 있으며, 변동 범위는 각각 0~10 (m/s), 0.0~1.0 (m) 정도로 파악되었으며, 일반적으로 제시되는 범위와는 그 차이를 보이는 것으로 파악되었다. 이러한 차이는 기존의 고도 분포함수가 대기 중립 조건을 가정하고 있는 영향으로 판단되며, 보다 정확한 추정을 위해서는 대기조건을 고려한 비선형 고도분포함수의 도입이 필요하다. A wind-speed estimation at the arbitrary elevations is key component for the design of the offshore wind energy structures and the computation of the wind-wave generation. However, the wind-speed estimation of the target elevation has been carried out by using the typical functions and their typical parameters, e.g., power and logarithmic functions because the available wind speed data is limited to the specific elevation, such as 2~3m , 10 m, and so on. In this study, the parameters of the vertical profile functions are estimated with optimal and analyzed the parameter ranges using the HeMOSU-1 platform wind data monitored at the eight different locations. The results show that the mean value of the exponent of the power function is 0.1, which is significantly lower than the typically recommended value, 0.14. The values of the exponent, the friction velocity, and the roughness parameters are in the ranges 0.0~0.3, 0~10 (m/s), and 0.0~1.0 (m), respectively. The parameter ranges differ from the typical ranges because the atmospheric stability condition is assumed as the neutral condition. To improve the estimation accuracy, the atmospheric condition should be considered, and a more general (non-linear) vertical profile functions should be introduced to fit the diverse profile patterns and parameters.

Two case studies on structural analysis of transmission towers under downburst

Yang, FengLi,Zhang, HongJie Techno-Press 2016 Wind and Structures, An International Journal (WAS Vol.22 No.6

Downbursts are of great harm to transmission lines and many towers can even be destroyed. The downburst wind field model by Chen and Letchford was applied, and the wind loads of two typical transmission towers in inland areas and littoral areas were calculated separately. Spatial finite element models of the transmission towers were established by elastic beam and link elements. The wind loads as well as the dead loads of conductors and insulators were simplified and applied on the suspension points by concentrated form. Structural analysis on two typical transmission towers under normal wind and downburst was completed. The bearing characteristics and the failure modes of the transmission towers under downburst were determined. The failure state of tower members can be judged by the calculated stress ratios. It shows that stress states of the tower members were mainly controlled by 45 degree wind load. For the inland areas with low deign wind velocity, though the structural height is not in the highest wind velocity zone of downburst, the wind load under downburst is much higher than that under normal wind. The main members above the transverse separator of the legs will be firstly destroyed. For the littoral areas with high deign wind velocity, the wind load under downburst is lower than under normal wind. Transmission towers are not controlled by the wind loads from downbursts in design process.

Two case studies on structural analysis of transmission towers under downburst

FengLi Yang,HongJie Zhang 한국풍공학회 2016 Wind and Structures, An International Journal (WAS Vol.22 No.6

Downbursts are of great harm to transmission lines and many towers can even be destroyed. The downburst wind field model by Chen and Letchford was applied, and the wind loads of two typical transmission towers in inland areas and littoral areas were calculated separately. Spatial finite element models of the transmission towers were established by elastic beam and link elements. The wind loads as well as the dead loads of conductors and insulators were simplified and applied on the suspension points by concentrated form. Structural analysis on two typical transmission towers under normal wind and downburst was completed. The bearing characteristics and the failure modes of the transmission towers under downburst were determined. The failure state of tower members can be judged by the calculated stress ratios. It shows that stress states of the tower members were mainly controlled by 45 degree wind load. For the inland areas with low deign wind velocity, though the structural height is not in the highest wind velocity zone of downburst, the wind load under downburst is much higher than that under normal wind. The main members above the transverse separator of the legs will be firstly destroyed. For the littoral areas with high deign wind velocity, the wind load under downburst is lower than under normal wind. Transmission towers are not controlled by the wind loads from downbursts in design process.

균일 흐름과 지상 전단 흐름에 놓인 수평축 풍력터빈 블레이드의 공력 하중 비교

김진(Kim, Jin),유기완(Ryu, Ki-Wahn) 한국신재생에너지학회 2007 한국신재생에너지학회 학술대회논문집 Vol.2007 No.11

Recently the diameter of the 5MW wind turbine reaches 126m, and the tower height is nearly the same with the wind turbine diameter. The blade will experience periodic inflow oscillation due to blade rotation inside the ground shear flow region, that is, the inflow velocity is maximum at uppermost position and minimum at lowermost position. In this study we compare the aerodynamic data between two inflow conditions, i.e, uniform flow and normal wind profile. From the computed results all of the relative errors for oscillating amplitudes increased due to the ground shear flow effect. Especially My at hub and F_x, M_y, M_z at LSS increased enormously. It turns out that the aerodynamic analysis including the ground shear flow effect must be considered for fatigue load analysis.

Prediction of typhoon design wind speed and profile over complex terrain

W.F. Huang,Y. L. Xu 국제구조공학회 2013 Structural Engineering and Mechanics, An Int'l Jou Vol.45 No.1

The typhoon wind characteristics designing for buildings or bridges located in complex terrain and typhoon prone region normally cannot be achieved by the very often few field measurement data, or by physical simulation in wind tunnel. This study proposes a numerical simulation procedure for predicting directional typhoon design wind speeds and profiles for sites over complex terrain by integrating typhoon wind field model, Monte Carlo simulation technique, CFD simulation and artificial neural networks (ANN). The site of Stonecutters Bridge in Hong Kong is chosen as a case study to examine the feasibility of the proposed numerical simulation procedure. Directional typhoon wind fields on the upstream of complex terrain are first generated by using typhoon wind field model together with Monte Carlo simulation method. Then, ANN for predicting directional typhoon wind field at the site are trained using representative directional typhoon wind fields for upstream and these at the site obtained from CFD simulation. Finally, based on the trained ANN model, thousands of directional typhoon wind fields for the site can be generated,and the directional design wind speeds by using extreme wind speed analysis and the directional averaged mean wind profiles can be produced for the site. The case study demonstrated that the proposed procedure is feasible and applicable, and that the effects of complex terrain on design typhoon wind speeds and wind profiles are significant.

Prediction of typhoon design wind speed and profile over complex terrain

Huang, W.F.,Xu, Y.L. Techno-Press 2013 Structural Engineering and Mechanics, An Int'l Jou Vol.45 No.1

The typhoon wind characteristics designing for buildings or bridges located in complex terrain and typhoon prone region normally cannot be achieved by the very often few field measurement data, or by physical simulation in wind tunnel. This study proposes a numerical simulation procedure for predicting directional typhoon design wind speeds and profiles for sites over complex terrain by integrating typhoon wind field model, Monte Carlo simulation technique, CFD simulation and artificial neural networks (ANN). The site of Stonecutters Bridge in Hong Kong is chosen as a case study to examine the feasibility of the proposed numerical simulation procedure. Directional typhoon wind fields on the upstream of complex terrain are first generated by using typhoon wind field model together with Monte Carlo simulation method. Then, ANN for predicting directional typhoon wind field at the site are trained using representative directional typhoon wind fields for upstream and these at the site obtained from CFD simulation. Finally, based on the trained ANN model, thousands of directional typhoon wind fields for the site can be generated, and the directional design wind speeds by using extreme wind speed analysis and the directional averaged mean wind profiles can be produced for the site. The case study demonstrated that the proposed procedure is feasible and applicable, and that the effects of complex terrain on design typhoon wind speeds and wind profiles are significant.