• 제목/요약/키워드: Runge-Kutta methods

검색결과 106건 처리시간 0.023초

PCV 밸브의 설계 프로그램 개발에 관한 연구 (A Study on Development of Design Program for PCV Valve)

  • 이종훈;;이연원;김영덕
    • 한국마린엔지니어링학회:학술대회논문집
    • /
    • 한국마린엔지니어링학회 2005년도 전기학술대회논문집
    • /
    • pp.228-232
    • /
    • 2005
  • Automobiles are very important as modern society is developed. Increase of the number of the automobiles cause environmental problem, that is, air pollution. So, many countries are adopting a environmental law. Automobile manufacturing companies have developing methods to prevent air pollution with increase of the efficiency of automotive engines. PCV(Positive Crankcase Ventilation) system which is one of them is made by the closed loop that consists of combustion chamber, crankcase, manifold suction tube and manifold. PCV valve is attached on manifold tube to control the flowrate of blowby gas. PCV valve is an important part in this system but it is difficult to design PCV valve which satisfies the required flowrate of blowby gas. In this study, our purpose is to help a PCV valve designer with the development of a design program. We used 4th order Runge-Kutta method and Bernoulli's equation to analyze the spool dynamic motion. By the comparison between our program and experiment, we think that a PCV designer can use our program in their work place.

  • PDF

일정체적 양단고정 기둥의 정·동적 최적형상 (Static and Dynamic Optimal Shapes of Both Clamped Columns with Constant Volume)

  • 이병구;김석기
    • 한국강구조학회 논문집
    • /
    • 제19권1호
    • /
    • pp.99-106
    • /
    • 2007
  • 이 논문은 일정체적 양단고정 기둥의 정 동적 최적형상에 관한 연구이다. 기둥의 단면은 정다각형이며, 단면깊이는 포물선으로 변화하는 변단면이다. 축방향 압축하중이 작용하는 기둥의 고유진동수 및 좌굴하중을 산정하는 수치해석 기법을 개발하였다. 그러한 기둥의 자유진동을 지배하는 미분방정식을 유도하고 Runge-Kutta법과 Regula-Falsi법을 이용하여 고유진동수를 산정하였다. 수치해석의 결과로부터 얻어진 하중-고유진동수 사이의 관계를 이용하여 기둥의 좌굴하중을 산정하였다. 기둥의 변수연구를 통하여 동적 안정영역, 동적 최적형상 및 최강기둥의 형상을 산출하였다.

Error Control Policy for Initial Value Problems with Discontinuities and Delays

  • Khader, Abdul Hadi Alim A.
    • Kyungpook Mathematical Journal
    • /
    • 제48권4호
    • /
    • pp.665-684
    • /
    • 2008
  • Runge-Kutta-Nystr$\"{o}$m (RKN) methods provide a popular way to solve the initial value problem (IVP) for a system of ordinary differential equations (ODEs). Users of software are typically asked to specify a tolerance ${\delta}$, that indicates in somewhat vague sense, the level of accuracy required. It is clearly important to understand the precise effect of changing ${\delta}$, and to derive the strongest possible results about the behaviour of the global error that will not have regular behaviour unless an appropriate stepsize selection formula and standard error control policy are used. Faced with this situation sufficient conditions on an algorithm that guarantee such behaviour for the global error to be asympotatically linear in ${\delta}$ as ${\delta}{\rightarrow}0$, that were first derived by Stetter. Here we extend the analysis to cover a certain class of ODEs with low-order derivative discontinuities, and the class of ODEs with constant delays. We show that standard error control techniques will be successful if discontinuities are handled correctly and delay terms are calculated with sufficient accurate interpolants. It is perhaps surprising that several delay ODE algorithms that have been proposed do not use sufficiently accurate interpolants to guarantee asymptotic proportionality. Our theoretical results are illustrated numerically.

Optimum Radius Size between Cylindrical Ion Trap and Quadrupole Ion Trap

  • Chaharborj, Sarkhosh Seddighi;Kiai, Seyyed Mahmod Sadat;Arifin, Norihan Md;Gheisari, Yousof
    • Mass Spectrometry Letters
    • /
    • 제6권3호
    • /
    • pp.59-64
    • /
    • 2015
  • Quadrupole ion trap mass analyzer with a simplified geometry, namely, the cylindrical ion trap (CIT), has been shown to be well-suited using in miniature mass spectrometry and even in mass spectrometer arrays. Computation of stability regions is of particular importance in designing and assembling an ion trap. However, solving CIT equations are rather more difficult and complex than QIT equations, so, analytical and matrix methods have been widely used to calculate the stability regions. In this article we present the results of numerical simulations of the physical properties and the fractional mass resolutions m/Δm of the confined ions in the first stability region was analyzed by the fifth order Runge-Kutta method (RKM5) at the optimum radius size for both ion traps. Because of similarity the both results, having determining the optimum radius, we can make much easier to design CIT. Also, the simulated results has been performed a high precision in the resolution of trapped ions at the optimum radius size.

불응축가스가 평판위 응축열전달에 미치는 영향에 관한 연구 (A study on effect of heat transfer of condensation including noncondensable gas over a flat plate)

  • 양대일;정형호
    • Journal of Advanced Marine Engineering and Technology
    • /
    • 제24권1호
    • /
    • pp.25-30
    • /
    • 2000
  • In present paper, mass transfer over a flat plate with film condensation including noncondesable gas is analyzed with the help of similarity methods. Couette flow was assumed in liquid film and boundary-layer approximation was used in the ambient flow. Governing equations were transformed into the ordinary differential equtions by the similarity methods. Runge-Kutta and shooting method were used in order to fine the effect of mass transfer on the velocity and concentrations at the liquid-vapor interface.

  • PDF

NUMERICAL METHODS FOR CAVITATING FLOW

  • SHIN Byeong Rog
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2001년도 추계 학술대회논문집
    • /
    • pp.1-9
    • /
    • 2001
  • In this paper, some numerical methods recently developed for gas-liquid two-phase flows are reviewed. And then, a preconditioning method to solve cavitating flow by the author is introduced. This method employs a finite-difference Runge-Kutta method combined with MUSCL TVD scheme, and a homogeneous equilibrium cavitation model. So that it permits to treat simply the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Finally, numerical results such as detailed observations of the unsteady cavity flows, a sheet cavitation break-off phenomena and some data related to performance characteristics of hydrofoils are shown.

  • PDF

Error Control Strategy in Error Correction Methods

  • KIM, PHILSU;BU, SUNYOUNG
    • Kyungpook Mathematical Journal
    • /
    • 제55권2호
    • /
    • pp.301-311
    • /
    • 2015
  • In this paper, we present the error control techniques for the error correction methods (ECM) which is recently developed by P. Kim et al. [8, 9]. We formulate the local truncation error at each time and calculate the approximated solution using the solution and the formulated truncation error at previous time for achieving uniform error bound which enables a long time simulation. Numerical results show that the error controlled ECM provides a clue to have uniform error bound for well conditioned problems [1].

ANALYTICAL AND NUMERICAL SOLUTIONS OF A CLASS OF GENERALISED LANE-EMDEN EQUATIONS

  • RICHARD OLU, AWONUSIKA;PETER OLUWAFEMI, OLATUNJI
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제26권4호
    • /
    • pp.185-223
    • /
    • 2022
  • The classical equation of Jonathan Homer Lane and Robert Emden, a nonlinear second-order ordinary differential equation, models the isothermal spherical clouded gases under the influence of the mutual attractive interaction between the gases' molecules. In this paper, the Adomian decomposition method (ADM) is presented to obtain highly accurate and reliable analytical solutions of a class of generalised Lane-Emden equations with strong nonlinearities. The nonlinear term f(y(x)) of the proposed problem is given by the integer powers of a continuous real-valued function h(y(x)), that is, f(y(x)) = hm(y(x)), for integer m ≥ 0, real x > 0. In the end, numerical comparisons are presented between the analytical results obtained using the ADM and numerical solutions using the eighth-order nested second derivative two-step Runge-Kutta method (NSDTSRKM) to illustrate the reliability, accuracy, effectiveness and convenience of the proposed methods. The special cases h(y) = sin y(x), cos y(x); h(y) = sinh y(x), cosh y(x) are considered explicitly using both methods. Interestingly, in each of these methods, a unified result is presented for an integer power of any continuous real-valued function - compared with the case by case computations for the nonlinear functions f(y). The results presented in this paper are a generalisation of several published results. Several examples are given to illustrate the proposed methods. Tables of expansion coefficients of the series solutions of some special Lane-Emden type equations are presented. Comparisons of the two results indicate that both methods are reliably and accurately efficient in solving a class of singular strongly nonlinear ordinary differential equations.

Nonlinear vibration of oscillatory systems using semi-analytical approach

  • Bayat, Mahmoud;Bayat, Mahdi;Pakar, Iman
    • Structural Engineering and Mechanics
    • /
    • 제65권4호
    • /
    • pp.409-413
    • /
    • 2018
  • In this paper, He's Variational Approach (VA) is used to solve high nonlinear vibration equations. The proposed approach leads us to high accurate solution compared with other numerical methods. It has been established that this method works very well for whole range of initial amplitudes. The method is sufficient for both linear and nonlinear engineering problems. The accuracy of this method is shown graphically and the results tabulated and results compared with numerical solutions.

다중격자 기법을 적용한 극초음속 평형 유동장 계산 (Application of Multigrid Method for Computing Hypersonic, Equilibrium Flows)

  • 김성수;김종암;노오현
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 1999년도 춘계 학술대회논문집
    • /
    • pp.23-28
    • /
    • 1999
  • A mutigrid convergence acceleration technique is presented for computing hypersonic inviscid and viscous flows in equilibrium state. The governing equations are solved using an explicit Runge-Kutta method. Curve fitting data in NASA Reference Publication 1181, 1260 are used to calculate equilibrium properties. In order to ensure stability, damped prolongation and modified implicit residual smoothing are proposed. Blunt body test cases are presented to demonstrate the robustness and the efficiency in performance of the proposed methods

  • PDF