• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.2971-2980
• /
• 1995

A numerical study of fluid flow in driven cavity was carried out using singular finite element method. The driven cavity problem is known to have infinite velocity gradients as well as dual velocity conditions at the singular points. To overcome such difficulties, a finite element method with singular shape functions was used and a special technique was employed to allow multiple values of velocities at the singular points. Application of singular elements in the driven cavity problem has a significant influence on the stability of solution. It was found the singular elements gave a stable solution, especially, for the pressure distribution of the entire flow field by keeping up a large pressure at the singular points. In the existing solutions of driven cavity problem, most efforts were focused on the study of streamlines and vorticities, and pressure were seldom mentioned. In this study, however, more attention was given to the pressure distribution. Computations showed that pressure decreased very rapidly as the distance from the singular point increased. Also, the pressure distribution along the vertical walls showed a smoother transition with singular elements compared to those of conventional method. At the singular point toward the flow direction showed more pressure increase compared with the other side as Reynolds number increased.

• Journal of the Korean Society of Visualization
• /
• v.19 no.3
• /
• pp.136-142
• /
• 2021

Submerged and semi-submerged vehicles expel cooling water through an outlet. In this process, induced noise and vibration by the flow around the outlet have been reported, and it may cause problems directly related to survivability of the navy vessels. The coolant outlet has a net-type structure and circular columns are mostly used. In this study, flow measurements using PIV and LDV were performed for different type outlets; conventional (flat plate with round bar) and improved (flat and flat plate) configurations. Experiments were conducted at a cavitation tunnel where pressure and steady flow rate conditions are ensured for sufficient time to measure the flow. The average velocity field of the outlets were measured and compared through LDV measurements, and instantaneous vorticities were evaluated through PIV measurements. The results show that the improved type of the outlet is advantageous in terms of flow stability compared to the conventional type of the outlet.

• Advances in aircraft and spacecraft science
• /
• v.8 no.4
• /
• pp.303-330
• /
• 2021

The paper discusses the effect of the winglets on the aerodynamic and aeroacoustic performance of Boeing 737-800 aircraft by numerical approach. For this purpose, computational fluid dynamics and fluent commercial software are used to solve the compressible flow governing equations. The RANS method and the K-ω SST turbulence model are selected to simulate the subsonic flow around the wing with acceptable accuracy and low computational cost. The main variables of steady flow around the simple and blended wing in constant atmospheric conditions are computed by numerical solution of governing equations. The solution of the acoustic field has also been accomplished by the broad-band acoustic source model. The results reveal that adding a blended winglet increases the pressure difference near the wingtip,which increases the lift force. Also, the blended winglet reduces the power and magnitude of vorticities around the wingtip, which reduces the wing's drag force. The effects of winglets on aerodynamic forces lead to a 3.8% increase in flight range and a 3.6% increase in the maximum payload of the aircraft. Also, the acoustic power level variables on the surfaces and fields around the wing have been investigated integrally and locally.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.8 no.3
• /
• pp.10-16
• /
• 2004

Untypical unstable burning with very low frequency was observed at firing test of a gas generator using bundle cylindrical grain. The pressure unbalance between inside and outside of cylindrical grain brought such a low unstable burning. The grains were radially holed so that the high pressure gas inside of grain could quickly moved outward of gain, resulting dissipation of the pressure unbalance However too many holes were required to let the burning be stable for all operation regime from low to high temperature of grain and resultantly deteriorate the Progressive increase of gas amount produced by a gas generator. So another idea using grids located both sides of a bundle grain was applied to dissipate actively large vorticities enhanced by unbalance pressure distribution in a combustor. Finally stable burning with progressively increase of gas was established by application of 5${\times}$5 grid slightly away bundle grain to move bundle gain freely in case pressure unbalance were occurred inside of combustor.

• Wind and Structures
• /
• v.26 no.4
• /
• pp.191-204
• /
• 2018

This article presents a study of the largest-ever (height = 220 m) cooling tower using the large eddy simulation (LES) method. Information about fluid fields around the tower and 3D aerodynamic time history in full construction process were obtained, and the wind pressure distribution along the entire tower predicted by the developed model was compared with standard curves and measured curves to validate the effectiveness of the simulating method. Based on that, average wind pressure distribution and characteristics of fluid fields in the construction process of ultra-large cooling tower were investigated. The characteristics of fluid fields in full construction process and their working principles were investigated based on wind speeds and vorticities under different construction conditions. Then, time domain characteristics of ultra-large cooling towers in full construction process, including fluctuating wind loads, extreme wind loads, lift and drag coefficients, and relationship of measuring points, were studied and fitting formula of extreme wind load as a function of height was developed based on the nonlinear least square method. Additionally, the frequency domain characteristics of wind loads on the constructing tower, including wind pressure power spectrum at typical measuring points, lift and drag power spectrum, circumferential correlations between typical measuring points, and vertical correlations of lift coefficient and drag coefficient, were analyzed. The results revealed that the random characteristics of fluctuating wind loads, as well as corresponding extreme wind pressure and power spectra curves, varied significantly and in real time with the height of the constructing tower. This study provides references for design of wind loads during construction period of ultra-large cooling towers.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.5
• /
• pp.819-830
• /
• 2016

Hydrodynamic characteristics of wave propagation through density stratification have not been identified in details. So this study conducted a numerical simulation using LES-WASS-3D ver. 2.0 for analysis of density current due to water temperature and salinity in order to analyze hydraulic characteristics under wave action in a two-layer density stratified fluid. For the validity and effectiveness of numerical wave tank used, it was compared and analyzed with the experiment to show waveform based on $3^{rd}$-order Stoke wave theory at the internal of a density stratification. Using the results obtained from numerical simulation, the surface and internal wave heights are reduced as the wave propagates in a two-layer density stratified water. And the surface or internal wave attenuation became more serious as the vorticities were increased by the velocity difference of wave propagation due to the upper-lower density difference around the interface of a density stratification. As well, the surface and internal wave attenuations became more serious with higher density difference and depth ratio between upper and lower layers when the wave propagates through a density stratification.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.449-452
• /
• 2008

Flow difference between reacting and non-reacting case in a swirl stabilized annular combustor is investigated using 3D Large Eddy Simulation with flamelet turbulent combustion model. The combustor of concern is the LM6000, lean premixed dry low-NOx annular combustor, developed by GEAE. Boundary conditions are based on experimental data. Heat release as a result of combustion put the dilatation of density in primary combustion zone highly increased so that the main swirl stream behind of a swirl cup stretched further downstream than that of non-reacting case. The oval shape of core flow in cross-section to flow direction, which clearly observed in non-reacting case, tends to be circle, and small vorticities in wide range in non-reacting case disappears, but the size of iso-vorticity increase in reacting case.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.4
• /
• pp.405-414
• /
• 2008

To better understand the physical processes that control the high-latitude lower thermospheric dynamics, we analyze the divergence and vorticity of the high-latitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions. For this study the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEG CM) is used. The analysis of the large-scale vorticity and divergence provides basic understanding flow configurations to help elucidate the momentum sources that ulti-mately determine the total wind field in the lower polar thermosphere and provides insight into the relative strengths of the different sources of momentum responsible for driving winds. The mean neutral wind pattern in the high-latitude lower thermosphere is dominated by rotational flow, imparted primarily through the ion drag force, rather than by divergent flow, imparted primarily through Joule and solar heating. The difference vorticity, obtained by subtracting values with zero IMF from those with non-zero IMF, in the high-latitude lower thermosphere is much larger than the difference divergence for all IMF conditions, indicating that a larger response of the thermospheric wind system to enhancement in the momentum input generating the rotational motion with elevated IMF than the corresponding energy input generating the divergent motion. the difference vorticity in the high-latitude lower thermosphere depends on the direction of the IMF. The difference vorticity for negative and positive $B_y$ shows positive and negative, respectively, at higher magnetic latitudes than $-70^{\circ}$. For negative $B_z$, the difference vorticities have positive in the dusk sector and negative in the dawn sector. The difference vorticities for positive $B_z$ have opposite sign. Negative IMF $B_z$ has a stronger effect on the vorticity than does positive $B_z$.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.408-417
• /
• 2007

Recently, an artificial beach has been constructed compensating for loss of the natural one caused by the development of coastal area, as well as serving as a location for recreational activities such as sea bathing. It is well known that some structure should be constructed to protect an artificial beach from the outflow due to wave action of the reclaimed sand. In general, dike is utilized as the structure to protect an artificial beach. And, one of the factors which may need to be taken into consideration for stability of dike on seabed foundation is the ground water behavior behind dike. However, the interrelated phenomena of nonlinear wave and ground water response have relatively little attention although these interactions are important for stability of structure and sand suction to the artificial beach. In this paper, the numerical wave tank was developed to clarify nonlinear wave, dike and ground water dynamic interaction, which can simulate the difference of ground water and mean water level. Using the developed numerical wave tank, the present study investigates how variation of ground water level influences hydrodynamic characteristics in seabed around dike and numerically simulates the wave fields, pore flow patterns, pore water pressures and vorticities according to variation of ground water level. Numerical results explain well how hydrodynamic characteristics in seabed around dike is affected by the variation of ground water level.

• Journal of Environmental Impact Assessment
• /
• v.24 no.5
• /
• pp.456-471
• /
• 2015

Emission source term is one of the strong controlling factors for the air quality simulation capability, particularly over the urban area. Ulsan is an industrial area and frequently required to simulate for environmental assessment. In this study, two CAPSS (Clean Air Policy Support System) emission data; CAPSS-2003 and CAPSS-2010 in Ulsan, were employed as an input data for WRF-CMAQ air quality model for emission assessment. The simulated results were compared with observations for the local emission dominant synoptic conditions which had negative vorticities and lower geostrophic wind speed at 850hPa weather maps. The measurements of CO, $NO_2$, $SO_2$ and $PM_{10}$ concentrations were compared with simulations and the 'scaling factors' of emissions for CO, $NO_2$, $SO_2$, and $PM_{10}$ were suggested in in aggregative and quantitative manner. The results showed that CAPSS-2003 showed no critical discrepancies of CO and $NO_2$ observations with simulations, while $SO_2$ was overestimated by a factor of more than 12, while $PM_{10}$ was underestimated by a factor of more than 20 times. However, CAPSS-2010 case showed that $SO_2$ and $PM_{10}$ emission were much more improved than CAPSS-2003. However, $SO_2$ was still overestimated by a factor of more than 2, and $PM_{10}$ underestimated by a factor of 5, while there was no significant improvement for CO and $NO_2$ emission. The estimated factors identified in this study can be used as'scaling factors'for optimizing the emissions of air pollutants, particularly $SO_2$ and $PM_{10}$ for the realistic air quality simulation in Ulsan.