• Journal of Ocean Engineering and Technology
• /
• v.23 no.6
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.181-188
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.6
• /
• pp.34-40
• /
• 2010

Bragg reflection of nonlinear waves is simulated by a 2D fully nonlinear numerical wave tank (NWT). The developed NWT was based on the Boundary Element Method (BEM) with potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme with Runge-Kutta 4th-order time integration. A spatial variation of wave elevations and their Fourier amplitudes of each component are compared to investigate the effect of sea bottom ripples and their relative heights. The incident waves over an undulated sea bottom are partially reflected and changed to partial standing waves due to Bragg reflection. The present results are verified with linear calculations and experimental data. It is found that the 1st-order wave component is mainly affected by Bragg reflection and its spatial modulation is significant in front of the bottom ripples.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.127-132
• /
• 2003

A study of high-pressure n-heptane droplet vaporization is conducted with emphasis placed on equilibrium at vapor-liquid interface. General frame of previous rigorous model[1] is retained but tailored for flash equilibrium calculation of vapor-liquid interfacial thermodynamics. The model is based on complete time-dependent conservation equations with a full account of variable properties and vapor-liquid interfacial thermodynamics. The influences of high-pressure phenomena, including ambient gas solubility, thermodynamic non-ideality, and property variation on the droplet evaporation are investigated. The governing equations and associated moving interfacial boundary conditions are solved numerically using a implicit scheme with the preconditioning method and the dual time integration technique. And a parametric study of entire droplet vaporization history as a function of ambient pressure, temperature has been conducted. Some computational results are compared with Sato's experimental data for the validation of calculations. For low ambient temperatures, the droplet lifetime first increases with pressures, then decreases for high pressures. For higher ambient temperatures, the droplet lifetime increase with less amplitude than that of low ambient temperatures, which then decreases with more amplitude than that of low temperatures. The solubility of nitrogen can not be neglected in the high pressure and it becomes higher as the pressure goes up.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.405-412
• /
• 2012

A flow control on airfoil installed a vortex cell for high efficiency wind turbine blade in stationary and dynamic stall conditions have been numerically investigated by solving the compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with Roe's flux-difference splitting and an implicit time-integration method coupled with dual time step sub-iteration. The computed result for the airfoil in the stationary showed that lift-drag ratio increases due to low pressure by the vortex cell. The oscillating airfoil with the vortex cell showed that the magnitude of hysteresis loop is reduced due to the enhanced vortex in the cell.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.375-382
• /
• 2008

The particular integral formulation for two(2D) and three(3D) dimensional inelastic transient dynamic stress analysis is presented. The elastostatic equation is used for the complementary solution. Using the concept of global shape function, the particular integrals for displacement and traction rates are obtained to approximate acceleration of the inhomogeneous equation. The Houbolt time integration scheme is used for the time-marching process. The Newton-Raphson algorithm for plastic multiplier is used to solve the system equation. Numerical results of four example problems are given to demonstrate the validity and accuracy of the present formulation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.75-84
• /
• 2002

A numerical analysis of unsteady motion in solid rocket motors with a nozzle has been conducted. The numerical formulation including modified $\kappa$-$\varepsilon$ turbulence model treats the complete conservation equation for the gas phase and the one-dimensional equations in the radial direction for the condensed phase. A fully coupled implicit scheme based on a dual time-stepping integration algorithm has been adopted to solve the governing equations. After obtaining a steady state solution, pulse and periodic oscillations of pressure are imposed at the head-end to simulate acoustic oscillations of a travelling-wave motion in the combustion chamber. Various steady and unsteady state features in the combustion chamber of a rocket motor has been analyzed as results of numerical calculations.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• Advances in Energy Research
• /
• v.3 no.2
• /
• pp.81-95
• /
• 2015

Microgrid, which can be considered as an integration of various dispersed resources (DRs), is characterized by number of DRs interfaced through the power electronics converters. The microgrid comprising these DRs is often operated in an islanded mode. To minimize the cost, reduce complexity and increase reliability, it is preferred to avoid any communication channel between them. Consequently, the droop control method is traditionally adopted to distribute active and reactive power among the DRs operating in parallel. However, the accuracy of distribution of active and reactive power among the DRs controlled by the conventional droop control approach is highly dependent on the value of line impedance, R/X i.e., resistance to reactance ratio of the line, voltage setting of inverters etc. The limitations of the conventional droop control approach are demonstrated and a modified droop control approach to reduce the effect of impedance mis-match and improve the time response is proposed. The error in reactive power sharing is minimized by inserting virtual impedance in line with the inverters to remove the mis-match in impedance. The improved time response is achieved by modifying the real-power frequency droop using arctan function. Simulations results are presented to validate the effectiveness of the control approach.