• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.9
• /
• pp.907-915
• /
• 2008

The computational fluid dynamic analysis has been conducted for the thermo-chemical flow field in an arcjet thruster with mono-propellant Hydrazine (N2H4) as a working fluid. The Reynolds Averaged Navier-Stokes (RANS) equations are modified to analyze compressible flows with the thermal radiation and electric field. the Maxwell equation, which is loosely coupled with the fluid dynamic equations through the Ohm heating and Lorentz forces, is adopted to analyze the electric field induced by the electric arc. The chemical reactions of Hydrazine were assumed to be infinitely fast due to the high temperature field inside the arcjet thruster. The chemical and the thermal radiation models for the nitrogen-hydrogen mixture and optical thick media respectively, were incorporated with the fluid dynamic equations. The results show that performance indices of the arcjet thruster with 1kW arc heating are improved by amount of 180% in thrust and 200% in specific impulse more than frozen flow. In addition thermo-physical process inside the arcjet thruster is understood from the flow field results.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.60-60
• /
• 2021

Investigating Backward-Facing Step(BFS) flow is important in that it is a representative case for separation flows in various engineering flow systems. There have been a wide range of experimental, theoretical, and numerical studies to investigate the flow characteristics over BFS, such as flow separation, reattachment length and recirculation zone. However, most of such previous studies were concentrated only on the perpendicular step angle. In this study, several numerical investigations on the flow pattern over BFS with various step angles (10° ~ 90°) and expansion ratios (1.48, 2 and 3.27) under different Reynolds numbers (5000 ~ 64000) were carried out, mainly focused on the reattachment length. The numerical simulations were performed using an open source 3D CFD software, OpenFOAM, in which the velocity profiles and turbulence intensities are calculated by RANS (Reynolds Averaged Navier-Stokes equation) and 3D LES (Large Eddy Simulation) turbulence models. Overall, it shows a good agreement between simulations and the experimental data by Ruck and Makiola (1993). In comparison with the results obtained from RANS and 3D LES, it was shown that 3D LES model can capture much better and more details on the velocity profiles, turbulence intensities, and reattachment length behind the step for relatively low Reynolds number(Re < 11000) cases. However, the simulation results by both of RANS and 3D LES showed very good agreement with the experimental data for the high Reynolds number cases(Re > 11000). For Re > 11000, the reattachment length is no longer dependent on the Reynolds number, and it tends to be nearly constant for the step angles larger than 30°.) Based on the calibrated and validated numerical simulations, several additional numerical simulations were also conducted with higher Reynolds number and another expansion ratio which were not considered in the experiments by Ruck and Makiola (1993).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.627-634
• /
• 2008

In a numerical simulation of open channel turbulent flows using RANS (Reynolds averaged Navier-Stokes) equations model equipped with VOF (Volume of Fluid) scheme, the determination of wall roughness for wall function was studied. The roughness constant, based on the law-of-the-wall for flow on rough walls, obtained by experimental works for pipe flows is employed in general wall functions. However, this constant of wall function is the function of Froude number in open channel flows. Thus, the wall roughness should be determined by taking into account the effect of Froude number. In addition, the wall roughness should be corresponding to Manning's roughness coefficient widely used for open channels. In this study, the relation between wall roughness height as an input condition and Manning's roughness coefficient was investigated, and an equation for effective wall roughness height considering the characteristics of numerical models was proposed as a function of Manning's roughness coefficient.

• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.181-188
• /
• 2005

Large eddy simulation(LES) methodology used to model the isothermal swirling flows in a dump combustor and the turbulent premixed flame in a model gas turbine combustor. The LES solver was implemented on parallel computer consisting 16 processors. In isothermal flow simulation, the results was compared with that of ${\kappa}-{\varepsilon}$ model as well as experimental data, in order to verify the capability of LES code. To model the turbulent premixed flame in a gas turbine, the G-equation flamelet model was used. The results showd that LES and RANS well predicted the mean velocity field of a non-swirling flow. However, in swirling flow, LES showed a better performance in predicting the mean axial and azimuthal velocities, and the central recirculation zone than those of RANS. In a model gas turbine combustor, the operation condition of high pressure and temperature induced the different phenomena, such as flame length and flow-field information, comparing with the condition of ambient pressure and temperature. Finally, it was identified that the flame and heat release oscillations are related to the vortex shedding generated by swirl flow and pressure wave propagation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.18-25
• /
• 2013

Assessment and validation of RANS turbulence models are conducted for the optimal analysis of supersonic converging-diverging nozzle through the comparison between computational results and experimental data. One/two equation turbulence closures such as Spalart-Allmaras, RNG k-${\varepsilon}$, and k-${\omega}$ SST are employed to simulate the two-dimensional nozzle flow. Computational results with the turbulence models mentioned fairly well predict shock structure of the nozzle-inside and pressure distribution along the wall. Especially, SST model among the employed ones shows the best agreement to experimental results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.201-209
• /
• 2011

An unsteady RANS analysis study of the 3-D flow around the NREL Phase VI horizontal axis wind turbine(HAWT) was performed using sliding mesh approach. Two different analysis models such as rotor-only and rotor with tower/nacelle were constructed to investigate the blade/tower interaction. Analysis results for the rotor with tower/nacelle were compared with the corresponding NREL's experimental data which produced fairly good validation of the present CFD model. Comparison of flows around those two models also clearly showed the blade/tower interaction even it was small for upwind configuration. Other visualization results and integrated aerodynamic loads including torque of the blade demonstrated the effective unsteady flow simulation capability of the present CFD model.

• International Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.67-72
• /
• 2015

The main concern of this study is to integrate the EARSM into an industrial RANS solver in conjunction with the $k-{\omega}$ model, as proposed by Hellsten (EARSMKO2005). In order to improve the robustness, particular limiters are introduced to turbulent conservative variables, and a suitable full-approximation storage (FAS) multi-grid (MG) strategy is designed to incorporate turbulence model equations. The present limiters and MG strategy improve both robustness and efficiency significantly but without degenerating accuracy. Two discretization approachs for velocity gradient on cell interfaces are implemented and compared with each other. Numerical results of a three-dimensional supersonic square duct flow show that the proper discretization of velocity gradient improves the accuracy essentially. To assess the capability of the resulting EARSM $k-{\omega}$ model to predict complex engineering flow, the case of Common Research Model (CRM, Wing-Body) is performed. All the numerical results demonstrate that the resulting model performs well and is comparable to the standard two-equation models such as SST $k-{\omega}$ model in terms of computational effort, thus it is suitable for industrial applications.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.331-334
• /
• 2011

The performance of ramjet engine is closely associated with a supersonic intake. In this study, experiments and computational simulations were conducted to observe the internal flow characteristics of the supersonic intake. The supersonic intake which have self-starting characteristics was designed and manufactured. The flow characteristics was analyzed from the experimental results using the supersonic wind tunnel testing and computational results using RANS equation and Menter's SST turbulence model. The detailed visualization results were suggested for the pseudo-shock wave of stable operations and for the inlet buzz phenomenon of unstable operations.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.25 no.6
• /
• pp.36-44
• /
• 2021

Numerical analysis was performed on the shock wave-boundary layer interaction generated in the internal flow path of the curved interstage of the scramjet engine flight test vehicle. For numerical analysis, the turbulence model k-ω SST was used in the compressibility Raynolds Averaged Navier Stokes(RANS) equation. Representatively, the separation bubbles on the upper wall of the nozzle, the interaction between the concave shock wave and the boundary layer, and the shock wave-shock wave interaction at the edge were captured. The analysis result visualizes the shock wave-boundary layer interaction of the curved internal flow path to enhance understanding and suggest design considerations.