• Journal of The Korean Astronomical Society
• /
• v.49 no.1
• /
• pp.1-8
• /
• 2016

The aim of this study is to describe the physical processes taking place in the solar photosphere. Based on 3D hydrodynamic simulations including a detailed radiation transfer scheme, we investigate thermodynamic structures and radiation fields in solar surface convection. As a starting model, the initial stratification in the outer envelope calculated using the solar calibrations in the context of the standard stellar theory. When the numerical fluid becomes thermally relaxed, the thermodynamic structure of the steady-state turbulent flow was explicitly collected. Particularly, a non-grey radiative transfer incorporating the opacity distribution function was considered in our calculations. In addition, we evaluate the classical approximations that are usually adopted in the onedimensional stellar structure models. We numerically reconfirm that radiation fields are well represented by the asymptotic characteristics of the Eddington approximation (the diffusion limit and the streaming limit). However, this classical approximation underestimates radiation energy in the shallow layers near the surface, which implies that a reliable treatment of the non-grey line opacities is crucial for the accurate description of the photospheric convection phenomenon.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1774-1783
• /
• 1998

The transition from steady laminar to chaotic convection in a glass melting furnace specified by upper surface temperature distribution has been studied by the direct numerical analysis of the two and three-dimensional time dependent Navier-Stokes equations. The thermal instability of convection roll may take place when modified Rayleigh number($Ra_m$) is larger than $9.71{\times}10^4$. It is shown that the basic flows in a glass melting furnace are steady laminar, unsteady periodic, quasi-periodic or chaotic flow. The dimensionless time scale of unsteady period is about the viscous diffusion time, ${\tau}_d=H^2/{\nu}_0$. Through primary and secondary instability analyses the fundamental unsteady feature in a glass melting furnace is well defined as the unsteady periodic or weak chaotic flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.382-390
• /
• 1992

The problem of natural convection from a vertical fin is solved by coupling the thermal diffusion equation in the fin to the constitutive equations of the ambient medium involving the radiation of the medium. The analysis is accomplished by employing an integral method. The governing equations for the problem are solved by shooting method based on the Runge-Kutta Scheme at Pr= 0.7. For the range of values of the fin parameter and the radiation-conduction parameter in the analysis, the numerical results show that the radiation effects play an important role in the heat transfer and enhance the heat transfer.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.48-60
• /
• 2016

Accurate simulation of free-surface wave flows around a ship is very important for better hull-form design. In this paper, a computational fluid dynamics (CFD) code which is based on the open source libraries, OpenFOAM, was developed to predict the wave patterns around a ship. Additional anti-diffusion source term for minimizing a numerical diffusion, which was caused by convection differencing scheme, was considered in the volume-fraction transport equation. The influence of the anti-diffusion source term was tested by applying it to free-surface wave flow around the Wigley and KCS model ships. In results, the wave patterns and hull wave profiles of the Wigley and KCS model ships for various anti-diffusion coefficients showed quite close patterns. While, the band width of the water volume-fraction values between 0.1 to 0.9 at the Wigley and KCS model hull surfaces was narrowed by considering the anti-diffusion term. From the results, anti-diffusion source term decreased free-surface smearing.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.26-31
• /
• 2007

A comparative study on the treatment of the turbulent heat flux with the elliptic blending second-moment closure for a natural convection flow is performed. Three cases of different treating the turbulent heat flux are considered. Those are the generalized gradient diffusion hypothesis (GGDH), the algebraic flux model (AFM) and the differential flux model (DFM). The constants in the models are adjusted with a primary emphasis placed on the accuracy of predicting the local Nusselt number. These models are implemented in a computer code specially designed for evaluation of turbulent models. Calculations are performed for a turbulent natural convection in the 1:5 rectangular cavity and the calculated results are compared with the available experimental data. The results show that the three models produce nearly the same accuracy of solutions. These results show that the GGDH, AFM and DFM models for treating the turbulent heat flux are sufficient for this simple shear flow where the shear production is dominant. It is observed that, in the weakly stratified region at the center zone of the cavity, the vertical velocity fluctuation is nearly zero in the GGDH solutions, which shows that the GGDH model may not be suitable for the strongly stratified flow. Thus, further study on the strongly stratified flow should be followed.

• Journal of ILASS-Korea
• /
• v.26 no.2
• /
• pp.88-95
• /
• 2021

This study aims to investigate the influence of the droplet volume on the evaporation characteristics of the sessile droplet. In particular, the effect of the free convection in the vapor domain on the evaporation rate was analyzed through the numerical simulation. The commercial code of the ANSYS Fluent (V.2020 R2) was used to simulate the heat transfer in the liquid-vapor domain. Moreover, we used the diffusion model to estimate the evaporation rate for the different droplet volume under the room temperature. It was found that the evaporation rate significantly increases with the droplet volume because of the larger surface area for the mass transfer. Also, the effect of free convection on the evaporation rate becomes significant with an increment of droplet volume owing to the increase in the droplet radius corresponding to the characteristic length of the free convection.

• Journal of the Korean Electrochemical Society
• /
• v.11 no.1
• /
• pp.6-10
• /
• 2008

The amount of methanol crossover was measured with changing the operating condition by using a liquid methanol concentration sensor. Appropriate operating condition was discussed in terms of methanol crossover. Mechanism of methanol crossover was classified into three items which are diffusion, convection and electro-osmosis. Contribution of each mechanism to methanol crossover and the effect of operating condition were analyzed with varying methanol concentration, pressure difference between anode and cathode, current, temperature, and stoichiometry of anode fuel. Among the three mechanisms diffusion affected mostly and electro-osmosis effect was observed only under high methanol concentration.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.7 no.2
• /
• pp.75-81
• /
• 2003

Compact schemes are shown to be effective for a class of problems including convection-diffusion equations when combined with multigrid algorithms [7, 8] and V-cycle convergence is proved[5]. We apply the multigrid algorithm for an semianalytic finite difference scheme, which is desinged to preserve high order accuracy despite of singularities.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.17 no.6
• /
• pp.256-263
• /
• 2007

For $P_B=50Torr,\;P_T=5401Torr,\;T_S=450^{\circ}C,\;{\Delta}T=20K$, Ar=5, Pr=3.34, Le=0.01, Pe=4.16, Cv=1.05, adiabatic and linear thermal profiles at walls, the intensity of solutal convection (solutal Grashof number $Grs=7.86{\times}10^6$) is greater than that of thermal convection (thermal Grashof number $Grt=4.83{\times}10^5$) by one order of magnitude, which is based on the solutally buoyancy-driven convection due to the disparity in the molecular weights of the component A ($Hg_2Cl_2$) and B (He). With increasing the partial pressure of component B from 20 up to 800 Torr, the rate is decreased exponentially. It is also interesting that as the partial pressure of component B is increased by a factor of 2, the rate is approximately reduced by a half. For systems under consideration, the rate increases linearly and directly with the dimensionless Peclet number which reflects the intensity of condensation and sublimation at the crystal and source region. The convective transport decreases with lower g level and is changed to the diffusive mode at $0.1g_0$. In other words, for regions in which the g level is $0.1g_0$ or less, the diffusion-driven convection results in a parabolic velocity profile and a recirculating cell is not likely to occur. Therefore a gravitational acceleration level of less than $0.1g_0$ can be adequate to ensure purely diffusive transport.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.25-34
• /
• 2013

A numerical study of a laminar natural convection of the CuO-water nanofluid in a square cavity using the Buongiorno's nonhomogeneous model is presented. All the governing equations including the volume fraction equation are discretized on a cell-centered, non-uniform grid employing the finite-volume method with a primitive variable formulation. Calculations are performed over a range of Rayleigh numbers and volume fractions of the nanopartile. From the computed results, it is shown that both the homogeneous and nonhomogeneous models predict the deterioration of the natural convection heat transfer well with an increase of the volume fraction of nanoparticle at the same Rayleigh number, which was observed in the previous experimental studies. It is also shown that the differences in the computed results of the average Nusselt number at the wall between the homogeneous and nonhomogeneous models are very small, and this indicates that the slip mechanism of the Brown diffusion and thermophoresis effects are negligible in the laminar natural convection of the nanofluid. The degradation of the heat transfer with an increase of the volume fraction of the nanoparticle in the natural convection of nanofluid is due to the increase of the viscosity and the decrease of the thermal expansion coefficient and the specific heat. It is clarified in the present study that the previous controversies between the numerical and experimental studies are owing to the different definitions of the Nusselt number.