• 대한조선학회논문집
• /
• 제39권3호
• /
• pp.28-40
• /
• 2002

지면효과를 받는 2차원 날개 주위의 난류유동을 비압축성 RANS(Reynolds Averaged Navier Stokes) 방정식과 유한차분법(Finite Difference Method)을 이용하여 해석하였다. 높은 레이놀즈수에 효과적인 Baldwin-Lomax 난류모델을 사용하였다. 본 연구의 목적은 지면효과를 받는 2차원 날개단면에서의 각기 다른 두 바닥 경계조건(이동지면, 고정지면)에 따른 유동의 특성을 파악하는 것이다. Clark-Y(t/C 11.7%)날개단면의 계산 결과와 발표된 계산결과 및 실험 값과의 비교를 통해 본 수치해석 프로그램의 정확성을 검증하였다. NACA4412 날개단면에 대해 지면과의 높이변화에 대해서 두 바닥 경계조건에 대해서 유동해석을 수행하였다 계산결과에 의하면 이동지면과 고정지면에 대해서 양력과 모멘트는 별 차이가 없으나 항력은 고정지면의 경우가 이동지면의 경우보다 다소 작았다. 따라서 풍동시험에서 고정지면의 결과는 이동지면에 비해 상대적으로 저항이 낮게 평가될 가능성이 있다고 본다.

• 대한토목학회논문집
• /
• 제28권5B호
• /
• pp.575-589
• /
• 2008

본 연구는 지진해일단파(tsunami bore) 혹은 조석단파(tidal bore)와 같은 단파의 동수역학적인 거동특성을 검토할 목적으로, 댐파괴류에서 단파의 형성과 동일한 방법, 즉 수조의 한쪽 끝단에 있는 고수위의 저수조(貯水槽) 게이트를 순간적으로 제거하는 방법으로 단파를 발생시킨다. 이러한 단파의 형성과 전파에 관한 수치시뮬레이션에 이상유(二相流)모델에 기초한 Navier-Stokes식을 적용하며, 이 때 비압축성 및 비혼합성의 액체와 기체흐름을 각각 고려한다. 기체와 액체의 접면을 VOF법으로 추적하고, Navier-Stokes방정식을 수치적으로 풀기 위하여 MCIP법을 적용한다. 1차원인 CIP법을 분할스텝기법을 사용하여 고차원으로 확장한 MCIP법은 수치확산이 매우 작고, 또한 안정된 스킴으로 알려져 있다. 게다가, 난류를 시뮬레이션하기 위하여 그의 유용성이 잘 알려져 있는 LES모델을 사용한다. 단파의 형성과 전파에 관한 수치해석결과를 검증하기 위하여 수리실험을 수행하였으며, 시간경과에 따른 수위변동과 평균유속변동에 대한 수치해석결과 및 실험결과를 비교하여 매우 양호한 상호대응관계를 확인할 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.225-234
• /
• 1995

터빈익렬내부의 유동해석을 위해 비압축성 점성유동해석을 이용한 수치 해석 프로그램을 개발하였다. 지배방정식으로는 2차원의 비정상 비압축성 Navier-Stokes 방정식을 일반화된 곡선좌표계로 전환하여 암시적으로(implicitly) 반복적인 시간진행방법을 이용하여 유동해석을 하였다. 지배방정식의 각항들은 시간에 대해 1차의 정확도 그리고 영역에 대해서는 2차의 정확도, 대류항에 대해서는 3차의 정확도를 가지는 Upwind기법을 적용하였다. 특히, 실험적 접근이 매우 어려운 터빈의 정익과 회전하고 있는 동익과의 상호운동을 멀티블럭기법과 데이터 interface를 통해 보다 쉽게 해석할 수 있었다. 본 연구결과는 정익만을 계산한 타 연구자의 결과와의 비교시 매우 일치하였으며 물리적인 유동을 잘 파악할 수 있었다. 난류유동 해석을 위해서 Baldwin-Lomax 모델을 적용하였다.

• 대한조선학회논문집
• /
• 제54권1호
• /
• pp.18-25
• /
• 2017

This study deals with numerically comparing performance according to rudder shape called 'Twisted rudder and Wavy twisted rudder'. In comparison with conventional rudder, rudder with wavy shape has showed a better performance at high angles of attack($30^{\circ}{\sim}40^{\circ}$) due to delaying stall. But most of study concerned with wavy shape had been performed in uniform flow condition. In order to identify the characteristics behind a rotating propeller, the present study numerically carries out an analysis of resistance and self-propulsion for KCS with twisted rudder and wavy twisted rudder. The turbulence closure model, Realizable $k-{\epsilon}$, is employed to simulate three-dimensional unsteady incompressible viscous turbulent and separation flow around the rudder. The simulation of self-propulsion analysis is performed in two step, because of finding optimization case of wavy shape. The first step presents there are little difference between twisted rudder and case of H_0.65 wavy twisted rudder in delivered power. So two kind of rudders are employed from first step to compare lift-to-drag ratio and torque at high angles of attack. Consequently, the wavy twisted rudder is presented as a possible way of delaying stall, allowing a rudder to have a better performance containing superior lift-to-drag ratio and torque than twisted rudder at high angles of attack. Also, as we indicate the flow visualization, check the quantity of separation flow around the rudder.

• 대한기계학회논문집B
• /
• 제37권11호
• /
• pp.999-1007
• /
• 2013

본 연구에서는 다양한 형상의 파이프에 대한 압력강하와 열전달 특성을 수치적으로 해석하였다. 원형 파이프에서부터 타원형, 톱니형, 비틀어진 형태와 같은 다양한 형상의 파이프를 3차원으로 수치해석을 통해 비교하였다. 수치해석은 층류에서 난류영역까지 계산을 수행하였다. 파이프 유동해석은 완전발달된 영역에서 정상상태, 비압축성 RANS수식을 이용하여 계산하였다. 유동의 손실은 friction factor를 통해 비교하였고, 열전달 성능은 각 파이프 표면에서의 Colburn factor를 통해 비교하였다. 종합적인 열유동 성능평가는 Volume and Area goodness factor를 통해 평가하였다. 열전달 성능을 향상시키고 유동의 손실은 최소화하는 최적의 형상을 연구하였다.

• 해양환경안전학회지
• /
• 제22권6호
• /
• pp.750-757
• /
• 2016

본 연구에서는 파 주기와 기울기에 따른 단일 실린더와 다중 실린더의 wave run-up을 추정하였다. 3차원 비압축성 점성 난류 유동이며 서로 다른 상을 가지는 이상유체에 대한 수치해석을 수행하기 위하여, 본 연구에서는 전산 유체 역학 상용코드인 "STAR-CCM+"을 기반으로 하여 VOF기법과 realizable $k-{\varepsilon}$ 난류 모델을 사용하였다. 모델스케일에 대한 파 주기는 단일 실린더의 경우 1.269초와 1.692초이고, 다중 실린더의 경우 1.716초이다. 각 케이스 별로 1/30, 1/16의 파 기울기를 가진다. 최종적으로 파 기울기와 주기에 대한 wave run-up 추정 결과는 관련 실험 결과와 비교하였다. 수치해석 결과는 실험과 비교하여 정성적으로 유사함을 확인하였다.

• 한국전산응용수학회:학술대회논문집
• /
• 한국전산응용수학회 2003년도 KSCAM 학술발표회 프로그램 및 초록집
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Wind and Structures
• /
• 제34권1호
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• Wind and Structures
• /
• 제20권4호
• /
• pp.549-564
• /
• 2015

The three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations and k-${\varepsilon}$ double equations turbulent model were used to investigate the effect on the measurements of anemometers due to a passing high-speed train. Sliding mesh technology in Fluent was utilized to treat the moving boundary problem. The high-speed train considered in this paper was with bogies and inter-carriage gaps. Combined with the results of the wind tunnel test in a published paper, the accuracy of the present numerical method was validated to be used for further study. In addition, the difference of slipstream between three-car and eight-car grouping models was analyzed, and a series of numerical simulations were carried out to study the influences of the anemometer heights, the train speeds, the crosswind speeds and the directions of the induced slipstream on the measurements of the anemometers. The results show that the influence factors of the train-induced slipstream are the passing head car and tail car. Using the three-car grouping model to analyze the train-induced flow is reasonable. The maxima of horizontal slipstream velocity tend to reduce as the height of the anemometer increases. With the train speed increasing, the relationship between $V_{train}$ and $V_{induced\;slipstream}$ can be expressed with linear increment. In the absence of natural wind conditions, from the head car arriving to the tail car leaving, the induced wind direction changes about $330^{\circ}$, while under the crosswind condition the wind direction fluctuates around $-90^{\circ}$. With the crosswind speed increasing, the peaks of $V_X,{\mid}V_{XY}-V_{wind}{\mid}$ of the head car and that of $V_X$ of the tail car tend to enlarge. Thus, when anemometers are installed along high-speed railways, it is important to study the effect on the measurements of anemometers due to the train-induced slipstream.