• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.120-126
• /
• 1991

An expanded scheme of a direct numerical solution method for solving the Navier-Stokes equation considering the modified boundary conditions for gas lubrication with ultra low clearance at high .LAMBDA. region is presented. Many examples are calculated by this scheme and their results are compared to the previous solutions using P$^{2}$H$^{[-992]}$ . This scheme has the advantages of fast calculation time and stable convergence in high .LAMBDA. region, and gives very good results in the case of fluid film thickness discontinuity.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.377-381
• /
• 2006

Recently, as economy grow and population increase we need to develop our coastal area and make good use of it for various purposes. That's why large structures are being installed on the sea. Some samples are petroleum storage tanks, pier of sea bridges. These are large structures which have been installed at coastal area. When we design such vertical cylinder, we should avoid too much construction expense caused by excessive designing or by lack of sufficient design. In order to prevent excessive expenditure, it is important to correctly calculate the force of waves acting on structures and predict the wave transformation. In this study, apply to VOF method based on Navier-Stokes equation and then discussed that nonlinear wave force and wave transformation. A comparison between the numerical model and existing experimental results showed nice agreement among them.

• Journal of Navigation and Port Research
• /
• v.28 no.8
• /
• pp.753-757
• /
• 2004

The free-surface flow is simulated to make clear the viscous interaction of stem waves and the sub-breaking phenomena around a high speed vehicle. The Navier-Stokes equation is solved by a finite difference method where the body-fitted coordinate system, the wall function and the triple-grid system are invoked They are applied to study precisely on the stem flow of S-103 as to which extensive experimental data are available. Computations are extended to the submerged revolutional body. The numerical result shows that the gradient of M/Us is greatly influenced by the submerged depth And the stem wave is influenced by the separation due to the bow wave.

• Journal of applied mathematics & informatics
• /
• v.16 no.1_2
• /
• pp.383-405
• /
• 2004

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregular-shapes.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.91-98
• /
• 2000

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the flow field and compression wave around the high speed train which Is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation around the nose region was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed and presented.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.171-176
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.182-187
• /
• 2003

Numerical analysis of three-dimensional vicous flow is used to compute the design speed operating line of a transonic axial-flow compressor. The Navier-Stokes equation was solved by an explicit finite-difference numerical scheme and the Baldwin-Lomax turbulence model was applied. A spatially-varying time-step and an implicit residual smoothing were used to improve convergence. Two-stage axial compressor of a turboshaft engine developed KARI was chosen for the analysis. Numerical results show reasonably good agreements with experimental measurements made by KARI. Numerical solutions indicate that there exist a strong shock-boundary layer interaction and a subsequent large flow separation. It is also observed that the shock is moved ahead of the blade passage at near-stall condition.