• 한국자동차공학회논문집
• /
• 제18권3호
• /
• pp.60-68
• /
• 2010

In order to understand the flow and filtration characteristics in a wall-flow type DPF(Diesel Particulate Filter), 0-D, 1-D, and 3-D simulations are preformed. In this paper, three model are explained and validated with each other. Based on the comparisons with 1-D and 3-D results for the steady state solution, 3-D CFD analysis is preferable to 1-D for the prediction of wall velocity at the inlet and exit plane. Because PM loading process is transient state phenomena, the combination of full 3-D and time dependent simulation is crucial for the configuration of wall channels. New coupling technique, which is the connection between calculated permeability from 0-D lumped parameter model and UDF(User Defined Functions) of main solver, is proposed for the realisti

• 대한기계학회논문집B
• /
• 제31권1호
• /
• pp.68-75
• /
• 2007

A three-dimensional computation is conducted to simulate a three-dimensional rotating stall in a low speed axial compressor. It is generally known that a tip leakage flow has an important role on a stall inception. However, almost of researchers have taken no interest in a role of the hub-comer-stall on the rotating stall even though it is a common feature of the flow in an axial compressor operating near stall and it has a large effect on the flows and loss characteristics. Using a time-accurate unsteady simulation, it is found that the hub-comer-stall may be a trigger to collapse the axisymmetric flows under high loads. An asymmetric disturbance is initially originated in the hub-comer-stall because separations are naturally unstable flow phenomena. Then this disturbance is transferred to the tip leakage flows from the hub-comer-stall and grows to be stationary stall cells, which adheres to blade passage and rotate at the same speed as the rotor. When stationary stall cells reach a critical size, these cells then move along the blade row and become a short-length-scale rotating stall. The rotational speed of stall cells quickly comes down to 79 percent of rotor so they rotate in the opposite direction to the rotor blades in the rotating frame.

• 대한기계학회논문집
• /
• 제18권5호
• /
• pp.1218-1226
• /
• 1994

The shock tube is a useful device for investigating shock phenomena, spray combustion, unsteady gas dynamics, etc. Therefore, it is necessary to analyze exactly the flow phenomena in shock tube. In this study, the mechanics of its reflected shock zone has been investigated by using of the one-dimensional gas dynamic theory in order to estimate the transition from initial reflection of shock wave region. Calulation for four kinds of reflected shock tube temperature (i.e. (a) 1388 K (b) 1276 K (c) 1168 K (d) 1073 K) corresponding to the experimental conditions have been carried out sumarized as follows. (1) The qualitative tendency is almost the same as in that conditions in region of reflected wave region. (2) High temperature period (reflected shock wave temperature) $T_{5}$, exists 0-2.65 ms. (3) Transition period from temperature of reflection shock wave is far longer than the calculated one. This principally attributed to the fact that the contact surface is accelerated, also, due to the release of energy by viscoity effect. This apparatus can advance the ignition process of a spray in a ideal condition that involved neither atomization nor turbulent mixing process, where, using a shock tube, a column of droplets freely from atomizer was ignited behind a reflected shock.

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

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2007년도 추계학술대회
• /
• pp.116-118
• /
• 2007

Three-dimensional (3D) velocity field information of a laminar flow in a curved micro tube of circular cross-section has been measured using a digital micro holographic particle tracking velocimetry (HPTV). The temporal evolution of instantaneous velocity field of a water flow in a curved micro tube of $100\;{\mu}\;m$ and $300\;{\mu}\;m$ in inner diameter was obtained. The 3D mean velocity field distribution was obtained quantitatively by statistical-averaging of instantaneous velocity fields. At low Dean number (De), a secondary flow was not generated in the curved tube. With increasing Dean number, the secondary flow constituted of two large-scale counter-rotating vortices was formed due to enhanced centrifugal force. To reveal the flow characteristics of high Dean numbers, trajectories of fluid particles were evaluated experimentally from the 3D velocity fields data measured by the HPTV technique. The present experimental results, especially the 3D particle trajectories, would be helpful to design and to understand the mixing phenomena in 3D curved passages of various curved micro-tubes or micro-channels.

• Nuclear Engineering and Technology
• /
• 제55권5호
• /
• pp.1651-1664
• /
• 2023

It is very difficult to capture the multi-dimensional phenomena such as asymmetric flow and temperature distributions with the one-dimensional (1D) model, obviously, due to its inherent limitation. In order to overcome such a limitation of the 1D representation, many state-of-the-art system codes have equipped a three-dimensional (3D) component for multi-dimensional analysis capability. In this study, a standard multi-dimensional analysis model of APR1400 (Advanced Power Reactor 1400) has been developed using TRACE (TRAC/RELAP Advanced Computational Engine). The entire reactor pressure vessel (RPV) of APR1400 has been modeled using a single 3D component. The fuels in the reactor core have been described with detailed and coarse representations, respectively, to figure out the impact of the fuel description. Using both 3D RPV models, a comparative analysis has been performed postulating a double-ended guillotine break at a cold leg. Based on the results of comparative analysis, it is revealed that both models show no significant difference in general plant behavior and the model with coarse fuel model could be used for faster transient analysis without reactor kinetics coupling. The analysis indicates that the asymmetric temperature and flow distributions are captured during the transient, and such nonuniform distributions contribute to asymmetric quenching behaviors during blowdown and reflood phases. Such asymmetries are directly connected to the figure of merits in the LBLOCA analysis. Therefore, it is recommended to employ a multi-dimensional RPV model with a detailed fuel description for a realistic safety analysis with the consideration of the spatial configuration of the reactor core.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

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

In the present study, a micro holographic PTV (HPTV) system was used to experimentally investigate the structure of 3D flow within a curved micro-tube with varying Dean number. The employed HPTV system incorporated a high-speed digital camera to measure the temporal evolution of the 3D velocity fields of micro-scale fluid flows. With increasing Dean number, flow in the curved tube is transformed from a steady flow to a secondary flow with two counter-rotating vortices. In this study, to analyze the 3D flow characteristics in the curved section of tube at a high Dean number, the trajectories of fluid particles were obtained experimentally using the whole 3D velocity field data obtained by the micro HPTV technique. The mean velocity field distribution was then obtained by ensemble averaging the instantaneous velocity fields. These results would be helpful in the design of various passages within micro-scale devices or micro-chips and in understanding the mixing phenomena that occur in curved conduits along the trajectories of fluid particles.

• 한국항공우주학회지
• /
• 제37권9호
• /
• pp.879-888
• /
• 2009

초음속 주 유동내 수소의 수직분사에 의한 비정상 반응 유동장에 대한 3차원 수치해석이 DES 난류 모델과 상세 화학반응 모델을 이용하여 수행되었다. 난류 반응 유동의 물리적 현상을 이해하기 위하여 해석 및 실험 결과를 비교하였다. 계산에 의해 구해진 OH 분포는 실험의 OH-PLIF 결과를 잘 모사하고 있다. 반면, 점화 지연 시간은 계산과 실험 사이에 차이를 보이고 있으며, 이는 실험적 계측의 한계에 기인하는 것으로 생각된다. RANS 및 DES 계산 결과의 비교로부터 간헐 현상을 확인하였으며, 유선을 따른 온도 분포 및 중첩된 OH 질량 분율을 통해 시 공간적 간헐 현상을 정량적으로 측정하였다.