• 제목/요약/키워드: Stokes equations

검색결과 1,394건 처리시간 0.248초

회전원판 근처에서 회전하는 유연디스크에 대한 실험 및 수치해석 (Experimental and Numerical Study on an Air-Stabilized Flexible Disk Rotating Close to a Rigid Rotating Disk)

  • 가드압델라솔;임윤철
    • 정보저장시스템학회논문집
    • /
    • 제5권1호
    • /
    • pp.19-35
    • /
    • 2009
  • The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

  • PDF

미세채널내 증발을 고려한 두 유체간 열전달현상에 대한 해석적인 연구 (An Analytical Study on a Heat Transfer Mechanism with Boiling Effect between Two Fluids in a Mini-channel)

  • 유영준;최상민
    • 한국추진공학회지
    • /
    • 제17권2호
    • /
    • pp.114-121
    • /
    • 2013
  • 미세채널을 갖는 증발형 열교환기의 효율을 평가하기 위하여, 공기의 온도와 물의 온도와 같은 열교환기의 상태값들을 계산하기위한 관계식들이 문제를 단순화하기 위한 몇 가지 가정을 적용한 Navier-Stokes 방정식으로부터 유도되었다. 미세채널내부는 물의 상태에 따라 3가지 영역으로 나누었다. 이 연구의 결과로써, 미세유로를 갖는 증발형 열교환기의 증발시작시점과 건조완료점을 계산하는 방정식이 제시되었다. 본 연구결과는 증발효과를 이용하는 미세채널형 열교환기의 설계, 성능예측 및 시험결과 분석 등에 효과적으로 활용될 수 있을 것으로 기대된다.

2차 정확도 VOF기법을 활용한 수중구조물에 의한 파랑변화 예측 (Numerical Simulation of Wave Deformation due to a Submerged Structure with a Second-order VOF Method)

  • 하태민;조용식
    • 한국방재학회 논문집
    • /
    • 제10권1호
    • /
    • pp.111-117
    • /
    • 2010
  • 수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험결과와 비교 및 분석하였다.

비압축성 Navier-Stokes 방정식의 해석을 위한 MQUICK 상류해법 (MQUICK Upwind Scheme for the Incompressible Navier-Stokes Equations)

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

  • PDF

진동하는 날개의 3차원 동실속에 관한 수치해석 (Numerical Analysis of Three Dimensional Dynamic Stall of an Oscillating Wing)

  • 고성호
    • 대한기계학회논문집B
    • /
    • 제25권1호
    • /
    • pp.62-70
    • /
    • 2001
  • Three-dimensional dynamic stall over an oscillating wing has been analyzed by using a compressible Navier-Stokes code. The code solved the thin-layer Navirer-Stokes equations with a second-order time accuracy for a semispan wing with 0.3048m chord, a NACA 0015 airfoil section, and zero twist Computations were made for a freestream Mach number of 0.29, a chord Reynolds number of 1.95$\times$10(sup)6 and a reduced frequency equal to 0.1. Numerical results were compared with experimental data which include the hysteresis of lift, drag and moment at various wing span. The comparison reveals the quantitative as well as qualitative nature of the three-dimensional dynamic stall.

Slip flow 영역에서 Navier Stokes 방정식의 해석 연구 (Solutions of the Navier-Stokes equation in slip flow region)

  • 박원희;김태국
    • 대한기계학회:학술대회논문집
    • /
    • 대한기계학회 2000년도 추계학술대회논문집B
    • /
    • pp.597-602
    • /
    • 2000
  • In a MEMS(micro-electro mechanical system), the fluid may slip near the surface of a solid and have a discontinuous temperature profile. A numerical prediction in this slip flow region can provide a reasonable guide for the design and fabrication of micro devices. The compressible Navier-Stokes equation with Maxwell/smoluchowski boundary condition is solved for two simple systems; couette flow and pressure driven flow in a long channel. We found that the couette flow could be regarded as an incompressible system in low speed regions. For the pressure driven flow system, we observed nonlinear distribution of pressure in the long channel and numerical results showed a good agreement with the experimental results.

  • PDF

Navier-Stokes 방정식을 이용한 천음속 익형의 설계최적화 연구 (Design Optimization of Transonic Airfoils Based on the Navier-Stokes Equation)

  • 이형민;조창열
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 1999년도 춘계 학술대회논문집
    • /
    • 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.

  • PDF

비대칭급확대채널의 층류유동 및 열전달 해석 (Analysis of Laminar Flow and Heat Transfer in Asymmetric, Sudden Expansion Channel)

  • 원승호;맹주성;손병진
    • 대한설비공학회지:설비저널
    • /
    • 제13권1호
    • /
    • pp.5-13
    • /
    • 1984
  • This analysis of numerical procedure is prediction of laminar flow and heat transfer at two dimension and steady flow in asymmetric sudden expansion channel. At former study, to analyse the flows with separation, the full Navier-Stokes equation is used, but there are many difficulties to analyse, and although significant progress has been made in the development of efficient computational methods for the Navier-Stokes equations, very large computation times are still required. In case of reward-facing flow, boundary-layer equation is used instead of full Navier-Stokes equation to analyse velocity fields, and result of this numerical analysis is good agreement with the given experimental study. In this case, since the computer time required for the boundary-layer calculation is an order of magnitude less than required for the solution of the full Navier-Stokes equation, this boundary-layer model provides a good approximate solution.

  • PDF

이산화된 Navier-Stokes 방정식의 영역분할법을 위한 병렬 예조건화 (Parallel Preconditioner for the Domain Decomposition Method of the Discretized Navier-Stokes Equation)

  • 최형권;유정열;강성우
    • 대한기계학회논문집B
    • /
    • 제27권6호
    • /
    • pp.753-765
    • /
    • 2003
  • A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

예조건화된 Navier-Stokes 방정식에서의 풍상차분법의 수치특성 (Numerical Characteristics of Upwind Schemes for Preconditioned Navier-Stokes Equations)

  • 길재흥;이두환;손덕영;최윤호;권장혁;이승수
    • 대한기계학회논문집B
    • /
    • 제27권8호
    • /
    • pp.1122-1133
    • /
    • 2003
  • Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.