• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2008년도 학술발표회 논문집
• /
• pp.1624-1629
• /
• 2008

The governing equations in generalized curvilinear coordinates for a 3D pulsatile flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. The computational technique implements the implicit approximate factorization method of the Beam and Warming method (1978), which is the extension of the Alternate Direction Implicit (ADI) method. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Rindt & Steenhoven, 1991).

• 대한기계학회논문집B
• /
• 제24권10호
• /
• pp.1335-1342
• /
• 2000

Three-dimensional flow analysis and numerical optimization methods are presented for the design of an axial-flow fan. Steady, incompressible, three-dimensional Reynolds-averaged Navier-Stokes equations are used as governing equations, and standard k- ${\varepsilon}$ turbulence model is chosen as a turbulence model. Governing equations are discretized using finite volume method. Steepest descent method, conjugate gradient method and BFGS method are compared to determine the searching directions. Golden section method and quadratic fit-sectioning method are tested for one dimensional search. Objective function is defined as a ratio of generation rate of the turbulent kinetic energy to pressure head. Two variables concerning sweep angle distribution are selected as the design variables. Performance of the final fan designed by the optimization was tested experimentally.

• 대한수학회지
• /
• 제57권4호
• /
• pp.915-933
• /
• 2020

A two-level finite element method based on the Newton iterative method is proposed for solving the Navier-Stokes/Darcy model. The algorithm solves a nonlinear system on a coarse mesh H and two linearized problems of different loads on a fine mesh h = O(H^4-𝜖). Compared with the common two-grid finite element methods for the considered problem, the presented two-level method allows for larger scaling between the coarse and fine meshes. Moreover, we prove the stability and convergence of the considered two-level method. Finally, we provide numerical experiment to exhibit the effectiveness of the presented method.

• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.34-40
• /
• 2005

본 연구에서는 비정상 유동해석을 위한 CFD 코드의 개발을 위해 대각화 ADI 기법을 적용한 정상 해석기법과 내재적 이중시간 전진기법을 도입하였다. 정상상태 Navier-Stokes 방정식의 Jacobian 행렬은 비점성항에 대해서만 적용하였고 여기에 내재적 인공점성 연산자를 첨가하여 블록 5대각 행렬을 유도하였다. 시간단축을 위해 스칼라 5대각 행렬로 대체하였다. 가상시간에 대한 정상상태기법에 실시간에 대한 미분항이 포함된 새로운 잔류항을 정의하였다. 가상시간에 대해 수렴된 해로부터 실시간 해를 구하고 시간에 대해 적분을 수행하는 내재적 이중시간 전진기법을 이용한 비정상 Navier-Stokes 코드를 개발하였다. 이에 대한 검증으로 정지한 유체속에 진동하는 평판문제, 원기둥 후방의 주기적인 Karman 와류생성, 이중원호 익형주위의 충격파 진동문제등을 수치해석하여 이론치, 실험치, 타연구자의 계산결과와 비교, 분석하였다.

• 한국구조물진단유지관리공학회 논문집
• /
• 제22권5호
• /
• pp.13-22
• /
• 2018

이 연구는 이방성 평판의 휨 해석을 위한 지배방정식을 유도하고 다양한 경계조건을 갖는 평판의 정확한 풀이과정을 제시하였다. 이 해법은 삼각급수를 이용하여 미분 방정식을 대수학적 방정식으로 변환시키는 전통적인 Navier와 Levy의 방법을 따랐다. Levy의 방법을 이용해 해를 구하려면 평판의 마주보는 두 끝단이 단순지지단인 경우에만 가능하다. Navier의 방법은 사각평판의 네 끝단이 모두 단순지지단 이어야 한다. 본 연구는 Navier와 Levy해법이 갖는 경계조건 한계를 극복하였다. 이 해법은 평판 네 끝단의 경계조건이 단순지지단과 고정단의 어떤 조합이라도 적용될 수 있다. 하중조건도 분포하중, 부분하중 그리고 선하중에 대해 적용할 수 있다. 이 해법의 장점은 Navier와 Levy해법이 갖는 경계조건 한계를 극복하였을 뿐만 아니라 정확한 해를 구할 수 있다. 비대칭 경계조건을 갖는 이방성평판에 대하여 이 해법을 이용한 계산결과를 나타냈다. 또한 Navier해법과 Levy해법 그리고 Szilard의 계산결과와 비교를 보여주었는데 계산된 처짐량이 잘 일치한다.

• 대한기계학회논문집B
• /
• 제25권11호
• /
• pp.1581-1593
• /
• 2001

An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제14권1호
• /
• pp.43-56
• /
• 2010

The representative numerical algorithms to solve the time dependent Navier-Stokes equations are projection type methods. Lots of projection schemes have been developed to find more accurate solutions. But most of projection methods [4, 11] suffer from inconsistency and requesting unknown datum. E and Liu in [5] constructed the gauge method which splits the velocity $u=a+{\nabla}{\phi}$ to make consistent and to replace requesting of the unknown values to known datum of non-physical variables a and ${\phi}$. The errors are evaluated in [9]. But gauge method is not still obvious to find out suitable combination of discrete finite element spaces and to compute boundary derivative of the gauge variable ${\phi}$. In this paper, we define 4 gauge algorithms via combining both 2 decomposition operators and 2 boundary conditions. And we derive variational derivative on boundary and analyze numerical results of 4 gauge algorithms in various discrete spaces combinations to search right discrete space relation.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.64-70
• /
• 1996

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.35-43
• /
• 1997

An efficient inverse method for 1.he supersonic/hypersonic axisymmetric body design is developed for the parabolized Navier-Stokes equations. The developed method is examined numerically for three extreme testcases in the supersonic(M/sub ∞/=3.0) and hypersonic(M/sub ∞/=6.28) speeds. The first one is a negative pressure distribution near a vacuum pressure and the second one is a positive pressure distribution over the whole region of the body. The last one is the case of abrupt change of pressure distribution to zero in the forward region of the body. These testcases show the robustness of the method. By introducing a regular-falsi method and by using a not-fully converged inverse solution, the convergence behavior was greatly improved.