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

• Journal of computational fluids engineering
• /
• v.5 no.3
• /
• pp.23-31
• /
• 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 viscous 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 owing to the viscous interaction around the train 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.541-550
• /
• 1997

Computational and experimental investigations on the three-dimensional flowfield through an automotive cooling fan are carried out in this work. Steady, incompressible, three-dimensional, turbulent flow through a rotating axial-flow fan is analyzed with Reynolds averaged Navier-Stokes equations and standard k-.epsilon. turbulence model. The governing equations are discretized with finite-volume approximations in non-orthogonal curvilinear coordinates. Computational static pressures on the casing wall agree well with the experimental data which are measured in this work. And, they are sensitive to the change of tip clearance. The flowfield is not significantly affected by the thickness of the blade. The k-.omega. model gives the static pressure rise on the casing wall which is similar to that with the k-.epsilon. model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.82-91
• /
• 1999

A numerical analysis based on the three-dimensional Reynolds-averaged Navier-Stokes equation has been conducted to investigate the flow within a NASA rotor 67 transonic fan. General coordinate transformations are used to represent the complex blade geometry and an H-type grid is used. The governing equations are solved using implicit LU-SGS scheme for the time-marching integration and a standard ${\kappa}-{\varepsilon}$ model is used with wall functions for the turbulence modeling. The computations are compared with the experimental data and a detailed study of the flow structures near peak efficiency and near stall is presented. The calculated overall aerodynamic efficiency and three-dimensional shock system agree well with the laser anemometer data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.9
• /
• pp.1325-1334
• /
• 1998

The objectives of the present study are to develop a systematic method rather than a conventional trial-and-error method for an optimal shape design using a mathematical theory, and to apply it to engineering problems. In the present study, an optimal condition for a minimum pressure loss in a two-dimensional curved duct flow is derived and then an optimal shape of the curved duct is designed from the optimal condition. In the design procedure, one needs to solve the adjoint Navier-Stokes equations which are derived from the Navier-Stokes equations and the cost function. Therefore, a computer code of solving both the Navier-Stokes and adjoint Navier-Stokes equations together with an automatic grid generation is developed. In a curved duct flow, flow separation occurs due to an adverse pressure gradient, resulting in an additional pressure loss. Optimal shapes of a curved duct are obtained at three different Reynolds numbers of 100, 300 and 800, respectively. In the optimally shaped curved ducts, the separation region does not exist or is significantly reduced, and thus the pressure loss along the curved duct is significantly reduced.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.1005-1012
• /
• 2000

Numerical optimization techniques combined with a three-dimensional thin-layer Navier-Stokes solver are presented to find an optimum shape of a stator blade in an axial compressor through calculations of single stage rotor-stator flow. Governing differential equations are discretized using an explicit finite difference method and solved by a multi-stage Runge-Kutta scheme. Baldwin-Lomax model is chosen to describe turbulence. A spatially-varying time-step and an implicit residual smoothing are used to accelerate convergence. A steady mixing approach is used to pass information between stator and rotor blades. For numerical optimization, searching direction is found by the steepest decent and conjugate direction methods, and the golden section method is used to determine optimum moving distance along the searching direction. The object of present optimization is to maximize efficiency. An optimum stacking line is found to design a custom-tailored 3-dimensional blade for maximum efficiency with the other parameters fixed.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.66-72
• /
• 1997

Circular-to-rectangular transition ducts are used as exhaust components of high performance fighter aircraft with vectored thrust nozzles. Three-dimensional incompressible Navier-Stokes solver is used to analyze the transition duct. Cross sections of transition duct are defined by superelliptic equation. The grid system is generated by Non-Uniform Rational B-Spline, after generating surface grid by blending the cross sections. Good agreement between the results of the computational simulation and the experimental data is observed.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.