• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.190-193
• /
• 2010

The turbulent mixing and combustion of a shear coaxial injector under supercritical pressures have been theoretically/numerically investigated. Turbulent numerical model is based on large eddy simulation with real-fluid transport and thermodynamics over the entire pressure; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The results are compared with previous researcher's. The large-scale vortices shedding from the outer rim into the recirculation region to react with gaseous oxygen was investigated.

• Nuclear Engineering and Technology
• /
• v.27 no.2
• /
• pp.189-202
• /
• 1995

An improved method is presented for the prediction of heat transfer coefficients in turbulent fall-ing liquid films with or without interfacial shear for both heating or condensation. A modified Mudawwar and El-Masri's semi-empirical turbulence model, particularly to extend its use for the turbulent falling film with high interfacial shear, is used to replace the eddy viscosity model incorporated in the unified approach unposed by Yih and Liu. The liquid film thickness and asymptotic heat transfer coefficients against the film Reynolds number for wide range of interfacial shear predicted by both present and existing methods are compared with experimental data. The results show that in general, predictions of the modified model agee more closely with experimental data than that of existing models.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.48-56
• /
• 2002

The flow characteristics around KRISO 3600TEU container ship model have been experimentally investigated in a circulating water channel. The instantaneous velocity vectors were measured using 2-frame PIV measurement system. The mean velocity fields and turbulent statistics including turbulent kinetic energy and vorticity were obtained by ensemble-averaging 400 instantaneous velocity fields. The free stream velocity was fixed at 0.6m/s and the corresponding Reynolds number was $9{\times}10^5$. The test sections were divided into two regions, three transverse sections of the wake region(Station -0.5767, -1, -3) and five longitudinal sections of the wake((Z/(B/2)=0, 0.1, 0.2, 0.4, 0.6). In the wake region, large-scale longitudinal vortices of nearly same strength are symmetric with respect to the wake centerline and a relatively weak secondary vortex is formed near the waterline. With going downstream, the strength of longitudinal vortex is decreased and the wake region expands.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.484-491
• /
• 2007

Numerical experiments using a numerical wave tank have been performed to verier the nonlinear interaction between monochromatic waves and a submerged horizontal plate. As a model for numerical wave tank, we used a higher-order Boundary Element Method(BEM) based on fully nonlinear potential flow theory and CADMAS-SURF for solving Navier Stokes equations and exact free surface conditions. Both nonlinear models are able to predict the higher harmonic generation in the shallow water region over a submerged horizontal plate. CADMAS-SURF, which involves the viscous effect, can evaluate the higher harmonic generation by flow separation and vortices at the each ends of plate. The comparison of reflection and transmission coefficients with experimental results(Patarapanich and Cheong, 1989) at different lengths and submergence depths of a horizontal plate are presented with a good agreement. It is found that the transfer of energy from the incident fundamental waves to higher harmonics becomes larger as the submergence depth ratio decreases and the length ratio increases.