• 한국해안·해양공학회논문집
• /
• 제32권1호
• /
• pp.69-84
• /
• 2020

자연 해빈은 해양환경에 따라 침·퇴적을 반복하며 고 파랑에 의해 해빈이 대규모로 침식되더라도 폭풍이 잦아들고 다시 너울이 우세한 해양환경이 회복되는 경우 점진적으로 복원되며, 이러한 해빈 복원은 경계층 streaming을 통해 이루어진다. 이처럼 경계층 streaming은 그 공학적 가치에도 불구하고 해안에서 가용한 표사의 대부분이 공급되는 쇄파 역에서의 경계층 streaming에 대한 우리의 이해는 아직 상당히 부족하다. 이러한 인식에 기초하여 본 연구에서는 쇄파역 경계층 streaming 수리특성을 살펴보기 위해 단조 해안과 사주를 포함한 해안에서의 천수 과정을 수치모의하였다. 수치 모의는 Spatially filtered Navier-Stokes Eq., LES(Large Eddy Simulation), Dynamic Smagorinsky 난류모형으로 구성된 정교한 수치모형에 기초하여 수행되었으며, 이 과정에서 k-ε 난류모형과 LES Turbulence Closure가 모의결과에 미치는 영향도 함께 살펴보았다. 모의결과 해안공학계에 잘 알려진 k-ε 난류모형의 한계로 인해 wall function에 기반한 k-ε 난류모형의 경우 LES와 비교하면 저면 인근 유속이 다소 과다하게 모의 되었다. 또한, 바닥과 가까운 해역에서의 유속이 바닥의 영향으로부터 비교적 자유로운 상층부에서의 유속보다 우월한 Longuet-Higgins(1957)가 이야기하는 전형적인 경계층 streaming이 천수 초입부에서부터 쇄파 역 깊숙이까지 존재하는 것을 확인하였다. 또한, 주기가 상대적으로 긴 경우 경계층 streaming의 세기와 생성범위는 해안 방향으로 확대되며 이러한 경향은 경계층 streaming이 바닥 인근에서 진행되는 마찰로 인한 파랑에너지손실로 결과되며 주기가 긴 경우 천수 과정이 일찍 시작된다는 사실을 상기하면 충분히 수용 가능해 보이며, Longuet-Higgins(1957)의 해석 해에서도 같은 경향을 확인할 수 있다.

• 대한기계학회논문집
• /
• 제19권7호
• /
• pp.1691-1701
• /
• 1995

Two aspects of Large-Eddy Simulation(LES) are investigated in order to improve its performance. The first one is on how to determine the model coefficient in conjunction with a dynamic subgrid-scale model, and the second one is on a wall-layer model(WLM) which allows one to skip near-wall regions to save a large number of grid points otherwise required. Especially, a WLM suitable for a separated flow is considered. Firstly, an averaging technique to calculate the model coefficient of dynamic subgrid-scale modeling(DSGSM) is introduced. The technique is based on the concept of local averaging, and useful to stabilize numerical solution in conjunction with LES of complex turbulent flows using DSGSM. It is relatively simple to implement, and takes very low overhead in CPU time. It is also able to detect the region of negative model coefficient where the "backscattering" of turbulence energy occurs. Secondly, a wall-layer model based on a local turbulence intensity is considered. It locally determines wall-shear stresses depending on the local flow situations including separation, and yields better predictions in separated regions than the conventional WLM. The two techniques are tested for a turbulent obstacle flow, and show the direction of further improvements.rovements.

• Nuclear Engineering and Technology
• /
• 제55권9호
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• 한국해양공학회지
• /
• 제18권2호
• /
• pp.1-9
• /
• 2004

As the first step to investigate the nonlinear interactions between turbulence and marine structures inside a viscous NWT, a LES technique was applied to solve the turbulent channel flow for =150. The employed turbulence models included 4 types: the Smagorinsky model, the Dynamic SGS model, the Structure Function model, and the Generalized Normal Stress model. The simulated data in time-series for the LESs were averaged in both time and space, and statistical analyses were performed. The results of the LESs were compared with those of a DNS, developed in the present study and two spectral methods by Yoon et al.(2003) and Kim et a1.(1987). Based on this research, the accuracy of LESs has been found to be still related to the number of grids for fine grid size).

• 한국해안·해양공학회논문집
• /
• 제32권1호
• /
• pp.55-68
• /
• 2020

너울이 우세한 온화한 해양환경에서 출현하는 beach cusp에서의 경계층 streaming 수리특성을 살펴보기 위해 edge waves의 천수 과정을 수치 모의하였다. Beach cusp을 유지하는 것으로 알려진 synchronous edge waves는 같은 주기와 파고를 지니는 두 개의 Cnoidal wave가 전면해역에서 비스듬히 조우 되도록 조파하여 재현하였다. Beach cusp의 진폭 A_B과 파장 L_B은 맹방 해변에서 수행된 관측결과를 토대로 각각 1.25 m, 18 m로 선정하였다. 모의결과 천수 각 단계에서 예외 없이 경계층 streaming을 관측할 수 있었으며 최대 경계층 streaming은 사주 정점에서 발생하였다. 주기가 가장 짧은 RUN 1의 경우 그 세기는 약 0.32 m/s 내외에 분포하며 이러한 수치는 free stream 유속 u_∞ 진폭의 두 배에 달하는 것으로 wave Reynolds 응력에 기반한 Longuet-Higgins(1957)의 해석 해와는 상당한 차이를 보였다. 수치 모의과정에서 온화한 해양환경에서 해빈이 복원되는 과정을 특정할 수 있었으며 이 과정을 정리하면 다음과 같다: 너울로 구성된 파랑 무리에서 성분 파랑 간의 공진성 상호작용으로 생성된 외 중력파가 쇄파선 인근에 도달하는 경우 중력으로 인한 가속이 더해진 Phase II 파랑 궤도 운동으로 수면 가까이 상승한 많은 모래가 쇄파 시 발생하는 파 마루로부터 시작된 up-rush에 의해 전 빈 정점 가까이 이동하며 이 과정에서 발생하는 침투로 인해 퇴적되는 것으로 모의 되었다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제30권3호
• /
• pp.369-375
• /
• 2006

Large eddy simulation has been applied to simulate turbulent flow around a cubic obstacle mounted on a channel surface for a Reynolds number of 40000(based on the incoming bulk velocity and the obstacle height) using a Smagorinsky model and a Lagrangian dynamic model. In order to develop the LES to the practical engineering application, the effect of upwind scheme, turbulent sub-grid scale model were investigated. The computed velocities. turbulence quantifies, separation and reattachment length were evaluated by compared with the previous experimental results.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.28-34
• /
• 2008

In engineering application of large eddy simulation, there are still questions as follows grid dependency on numerical results, the effect of upwind scheme against a calculation instability, appropriate boundary conditions dealing with turbulence fluctuation and the performance of SGS models. In this study, in order to develop the LES to the engineering application, large eddy simulation was carried out to investigate the effect of upwind scheme, turbulent subgrid model and the grid dependancy of the flow around a wall-mounted cube in a channel at Re=40,000 based on cubic height and bulk mean velocity. The computed velocities, turbulence quantities, separation and reattachment length were evaluated compared with the experimental results of R. Matinuzzi and C. Tropea.

• 대한기계학회논문집B
• /
• 제32권9호
• /
• pp.720-727
• /
• 2008

Time delay effects on near-wall turbulent structures are investigated by performing a large-eddy simulation of a transient turbulent flow in a pipe. To elucidate the time delay effects on the near-wall turbulence, we selected the dimensionless acceleration parameter which was used in the previous study. Various turbulent statistics revealed the distinctive features of the delay. It was shown that the dynamic Smagorinsky model is valid to capture the alterations of the turbulence physics well. A dimensionless time for the responses of the flow quantities was introduced to give the detailed information on the delay of the nearwall turbulence. The conditionally-averaged flow fields associated with Reynolds shear stress producing events show that sweep and ejections are closely related to the delays of the turbulence production and the turbulence propagation toward the pipe center. The present study suggested that the enhanced anisotropy of the turbulence in the initial and transient stages would be a challenging problem to standard turbulence models.

• Journal of Mechanical Science and Technology
• /
• 제15권12호
• /
• pp.1835-1844
• /
• 2001

The turbulent flow with wake, reattachment and recirculation is a very important problem that is related to vehicle dynamics and aerodynamics. The Smagorinsky Model (SM), the Dynamics Subgrid Scale Model (DSM), and the Lagrangian Dynamic Subgrid Scale Model (LDSM) are used to predict the three-dimensional flow field around a bluff body model. The Reynolds number used is 45,000 based on the bulk velocity and the height of the bluff body. The fully developed turbulent flow, which is generated by the driver part, is used for the inlet boundary condition. The Convective boundary condition is imposed on the outlet boundary condition, and the Spalding wall function is used for the wall boundary condition. We compare the results of each model with the results of the PIV measurement. First of all, the LES predicts flow behavior better than the k-$\xi$ turbulence model. When ew compare various LES models, the DSM and the LDSM agree with the PIV experimental data better than the SM in the complex flow, with the separation and the reattachment at the upper front part of th bluff body. But in the rear part of the bluff body, the SM agrees with the PIV experimental results better than them. In this case, the SM predicts overall flow behavior better than the DSM nd the LDSM.

• Wind and Structures
• /
• 제39권4호
• /
• pp.287-304
• /
• 2024

An accurate assessment of aerodynamic effects on structures is essential for a reliable wind design for high-rise buildings. Turbulence model is a key ingredient of computational fluid dynamics (CFD) in calculating the wind flow fields. This paper aims to identify the properties of representative RANS and LES models particularly for wind load determination. The models investigated are the realizable k-ε model for RANS and the dynamic Smagorinsky model for LES. In this study, their application aspects are discussed to provide enhanced reproducibility and reliability. The airflow around a building at Reynolds number 76,000 is simulated and the numerical results are also compared with wind tunnel experiment data. The wind design loads, such as story shear forces and overturning or torsional moments, are calculated based on the numerical results. Both RANS and LES models accurately capture surface pressure profiles, while LES results demonstrate proper energy decay in the power spectra. The numerical results highlight the effects of aspect ratio of building and the attack angle on the wind loads. This information would be of great help in designing tall buildings resilient to wind environments using CFD models.