• Journal of The Korean Astronomical Society
• /
• v.47 no.3
• /
• pp.87-98
• /
• 2014

We develop a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is implemented based on the smoothed particle hydrodynamics. The initial displacement and velocity of simulation particles are calculated according to second-order Lagrangian perturbation theory using the power spectra of dark matter and baryonic matter. The initial background temperature is given by Recfast and the temperature uctuations at the initial particle position are assigned according to the adiabatic model. We use a time-limiter scheme over the individual time steps to capture shock-fronts and to ease the time-step tension between the shock and preshock particles. We also include the astrophysical gas processes of radiative heating/cooling, star formation, metal enrichment, and supernova feedback. We test the code in several standard cases such as one-dimensional Riemann problems, Kelvin-Helmholtz, and Sedov blast wave instability. Star formation on the galactic disk is investigated to check whether the Schmidt-Kennicutt relation is properly recovered. We also study global star formation history at different simulation resolutions and compare them with observations.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.1
• /
• pp.100-113
• /
• 2017

Snap rolling during hard turning and instability during emergency rising are important features of submarine operation. Hydrodynamics modeling using a high incidence flow angle is required to predict these phenomena. In the present study, a quasi-steady dynamics model of a submarine suitable for high-incidence-angle maneuvering applications is developed. To determine the hydrodynamic coefficients of the model, static tests, dynamic tests, and control surface tests were conducted in a towing tank and wind tunnel. The towing tank test is conducted utilizing a Reynolds number of $3.12{\times}10^6$, and the wind tunnel test is performed utilizing a Reynolds number of $5.11{\times}10^6$. In addition, least squares, golden section search, and surface fitting using polynomial models were used to analyze the experimental results. The obtained coefficients are presented in tabular form and can be used for various purposes such as hard turning simulation, emergency rising simulation, and controller design.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.644-653
• /
• 1991

A wave instability problem is formulated for natural convection flows adjacent to a inclined isothermal surface in pure water near the density extremum. It accounts for the nonparallelism of the basic flow and temperature fields. Numerical solutions of the hydrodynamic stability equations constitute a two-point boundary value problem which are accurately solved using a computer code COLSYS. Neutral stability results for Prandtl number of 11.6 are obtained for various angles of inclination of a surface in the range from-10 to 30 deg. The neutral stability curves are systematically shifted toward modified Grashof number G=0 as one proceeds from downward-facing inclined plate(.gamma.<0.deg.) to upward-facing inclined plate (.gamma.>0.deg.). Namely, an increase in the positive angle of inclination always cause the flows to be significantly more unstable. The present results are compared with the results for the parallel flow model. The nonparallel flow model has, in general, a higher critical Grashof number than does the parallel flow model. But the neutral stability curves retain their characteristic shapes.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.679-682
• /
• 2002

Energy dissipations in a general PHE flow are the compounded effects of the piled corrugate geometries and its wall pressure and temperature distributions. In addition, although the exchangers are substantial pieces of engineering equipment, they are composed of a very large number of nominally identical and small geometrical elements. In the present numerical study, the three-dimensionally complicated energy dissipation fields and those wall-shape-induced flow destabilization are investigated in the cross-corrugated passages, which result in high energy transports with comparatively low pressure drop. We revealed the critical conditions as $Re=157.3 for the wall-shape-induced flow destabilization in a general PHE element by initial value method, or shooting method, and compare its value to that of analytical solution of plane Poiseille flow, two-dimensional grooved flow and so on. We also observed the detailed variation of flow field and energy transportation with changes in time and flow variables such as Reynolds number. Lastly, we considered the flow natural frequency, or Strouhal number, with variation of hydrodynamic conditions for the best use of active control, such as forced mass flow rate pulsative flow, to enhance energy transportation.

• Journal of the KSME
• /
• v.16 no.2
• /
• pp.84-91
• /
• 1976

An experimental investigation of the natural circulating boiling flow characteristic in three cases of annulus with different outer diameter, and the effect of annular gap size on the burnout behavior is presented. The experimental work was conducted for each case of test section at system pressure of $1kg/cm^2$ and inlet subcooling $0-20^{\circ}C$ in the full range of throttling ratio. As the result, the following facts were found. 1) With the increase of ${\Delta}T_{sub}$, $D_{2}$ and A/A_{o}$, $q_{BO}$ increases on the whole, and with the decrease of ${\Delta}T_{sub}$ and $D_{2}$, hydrodynamic instability is accelerated to happen prematually. 2)With the increase of ${\Delta}T_{sub}$, $D_{2}$ and A/A_{o}$ burnout characteristic shows the high velocitylow quality burnout, and with the decrease, low-velocity-high quality burnout. 3)With the decrease of A/A_{o}$, hyddrodynamic instability is singnificantly restrained and the difference of $q_{BO}$ in each $D_2$ under same condition is gradually reduced, finally converging into $1.9{\times}10^{5}kcal/m^{2}-hr$.

• Kyungpook Mathematical Journal
• /
• v.47 no.4
• /
• pp.455-471
• /
• 2007

The magnetohydrodynamic axisymmetric instability of a streaming gas jet surrounded by bounded fluid with non-uniform field has been developed. The problem is formulated, solved and the boundary conditions are applied across the interfaces. The eigenvalue relation is derived and discussed analytically and the results are confirmed numerically. Some reported works are recovered as limiting cases from the present general results. The streaming has a destabilizing effect for all short and long wavelengths. The capillary force is stabilizing for short wavelengths but it is destabilizing for long wavelengths. The axial magnetic fields interior the gas and fluid media are stabilizing. The transverse field is destabilizing for all wavelengths. The radii ratio of the gas and fluid cylinders plays an important role for stabilizing the model and made it more realistic one than the full liquid jet or/and the ordinary hollow jet. The numerical analysis clarify the stable and unstable domains based on different values of the various parameters of the problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.4
• /
• pp.399-408
• /
• 2010

The characteristics of lifted laminar propane flames diluted with nitrogen have been investigated experimentally to elucidate self-excitation and the effects of flame curvature. Flame oscillation modes are classified as follows: oscillation induced by heat loss, a combination of oscillations induced by heat loss and buoyancy, and a combination of the oscillations induced by heat loss and diffusive thermal instability. It is shown that the oscillation induced only by heat loss is not relevant to the diffusive thermal instability and hydrodynamic instability caused by buoyancy; this oscillation is observed under all lift-off flame conditions irrespective of the fuel Lewis number. These experimental evidences are displayed through the analysis of the power spectrum for the temporal variation of lift-off height. The possible mechanism of the oscillation induced by heat loss is also discussed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.1
• /
• pp.18-25
• /
• 1984

The transition boiling of a liquid droplet on a heated flat surface was studied utilizing Kotake's model with the effects of viscosity of a thin vapor layer between the droplet and the hot plate taken into account. This problem was analyzed considering the process of the droplet evaporation which resulted in hydrodynamic instability at the liquid-vapor interface. The results of the study are as follows; (1)The effect of the viscosity in the vapor layer at the interface appears as a dimensionless number N, namely .sigma. .delta.$_{0}$ /.rho.nu.$^{2}$ (2)The time required for evaporation at the transitional region increases with the temperature difference ratio .DELTA. T$_{r}$. The rate of increase of the total evaporation time becomes larger as increasing of N$_{m}$(N number at maximum heat flux) increases.s.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.6 no.4
• /
• pp.89-99
• /
• 2003

Smoothed particle hydrodynamics, SPH, is a gridless Lagrangian technique which is a useful alternative numerical analysis method to simulate high velocity deformation problems as well as astrophysical and cosmological problems. The SPH method brings about some difficulties such as tensile Instability and stress oscillation. A new SPH method, so called normalized algorithm, was introduced to overcome these difficulties. In this paper we aimed to estimate this method and have developed an one-dimensional normalized SPH program. The high velocity impact model of an aluminum bar has been analysed by using the developed program and a commercial hydrocode, LS-DYNA. The obtained numerical results showed good agreement with the results of the same model in reference. The program also showed more stable results than those of LS-DYNA in stress oscillation. We hopefully expect that the developed one-dimensional normalized SPH program can be used to solve hydrodynamic problems especially for explosive detonation analysis.