• 대한기계학회논문집B
• /
• 제39권5호
• /
• pp.415-423
• /
• 2015

본 연구에서는 2차원 비압축성 유체 유동을 해석하기 위하여 격자 세밀화 모델을 적용한 격자볼츠만 방법(LBM)을 수치적으로 연구하였다. 일반적으로 유동해석에서 충분한 정확도를 얻기 위해서는 관심 있는 영역에서 격자가 세밀하게 구성되어야 한다. 그래서 본 연구에서는 유동함수-와도 공식을 적용한 LBM에 격자 세밀화 모델을 적용하여 유동해석을 수행하였다. 공동형상 유동에서의 기존의 신뢰성 있는 유동장 결과와의 비교를 통해 본 연구의 격자 세밀화 모델을 적용한 격자볼츠만 방법의 신뢰성과 유용성을 검토하였다.

• Wind and Structures
• /
• 제14권5호
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.

• Ocean and Polar Research
• /
• 제44권4호
• /
• pp.269-285
• /
• 2022

The interaction between ocean and ice shelf is a critical physical process in relation to water mass transformations and ice shelf melting/freezing at the ocean-ice interface. However, it remains challenging to thoroughly understand the process due to a lack of observational data with respect to ice shelf cavities. This is the first study to simulate the variability and circulation of water mass both overlying the continental shelf and underneath an ice shelf and an ice tongue in the Terra Nova Bay (TNB), East Antarctica. To explore the properties of water mass and circulation patterns in the TNB and the corresponding effects on sub ice shelf basal melting, we explicitly incorporate the dynamic-thermodynamic processes acting on the ice shelf in the Regional Ocean Modeling System. The simulated water mass formation and circulation in the TNB region agree well with previous studies. The model results show that the TNB circulation is dominated by the geostrophic currents driven by lateral density gradients induced by the releasing of brine or freshwater at the polynya of the TNB. Meanwhile, the circulation dynamics in the cavity under the Nansen Ice shelf (NIS) are different from those in the TNB. The gravity-driven bottom current induced by High Salinity Shelf Water (HSSW) formed at the TNB polynya flows towards the grounding line, and the buoyance-driven flow associated with glacial meltwater generated by the HSSW emerges from the cavity along the ice base. Both current systems compose the thermohaline overturning circulation in the NIS cavity. This study estimates the NIS basal melting rate to be 0.98 m/a, which is comparable to the previously observed melt rate. However, the melting rate shows a significant variation in space and time.

• 설비공학논문집
• /
• 제6권1호
• /
• pp.20-28
• /
• 1994

A numerical study has been performed on enhancement of heat transfer in a forced convection of the modified driven-cavity flow which was previously found by the author to give a regular or chaotic stirring depending on the parameter value. It is found that for the present case wherein heat is transmitted between fluid and the surrounding walls, the chaotic stirring enhances the heat transfer at high Peclet numbers. The optimal condition of the flow modulation for the best heat transfer can be predicted by purely investigating the hydrodynamic facet, i.e. the stirring effect.

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

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.169-174
• /
• 1995

A finite element scheme using the concept of PISO method has been developed to solve the viscous flow problems in all speed range. In this study, new pressure equation is proposed such that both the hyperbolic term related with the density variations and elliptic term reflecting the incompressibility constraint are included. Present method has been applied to incompressible flow in two-dimensional driven cavity(Re=100, 400 and 1,000), and its computed results are compared with other's. Also, Carter plate problem(M=3 and Re=1,000) is computed and the comparison is made with Carter's results. Finally, we simulate a shock-boundary layer interaction problem(M=2 and $Re=2.96{\times}10^5$) to illustrate the shock capturing capability of the present solution algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.61-62
• /
• 2003

In this paper, the new higher order wall boundary conditions are proposed for solving the incompressible flows. The square driven cavity flows are simulated by using the variable order method of lines with the present wall boundary conditions. The variable order method of lines is constructed by the spatial discretization, i.e., the variable order proper convective scheme for convective terms and the modified differential quadrature method for diffusive terms, and time integration. The 2nd, 4th, and 6th order solutions are presented and these results show this higher order boundary conditions are very promising for the incompressible flow simulations.

• 한국자동차공학회논문집
• /
• 제11권2호
• /
• pp.157-163
• /
• 2003

This work was undertaken for the numerical analysis of defrosting phenomena of automobile windshield. To analyze the defrost, the flow and temperature field of cabin interior, heat transfer through the windshield glass, and phase change of the frost should be analyzed simultaneously. The flow field was obtained by solving the 3-D unsteady Navier-Stokes equation and the temperature field was computed by energy equation. The phase-change process of Stefan problem was solved by enthalpy method. For code validation, the temperature field of the driven cavity was calculated. The result of calculation shows a good agreement with the other numerical results. Then, the present code was applied to the defrosting analysis of a real automobile and, also, a good agreement with experiment was obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.86-94
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation or condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure distribution that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new reconstruction method to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function, a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the CUPID code.