• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.13-17
• /
• 2006

RoeM and AUSMPW+ schemes are two of the most accurate and efficient schemes which are recently developed for the analysis of single phase gas dynamics. In this paper, we developed two-phase versions of these schemes for the analysis of gas-liquid large density ratio two-phase flow. We adopt homogeneous equilibrium model (HEM) using mass fraction to describe different two phases. In the Eulerian-Eulerian framework, HEM assumes dynamic and thermal equilibrium of the two phases in the same computational mesh. From the mixture equation of state (EOS), we derived new shock-discontinuity sensing term (SDST), which is commonly used in RoeM and AUSMPW+ for the stable numerical flux calculation. The proposed two-phase versions of RoeM and AUSMPW+ schemes are applied on several air-water two-phase test problems. In spite of the large discrepancy of material properties such as density, enthalpy, and speed of sound, the numerical results show that both schemes provide very satisfactory solutions.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.11-19
• /
• 2010

A code is developed to simulate incompressible free surface flows using the Roe's flux-difference splitting scheme. An interface of two fluids is considered as a moving contact discontinuity. The continuities of pressure and normal velocity across the interface are enforced by the conservation law in the integral sense. The fluxes are computed using the Roe's flux-difference splitting scheme for two incompressible fluids. The interface can be identified based on the computed density distribution. However, no additional treatment is required along the interface during the whole computations. Complicated time evolution of the interface including topological change can be captured without any difficulties. The developed code is applied to simulate the Rayleigh-Taylor instability of two incompressible fluids in the density ratio of 7.2:1 and the broken dam problem of water-air. The present results are compared with other available results and good agreements are achieved for the both cases.

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

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.48-58
• /
• 1999

The Characteristic Flux Difference Splitting (CFDS) scheme designed to adapt the characteristic boundary conditions at the wall and inflow/outflow boundary planes satisfies Roe's property U, although the CFDS Jacobian matrix is decomposed by a product of elaborate transformation matrices and explicit eigenvalue matrix. When the CFDS algorithm, thus a variant of Roe's scheme, is applied straightforwardly to hypersonic flows over a blunt body, the strong bow shock gradually breaks down near the stagnation point. This numerical instability is widely observed by many researchers employing flux-difference method, known in the literature as the carbuncle phenomenon. Many remedies have been proposed and resulted in partial cures. When the idea of Sanders et al. which identifies the minimum eigenvalues near the discontinuity present is applied to CFDS method, it is shown that the instability problem can be controlled successfully. A few flux splitting methods have also been tested and results are compared against the Nakamori's Mach 8 blunt body flow.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.569-574
• /
• 2020

Though the discontinuous Galerkin (DG) method has been developed and applied to shallow water equations mainly in explicit schemes, they have been criticized for the limitation in treatment of bottom friction terms and severe CFL conditions. In this study, an implicit scheme is devised and applied to some representative benchmark problems. The linear triangular elements were employed and the Roe numerical fluxes were adopted for convective fluxes. To preserve TVD property, the slope limiter was employed. As the case studies, the model is applied to the flow around the cylinders and the dam-break flow. Then, the results are compared with the experimental and numerical data of previous studies and good agreements were observed.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.1
• /
• pp.1-16
• /
• 2009

In this paper, we introduce two-phase versions of RoeM and AUSMPW+ schemes. Both schemes are originally developed for the gas dynamic problems, and have shown superior accuracy, efficiency and robustness. A new shock discontinuity sensing term is derived from the mixture equation of state, which is commonly used in the RoeM and AUSMPW+ schemes for the stable numerical flux calculation. The developed two-phase versions of the schemes are applied to several liquid-gas, large property discrepancy two-phase test problems, including several shock stability test problems. The results show that both schemes maintain the merits exhibited in gas dynamic problems even in two-phase flows.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.43-59
• /
• 2000

Characteristic boundary conditions are discussed in conjunction with a flux-difference splitting formulation as modified from Roe's linearization. Details of how one can implement the characteristic boundary conditions which are made compatible with the interior point formulation are described for different types of boundaries including subsonic outflow and adiabatic wall. The validity of boundary conditions are demonstrated through computation of transonic airfoil, supersonic ogive-cylinder, hypersonic cylinder, and S-duct internal flows. The computed wall pressure distributions are compared with published experimental and computed data. Objectives of this paper are thus to give insight of formulation procedure of a flux-difference splitting method and to pave ways for other users to adopt present boundary procedure on their numerical methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.