• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.61-97
• /
• 2016

This is the second of two papers on the 3D numerical modeling of nearshore hydro- and morphodynamics. In Part I, the focus was on surf and swash zone hydrodynamics in the cross-shore and longshore directions. Here, we consider nearshore processes with an emphasis on the effects of oceanic forcing and beach characteristics on sediment transport in the cross- and longshore directions, as well as on foreshore bathymetry changes. The Delft3D and XBeach models were used with four turbulence closures (viz., ${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES) to solve the 3D Navier-Stokes equations for incompressible flow as well as the beach morphology. The sediment transport module simulates both bed load and suspended load transport of non-cohesive sediments. Twenty sets of numerical experiments combining nine control parameters under a range of bed characteristics and incident wave and tidal conditions were simulated. For each case, the general morphological response in shore-normal and shore-parallel directions was presented. Numerical results showed that the ${\kappa}-{\varepsilon}$ and H-LES closure models yield similar results that are in better agreement with existing morphodynamic observations than the results of the other turbulence models. The simulations showed that wave forcing drives a sediment circulation pattern that results in bar and berm formation. However, together with wave forcing, tides modulate the predicted nearshore sediment dynamics. The combination of tides and wave action has a notable effect on longshore suspended sediment transport fluxes, relative to wave action alone. The model's ability to predict sediment transport under propagation of obliquely incident wave conditions underscores its potential for understanding the evolution of beach morphology at field scale. For example, the results of the model confirmed that the wave characteristics have a considerable effect on the cumulative erosion/deposition, cross-shore distribution of longshore sediment transport and transport rate across and along the beach face. In addition, for the same type of oceanic forcing, the beach morphology exhibits different erosive characteristics depending on grain size (e.g., foreshore profile evolution is erosive or accretive on fine or coarse sand beaches, respectively). Decreasing wave height increases the proportion of onshore to offshore fluxes, almost reaching a neutral net balance. The sediment movement increases with wave height, which is the dominant factor controlling the beach face shape.

• Proceedings of the Acoustical Society of Korea Conference
• /
• spring
• /
• pp.469-474
• /
• 2002

In recent years, modularization of engine parts has increased the application of plastic products in air intake systems. Plastic intake manifolds provide many advantages including reduced weight, contracted cost, and lower intake air temperatures. These manifolds, however, have some weakness when compared with customary aluminium intake manifolds, in that they have low sound transmission loss because of their lower material density. This low transmission loss of plastic intake manifolds causes several problems related to flow noise, especially when the throttle is opened quickly. The physical processes, responsible for this flow noise, include turbulent fluid motion and relative motion of the throttle to the airflow. The former is generated by high-speed airflow in the splits between the throttle valve and the inner-surface of the throttle body and surge-tank, which can be categorized into the quadrupole source. The latter induces the unsteady force on the flow, which can be classified into the dipole source. In this paper, the mechanism of noise generation from the turbulence is only investigated as a preliminary study. Stochastic noise source synthesis method is adopted for the analysis of turbulence-induced, i.e. quadrupole noise by throttle at quick opening state. The method consists of three procedures. The first step corresponds to the preliminary time-averaged Navier-Stokes computation with a $k-\varepsilon$ turbulence model providing mean flow field characteristics. The second step is the synthesis of time-dependent turbulent velocity field associated with quadrupole noise sources. The final step is devoted to the determination of acoustic source terms associated with turbulent velocity. For the first step, we used market available analysis tools such as STAR-CD, the trade names of fluid analysis tools available on the market. The steady state flows at three open angle of throttle valve, i.e. 20, 35 and 60 degree, are numerically analyzed. Then, time-dependent turbulent velocity fields are produced by using the stochastic model and the flow analysis results. Using this turbulent velocity field, the turbulence-originated noise sources, i.e. the self-noise and shear-noise sources are synthesized. Based on these numerical results, it is found that the origin of the turbulent flow and noise might be attributed to the process of formulation and the interaction of two vortex lines formed in the downstream of the throttle valve. These vortex lines are produced by the non-uniform splits between the throttle valve and inner cylinder surface. Based on the analysis, we present the low-noise design of the inner geometry of throttle body.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.27 no.6
• /
• pp.699-710
• /
• 2009

If a country has her geoid model, it could be determine accurate orthometric height because the geoid model could provide continuous equi-gravity potential surface. And it is possible to improve the coordinates accuracy of national control points through geodetic network adjustment considering geoidal heights. This study aims to find the best way to develop geoid model in Philippines which have similar topographic conditions as like Malaysia and Indonesia in Eastsouth asia. So, in this study, it is surveyed the general theories of geoid determination and development cases of geoid model in Asia and it is computed that the geoidal heights and gravity anomalies by spherical harmonic analysis using EGM2008, the latest earth geopotential model. The results show that first, the development of gravimetric geoid model based on airborne gravimetry is needed and second, about 200 GPS surveying data at national benchmark is needed. It is concluded that it is the most reasonable way to develop the hybrid geoid model through fitting geometric geoid by GPS/leveling data to gravimetric geoid. Also, it is proposed that four band spherical Fast fourier transformation(FFT) method for evaluation of Stokes integration and remove and restore technique using EGM2008 and SRTM for calculation of gravimetric geoid model and least square collocation algorithm for calculation of hybrid geoid model.

• Korean Journal of Food Science and Technology
• /
• v.47 no.2
• /
• pp.158-163
• /
• 2015

Fustin comprised >98 and 73.0-86.7% of the total flavonoid content in the bark, and the stems and lignum of Rhus verniciflua, respectively. The butein, fisetin, and sulfuretin content varied between 0.31-2.17, 0.27-3.32, and 0.15-0.80 mg/g on a dry weight basis, respectively, in different parts of Rhus verniciflua. The urushiol content was 5.09-6.29, 55.05-56.30, and 0.38-0.39 mg/100 g on a dry weight basis in the stems, bark, and lignum, respectively. This showed that the bark of the tree had the highest urushiol content. C15:3 (pentadecatrienyl catechol), C15:1 (pentadecenyl catechol), and C15:2 (pentadecadienyl catechol) comprised 63, 33-35, and 2-3% of urushiol congeners in the tree bark, respectively.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.214-214
• /
• 2003

During a model test of Hutton TLP, a "ringing" response was first observed about 20 years ago. This phenomenon is a resonant build up over the time of wave period and this burst-like motion can cause the extreme load on the TLP tether. It is often detected in the large and steep irregular waves but the generation mechanism leading to the "ringing" is not yet well understood. According to the research since then, the higher order harmonic components may account for the "ringing" on the floating offshore structures. The main purpose of the present research is, thus, to measure the higher harmonic forces exerted on a vertical truncated circular column and to compare them with available data. A vertical truncated cylinder with a diameter of 3.5inch and a draft of 10.5inch is used as a test structure, which is a scaled model of ISSC TLP column. The cylinder is installed at a distance of 45ft from the wave maker in order to avoid parasitic waves created in the wave flap. Attached to the upper part of the cylinder are two force gages to measure the horizontal (surge) and vertical (heave) forces on the cylinder. The incoming waves are Stokes waves with a slope ranging from 0.06 to 0.24. The forces and waves are measured for 60 seconds with a sampling rate of 50 Hz. Among the recorded data, the first 10 waves are excluded because of transient behavior of the waves and the next The horizontal and vertical forces are analyzed up to 5th order harmonics. The horizontal forces are then compared to the values from the theoretical model called "FNV model". In addition, force transfer functions are also investigated. Major findings in this research are below. 1) The first order forces measured are slightly larger than the theoretical values of "FNV model" 2) The "FNV model" considerably overpredicts the second order forces. 3) The larger the amplitude and more extreme the wave slope, the smaller the predictions are compared to the experimental. 4) The higher harmonic forces are significantly smaller than the first harmonic force for all wave parameters. 5) The normalized forces vs. waves slopes are almost constant in the lower harmonics but vary a lot in the higher harmonics. 6) The trend of forces is more nonlinear in the horizontal forces than in the vertical forces as the wave slope increases. 7) The part of the results above is also observed by other researchers and confirmed again through the present work.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.169-169
• /
• 2015

유사 입자의 크기는 유사의 특성 및 그에 따른 거동 변화에 중요한 영향을 미친다. Stokes 침강 속도 모형에서 유사의 침강 속도에 가장 많은 영향을 주는 인자는 유사의 크기인 것 또한 확인된다. 유사 입자의 크기가 약 $60{\mu}m$보다 작은 유사들은 알갱이 사이의 점착력을 무시할 수 없다. 이로 인해 유사들은 응집 현상을 겪으며 입자 본래의 크기보다 크기가 큰 플럭을 형성하는 점착성 유사로 분류된다. 응집 현상이란, 흐름 내 점착성을 띠는 일차입자(Primary Particle)가 응집과 파괴를 반복하며 플럭을 형성하는 현상을 뜻한다. 입자 간의 충돌을 통해 응집이 진행되며 난류 전단으로 인해 형성된 플럭의 파괴가 발생한다. 많은 연구에서 점착성 유사의 충돌을 야기하는 가장 지배적인 원리는 난류라 알려져 있다. 이러한 응집 현상으로 인하여 플럭의 크기와 밀도는 지속적으로 변화를 겪으며 비점착성 유사와 다른 특징들을 보인다. 흐름에 존재하는 유사의 이동은 이송-확산 방정식을 통해 표현된다. 이송-확산 방정식은 시간 변화에 따른 농도의 변화를 입자의 침강과 난류 및 유사 자체의 특징에 의한 확산으로 해석한다. 침강속도로 대변되는 이송과 달리, 확산은 난류흐름 내에서 유사가 확산되는 정도를 정량화하기 위한 인자가 요구된다. 난류에 의한 유사의 확산은 유사 자체 특성에 따른 물질 확산에 비하여 매우 큰 값을 가지며, 이를 확산 계수로 개념화 한다. 확산계수는 와점성계수와 Schmidt 수(${\sigma}_c$)의 비로 정의된다. ${\sigma}_c$는 난류의 점성과 난류로 인한 부유과정에 의해 유사가 확산되는 정도를 나타낸다. 이에 따라 ${\sigma}_c$의 변화가 유사의 부유 및 침강거동에 많은 영향을 미칠 것이라 판단되나, 국내외에서 수행된 연구 동향에서는 ${\sigma}_c$를 0.5부터 1.0 사이의 상수를 적용하여 수행되었다. 이에 본 연구에서는 ${\sigma}_c$의 크기에 따라 달라지는 유사의 부유 및 침강 변화에 의한 총 부유량을 살펴보고자 한다. 유사의 점착성을 고려할 수 있는 1DV 수치 모형을 이용하여 비점착성 유사와 점착성 유사를 대상으로 수치연구를 수행하며, 유사의 크기 및 ${\sigma}_c$의 변화에 따른 총 부유량 경향을 살펴본다. 그 결과, 점착성 유사는 ${\sigma}_c$의 증가에 따라서 유사의 총 부유량이 증가하는 현상이 나타난 반면 비점착성 유사는 ${\sigma}_c$의 증가에 따라 유사의 총 부유량이 감소하는 경향이 나타났다. 그러나 크기가 아주 작은 비점착성 유사를 대상으로 수치 연구를 수행한 결과, ??에 따른 총 부유량의 경향은 유사의 점착성에서 기인하는 것이 아닌 입자의 크기로부터 야기되는 특성이라는 결론이 도출되었다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.25-34
• /
• 2008

In this paper, three-dimensional multiblade row unsteady Navier-Stokes simulations at a hot streak temperature ratio of 2.0 have been performed to reveal the effects of rotor tip clearance on the inlet hot streak migration characteristics in low pressure stage of a Vaneless Counter-Rotating Turbine. The hot streak is circular in shape with a diameter equal to 25% of the high pressure turbine stator span. The hot streak center is located at 50% of the span and the leading edge of the high pressure turbine stator. The tip clearance size studied in this paper is 2.0mm(2.59% high pressure turbine rotor height, and 2.09% low pressure turbine rotor height). The numerical results show that the hot streak is not mixed out by the time it reaches the exit of high pressure turbine rotor. The separation of colder and hotter fluid is observed at the inlet of low pressure turbine rotor. Most of hotter fluid migrates towards the rotor pressure surface, and only little hotter fluid migrates to the rotor suction surface when it convects into the low pressure turbine rotor. And the hotter fluid migrated to the tip region of the high pressure turbine rotor impinges on the leading edge of the low pressure turbine rotor after it goes through the high pressure turbine rotor. The migration of the hotter fluid directly results in very high heat load at the leading edge of the low pressure turbine rotor. The migration characteristics of the hot streak in the low pressure turbine rotor are dominated by the combined effects of secondary flow and leakage flow at the tip clearance. The leakage flow trends to drive the hotter fluid towards the blade tip on the pressure surface and to the hub on the suction surface, even partial hotter fluid near the pressure surface is also driven to the rotor suction surface through the tip clearance. Compared with the case without rotor tip clearance, the heat load of the low pressure turbine rotor is intensified due to the effects of the leakage flow. And the numerical results also indicate that the leakage flow effect trends to increase the low pressure turbine rotor outlet temperature at the tip region.

• Solar Energy
• /
• v.14 no.2
• /
• pp.75-90
• /
• 1994

Thermal energy transport in a two-dimensional horizontal and vertical channel with an isothermal rectangular beam attached to one adiabatic wall is investigated from the numerical solution of Navier-Stokes and energy equations. The solutions have been obtained for dimensionless aspect equations. The solutions have been obtained for dimensionless aspect ratios of beam, H/B=$0.25{sim}4$, Reynolds numbers, Re=$50{\sim}500$ and Grashof numbers, Gr=$0{\sim}5{\times}10^4$. The mean Nusselt number, $\overline{Nu}$ for horizontal and vertical channels shows same value at Gr=0 and increases as Gr increases and decreases as H/B increases at Re=100. $\overline{Nu}$ of vertical channel shows higher in $0.25{\leq}H/B<1.1$ and lower in $1.1{\leq}H/B{\leq}4.0$ than that of horizontal channel at $Gr=10^4$, Re=100. $\overline{Nu}$ of vertical channel shows higher in $0.25{\leq}H/B<1.1$ and lower in $1.1{\leq}H/B=1.0$ than that of horizontal channel at Re=100, $0<Gr{\leq}5{\times}10^4$. A comparison between the experimental and numerical results shows good agreement.