• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Korean Journal of Optics and Photonics
• /
• v.13 no.2
• /
• pp.166-170
• /
• 2002

Reflectors for reflective LCDs must be designed to be able to control the direction of the image and the viewing angle range for optimal display. In this paper, we tried to optimize the reflector for reflective LCDs. The conventional surface structures showed many defects and problems. Therefore we made efforts to improve the surface structures of the conventional reflectors through both a theoretical approach and analytic work. As a result, we proposed surface structures which showed almost uniform reflectance within the viewing angle range specified by a user even at any oblique incidence. By using Beckmann's equations, we could check the reflective property of the proposed surface structures.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.634-639
• /
• 2003

It has been studied the flow near a rotating disk with surface topography. The system Ekman number is assumed very small, i.e., $E[{\equiv}\frac{\nu}{{\Omega}^{\ast}L^{\ast2}}]<<1$ in which $L^{\ast}$ denotes a disk radius, ${\nu}$ kinematic viscosity of the fluid and ${\Omega}^{\ast}$ angular velocity of the basic state. Disk surface has a sinusoidal topographic variation along radial coordinate, i.e., $z={\delta}cos(2{\pi}{\omega}r)$, where ${\delta}$ and ${\omega}$ are, respectively, nondimensional amplitude and wave number of the disk surface. Analytic solutions, being useful over the parametric ranges of ${\delta}{\sim}O$( $E^{1/2}$ ) and ${\omega}{\leq}O$ ( $E^{1/2}$ ), are secured in a series-function form of Fourier-Bessel type. An asymptotic behavior, when $E{\rightarrow}0$, is clarified as : for a disk with surface roughness, in contrast to the case of a flat disk, the azimuthal velocity increases in magnitude, together with the thickening boundary layer. The radial velocity, however, decreases in magnitude as the amplitude of surface waviness increases. Consequently, the overall Ekman pumping at the edge of the boundary layer remains unchanged, maintaining the constant value equal to that of the flat disk.

• Journal of Energy Engineering
• /
• v.28 no.3
• /
• pp.1-9
• /
• 2019

A major stress determining the remaining life of the tube in feedwater heater of fossil fuel power plant is hoop stress by the internal pressure. However, thermal stress due to temperature difference across the wall thickness also contributed to reduce the remaining life of the tube. Therefore, thermal loading must be considered even though the contribution of internal pressure loading to the stresses of the tube was known to be much higher than that of the thermal loading. In this study, thermal stress of the tubes in the de-superheating zone was estimated, which was generated due to the temperature difference across the tube thickness. Analytic equations were shown for determining the hoop stress and the radial stress of the tube with uniform thinning and for the temperature across the tube thickness. Accuracy and effectiveness of the analytic equations for the stresses were verified by comparing the results obtained by the analytic equations with those obtained from finite element analysis. Using finite element analysis, the stresses for eccentric thinning were also determined. The effect of heat transfer coefficient on thermal stress was investigated using series of finite element analyses with various values of heat transfer coefficient for both inner and outer surface of the tube. It was shown that the effect of heat transfer coefficient at outer surface was larger than that of heat transfer coefficient at inner surface on the thermal stress of the tube. Also, the hoop stress was larger than the radial stress for both cases of uniformly and eccentrically thinned tubes when the thermal loading was only considered without internal pressure loading.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.10
• /
• pp.3126-3134
• /
• 1996

In this paper, to examine the propagation of inclined surface crack due to frictional heating, analytic model is considered as the semi-infinite elastic body subjected to the thermo-mechanical loading of an asperity moving with a high speed. Considering the moving of frictional heat source and convection on a semi-infinite surface having inclined crack, theoretical analysis was carried out to estimate the propagation characteristics of thermo-mechanical crack. Numerical results showed that stress intensity factor $K_\prod/P_0\sqrt{c}$ is increasing with increasing velocity and frictional coefficient, inclined degree, decreasing crack length and the maximum value of it is positioned at the trailing edge. So it is shown that the propagation probability of surface crack is high at the trailing edge of contact area as increasing velocity and frictional coefficient, inclined degree, as decreasing crack length.

• Journal of The Korean Astronomical Society
• /
• v.48 no.2
• /
• pp.139-154
• /
• 2015

We carry out three-dimensional hydrodynamic simulations of the supernova remnants (SNRs) produced inside molecular clouds (MCs) near their surface using the HLL code (Harten et al. 1983). We explore the dynamical evolution and the X-ray morphology of SNRs after breaking through the MC surface for ranges of the explosion depths below the surface and the density ratios of the clouds to the intercloud media (ICM). We find that if an SNR breaks out through an MC surface in its Sedov stage, the outermost dense shell of the remnant is divided into several layers. The divided layers are subject to the Rayleigh-Taylor instability and fragmented. On the other hand, if an SNR breaks through an MC after the remnant enters the snowplow phase, the radiative shell is not divided to layers. We also compare the predictions of previous analytic solutions for the expansion of SNRs in stratified media with our onedimensional simulations. Moreover, we produce synthetic X-ray surface brightness in order to research the center-bright X-ray morphology shown in thermal composite SNRs. In the late stages, a breakout SNR shows the center-bright X-ray morphology inside an MC in our results. We apply our model to the observational results of the X-ray morphology of the thermal composite SNR 3C 391.

• Proceedings of the Korea Air Pollution Research Association Conference
• /
• 2002.04a
• /
• pp.297-298
• /
• 2002

전기집진기(ESP)는 보일러, 소각로등 많은 산업 공정에서 발생되는 입자상 물질을 제거하는데 일반적으로 사용되어 왔다. 가장 널리 쓰이는 ESP의 집진효율을 예측하기 위한 수식으로는 Deutsch-Anderson식으로 다음과 같다.(equation omitted) 여기서 η는 포집효율, Ac는 집진판의 표면적(surface area), Ve는 전기집진기 내 입자의 유효 이동속도(effective migration velocity), 그리고 Q는 단위면적을 통과하는 가스의 부피 유량이다. (중략)

• Proceedings of the KIEE Conference
• /
• 2000.07b
• /
• pp.828-830
• /
• 2000

This paper describes BIEM(Boundary Integral Equation Method) for computation of three dimensional electric field distribution and numerical method that an equivalent charge density is unknown variable. After computing numerically the surface charge distribution. the distribution of both potential and electric field are obtained. Finally, this numerical method is applied to the concentric sphere and the coaxial cylindrical model and numerical result is compared to the analytic solution.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1203-1216
• /
• 2010

The homotopy analysis method (HAM) has been applied to develop an analytic solution for the coupled radiation-convection dissipative non-gray gas flow in a non-Darcy porous medium. Results are presented for the surface shear and temperature profiles are presented to illustrate the effect of various parameters appearing in the analytical formulation. The accuracy and convergence of the method is also discussed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.8
• /
• pp.71-77
• /
• 2018

In order to estimate behavior of soil mass which is located straight up of reinforced concrete culvert, Ritz method and FEM were applied and arching effects between the soil mass and adjacent soil were considered for the analyses. Analysis results obtained from the Ritz method and finite element method were compared with analytical solution. In the case of estimating nodal forces considered in FEM, caution is needed that shear stress depending on depth from ground surface should be reflected regardless of local coordinate system. Comparing the displacements computed from Ritz method with those of the analytic solution, it is seen that as the power of assumed displacement function increases, differences between the computed displacements and those of analytic solution decreases. It seems that displacements of FEM becomes closer to those of analytical solution as the number of elements are increased. It is seen that stresses computed from the Ritz method don't get closer to those of the analytic solution as the power of assumed displacement function. Stresses from FEM become closer to those of analytic solution as the number of elements are increased. Comparing the analysis results from the Ritz method and FEM with those of analytic solution, it can be seen that FEM is more reliable than Ritz method.