• International Journal of Aeronautical and Space Sciences
• /
• 제13권3호
• /
• pp.307-316
• /
• 2012

Numerical simulations of 3D aircraft configurations are performed in order to understand the effects of turbulence models on the prediction of aircraft's aerodynamic characteristics. An in-house CFD code that solves 3D RANS equations and two-equation turbulence model equations are used. The code applies Roe's approximated Riemann solver and an AF-ADI scheme. Van Leer's MUSCL extrapolation with van Albada's limiter is also adopted. Various versions of Menter's $k-{\omega}$ SST turbulence models as well as Coakley's $q-{\omega}$ model are incorporated into the CFD code. Menter's $k-{\omega}$ SST models include the standard model, the 2003 model, the model incorporating the vorticity source term, and the model containing controlled decay. Turbulent flows over a wing are simulated in order to validate the turbulence models contained in the CFD code. The results from these simulations are then compared with computational results from the $3^{rd}$ AIAA CFD Drag Prediction Workshop. Numerical simulations of the DLR-F6 wing-body and wing-body-nacelle-pylon configurations are conducted and compared with computational results of the $2^{nd}$ AIAA CFD Drag Prediction Workshop. Aerodynamic characteristics as well as flow features are scrutinized with respect to the turbulence models. The results obtained from each simulation incorporating Menter's $k-{\omega}$ SST turbulence model variations are compared with one another.

• 대한기계학회논문집B
• /
• 제20권8호
• /
• pp.2561-2571
• /
• 1996

The flow in the turbomachinery is very unsteady due to the stator-rotor interaction. It has been indicated that the stator-rotor interaction has three distinct causes of unsteadiness: that is, the viscous vortex shedding, wake rotor interaction and potential stator-rotor interaction. In this paper, the mechanism of unsteady potential interaction and wake interaction in the stator-rotor stage flow is numerically investigated in two-dimensional view point. The numerical technique used is the upwind scheme of Van Leer's Flux Vector Splitting(FVS) and cubic spline interpolation is applied on zonal interface. Then, the flow field of a compressor stage composed of NACA 65410 is analyzed. Flow fields are found to be simulated reasonably by this method and the sensitivity due to back-pressure variation is more stronger than rotor-velocity variation.

• International Journal of Aeronautical and Space Sciences
• /
• 제12권2호
• /
• pp.149-155
• /
• 2011

The viscous flow in an axisymmetric nozzle was analyzed while accounting for the mesh sizes in both in the free stream and the boundary layer. The Navier-Stokes equations were resolved using the finite volume method in order to determine the supersonic flow parameters at the exit of the converging-diverging nozzle. The numerical technique in the aforementioned method uses the flux vector splitting of Van Leer. An adequate time stepping parameter, along with the Courant, Friedrich, Lewis coefficient and mesh size level, was selected to ensure numerical convergence. The boundary layer thickness significantly affected the viscous flow parameters at the exit of the nozzle. The best solution was obtained using a very fine grid, especially near the wall at which a strong variation of velocity, temperature and shear stress was observed. This study confirmed that the boundary layer thickness can be obtained only if the size of the mesh is lower than a certain value. The nozzles are used at the exit of the shock tube in order to obtain supersonic flows for various tests. They also used in propulsion to obtain the thrust necessary to the displacement of the vehicles.

• 한국전산유체공학회지
• /
• 제8권4호
• /
• pp.16-25
• /
• 2003

The design method for transonic turbine blades has been developed based on Wavier-Stokes equations. The present computing process is done on the four separate steps, i.e., determination of the blade profile, generation of the computational grids, cascade flow simulation and analysis of the computed results in the sense of the aerodynamic performance. The blade shapes are designed using the cubic polynomials under the control of the design parameters. Numerical methods for the flow equations are based on Van-Leer's FVS with an upwind TVD scheme on the finite volume. In the present study, numerical simulation has been done to investigate the effects of the design parameters on the aerodynamic peformance of the axial-flow turbine blades. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to four parameters and compared with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

• 한국전산유체공학회지
• /
• 제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.

• 대한기계학회논문집B
• /
• 제23권11호
• /
• pp.1371-1378
• /
• 1999

Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

• Nuclear Engineering and Technology
• /
• 제42권2호
• /
• pp.211-218
• /
• 2010

An analysis model on debris transport in the containment floor of pressurized water reactors is developed in which the flow field is calculated by Eulerian conservation equations of mass and momentum and the debris particles are traced by Lagrange equations of motion using the pre-determined flow field data. For the flow field calculation, two-dimensional Shallow Water Equations derived from Navier Stokes equations are solved using the Finite Volume Method, and the Harten-Lax-van Leer scheme is used for accuracy to capture the dry-to-wet interface. For the debris tracing, a simplified two-dimensional Lagrangian particle tracking model including drag force is developed. Advanced schemes to find the positions of particles over the containment floor and to determine the position of particles reflected from the solid wall are implemented. The present model is applied to calculate the transport fraction to the Hold-up Volume Tank in Advanced Power Reactors 1400. By the present model, the debris transport fraction is predicted, and the effect of particle density and particle size on transport is investigated.

• Nuclear Engineering and Technology
• /
• 제41권8호
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• 한국전산유체공학회지
• /
• 제5권2호
• /
• pp.38-46
• /
• 2000

It has been indicated that the rotor/stator interaction has distinct causes of unsteadiness, such as the viscous vortex shedding, wake/stator interaction and potential rotor/stator interaction. In this paper, the mechanism of unsteady potential interaction in one stage axial compressor is numerically investigated for blade row ratio 1:1 and 2:3 at design point and for blade row ratio 2:3 at off-design point in two-dimensional view point. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting(FVS) and Cubic spline interpolation is applied on zonal interface. In this study the flow unsteadiness due to potential interaction are found to be larger in blade row ratio 2:3 than in 1:1. The total pressure rise in blade row ratio 2:3 is closer to the real value in design point than that in 1:1. The change of unsteady pressure amplitude according to the variation of stator exit pressure is very small.