• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.2B
• /
• pp.123-130
• /
• 2012

MLP (Multi dimensional Limiting Process) is implemented to simulate shallow water flows, and its performance over conventional TVD limiters in multidimensional flows is verified through several numerical simulations. MLP was developed to control oscillations for multi-dimensional compressible flows and proved to improve accuracy, efficiency and robustness in compressible flows. In this study, we applies MLP to modeling shallow water equations(SWEs) given that the SWEs are amenable to be solved using the large range of numerical methods developed to deal with compressible flows and MLP has been yet used for SWEs. Simulation results through the benchmark tests show that MLP has favorable features such as numerical oscillation control and convergence behaviors comparable to the conventional limiters. Both numerical accuracy and stability are improved in multi-dimensional discontinuous flows.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.22-27
• /
• 2006

This paper describes the development of the stream-tube based dispersion model for modeling contaminant transport in open channels. The operator-splitting technique is employed to separate the 2D contaminant transport equation into the pure advection and pure dispersion equations. Then the total variation diminishing (TVD) schemes are combined with the second-order Lax-Wendroff and third-order QUICKEST explicit finite difference schemes respectively to solve the pure advection equation in order to prevent the occurrence of numerical oscillations. Due to various limiters owning different features, the numerical tests for 1D pure advection and 2D dispersion are conducted to evaluate the performance of different TVD schemes firstly, then the TVD schemes are applied to experimental data for simulating the 2D mixing in a straight trapezoidal channel to test the model capability. Both the numerical tests and model application show that the TVD schemes are very competent for solving the advection-dominated transport problems.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.

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

A new treatment of cell-interface flux in AUSM-type methods is introduced to reduce the numerical dissipation. Through analysis of TVD limiters, a criterion for the more accurate prediction of cell-interface state is found out and M-AUSMPW+ is developed by determining the transferred property newly and appropriately within the criterion. The superiority of M-AUSMPW+ is clearly revealed in multi-dimensional flow problems. It can eliminate numerical dissipation effectively in a non-flow aligned grid system. As a result, M-AUSMPW+ is shown to be much more accurate and effective than other previous schemes in multi-dimensional problems. Through a stationary contact discontinuity, a vortex flow, a shock wave/boundary layer interactions and viscous shock tube problems, it is verified that accuracy of M-AUSMPW+ is improved.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.57-66
• /
• 1999

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit. The waves reflect repeatedly on the center axis and the sonic surface in the shear layer. The pressure difference is resolved across these reflected waves. In this paper, the axi-symmetric Navier-Stokes equation has been used with the κ-ε turbulence model. The second order TVD scheme with flux limiters, based on the flux vector split with the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

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

We performed numerical analysis of base drag phenomenon, when a projectile with backward step flies into atmosphere at supersonic speed. We compared with other researchers. From our previous studies that were 2-dimensional simulation, we found out from sophisticated simulations that need dense mesh points to compare base pressure and velocity profile after from base with experimental data. Therefore, we focus on high order spatial disceretization over 3rd order with TVD such as MUSCL TVD 3rd, 5th, and WENO 5th order, and Limiters such as minmod, Triad. Moreover, we enforce to flux averaging schemes such as Roe, RoeM, HLLE, AUSMDV. In present, one dimensional result of Euler tests, there are Sod, Lax, Shu-Osher and interacting blast wave problems. AUSMDV as a flux averaging scheme with MUSCL TVD 5th order as spatial resolution is good agreement with exact solutions than other combinations. We are carrying out the same approaches into 3-dimensional base flow only candidate flux schemes that are Roe, AUSMDV. Additionally, turbulence models are used in 3-dimensional flow, one is Menter s SST DES model and another is Sparlat-Allmaras DES/DDES model in Navier-Stokes equations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.2
• /
• pp.10-16
• /
• 2003

The unsteady supersonic flow over tandem cavities has been analyzed by the integration of Navier-Stokes equations with the k-$\varepsilon$ 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. The upwind TVD scheme based on the flux vector split with the van Leer limiters is used. The results show the principal frequency is very reasonable. The principal frequency of the rear cavity due to the front cavity has been analyzed by the combination of the several aspect ratios of cavities. In the case of the front cavity of low aspect ratio, the frequencies of tandem cavities are almost same, because two shear layers developed from each cavity are mixed and developed to one shear layer. However, in the case of the front cavity of high aspect ratio, the characteristis of frequency are very different, because the second shear layer is developed in the diffused first shear layer.

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

• Journal of computational fluids engineering
• /
• v.10 no.4 s.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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.7
• /
• pp.1-11
• /
• 2006

Through the analysis of conventional TVD limiters, a new multi-dimensional limiting function is derived for an oscillation control in multi-dimensional flows. Then, Multi-dimensional Limiting Process (MLP) is developed with the multi-dimensional limiting function. The major advantage of MLP is to prevent oscillations across a multi-dimensional discontinuity, and it is readily compatible with more than 3rd order spatial interpolation. Moreover, MLP shows a good convergence characteristic in a steady problem and it is very simple to be implemented. Through numerical test cases, it is verified that MLP substantially improves accuracy, efficiency and robustness both in continuous and discontinuous flows.