• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.5
• /
• pp.652-659
• /
• 2000

The purpose of this study is to present the most effective numerical scheme to calculate the unsteady flows. In order to calculate the flow quantities of flow past a circular cylinder, Three-time level and five convective schemes are applied to unsteady and convective terms, respectively. The values obtained are compared with those from the existing experimental and numerical studies. At Reynolds numbers up to 160, time intervals can be expanded 10 times of Implicit Euler scheme using Three-time level method, and it is found that QUICK and CUI schemes work much stable than others even if less grid density conditions. The combination of Three-time level and QUICK scheme gives high resolutions for laminar unsteady problems with PC level.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1179-1194
• /
• 1992

QUICKER scheme has several attractive properties. However, under highly convective conditions, it produces overshoots and possibly some oscillations on each side of steps in the dependent variable when the flow is convected at an angle oblique to the grid line. Fortunately, it is possible to modify the QUICKER scheme using non-linear and linear functional relationship. Details of the development of polynomial upwinding scheme are given in this paper, where it is seen that this non-linear scheme has also third order accuracy. This polynomial upwinding scheme is used as the basis for the SHARPER and SMARTER schemes. Another revised scheme was developed by partial modification of QUICKER scheme using CDS and UPWIND schemes(QUICKUP). These revised schemes are tested at the well known bench mark flows, Two-Dimensional Pure Convection Flows in Oblique-Step, Lid Driven Cavity Flows and Buoyancy Driven Cavity Flows. For pure convection oblique step flow test problem, QUICKUP, SMARTER and SHARPER schemes remain absolutely monotonic without overshoot and oscillation. QUICKUP scheme is more accurate than any other scheme in their relative accuracy. In high Reynolds number Lid Driven Cavity Flow, SMARTER and SHARPER schemes retain lower computational cost than QUICKER and QUICKUP schemes, but computed velocity values in the revised schemes produced less predicted values than QUICKER scheme which is strongly effected by overshoot and undershoot values. Also, in Buoyancy Driven Cavity Flow, SMARTER, SHARPER and QUICKUP schemes give acceptable results.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• Atmosphere
• /
• v.21 no.4
• /
• pp.391-404
• /
• 2011

This study investigates the impact of cumulus parameterization scheme (CPS) with different horizontal grid sizes on the simulation of the local heavy rainfall case over the Korean Peninsula. The Weather Research and Forecasting (WRF)-based real-time forecast system of the Joint Center for High-impact Weather and Climate Research (JHWC) is used. Three CPSs are used for sensitivity experiments: the BMJ (Betts-Miller-Janjic), GD (Grell-Devenyi ensemble), and KF (Kain-Fritsch) CPSs. The heavy rainfall case selected in this study is characterized by low-level jet and low-level transport of warm and moist air. In 27-km simulations (DM1), simulated precipitation is overestimated in the experiment with BMJ scheme, and it is underestimated with GD scheme. The experiment with KF scheme shows well-developed precipitation cells in the southern and the central region of the Korean Peninsula, which are similar to the observations. All schemes show wet bias and cold bias in the lower troposphere. The simulated rainfall in 27-km horizontal resolution has influence on rainfall forecast in 9-km horizontal resolution, so the statements on 27-km horizontal resolution can be applied to 9-km horizontal resolution. In the sensitivity experiments of CPS for DM3 (3-km resolution), the experiment with BMJ scheme shows better heavy rainfall forecast than the other experiments. The experiments with CPS in 3-km horizontal resolution improve rainfall forecasts compared to the experiments without CPS, especially in rainfall distribution. The experiments with CPS show lower LCL(Lifted Condensation Level) than those without CPS at the maximum rainfall point, and weaker vertical velocity is simulated in the experiments with CPS compared to the experiments without CPS. It means that CPS suppresses convective instability and influences mainly convective rainfall. Consequently, heavy rainfall simulation with BMJ CPS is better than the other CPSs, and even in 3-km horizontal resolution, CPS should be applied to control convective instability. This conclusion can be generalized by conducting more experiments for a variety of cases over the Korean Peninsula.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.4
• /
• pp.147-157
• /
• 1987

Five finite difference schemes (explicit, Bresler, implicit, upstream, and Chaudhari scheme) for the convective dispersion model are analyzed numerically to investigate their characteristics and applicabilities. Camparative study results show that the conditionally stable Chaudhari scheme has the smallest numerical dispersion and that the unconditionally stable Bresler scheme has the overshooting in regions of oscillation. Explicit scheme is the most accurate for a dispersion-dominated flow whereas Chaudhari scheme is for a convection-dominated flow. The computation time (CPU) is the shortest for the explicit or Chaudhari scheme with the same order of magnitude and is always the longest for the Bresler scheme.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.57-67
• /
• 1996

The formulation for hybrid-QUICK scheme of convective transport terms in finite-volume calculation procedure is presented. Source terms are modified to apply the hybrid-QUICK scheme. Test calculations are performed for wall-driven cavity flow at Re=$10_2$, $10_3$, and $10_4$. These include the evaluation of boundary conditions approximated by third-order finite difference scheme. The stable and converged solutions are obtained without unsteady terms in the momentum equations. The results using hybrid-QUICK scheme show no difference with those using hybrid scheme at low Re ($=10_2$) and are better at higher Re ($10_3$, and $10_4$).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.5
• /
• pp.1661-1670
• /
• 1996

Incompressible flow over a backward-facing step is computed by low Reynolds number turbulence models in order to compare with direct simulation results. In this study, selected low Reynolds number 1st and 2nd (Algebraic Stress Model : ASM) moment closure turbulence models are adopted and compared with each other. Each turbulence model predicts different flow characteristics, different re-attachment point, velocity profiles and Reynolds stress distribution etc. Results by .kappa.-.epsilon. turbulence models indicate that predicted re-attachment lengths are shorter than those by standard model. Turbulent intensity and eddy viscosity by low Reynolds number .kappa.-.epsilon. models are still greater than DNS results. The results by algebraic stress model (ASM) are more reasonable than those by .kappa.-.epsilon. models. The convective scheme is QUICK (Quadratic Upstream Interpolation for Convective Kinematics) and SIMPLE algorithm is adopted. Reynolds number based on step height and inlet free stream velocity is 5100.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.49-56
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Navier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. The endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency. The effects of both turbulence model and convective differencing scheme on the prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence model on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.1-6
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Wavier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. the endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency The effects of both turbulence model and convective differencing scheme on the Prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence modei on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.