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

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.96-101
• /
• 1998

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.141-146
• /
• 1998

An Euler/Navier Stokes solver has been developed for the analysis of steady and unsteady flows. The $q-{\omega}$ turbulent model has been incorporated into the solver in strongly coupled manner for stability and robustness. A new Chimera hole cutting algorithm, Cut-paste algorithm, has been devised for automatic Chimera hole cutting. Number of viscous/inviscid numerical computations demonstrate the accuracy and the versatility of the solver.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.59-64
• /
• 2002

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.54-61
• /
• 2015

In the present study, a flow solver using a kinetic BGK scheme was developed for the compressible Navier-Stokes equation. The kinetic BGK scheme was used to simulate flow field from the continuum up to the transitional regime, because the kinetic BGK scheme can take into account the statistical properties of the gas particles in a non-equilibrium state. Various numerical simulations were conducted by the present flow solver. The laminar flow around flat plate and the hypersonic flow around hollow cylinder of flare shape in the continuum regime were numerically simulated. The numerical results showed that the flow solver using the kinetic BGK scheme can obtain accurate and robust numerical solutions. Also, the present flow solver was applied to the hypersonic flow problems around circular cylinder in the transitional regime and the results were validated against available numerical results of other researchers. It was found that the kinetic BGK scheme can similarly predict a tendency of the flow variables in the transitional regime.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.41-52
• /
• 1999

이 논문에서는, QUICK해법의 불안정성을 개량하므로써, 수치계산에 있어서 수렴이 빠르고, 수치적으로 안정한 계산을 할 수 있는 새로운 MQUICK 상류해법을 제안하고, 이를 비압축성 층류유동의 계산에 적용하였다. 또한, 해법의 정확성, 안정성, 수렴속도에 대한 검토를 통하여 본 MQUICK 상류해법의 유효성과 타당성이 평가되었다. 이 해법에서는 인공산일의 가감을 조절하기 위하여 가중계수 α를 써서 정식화 하였고, 위의 검토를 통하여 α의 최적값을 조사하였다. 이 해법을 SMAC 음해법에 적용하여 2 차원 공동유동, 3 차원 덕트유동과 같은 몇몇 표준문제를 계산하고, 계산된 결과를 실험값 또는, 3 차 정확도의 상류해법 및 QUICK해법에 의한 결과 들과 비교 하므로써, 본 MQUICK 상류해법이 위의 다른 해법에 비하여 안정하고, 유효성이 높은 해법임을 확인 하였다.

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

A multiblock grid generation has been applied to a Navier-Stokes design procedure of a axial-flow compressors. A multiblock structure simplifies the creation of structured H-grids about complex turbomachinery geometries and facilitate the creation of a grid in the tip flow region. The numerical algorithm adopts the combination of the algebraic and elliptic method to create the internal grids efficiently and quickly. The input module is made of the results of the preliminary design, i.e., flow-path, aerodynamic conditions along the spanwise direction, and the blade profile data. The final grids generated from each module of the system are used as the preprocessor for the performance prediction of the sectional blade, the blade-stacking process and the three-dimensional flow simulation inside the blade passage. Application to the blade design of the LP compressor was demonstrated to be very reliable and practical in support of design activities. This customized system are coupled strongly with the design procedure of the turbomachinery cascades using the Navier-Stokes technique.

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

The airfoil design optimization procedures based on the Navier-Stokes equations were developed, This procedure enables more realistic and practical transonic airfoil designs. The modified Hicks-Henne functions were used to generate the shape of airfoils. Five Hick-Henne functions were used to design upper surface of airfoil only. To enhance the ability of Hick-Henne function to generate various airfoil shape with limited number of functions, the positions of control points were adjusted through optimization procedure. The design procedure was applied to the single-point design for the drag minimization problem with lift and area constraints. The result shows the capability of the procedure to generate much realistic airfoils with very small drag-creep in the low transonic regime. This is mainly due to the viscosity effect of Navier-Stokes flow analysis. However, in the higher transonic range tile drag-creep appears. The multi-point design is shown to be an effective way to avoid the drag-creep and improve off-design performance which is very similar in the Euler design.