• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.5
• /
• pp.613-620
• /
• 1985

An analysis on the thermal instability of horizontal Blasius flow in the form of longitudinal vortices has been carried out by introducing the 3-dimensional spatial dependence of the disturbance quantities. The stability problem has been simplified significantly by considering the limiting case of infinite Prandtl number and by skilfully replacing the boundary conditions at infinity with the interface conditions at the edge of the thermal boundary layer (or by simply confining the thermal disturbances in the thermal boundary layer). The advantage of this approach is that the critical values marking the onset of thermal instability can be readily obtained as solutions of the eigenvalues problems formulated by a 6*6(or a 5*5) determinant. Present analysis provides reasonable explanations on the existing experimental and theoretical data. Especially, the heat transfer correlation based on the present analysis agrees well with the existing experimental data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.4
• /
• pp.390-396
• /
• 1982

The onset of longitudinal vortices in horizontal Blasius flow isothermally heated from below is studied analytically. The assumption that at the onset of thermal instability the thermal disturbances are confined within the thermal boundary layer is employed for the limiting case of large Prandtl number. Polynomial representations for the basic quantities obtained by the integral method of the boundary layer analysis have been used. Then the system of differential equations and boundary conditions for disturbance quantities is reformulated in a convenient form so that the solutions may be constructed as rapidly convergent power series. The critical buoyancy parameter G $r_{x}$ $^{*}$ /R $e^{*1.5}$ falls between 2 and 6, which is about one order of magnitude lower than the existing experimental values. It is also shown that the positions of the onset of instability can be closely predicted by the present theory.y.y.

• 한국가시화정보학회:학술대회논문집
• /
• 2005.12a
• /
• pp.117-120
• /
• 2005

This study of thermal convection uses the following geometry: a horizontal layer of fluid heated from below of solid lid at bottom and cooled from above. A variety range of thermal conductivity ratio, $\kappa$ is considered to investigate the interface temperature, $\theta_{i}$ between solid and fluid region. Periodic boundary conditions are employed in the horizontal direction to allow for lateral freedom for the convection cells. A two-dimensional solution for unsteady natural convection is obtained, using an accurate and efficient Chebyshev spectral multi-domain methodology, for different effective Rayleigh numbers, $Ra_{eff}$ varying over the range of $10^{4}\;to\;10^{7}$ in which the effective Rayleigh number is defined as $Ra{\times}<\overline{T}_{i}>$.

• Journal of Environmental Science International
• /
• v.10 no.5
• /
• pp.329-336
• /
• 2001

Characgcteristics of nocturnal boundary layer(NBL) were analyzed by the upper-air observations data using with the airsonde and pilot balloons from 1994 to 1999 in Kyungpook province. The automate weather boundary layer can become stably stratified when the surface is cooler than the air. Stable nocturnal boundary layer height were estimated from the top of surface stable layer where the vertical gradient of temperature and mixing ratio tend to zero or negative. The depth of the stable nocturnal boundary layer depended largely on the thermal effect rather than the wind effect at nighttime. The NBL was more developed on the land than on the coastal region. The stability index (bulk Richardson number) showed that the NBL was stable when the wind was weak and the vertical gradient of the temperature was strong. The heat budget in the NBL was studied by considering the effect of the radiative and the cooled by both the longwave radiative flux and the divergence of the heat flux, while NBL under the cloudy sky the longwave radiative flux played a role of the warming. It was noted that the heat was not conserved in both cases. To complete the heat budget in the NBL the warming/cooling by advection and subsidence must be considered.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.10
• /
• pp.795-803
• /
• 2015

The boundary layer developing over the suction surface of a first-stage turbine blade for power generation has been investigated in this study. For three locations selected in the region where local thermal load changes dramatically, mean velocity, turbulence intensity, and one-dimensional energy spectrum are measured with a hot-wire anemometer. The results show that the suction-surface boundary layer suffers a transition from a laminar flow to a turbulent one. This transition is confirmed to be a "separated-flow transition", which usually occurs in the shear layer over a separation bubble. The local minimum thermal load on the suction surface is found at the initiation point of the transition, whereas the local maximum thermal load is observed at the location of very high near-wall turbulence intensity after the transition process. Frequency characteristics of turbulent kinetic energy before and after the transition are understood clearly from the energy spectrum data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.80-86
• /
• 2005

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with heat addition. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, start and unstart phenomena, fluctuations of specific thrust caused by thermal choking and viscous effects. We prove that the control of separation of boundary layer is the most important key problem to prevent the thermal choking.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1189-1196
• /
• 1988

Thermal entrance lengths of turbulent tube flow for viscoelastic polymer solutions are investigated experimentally in the recirculating flow system with tubes of inside diameters 8.5mm(L/D=710) and 10.3mm(L/D=1158), respectively. In the present system, the hydrodynamic and thermal boundary layers develop simultaneously from the beginning of the test section. To provide the boundary condition of constant heat flux at the wall, the test tubes are heated directly by electricity. The polymer solution used in the current study is 1000 wppm aqueous solution of polyacrylamide(Separan AP-273). The apparent viscosity of the polymer solutions circulating in the flow system are measured by the capillary tube viscometer at regular time intervals. Thermal entrance lengths vary due to the rate of degradation. The entrance lengths of degraded polymer solutions are about 500~600 times the diameter. However, the entrance lengths of fresh polymer solutions are greater than the lengths of the test tubes used in this study suggesting that thermal entrance lengths for viscoelastic polymer solutions are greater than 1100 tube times the diameters. Friction factor is almost insensitive to the degradation, but the heat transfer $j_{H}$-factor is affected seriously by degradation. Based on the present experimental data of fresh solutions a correlation for the heat transfer $j_{H}$-factor is presented.ted.

• Journal of Environmental Science International
• /
• v.13 no.8
• /
• pp.721-731
• /
• 2004

The dispersion of suspended particulates in the coastal complex terrain of mountain-inland basin (city)-sea, considering their recycling was investigated using three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). Convective boundary layer under synoptic scale westerly wind is developed with a thickness of about I km over the ground in the west of the mountain, while a thickness of thermal internal boundary layer (TIBL) is only confined to less than 200m along the eastern slope of the mountain, below an easterly sea breeze circulation. At the mid of the eastern slop of the mountain, westerly wind confronts easterly sea breeze, which goes to the height of 1700 m above sea level and is finally eastward return flow toward the sea. At this time, particulates floated from the ground surface of the city to the top of TIBL go along the eastern slope of the mountain in the passage of sea breeze, being away the TIBL and reach near the top of the mountain. Then those particulates disperse eastward below the height of sea-breeze circulation and widely spread out over the coastal sea. Total suspended particulate concentration near the ground surface of the city is very low. On the other hand, nighttime radiative cooling produces a shallow nocturnal surface inversion layer (NSIL) of 200 m thickness over the inland surface, but relatively thin thickness less than 100m is found near the mountain surface. As synoptic scale westerly wind should be intensified under the association of mountain wind along the eastern slope of mountain to inland plain and further combine with land-breeze from inland plain toward sea, resulting in strong wind as internal gravity waves with a hydraulic jump motion bounding up to about 1km upper level in the atmosphere in the west of the city and becoming a eastward return flow. Simultaneously, wind near the eastern coastal side of the city was moderate. Since the downward strong wind penetrated into the city, the particulate matters floated near the top of the mountain in the day also moved down along the eastern slope of the mountain, reaching the. downtown and merging in the ground surface inside the NSIL with a maximum ground level concentration of total suspended particulates (TSP) at 0300 LST. Some of them were bounded up from the ground surface to the 1km upper level and the others were forward to the coastal sea surface, showing their dispersions from the coastal NSIL toward the propagation area of internal gravity waves. On the next day at 0600 LST and 0900 LST, the dispersed particulates into the coastal sea could return to the coastal inland area under the influence of sea breeze and the recycled particulates combine with emitted ones from the ground surface, resulting in relatively high TSP concentration. Later, they float again up to the thermal internal boundary layer, following sea breeze circulation.