• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.123-141
• /
• 1996

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone(5°) geometry. The effective gamma(γ), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about 3 ∼ 5 percent. The heat transfer coefficient were also calculated. The results were compared with VSL results in order to validate the current numerical analysis. The results from current method were compared well VSL results ; however, not well at near nose. The proper boundary condition and grid system will be studied to improve the solution quality.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.41-48
• /
• 2013

Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.

• Journal of KSNVE
• /
• v.8 no.1
• /
• pp.122-131
• /
• 1998

Internal and external acoustic fields of the engine inlet are calculated by using a finite element method. The far fields non reflecting boundary condition is enforced by using a wave envelope element, which is a kind of infinite element. The geometry is assumed an axisymetric duct. Sources of the fan are modeled by the Tyler and Sofrin's theory. Effects of uniformly moving medium are considered. A pulsating sphere and an oscillating piston problem are calculated to verify the external problems, and compared with exact solutions. When the wave envelope element is applied at the far boundary, the calculated finite element solutions show good agreements with the exact solutions. The engine inlet is solved with the combined internal and external grid. The cut-off phenomena on engine inlet duct are observed.

• Journal of Ship and Ocean Technology
• /
• v.3 no.1
• /
• pp.1-11
• /
• 1999

Numerical simulation for the axisymmetric body with contractive and dilative periodic motion is carried out. The present analysis shows that a propulsive force can be obtained in highly viscous fluid by the contractive and dilative motion of axisymmetric body. An axisymmetric code is developed with unstructured grid system for the simulation of complicated motion and geometry. It is validated by comparing with the results of Stokes approximation with the problem of uniform flow past a sphere in low Reynolds number($R_n$ = 1). The validated code is applied to the simulation of contractive and dilative periodic motion of body whose results are quantitatively compared with the two dimensional case. The simulation is extended to the analysis of waving surface with projecting part for finding out the difference of hydrodynamics performance according to variation of waving surface configuration. The present study will be the basic research for the development of the propulsor of an axisymmetric micro-hydro-machine.

• Geomechanics and Engineering
• /
• v.18 no.6
• /
• pp.651-659
• /
• 2019

The IB-SEM numerical method combines the spectral/hp element method and the rigid immersed boundary method. This method avoids the problems of low computational efficiency and errors that are caused by the re-division of the grid when the solids move. Based on the Fourier transformation and the 3D immersed boundary method, the 3D IB-SEM system was established. Then, using the open MPI and the Hamilton HPC service, the computational efficiency was increased substantially. The flows around a cylinder and a sphere were simulated by the system. The surface of the cylinder generates vortices with alternating shedding, and these vortices result in a periodic force acting on the surface of the cylinder. When the shedding vortices enter the flow field behind the cylinder, a recirculation zone is formed. Finally, the three-dimensional pore flow was successfully investigated.

• Atmosphere
• /
• v.25 no.3
• /
• pp.529-541
• /
• 2015

Korea Institute of Atmospheric Prediction System (KIAPS) has been developing a global numerial prediction model and data assimilation system. We has implemented LETKF (Local Ensemble Transform Kalman Filter, Hunt et al., 2007) data assimilation system to NCAR CAM-SE (National Center for Atmospheric Research Community Atmosphere Model with Spectral Element dynamical core, Dennis et al., 2012) that has cubed-sphere grid, known as the same grid system of KIAPS Integrated Model (KIM) now developing. In this study, we have assimilated Global Positioning System Radio Occultation (GPS-RO) bending angle measurements in addition to conventional data within ensemble-based data assimilation system. Before assimilating bending angle data, we performed a vertical unit conversion. The information of vertical localization for GPS-RO data is given by the unit of meter, but the vertical localization method in the LETKF system is based on pressure unit. Therefore, with a clever conversion of the vertical information, we have conducted experiments to search for the best vertical localization scale on GPS-RO data under the Observing System Simulation Experiments (OSSEs). As a result, we found the optimal setting of vertical localization for the GPS-RO bending angle data assimilation. We plan to apply the selected localization strategy to the LETKF system implemented to KIM which is expected to give better analysis of GPS-RO data assimilation due to much higher model top.

• Journal of the Korea Society of Computer and Information
• /
• v.29 no.3
• /
• pp.91-98
• /
• 2024

In this paper, we present a particle-grid blending framework based on a geometric approach to efficiently represent opaque ice spheres with air bubbles. The water temperature is diffused through the grid and the air bubbles represented inside the ice through the particles. To solve the problem of previous methods that generate noisy dissolved air fields, we use levelsets to lighten the algorithm, i.e., the number of active particles and the initial amount of dissolved oxygen can be used to efficiently control the termination conditions of heat diffusion. We also extend the previous dissolved air field method, which only computes near air bubbles, to transparent regions to represent realistic ice spheres, and introduce a levelset-based approach to accurately compute the orientation of particles. As a result, the method presented in this paper is about three times faster than the existing methods and shows visually improved visualization of opaque ice spheres, which can be used in the field of representing physical virtual ice forms.

• Smart Media Journal
• /
• v.8 no.4
• /
• pp.72-79
• /
• 2019

The automated demand response (DR) program encourages consumers to participate in grid operation by reducing power consumption or deferring electricity usage at peak time automatically. However, successful deployment of the automated DR program sphere needs careful assessment of appliances load profile (ALP). To this end, the recent method estimates frequency, consistency, and peak time consumption parameters of the daily ALP to compute their potential score to be involved in the DR event. Nonetheless, as the daily ALP is subject to varying with respect to the DR time ALP, the existing method could lead to an inappropriate estimation; in such a case, inappropriate appliances would be selected at the automated DR operation that effected a consumer comfort level. To address this challenge, we propose a more proper method, in which all the three parameters are calculated using ALP that overlaps with DR time, not the total daily profile. Furthermore, evaluation of our method using two public residential electricity consumption data sets, i.e., REDD and REFIT, shows that our energy management systems (EMS) could properly match a DR target. A more optimal selection of appliances for the DR event achieves a power consumption decreasing target with minimum comfort level reduction. We believe that our approach could prevent the loss of both utility and consumers. It helps the successful automated DR deployment by maintaining the consumers' willingness to participate in the program.

• Journal of KIISE:Software and Applications
• /
• v.36 no.11
• /
• pp.936-946
• /
• 2009

A graph cuts method has recently attracted a lot of attentions for image segmentation, as it can globally minimize energy functions composed of data term that reflects how each pixel fits into prior information for each class and smoothness term that penalizes discontinuities between neighboring pixels. In previous approaches to graph cuts-based automatic image segmentation, GMM(Gaussian mixture models) is generally used, and means and covariance matrixes calculated by EM algorithm were used as prior information for each cluster. However, it is practicable only for clusters with a hyper-spherical or hyper-ellipsoidal shape, as the cluster was represented based on the covariance matrix centered on the mean. For arbitrary-shaped clusters, this paper proposes graph cuts-based image segmentation using mean shift analysis. As a prior information to estimate the data term, we use the set of mean trajectories toward each mode from initial means randomly selected in $L^*u^*{\upsilon}^*$ color space. Since the mean shift procedure requires many computational times, we transform features in continuous feature space into 3D discrete grid, and use 3D kernel based on the first moment in the grid, which are needed to move the means to modes. In the experiments, we investigate the problems of mean shift-based and normalized cuts-based image segmentation methods that are recently popular methods, and the proposed method showed better performance than previous two methods and graph cuts-based automatic image segmentation using GMM on Berkeley segmentation dataset.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.