• Solar Energy
• /
• v.17 no.2
• /
• pp.35-42
• /
• 1997

Thermal diode is a device which allows heat to be transferred in one direction by convection due to difference of density of fluid, and blocks heat flow in the opposite direction. Vertical plate for heat collection and radiation are of utility for design of thermal diode. It was considered the transient and steady state of air filled thermal diode with guide vane which combined rectangular and parallelogrammic shape enclosures. $Gr^*$ was kept constantly on $1.11{\times}10^{10}{\sim}1.4{\times}10^{10}$ and error range was ${\pm}3%$ during the experiment. Nu was examined when inclined angle are $15^{\cir\c}\;and\;45^{\circ}$ and, also the experiments was carried out with and without guide vane as well. Specially, Nu was linearly increased due to increase of $Gr^*$, and the effect of guide vane and dimensionless channel depth was sensitive. Developed state of temperature began at dimensionless time $0.5{\sim}0.6$ due to variation of inclined angle, which is characteristic of system.

• Journal of Advanced Marine Engineering and Technology
• /
• v.35 no.4
• /
• pp.414-420
• /
• 2011

Thermal stress analysis of a planar anode-supported SOFC considering electrochemical reactions has been performed under operating conditions where average current density varies from 0 to 2000 $A/m^2$. For the case of the 2000 $A/m^2$ operating condition, Structural stress analysis based on the temperature distributions obtained from the CFD analysis of the unit cell has also been done. From this one way Fluid-Structure Interaction(FSI) analysis, Maximum Von-Mises stress under negligible temperature gradient fields occurs when cell components are perfectly bonded. The maximum stress of the electrolyte, cathode and anode in a unit cell SOFC is 262.58MPa, 28.55MPa and 15.1MPa respectively. The maximum thermal stress is critically dependent on static friction coefficient.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.12
• /
• pp.5676-5683
• /
• 2012

In this study, numerical analysis was performed on three shapes of heat transfer plates (chevron, wave and dimple type), which are currently used as fillers of cooling towers. Results show that heat transfer rates per consumed power were larger for enhanced plates as compared with that of plain plate. Highest heat transfer coefficient was obtained for wave shape followed by chevron and dimple shape. For wave shape, cross corrugations induced significant mixing of fluids, which enhanced the heat transfer. Friction factor yielded a similar trend with the heat transfer coefficient. However, heat transfer rate and pressure drop per sheet was the largest for chevron shape, due to the largest heat transfer area per sheet.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.18-28
• /
• 2001

The present study investigates flow dynamics between two dimensional compliant plates under sinusoidal flow conditions in order to understand influence of wall motion, impedance phase angle (time delay between pressure and flow waveforms), and non-Newtonian fluid on wall shear stress using computational fluid dynamics. The results showed that wall motion induced additional terms in the streamwise velocity profile and the pressure gradient. These additional terms due to wall motion reduced the amplitude of wall shear stress and also changed the mean wall shear stress. The trend of the changes was very different depending on the impedance phase angle. As the impedance phase angle was changed to more negative values, the mean wall shear stress decreased while the amplitude of wall shear stress increased. As the phase angle was reduced from 0°to -90°under $\pm$4% wall motion, the mean wall shear stress decreased by 12% and the amplitude of wall shear stress increased by 9%. Therefore, for hypertensive patients who have large negative phase angles, the ratio of amplitude and mean of the wall shear stress is raised resulting in a more vulnerable state to atherosclerosis according to the low and oscillatory shear stress theory. We also found that non-Newtonian characteristics of the blood protect atherosclerosis by decreasing the oscillatory shear index.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.9
• /
• pp.815-821
• /
• 2013

The efficiency and lifetime of a polymer electrolyte membrane fuel cell (PEMFC) system is critically affected by the humidity of the incoming gas, which should be maintained properly under normal operating conditions. Typically, the incoming gas of a fuel cell is humidified by an external humidifier, but few studies have reported on the device characteristics. In this study, a laboratory-scale planar membrane humidifier is designed to investigate the characteristics of water transport through a hydrophilic membrane. The planar membrane humidifier is immersed in a constant temperature bath to isolate the humidifier from the effect of temperature variations. The mass transfer capability of the hydrophilic membrane is first examined under isothermal conditions. Then, the mass transfer capability is investigated under various conditions. The results show that water transport in the hydrophilic membrane is significantly affected by the flow rate, operating temperature, operating pressure, and flow arrangement.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.11 s.242
• /
• pp.1265-1276
• /
• 2005

The elliptic conceptual second moment models for turbulent heat fluxes, which are proposed on the basis of elliptic-blending and elliptic-relaxation equations, are applied to calculate the combined forced and natural turbulent convection in a vertical plane channel. The models satisfy the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also have the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. Also the models are closely linked to the elliptic blending model which is used for the prediction of Reynolds stress. In order to calibrate the heat flux models, firstly, the distributions of mean temperature and scala flux in fully developed channel flow with constant wall difference temperature are solved by the present models. The buoyancy effect on the turbulent characteristics including the mean velocity and temperature, the Reynolds stress tensor, and the turbulent heat flux vector are examined. In the opposing flow, the turbulent transport is greatly enhanced with both the Reynolds stresses and the turbulent heat fluxes being remarkably increased; whereas, in the aiding flow, the opposite change is observed. The results of prediction are directly compared to the DNS to assess the performance of the model predictions and show that the behaviors of the turbulent heat transfer in the whole flow region are well captured by the present models.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.217-228
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of the p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity using the constitutive equation for work-hardening materials, and the associated flow rule. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the circular plate with uniformly distributed load. Those results are compared with the theoretical solutions and the numerical solutions of ADINA

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• Solar Energy
• /
• v.15 no.1
• /
• pp.85-94
• /
• 1995

This research has been proceeded over the potential core region(H/B=2) of two-dimensional impinging air jet sytem, in which square rods(width of 6 mm) has been set up in front of heating surface in order to increase heat transfer. The objective of this research was to investigate the characteristics of heat transfer and air flow, in eases of the clearance from rods to heating surface(C=1, 2, 4, 6 mm) and the pitch between each rods(P=30, 40, 50 mm) changed. And this research compared the above with the experimentation without rods. As result, heat transfer performance was best under the condition of C=1mm, in case clearance changed, and there was no serious difference in the effect of heat transfer augmentation in the case of pitch of rods changed.