• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.

• Structural Engineering and Mechanics
• /
• v.65 no.4
• /
• pp.453-464
• /
• 2018

In this paper, an efficient higher-order shear deformation theory is presented to analyze thermomechanical bending of temperature-dependent functionally graded (FG) plates resting on an elastic foundation. Further simplifying supposition are made to the conventional HSDT so that the number of unknowns is reduced, significantly facilitating engineering analysis. These theory account for hyperbolic distributions of the transverse shear strains and satisfy the zero traction boundary conditions on the surfaces of the plate without using shear correction factors. Power law material properties and linear steady-state thermal loads are assumed to be graded along the thickness. Nonlinear thermal conditions are imposed at the upper and lower surface for simply supported FG plates. Equations of motion are derived from the principle of virtual displacements. Analytical solutions for the thermomechanical bending analysis are obtained based on Fourier series that satisfy the boundary conditions (Navier's method). Non-dimensional results are compared for temperature-dependent FG plates and validated with those of other shear deformation theories. Numerical investigation is conducted to show the effect of material composition, plate geometry, and temperature field on the thermomechanical bending characteristics. It can be concluded that the present theory is not only accurate but also simple in predicting the thermomechanical bending responses of temperature-dependent FG plates.

• Geomechanics and Engineering
• /
• v.22 no.4
• /
• pp.361-374
• /
• 2020

This paper is concerned with the thermoelastic bending of FG beams resting on two-layer elastic foundations. One of these layers is Winkler springs with a variable modulus while the other is considered as a shear layer with a constant modulus. The beams are considered simply supported and subjected to thermo-mechanical loading. Temperature-dependent material properties are considered for the FG beams, which are assumed to be graded continuously across the panel thickness. The used theories contain undetermined integral terms which lead to a reduction of unknowns functions. Several micromechanical models are used to estimate the effective two-phase FG material properties as a function of the particles' volume fraction considering thermal effects. Analytical solutions for the thermo-mechanical bending analysis are obtained based on Navier's method that satisfies the boundary conditions. Finally, the numerical results are provided to reveal the effect of explicit micromechanical models, geometric parameters, temperature distribution and elastic foundation parameters on the thermoelastic response of FG beams.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1675-1683
• /
• 1995

A numerical study of the fluid flow and heat transfer characteristics of the two-dimensional impingement jet with a confinement plate has been carried out. The fluid flow was calculated by solving the full Navier-Stokes equation. In doing that, the well known SIMPLER algorithm was used and the trouble making convection term was discretized according to QUICKER scheme. The energy equation was simply solved by using the SOR method. For the Reynolds number of 10000, two channel heights, say 1.5 and 3.0 times the jet exit width, and two thermal boundary conditions constant wall temperature and constant wall heat flux were considered. Discrete heat sources were flush mounted along the impingement plate at a distance of 0, 2, 3, 4, 5, 6, 10, 12, times the jet exit width from the stagnation point. The length of each heat source is 4 times the jet exit width long. The Nusselt number averaged over each heat source was compared with experiment. Comparison shows that both calculations and experiment have the secondary peak of Nusselt number at downstream of stagnation point, even though there is a little quantitative difference in between. The difference is believed due to abscure thermal boundary condition in experiment and also accuracy of turbulence model used. The secondary peak is shown to be caused by rigorous turbulent flow motion generated as the wall jet flow is retarded and developes into the channel flow without flow reversal.

• 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 the Korean Society for Marine Environment & Energy
• /
• v.12 no.1
• /
• pp.9-14
• /
• 2009

A numerical investigation was performed based on the Reynolds-Averaged Navier-Stokes(RANS) equations for the two-dimensional unsteady flow around a vertical axis turbine(VAT) with three or four blades. VAT is one of the promising devices for tidal current energy conversion. The geometry of the turbine blade was $NACA65_3$-018 airfoil, for which CFD analysis using Fluent was carried out at several angles of attack and the results were compared with the corresponding experimental data for validation and calibration. Then CFD simulations were carried out for the whole vertical axis turbine with a two-dimensional setup. The CFD simulation demonstrated the usefulness of the method to study the typical unsteady flows around VATs and the results showed that the optimum turbine efficiency could be achieved for carefully selected combinations of the number of blade and Tip-Speed Ratio(TSR).

• Tribology and Lubricants
• /
• v.33 no.6
• /
• pp.288-295
• /
• 2017

In order to improve the efficiency and reliability of the machine, the friction should be minimized. The most widely used method to minimize friction is to maintain the fluid lubrication state. However, we can reduce friction only up to a certain limit because of viscosity. As a result of several recent studies, surface texturing has significantly reduced the friction in highly sliding machine elements, such as mechanical seals and thrust bearings. Thus far, theoretical studies have mainly focused on isothermal/iso-viscous conditions and have not taken into account the heat generation, caused by high viscous shear, and the temperature conditions on the bearing surface. In this study, we investigate the effect of dimple depth and film-temperature boundary conditions on the thermohydrodynamic (THD) lubrication of textured parallel slider bearings. We analyzed the continuity equation, the Navier-Stokes equation, the energy equation, and the temperature-viscosity and temperature-density relations using a computational fluid dynamics (CFD) code, FLUENT. We compare the temperature and pressure distributions at various dimple depths. The increase in oil temperature caused by viscous shear was higher in the dimple than in the bearing outlet because of the action of the strong vortex generated in the dimple. The lubrication characteristics significantly change with variations in the dimple depths and film-temperature boundary conditions. We can use the current results as basic data for optimum surface texturing; however, further studies are required for various temperature boundary conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

• 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 the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.281-289
• /
• 2016

Accurate prediction of the trajectory and time of a time-varying mass parachute system remains essential in the mission requiring a precision airdrop to the ground. In this study, we investigate the altitude-varying behavior of a cross-type parachute system designed to deliver a time-varying mass object like flare. The dynamics of the descending parachute system was analyzed based on the Runge-Kutta method of the ordinary differential system. The drag coefficients of the cross-type parachute and flare were calculated by a CFD code based on the incompressible Navier-Stokes equation. Finally, by using a simplified gust wind model in troposphere, the combined effects of gust wind and time-varying mass were examined in detail.