• 한국가시화정보학회:학술대회논문집
• /
• 2006.12a
• /
• pp.149-152
• /
• 2006

This paper describes the flow visualization around the model of a commercial passenger airplane, Boeing 747-400, which cruises in the transonic speed. The geometry was realized through the reverse engineering based on the photogrammetry. The three-dimensional inviscid steady compressible governing equations are solved in the unstructured grid system under the cruise condition and in a finite volume method. The convective term is processed by the Crank-Nicholson scheme and first order upwind scheme is applied. The lift and drag forces in the wing with engines increase by 1.49% End 3.9%, respectively compared with the wing without engines.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.8 no.3
• /
• pp.53-60
• /
• 2004

In order to investigate the nozzle - rotor interactions and the effect of partial admission, the flows in supersonic turbine rotor cascades with a nozzle have been computed. Extensive computations of partial admission supersonic turbines provide the shock structures and flow patterns in the nozzle and rotor. The governing equations were discretized with Euler implicit method in time and 2nd-order upwind scheme of FVM in space. The $\kappa$-$\varepsilon$ turbulence model was utilized to describe the turbulent flow field. It is clearly shown that the nozzle flow is highly affected by the shocks or expansion waves propagated from the rotor leading edge. And the rotor flow is also affected by the shocks or wakes originated from the nozzle.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.115-126
• /
• 2002

A mesoscale meteorological model is applied to simulate turbulent airflow and eddy shedding over the Isle of Arran, SW Scotland, UK. Under conditions of NW flow, the mountain ridge of Kintyre, located upwind of Arran, induces gravity waves that also affect the airflow over the island. The possibility to nest domains allows description of the airflow over Arran with a very high resolution grid, while also including the effects of the surrounding mainland of Scotland, in particular of the mountain ridge of Kintyre. Initialised with a stably stratified NW flow, the mesoscale model simulates quasi-stationary gravity waves over the island induced by Kintyre. Embedded in the larger scale wave trains there is continuous development of small-scale transient eddies, created at the Arran hill tops, that move downstream through the stationary wave field. Although the transient eddies are more frequently simulated on the northern island where the terrain is more pronounced, they are also produced over Tighvein, a hill of 458 m on the southern island where measurements of surface pressure and 2 m meteorological variables have been recorded at intermittent intervals between 1996 and 2000. Comparison between early observations and simulations so far show qualitatively good agreement. Overall the computations demonstrate that turbulent flow can be modelled with a horizontal resolution of 70 m, and describe turbulent eddy structure on wavelength of only a few hundred metres.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• Journal of Mechanical Science and Technology
• /
• v.15 no.3
• /
• pp.366-373
• /
• 2001

A numerical analysis of shock wave/boundary layer interaction in transonic/supersonic axial flow compressor cascade has been performed by using a characteristics upwind Navier-Stokes method with various turbulence models. Two equation turbulence models were applied to transonic/supersonic flows over a NACA 0012 airfoil. The results are superion to those from an algebraic turbulence model. High order TVD schemes predicted shock wave/boundary layer interactions reasonably well. However, the prediction of SWBLI depends more on turbulence models than high order schemes. In a supersonic axial flow cascade at M=1.59 and exit/inlet static pressure ratio of 2.21, k-$\omega$ and Shear Stress Transport (SST) models were numerically stables. However, the k-$\omega$ model predicted thicker shock waves in the flow passage. Losses due to shock/shock and shock/boundary layer interactions in transonic/supersonic compressor flowfields can be higher losses than viscous losses due to flow separation and viscous dissipation.