• 한국전산유체공학회지
• /
• 제7권4호
• /
• pp.35-41
• /
• 2002

The projectile afterbodies for zero-lift drag reduction has been analyzed using the Navier-Stokes equations with the κ-εturbidence model. The numerical method of a second order upwind scheme has been used on an unstructured adaptive grid system. Base drag reduction methods that have been found effective on axisymmetric bodies are boattailing, base bleed, base combustion, locked vortex afterbodies and multistep afterbodies. In this paper, turbulence flow and pressure charateristics have been studied for geometries of multistep afterbodies. The important geometrical and flow parameters relevant to the design of such afterbodies have been identified by step number, length and height. The flow over multistep aftoerbodies or base have many kinds of compressible flow characteristics including expansion waves at the trailing edge, recompression waves, separation and recirculating flow in the base region, shear flow and wake flow. The numerical results have been compared and analyzed with the experimental data. The flow characteristics have been clearly shown.

• 한국전산유체공학회지
• /
• 제10권4호통권31호
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.

• 대기
• /
• 제25권3호
• /
• pp.473-482
• /
• 2015

In this study, street canyons with a higher downwind building (so called, step-up street canyons) are considered for understanding characteristics of flow and reactive pollutants' dispersion as a basic step to understand the characteristics in wider urban areas. This study used a CFD_NIMR_SNU coupled to a chemistry module just including simple $NO_X-O_3$ photochemical reactions. First, flow characteristics are analyzed in step-up street canyons with four aspect ratios (0.33, 0.47, 0.6, 0.73) defined as ratios of upwind building heights to downwind building height. The CFD_NIMR_SNU reproduced very well the main features (that is, vortices in the street canyons) which appeared in the wind-tunnel experiment. Wind speed within the street canyons became weak as the aspect ratio increased, because volume of flow incoming over the upwind building decreased. For each step-up street canyon, chemistry transport model was integrated up to 3600 s with the time step of 0.5 s. The distribution patterns of $NO_X$ and $O_3$ were largely dependent on the mean flow patterns, however, $NO_X$ and $O_3$ concentrations were partly affected by photochemical reactions. $O_3$ concentration near the upwind lower region of the street canyons was much lower than background concentration, because there was much reduction in $O_3$ concentration due to NO titration there. Total amount of $NO_X$ in the street canyons increased with the aspect ratio, resulting from the decrease of mean wind intensity.

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

A numerical study of axisymmetric rocket main nozzle flow has been accomplished. The CSCM upwind flux difference splitting method with an iterative time marching scheme having second order accuracy in time and space has been used to simulate unsteady flow characteristics in an axisymmetric rocket main nozzle. Though the pressure vary at nozzle inlet with the lapse of time, Mach No. and the density were not changed significontly compared with the temperature. Specific heat ratio $\gamma$=1.134 predicted higher temperature at nozzle throat and exit and nondimensional thrust coefficients at exit than specific heat ratio $\gamma$=1.4 did.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2001년도 추계학술대회 논문집
• /
• pp.202-206
• /
• 2001

A three dimensional incompressible Navier-Stokes code based on the third derivative upwind is employed to simulate the flow around the underwater vehicle advancing on the calm water. Computations are carried out in the range of Froude numbers 0.4 to 0.7. The wave resistance, lift, moment and the pressure distribution on the body are calculated. Computations are performed in a rectangular grid system based on the Marker & Cell method. For validation, computation results are compared with existing experimental results.

• 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.

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

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

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

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

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

A two-dimensional/axisymmetric CSCM upwind flux difference splitting Wavier-Stokes method has been developed to study the finite rate chemically react-ing invisicd and viscous hypersonic flows over blunt-body. A upwind method was chosen due to its robustness in capturing the strong bow shock waves. For the nonequilibrium chemically reacting air, NS-I species conservation equations were strongly coupled with flowfield equations through convection and species production terms. The nonequilibrium wall pressure and heat transfer rate distributions along the vehicle were compared with those from equilibrium and perfect gas calculations. The nonequilibrium species distribution shows the reduced concentrations of O and N species when compared with equilibrium species distribution. The solutions resolved strong bow shock waves md heat transfer rate very accurately when compared with central difference schemes.

• 대한기계학회논문집B
• /
• 제27권8호
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.