• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.2
• /
• pp.83-92
• /
• 1998

To evaluate the hydraulic resistance behind bodies in a large scale grid numerical model, a drag stress term which is formulated by the drag force is introduced in the depth-integrated Reynolds equations. And also, the applicability and problems of this model are discussed through various numerical experiments where the analytical solutions exist. In the case of a single body, the error range of velocity difference between analytical and numerical solutions is within $\pm$10% and the wake width behind the body shows a good agreement with the analytical solution. When the drag coefficient and the eddy viscosity are precisely decided, the numerical solutions behind a row of bodies will be efficiently used in real situations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.5
• /
• pp.1-9
• /
• 2017

The cryogenic nitrogen spray of a swirl injector has been numerically investigated using three dimensional LES turbulence model to analyze the dynamic characteristics under supercritical condition. To predict the precise nitrogen properties under supercritical condition, SRK equation of state, Chung's method for viscosity and thermal conductivity and Takahashi's correlation based on Fuller's theory for diffusion coefficient are implemented. The complex flow structures due to interaction between flow field and acoustic field are observed inside and outside the injector under supercritical condition. FFT, POD, and DMD techniques are employed to understand the coherent structures. By implementing the FFT, the dominant frequencies are identified inside and outside the injector. The coherent flow structures related to the dominant frequencies are visualized using the POD and DMD techniques. In addition, the DMD provides the damping coefficient which is related with the instability prediction.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.115-125
• /
• 2022

The high-Reynolds number flow around a rectangular cylinder, having streamwise to crossflow length ratio equal to 5 is analyzed in the present paper. The flow is characterized by shear-layer separation from the upstream edges. Vortical structures of different size form from the roll-up of these shear layers, move downstream and interact with the classical vortex shedding further downstream in the wake. The corresponding mean flow is characterized by a recirculation region along the lateral surface of the cylinder, ending by mean flow reattachment close to the trailing edge. The mean flow features on the cylinder side have been shown to be highly sensitive to set-up parameters both in numerical simulations and in experiments. The results of 21 Large Eddy Simulations (LES) are analyzed herein to highlight the impact of the lateral mean recirculation characteristics on the near-wake flow features and on some bulk quantities. The considered simulations have been carried out at Reynolds number Re=DU_∞/ν=40 000, being D the crossflow dimension, U_∞ the freestream velocity and ν the kinematic viscosity of air; the flow is set to have zero angle of attack. Some simulations are carried out with sharp edges (Mariotti et al. 2017), others with different values of the rounding of the upstream edges (Rocchio et al. 2020) and an additional LES is carried out to match the value of the roundness of the upstream edges in the experiments in Pasqualetto et al. (2022). The dimensions of the mean recirculation zone vary considerably in these simulations, allowing us to single out meaningful trends. The streamwise length of the lateral mean recirculation and the streamwise distance from the upstream edge of its center are the parameters controlling the considered quantities. The wake width increases linearly with these parameters, while the vortex-shedding non-dimensional frequency shows a linear decrease. The drag coefficient also linearly decreases with increasing the recirculation length and this is due to a reduction of the suctions on the base. However, the overall variation of C_D is small. Finally, a significant, and once again linear, increase of the fluctuations of the lift coefficient is found for increasing the mean recirculation streamwise length.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.3
• /
• pp.132-139
• /
• 1997

CTD and current observation were taken to investigate the structure of the cold water in the Western Channel of the Korea Strait in October 1993. Thickness of the cold water in the deep trough of the strait changes from 20 m to 70 m according to the water depth. Thermocline between the Tsushima Warm Water and the cold water deepens from north to south with 0.00057 in slope. Temporal variation of the thickness appears to be related with the tidal current. The maximum variation is 20 m for 48 hours. Mean velocity of the cold water for 72 hours is 17 cm/sec southward. A simple model was used to understand dynamically the southward flow of the cold water and the return flow at the upper part in the lower layer. Calculated maximum southward flow and eddy viscosity coefficient are 7 cm/sec and 0.038 $m^2$/sec respectively in the model. Southward transport is $0.032$\times$10^6㎥/sec$ at the northern part in the trough and decreases from north to south due to the presence of the return flow. Southward transport increases with the increase in the upper layer transport but is not affected by the density of the upper layer or the interface slope.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.3
• /
• pp.357-366
• /
• 1986

The flow characteristics of two-dimensional turbulent offset jet which is discharged parallel to a solid wall has been studied experimentally and numerically. In the experiment, 3-hole pitot tube and 2 channel constant temperature hot-wire anemometer are used to measure local mean velocity, turbulence intensity and Reynolds stress while scannivalve is used to measure the wall pressure distribution. It is confirmed experimentally that local mean velocity is closely related to wall pressure distribution. It is also verified that for large Reynolds numbers and fixed step height there exists a similarity in the distribution of wall pressure coefficient. The maximum values of turbulence intensity occur in the top and bottom mixing layers and the magnitude of Reynolds stress becomes large in the lower mixing layer than in the top mixing layer due to the effect of streamline curvature and entrainment. In the numerical analysis, standard k-.epsilon. model based on eddy viscosity model and Leschziner and Rodi model based on algebraic stress model are adopted. The numerical analyses predict shorter reattachment lengths than the experiment, and this difference is judged to be due mainly to the problem of turbulence model constants and numerical algorithm. This also causes the inconsistency between the two results for other turbulence quantities in the recirculation region and impingement region, which constitutes a subject of a continued future study.

• Journal of Environmental Science International
• /
• v.18 no.9
• /
• pp.971-981
• /
• 2009

An objective of this study is as follows: 1) performing sensitivity analysis and parameter estimation of RMA2 and RMA4 models for the West-Nakdong River, 2) drawing up alternatives of gates-operation for water-quality enhancement, and 3) quantitative evaluation of methodology of 'flow-restoration by gates-operation' among 'Comprehensive Plan Improving Water-Quality in the West-Nakdong River(WNR)' with the target water-quality(BOD at Nakbon-N point: below 4.3 mg/L). The parameters for the RMA2 (depth-averaged two-dimensional flow model) and RMA4 (depth-averaged two-dimensional water-quality model) were determined by sensitivity analysis. Result of parameter estimation for RMA2 and RMA4 models is $1,000\;Pa{\cdot}s$ of the eddy viscosity, 20 of the Peclet number, 0.025 of the Manning coefficient, and $1.0\;m^2/s$ of the diffusion coefficient. We have evaluated the effects of water-quality enhancement of the selected alternatives by numerical simulation technique with the models under the steady-state flow condition and the time-variant transport condition. Because of no-resuspension from river bottom and considering BOD as conservative matter, these simulation results slightly differ from real phenomena. In the case of $50\;m^3/s$ of Daejeo-gate inflow, two-dimensional flow pn results result represents that small velocity occurs in the Pyungkang Stream and no flow in the Maekdo River. In the WNR, there occurs the most rapid flow near timhae-bridge. In the WNR, changes of water-quality for the four selected simulation cases(6, 10, 30, $50\;m^3/s$ of the Daejeo-gate inflow) were predicted. Since the Daejeo-Gate and the Noksan-Gate can be opened up to 7 days, it would be found that sustainable inflow of $30\;m^3/s$ at the Daejeo-gate makes BOD in the WNR to be under the target of water-quality.

• Fisheries and Aquatic Sciences
• /
• v.13 no.2
• /
• pp.167-181
• /
• 2010

To predict changes in the marine environment of the Beolgyo Stream Estuary in Jeonnam Province, South Korea, where cohesive tidal flats cover a broad area and a large bridge is under construction, this study conducted numerical simulations involving tidal flow and cohesive sediment transport. A wetting and drying (WAD) technique for tidal flats from the Princeton Ocean Model (POM) was applied to a large-scale-grid hydrodynamic module capable of evaluating the flow resistance of structures. Derivation of the eddy viscosity coefficient for wakes created by structures was accomplished through the explicit use of shear velocity and Chezy's average velocity. Furthermore, various field observations, including of tide, tidal flow, suspended sediment concentrations, bottom sediments, and water depth, were performed to verify the model and obtain input data for it. In particular, geologic parameters related to the evaluation of settling velocity and critical shear stresses for erosion and deposition were observed, and numerical tests for the representation of suspended sediment concentrations were performed to determine proper values for the empirical coefficients in the sediment transport module. According to the simulation results, the velocity variation was particularly prominent around the piers in the tidal channel. Erosion occurred mainly along the tidal channels near the piers, where bridge structures reduced the flow cross section, creating strong flow. In contrast, in the rear area of the structure, where the flow was relatively weak due to the formation of eddies, deposition and moderated erosion were predicted. In estuaries and coastal waters, changes in the flow environment caused by artificial structures can produce changes in the sedimentary environment, which in turn can affect the local marine ecosystem. The numerical model proposed in this study will enable systematic prediction of changes to flow and sedimentary environments caused by the construction of artificial structures.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.