• Journal of Ocean Engineering and Technology
• /
• v.29 no.1
• /
• pp.9-15
• /
• 2015

We investigated the Reynolds number effects on the flow over a twisted offshore structure in the range of 3×10³≤ Re ≤ 1 × 10⁴. To analyze the effect of the twisted surface treatment, a large eddy simulation (LES) with a dynamic subgrid model was employed. A simulation of the cylindrical structure was also carried out to compare the results with those of the twisted offshore structure. As Re increased, the mean drag and lift coefficient of the twisted offshore structure increased with the same tendency as those of the cylindrical structure. However, the increases in the mean drag and lift coefficient of the twisted offshore structure were much smaller than those of the cylindrical structure. Furthermore, elongated shear layer and suppressed vortex shedding from the twisted offshore structure occurred compared to those of the cylindrical cylinder, resulting in a drag reduction and suppression of the vortex-induced vibration (VIV). In particular, the twisted offshore structure achieved a significant reduction of over 96% in VIV compared with that of the cylindrical structure, regardless of increasing Re. As a result, we concluded that the twisted offshore structure effectively controlled the flow structures with reductions in the drag and VIV compared with the cylindrical structure, irrespective of increasing Re.

• Journal of Power System Engineering
• /
• v.13 no.6
• /
• pp.13-21
• /
• 2009

This basic study is required to examine spray or jet behavior depending on fuel phase. In this study, analyses of diesel fuel(n-Tridecane, $C_{13}H_{28}$) spray and natural gas fuel(Methane, $CH_4$) jet under high temperature and pressure are performed by a general-purpose program, ANSYS CFX release 11.0, and the results of these are compared with experimental results of diesel fuel spray using the exciplex fluorescence method. The simulation results of diesel spray is analyzed by using the combination of Large-Eddy Simulation(LES) and Lagrangian Particle Tracking(LPT) and of a natural gas jet is analyzed by using Multi-Component Model(MCM). There are two study variables considered, that is, ambient pressure and injection pressure. In a macroscopic analysis, the higher ambient pressure is, the shorter spray or jet tip penetration is at each time after start of injection. And the higher injection pressure is, the longer spray or jet tip penetration is at each time after start of injection. When liquid fuel is injected, droplets of the fuel need some time to evaporate. However, when natural gas fuel is injected, the fuel does not need time to evaporate. Gas fuel consists of minute particles. Therefore, the gas fuel is mixed with the ambient gas more quickly at the initial time of injection than the liquid fuel is done. The experimental results also validate the usefulness of this analysis.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.2
• /
• pp.113-121
• /
• 2010

LES(Large Eddy Simulation) with a cavitation model was performed to calculate an unsteady flow for a mixed flow pump with a closed type impeller. First, the comparison between the numerical and experimental results was done to evaluate a computational accuracy. Second, the torque acting on the blade was calculated by simulation to investigate how the cavitation caused the fluctuation of torque. The absolute pressure around the leading edge on the suction side of blade surface had positive impulsive peaks in both the numerical and experimental results. The simulation showed that those peaks were caused by the cavitaion which contracted and vanished around the leading edge. The absolute pressure was predicted by simulation with -10% error. The absolute pressure around the trailing edge on the suction side of blade surface had no impulsive peaks in both the numerical and experimental results, because the absolute pressure was 100 times higher than the saturated vapor pressure. The simulation results showed that the cavitation was generated around the throat, then contracted and finally vanished. The simulated pump had five throats and cavitation behaviors such as contraction and vanishing around five throats were different from each other. For instance, the cavitations around those five throats were not vanished at the same time. When the cavitation was contracted and finally vanished, the absolute pressure on the blade surface was increased. When the cavitation was contracted around the throat located on the pressure side of blade surface, the pressure became high on the pressure side of blade surface. It caused the 1.4 times higher impulsive peak in the torque than the averaged value. On the other hand, when the cavitation was contracted around the throat located on the suction side of blade surface, the pressure became high on the suction side of blade surface. It caused the 0.4 times lower impulsive peak in the torque than the averaged value. The cavitation around the throat caused the large fluctuation in torque acting on the blade.

• 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.831-838
• /
• 2016

The disaster preventing system using vegetation has been growing in the field of coastal engineering in recent years. To analyze wave and flow fields under nonlinear interactions between tsunami and vegetation, the purpose of this study is to evaluate newly-developed 3-D numerical wave tank including energy dissipation by tsunami-vegetation interaction based on existing N-S solver with porous body model. Comparing numerical results using mean drag coefficient and dynamic drag coefficient due to Reynolds number to existing experimental results it is revealed that computed results considering the dynamic drag coefficient are in good agreement with the laboratory test results for time-domain waveform. In addition, the calculated transmission coefficients of solitary waves in various vegetation densities and incident wave heights are also in good agreement with the experimental values. This confirms the validity and effectiveness of the developed 3-D numerical wave tank with the fluid resistance by vegetation.

• Nuclear Engineering and Technology
• /
• v.38 no.2
• /
• pp.147-154
• /
• 2006

Thermal fatigue of the coolant circuits of PWR plants is a major issue for nuclear safety. The problem is especially accute in mixing zones, like T-junctions, where large differences in water temperature between the two inlets and high levels of turbulence can lead to large temperature fluctuations at the wall. Until recently, studies on the matter had been tackled at EDF using steady methods: the fluid flow was solved with a CFD code using an averaged turbulence model, which led to the knowledge of the mean temperature and temperature variance at each point of the wall. But, being based on averaged quantities, this method could not reproduce the unsteady and 3D effects of the problem, like phase lag in temperature oscillations between two points, which can generate important stresses. Benefiting from advances in computer power and turbulence modelling, a new methodology is now applied, that allows to take these effects into account. The CFD tool Code_Saturne, developped at EDF, is used to solve the fluid flow using an unsteady L.E.S. approach. It is coupled with the thermal code Syrthes, which propagates the temperature fluctuations into the wall thickness. The instantaneous temperature field inside the wall can then be extracted and used for structure mechanics computations (mainly with EDF thermomechanics tool Code_Aster). The purpose of this paper is to present the application of this methodology to the simulation of a straight T-junction mock-up, similar to the Residual Heat Remover (RHR) junction found in N4 type PWR nuclear plants, and designed to study thermal striping and cracks propagation. The results are generally in good agreement with the measurements; yet, in certain areas of the flow, progress is still needed in L.E.S. modelling and in the treatment of instantaneous heat transfer at the wall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.3
• /
• pp.8-14
• /
• 2003

The particle dispersion in the turbulent mixing layer has been numerically investigated to clarify the effect of the velocity ratio in the large-scale vortical structures. In this study the LES with subgrid-scale model is employed. The Lagrangian method to predict the particle motion is applied. The particles of 10, 50, 150, 200${\mu}m$ in mean diameter were loaded into the origin of the mixing layer. It is shown that the characteristics of flow and growth rate are strongly dependent on the variation of the velocity ratio. It is also shown the relationship between the Stokes number and the particle dispersion. As a result, in the case of St~1 the particle dispersion is faster than the diffustion of the flow field while in the cases of both St<<1 and St>>1 it is shown that the particle dispersion in lower than the diffusion of the flow filed.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.5
• /
• pp.484-491
• /
• 2017

Recent research has shown that quadrupole noise has a significant influence on the overall characteristics of flow-induced noise and on the performance of underwater appendages such as sonar domes. However, advanced research generally uses the Ffowcs Williams-Hawkings analogy without considering the quadrupole source to reduce computational cost. In this study, flow-induced noise is predicted by using an LES turbulence model and a developed formulation, called the formulation Q1As method to properly take into account the quadrupole source. The noise around a circular cylinder in an underwater environment is examined for two cases with different velocities. The results from the method are compared to those obtained from the experiments and the permeable FW-H method. The results are in good agreement with the experimental data, with a difference of less than 1 dB, which indicates that the formulation Q1As method is suitable for use in predicting quadrupole noise around underwater appendages.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.3
• /
• pp.227-235
• /
• 2014

Using a 3-D numerical scheme (LES-WASS-3D) that considered wave-structure-sandy seabed interactions in a 3-D wave field, we analyzed the wave reflection and hydraulic characteristics inside a slit caisson with two chambers in a 3-D oblique wave field. To verify the 3-D numerical analysis method suggested in this study, we compared the numerical results with existing experimental results and found good agreement. The numerical analysis revealed that a standing wave field is generated on the front side of the slit caisson due to the effect of wave reflection. For incident waves propagating perpendicular to the slit caisson, the nodes and anti-nodes of the standing wave are apparent and symmetrical. However, in an oblique wave field, as the incident wave angle decreases, the nodes and anti-nodes of the standing wave become ambiguous and unsymmetrical. It was also found that the wave reflection coefficient decreases as the incident wave angle decreases. It can be pointed out that as the incident wave angle decreases, the turbulent intensity in the chamber increases. Thereby, the increased wave energy dissipation by the increased turbulent intensity reduces the rate of wave reflection. In addition, a strong turbulent intensity generally occurs in the first chamber.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.1
• /
• pp.76-84
• /
• 2019

For decades, solitary waves have commonly been used to simulate tsunami conditions in numerical studies. However, the main component of a tsunami waveform acts at completely different spatial and temporal distributions than a solitary waveform. Thus, this study applied a 2-D numerical wave tank that included a non-reflected tsunami generation system based on Navier-Stokes equations (LES-WASS-2D) to directly simulate the run-up of a tsunami-like solitary wave on a slope. First, the waveform and velocity due to the virtual depth factor were applied to the numerical wave tank to generate a tsunami, which made it possible to generate the wide waveform of a tsunami, which was not reproduced with the existing solitary wave approximation theory. Then, to validate the applied numerical model, the validity and effectiveness of the numerical wave tank were verified by comparing the results with the results of a laboratory experiment on a tsunami run-up on a smooth impermeable 1:19.85 slope. Using the numerical results, the run-up characteristics due to a tsunami-like solitary wave on an impermeable slope were also discussed in relation to the volume ratio. The maximum run-up heights increased with the ratio of the tsunami waveform. Therefore, the tsunami run-up is highly likely to be underestimated compared to a real tsunami if the solitary wave of the approximation theory is applied in a tsunami simulation in a coastal region.