• Title/Summary/Keyword: unstructured finite volume model

Search Result 66, Processing Time 0.02 seconds

NUMERICAL METHOD FOR TWO-PHASE FLOW ANALYSIS USING SIMPLE-ALGORITHM ON AN UNSTRUCTURED MESH (비정렬격자 SIMPLE 알고리즘기반 이상유동 수치해석 기법)

  • Kim, Jong-Tae;Park, Ik-Kyu;Cho, Hyung-Kyu;Kim, Kyung-Doo;Jeong, Jae-Jun
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2008.10a
    • /
    • pp.71-78
    • /
    • 2008
  • For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

  • PDF

NUMERICAL METHOD FOR THE TWO-FLUID THREE-FIELD MODEL ON AN UNSTRUCTURED MESH (비정렬격자 2-유체 3-상 유동 해석 기법)

  • Kim, J.;Park, I.K.;Cho, H.K.;Yoon, H.Y.;Jeong, J.J.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2007.10a
    • /
    • pp.243-248
    • /
    • 2007
  • A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

  • PDF

Development of a 3-D Viscous Flow Solver Based on Unstructured Hybrid Meshes (비정렬 혼합 격자계 기반의 삼차원 점성 유동해석코드 개발)

  • Jung, Mun-Seung;Kwon, Oh-Joon
    • Journal of the Korean Society for Aeronautical & Space Sciences
    • /
    • v.35 no.8
    • /
    • pp.677-684
    • /
    • 2007
  • In the Present Study, a 3-D viscous flow solver, based on unstructured hybrid meshses containing tetrahedra, prisms and pyramids, has been developed. A finite-volume discretization scheme is used for solving the compressible Navier-Stokes equations. A cell-vertex median dual volume is used for spatial discretization. The one-equation Spalart-Allmaras turbulence model has been adopted to evaluate the eddy viscosity. Validation were made by computing laminar and turbulent flows around a 3-D wing for steady flows and turbulent flows around an oscillating 3-D wing in harmonic motion for unsteady flows.

Flow Visualization and Unstructured Grid Computation of Flow over a High-Speed Projectile (고속탄자 유동의 가시화 실험 및 비정렬격자 계산)

  • 이상길;최서원;강준구;임홍규;백영호;김두연;강호철
    • Transactions of the Korean Society of Automotive Engineers
    • /
    • v.6 no.2
    • /
    • pp.12-20
    • /
    • 1998
  • Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

  • PDF

A Numerical Study on Characteristics of Flood Wave Passing through Urban Areas (1) : Development and Verification of a Numerical Model (도시지역을 관통하는 홍수파의 특성에 관한 수치적 연구 (1) : 수치모형의 개발 및 검증)

  • Jeong, Woo-Chang;Lee, Jin-Woo;Cho, Yong-Sik
    • Journal of the Korean Society of Hazard Mitigation
    • /
    • v.9 no.6
    • /
    • pp.89-97
    • /
    • 2009
  • In this study, a two-dimensional unstructured finite volume model based on the shallow-water equations and well-balanced HLLC scheme is developed. The model is verified by applying to various one- and two-dimensional problems related to the analyses of dam-break wave. The predicted numerical results agree very well with available analytical solutions and laboratory measurements. The model provides slightly more accurate results compared with the existing models.

Assessment of Rotor Hover Performance Using a Node-based Flow Solver

  • Jung, Mun-Seung;Kwon, Oh-Joon;Kang, Hee-Jung
    • International Journal of Aeronautical and Space Sciences
    • /
    • v.8 no.2
    • /
    • pp.44-53
    • /
    • 2007
  • A three-dimensional viscous flow solver has been developed for the prediction of the aerodynamic performance of hovering helicopter rotor blades using unstructured hybrid meshes. The flow solver utilized a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart- Allmaras one-equation turbulence model. Calculations were performed at three operating conditions with varying tip Mach number and collective pitch setting for the Caradonna-Tung rotor in hover. Additional computations are made for the UH-60A rotor in hover. Reasonable agreements were obtained between the present results and the experiment in both blade loading and overall rotor performance. It was demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

CAVITATION FLOW ANALYSIS OF 2-D HYDROFOIL USING A HOMOGENEOUS MIXTURE MODEL ON UNSTRUCTURED MESHES (비정렬 격자계에서 균질혼합 모델을 이용한 2차원 수중익형 주위의 캐비테이션 유동 해석)

  • An, S.J.;Kwon, O.J.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2011.05a
    • /
    • pp.20-24
    • /
    • 2011
  • In this paper, numerical simulation of cavitation flow for modified NACA66 hydrofoil was made by using the multi-phase RANS equation based on pseudo-compressibility. The Homogeneous mixture model comprised of the mixture continuity, mixture momentum and liquid volume fraction equations was utilized. A vertex-centered finite-volume method was used in conjunction 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing The Spalart-Allmaras one equation model was employed for the closure of turbulence. Reasonable agreements were obtained between the calculation results and the experiment for pressure coefficients on the hydrofoil surface.

  • PDF

A Study on Simulation of Dam-break Wave Using Two-dimensional Finite Volume Model (2차원 유한체적모형을 이용한 댐 붕괴파 모의에 관한 연구)

  • Jeong, Woo-Chang;Park, Young-Jin
    • Journal of Korea Water Resources Association
    • /
    • v.44 no.3
    • /
    • pp.249-262
    • /
    • 2011
  • In this study, in order to reduce the numerical oscillation due to the unbalance between source and flux terms as the HLLC scheme is applied to the flow analysis on the irregular bed topography, a unstructured finite volume model based on the well-balanced HLLC scheme and the shallow water equations is developed and applied to problems of dam-break waves. The well-balanced HLLC scheme considers directly the gradient of bed topography as the flux terms is calculated. This scheme provides the good numerical balance between the source and flux terms in the case of the application to the steady-state transcritical flow. To verify the numerical model developed in this study, it is applied to three cases of hydraulic model experiments and a field case study of Mapasset dam failure (France). As a result of the verification, the predicted numerical results agree relatively well with available laboratory and field measurements. The model provides slightly more accurate results compared with the existing models.

CAVITATION FLOW SIMULATION FOR A 2-D HYDROFOIL USING A HOMOGENEOUS MIXTURE MODEL ON UNSTRUCTURED MESHES (비정렬 격자계에서 균질혼합 모델을 이용한 2차원 수중익형 주위의 캐비테이션 유동 해석)

  • Ahn, S.J.;Kwon, O.J.
    • Journal of computational fluids engineering
    • /
    • v.17 no.1
    • /
    • pp.94-100
    • /
    • 2012
  • In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

Numerical Simulation of MIT Flapping Foil Experiment : Unsteady Flow Characteristics (MIT 요동 익형의 수치해석 : 비정상 유동 특성)

  • Bae Sang Su;Kang Dong Jin;Kim Jae Won
    • 한국전산유체공학회:학술대회논문집
    • /
    • 1998.11a
    • /
    • pp.133-140
    • /
    • 1998
  • A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

  • PDF