• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.6
• /
• pp.1-7
• /
• 2017

The atomization of a liquid into multiple droplets has many important industrial applications, including the atomization of fuels in combustion processes and coating of surfaces and particles. Ultrasonic atomizing nozzle has a transducer that receives electrical input in the form of a high frequency signal from a power generator and converts that into mechanical energy at the same frequency. Liquid is atomized into a fine mist spray using high frequency sound vibrations. In coating applications, the unpressurized, low-velocity spray reduces the amount of overspray significantly because the droplets tend to settle on the substrate, rather than bouncing off it. The spray can be controlled and shaped precisely by entraining the slow-moving spray in an ancillary air stream using specialized types of spray-shaping equipment. The desired patterns of spray can be obtained using an air stream. To simulate the water mist behavior of an ultrasonic atomizing nozzle using an air stream, the Lagrangian dispersed phase model was employed using the commercial code FLUENT. The effects of the nozzle contraction shape, water droplet size and the pneumatic pressure drop on the spray characteristics were investigated to obtain the optimal condition for coating applications.

• Nuclear Engineering and Technology
• /
• v.52 no.6
• /
• pp.1172-1179
• /
• 2020

As a Lagrangian particle method, Moving Particle Semi-implicit (MPS) method has great capability to capture interface/surface. In recent years, the multiphase flow simulation using MPS method has become one of the important directions of its developments. In this study, some key methods for multiphase flow have been introduced. The interface tension model in multiphase flow is modified to maintain the smooth of the interface and suitable for the three-phase flow. The mass transfer at immiscible liquid interface entrained by single bubble which could occur in Molten Core-Concrete Interaction (MCCI) has been investigated using this particle method. With the increase of bubble size, the height of entrainment column also increases, but the time of film rupture is slightly different. With the increase of density ratio between the two liquids, the height of entrained column decreases significantly due to the decreasing buoyancy of the denser liquid in the lighter liquid. In addition, the larger the interface tension coefficient is, the more rapidly the entrained denser liquid falls. This study validates that the MPS method has shown great performance for multiphase flow simulation. Besides, the influence of physical parameters on the mass transfer at immiscible interface has also been investigated in this study.

• Journal of Environmental Science International
• /
• v.26 no.2
• /
• pp.183-200
• /
• 2017

The classification of airflow patterns during high ozone ($O_3$) and $PM_{10}$ episodes on Jeju Island in recent years (2009-2015), as well as their correlation with meteorological conditions according to classified airflow patterns were investigated in this study. The airflow patterns for $O_3$ and $PM_{10}$ were classified into four types (Types A-D) and three types (Types E-G), respectively, using the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model and synoptic weather charts. Type A was the most dominant airflow pattern for $O_3$ episodes, being characterized by the transport of airflows from urban and industrial areas in China with the highest frequency (about 69%, with a mean of 67 ppb). With regard to the $PM_{10}$ episodes, Type E was the most dominant airflow pattern, and was mostly associated with long distance transport from Asian dust source regions along northwesterly winds, having the highest frequency (about 92%, with a mean of $136{\mu}g/m^3$). The variations in the concentration of $O_3$ and $PM_{10}$ during the study period were clarified in correlation with two pollutant and meteorological variables; for example, the high (low) $O_3$ and $PM_{10}$ concentrations with high (low) air temperature and/or wind speed and vice versa for precipitation. The contribution of long-range transport to the observed $PM_{10}$ levels in urban sites for different airflow patterns (Types E-F), if estimated in comparison to the data from the Gosan background site, was found to account for approximately 87-93% (on average) of its input. The overall results of the present study suggest that the variations in $O_3$ and $PM_{10}$ concentrations on Jeju Island are mainly influenced by the transport effect, as well as the contribution of local emissions.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Economic and Environmental Geology
• /
• v.14 no.2
• /
• pp.93-102
• /
• 1981

Studies of model potential fields continued upward and downward show differences depending on the method of continuation. Beginning with a magnetic field computed over a buried vertical cylinder, the field was continued to various levels by a method introduced by Henderson (Lagrangian interpolation) and by a spectral method (frequency domain analysis). Resultant fields show (1) no significant differences in upward continued values, (2) in downward continuation, accurate values are obtained with the spectral method over the central part of the anomaly, and (3) accurate values are obtained with Henderson's method on the flanks of the anomaly, while oscillations usually characterize the spectral method in this region. Essentially the same observations are made for derivative calculations. Field oscillations are empirically predicted at levels continued to approximately two-thirds of the depth of the source. Our spectral computer program output yields marked oscillations at one-half of the depth of the source. Henderson's method shows no oscillations at this depth and only minor oscillations at the top of the body (some negative values appear on the flanks of the anomaly). The Henderson output is a smooth field even if continued below the top of the body. These results suggest that the presence of oscillations cannot be used to identify the top of a buried source without careful consideration of the method used to continue the field. Use of the derivative to outline and isolate anomalies must similarly include consideration of the method of calculation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.731-736
• /
• 2004

강우-유출과정의 수문학적 현상을 보다 정확히 분석하고 예측하는 기법으로 강우에 의한 유출의 반응을 나타내는 지체시간, 도달시간 등 수문학적 반응시간을 유역의 지형형태학적 인자들과 연계하는 방법이 많이 이용되고 있다. 이에 본 연구에서는 Clark방법과 지형형태학적 순간단위도(GIUH)를 이용하여 계측유역의 강우-유출반응을 모의하였고, 이를 관측된 값과 비교하여 미계측유역의 적용성 여부를 검토해보았다. 대상 유역의 하상지형인자 및 지형형태학적 특성은 Arc-View를 이용하여 구하였으며, 이를 기존의 문헌자료와 비교해보았다. Clark방법의 매개변수의 결정에 있어서 시간-면적곡선은 HEC-1의 무차원 식을 이용하였고, 도달시간은 Kirpich 공식을 이용하여 구하였으며, 저류상수는 Clark방법에 의해 추정된 순간단위도의 첨두유량이 Horton의 차수비의 함수로 구한 철두유량과 같아지는 값으로 결정하였다. 본 연구는 전적비교를 출구점으로하는 유역면적 $8.5km^2$인 설마천을 대상유역으로 하였으며, 모의된 강우-유출반응과 비교하기 위해 사용된 강우사상은 2002년의 8월 4일과 2002년 10월 6일의 10분 단위 우량이다. Clark방법과 GIUH를 이용하여 모의한 유출곡선과 관측된 유출곡선을 비교해본 결과 첨두유량은 8월의 강우사상 때는 $21\%$크게, 10월의 강우사상 때는 $35\%$작게 나타났다. 첨두시간은 모의된 경우가 각각 10분, 20분 빨리 도달하였다. 또한 이러한 결과는 유역의 도달시긴에 가장 민감하게 반응함을 알 수 있었다. 따라서, 유역의 도달시간 산정에 주의를 요한다면 프랙탈 차원이 유사한 미계측유역의 수문곡선 산정에 있어서 Clark방법과 GIUH를 이용하는 방법도 유용하다고 사료된다. 주는 각 수문인자 중 강우시간분포와 유효우량 산정방법 그리고 유출모형에 대해 자자 검토하였으며, 최종적으로 면적에 따른 임계지속기간과 유출량의 변화를 검토해 보았다.이를 각각의 경우의 해석해 결과와 비교${\cdot}$분석하였다. 후방추적 퍼프모형은 전방추적 퍼프모형에 비하여 사용된 퍼프수와 관계없이 작은 오차를 발생하였으며, 전체적으로 퍼프 모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 향후, 보다 다양한 흐름영역에서 장${\cdot}$단점 분석 및 오차해석을 수행한 후에 각각의 Lagrangian 모형의 장점만을 갖는 모형결합 방법을 제시할 수 있을 것으로 판단된다.mm/$m^{2}$로 감소한 소견을 보였다. 승모판 성형술은 전 승모판엽 탈출증이 있는 두 환아에서 동시에 시행하였다. 수술 후 1년 내 시행한 심초음파에서 모든 환아에서 단지 경등도 이하의 승모판 폐쇄 부전 소견을 보였다. 수술 후 조기 사망은 없었으며, 합병증으로는 유미흉이 한 명에서 있었다. 술 후 10개월째 허혈성 확장성 심근증이 호전되지 않아 Dor 술식을 시행한 후 사망한 예를 제외한 나머지 6명은 특이 증상 없이 정상 생활 중이다 결론: 좌관상동맥 페동맥이상 기시증은 드물기는 하나, 영유아기에 심근경색 및 허혈성 심근증 또는 선천성 승모판 폐쇄 부전등을 초래하는 심각한 선천성 심질환이다. 그러나 진단 즉시 직접 좌관상동맥-대동맥 이식술로 수술적 교정을 해줌으로써 좋은 성적을 기대할 수 있음을 보여주었다.특히 교사들이 중요하

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.537-542
• /
• 2004

오염물질의 이송${\cdot}$확산 과정을 해석하기 위하여 계산효율이 높은 3차원 퍼프모형을 개발하였다. 본 연구에서 개발된 퍼프모형은 추적방법에 따라 전방추적모형과 후방추적모형으로 나눌 수 있으며, 이 두 가지 추적방법은 계산효율과 수치오차에 있어서 상이한 특성을 나타내었다. 전방추적 퍼프모형은 일정한 시간간격을 가지므로 정상상태의 연속오염원의 경우에 각각의 퍼프들이 동일한 질량을 갖는다. 그러므로 전방추적 퍼프 모형은 각 퍼프들간의 중첩정도가 일정하지 않다. 이에 관한 오차분석을 수행하기 위하여 본 연구에서는 무차원 퍼프중첩계수를 정의하였다. 퍼프중첩계수란 퍼프의 크기에 대하여 퍼프중심간의 거리가 떨어진 정도를 나타내는 무차원수로서 너무 작은 경우에는 정확도가 떨어지고 너무 큰 경우에는 계산효율이 감소한다. 전방 추적 퍼프모형의 경우, 중첩계수가 작은 초기구간에는 정확도가 떨어지며, 시간이 지남에 따라 중첩계수가 필요이상으로 증가하여 계산효율이 떨어지는 것으로 나타났다. 이에 비해 일정한 중첩정도를 갖는 후방추적 퍼프모형의 경우에는 전 영역에 걸쳐서 정확도와 계산효율이 높은 것으로 나타났다. 하지만 일정한 시간간격 마다 농도장을 계산하고자 할 때, 전방추적법은 단 한번의 전체계산을 통하여 수행가능하지만 후방추적법의 경우에는 매 출력시간마다 초기시점까지 반복해서 계산해야하는 단점이 있다. 경계처리에 있어서 입자추적모형과 상이한 방법을 사용하는 퍼프모형은 폐경계에서는 입자추적모형과 동일한 결과를 나타내지만 개경계에서는 다른 결과를 나타내었다. 또한 오염원이 임의의 공간적 분포를 갖는 경우, 퍼프모형은 입자추적모형보다는 적은 수의 퍼프를 사용할 수 있지만 이에 따른 경계면에서의 수치오차를 발생하였다. 본 연구에서는 개발된 퍼프모형을 검증하고 장${\cdot}$단점을 분석하기 위하여 흐름이 일정한 경우와 전단흐름의 경우에 대하여 이송${\cdot}$확산 수치모의를 수행하였으면, 이를 각각의 경우의 해석해 결과와 비교${\cdot}$분석하였다. 후방추적 퍼프모형은 전방추적 퍼프모형에 비하여 사용된 퍼프수와 관계없이 작은 오차를 발생하였으며, 전체적으로 퍼프 모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 향후, 보다 다양한 흐름영역에서 장${\cdot}$단점 분석 및 오차해석을 수행한 후에 각각의 Lagrangian 모형의 장점만을 갖는 모형결합 방법을 제시할 수 있을 것으로 판단된다.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.40 no.5
• /
• pp.91-99
• /
• 1998

The widespread use of thin shell structures has created a need for a systematic method of analysis which can adequately account for arbitrary geometric form and boundary conditions as well as arbitrary general type of loading. Therefore, the stress and analysis of thin shell has been one of the more challenging areas of structural mechanics. A wide variety of numerical methods have been applied to the governing differential equations for spherical and cylindrical structures with a few results applicable to practice. The analysis of axisymmetric spherical shell is almost an every day occurrence in many industrial applications. A reliable and accurate finite element analysis procedure for such structures was needed. Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometry changes on the response is also significant in many cases. Therefore both material and geometric nonlinear effects should be considered. In general, the shell structures designed according to quasi-static analysis may fail under conditions of dynamic loading. For a more realistic prediction on the load carrying capacity of these shell, in addition to the dynamic effect, consideration should also include other factors such as nonlinearities in both material and geometry since these factors, in different manner, may also affect the magnitude of this capacity. The objective of this paper is to demonstrate the dynamic characteristics of spherical shell. For these purposes, the spherical shell subjected to uniformly distributed step load was analyzed for its large displacements elasto-viscoplastic static and dynamic response. Geometrically nonlinear behaviour is taken into account using a Total Lagrangian formulation and the material behaviour is assumed to elasto-viscoplastic model highly corresponding to the real behaviour of the material. The results for the dynamic characteristics of spherical shell in the cases under various conditions of base-radius/central height(a/H) and thickness/shell radius(t/R) were summarized as follows : The dynamic characteristics with a/H. 1) AS the a/H increases, the amplitude of displacement in creased. 2) The values of displacement dynamic magnification factor (DMF) were ranges from 2.9 to 6.3 in the crown of shell and the values of factor in the mid-point of shell were ranged from 1.8 to 2.6. 3) As the a/H increases, the values of DMF in the crown of shell is decreased rapidly but the values of DMF in mid-point shell is increased gradually. 4) The values of DMF of hoop-stresses were range from 3.6 to 6.8 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.3 to 2.6, and the values of DMF of stress were larger than that of displacement. The dynamic characteristics with t/R. 5) With the thickness of shell decreases, the amplitude of the displacement and the period increased. 6) The values of DMF of the displacement were ranged from 2.8 to 3.6 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.1 to 2.2.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Nuclear Engineering and Technology
• /
• v.50 no.5
• /
• pp.665-672
• /
• 2018

A venturi scrubber is an important element of Filtered Containment Venting System (FCVS) for the removal of aerosols in contaminated air. The present work involves computational fluid dynamics (CFD) study of dust particle removal efficiency of a venturi scrubber operating in self-priming mode using ANSYS CFX. Titanium oxide ($TiO_2$) particles having sizes of 1 micron have been taken as dust particles. CFD methodology to simulate the venturi scrubber has been first developed. The cascade atomization and breakup (CAB) model has been used to predict deformation of water droplets, whereas the Eulerian-Lagrangian approach has been used to handle multiphase flow involving air, dust, and water. The developed methodology has been applied to simulate venturi scrubber geometry taken from the literature. Dust particle removal efficiency has been calculated for forced feed operation of venturi scrubber and found to be in good agreement with the results available in the literature. In the second part, venturi scrubber along with a tank has been modeled in CFX, and transient simulations have been performed to study self-priming phenomenon. Self-priming has been observed by plotting the velocity vector fields of water. Suction of water in the venturi scrubber occurred due to the difference between static pressure in the venturi scrubber and the hydrostatic pressure of water inside the tank. Dust particle removal efficiency has been calculated for inlet air velocities of 1 m/s and 3 m/s. It has been observed that removal efficiency is higher in case of higher inlet air velocity.