• Title/Summary/Keyword: Flux-difference splitting

Search Result 67, Processing Time 0.02 seconds

Numerical Simulation of Free Surface Flows Using the Roe's Flux-difference Splitting Scheme (Roe의 Flux-difference Splitting 기법을 이용한 자유표면 유동 모사)

  • Shin, Sang-Mook;Kim, In-Chul;Kim, Yong-Jig
    • 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.

Compressor Cascade Flow Analysis by Using Upwind Flux Difference Splitting Method (풍상차분법을 이용한 압축기 익렬유동 해석)

  • 권창오;송동주;강신형
    • Transactions of the Korean Society of Mechanical Engineers
    • /
    • v.18 no.3
    • /
    • pp.653-661
    • /
    • 1994
  • In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 supersonic/transonic compressor cascade flow. H-type grid was chosen for its simplicity in applying cyclic tridiagonal matrix algorithm along with conventional slip/no-slip boundary conditions. The thin-layer algebraic model of Baldwin-Lomax was employed for the calculation of turbulent flows. The test case inlet Mach No. was 1.612 and inlet/exit pressure ratio($P_2/P_1$) was 2.15. The results were compared with experimental results from current method were compared well in suction surface with the experiments and other computational results; however, not well in pressure surface. It might be due to the complex flowfields such as shock/boundary layer interaction, turbulence, and flow separation, etc. In the future, a proper turbulence modelling and adaptive grid system will be studied to improve the solution quality.

Applications of Characteristic Boundary Conditions within CFDS Numerical Framework (CFDS기법에 연계된 특성경계조건에 응용성에 대한 소개)

  • Hong S. K.;Lee K. S.
    • 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.

  • PDF

Characteristic Flux-Difference Improvement for Inviscid and Viscous Hypersonic Blunt Body Flows

  • Lee Gwang-Seop;Hong Seung-Gyu
    • 한국전산유체공학회:학술대회논문집
    • /
    • 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.

  • PDF

Numerical Analysis of Nonequilibrium Chemically Reacting Inviscid flow over Blunt-bodies Using Upwind Method (Upwind 방법을 이용한 무딘물체 주위의 화학적 비평형 비점성 유동장의 수치 해석)

  • Seo Jeong Il;Song Dong Joo
    • 한국전산유체공학회:학술대회논문집
    • /
    • 1997.10a
    • /
    • pp.99-105
    • /
    • 1997
  • A finite-difference method based on conservative supra characteristic method type upwind flux difference splitting has been developed to study the nonequilibrium chemically reacting inviscid flow. For nonequilibrium air, NS-1 species equations were strongly coupled with flowfield equations through convection and species production terms. Inviscid nonequilibrium chemically reacting air mixture flows over Blunt-body were solved to demonstrate the capability of the current method. At low altitude flight conditions the nonequilibrium air models predicted almost the same temperature, density and pressure behind the shock as equilibrium flow: however, at high altitudes they showed substantial differences due to nonequilibrium chemistry effect. The new nonequilibrium chemically reacting upwind flux difference splitting mettled can be extended to viscous flow and multi-dimensional flow conditions.

  • PDF

The Application of Multigrid Algorithm to Low-Speed Precondition

  • Yang, Zhong;Xu, Jianzhong
    • Proceedings of the Korean Society of Propulsion Engineers Conference
    • /
    • 2008.03a
    • /
    • pp.666-670
    • /
    • 2008
  • The low-speed preconditioning technique is applied to solve the compressible Reynolds averaged Navier-Stokes equations for low-speed flows. The space discretization is based on Roe's flux-difference splitting with third-order-accurate MUSCL extrapolation. Time integration is performed employing a diagonal approximate factorization algorithm. The dual-time stepping has been incorporated to solve the unsteady flows. Full multigrid method is implemented to accelerate the convergence rate. To verify the algorithms several cases have been tested. Demonstrated the improvement on convergence and quality of the solution.

  • PDF

Performance Enhancement Study Using Passive Control of Shock-Boundary Layer Interaction in a Transonic/Supersonic Compressor Cascade (천음속/초음속 압축기 익렬에서 Shock-Boundary Layer 상호작용의 수동적 제어에 의한 성능 향상 연구)

  • Kim, Sang-Deok;Gwon, Chang-O;Sa, Jong-Yeop
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.20 no.9
    • /
    • pp.2944-2952
    • /
    • 1996
  • In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 transoni $c^ersonic compressor cascade flow. First, the general characteristics of baseline cascade flow were analyzed. At freestream Mach n.1.612 and exit/inlet pressure ratio 2.15, the results from current laminar flow were compared well in suction surface with the experiment; however, not well in pressure surface. Second, numerical study of the transoni $c^ersonic compressor cascade flow demonstrated the effectiveness of a passive control by the various size cavities. A cavity under the shock foot point at the suction surface of the blades was used as a passive control. The passive control of shock-boundary layer interaction by a cavity reduced total pressure losses. The effect of cavity length and depth was studied. The total pressure loss was reduced by about 10% and the isentropic efficiency was improved slightly. The effect of cavity depth in current study(d/l = 0.05, 0.02) was not found strong. Further adequate turbulence modeling and TVD schemes would help to capture the shock more accurately and increase the effectiveness of the current shock-boundary layer interaction study using upwind flux difference splitting computational methods.thods.

An Implementation of the Robust Inviscid Wall Boundary Condition in High-Speed Flow Calculations

  • Kim, Moon-Sang;Jeon, Byung-Woo;Kim, Yong-Nyun;Kwon, Hyeok-Bin;Lee, Dong-Ho
    • Journal of Mechanical Science and Technology
    • /
    • v.15 no.5
    • /
    • pp.671-680
    • /
    • 2001
  • Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

  • PDF

SIMULATION OF FLUID-STRUCTURE INTERACTION OF A TOWED BODY USING AN ASYMMETRIC TENSION MODEL (비대칭 장력 모델을 이용한 예인 물체의 유체-구조 상호작용 모사)

  • Shin, Sang-Mook
    • Journal of computational fluids engineering
    • /
    • v.16 no.1
    • /
    • pp.7-13
    • /
    • 2011
  • The fluid-structure interaction of a towed body is simulated using a developed code, which is based on the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method. To improve the stability in the coupling between the fluid and structure domains, a scheme is used, in which the effects of structure deformation are treated implicitly. The developed code is validated for the fluid-structure interaction problem through comparisons with other results on the vortex-induced vibration of elastically mounted cylinders. To simulate behavior of a towed body, an asymmetric tension modelling for a towing cable is suggested. In the suggested model, the tension is proportional to the elongation of the cable, but the cable has no effect on the body motion whenever the distance between the endpoints of the cable is smaller than the original length of the cable. The fluid-structure interactions of a towed body are simulated on the basis of different parameters of the towing cables. It is observed that the suggested tension model predicts the snapping for a shorter towing cable, which is in accordance with the reported results.

NUMERICAL SIMULATION OF THREE-DIMENSIONAL INTERNAL WAVES USING THE FDS SCHEME ON THE HCIB METHOD (FDS 기법과 HCIB법을 이용한 3차원 내면파 수치 모사)

  • Shin, Sang-Mook
    • Journal of computational fluids engineering
    • /
    • v.17 no.1
    • /
    • pp.8-15
    • /
    • 2012
  • A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.