• Title/Summary/Keyword: fourth-order Runge-Kutta method

Search Result 69, Processing Time 0.028 seconds

Computational Study on the Hemodynamic Behaviors of the Human Cardiovascular System with an Acute Arteriovenous Fistula (급성 동정맥루를 포함하는 인체 심혈관계의 혈류역학적 거동에 관한 수치 해석적 연구)

  • 변수영;손정락;심은보;노승탁
    • Journal of Biomedical Engineering Research
    • /
    • v.24 no.4
    • /
    • pp.329-337
    • /
    • 2003
  • Blood in congenital or acquired AY fistula(arteriovenous fistula) flows from arteries directly to veins. detouring peripheral micro-circulation. This makes a great effect on the hemodynamics of human cardiovascular system. In this study, a computational method using lumped parameter mode) was proposed to simulate the cardiovascular hemodynamics of patients with acute AV fistula The cardiovascular system model with a fistula compartment in left lower limb was built using 17 standard lumped compartments. Using fourth order Runge-Kutta method. we solved numerically the unsteady linear set of the ordinary differential equations resulting from application of Kirchhoff's law to the lumped parameter hemodynamic model. The baroreceptor reflex system was implemented to explain the auto-regulation effect of the cardiovascular system with acute AV fistula.

NUMERICAL ANALYSIS OF AXISYMMETIC SCREECH TONE FROM SUPERSONIC JET USING HIGH-ORDER HIGH-RESOLUTION COMPACT SCHEME (고차고해상도 수치기법을 이용한 초음속 제트 screech tone의 axisymmetric mode 해석)

  • Lee, I.C.;Lee, D.J.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2007.10a
    • /
    • pp.56-59
    • /
    • 2007
  • The screech tone of an underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. A fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. The Fourier transform of pressure signal at upstream shows the directivity pattern of the screech tone very clearly. Pressure signal is shown to observe the generation of the screech tone. Most importantly, we can simulate the axisymmetric mode change of the screech tone very precisely with the proposed method. It can be concluded that the basic phenomenon of the screech tone including its frequency can be calculated and its mode change can be simulated with inviscid Euler equations.

  • PDF

Numerical Simulation of Periodic and Oscillatory Problems by Using RK-Butcher Algorithms (RK-Butcher알고리듬의 사용에 의한 주기적 진동 문제의 수치적 시뮬레이션)

  • Park, Dae-Chul;Gopal, Devarajan;Murugesh, V.
    • Journal of the Institute of Convergence Signal Processing
    • /
    • v.9 no.1
    • /
    • pp.82-88
    • /
    • 2008
  • In this paper, Runge-Kutta (RK)-Butcher algorithm is proposed to study the periodic and oscillatory problems. Simulation results obtained using RK-Butcher algorithms and the classical fourth order Runge-Kutta (RK(4)) methods are compared with the exact solutions of a few periodic and oscillatory problems to confirm the performance of the proposed algorithm. The simulation results from RK-Butcher algorithms are always found to be very close to the exact solutions of these problems. Further, it is found that the RK-Butcher algorithm is superior when compared to RK(4) methods in terms of accuracy. The RK-Butcher algorithm can be easily implemented in a programming language and a more accurate solution may be obtained for any length of time. RK-Butcher algorithm is applicable as a good numerical algorithm for studying the problems of orbit and two body as it gives the nearly identical solutions.

  • PDF

Entropy Generation Minimization in MHD Boundary Layer Flow over a Slendering Stretching Sheet in the Presence of Frictional and Joule Heating

  • Afridi, Muhammad Idrees;Qasim, Muhammad;Khan, Ilyas
    • Journal of the Korean Physical Society
    • /
    • v.73 no.9
    • /
    • pp.1303-1309
    • /
    • 2018
  • In the present paper, we study the entropy analysis of boundary layer flow over a slender stretching sheet under the action of a non uniform magnetic field that is acting perpendicular to the flow direction. The effects of viscous dissipation and Joule heating are included in the energy equation. Using similarity transformation technique the momentum and thermal boundary layer equations to a system of nonlinear differential equations. Numerical solutions are obtained using the shooting and fourth-order Runge-Kutta method. The expressions for the entropy generation number and Bejan number are also obtained using a suggested similarity transformation. The main objective of this article is to investigate the effects of different governing parameters such as the magnetic parameter ($M^2$), Prandtl number (Pr), Eckert number (Ec), velocity index parameter (m), wall thickness parameter (${\alpha}$), temperature difference parameter (${\Omega}$), entropy generation number (Ns) and Bejan number (Be). All these effects are portrayed graphically and discussed in detail. The analysis reveals that entropy generation reduces with decreasing wall thickness parameter and increasing temperature difference between the stretching sheet and the fluid outside the boundary layer. The viscous and magnetic irreversibilities are dominant in the vicinity of the stretching surface.

A Hydrodynamic-Acoustic Splitting Method for Aeroacoustic Noise Prediction of Wall-bounded Shear Flow (벽면 전단 유동의 공력 소음 해석을 위한 유동-음향 분리 기법)

  • Seo, J.H.;Moon, Young-J.
    • Proceedings of the KSME Conference
    • /
    • 2003.04a
    • /
    • pp.2060-2065
    • /
    • 2003
  • Aeolian tone generation from a two dimensional circular cylinder is numerically investigated via direct numerical simulation and hydrodynamic-acoustic splitting method. All governing equation are spatially discretized with the sixth-order compact scheme and fourth-order Runge-Kutta method to avoid excessive numerical dissipations and dispersions of acoustic quantities. Comparisons of two results show that the previous splitting method can not accurately predict the aeroacoustic noise of wall bounded shear flow. In this study, a perturbation viscous term and a new energy equation have been developed. This modified splitting method accurately predicts aeroacoustic noise from wall-bounded shear flow. The present results agree very well with the direct numerical simulation solution.

  • PDF

Resonant response of spar-type floating platform in coupled heave and pitch motion

  • Choi, E.Y.;Cho, J.R.;Jeong, W.B.
    • Structural Engineering and Mechanics
    • /
    • v.65 no.5
    • /
    • pp.513-521
    • /
    • 2018
  • In this paper, the resonance response of spar-type floating platform in coupled heave and pitch motion is investigated using a CPU time-effective numerical method. A coupled nonlinear 2-DOF equation of motion is derived based on the potential wave theory and the rigid-body hydrodynamics. The transient responses are solved by the fourth-order Runge-Kutta (RK4) method and transformed to the frequency responses by the digital Fourier transform (DFT), and the first-order approximation of heave response is analytically derived. Through the numerical experiments, the theoretical derivation and the numerical formulation are verified from the comparison with the commercial software AQWA. And, the frequencies of resonance arising from the nonlinear coupling between heave and pitch motions are investigated and justified from the comparison with the analytically derived first-order approximation of heave response.

Study of Greitzer's B-Parameter Model Using ANOVA & Taguchi Method

  • Ng E. Y-K;Liu N.;Tan S. Y.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2003.10a
    • /
    • pp.197-199
    • /
    • 2003
  • In this work, the Greitzer's B-parameter model is applied for analyzing the stall and surge characteristics. The four parameters in the model are highlighted in order to establish the influence of each parameter on the system. First of all, the governing equations of stall and surge behavior are solved numerically using fourth-order Runge-Kutta method. The Taguchi method is then used to analyze the results generated to obtain the extent of effects of the parameters on the system by varying the parameters in a series of combinations. Finally, a thorough analysis is carried out on the results generated from the Taguchi method and the graphs.

  • PDF

Risk Model for the Safety Evaluation of Dam and Levee : I. Theory and Model (댐 및 하천제방에 대한 위험도 해석기법의 개발 : I. 이론 및 모형)

  • Han, Geon-Yeon;Lee, Jong-Seok;Kim, Sang-Ho
    • Journal of Korea Water Resources Association
    • /
    • v.30 no.6
    • /
    • pp.679-690
    • /
    • 1997
  • The risk assessment model for hydrlolgic safety analysis of dam and levee in developed by using Monte-Carlo and AFOSM (Advanced First-Order Second-Moment) method. The fault tree analysis and four phases approach are presented for the safety eveluation of risk of dam and levee. The risk model consists of rainfall-runoff analysis, reservoir routing and channel routing considering the variations in the model parameter. For the rainfall-runoff analysis, KRRL method is adopted with 200-year precipitation and PMP (Probable Maximum Precipitation). Reservoir routing is performed by fourth order Runge-Kutta method and channel routing by standard step method. The suggested model will contribute to safety evaluation of dam and levee and their rehabilitation decision problem.

  • PDF

Highly accurate family of time integration method

  • Rezaiee-Pajand, Mohammad;Esfehani, S.A.H.;Karimi-Rad, Mahdi
    • Structural Engineering and Mechanics
    • /
    • v.67 no.6
    • /
    • pp.603-616
    • /
    • 2018
  • In this study, the acceleration vector in each time step is assumed to be a mth order time polynomial. By using the initial conditions, satisfying the equation of motion at both ends of the time step and minimizing the square of the residual vector, the m+3 unknown coefficients are determined. The order of accuracy for this approach is m+1, and it has a very low dispersion error. Moreover, the period error of the new technique is almost zero, and it is considerably smaller than the members of the Newmark method. The proposed scheme has an appropriate domain of stability, which is greater than that of the central difference and linear acceleration techniques. The numerical tests highlight the improved performance of the new algorithm over the fourth-order Runge-Kutta, central difference, linear and average acceleration methods.

Nonlinear dynamic stability and vibration analysis of sandwich FG-CNTRC shallow spherical shell

  • Kamran Foroutan;Akin Atas;Habib Ahmadi
    • Advances in nano research
    • /
    • v.17 no.2
    • /
    • pp.95-107
    • /
    • 2024
  • In this article, the semi-analytical method was used to analyze the nonlinear dynamic stability and vibration analysis of sandwich shallow spherical shells (SSSS). The SSSS was considered as functionally graded carbon nanotube-reinforced composites (FG-CNTRC) with three new patterns of FG-CNTRC. The governing equation was obtained and discretized utilizing the Galerkin method by implementing the von Kármán-Donnell nonlinear strain-displacement relations. The nonlinear dynamic stability was analyzed by means of the fourth-order Runge-Kutta method. Then the Budiansky-Roth criterion was employed to obtain the critical load for the dynamic post-buckling. The approximate solution for the deflection was represented by suitable mode functions, which consisted of the three modes of transverse nonlinear oscillations, including one symmetrically and two asymmetrical mode shapes. The influences of various geometrical characteristics and material parameters were studied on the nonlinear dynamic stability and vibration response. The results showed that the order of layers had a significant influence on the amplitude of vibration and critical dynamic buckling load.