• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.

• 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.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4357-4366
• /
• 2023

MARS-KS, a domestic regulatory confirmatory code of Republic of Korea, had been developed by integrating RELAP5/MOD2 and COBRA-TF. The integration of COBRA-TF allowed to extend the capability of MARS-KS, limited to one-dimensional analysis, to multi-dimensional analysis. The use of COBRA-TF was mainly focused on subchannel analyses for simulating multi-dimensional behavior within the reactor core. However, this feature has been remained as a legacy without ongoing maintenance. Meanwhile, MARS-KS also includes its own multidimensional component, namely MULTID, which is also feasible to simulate three-dimensional convection and diffusion. The MULTID is capable of modeling the turbulent diffusion using simple mixing length model. The implementation of the turbulent mixing is of importance for analyzing the reactor core where a disturbing cross-sectional structure of rod bundle makes the flow perturbation and corresponding mixing stronger. In addition, the presence of this turbulent behavior allows the secondary transports with net mass exchange between subchannels. However, a series of assessments performed in previous studies revealed that the turbulence model of the MULTID could not simulate the aforementioned effective mixing occurred in the subchannel-scale problems. This is obvious consequence since the physical models of the MULTID neglect the effect of mass transport and thereby, it cannot model the void drift effect and resulting phasic distribution within a bundle. Thus, in this study, the turbulence mixing model of the MULTID has been improved by means of the inter-channel mixing model, widely utilized in subchannel analysis, in order to extend the application of the MULTID to small-scale problems. A series of assessments has been performed against rod bundle experiments, namely GE 3X3 and PSBT, to evaluate the performance of the introduced mixing model. The assessment results revealed that the application of the inter-channel mixing model allowed to enhance the prediction of the MULTID in subchannel scale problems. In addition, it was indicated that the code could not predict appropriate phasic distribution in the rod bundle without the model. Considering that the proper prediction of the phasic distribution is important when considering pin-based and/or assembly-based expressions of the reactor core, the results of this study clearly indicate that the inter-channel mixing model is required for analyzing the rod bundle, appropriately.

• Journal of Korea Water Resources Association
• /
• v.57 no.5
• /
• pp.347-357
• /
• 2024

Numerous prior studies have delineated the size distribution of noncohesive sediment in suspension, focusing on mean size and standard deviation. However, suspensions comprise a heterogeneous mixture of sediment particles of varying sizes. The transport dynamics of suspended sediment in turbulent flow are intimately tied to settling velocities calculated based on size and density. Consequently, understanding the grain size distribution becomes paramount in comprehending sediment transport phenomena for noncohesive sediment. This study aims to introduce a straightforward modeling approach for simulating the grain size distribution of suspended sediment amidst turbulence. Leveraging insights into the contrast between cohesive and noncohesive sediment, we have meticulously revised a stochastic flocculation model originally designed for cohesive sediment to aptly simulate the grain size distribution of noncohesive sediment in suspension. The efficacy of our approach is corroborated through a meticulous comparison between experimental data and the grain size distribution simulated by our newly proposed model. Through numerical simulations, we unveil that the modulation of grain size distribution of suspended sediment is contingent upon the sediment transport capacity of the carrier fluid. Hence, we deduce that our simplified approach to simulating the grain size distribution of suspended sediment, integrated with a sediment transport model, serves as a robust framework for elucidating the pivotal bulk properties of sediment transport.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.849-861
• /
• 1997

In this study, the effect of hole geometry of the cooling system on the flow and temperature field was numerically calculated. The finite volume method was employed to discretize the governing equation based on the non-orthogonal coordinate with non-staggered variable arrangement. The standard k-.epsilon. turbulence model was used and also the predicted results were compared with the experimental data to validate numerical modeling. The predicted results showed good agreement in all cases. To analyze the effect of the discharge coefficient for slots of different length to width, the inlet chamfering and radiusing holes were considered. The discharge coefficient was increased with increment of the chamfering ratio, radiusing ratio and slot length to width and also the effect of radiusing showed better result than chamfering in all cases. In order to analyze the difference between the predicted results with plenum region and without plenum region, the velocity profiles of jet exit region for a various flow conditions were calculated. The normal velocity components of jet exit showed big difference for the low slot length to width and high blowing rate cases. To analyze the flow phenomena injected from a row of inclined holes in a real turbine blade, three dimensional flow and temperature distribution of the region including plenum, hole and cross stream with flow conditions were numerically calculated. The results have shown three-dimensional flow characteristics, such as the development of counter rotating vortices, jetting effect and low momentum region within the hole in addition to counter rotating vortex structure in the cross stream.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.40 no.4
• /
• pp.407-414
• /
• 2016

As wind turbines are getting larger in size with multi-MW capacity, the blades are getting longer, over 40 m, and hence the asymmetric loads produced during the rotation of the rotor blades are increasing. Some factors such as wind shear, tower shadow, and turbulence have an effect on the asymmetric loads on the blades. This paper focuses on a method of modeling the dynamic load acting on a blade because of thrust variation under wind shear. A method that uses thrust coefficient is presented. For this purpose, "wind shear coefficient of thrust variation" is defined and introduced. Further, we calculate the values of the "wind shear coefficient of thrust variation" for a 2 MW on-shore wind turbine, and analyze them for speeds below the rated wind speed. Then, we implement a dynamic model that represents the thrust variation under wind shear on a blade, using MATLAB/Simulink. It is shown that it is possible to express thrust variations on three blades under wind shear by using both thrust coefficient and "wind shear coefficient of thrust variation."

• Journal of Korean Society of Occupational and Environmental Hygiene
• /
• v.18 no.2
• /
• pp.161-170
• /
• 2008

푸쉬-풀 환기시스템은 도금조와 같이 흡인해야 할 거리가 상대적으로 긴 경우에 많이 사용되고 있다. 그러나, 창문이나 출입문을 통한 방해기류가 푸쉬-풀 환기시스템의 오염물질 제어효율을 심각하게 훼손시키고 있다고 추측하고 있으나 이에 대한 세부적인 연구가 부족한 상태에 있다. 따라서, 본 연구에서는 전산유체역학(Computational fluid dynamics)을 이용하여 푸쉬-풀 환기시스템에서의 방해기류의 방향과 세기가 흡인효율에 어떠한 영향을 미치는지에 대해 평가해 보았다. 선형흡인효율(Linear capture efficiency) 방법을 이용하여 푸쉬-풀 환기시스템에서 가상의 개방조에서 발생한 오염물질이 푸쉬-풀 시스템에 의하여 포집되지 못하고 누출되는 구역이 어딘지를 찾아낼 수 있었다. 전산유체역학 컴퓨터시뮬레이션은 AIRPAK2.1 (FLUENT CODE) 소프트웨어를 사용하였다. 푸쉬-풀 후드시스템에 방해기류가 강하게 작용하면 상대적으로 강한 와류가 발생하는데, 일반적인 난류모델인 ${\kappa}-{\varepsilon}$모델은 와류현상을 충분히 보여주지 못한 반면에 RNG 모델을 사용했을 때 실험결과를 적절히 모사해낼 수 있었다. RNG 모델을 이용하여 세가지 방향, 즉 푸쉬에서 풀 방향으로, 풀에서 푸쉬 방향으로 그리고 그에 수직되는 방향으로 방해기류가 있을 때의 푸쉬-풀 환기시스템의 흡인효율을 분석하였다. 방해기류가 0.25m/s이하일 때에는 흡인효율이 거의 떨어지지 않았으나, 방해기류가 0.6m/s에서 흡인효율이 40-70%로 떨어짐을 알 수 있었다. 따라서, 방해기류를 감소시킬 수 있는 방안에 대해서도 연구를 해야 되겠지만, 방해기류 존재 하에서 충분한 흡인 효율을 유지할 수 있는 푸쉬-풀 후드 설계기준에 대한 연구도 필요할 것으로 판단된다.

• Journal of Korean Society for Atmospheric Environment
• /
• v.22 no.1
• /
• pp.73-84
• /
• 2006

Hazardous Air Pollutants (HAPs) are characterized by being relatively heavier and denser than that of ambient air due to the various reasons such as higher molecular weight, low temperature and other complicated chemical transformations (Witlox, 1994). In an effort to investigate transport and diffusion from instantaneous emission of heavy gas, Lagrangian Particle Dispersion Model (LPDM) coupled with the RAMS output was employed. Both deposition process and buoyancy term were added on the atmospheric diffusion equations of LPDM, and the locations and concentrations of dense gas particle released from instantaneous single point source (emitting initially for 10 minutes only) were analyzed. The result overall shows that adding deposition process and buoyancy terms on the diffusion equation of LPDM has very small but detectable effect on the vertical and horizontal distribution of Lagrangian particles that especially transported for a fairly long traveling time. Also the slumping of dense gas can be found to be ignored horizontally compared to the advection by the horizontal wind suggesting that it was essential to couple the Lagrangian particle dispersion model coupled with the RAMS model in order to explain the dispersion of HAPs more accurately. However, during the initial time of instantaneous emission, buoyancy term play an important role on the vertical locations of dense particles for near surface atmosphere and around source area, indicating the importance of densities of HAPs in the beginning stage or short duration for the risk assessment of HAPs or management of heavy vapors during the explosive accidents.

• International Journal of Fluid Machinery and Systems
• /
• v.5 no.3
• /
• pp.134-142
• /
• 2012

The competition to deliver ultra-low emitting vehicles at a reasonable cost is driving the automotive industry to invest significant manpower and test laboratory resources in the design optimization of increasingly complex exhaust after-treatment systems. Optimization can no longer be based on traditional approaches, which are intensive in hardware use and laboratory testing. The CFD is in high demand for the analysis and design in order to reduce developing cost and time consuming in experiments. This paper describes the development of a comprehensive practical model based on experiments for simulating the performance of automotive three-way catalytic converters, which are employed to reduce engine exhaust emissions. An experiment is conducted to measure species concentrations before and after catalytic converter for different loads on engine. The model simulates the emission system behavior by using an exhaust system heat conservation and catalyst chemical kinetic sub-model. CFD simulation is used to study the performance of automotive catalytic converter. The substrate is modeled as a porous media in FLUENT and the standard k-e model is used for turbulence. The flow pattern is changed from axial to radial by changing the substrate model inside the catalytic converter and the flow distribution and the conversion efficiency of CO, HC and NOx are achieved first, and the predictions are in good agreement with the experimental measurements. It is found that the conversion from axial to radial flow makes the catalytic converter more efficient. These studies help to understand better the performance of the catalytic converter in order to optimize the converter design.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.9
• /
• pp.2944-2952
• /
• 1996

In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 transoni $c^ersonic compressor cascade flow. First, the general characteristics of baseline cascade flow were analyzed. At freestream Mach n.1.612 and exit/inlet pressure ratio 2.15, the results from current laminar flow were compared well in suction surface with the experiment; however, not well in pressure surface. Second, numerical study of the transoni $c^ersonic compressor cascade flow demonstrated the effectiveness of a passive control by the various size cavities. A cavity under the shock foot point at the suction surface of the blades was used as a passive control. The passive control of shock-boundary layer interaction by a cavity reduced total pressure losses. The effect of cavity length and depth was studied. The total pressure loss was reduced by about 10% and the isentropic efficiency was improved slightly. The effect of cavity depth in current study(d/l = 0.05, 0.02) was not found strong. Further adequate turbulence modeling and TVD schemes would help to capture the shock more accurately and increase the effectiveness of the current shock-boundary layer interaction study using upwind flux difference splitting computational methods.thods.