• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.1
• /
• pp.36-45
• /
• 2015

Recently, a lot of studies have been made about the methods used to generate turbulent velocity fields stochastically in order to effectively predict broadband flow noise. Among them, the FRPM (Fast Random Particle Mesh) method which generates turbulence with specific statistical properties using turbulence kinetic energy and dissipation obtained from the steady solution of the RANS (Reynolds Averaged Navier-Stokes) equations has been successfully applied. However, the FRPM method cannot be applied to the flow noise problems involving intrinsic unsteady characteristics such as centrifugal fan. In this paper, to effectively predict the broadband noise generated by centrifugal fan, U-FRPM (unsteady FRPM) method is developed by extending the FRPM method to be combined with the unsteady numerical solutions of the unsteady RANS equations to generate the turbulence considered as broadband noise sources. Firstly, an unsteady flow field is obtained from the unsteady RANS equations through CFD (Computational Fluid Dynamics). Then, noise sources are generated using the U-FRPM method combined with acoustic analogy. Finally, the linear propagation model which is realized through BEM (Boundary Element Method) is combined with the generated sources to predict broadband noise at the listeners' position. The proposed technique is validated to compare its prediction result with the measured data.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.4
• /
• pp.351-363
• /
• 2020

Launch vehicles are subject to airborne acoustic loads during atmospheric flight and these effects become pronounced especially in transonic region. As the vibration due to the acoustic loads can cause malfunction of payloads, it is essential to predict and reduce the acoustic loads. In this study, a complete process has been developed for predicting airborne vibro-acoustic environment inside the payload pairing and subsequent noise reduction procedure employing acoustic blankets and Helmholtz resonators. Acoustic loads were predicted by Reynolds-Averaged Navier-Stokes (RANS) analysis and a semi-empirical model for pressure fluctuation inside turbulent boundary layer. Coupled vibro-acoustic analysis was performed using VA One SEA's Finite Element Statistical Energy Analysis (FE-SEA) hybrid module and ANSYS APDL. The process has been applied to a hammerhead launch vehicle to evaluate the effect of acoustic load reduction and accordingly to verify the effectiveness of the process. The presently developed process enables to obtain quick analysis result with reasonable accuracy and thus is expected to be useful in the initial design phase of a launch vehicle.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.10
• /
• pp.618-627
• /
• 2016

This study numerically analyzes the characteristics of the velocity distribution for each location of a square-sectional $180^{\circ}$ bent duct using a Reynolds Stress Turbulent model. The flow parameters were varied, including the working fluids, inlet velocity, surface roughness, radius of curvature, and hydraulic diameter. The boundary conditions for computational fluid dynamics analysis were inlet temperatures of air and water of 288 K and 293 K, inlet air velocity of 3-15 m/s, inner surface roughness of 0-0.001 mm, radius of curvature of 2.5-4.5 D, and hydraulic diameter of 70-100 mm. The working fluid characteristics were highly affected by changes in the viscous force. The maximum velocity profiles in the bent duct were indicated when the $90^{\circ}$ section was in the region of X/D=0.8 and the $180^{\circ}$ section was in the region of Y/D=0.8. Lower surface roughness and higher radius of curvature resulted in a higher rate of velocity change. Also, an efficient measuring location downstream of the bent duct is suggested since the flow deviations were the most stable when the straight duct length was in the region of L/D=30. The minimum deviations at the same velocity conditions according to the hydraulic diameter were mostly indicated in the range of L/D=15-30 based on the standard deviation characteristics.

• Journal of the Korean Institute of Gas
• /
• v.2 no.1
• /
• pp.28-39
• /
• 1998

This study investigated the natural convection and radiation in a rectangular enclosure with ceiling vent experimentally and numerically. A heat source is located on the center of the bottom surface. The analysis was peformed a pure convection and is combination of natural convection and radiation. The shape of the considered two dimensional model is a square whose center of ceiling($30\%$) is opened. The numerical simulations are carried out for the pure natural convection case and the combined heat transfer case by using the SIMPLE algorithm. For the turbulent flow, Reynolds stresses are closed by the standard $k-{\epsilon}$ model and the wall function is used to determine the wall boundary conditions. The experiment was performed on the same geometrical shape as the computations. The radiative heat transfer is analized by the S-N discrete ordinates method. The results of pure natural convection are compared with those of combined heat transfer by the velocity vectors, stream lines, isothermal lines. The results obtained are as follows 1. Comparing the results of pure convection with those of the combined convection-radiation through the shape of stream lines, isothermal lines are similar to each other. 2. The temperature fields obtained by numerical method are compared to those obtained by experimental one, and it is found that they are showed mean relative error $8.5\%$. 3. Visualization bt smoke is similar to computational results.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.20-26
• /
• 2019

Where hydraulic structures are to be installed over the entire width of a river or stream, usually a bed protection structure is to be installed. However, a local scour occurs in which the river bed downstream of the river protection system is eroded due to the influence of the upstream flow characteristics. This local scour is dominant in the flow and turbulence characteristics at the boundary of the flow direction and in the material of the bed materials, and may gradually become dangerous over time. Therefore, in this study, we compared the turbulence patterns in the local scour hole at the downstream of the river bed protection with the results of the analysis of the mobile bed experiment, and compared with the application of OpenFoam, a three dimensional numerical analysis model. The distribution of depth-averaged relative turbulence intensities along the flow direction was analyzed. In addition to this result, the stabilization of scour hole was compared with the bed shear stress and Shields parameter, and the results were compared by changing the initial turbulent flow conditions. From the results, it was confirmed that the maximum depth of generation of the three-stage was dominantly developed by the magnitude of depth-averaged relative turbulence intensity rather than the mean flow velocity. This result also suggests that design, construction or gate control are needed to control the depth-averaged relative turbulence intensities in order to reduce or prevent the local scour faults that may occur in the downstream part of the bed protection.

• Journal of the Korean earth science society
• /
• v.23 no.4
• /
• pp.357-368
• /
• 2002

In this study we have tried to explain the barrel-shaped morphology for young supernova remnants considering the dynamical effects of the ejecta. We consider the magnetic field amplification resulting from the Rayleigh-Taylor instability near the contact discontinuity. We can generate the synthetic radio image assuming the cosmic-ray pressure and calculate the azimuthal intensity ratio (A) to enable a quantitative comparison with observations. The postshock magnetic field are amplified by shearing, stretching, and compressing at the R-T finger boundary. The evolution of the instability strongly depends on the deceleration of the ejecta and the evolutionary stage of the remnant. the strength of the magnetic field increases in the initial phase and decreases after the reverse shock passes the constant density region of the ejecta. However, some memory of the earlier phases of amplification is retained in the interior even when the outer regions turn into a blast wave. The ratio of the averaged magnetic field strength at the equator to the one at the pole in the turbulent region can amount to 7.5 at the peak. The magnetic field amplification can make the large azimuthal intensity ratio (A=15). The magnitude of the amplification is sensitive to numerical resolution. This mens the magnetic field amplification can explain the barrel-shaped morphology of young supernova remnant without the dependence of the efficiency of the cosmic-ray acceleration on the magnetic field configuration. In order for this mechanism to be effective, the surrounding magnetic field must be well-ordered. The small number of barrel-shaped remnants may indicate that this condition rarely occurs.