• Title/Summary/Keyword: Numerical schemes

Search Result 758, Processing Time 0.025 seconds

Comparative Analysis of Finitc Difference Methods for the St, Venant Equation (St. Venant식에 관한 유한 차분법의 비교 분석)

  • 이상호;이길성
    • Water for future
    • /
    • v.21 no.2
    • /
    • pp.173-182
    • /
    • 1988
  • Characteristics of finite difference schemes for St. Venant equation were compared with two input cases. One is the monoclinal wave which has large friction slope without discontinuity and the other is the shock wave with discontinuity. For monoclinal wave, Keller Box scheme is the best in terms of accuracy, efficiency and stability when two parameters were selected with a rele : $0.5{\leq}{\theta}{\leq}1.0$, ${\theta}+{\psi}$=1, But for shock wave only the Preissmann type of parameter ${\psi}$(=0.5) showed stable results. Numerical experiments of monoclinal wave showed that Lax-Wendroff and Richtmyer schemes were more stable than leap Frog and more accurate than Lax-Fredrich scheme. For shock wave Lax-Fredrich showed larger numerical dissipation than other explicit schemes and Leap Frog produced slower mass transport.

  • PDF

A Numerical Study of Shock Wave/Boundary Layer Interaction in a Supersonic Compressor Cascade

  • Song, Dong-Joo;Hwang, Hyun-Chul;Kim, Young-In
    • Journal of Mechanical Science and Technology
    • /
    • v.15 no.3
    • /
    • pp.366-373
    • /
    • 2001
  • A numerical analysis of shock wave/boundary layer interaction in transonic/supersonic axial flow compressor cascade has been performed by using a characteristics upwind Navier-Stokes method with various turbulence models. Two equation turbulence models were applied to transonic/supersonic flows over a NACA 0012 airfoil. The results are superion to those from an algebraic turbulence model. High order TVD schemes predicted shock wave/boundary layer interactions reasonably well. However, the prediction of SWBLI depends more on turbulence models than high order schemes. In a supersonic axial flow cascade at M=1.59 and exit/inlet static pressure ratio of 2.21, k-$\omega$ and Shear Stress Transport (SST) models were numerically stables. However, the k-$\omega$ model predicted thicker shock waves in the flow passage. Losses due to shock/shock and shock/boundary layer interactions in transonic/supersonic compressor flowfields can be higher losses than viscous losses due to flow separation and viscous dissipation.

  • PDF

Comparison of Semi-Implicit Integration Schemes for Rate-Dependent Plasticity (점소성 구성식의 적분에 미치는 선형화 방법의 영향)

  • Yoon, Sam-Son;Lee, Soon-Bok
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.27 no.11
    • /
    • pp.1907-1916
    • /
    • 2003
  • During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

Localized Algorithm to Improve Connectivity and Topological Resilience of Multi-hop Wireless Networks

  • Kim, Tae-Hoon;Tipper, David;Krishnamurthy, Prashant
    • Journal of information and communication convergence engineering
    • /
    • v.11 no.2
    • /
    • pp.69-81
    • /
    • 2013
  • Maintaining connectivity is essential in multi-hop wireless networks since the network topology cannot be pre-determined due to mobility and environmental effects. To maintain the connectivity, a critical point in the network topology should be identified where the critical point is the link or node that partitions the network when it fails. In this paper, we propose a new critical point identification algorithm and also present numerical results that compare the critical points of the network and H-hop sub-network illustrating how effectively sub-network information can detect the network-wide critical points. Then, we propose two localized topological control resilient schemes that can be applied to both global and local H-hop sub-network critical points to improve the network connectivity and the network resilience. Numerical studies to evaluate the proposed schemes under node and link failure network conditions show that our proposed resilient schemes increase the probability of the network being connected in variety of link and node failure conditions.

Delta-form-based method of solving high order spatial discretization schemes for neutron transport

  • Zhou, Xiafeng;Zhong, Changming;Li, Fu
    • Nuclear Engineering and Technology
    • /
    • v.53 no.7
    • /
    • pp.2084-2094
    • /
    • 2021
  • Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor SN transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional SN codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the SN transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.

ON EXACT CONVERGENCE RATE OF STRONG NUMERICAL SCHEMES FOR STOCHASTIC DIFFERENTIAL EQUATIONS

  • Nam, Dou-Gu
    • Bulletin of the Korean Mathematical Society
    • /
    • v.44 no.1
    • /
    • pp.125-130
    • /
    • 2007
  • We propose a simple and intuitive method to derive the exact convergence rate of global $L_{2}-norm$ error for strong numerical approximation of stochastic differential equations the result of which has been reported by Hofmann and $M{\"u}ller-Gronbach\;(2004)$. We conclude that any strong numerical scheme of order ${\gamma}\;>\;1/2$ has the same optimal convergence rate for this error. The method clearly reveals the structure of global $L_{2}-norm$ error and is similarly applicable for evaluating the convergence rate of global uniform approximations.

A TREATMENT OF CONTACT DISCONTINUITY FOR CENTRAL UPWIND SCHEME BY CHANGING FLUX FUNCTIONS

  • Shin, Moungin;Shin, Suyeon;Hwang, Woonjae
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.17 no.1
    • /
    • pp.29-45
    • /
    • 2013
  • Central schemes offer a simple and versatile approach for computing approximate solutions of nonlinear systems of hyperbolic conservation laws. However, there are large numerical dissipation in case of contact discontinuity. We study semi-discrete central upwind scheme by changing flux functions to reduce the numerical dissipation and we perform numerical computations for various problems in case of contact discontinuity.

NUMERICAL MODELING OF NON-CAPACITY MODEL FOR SEDIMENT TRANSPORT BY CENTRAL UPWIND SCHEME

  • S. JELTI;A. CHARHABIL;J. EL GHORDAF
    • 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.

Multivariate Control Charts for Means and Variances with Variable Sampling Intervals

  • Kim, Jae-Joo;Cho, Gyo-Young;Chang, Duk-Joon
    • Journal of Korean Society for Quality Management
    • /
    • v.22 no.1
    • /
    • pp.66-81
    • /
    • 1994
  • Several sample statistics to simultaneously monitor both means and variances for multivariate quality characteristics under multivariate normal process are proposed. Performances of multivariate Shewhart schemes and cumulative sum(CUSUM) schemes are evaluated for matched fixed sampling interval(FSI) and variable sampling interval(VSI) feature. Numerical results show that multivariate CUSUM charts are more efficient than Shewhart charts for small or moderate shifts and VSI feature is more efficient than FSI feature.

  • PDF

CENTRAL SCHEMES WITH LAX-WENDROFF TYPE TIME DISCRETIZATIONS

  • Shin, Su-Yeon;Hwang, Woon-Jae
    • Bulletin of the Korean Mathematical Society
    • /
    • v.48 no.4
    • /
    • pp.873-896
    • /
    • 2011
  • The semi-discrete central scheme and central upwind scheme use Runge-Kutta (RK) time discretization. We do the Lax-Wendroff (LW) type time discretization for both schemes. We perform numerical experiments for various problems including two dimensional Riemann problems for Burgers' equation and Euler equations. The results show that the LW time discretization is more efficient in CPU time than the RK time discretization while maintaining the same order of accuracy.