• Title/Summary/Keyword: Runge Kutta Method

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

Three-dimensional dynamics of vortex-induced vibration of a pipe with internal flow in the subcritical and supercritical regimes

  • Duan, Jinlong;Chen, Ke;You, Yunxiang;Wang, Renfeng;Li, Jinlong
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • 제10권6호
    • /
    • pp.692-710
    • /
    • 2018
  • The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.

배급수계통에서 잔류염소 감소 특성 및 적용연구 (Modeling and Application of Chlorine Bulk Decay in Drinking Water Distribution System)

  • 안재찬;박창민;구자용
    • 상하수도학회지
    • /
    • 제19권4호
    • /
    • pp.487-496
    • /
    • 2005
  • Chlorine bulk decay tests were carried out by bottle test under controlled conditions in a laboratory. Experiments were performed at different temperatures: $5^{\circ}C$, $15^{\circ}C$, $25^{\circ}C$, and the water temperatures when samples were taken from the effluent just before entering to its distribution system. 38 bulk tests were performed for water of Al (water treatment plant), 4 bulk tests for A2 (large service reservoir), and A3(pumping station). Residual chlorine concentrations in the amber bottles were measured over time till about 100 hours and bulk decay coefficients were evaluated by assuming first-order, parallel first-order, second-order. and $n^{th}-order$ reaction. The $n^{th}-order$ coefficients were obtained using Fourth-order Runge-Kutta Method. A good-fit by the average coefficient of determination ($R^2$) was first-order ($R^2=0.90$) < parallel first-order ($R^2{_{fast}}=0.92$, $R^2{_{slow}}=0.95$) < second-order ($R^2=0.95$) < $n^{th}-order$ ($R^2=0.99$). But if fast reaction of parallel first-order bulk decay were applied to the effluent of large service reservoir with ca. 20 hours of travel time and slow reaction in the water distribution system following the first 20 hours, parallel first-order bulk decay would be best and easy for application of water quality modeling technique.

CWM 연료의 연소시 입자 가열속도와 입자 크기가 CWM 응집물 형성에 미치는 영향에 관한 전산모사 연구 (Simulation Study on the Effects of Heating Rate and Particle Size Distribution for the Formation of the Agglomerate During CWM Combustion)

  • 김수호;김영환;황갑성;홍성선
    • 공업화학
    • /
    • 제9권3호
    • /
    • pp.364-371
    • /
    • 1998
  • 이론적인 입자 응집모델이 CWM 응집물에서 인접한 석탄입자들 사이의 응집력을 예측하기 위해서 전개되었다. 유연탄은 약 $400^{\circ}C$ 또는 그 이상 가열될 때 석탄 입자들은 플라스틱성 거동으로의 전이 현상이 관찰될 수 있다. 석탄 입자들의 응집과정에서의 응집력은 입자들의 플라스틱성의 지속기간에 비례하는 것으로 알려져 있다. 본 연구에서는 가열단계에서 CWM 연료의 응집물 형성에 미치는 영향을 입자 응집모델에 의하여 고찰하였다. 실험에 사용한 전산모사 프로그램은 RKG(Runge-Kutta-Gill) 방법을 이용하여 포트란으로 작성을 하였다. 입자 응집 모델에 의한 CWM 연료의 입자 응집과정에서의 응집력은 입자 가열속도에 반비례하며, 입자크기에 비례하는 것으로 나타났다.

  • PDF

Basic Study on the Heat Transfer During Rapid Freezing of Syobean Seed by Liquid Nitrogen

  • Kawano, Toshio;Nakano, Kohei;Murata, Satoshi
    • 한국농업기계학회:학술대회논문집
    • /
    • 한국농업기계학회 1993년도 Proceedings of International Conference for Agricultural Machinery and Process Engineering
    • /
    • pp.443-452
    • /
    • 1993
  • Direct freezing tests of soybean seed by liquid nitrogen were carried out at various moisture contents and the following important conclusions were drawn from the results of temperature measurements of soybean seed and photographs of bubbles generated on its surface : 1) Assuming that the temperature gradient in a soybean seed is negligible because of its small seed size and the freezing ratio is followed the Heiss's formula, and a differential equation based on the heat energy balance was introduced . The equation was easily solved by the Runge-Kutta-Gill method and the predicted values of the temperature were in good agreement with the observed data. 2) The photographs of bubble generation during freezing showed the boiling mode was nucleate, and then the most suitable formula on the nucleate boiling heat transfer was introduced from many formulate proposed up to now by fitting the calculated values based on the formula to the observed data. The formula used for the predict on of the seed temperature was as follows: $\frac{{\partial}T_s}{\partial\theta}\;=\;-\frac{{\alpha}(T_s\;-\;T_L)^{3.3}}{W(C_s\;-\;\frac{{\delta}m(CT_s\;+\;{\sigma})}{T_s^2})}$ where C = difference of the specific heat between pure ice and water m=moisture content of soybean seed $T_s$ = seed temperature $T_L$ = Temperature of liquid nitrogen W = mass of soybean seed $\alpha$ = proportional constant $\delta$ = constant depends on variety or the type of seed $\theta$ = time $\sigma$ = latent heat of melting of pure ice This study will give important information in the hydro-freezing technique by liquid nitrogen, available as a new technique of processing agricultural products in the near future.

  • PDF

3D 스캔을 이용한 사이클 동작 전후 체표 변화 고찰 및 2D 전개 패턴의 비교 (2D Pattern Development of Body Surface from 3D Human Scan Data Using Standing and Cycling Postures)

  • 정연희;이예진
    • 한국생활과학회지
    • /
    • 제21권5호
    • /
    • pp.975-988
    • /
    • 2012
  • Although the pattern development for tight-fitting clothing has been carried out using 3D data on humans, the pattern development using 3D scan data obtained for various postures still remains an interesting subject. In this study, we have developed the 2D pattern using the 3D human body reflecting standing and cycling postures. The 3D scan data of a subject was obtained using Cyberware. 2C-AN program(Triangle simplification and the Runge-Kutta method) was used in the system to reduce the 3D scan data points and to make segmented triangular patches in a plane from 3D data. As results, surface distance and area of each body part of standing and cycling postures were also provided for the future application of the functional clothing construction. The area of center piece on the front (c.front) decreased by $106.45cm^2$(-13.08%) and that of lateral piece(s.back) on the back increased by $144.96cm^2$(18.69%) in the patterns of cycling posture. The girth of neck and waist for the cycling posture increased by 0.88cm (3.92%) and 1.56cm(4.40%) respectively, and the that of thigh decreased by 1.01cm(-2.24%). The differences between the area in the 2D pattern obtained from the 3D scan data and that in the 3D scan surface data for standing and cycling postures were very small($-10.34cm^2$(-0.32%) and $-44.33cm^2$(-1.32%)).

A Study on the Influence of Nonlinearity Coefficients in Air-Bearing Spindle Parametric Vibration

  • Chernopyatov, Y.A.;Lee, C.M.;Chung, W.J.;Dolotov, K.S.
    • International Journal of Precision Engineering and Manufacturing
    • /
    • 제6권1호
    • /
    • pp.51-58
    • /
    • 2005
  • The development of the high-efficiency machine-tools equipment and new cutting tool materials with high hardness, heat- and wear-resistance has opened the way to application of high-speed cutting process. The basic argument of using of high-speed cutting processes is the reduction of time and the respective increase of machining productivity. In this sense, the spindle units may be regarded as one of the most important units, directly affecting many parameters of high-speed machining efficiency. One of the possible types of spindle units for high-speed cutting is the air-bearing type. In this paper, we propose the mathematical model of the dynamic behavior of the air-bearing spindle. To provide the high-level of speed capacity and spindle rotation accuracy we need the adequate model of "spindle-bearings" system. This model should consider characteristics of the interactions between system components and environment. To find the working characteristics of spindle unit we should derive the equations of spindle axis movement under the affecting factors, and solve these equations together with equations which describe the behavior of lubricant layer in bearing (bearing stiffness equations). In this paper, the three influence coefficients are introduced, which describe the center of spindle mass displacement, angle of shaft rotation around the axes under the unit force application and that under the unit torque application. These coefficients are operated in the system of differential equations, which describes the spindle axis spatial movement. This system is solved by Runge-Kutta method. Obtained trajectories and amplitude-frequency characteristics were then compared to experimental ones. The analysis shows good agreement between theoretical and experimental results, which confirms that the proposed model of air-bearing spindle is correctis correct

인체의 3차원 스캔 데이터를 이용한 밀착 바디 슈트 개발 (2D Pattern Development of Tight-fitting Bodysuit from 3D Body Scan Data for Comfortable Pressure Sensation)

  • 정연희
    • 한국생활과학회지
    • /
    • 제15권3호
    • /
    • pp.481-490
    • /
    • 2006
  • Adjusting pressure level in the construction of athletes' tight-fitting garments by reducing the elastic knit pattern is a challenging subject, which influences the performance of the wearer directly. Therefore, in this study, relationship between the reduction rates of the basic pattern obtained from 3D human scan data and resultant clothing pressure was explored to improve the fit and pressure exerted by clothing. 3D scan data were obtained using Cyberware and they were transformed into a flat pattern using software based on Runge-Kutta method. Reduction rate was examined by subjective wear test as well as objective pressure measurement. As a result, difference in the length between the original 3D body scan data and the 2D tight-fitting pattern was 0.02$\sim$0.50cm (0.05$\sim$1.06%), which was within the range of tolerable limits in making clothes. Among the five garments, the 3T-pattern was superior in terms of subjective sensation and fit. The pressure of the 3T pattern was 2$\sim$4 gf/cm2 at five locations on the body, which is almost the same or a bit higher than that of Z-pattern. In the case of tight-fitting overall garment, the reduction rate of the pattern in the wale direction is more critical to the subjective sensation than the course direction. It is recommended that the reduction grading rules of course direction should be larger than that of Ziegert for a better fit of tight-fitting garments. In the case of wale direction, however, reduction grading rule should be kept the same as suggested earlier by Ziegert (1988).

  • PDF

심혈관계를 포함한 인공심장의 모델링 및 컴퓨터 시뮬레이션 (Modeling and Simulation of the Total Artificial Heart with Cardiovascular System)

  • 박준우;박성근;최종훈;조영호;최재순;안재목;민병구
    • 대한의용생체공학회:학술대회논문집
    • /
    • 대한의용생체공학회 1998년도 추계학술대회
    • /
    • pp.249-250
    • /
    • 1998
  • In this study, we modeled moving-actuator type Total Artificial Heart (TAH) with cardiovascular system as a form of electric circuit. The bronchial circulation, important for the imbalance between the left cardiac output and the right one, was considered and added to the model. In the model, the relations of hemodynamic variables, just as blood pressures, volumes, or flow rates of each part of body, can be expressed as simultaneous first order ordinary differential equations. To solve the equations by the numerical analysis, Runge-Kutta forth order approximation method was adopted. The simulation software (SimTAH), implemented in C++ as a window-based application program, was developed to display the hemodynamic variables and to receive control inputs from users. SimTAH was evaluated by comparison of the simulation results with the results of mock-circulation tests, in vitro.

  • PDF

Large deflection behavior and stability of slender bars under self weight

  • Goncalves, Paulo B.;Jurjo, Daniel Leonardo B.R.;Magluta, Carlos;Roitman, Ney;Pamplona, Djenane
    • Structural Engineering and Mechanics
    • /
    • 제24권6호
    • /
    • pp.709-725
    • /
    • 2006
  • In this paper the buckling and post-buckling behavior of slender bars under self-weight are studied. In order to study the post-buckling behavior of the bar, a geometrically exact formulation for the non-linear analysis of uni-directional structural elements is presented, considering arbitrary load distribution and boundary conditions. From this formulation one obtains a set of first-order coupled nonlinear equations which, together with the boundary conditions at the bar ends, form a two-point boundary value problem. This problem is solved by the simultaneous use of the Runge-Kutta integration scheme and the Newton-Raphson method. By virtue of a continuation algorithm, accurate solutions can be obtained for a variety of stability problems exhibiting either limit point or bifurcational-type buckling. Using this formulation, a detailed parametric analysis is conducted in order to study the buckling and post-buckling behavior of slender bars under self-weight, including the influence of boundary conditions on the stability and large deflection behavior of the bar. In order to evaluate the quality and accuracy of the results, an experimental analysis was conducted considering a clamped-free thin-walled metal bar. As this kind of structure presents a high index of slenderness, its answers could be affected by the introduction of conventional sensors. In this paper, an experimental methodology was developed, allowing the measurement of static or dynamic displacements without making contact with the structure, using digital image processing techniques. The proposed experimental procedure can be used to a wide class of problems involving large deflections and deformations. The experimental buckling and post-buckling behavior compared favorably with the theoretical and numerical results.

열차-교량의 동적 상호작용을 고려한 중·저속 및 초고속 자기부상열차와 가이드웨이의 동특성 해석 (Dynamic Interaction Analysis of Low, Medium and Super-high Speed Maglev and Guideways)

  • 민동주;정명락;이준석;김이현;김문영
    • 한국지진공학회논문집
    • /
    • 제15권3호
    • /
    • pp.1-9
    • /
    • 2011
  • 자기부상열차가 중 저속 및 초고속 주행 시 차량의 주행특성 및 교량의 동적 응답 결과를 제시하고자 한다. 수직 자유도 및 회전 자유도를 포함한 10자유도 자기부상열차에 대한 운동방정식을 유도하고, 모드 중첩법을 이용하여 교량의 운동방정식을 구성하였다. 또한 제어 방법으로 UTM01제어기법을 적용하여 수치해석을 수행하였다. 해석 예제로 노면조도, 가이드웨이의 처짐비, 차량 속도 등이 교량의 처짐과 차량의 부상공극 및 여러 가지 변수에 미치는 영향을 파악하였다. 부상공극은 조도의 조건에 따라 그 차이가 확연히 드러나고 또한 자기부상열차의 속도가 증가함에 따라 부상공극이 증가함을 알 수 있었다. 그리고 자기부상열차가 중, 저속 주행 시에는 교량에 대한 영향이 미비하지만 초고속 주행 시 교량에 대한 동적확대계수가 큰 값을 보여주었다.