• 한국항공우주학회지
• /
• 제45권6호
• /
• pp.447-454
• /
• 2017

난류 유동장으로 분사되는 두가지 형태의 액체 제트 (수평분사는 디젤연료와 수직분사는 물)의 액주 분열과 미립화 현상에 관한 이상(Two-phase) 유동에 대해 3차원 LES 수치해석을 수행하였다. 기체상태의 공기 유동은 오일러리안 해법을 사용하고, 액체 제트의 액적 추적은 라그랑지안 해법을 사용하여 기체-액체간 이상유동 해석을 수행하였다. 두 종류의 확률론적 분열 모델(Stochastic breakup model)을 사용하여 액적 분열을 모사하였으며, 액체제트의 침투깊이와 액적 분포(Sauter Mean Diameter)를 실험결과와 비교하여 미세하게 분열되는 액체 제트의 분열 현상에 대해 확률론적 분열 모델링의 적합성을 제시하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2011년도 제37회 추계학술대회논문집
• /
• pp.274-278
• /
• 2011

다양한 고에너지 물질이 포함된 폭발 장치를 정확하게 해석하기 위해서는 고에너지 물질의 다양한 연소 현상과 주변을 둘러싸고 있는 고체물질에 대한 대변형 현상을 정확하게 모사하는 것이 필수적이다. 이를 위하여 본 연구에서는 다중물질 각각의 경계면을 level set 함수를 이용하여 정확하게 표현하였으며, 경계면은 ghost fluid 기법을 사용하여 나타내었다. 각각의 물질에 대하여 대변형에 적합한 구성방정식을 사용하였으며, 지배방정식을 고차의 수치해석 기법을 사용하여 해석하였다. 다양한 폭발장치 중 실험적인 데이터를 이용하여 검증이 가능한 rate stick 문제를 해석하였으며, 실험과 유사한 결과를 획득하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2009년도 춘계학술대회 논문집
• /
• pp.99-102
• /
• 2009

비정상 난류 유동장으로 분사되는 액체 제트의 액주 분열과 미립화 현상에 관한 LES를 수행하였다. 기체상태의 공기 유동 해석에 오일러리안 해법을 사용하고, 액적 추적을 위하여 라그랑지안 해법을 사용하여 기체-액체간 이상유동(two phase flow) 해석을 수행하였다. 액주의 1차 및 2차 분열이 관찰되었다. 일정한 속도로 유입되는 공기유동 중에 액체 분사 속도를 달리하여, 액체-기체 운동량 플럭스 비의 변화를 고려하여 액체 제트의 침투깊이를 조사하였으며 실험결과와 유사함을 알 수 있었다. 제트 후류에서 입자 평균직경에 대한 분석을 수행하였다.

• 설비공학논문집
• /
• 제26권1호
• /
• pp.42-47
• /
• 2014

Adequate relative humidity is needed in the common living space, for human health, and industrial space, for maintenance and efficiency. Evaporative humidifier systems generally have slow response tendencies of air loads, but relatively low initial investment and maintenance costs. Humidification efficiency in the evaporative humidifier is dependent mostly on the shapes and arrays of humidifying modules. So, we tried to apply the computational fluid dynamics to the evaporative humidifier systems, and studied the humidification efficiency of evaporative humidifiers, through the outlet air conditions.

• 한국항공우주학회지
• /
• 제34권8호
• /
• pp.47-54
• /
• 2006

본 연구에서는 Lagrangian 단계와 Remap 단계로 이루어진 Two-step Eulerian 기법에 적용되는 수치기법에 대해 살펴보고, 이를 구현하였다. 요소 값은 donor cell 과 van Leer 기법으로 이류처리 하였고, 절점 값은 Simplifed ALE [SALE] 기법으로 이류처리 하였으며, 재료 경계면은 Simple Line Interface Calculation [SLIC] 기법으로 결정하였다. 또한, 프로그램 구현 과정과 Lagrangian 단계와 Remap 단계의 연결 시 주의 깊게 고려해야할 사항들에 대해 자세히 기술하였다. 본 연구를 통해 구현된 Two-step Eulerian 코드의 검증을 위해, 순수 Lagrangian 기법, 그리고 상용코드의 결과와 비교하여 해석 결과의 정확성을 검증하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2017년도 제48회 춘계학술대회논문집
• /
• pp.451-455
• /
• 2017

ZPP와 THPP 화약의 압력 카트리지가 밀폐용기에 장착되어 연소될 때의 현상을 반응성 오일러리안-라그랑지안 이상 유동 해석 코드를 통해 모사 하였다. ZPP와 THPP는 압력 카트리지 내에서 boron nitride 판으로 격리되어있고, ZPP만 열선에 의해 직접 점화되기 때문에 THPP의 연소지연효과가 발생할 가능성이 높다. 실험을 통한 THPP의 점화지연 측정은 힘들기 때문에 기존의 연구를 통해 검증된 수치해석 코드를 통해 점화지연에 대한 케이스 스터디를 수행하고 현상학적 분석을 수행하였다. 해석 결과 THPP의 점화지연 정도에 따라 초기 충격파의 강도가 변하여 압력선도의 초기 피크특성 뿐만 아니라 주파수에도 영향을 미친다는 것을 알 수 있었다.

• 한국항공우주학회지
• /
• 제37권6호
• /
• pp.571-579
• /
• 2009

에너지 물질과 같이 연소 반응을 하는 압축성 물질을 해석하기 위하여 Hydro-SCCM (Shock Compression of Condensed Matter)이라는 에너지 물질과 비반응 물질을 포함한다중 물질 해석툴을 개발하였다. 고에너지 물질은 강한 충격파와 고온과 고압을 가진 물질경계면에서 높은 변형률을 발생시킨다. 이러한 큰 구배를 가진 현상을 해석하기 위하여 새로운 오일러리안 기법을 사용하였다. 본 논문에서는 현상을 해석하기 위한 수학적 방법과 해석결과를 소개하였다.

• 한국전산유체공학회지
• /
• 제20권2호
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.

• 한국소음진동공학회논문집
• /
• 제25권12호
• /
• pp.866-873
• /
• 2015

The transverse vibration of the deploying beam from rigid hub was analyzed by simulation and experiment. The linear governing equation of the deploying beam was obtained using the Euler-Bernoulli beam theory. To discretize the governing equation, the Galerkin method was used. After transforming the governing equation into the weak form, the weak form was discretized. The discretized equation was expressed by the matrix-vector form, and then the Newmark method was applied to simulate. To consider the damping effect of the beam, we conducted the modal test with various beam length. The mass proportional damping was selected by the relation of the first and second damping ratio. The proportional damping coefficient was calculated using the acquired natural frequency and damping ratio through the modal test. The experiment was set up to measure the transverse vibration of the deploying beam. The fixed beam at the carriage of the linear actuator was moved by moving the carriage. The transverse vibration of the deploying beam was observed by the Eulerian description near the hub. The deploying or retraction motion of the beam had the constant velocity and the velocity profile with acceleration and deceleration. We compared the transverse vibration results by the simulation and experiment. The observed response by the Eulerian description were analyzed.

• 한국환경과학회지
• /
• 제15권5호
• /
• pp.387-396
• /
• 2006

The numerical modeling and comparison with observations are performed to find out the detailed structure of meteorology and the characteristic of related dispersion phenomena of the non-reactive air pollutant at Kyoungin region, South Korea, where several industrial complex including Siwha, Banwol and Namdong is located. MM5 (Fifth Generation NCAR/Penn State Mesoscale Model), 3-D Land/sea breeze model and 3-D diagnostic meteorological model have been utilized for the meteorological simulation for September, 2002 with each different spatial resolution, while 3-D Eulerian air dispersion model for the air quality study. We can see the simulated wind field shows the very local circulation quitely well compared with in-site observations in shoreline area with complex terrains, at which the circulation of Land/sea breeze has developed and merged with the mountain and valley breeze eventually. Also it is shown in the result of the dispersion model that the diurnal variation and absolute value of daily mean $SO_2$ concentrations have good agreement with observations, even though the instant concentration of $SO_2$ simulated overestimates around 1.5 times rather than that of observation due to neglecting the deposition process and roughly estimated emission rate. This results may indicate that it is important for the air quality study at shoreline region with the complex terrain to implement the high resolution meteorological model which is able to handle with the complicate local circulation.