• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.2 s.233
• /
• pp.189-196
• /
• 2005

The physical model considered here is a horizontal layer of fluid heated below and cold above with a conducting body placed at the center of the layer. The body has dimensionless thermal conductivities to the fluid of 0.1, 1 and 50. Two-dimensional solution for unsteady natural convection is obtained using an accurate and efficient Chebyshev spectral methodology for different Rayleigh numbers. Multi-domain technique is used to handle a square-shaped conducting body. The results for the case of a conducting body are also compared to those of adiabatic and neutral isothermal bodies. When the dimensionless thermal conductivity is 0.1, a pattern of fluid flow and isotherms and the corresponding time-averaged surface Nusselt number are almost the same as the case of an adiabatic body. When the dimensionless thermal conductivity is 50, a pattern of flow and isotherm and the corresponding surface and time-averaged Nusselt number are similar to those of neutral body. The results for the case of dimensionless thermal conductivity of unity are also compared to those of pure natural convection.

• Journal of applied mathematics & informatics
• /
• v.42 no.3
• /
• pp.607-623
• /
• 2024

This article studies the effects of heat generation/absorption and thermal radiation on the unsteady magnetohydrodynamic (MHD) Casson fluid flow past a vertical plate through rotating porous medium with constant temperature and mass diffusion. It is assumed that the plate temperature and concentration level are raised uniformly. For finding the exact solution, a set of non-dimensional partial differential equations is solved analytically using the Laplace transform technique. The influence of various non-dimensional parameters on the velocity are discussed, including the effects of the magnetic parameter M, heat generation/absorption Q, thermal radiation parameter R, Prandtl number Pr, Schmidt number Sc, permeability of porous medium parameter, Casson fluid parameter γ, on velocity, temperature, and concentration profiles, which are discussed through several figures. It is found that velocity, temperature, and concentration profiles in the case of heat generation parameter Q, Casson fluid parameter γ, thermal Grashof number Gr, mass Grashof number Gc, Permeability Porous medium parameter K, and time t have retarding effects. It is also seen that the magnetic field M, Thermal Radiation parameter R, Prandtl field Pr, Schmidt number Sc have reverse effects on it.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.25-34
• /
• 2008

In this paper, three-dimensional multiblade row unsteady Navier-Stokes simulations at a hot streak temperature ratio of 2.0 have been performed to reveal the effects of rotor tip clearance on the inlet hot streak migration characteristics in low pressure stage of a Vaneless Counter-Rotating Turbine. The hot streak is circular in shape with a diameter equal to 25% of the high pressure turbine stator span. The hot streak center is located at 50% of the span and the leading edge of the high pressure turbine stator. The tip clearance size studied in this paper is 2.0mm(2.59% high pressure turbine rotor height, and 2.09% low pressure turbine rotor height). The numerical results show that the hot streak is not mixed out by the time it reaches the exit of high pressure turbine rotor. The separation of colder and hotter fluid is observed at the inlet of low pressure turbine rotor. Most of hotter fluid migrates towards the rotor pressure surface, and only little hotter fluid migrates to the rotor suction surface when it convects into the low pressure turbine rotor. And the hotter fluid migrated to the tip region of the high pressure turbine rotor impinges on the leading edge of the low pressure turbine rotor after it goes through the high pressure turbine rotor. The migration of the hotter fluid directly results in very high heat load at the leading edge of the low pressure turbine rotor. The migration characteristics of the hot streak in the low pressure turbine rotor are dominated by the combined effects of secondary flow and leakage flow at the tip clearance. The leakage flow trends to drive the hotter fluid towards the blade tip on the pressure surface and to the hub on the suction surface, even partial hotter fluid near the pressure surface is also driven to the rotor suction surface through the tip clearance. Compared with the case without rotor tip clearance, the heat load of the low pressure turbine rotor is intensified due to the effects of the leakage flow. And the numerical results also indicate that the leakage flow effect trends to increase the low pressure turbine rotor outlet temperature at the tip region.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.257-268
• /
• 2014

The objective of the present study is to investigate thermal diffusion and radiation effects on unsteady MHD flow past a linearly accelerated vertical porous plate with variable temperature and also with variable mass diffusion in presence of heat source or sink under the influence of applied transverse magnetic field. The fluid considered here is a gray, absorbing/emitting radiation but a non-scattering medium. At time t > 0, the plate is linearly accelerated with a velocity $u=u_0t$ in its own plane. And at the same time, plate temperature and concentration levels near the plate raised linearly with time t. The dimensionless governing equations involved in the present analysis are solved using the closed analytical method. The velocity, temperature, concentration, skin-friction, the rate or heat transfer and the rate of mass transfer are studied through graphs in terms of different physical parameters like magnetic field parameter (M), radiation parameter (R), Schmidt parameter (Sc), Soret number (So), Heat source parameter (S), Prandtl number (Pr), thermal Grashof number (Gr), mass Grashof number (Gm) and time (t).

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.95-105
• /
• 2009

The flow fields around the HARTII rotor were numerically investigated using a viscous flow solver on adaptive unstructured meshes. An overset mesh and a deforming mesh technique were used to handle the blade motion including blade deflection, which was obtain from the HARTII experimental data. A solution-adaptive mesh refinement technique was also used to capture the rotor wake effectively. Comparison of the sectional normal force and pitching moment at 87% radial station between the two cases, with and without the blade deflection, showed that the blade loading is significantly affected by blade torsion. It was found that as the mesh was refined, the strength of tip vortex is better preserved, and the magnitude of high frequency blade loading, caused by blade-vortex interaction (BVI), is further magnified. It was also found that a proper time step size, which corresponds to the cell size, should be used to predict unsteady solutions accurately. In general, the numerical results in terms of the unsteady blade loading and the rotor wake show good agreement with the experimental data.

• International Journal of Aeronautical and Space Sciences
• /
• v.4 no.1
• /
• pp.34-44
• /
• 2003

An advanced aeroelastic analysis using a computational structural dynamics (CSD), finite element method (FEM) and computational fluid dynamics (CFD) is presented in this Paper. A general aeroelastic analysis system is originally developed and applied to realistic design problems in the transonic flow region, where strong shock wave interactions exist. The present computational approach is based on the modal-based coupled nonlinear analysis with the matched-point concept and adopts the high-speed parallel processing technique on the low-cost network based PC-clustered machines. It can give very accurate and useful engineering data on the structural dynamic design of advanced flight vehicles. For the nonlinear unsteady aerodynamics in high transonic flow region, Euler equations using the unstructured grid system have been applied to easily consider complex configurations. It is typically shown that the advanced numerical approach can give very realistic and practical results for design engineers and safe flight tests. One can find that the present study conducts a virtual flutter flight test which are usually very dangerous in reality.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.295-302
• /
• 2009

The flow in the draft tube cone of Francis turbines operated at partial discharge is a complex hydrodynamic phenomenon where an incoming steady axisymmetric swirling flow evolves into a three-dimensional unsteady flow field with precessing helical vortex (also called vortex rope) and associated pressure fluctuations. The paper addresses the following fundamental question: is it possible to compute the circumferentially averaged flow field induced by the precessing vortex rope by using an axisymmetric turbulent swirling flow model? In other words, instead of averaging the measured or computed 3D velocity and pressure fields we would like to solve directly the circumferentially averaged governing equations. As a result, one could use a 2D axi-symmetric model instead of the full 3D flow simulation, with huge savings in both computing time and resources. In order to answer this question we first compute the axisymmetric turbulent swirling flow using available solvers by introducing a stagnant region model (SRM), essentially enforcing a unidirectional circumferentially averaged meridian flow as suggested by the experimental data. Numerical results obtained with both models are compared against measured axial and circumferential velocity profiles, as well as for the vortex rope location. Although the circumferentially averaged flow field cannot capture the unsteadiness of the 3D flow, it can be reliably used for further stability analysis, as well as for assessing and optimizing various techniques to stabilize the swirling flow. In particular, the methodology presented and validated in this paper is particularly useful in optimizing the blade design in order to reduce the stagnant region extent, thus mitigating the vortex rope and expending the operating range for Francis turbines.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.167-174
• /
• 2008

In this study, flutter analyses for supercritical airfoil have been conducted in transonic region. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed static and dynamic responses of supercritical airfoil. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of cascades for fluid-structure interaction (FSI) problems. Also, flow-induced vibration (FIV) analyses for various supercritical airfoil models have been conducted. Detailed flutter responses for supercritical are presented to show the physical performance and vibration characteristics in various angle of attack.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.04a
• /
• pp.68-75
• /
• 2009

In this study, aeroelastic response analyses have been conducted for a 3D wind turbine blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responsed of wind turbine blade. Vibration analyses of rotating wind-turbine blade have been conducted using the general nonlinear finite element program, SAMCEF (Ver.6.3). Reynolds-averaged Navier-Stokes (RANS)equations with spalart-allmaras turbulence model are solved for unsteady flow problems of the rotating turbine blade model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D turbine blade for fluid-structure interaction (FSI) problems. Detailed dynamic responses and instantaneous Mach contour on the blade surfaces considering flow-separation effects are presented to show the multi-physical phenomenon of the rotating wind-turbine blade model.

• Tribology and Lubricants
• /
• v.21 no.4
• /
• pp.165-170
• /
• 2005

Engine lubrication system has very complex oil flow networks. Therefore it is difficult to evaluate and optimize the system only depending on experimental results. When we plan upgrading the engine performance the engine lubrication system must be considered in the plan. In this paper, engine lubrication system design using computational fluid dynamic analysis was studied. To do this, unsteady transient flow network analysis on the engine oil circuit system was carried out. Finally we discussed the design process in the modified engine lubrication system.