• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.12
• /
• pp.1075-1081
• /
• 2007

The temperature preconditioning is applied to the Navier-Stokes equations. Also, a new concept of diffusion Mach numbers is introduced to modify the reference Mach number for the Navier-Stokes equations. Flows over a circular cylinder were calculated at different Reynolds numbers. It is shown that the temperature preconditioning improves the convergence characteristics of Navier-Stokes equations. Also, it is shown that the modified reference Mach number alleviates the convergence problems at locally low speed regions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.5
• /
• pp.755-761
• /
• 1987

Two-Dimensional incompressible laminar boundary layer with the reversed flow region is computed using the parially parabolized Navier-Stokes equations in primitive variables. The velocities and the pressure are explicity coupled in the difference equation and the resulting penta-diagonal matrix equations are solved by a streamwise marching technique. The test calculations for the trailing edge region of a finite flat plate and Howarth's linearly retarding flows demonstrate that the method is accurate, efficient and capable of predicting the reversed flow region.

• Tunnel and Underground Space
• /
• v.13 no.5
• /
• pp.389-396
• /
• 2003

It is shown that the cubic law can be modified regarding the steady-state Navier-Stokes equations by using perturbation approximation method for a sinusoidal aperture variation. In order to adopt the perturbation theory, the sinusoidal function needs to be non-dimensionalized for the amplitude and wavelength. Then, the steady-state Navier-Stokes equations can be solved by expanding the non-dimensionalized stream function with respect to the small value of the parameter (the ratio of the mean aperture to the wavelength), together with the continuity equation. From the approximate solution of the Navier-Stokes equations, the basic cubic law is successfully modified for the steady-state condition and a sinusoidal aperture variation. A finite difference method is adopted to calculate the pressure within a fracture model, and the results of numerical experiments show the accuracy and applicability of the modified cubic law. As a result, it is noted that the modified cubic law, suggested in this study, will be used for the analysis of fluid flow through aperture geometry of sinusoidal distributions.

• Economic and Environmental Geology
• /
• v.52 no.4
• /
• pp.291-297
• /
• 2019

Fluid flow through rock fractures has been quantified using equations such as Stokes equations, Reynolds equation (or local cubic law), cubic law, etc. derived from the Navier-Stokes equations under the assumption that linear flow prevails. Therefore, these simplified equations are limited to linear flow regime, and cause errors in nonlinear flow regime. In this study, causal mechanism of nonlinear flow and critical Reynolds number were presented by carrying out fluid flow modeling with both the Navier-Stokes equations and the Stokes equations for a three-dimensional rough-walled rock fracture. This study showed that flow regimes changed from linear to nonlinear at the Reynolds number greater than 10. This is because the inertial forces, proportional to the square of the fluid velocity, increased enough to overwhelm the viscous forces. This tendency was also shown for the unmated (slightly sheared) rock fracture. It was found that nonlinear flow was caused by the rapid increase in the inertial forces with increasing fluid velocity, not by the growing eddies that have been ascribed to nonlinear flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.111-116
• /
• 2005

Two coupling methods for the Navier-Stokes equations and a two-equation turbulence model equations are compared. They are the strongly coupled method and the loosely coupled method. The strongly coupled method solves the Navier-Stokes equations and the two-equation turbulence model equations simultaneously, while the loosely coupled method solves the Navier-Stokes equation with the turbulence viscosity fixed and subsequently solves the turbulence model equations with all the flow quantities fixed. In this paper, performances of two coupling methods are compared for two and three-dimensional problems.

• Journal of the Korean Chemical Society
• /
• v.40 no.6
• /
• pp.409-414
• /
• 1996

The limit of applicability of Navier-Stokes equation to stationary plane shock-waves is examined by using the principle of minimum entropy production of linear irreversible thermodynamics. In order to obtain analytic results, the equation is linearized near the equilibrium of downstream. Results show that the solution of Navier-Stokes equation which fits the boundary condition of far downstream flow is consistent with the thermodynamic requirement within the first order when the solution is expanded around the M=1, where M is the Mach number of upstream speed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.2
• /
• pp.41-49
• /
• 2019

We approximately obtain heights of cnoidal waves overtopping on a porous breakwater using both the one-layer Boussinesq equations (Vu et al., 2018) and the two-layer Boussinesq equations (Huynh et al., 2017). For cnoidal waves overtopping on a porous breakwater, we find through numerical experiments that the heights of cnoidal waves overtopping on a low-crested breakwater (obtained by the Navier-Stokes equations) are smaller than the heights of waves passing through a high-crested breakwater (obtained by the one-layer Boussinesq equations) and larger than the heights of waves passing through a submerged breakwater (obtained by the two-layer Boussinesq equations). As the cnoidal wave nonlinearity becomes smaller or the porous breakwater width becomes narrower, the heights of transmitting waves obtained by the one-layer and two-layer Boussinesq equations become closer to the height of overtopping waves obtained by the Navier-Stokes equations.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.59-59
• /
• 2011

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.1-11
• /
• 2004

A comprehensive study has been made for the investigation of the convergence and stability characteristics of the LU scheme for solving the Navier-Stokes equations on unstructured meshes. For this purpose the characteristics of the LU scheme was initially studied for a scalar model equation. Then the analysis was extended to the Navier-Stokes equations. It was shown that the LU scheme has an inherent stiffness in the streamwise direction. This stiffness increases when the grid aspect ratio becomes high and the cell Reynolds number becomes small. It was also shown that the stiffness related to the grid aspect ratio can be effectively eliminated by performing proper subiteration. The results were validated for a flat-plate turbulent flow.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.326-334
• /
• 2017

We obtain height of waves overtopping on a porous breakwater using both the one-layer and two-layer Boussinesq equations. The one-layer Boussinesq equations of Lee et al. (2014) are used and the two-layer Boussinesq equations are derived following Cruz et al. (1997). For solitary waves overtopping on a porous breakwater, we find through numerical experiments that the height of waves overtopping on a low-crested breakwater (obtained by the Navier-Stokes equations) are smaller than the height of waves passing through a high-crest breakwater (obtained by the one-layer Boussinesq equations) and larger than the height of waves passing through a submerged breakwater (obtained by the two-layer Boussinesq equations). As the wave nonlinearity becomes smaller or the porous breakwater width becomes narrower, the heights of transmitting waves obtained by the one-layer and two-layer Boussinesq equations become closer to the height of overtopping waves obtained by the Navier-Stokes equations.