• 유체기계공업학회:학술대회논문집
• /
• 2004.12a
• /
• pp.396-401
• /
• 2004

In this paper, a 1MW HAWT(FIL-1000) rotor blade has been designed by BEMT(Blade Element Momentum Theory) with Prandtl's tip loss. Also, a 3-D flow and performance analysis on the FIL-1000 rotor blade has been carried out by using the 3-D Navier-Stokes commercial solver (CFX-5.7) to provide more efficient design techniques to the large-scale HAWT engineers. The rated power and itsapproaching wind velocity at design point (TSR=7.5) are 1MW and 9.99m/s respectively. The rotor diameter is 54.5m and the rotating speed is 26.28rpm. Airfoils such as FFA W-301, DU91-W-250, DU93-W-210, NACA 63418, NACA 63415 consist of the rotor blade from hub to tip. Recent CFX version, 5.7 was adopted to simulate 3-D flow field and to analyze the performance characteristics of the rotor blade. Entire mesh node number is about 730,000 and it is generated by ICEM-CFD to achieve better mesh quality The predicted maximum power occurringat the design tip speed ratio is 931.45kW. Approaching to the root, the inflow angle becomes large, which causesthe blade to be stalled in the region. Therefore, k-$\omega$ SST turbulence model was used to predict the quantitative flow information more accurately. Application of commercial CFD code to optimum blade design and performance analysis was proved to be more effective environment to HAWT blade designers.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.435-449
• /
• 2014

The effects of the gravity field and the free surface on the cavity shape and the drag are investigated through a numerical analysis for the steady supercavitating flow past a simple two-dimensional body underneath the free surface. The continuity and the RANS equations are numerically solved for an incompressible fluid using a $k-{\epsilon}$ turbulence model and a mixture fluid model has been applied for calculating the multiphase flow of air, water and vapor using the method of volume of fluid and the Schnerr-Sauer cavitation model. Numerical solutions have been obtained for the supercavitating flow about a two-dimensional $30^{\circ}$ wedge in wide range of depths of submergence and inflow velocities. The results are presented for the cavity shape, especially the length and the width, and the drag of the wedge in comparison with those of the case for the infinite fluid flow neglecting the gravity and the free surface. The influences of the gravity field and the free surface on the aforementioned quantities are discussed. The length and the width of the supercavity are reduced and the centerline of the cavity rises toward the free surface due to the effects of the gravity field and the free surface. The drag coefficient of the wedge, however, is about the same except for shallow depths of submergence. As the supercavitating wedge is approaching very close to the free surface, it is found the length and the width of a cavity are shorten even though the cavitation number is reduced. Also the present result suggests that, under the influence of the gravity field and the free surface, the length of the supercavity for a certain cavitation number varies and moreover is proportional to the inverse of the submergence depth Froude number.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.2
• /
• pp.14-21
• /
• 2015

Design lifetime of a wind turbine is required to be at least 20 years. The most important step to ensure the deign is to evaluate the loads on the wind turbine as accurately as possible. In this study, extreme design load of a offshore wind turbine using Garrad Hassan (GH) Bladed and National Renewable Energy Laboratory (NREL) FAST codes are calculated considering structural dynamic loads. These wind turbine aeroelastic analysis codes are high efficiency for the rapid numerical analysis scheme. But, these codes are mainly based on the mathematical and semi-empirical theories such as unsteady blade element momentum (UBEM) theory, generalized dynamic wake (GDW), dynamic inflow model, dynamic stall model, and tower influence model. Thus, advanced CFD-dynamic coupling method is also applied to conduct cross verification with FAST and GH Bladed codes. If the unsteady characteristics of wind condition are strong, such as extreme design wind condition, it is possible to occur the error in analysis results. The NREL 5 MW offshore wind turbine model as a benchmark case is practically considered for the comparison of calculated designed loads. Computational analyses for typical design load conditions such as normal turbulence model (NTM), normal wind profile (NWP), extreme operation gust (EOG), and extreme direction change (EDC) have been conducted and those results are quantitatively compared with each other. It is importantly shown that there are somewhat differences as maximum amount of 18% among numerical tools depending on the design load cases.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2005.10a
• /
• pp.135-140
• /
• 2005

There are cofferdam and watertight wall to prevent of circulation or pollution during building of ocean structures like a dam and bridge in the harbors area and the sea. Inflow fluid and base of structure is important thing as one of the structural design factors for this interception wall like a cofferdam and watertight wall. In this study, it is revealed that at least 500 instantaneous velocity field data are required for ensemble average to get reliable turbulence statistics. The turbulent shear flow around a surface-mounted vertical wall was investigated by using the two-frame PIV(CACTUS 3.1) system and Mean velocity distributions have also been measured in the whole flow field.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.6
• /
• pp.494-502
• /
• 2016

The accurate assessment of hull-appendage interaction in the early design stage is important to control the inflow to the propeller plane, which can cause undesirable hydrodynamic effects in terms of cavitation phenomenon. This paper describes a numerical analysis for the flow around a fully appended surface ship model for which KRISO has carried out a model test in the Large Cavitation Tunnel(LCT). This numerical study was performed with the LCT model test in a complementary manner for a good reproduction of the wake distribution of surface ships. A second order accurate finite volume method provided by a commercial computational fluid dynamics(CFD) program was used to solve the governing Reynolds Averaged Navier-Stokes(RANS) equations, where the SST $k-{\omega}$ model was used for turbulence closure. The numerical results were compared to available LCT experimental data for validation. The calculations gave good predictions for the boundary layer profiles on the walls of the empty cavitation tunnel and the wake at the propeller plane of the fully appended hull model in the LCT.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.102-115
• /
• 2020

In this study, the SST k - ω turbulence model and the sliding mesh technology based on RANS method have been adopted to simulate the exciting force and hydrodynamic of a pump-jet propulsor in different oblique inflow angle (0°, 10°, 20°, 30°) and different advance ratio (J = 0.95, J = 1.18, J = 1.58).The fully structured grid and full channel model have been adopted to improved computational accuracy. The classical skewed marine propeller E779A with different advance ratio was carried out to verify the accuracy of the numerical simulation method. The grid independence was verified. The time-domain data of pump-jet propulsor exciting force including bearing force and fluctuating pressure in different working conditions was monitored, and then which was converted to frequency domain data by fast Fourier transform (FFT). The variation laws of bearing force and fluctuating pressure in different advance ratio and different oblique flow angle has been presented. The influence of the peak of pulsation pressure in different oblique flow angle and different advance ratio has been presented. The results show that the exciting force increases with the increase of the advance ratio, the closer which is to the rotor domain and the closer to the blades tip, the greater the variation of the pulsating pressure. At the same time, the exciting force decrease with the oblique flow angle increases. And the vertical and transverse forces will change more obviously, which is the main cause of the exciting force. In addition, the pressure distribution and the velocity distribution of rotor blades tip in different oblique flow angles has been investigated.

• Journal of Environmental Science International
• /
• v.19 no.11
• /
• pp.1423-1436
• /
• 2010

The effects of meteorological and reclaiming conditions on the reduction of suspended particles are investigated using a computational fluid dynamics (CFD) model with the k-$\varepsilon$ turbulence closure scheme based on the renormalization group (RNG) theory. Twelve numerical experiments with different meteorological and reclaiming conditions are performed. For identifying the meteorological characteristics of the target area and providing the inflow conditions of the CFD model, the observed data from the automatic weather station (AWS) near the target area is analyzed. Complicated flow patterns such as flow distortion, horse-shoe vortex, recirculation zone, and channeling flow appeared due to the topography and buildings in the domain. Specially, the flow characteristics around the reclamation area are affected by the reclaiming height, reclaiming size and windbreak height. Reclaiming height affected the wind speed above the reclaiming area. Windbreak induces more complicated flow patterns around the reclaiming area as well as within the reclaiming area. In front of the windbreak, flow is distorted as it impinges on the windbreak. As a result, upward flow is generated there. Behind the windbreak, a secondary circulation, so called, a recirculation zone is generated and flow is reattached at the end of the recirculation zone (reattachment point). At the lower part of the recirculation zone, there is a reverse flow toward the windbreak. Flow passing to the reattachment point starts to be recovered. Total amounts of suspended particles are calculated using the frictional and threshold frictional velocities, erosion potential function, and the number of surface disturbance. In the case of a 10 m-reclaiming and northerly wind, the amount of suspended particles is largest. In the presence of 5 m windbreak, the friction velocity above the reclaiming area is largely reduced. As a result, the total amount of the suspended particles largely decreases, compared to the case with the same reclaiming and meteorological conditions except for the windbreak The calculated suspended particle amounts are used as the emission rate of the dispersion model simulations and the dispersion characteristics of the suspended particles are analyzed.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.8 no.4
• /
• pp.357-368
• /
• 2003

Time-series CTDT (Conductivity/Temperature/Depth/Transmissivity) were obtained at one point near tidewater glacier of Marian Cove (King George Islands, Antarctica) to present water column properties and SPM (suspended particulate matter) dispersal pattern in relation with tide, current, meteorological data, and SPM concentration. Four layers were divided from the water column characteristics measured in the interval of an hour for about 2 days: 1) cold, fresh, and turbid surface mixed layer between 0-20 m in water depth, 2) warm, saline, and relatively clean Maxwell Bay inflow between 20-40 m in water depth, 3) turbid/cold tongue of subglacial discharges compared with the ambient waters between 40-70 m in water depth, and 4) cold, saline, and clean bottom water beneath 70 m in water depth. Surface plume, turbid freshwater at coastal/cliff area in late summer (early February), had the characteristic temperature and SPM concentration according to morphology, glacial condition, and composition of sediments. The restrict dispersion only over the input source of meltwater discharges was due to calm wether condition. Due to strong wind-induced surface turbulence, fresh and turbid surface plume, englacial upwelling cold water, glacier-contact meltwater, and Maxwell Bay inflow was mixing at ice-proximal zone and the consequent mixed layer deepened at the surface. Large amount of precipitation, the major controlling factor for increasing short-term glacial discharges, was accompanied by the apparent development of subglacial discharge that resulted in the rapid drop of salinity below the mid depth. Although amount of subglacial discharge and englacial upwelling may be large, however, their low SPM concentration would have small influence on bottom deposition of terrigenous sediments.

• Journal of the Institute of Electronics and Information Engineers
• /
• v.51 no.12
• /
• pp.180-188
• /
• 2014

It is very important accurate diagnosis and quick treatment in cerebrovascular disease, i.e. stenosis or occlusion that could be caused by risk factors such as poor dietary habits, insufficient exercise, and obesity. Time-of-flight magnetic resonance angiography (TOF-MRA), it is well known as diagnostic method without using contrast agent for cerebrovascular disease, is the most representative and reliable technique. Nevertheless, it still has measurement errors (also known as overestimation) for length of stenosis and area of occlusion in celebral infarction that is built by accumulation and rupture of plaques generated by hemodynamic turbulence. The purpose of this study is to show clinical trial feasibility for 3D-SPACE T2, which is improved by using signal attenuation effects of fluid velocity, in diagnosis of cerebrovascular disease. To model angiostenosis, strictures of different proportions (40%, 50%, 60%, and 70%) and virtual blood stream (normal saline) of different velocities (0.19 ml/sec, 1.5 ml/sec, 2.1 ml/sec, and 2.6 ml/sec) by using dialysis were made. Cross-examinations were performed for 3D-SPACE T2 and TOF-MRA (16 times each). The accuracy of measurement for length of stenosis was compared in all experimental conditions. 3D-SPACE 2T has superiority in terms of accuracy for measurements of the length of stenosis, compared with TOF-MRA. Also, it is robust in fast blood stream and large stenosis than TOF-MRA. 3D-SPACE 2T will be promising technique to increase diagnosis accuracy in narrow complex lesions as like two cerebral small vessels with stenosis, created by hemodynamic turbulence.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.