• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.13 no.6
• /
• pp.88-98
• /
• 2014

A moving mass on a skewed linear guideway model to analyze the friction-induced stick-slip behavior of ball-screw-driven slides is proposed. To describe the friction force, a friction coefficient function is modelled as a third-order polynomial of the relative velocity between the slide mass and a guideway. A nonlinear differential equation of motion is derived and an approximate solution is obtained using a perturbation method for the amplitudes and base frequencies of both pure-slip and stick-slip oscillations. The results are presented with time responses, phase plots, and amplitude plots, which are compared adequately with those obtained by Runge Kutta 4th-order numerical integration, as long as the difference between the static and kinematic friction coefficients is small. However, errors in the results by the approximate solution increase and are not negligible if the difference between the friction coefficients exceeds approximately 40% of the static friction coefficient.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.234-239
• /
• 2002

Acoustic sounds generated by a uniform flow around a two-dimensional circular cylinder at Re=200 are simulated using finite difference lattice Boltzmann method. A third-order-accurate up-wind scheme is used for spartial derivatives, and a second-order-accurate Runge-Kutta scheme is used for time marching. The results show that in capturing very small acoustic pressure fluctuation with same frequency of Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of acoustic sound is presented that acoustic which approaches tire upstream due to Doppler effect in the uniform flow slowly propagates, while that for the downstream quickly propagates. It is also apparent that the size of sound pressure is proportional for central distance $r^{-1/2}$ of the cylinder.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3653-3664
• /
• 2021

In PWR three-dimensional transient coupling calculation software CORCA-K, the nodal Green's function method and diagonal implicit Runge Kutta method are used to solve the spatiotemporal neutron dynamic diffusion equation, and the single-phase closed channel model and one-dimensional cylindrical heat conduction transient model are used to calculate the coolant temperature and fuel temperature. The LMW, NEACRP and PWR MOX/UO₂ benchmarks and FangJiaShan (FJS) nuclear power plant (NPP) transient control rod move cases are used to verify the CORCA-K. The effects of burnup, fuel effective temperature and ejection rate on the control rod ejection process of PWR are analyzed. The conclusions are as follows: (1) core relative power and fuel Doppler temperature are in good agreement with the results of benchmark and ADPRES, and the deviation between with the reference results is within 3.0% in LMW and NEACRP benchmarks; 2) the variation trend of FJS NPP core transient parameters is consistent with the results of SMART and ADPRES. And the core relative power is in better agreement with the SMART when weighting coefficient is 0.7. Compared with SMART, the maximum deviation is -5.08% in the rod ejection condition and while -5.09% in the control rod complex movement condition.

• Advances in nano research
• /
• v.10 no.5
• /
• pp.481-491
• /
• 2021

Graphene Nanosheets play an important role in nanosensors due to their proper surface to volume ratio. Therefore, the main purpose of this paper is to consider the nonlinear vibration behavior of graphene nanosheets (GSs) under the influence of electromagnetic fields and electrical current create forces. Considering more realistic assumptions, new equations have been proposed to study the nonlinear vibration behavior of the GSs carrying electrical current and placed in magnetic field. For this purpose, considering the influences of the magnetic tractions created by electrical and eddy currents, new relationships for electromagnetic interaction forces with these nanosheets have been proposed. Nonlinear coupled equations are discretized by Galerkin method, and then solved via Runge-Kutta method. The effect of different parameters such as size effect, electrical current magnitude and magnetic field intensity on the vibration characteristics of GSs is investigated. The results show that the magnetic field increases the linear natural frequency, and decreases the nonlinear natural frequency of the GSs. Excessive increase of the magnetic field causes instability in the GSs.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.19 no.11
• /
• pp.49-57
• /
• 2020

The forced vibration characteristics for two impulse forces with time lag was discussed in the vehicle dull progress model. Detailed numerical analyses of the time domain were performed systematically. By the two exciting impulse forces, the responses of displacement, the velocity, and the acceleration were investigated in detail for the vehicle's vibration. Notably, the forced vibration responses in the time domain can be used to identify and monitor several vehicle vibration models. In order to define the responses of displacement, the velocity, and the acceleration, we applied a numerical technique (i.e., the Runge-Kutta-Gill method[1,2]). These variables were subsequently used to analyze the vehicle's vibration according to the time lapse and while it passed over a bump stock; moreover, the characteristics of the variables were analyzed in detail according to their force conditions. Finally, the intrinsic characteristics of the forced vibration were discussed in the context of the automobile model. Overall, our results indicate that the tested method can be successfully applied under different damped conditions.

• Structural Engineering and Mechanics
• /
• v.90 no.4
• /
• pp.357-370
• /
• 2024

Due to the fact that the mechanism of the effects of temperature and initial geometric imperfection on low-velocity impact problem of axially moving plates is not yet clear, the present paper is to fill the gap. In the present paper, the nonlinear dynamic behavior of axially moving imperfect graphene platelet reinforced metal foams (GPLRMF) plates subjected to lowvelocity impact in thermal environment is analyzed. The equivalent physical parameters of GPLRMF plates are estimated based on the Halpin-Tsai equation and the mixing rule. Combining Kirchhoff plate theory and the modified nonlinear Hertz contact theory, the nonlinear governing equations of GPLRMF plates are derived. Under the condition of simply supported boundary, the nonlinear control equation is discretized with the help of Gallekin method. The correctness of the proposed model is verified by comparison with the existing results. Finally, the time history curves of contact force and transverse center displacement are obtained by using the fourth order Runge-Kutta method. Through detailed parameter research, the effects of graphene platelet (GPL) distribution mode, foam distribution mode, GPL weight fraction, foam coefficient, axial moving speed, prestressing force, temperature changes, damping coefficient, initial geometric defect, radius and initial velocity of the impactor on the nonlinear impact problem are explored. The results indicate that temperature changes and initial geometric imperfections have significant impacts.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.372-378
• /
• 1995

This study describes an analysis on the sinking characteristics in purse seine. The experiment was carried out using three simplified model seines in a flume tank under still water condition. The densities($\rho$) of netting materials were 0.91g/$cm^3$ for PP seine, 1.14g/$cm^3$ for PA seine and 1.38g/$cm^3$ for PES seine. Differential equations were derived from the conservation of momenta of the model seines and used to determine the sinking speeds of the depths of leadline and the other portions of seines. An analysis carried out by simultaneous differential equations for numerical method by subroutine Runge-Kutta-Gill. The results obtained were as follows: 1) Sinking speed of net margin was fastest for PP seine, followed by PA and PES seines. 2) The coefficient of resistance for netting of seine was estimated to be $K_D=0.061({\frac{\rho}{{\rho}_w}})^4$. 3) The coefficient of resistance for netting bundle of seine was estimated to be $C_R=0.91({\frac{\rho}{{\rho}_w}}$. 4) In all seines, the calculated depths of leadline closely agreed with the measured ones(meas.=0.99cal.).

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.34 no.3
• /
• pp.274-282
• /
• 1998

This study deals with sinking resistance for purse seine, in the case of different d/l, ratio of twine diameter and leg length. Experiments carried out on the six types simplified reduced model seines which were made of knotless netting. The nettings were woven in different leg length 4.3, 5.0, 5.5, 6.0, 6.6 and 7.7mm of polyester 28 tex two threads two -ply twine, and each of the seines were named I, II, III, IV, V and Ⅵ seine. Dimension of seine models were 450cm for corkline and 85cm for seine depth, each seines rigged up 160g of float for a floatline and 50g (underwater weight) of lead for a leadline. Experiments were measured in the observation channel of a flume tank at the static conditions Sinking motion was recorded by the two sets TV-camera for VTR which were placed in top and side of the model seine, and reading coordinate carried out by the video digitization system. An analysis were calculated out by simultaneous differential equations for numerical method by Runge - Kutta - Gill sub - routine. The results obtained were as follows: 1. Average sinking speed of seine of seine margin was fastest for Ⅵ seine followed by V, IV, III, II and I seines. 2. The coefficient of resistance for a seine wall was depended upon the ratio of d/l : KD =0.081 (d/l )-0.5 3. The coefficient of resistance for netting bundle was not depended upon the ratio of d/l :CR = 0.91 (), d : Twine diameter, l : Leg length, : Density of netting materals, $\omega$ : Density of water

• Membrane Journal
• /
• v.17 no.4
• /
• pp.352-358
• /
• 2007

A nonporous hollow fiber membrane contactor was used to control the concentration of oxygen dissolved in an aqueous solution, which was predicted along the hollow fiber membrane using a computer simulation. The governing ordinary differential equations were derived for the occurrent flows of the feed aqueous solution and the feed gas mixture in a membrane contactor and they were numerically solved using the 5th Runge-Kutta-Verner method with a personal computer, where the program was coded utilizing a software of the Compaq Visual Fortran 6.6. It is found that the concentration of oxygen dissolved in water increases from 30 to 64 ppm as the length of the hollow fiber increases from 0.4 to 1.2 m when the membrane of fibers are equal to be 16,000; the flow rate of the feed gas is kept to be 0.536 mol/sec; its pressure is maintained to be 486 kPa; the flow rate of the water is 16.69 mol/sec. As the flow rate of the water increases from 9.26 to 26.85 mol/sec, the concentration of oxygen decreases from 40 to 20 ppm with the constant fiber length of 0.4 m. Finally, it is observed that the concentration of oxygen increases from 33 to 69 ppm as the pressure of the feed gas increases from 298 to 847 kPa.