• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.127-138
• /
• 1999

In this study steady state solutions of cavity flows driven by two moving walls are studied. The north and east walls of the cavity are movable where as the remaining two walls are fixed in space. Numerical experiments for three different driving schemes for moving walls are done at two different Reynolds numbers of Re=40 and 400. The first scheme is to accelerate north and east walls simultaneously. In the second one, the north wall is started first and the east wall is accelerated later. In the third one the east wall starts first. It is usually expected that all these three cases yield the same steady state solution after sufficiently long time. However, present numerical experiments show that such a usual belief is valid only when the Reynolds number is low enough (Re=40). At higher Reynolds number (Re=400), the flow develops to three different steady states depending on the history of the boundary condition change.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.7
• /
• pp.683-693
• /
• 2013

In this study, the mixing characteristics in lid-driven cavity flows were studied numerically by using a hybrid lattice Boltzmann method (HLBM). First, we compared the numerical results from single-relaxation-time (LB-SRT) and multi-relaxation-time (LB-MRT) models to examine their reliability. In most of the cavity flow, the results from both the LB-SRT and the LB-MRT models were in good agreement with those using a Navier-Stokes solver for Re=100-5000. However, the LB-MRT model was superior to the LB-SRT model for the simulation of higher Reynolds number flows having a geometrical singularity with much lesser spatial oscillations. For this reason, the LB-MRT model was selected to study the mass transport in lid-driven cavity flows, and it was demonstrated that mass transport in the fluid was activated by a recirculation zone in the cavity, which is connected from the top to the bottom surfaces through two boundary layers. Various mixing characteristics such as the concentration profiles, mean Sherwood (Sh) numbers, and velocity were computed. Finally, the detailed transport mechanism and solutions for the concentration profile in the cavity were presented.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.10
• /
• pp.779-787
• /
• 2013

The effects of various nozzle configurations on infrared signature are investigated for the purpose of analysing the infrared signature level of aircraft propulsion system. A virtual subsonic aircraft is selected and then a circular convergent nozzle, which meets the mission requirements, is designed. Convergent nozzles of different configurations are designed with different geometric profiles. Using a compressible Navier-Stokes-Fourier CFD code, an analysis of thermal flow field and nozzle surface temperature distribution is conducted. From the information of plume flow field and nozzle surface temperature distribution, IR signature of plume and nozzle surface is calculated through the narrow-band model and the RadThermIR code. Finally, qualitative information for IR signature reduction is obtained through the analysis of the effects of various nozzle configurations on IR signature.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.5
• /
• pp.43-53
• /
• 2010

This study investigates a design optimization of a rotating two-pass rectangular cooling channel with staggered arrays of pin-fins. The radial basis neural network method is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis of fluid flow and heat transfer with shear stress transport turbulent model. The ratio of the diameter to height of the pin-fins and the ratio of the streamwise spacing between the pin-fins to height of the pin-fin are selected as design variables. The optimization problem has been defined as a minimization of the objective function, which is defined as a linear combination of heat transfer related term and friction loss related term with a weighting factor. Results are presented for streamlines, velocity vector fields, and contours of Nusselt numbers, friction coefficients, and turbulent kinetic energy. These results show how fluid flow in a two-pass square cooling channel evolves a converted secondary flows due to Coriolis force, staggered arrays of pin-fins, and a $180^{\circ}$ turn region. These results describe how the fluid flow affects surface heat transfer. The Coriolis force induces heat transfer discrepancy between leading and trailing surfaces, having higher Nusselt number on the leading surface in the second pass while having lower Nusselt number on the trailing surface. Dean vortices generated in $180^{\circ}$ turn region augment heat transfer in the turning region and in the upstream region of the second pass. As the result of optimization, in comparison with the reference geometry, thermal performance of the optimum geometry shows the improvement by 30.5%. Through the optimization, the diameter of pin-fin increased by 14.9% and the streamwise distance between pin-fins increased by 32.1%. And, the value of objective function decreased by 18.1%.

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

One of the main function of Ieodo Ocean Research Station is to service the information about the weather and fishing grounds condition which are collected through calibrating convection flow and ocean current around the station. However, due to the influence of the station's structure below sea level, it is difficult to obtain the exact flow data. Therefore, it is required to research on the effect of the structure and the method to evaluate and revise the observed data. In this paper, as a basic study, it deals with the algorithm that simulate the interaction between ocean current and the station structure, followed by discussions about the way to applicate the algorithm. Through 3-dimensional computational fluid dynamics analyses (using Navier-Stokes equations with K-turbulence model), the influence of the station and submerged rocks are quantitatively evaluated, and we would suggest methods how to obtain accurate flow information from the measured rough data.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.222-232
• /
• 2011

This article examines the scale effect of the flow characteristics, resistance and propulsion performance on a 317k VLCC. The turbulent flows around a ship in both towing and self-propulsion conditions are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out in both model- and full-scale. A double-body model is applied for the treatment of free surface. An asymmetric body-force propeller is used. The speed performances including resistance and propulsion factors are obtained from two kinds of methods. One is to analyze the computational results in model scale through the revised ITTC' 78 method. The other is directly to analyze the computational results in full scale. Based on the computational predictions, scale effects of the resistance and the self-propulsion factors including form factor, thrust deduction fraction, effective wake fraction and various efficiencies are investigated. Scale effects of the streamline pattern, hull pressure and local flow characteristics including x-constant sections, propeller and center plane, and transom region are also investigated. This study presents a useful tool to hull-form and propeller designers, and towing-tank experimenters to take the scale effect into consideration.

• Journal of Korea Water Resources Association
• /
• v.38 no.3 s.152
• /
• pp.189-198
• /
• 2005

This study showed that numerical simulation can be effectively used to analyze discharge capacity according to the form and arrangement of the lock gate of a tidal power plant. For the numerical simulation, FLOW-3D with Reynolds-averaged Navier-Stokes equation as a governing equation was used. This study found that improvement of apron length and approach angle of guide wall of the lock gate causes differences in discharge capacity of $10\%$ or more. In addition, there was a difference of discharge capacity caused by the connecting structures of the drainage gate and hydraulic turbine structure and the side slope at the end of apron. This study also showed that hydraulic investigation to enhance a discharge capacity is needed when the lock gate is designed and that numerical model experiments can be a useful analysis tool to design the drainage structure, as well as the hydraulic model experiment.

• 유체기계공업학회:학술대회논문집
• /
• 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.