• 한국수자원학회논문집
• /
• 제41권12호
• /
• pp.1187-1198
• /
• 2008

국내 많은 저수지들과 그 하류하천은 높은 탁도를 가진 물의 장기적인 방류로 인해 수자원 이용과 수생태계 관리에 많은 어려움을 겪고 있다. 탁도($C_T$)는 물의 탁한 정도를 나타내는 척도로써, 수질과 수환경의 건강을 평가하는 매우 중요한 지표로써 광범위하게 사용되어 왔지만, 모형의 검증에 필요한 실험 자료의 부족으로 인해 탁도 모델링에 대한 연구는 지금까지 매우 부족하였다. 본 연구의 목적은 광범위한 현장 실측자료를 이용하여 성층화된 대청호로 유입한 탁수의 밀도류 거동 모의를 위한 3차원 수리-입자동력학 연동 수치모형인 ELCOM-CAEDYM의 적용성을 검증하는데 있다. 입자크기에 따라 구분된 3개 그룹의 부유물질 (SS) 농도가 모형의 모의변수로 사용되었으며, 모형 변수인 SS와 저수지내 실측값의 $C_T$의 변환을 위해 저수지 지점별로 측정한 SS-$C_T$ 상관관계를 사용하였다. 모의결과는 2004년에 대청호의 회남과 댐앞 지점에서 수심별로 실측한 수온과 탁도 자료와 비교함으로써 검증하였다. 모형은 저수지의 성층구조, 탁수를 포함한 하천 밀도류의 시간에 따른 진행과정을 잘 재현하였으며, 탁도의 수직분포와 크기도 실측값과 부합하였다. 본 연구에서 제시한 3차원 수치모형과 탁도 모델링 방법론은 유사한 탁도 문제를 가지고 있는 다른 저수지에서도 탁수의 최적관리를 위한 지원 도구로써 사용 가능하다고 사료된다.

• 한국전산구조공학회논문집
• /
• 제25권4호
• /
• pp.315-321
• /
• 2012

본 연구에서는 분자동역학 전산모사와 유한요소해석 기반의 균질화 기법을 통해 나노복합재의 열전도 특성을 정확하고 효율적으로 예측할 수 있는 순차적 멀티스케일 균질화 해석기법을 제안하였다. 나노입자의 크기효과가 나노복합재의 유효 열전도 특성에 미치는 영향을 조사하기 위해 크기가 다른 구형 나노입자가 첨가된 나노복합재의 열전도 계수를 분자동역학 전산모사를 통해 예측했고, 그 결과 나노입자의 크기가 작아질수록 계면에서의 Kapitza열저항에 의해 나노복합재의 열전도 계수가 점차 감소하는 것으로 나타났다. 이러한 나노입자의 크기효과를 균질화 해석모델을 통해 정확하게 묘사하기 위해 Kapitza 열저항에 의한 계면에서의 온도 불연속 구간과 고분자 기지가 높은 밀도를 가지며 흡착되는 유효계면을 추가적인 상으로 도입하여 나노복합재를 입자, Kapitza 계면, 유효계면, 기지로 구성된 4상의 연속체 구조로 모델링하였다. 이후 순차적 멀티스케일 균질화 해석기법을 통해 유효계면의 열전도 계수를 나노복합재의 열전도 계수로부터 역으로 예측했으며, 이를 입자의 반경에 대한 함수로 근사하였다. 근사 함수를 토대로 다양한 입자 체적분율과 반경에 대한 나노복합재의 유효 열전도 특성을 예측하였으며, 유효계면에 대한 매개변수 연구를 수행하였다.

• 한국음향학회지
• /
• 제34권1호
• /
• pp.36-45
• /
• 2015

유동광대역소음을 효율적으로 예측하기 위하여 통계적으로 난류를 재생하는 방법에 대한 많은 연구들이 최근에 진행되고 있다. 그 중에서도, FRPM(Fast Random Particle Mesh) 기법은 RANS(Reynolds-Averaged Navier-Stokes) 방정식 해석을 통해 도출된 정상상태 유동장의 난류 운동에너지와 소산 값을 이용하여 특정한 통계적 특성을 가지는 난류를 재생하는 기법으로서 유동광대역소음 문제 등에 성공적인 적용 예에 대해서 보고되고 있다. 하지만 기존의 FRPM 방법은 축류팬과 같이 축 대칭 특성을 갖는 기계의 경우 정상상태의 유동장을 기초로 광대역소음을 예측하는 문제에는 적용할 수 있으나, 원심팬과 같이 볼루트 영역으로 인하여 축 대칭이 성립되지 않는 기계류의 유동광대역소음에는 적용할 수 없다. 본 연구에서는 이러한 FRPM 기법을 확장하여, 원심팬에서 발생하는 광대역소음을 효율적으로 예측하기 위하여 비정상 RANS 방정식의 수치해와 연계하여 광대역소음원으로 고려되는 난류를 특정한 통계적 특성을 가지도록 재생할 수 있는 U-FRPM(Unsteady-FRPM) 기법을 제안하였다. 먼저 전산유체역학을 사용하여 RANS 방정식을 해석함으로써, 원심팬 주위의 비정상상태 유동장 정보를 도출하고, 음향상사법(Acoustic Analogy)을 기초로 도출된 유동소음원을 U-FRPM을 이용하여 모델링하였다. 모델링된 소음원은 경계요소법을 통해 구현되는 선형음향전파모델과 연계하여 수음점에서 광대역소음을 예측하는데 이용되었다. 예측된 결과와 실험결과의 비교를 통해 본 논문에서 제시한 방법의 유효성을 확인하였다.

• 한국산학기술학회논문지
• /
• 제20권8호
• /
• pp.449-455
• /
• 2019

최근 도시화로 인한 불투수면의 증가는 지표 유출수를 집수하여 배수하는 기존의 배수시스템의 부담을 증가시킨다. 이러한 방식의 우수배수시스템은 표면유출수와 함께 이송되는 각종 쓰레기, 낙엽, 등의 부유물질에 의해 배수면적이 감소하는 구조적 한계를 가지고 있다. 이러한 한계를 극복하기 위해 최근 새로운 형태의 배수시스템이 개발 및 적용되고 있다. 본 연구는 3차원 전산 유체역학 프로그램 중 하나인 ANSYS CFX를 이용하여 투수성 포장 하부에 위치한 배수층의 체적 공극과 투수도 결정을 위한 모의를 각각 수행하였다. 모의결과 35% 체적공극을 가진 배수층의 배수구 유속이 20%, 50%에 해당하는 배수층보다 큰 값을 보여 체적공극과 배수성능 사이의 상관관계는 없는 것으로 파악되었다. 투수도는 구성물질의 입경에 따라 결정되며, 5가지 조건을 모의하여 배수구 유속을 분석한 결과 입경 2 mm의 사질토가 사용성과 시공성 측면에서 가장 적절하다고 분석되었다. 본 연구는 배수층의 적절한 체적 공극과 구성물질의 입경을 제시하였고 이러한 조건을 갖는 배수층이 침수피해를 저감 및 방지의 측면에서 유리할 것으로 판단된다.

• Composites Research
• /
• 제35권6호
• /
• pp.447-455
• /
• 2022

현재 해양 미세플라스틱에 의한 해양 환경 오염문제가 국제적으로 심각하게 대두되고 있다. 이와 관련하여 본 연구에서는 휴대용 미세플라스틱 수거 장비의 경량화 개발을 위해 대한민국 전국 21개소 해안가에서의 미세플라스틱종류와개수를조사하였고, 미세플라스틱수거장비진공장비의내구성, 내부식성, 경량화를위해 CFRP(Carbon Fiber Reinforced Plastic), GFRP(Glass Fiber Reinforced Plastic), ABS(Acrylonitrile Butadiene Styrene copolymer), aluminum과 같은 다양한 소재를 적용하여 설계 및 해석을 수행하였다. 해안가에서의 시료 채취 및 분류결과 함덕해수욕장에서 미세플라스틱이 가장 많이 분포되어 있는 것을 확인하였고, 주로 폴리스타이렌(Polystyrene)의 미세플라스틱이 분포되어 있는 것으로 나타났다. 해안가 현장 조사를 통한 입자정보 분석과 전산유체해석을 통해 모래 및 불순물 등의 입자 유동속도 및 분포를 분석하였고, 이를 진공 본체 내부의 중요 부품인 싸이클론 장치의 구조해석에 적용하였다. 싸이클론 장치의 모래 및 불순물 흡입 시 발생되는 입자 충격을 고려한 구조해석 결과 CFRP, GFRP, aluminum, ABS의 순으로 구조적 안전성이 우수한 것으로 평가되었다. 경량화 측면으로는 싸이클론 장치에 aluminum 소재를 적용했을 때와 비교하여 CFRP는 53%, GFRP는 47%, ABS는 61%의 경량화가 가능한 것으로 나타났다.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2010년도 춘계학술대회 초록집
• /
• pp.106.1-106.1
• /
• 2010

석탄가스화기술은 매장량이 풍부하여 안정적인 공급이 보장되는 석탄을 이용함과 동시에 환경오염물질 감소라는 사회적 요구조건을 충족시키면서 화학제품, 석탄-가스화, 석탄-디젤화, 연료전지, 복합발전 등 다양한 분야에 응용이 가능한 장점이 있다. 특히 석탄가스화복합기술(Intergrated Coal Gasification Combined Cycle, IGCC)은 석탄을 고온, 고압하에서 가스화시켜 일산화탄소(CO), 수소($H_2$)가 주성분인 합성가스를 제조, 정제 후 가스터빈 및 증기터빈을 복합으로 구동하여 전기를 생산하는 친환경 차세대 발전기술로 주목을 받고 있다. 현재 IGCC 기술은 세계적으로 볼 때 상용화단계에 있고, 우리나라의 경우 한국형 IGCC 기술의 확보를 위한 연구사업이 진행중에 있다. 본 연구는 IGCC 발전플랜트의 발전효율을 결정하는 가장 중요한 부분이라 할 수 있는 가스화반응기의 모델링 기술을 개발하는 목적으로 진행되었다. 본 연구에서는 석탄가스화 반응기에서 발생하는 석탄의 휘발화와 Char의 표면반응 그리고 기상에서의 가스화반응등의 현상을 전산유체역학(Computational Fluid Dynamics)을 이용하여 모델링하는 방법론이 연구되었다. 해석을 위한 형상은 해석에 소요되는 시간을 줄이고, 형상이 해석결과에 미치는 영향을 줄이고자 2차원으로 구성하였다. 해석을 위한 수학적모델으로는 난류모델, 가스화반응모델, Lagrangian particle tracking, Char reaction 등을 포함하였고, 해석을 위한 Solver는 Fluent를 이용하였다. 모델링결과에 의해 예측되는 합성가스의 조성을 상용급 IGCC 가스화기의 운전결과와 비교해 본 결과 본 연구에서 설정한 모델로 예측되는 온도 및 가스농도가 실험치와 유사하게 나타남을 알 수 있었고 이를 통하여 본 연구에서 설정한 모델링방법이 적절함을 알 수 있었다.

• 한국농업기계학회:학술대회논문집
• /
• 한국농업기계학회 2017년도 춘계공동학술대회
• /
• pp.170-170
• /
• 2017

기후변화 따른 스마트팜 돈사 외부 환경의 변화에 대응하고, 사육 환경을 능동적으로 개선하기 위한 연구가 수행 중이다. 돈사 내 열전달 요소 간 상호 역학성 분석을 위해서 고려해야할 사항은 입기구, 보온 등, 열풍기, 단열제, 위치, 방향, 돈사의 연평균 온도, 습도, 연중 일사량, 가축의 열복사 등 상호 복잡하게 연관되어 있는 물리량이다. 돈사 전체 열손실, 자연발생 에너지량, 강제발생 에너지량, 난방용량 등을 고려한 순간 열부하 산정을 위한 여러 방법 중 우선적으로 CFD(Computational Fluid Dynamics)를 이용하였다. 순간 열부하 산정을 위한 해석 도구 선정에 있어서 다양한 유체 및 기체 전산 유체역학 Solver(Fluent, Open-FOAM, Blender)를 고려하였다. 공간 Mech를 수행하기 위한 도구로는 공개 소프트웨어 인 FreeFem++ 3.51-4 (http://www.freefem.org)를 이용하였다. 이 과정에서 일부 기체 (암모니아)의 농도를 난수로 변화시키는 기법을 적용하여 가상적으로 돈사의 환경을 Pseudo 시뮬레이션 하였다. 결과적으로 Fluent에 비하여 OpenFOAM을 이용하여 얻은 열유동의 방향(속도)과 크기 백터가 상대적으로 크게 나타났다. Fluent가 시계열 상에서 혼합 기체 물리량 변화를 무시할 수 있는 안정되고 균일한 환경에 적합하기 때문인 것으로 판단되었다. Blender의 경우 Lattice Boltzmann methods 과 Smoothed-particle hydrodynamics 방법을 이용한 유체/입자 동력학 모델링을 제공함에 있어 시각적 효과를 강조하는 기능에 중점을 두었다. Fluent와 Blender에서 제공하는 해석 연산 모듈의 정확성 검증을 위해선 공간 분해능을 높인 정밀 계측 시스템을 이용하여 검증할 필요가 있다. Open-FOAM를 이용한 열부하 분석 수행이 상대적으로 높은 절대값을 보이는 특성은 열부하 제어 시스템의 Overshoot를 유발할 가능성이 있으므로 이에 대한 해석 모델의 보정이 추가적으로 필요할 것이다. CFD의 한계인 시간 복잡도를 낮추고 상대적으로 높은 시계열 분해능을 확보할 경우 돈사 내 환기시스템에 맞는 소요 환기량 실시간 산정이 가능해지고 외부기상 및 돈사내부 복사열을 활용함과 동시에 돈군 순환에 상응하는 실시간 열부하 관리 시스템 도출이 가능할 것이다.

• 천문학회보
• /
• 제38권2호
• /
• pp.93.1-93.1
• /
• 2013

The Earth's radiation belts consist of an inner belt and an outer belt, being separated by the slot region. It is well known that the variations of the inner edge of the outer belt and the location of the plasmapause (Lpp) are closely related to each other. Different waves exist inside and outside the plasmasphere, playing different roles in the particle dynamics. The plasmapause is well known to be influenced by solar wind conditions and geomagnetic disturbances. Therefore, it is important to precisely determine the location of the plasmapause and develop a prediction scheme. In this study, we identified the location of the plasmapause using the plasma density data from the Time History of Events and Macroscale Interactions During Substorms (THEMIS). The plasmapause is determined by requiring density gradient of a factor of 15 within L-change = 0.5. We statistically determined Lpp as a function of preceding geomagnetic indices. Also, we determined the relations between Lpp and preceding solar wind conditions by estimating correlation coefficients. These relations give us predicting models of Lpp as a function of preceding solar wind parameters and geomagnetic indices. As our database covers a period over the ascending phase from near-sunspot minimum, our statistical results differ somewhat from previous works that cover near-sunspot maximum. Finally, we give some comparative examples obtained from the Van Allen Probes data.

• 한국전산유체공학회지
• /
• 제20권3호
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Computers and Concrete
• /
• 제8권6호
• /
• pp.647-675
• /
• 2011

Our ability to predict hydration behavior is becoming increasingly relevant to the concrete community as modelers begin to link material performance to the dynamics of material properties and chemistry. At early ages, the properties of concrete are changing rapidly due to chemical transformations that affect mechanical, thermal and transport responses of the composite. At later ages, the resulting, nano-, micro-, meso- and macroscopic structure generated by hydration will control the life-cycle performance of the material in the field. Ultimately, creep, shrinkage, chemical and physical durability, and all manner of mechanical response are linked to hydration. As a way to enable the modeling community to better understand hydration, a review of hydration models is presented offering insights into their mathematical origins and relationships one-to-the-other. The quest for a universal model begins in the 1920's and continues to the present, and is marked by a number of critical milestones. Unfortunately, the origins and physical interpretation of many of the most commonly used models have been lost in their overuse and the trail of citations that vaguely lead to the original manuscripts. To help restore some organization, models were sorted into four categories based primarily on their mathematical and theoretical basis: (1) mass continuity-based, (2) nucleation-based, (3) particle ensembles, and (4) complex multi-physical and simulation environments. This review provides a concise catalogue of models and in most cases enough detail to derive their mathematical form. Furthermore, classes of models are unified by linking them to their theoretical origins, thereby making their derivations and physical interpretations more transparent. Models are also used to fit experimental data so that their characteristics and ability to predict hydration calorimetry curves can be compared. A sort of evolutionary tree showing the progression of models is given along with some insights into the nature of future work yet needed to develop the next generation of cement hydration models.