• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

• Journal of Korea Water Resources Association
• /
• v.38 no.2
• /
• pp.155-166
• /
• 2005

Upwind schemes are very well adapted to the discontinuous flow and have become popular for applications Involving dam break flow, transcritical Slow, etc. However, upwind schemes have been applied mainly to the idealized problems not to the natural channels with irregular geometry so far because of the error due to source terms. In this paper, the new type of upwind discretization of source terms, which uses the normalized Jacobian to discretize the source terms, is proposed. As results of tests to flows with source terms by the upwind models, the method proposed in this paper is proved as efficient and accurate. This generalized method for differencing source terms is simple and might beapplicable to diverse type of flux upwind discretization scheme in finite difference method.

• Journal of applied mathematics & informatics
• /
• v.41 no.1
• /
• pp.181-192
• /
• 2023

This work deals with the numerical modeling of dam-break flow over erodible bed. The mathematical model consists of the shallow water equations, the transport diffusion and the bed morphology change equations. The system is solved by central upwind scheme. The obtained results of the resolution of dam-beak problem is presented in order to show the performance of the numerical scheme. Also a comparison of central upwind and Roe schemes is presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.224-232
• /
• 2000

Four turbulent $k-{\varepsilon}$ models (i.e., standard model, modified models with streamline curvature modification and/or preferential dissipation modification) are applied in order to analyze the turbulent flow of wall-attaching offset jet. For numerical convergence, this paper develops a method of slowly increasing the convective effect induced by skew-velocity in skew-upwind scheme (hereafter called Partial Skewupwind Scheme). Even though the method was simple, it was efficient in view of convergent speed, computer memory storage, programming, etc. The numerical results of all models show good prediction in first order calculations (i.e., reattachment length, mean velocity, pressure), while they show some deviations in ·second order (i.e., kinetic energy and its dissipation rate). Like the previous results obtained by upwind scheme, the streamline curvature modification results in better prediction, while the preferential dissipation modification does not.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.2
• /
• pp.119.2-119.2
• /
• 2014

Many astrophysical phenomena involve processes of magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD). A number of numerical schemes have been developed to solve the equations of ideal MHD and RMHD. Recent codes are based on upwind schemes which solve hyperbolic systems of equations following the characteristics of the systems. Upwind schemes stand out by their robustness, clarity of the underlying physical model, and ability of achieving high resolution. We present MHD and RMHD codes based on the total variation diminishing (TVD) and weighted essentially non-oscillatory (WENO) schemes, which are second and higher order accurate extensions of upwind schemes. We demonstrate the ability and limitation of codes based on upwind schemes through a series of tests.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.865-875
• /
• 1988

A two-dimensional stepped wall jet was numerically investigated by applying three different models : One is the standard k-.epsilon. and the other is the modified k-.epsilon. model which takes account of the streamline curvature effect by modifying the Reynolds shear stress and a source term in the dissipation equation, and a third is curvature dependent third-order correlation model. In order to test the influences of the numerical result, both the upwind scheme and the skew-upwind scheme were sued for the computations. By comparing the numerical results with available experiments, it was found that the modified k-.epsilon. model gives best overall prediction accuracy only when the numerical diffusion is eliminated by using the skew-upwind scheme. The numerical scheme was found to have more pronounced effect on the accuracy of the turbulence computation than the turbulence models.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.3
• /
• pp.1-9
• /
• 2004

This study concentrates on the upwind schemes for convergence acceleration of the multigrid method for the Navier-Stokes equations. Comparative study of the upwind schemes in the Fourier space has been performed to identify why the second-order upwind scheme with enlarged stencil can be preconditioned better than the classical second-order upwind scheme. The full-coarsening multigrid method with implicit preconditioned multistage scheme has been implemented for verification of analysis. Numerical simulations on the inviscid and turbulent flows with the Spalart-Allmaras turbulent model have been performed. The results showed consistent trend with the analysis.

• 한국해양학회지
• /
• v.21 no.4
• /
• pp.203-210
• /
• 1986

A linear parallel transport model is formulated and applied to an idealized Yellow Sea, With this simple analytical model, the hither-to suspected upwind flow phenomena in the southern Yellow Sea can be reasonably explained. In deep waters where the local depth exceeds a critical depth (Hc=53m in the present model sea), pressure gradient force dominates over wind stress and contributes to an upwind flow. The estimated upwind flow velocity increases with wind speed and a maximum upwind flow occurs along the axis of the Yellow Sea embayment. For the typical south wind of 5-10 knots in summer, the upwind (southward) flow velocity along the axis of the Yellow Sea is estimated to be 1-5cm s$\^$-1/. While, for the typical north wind of 10-15 knots in winter, the upwind (northward) flow velocity is 5-12cm s$\^$-1/. These velocity ranges can be served as rough estimates for the intrusion velocity of the Yellow Sea Bottom Cold Water in summer and the Yellow Sea Warm Current in winter, respectively.

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

• Asian Journal of Atmospheric Environment
• /
• v.4 no.2
• /
• pp.89-96
• /
• 2010

To identify the domains that have the greatest impacts on air quality at the surface, both the upwind and downwind areas of Seoul were determined by season using refined wind fields. Four consecutive days were selected as the study period typical of each season. The mesoscale meteorology of the study period was reproduced by using the MM5 prognostic meteorological model (PSU/NCAR Mesoscale Model) with horizontally nested grids. The gridded meteorological field, which was used on the study area of $242\;km{\times}226\;km$ with grid spacing of 2 km, was generated by using the CALMET diagnostic meteorological model. Upwind and downwind areas of Seoul were determined by calculating 24-hour backward and forward air parcel trajectories, respectively, with u, v, and w velocity vectors. The results showed that the upwind and downwind areas were extended far to the northwest and the southeast as a result of high wind speeds in the spring and winter, while they were restricted on the fringe of Seoul in the summer and fall.