• 한국해양공학회지
• /
• 제34권6호
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Structural Engineering and Mechanics
• /
• 제47권5호
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Structural Engineering and Mechanics
• /
• 제90권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.

• 한국공간구조학회논문집
• /
• 제20권2호
• /
• pp.59-68
• /
• 2020

In this paper, the dynamic snapping of the 3-free-nodes spatial truss model was studied. A governing equation was derived considering geometric nonlinearity, and a model with various conditions was analyzed using the fourth order Runge-Kutta method. The dynamic buckling phenomenon was observed in consideration of sensitive changes to the force mode and the initial condition. In addition, the critical load level was analyzed. According to the results of the study, the level of critical buckling load elevated when the shape parameter was high. Parallelly, the same result was caused by the damping term. The sensitive asymmetrical changes showed complex orbits in the phase space, and the critical load level was also becoming lowly. In addition, as the value of damping constant was high, the level of critical load also increases. In particular, the larger the damping constant, the faster it converges to the equilibrium point, and the occurrence of snapping was suppressed.

• Journal of Mechanical Science and Technology
• /
• 제14권6호
• /
• pp.690-698
• /
• 2000

The laminar boundary layer flow and heat transfer of anisotropic fluids in the vicinity of a wedge have been examined with constant surface temperature. The similarity variables found by Falkner and Skan are employed to reduce the stream wise-dependence in the coupled nonlinear boundary layer equations. The numerical solutions are presented using the fourth-order Runge - Kutta method and the distribution of velocity, micro-rotation, shear and couple stresses and temperature across the boundary layer are plotted. These results are also compared with the corresponding flow problems for Newtonian fluid over wedges. It is found that for a constant wedge angle, the skin friction coefficient is lower for micropolar fluid, as compared to Newtonian fluid. For the case of the constant material parameter K, however, the magnitude of velocity for anisotropic fluid is greater than that of Newtonian fluid. The numerical results also show that for a constant wedge angle with a given Prandtl number, Pr = I, the effect of increasing values of K results in increasing thermal boundary layer thickness for anisotropic fluid, as compared with Newtonian fluid. For the case of the constant material parameter K, however, the heat transfer rate for anisotropic fluid is lower than that of Newtonian fluid.

• 대한기계학회논문집B
• /
• 제22권11호
• /
• pp.1635-1646
• /
• 1998

The present study investigated the effects of geometric parameters considered on the performance of an automotive scroll compressor by assuming ideal, semi-real and real gases for R-134a. The geometric parameters were center thickness of scroll, height of scroll and the size of discharge port. Fourth-order Runge-Kutta method was applied to solve the thermodynamic equations, leaking rate equation and the equation of motion of discharge valve for ideal, semi-real, and real gases. The volumetric and adiabatic efficiencies for semi-real and real gases differed little, but those for ideal gas differed by 18% and 25% compared with those for real gas at 2,000rpm. The volumetric efficiency changed little as the design angle of scroll (${\gamma}$) changed, but the adiabatic efficiency at ${\gamma}$ of $34^{\circ}$ was higher by 2.4% than that of $147^{\circ}$ for 2,000rpm. The volumetric and adiabatic efficiencies at scroll height of 29.8mm were higher by 1.7% and 2.8% than those of 65.8mm. The volumetric efficiency changed little as the size of discharge port changed, but the adiabatic efficiency increased a little as the size of discharge port decreased.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2009년도 춘계학술대회 논문집
• /
• pp.334-336
• /
• 2009

An elastic-plasticity model during the austenitic decomposition was derived and implemented to incorporate the two important deformation behaviors observed during the phase transformations: the volumetric strain and transformation induced plasticity due to the temperature change and phase transformation. To obtain transformed phase volume fractions during cooling, the fourth order Runge-Kutta method was used to solve the Kirkaldy's phase kinetics model which is function of temperature, austenitic grain size and chemical composition. The volumetric strain was calculated by considering the densities of constituent phases, while the transformation induced plasticity was based on the micro-plasticity due to the volume mismatch between soft austenitic phase and other harder phases. The constitutive equations were implemented into the implicit finite element software and a simple boundary value problem was chosen as a model problem to validate the effect of transformation plasticity on the deformation behavior of steel under cooling from high temperature. It was preliminary concluded that the transformation plasticity plays a critical role in relaxing the developed stress during forming and thus reducing the magnitude of springback.

• 한국항해학회지
• /
• 제4권1호
• /
• pp.31-50
• /
• 1980

It does not seem necessarily practicable to keep the system always in optimal condition, athough the control system of the follow-up mechanism on the most marine gyro compasses is to be adjusted by the operator through the gain adjustment. Sometimes a sustained oscillation or an incorrect gyro reading occurs to the system. For such a system any systematical research or theoretical basis of the guide for the optimal gain adjustment has not been reported yet. As a basic investigation of the theoretical system analysis to solve the problems concerned, the author attempts in this paper to express the system in a mathematical model deduced from the results of the theoretical approach and the experimental observation of each element contained in the follow-up mechanism of Hokshin D-1 gyro compass, and to constitute an over-all closed loop transfer function. This funciton being reverted to a fourth orderlinear differential equation, the first order simultaneous differential equations are obtained by means of the state-variables. The latter equations are solved by the Runge-Kutta method with digital computer. By comparing the characteristic of the simulated over-all output with that of the experimental result, it is shown that both outputs are nearly consistent with each other. It is also expected that the system representation proposed by this paper is valid and will be a prospective means in a further study on the design and optimal adjustment of the system.

• Structural Engineering and Mechanics
• /
• 제70권5호
• /
• pp.623-637
• /
• 2019

This paper analyzes the non-stationary vibration and super-harmonic resonances in nonlinear dynamic motion of viscoelastic nano-resonators. For this purpose, a new coupled size-dependent model is developed for a plate-shape nano-resonator made of nonlinear viscoelastic material based on modified coupled stress theory. The virtual work induced by viscous forces obtained in the framework of the Leaderman integral for the size-independent and size-dependent stress tensors. With incorporating the size-dependent potential energy, kinetic energy, and an external excitation force work based on Hamilton's principle, the viscous work equation is balanced. The resulting size-dependent viscoelastically coupled equations are solved using the expansion theory, Galerkin method and the fourth-order Runge-Kutta technique. The Hilbert-Huang transform is performed to examine the effects of the viscoelastic parameter and initial excitation values on the nanosystem free vibration. Furthermore, the secondary resonance due to the super-harmonic motions are examined in the form of frequency response, force response, Poincare map, phase portrait and fast Fourier transforms. The results show that the vibration of viscoelastic nanosystem is non-stationary at higher excitation values unlike the elastic ones. In addition, ignoring the small-size effects shifts the secondary resonance, significantly.

• 상하수도학회지
• /
• 제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.