• Ocean Systems Engineering
• /
• 제6권1호
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• 해양환경안전학회지
• /
• 제26권7호
• /
• pp.922-930
• /
• 2020

지난 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 % 감소하였다. 결론적으로 본 연구에서는 고도에 따른 풍속 변화 모델을 통해 기존보다 정밀한 수준의 풍하중 추정이 가능한 것을 확인하였으며, 그 결과가 선박의 풍하중 추정 평가기법 발전에 이바지하길 기대한다.

• 한국항해항만학회지
• /
• 제46권6호
• /
• pp.491-500
• /
• 2022

큰 상부 형상을 가지는 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 모델에서의 풍하중은 균일 풍속(10m/s) 조건과 비교하여, Surge와 Yaw 하중이 최대 20.35 % 와 34.27% 증가하였다. 특히, 특정 일부 구간에서만 큰 하중의 차이를 보인 Sway (45°< α < 135°, 225°< α < 315°)와 Roll (60° < α < 135°, 225° < α < 270°)은 구간별 평균 증가율이 15.60%와 10.89% 수준으로 나타났다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제20회 춘계학술대회 논문집
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.