• 제목/요약/키워드: Thermodynamic modeling

검색결과 100건 처리시간 0.024초

수소불화탄소 및 수소염화불화탄소 냉매(R-134a, R-227ea, R-236fa, R-141b)를 이용한 가스 하이드레이트 형성에 관한 계산화학적 해석 (Computational Chemistry Study on Gas Hydrate Formation Using HFC & HCFC Refrigerants (R-134a, R-227ea, R-236fa, R-141b))

  • 김경민;안혜영;임준혁;이제근;원용선
    • Korean Chemical Engineering Research
    • /
    • 제55권5호
    • /
    • pp.704-710
    • /
    • 2017
  • 가스 하이드레이트 형성원리를 이용한 해수담수화는 이미 상용화된 역삼투 방식에 비하여 아직 실증화 단계이지만 그 공정이 비교적 단순하고 특히 냉매를 객체가스로 사용할 경우 아주 낮은 공정 온도가 필요하지 않아 에너지 소비량(thermal budget)이 향상될 가능성이 있기 때문에 여전히 많은 관심을 받고 있다. 따라서 본 연구에서는 수소불화탄소(HFC, hydrofluorocarbon) 및 수소염화불화탄소(HCFC, hydrochlorofluorocarbon) 계열의 냉매들을 객체가스로 한 가스 하이드레이트 형성 거동을 에너지적인 관점에서 해석하고자 하였고 이를 위해 밀도 범함수(DFT, density functional method) 이론을 기반으로 한 분자모델링을 도입하였다. 객체가스(guest gas)로 R-134a, R-227ea, R-236fa, R-141b를 선정하였으며 계산을 위하여 물 분자로 이루어진 $5^{12}$, $5^{12}6^2$, $5^{12}6^4$의 세가지 구조의 동공들(cavities)을 구성하였다. 동공, 객체가스, 그리고 객체가스가 삽입된 동공의 구조를 분자모델링을 이용하여 각각 최적화하였고 계산된 각 구조의 에너지로부터 동공과 객체가스의 결합에너지(binding energy)를 계산하였다. 마지막으로 결합에너지를 비교함으로써 어느 냉매가 가장 유리한 조건에서 가스하이드레이트를 형성할 지를 판단하였다. 결과적으로 R-236fa가 가장 자발적(spontaneous)으로 가스 하이드레이트를 형성할 것으로 예상되었고 사람에 대한 낮은 독성과 물에 대한 작은 용해도 측면에서도 가장 적절한 선택으로 평가되었다.

반응표면분석법을 이용한 폐감귤박 활성탄에 의한 수중의 2,4-Dichlorophenol 흡착특성 해석 (Adsorption Characteristics Analysis of 2,4-Dichlorophenol in Aqueous Solution with Activated Carbon Prepared from Waste Citrus Peel using Response Surface Modeling Approach)

  • 이창한;감상규;이민규
    • Korean Chemical Engineering Research
    • /
    • 제55권5호
    • /
    • pp.723-730
    • /
    • 2017
  • 폐감귤박으로 제조한 활성탄(WCAC)에 의한 2,4-디클로로페놀(2,4-DCP) 흡착에서 온도, 초기농도, 접촉시간 및 흡착제 투여량과 같은 운전변수의 영향을 조사하기 위해 회분식 실험 및 반응표면분석법(Response Surface Methodology: RSM)을 적용하였다. 2,4-DCP 흡착부터 도출된 회귀식은 반응변수의 함수로 나타낼 수 있었다. 이 모델의 적합성은 응답에 대한 실험값과 예측값 간의 상관관계에 의해 평가되었다. $R^2$ 값은 0.9921로서 높은 상관성을 가지며, 회귀 모델에 의해 대부분의 데이터 변동을 설명할 수 있었다. 독립변수 및 그 상호작용의 유의성은 분산분석(ANOVA)과 t-검정 통계 기법으로 평가하였다. 이들 결과는 사용된 모델이 응답변수를 유의미하게 잘 부합되며, 응답과 독립 변수 간의 관계를 적합하게 잘 설명한다는 것을 보여 주었다. 흡착 속도 및 등온 실험결과는 각각 유사 2차 속도식 및 Langmuir 등온 모델에 의해 잘 설명될 수 있었다. Langmuir 등온 모델로부터 계산된 WCAC에 의한 2,4-DCP의 최대 흡착량은 345.49 mg/g이었다. 흡착과정에서 막확산과 입자내부확산이 동시에 일어나는 것을 흡착 메커니즘 연구로부터 확인하였다. 열역학적 파라미터는 WCAC에서 2,4-DCP의 흡착 반응이 흡열반응이고 자발적인 과정임을 나타내었다.

열화학 황-요오드 수소 생산 공정의 요오드 결정화기 설계를 위한 결정 침강 모델링 (Crystal Sinking Modeling for Designing Iodine Crystallizer in Thermochemical Sulfur-Iodine Hydrogen Production Process)

  • 박병흥;정성욱;강정원
    • Korean Chemical Engineering Research
    • /
    • 제52권6호
    • /
    • pp.768-774
    • /
    • 2014
  • 황-요오드(Sulfur-Iodine, SI) 공정은 물을 분해시켜 수소를 생산하는 열화학 공정으로 공정에 사용되는 황과 요오드는 재순환된다. SI 공정 중 요오드가 분리 순환되는 Section III에서는 공정 효율 개선을 위해 다양한 방법이 개발되고 있다. EED(electro-electrodialysis)를 이용한 방법은 추가적인 화합물이 필요하지 않는 공정으로 Section III의 효율을 높일 수 있으나 공정 흐름에 포함된 요오드에 의해 상당한 부하가 걸리게 된다. 이를 해결하기 설계를 위한 기초 자료 제거 공정으로 결정화 방법이 고려되고 있다. 본 연구에서는 요오드 결정화 반응기 설계를 위한 기초 자료 확보를 위해 $I_2$ 포화 $HI_x$ 용액에서 요오드 결정의 침강 속도를 모델링 하였다. $HI_x$ 용액 조성은 열역학 모델인 UVa를 이용하여 결정하였으며 용액 물성은 순수한 물성들과 상관관계식을 활용하여 추산하였다. Multiphysics 전산툴을 이용하여 침강에 따른 속도 변화를 계산하였으며 요오드 직경과 온도에 따른 변화를 추산하였다. 직경(1.0~2.5 mm)과 온도($10{\sim}50^{\circ}C$) 범위에서 요오드는 0.5 m/s 내외의 종말 속도를 보이며 이 속도는 용액의 점도 보다 밀도에 더 크게 영향을 받는 것으로 나타났다.

표면착물 모델을 이용한 합성 침철광과 석영의 카드뮴 흡착 모사 (Surface Complexation Modeling of Cadmium Sorption onto Synthetic Goethite and Quartz)

  • 옥용식;정진호;이옥민;임수길;김정규
    • 한국토양비료학회지
    • /
    • 제36권4호
    • /
    • pp.210-217
    • /
    • 2003
  • 산화물에 대한 카드뮴의 흡착특성을 모사하기 위하여 토양에 존재하는 대표적 산화광물인 침철광과 석영을 대상으로 회분형 카드뮴 흡착실험을 수행하였다. 토양 광물에 대한 카드뮴의 흡착을 예측하기 위해 기존의 경험적 흡착모델 대신 표면착물 모델을 적용하였다. 모델을 검증하기 위하여 카드뮴의 흡착에 영향을 미치는 pH, 중금속 농도, 산화물 농도 및 이온강도 등의 인자를 조절하며, 다양한 조건하에서 회분형 흡착실험을 수행하였다. 실험결과, 카드뮴의 농도가 증가함에 따라 pH에 따른 흡착 곡선은 알칼리 쪽으로 평행하게 이동하였다. 흡착 실험으로부터 계산한 50%의 카드뮴이 흡착되는 pH인 $pH_{50}$은 합성 침철광이 5.25로 석영의 7.83에 비해 매우 낮아 카드뮴에 대한 표면 결합력이 석영에 비해 합성 침철광이 우세함을 시사하였고 이는 모델의 모사 결과와도 잘 부합하였다. 또한 배경전해질로 사용한 $KNO_3$의 이온강도는 카드뮴의 흡착에 영향을 미치지 않았다. 따라서 카드뮴은 산화광물 표면에 내부이온흡착 메커니즘에 의하여 흡착되는 것으로 판단할 수 있었다.

1-Butyl-3-methylpiperidinium Bis(trifluoromethylsulfonyl)imide 이온성 액체에 대한 이산화탄소의 고압 용해도 (High-Pressure Solubility of Carbon Dioxide in 1-Butyl-3-methylpiperidinium Bis(trifluoromethylsulfonyl)imide Ionic Liquid)

  • 남상규;이병철
    • 분석과학
    • /
    • 제27권2호
    • /
    • pp.79-91
    • /
    • 2014
  • 1-Butyl-3-methylpiperidinium bis(trifluoromethylsulfonyl)imide ($[bmpip][Tf_2N]$) 이온성 액체를 대상으로 약 303 K로부터 약 343 K의 온도 범위와 약 30 MPa까지의 압력 범위에서 이온성 액체에 녹는 이산화탄소 ($CO_2$)의 용해도를 측정하였다. 우리가 아는 한, $[bmpip][Tf_2N]$에 대한 $CO_2$의 용해도 데이터는 다른 연구자들에 의해 지금까지 문헌에 발표된 바가 없다. 가변부피투시창 (variable-volume view cell)이 장착된 고압용 상평형 장치를 사용하여 온도를 변화시키면서 여러 가지 조성을 갖는 $CO_2+[bmpip][Tf_2N]$ 혼합물의 기포점 또는 구름점 압력을 측정함으로써 $[bmpip][Tf_2N]$에서의 고압 $CO_2$의 용해도를 결정하였다. 이온성 액체가 가지고 있는 양이온이 $CO_2$ 용해도에 미치는 영향을 관찰하기 위하여, 본 연구에서 사용한 $[bmpip][Tf_2N]$에 대한 $CO_2$ 용해도 데이터를 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide($[bmim][Tf_2N]$)에 대한 $CO_2$ 용해도 데이터와 비교하였다. 압력이 증가함에 따라 $[bmpip][Tf_2N]$에 대한 $CO_2$ 용해도는 급격히 증가하였으며 온도가 증가함에 따라 용해도는 감소하였다. $[bmpip][Tf_2N]$$[bmim][Tf_2N]$ 두 이온성 액체에 대하여 $CO_2$의 용해도는 몰분율 기준으로 온도 및 압력 조건에 관계없이 거의 같았다. Peng-Robinson 상태 방정식을 사용하여 $CO_2+[bmpip][Tf_2N]$ 혼합물 시스템에 대한 상평형 모델링을 수행하였다.

전해질 수용액에서 L-Alanine의 활동도계수와 용해도의 측정 및 모델링 (Measurements and Modeling of the Activity Coefficients and Solubilities of L-alanine in Aqueous Electrolyte Solutions)

  • 이봉섭;김기창
    • Korean Chemical Engineering Research
    • /
    • 제48권4호
    • /
    • pp.519-533
    • /
    • 2010
  • 본 연구에서는 L-형 아미노산인 L-Alanine과 무기염인 NaCl, KCl, $NaNO_3$$KNO_3$의 각 전해질로 이루어진 L-Alanine/전해질 수용액 계에서 L-Alanine의 활동도계수와 용해도를 298.15 K에서 측정하였다. L-Alanine의 활동도계수는 양이온 및 음이온의 선택성 전극으로 이루어진 화학전지에서 두 전극간의 기전력을 측정하는 전기화학 법으로 측정하였으며, 용해도는 L-Alanine의 고체상과 상평형을 이루고 있는 포화용액을 중량 분석하여 측정하였다. 한편 본 연구에서는 아미노산(L-Alanine)/전해질 수용액 계의 잔류(residual) Helmholtz 자유에너지를 섭동사슬-통계역학적 회합성유체이론(perturbed chain-statistical associating fluid theory)과 단순-평균구근사(primitive-mean spherical approximation)이론을 결합한 관계로 모델링 하였으며, 이로부터 아미노산의 활동도계수 및 용해도에 대한 열역학적 관계식을 얻었다. Helmholtz 자유에너지를 모델링 하는 과정에서는 아미노산은 양쪽성 이온(zwitterion) 형태로 존재하며 아미노산의 양쪽성 이온은 같은 이온끼리 자기-회합(self-association)하며 동시에 물분자와 교차-회합(cross-association)하는 회합체로 가정하였으며, 또한 아미노산의 양쪽성 이온이 전해질(무기염)로부터 해리된 양이온 및 음이온과 상호작용하여 이온복합체(ion complex)를 형성하는 과정을 회합현상으로 가정하였다. 본 연구에서 제안된 이론적 모델로부터 L-Alanine/전해질 수용액 계에서 계산되는 L-Alanine의 활동도계수 및 용해도 값은 본 연구의 실험값과 일치하는 경향을 보였다.

촉매를 이용한 이산화탄소 고정화 및 고리형 카보네이트 합성반응에 대한 계산화학적 해석 (Computational Chemistry Study of CO2 Fixation and Cyclic Carbonate Synthesis Using Various Catalysts)

  • 안혜영;김민경;정희철;엄기헌;원용선
    • 청정기술
    • /
    • 제22권1호
    • /
    • pp.35-44
    • /
    • 2016
  • 본 연구에서는 이론화학 또는 분자모델링이라 불리는 계산화학을 기존에 연구된 촉매반응 실험결과를 해석하는 데 접목시켜 보았다. 온실효과의 주범인 이산화탄소를 에폭사이드와 반응을 통해 고정화하고 카보네이트를 생성하는 반응을 선택하였는데 이 반응은 활성화 에너지(55~59 kcal/mol)가 높아 촉매의 사용이 불가피하다. 많은 기존 연구들 중에 ZIF-90/ionic liquid immobilized ZIF-90 (IL-ZIF-90), polystyrene-supported quaternized ammonium salt, KI/KI-glycine, dimethylethanolamine (DMEA)를 촉매로 사용한 경우에 대해 반응의 각 단계를 계산하여 반응의 경로를 예측하고 단계별로 구한 에너지를 바탕으로 에너지도를 구성함으로써 실험결과를 열역학적으로 해석하였다. ZIF-90/IL-ZIF-90과 KI/KI-glycine의 경우는 실험적으로 후자들이 상대적으로 높은 수율을 얻었는데 계산 결과 활성화에너지가 낮아진 이유가 아니라 전자들의 경우 반응 중간체가 높은 에너지를 가져 반응물로 되돌아가는 역반응에 의해 정반응의 진행이 방해를 받은 것으로 밝혀졌다. DMEA를 촉매로 사용하였을 경우는 활성화 에너지를 ~42 kcal/mol로 낮춰줌으로써 금속이나 할로겐염 없이도 촉매의 활성이 잘 일어남을 증명하였다. 폴리스티렌(polystyrene)으로 지지된 quaternized 암모늄염 촉매의 경우 클러스터 가정을 사용하여 계산을 진행하였으며 암모늄염의 -NH기와 에폭사이드의 O 원자 사이의 상호작용의 가능성을 확인할 수 있는 반응경로를 제시하였다.

열-수리-역학적 연계해석을 이용한 복공식 지하 압축공기에너지 저장공동의 역학적 안정성 평가 (Geomechanical Stability of Underground Lined Rock Caverns (LRC) for Compressed Air Energy Storage (CAES) using Coupled Thermal-Hydraulic-Mechanical Analysis)

  • 김형목;;류동우;신중호;송원경
    • 터널과지하공간
    • /
    • 제21권5호
    • /
    • pp.394-405
    • /
    • 2011
  • 본 연구에서는 복공식 지하 압축공기에너지 저장공동의 역학적 변형 및 누출 거동의 복합거동을 파악할 목적으로 비등온 다상다성분 유체유동 및 역학적 거동의 연계해석이 가능한 TOUGH-FLAC 해석을 실시하였다. 지하압축공기에너지 저장 공동의 초기 및 장기 운영 과정에서 고압 압축공기 인입 입출에 따른 콘크리트 라이닝 내부에 발생하는 응력 양상을 살펴보고 저장공동 내부 압력 및 온도 변화를 파악함으로써 기밀성능을 평가하였다. 최대 저장공동 운영압력 8 MPa 조건에서 콘크리트 라이닝 내부에서는 공기침투압에 의한 유효응력의 감소와 접선방향의 인장응력의 증가에 따라 인장균열이 발생할 수 있음을 확인하였다. 콘크리트 라이닝 내부의 인장균열 발생에 따른 투과특성 증가 모델을 이용한 해석 결과, 저장공동 천정부 및 측벽부 일부에서 인장파괴가 발생하여 이들 영역에서의 투과계수는 초기 $10{\times}10^{-20}m^2$에서 $5.0{\times}10^{-13}m^2$까지 증가하였다. 한편, 콘크리트 라이닝 내부 인장균열 발생 및 투과특성 증가에도 불구하고 저장공동 내부 압축공기 압력은 주변 암반의 기밀성능으로 인해 일정하게 유지되고 공기누출량은 일일주입량의 0.02%에도 못 미쳐 복공식 지하 압축공기에너지 저장공동의 유효성을 확인할 수 있었다.

해양-빙붕 상호작용을 고려한 남극 테라노바 만에서 수괴 형성과 순환의 수치 시뮬레이션 (Numerical Simulation of Ocean - Ice Shelf Interaction: Water Mass Circulation in the Terra Nova Bay, Antarctica)

  • 김태균;진경;나지성;이춘기;이원상;문재홍
    • Ocean and Polar Research
    • /
    • 제44권4호
    • /
    • pp.269-285
    • /
    • 2022
  • The interaction between ocean and ice shelf is a critical physical process in relation to water mass transformations and ice shelf melting/freezing at the ocean-ice interface. However, it remains challenging to thoroughly understand the process due to a lack of observational data with respect to ice shelf cavities. This is the first study to simulate the variability and circulation of water mass both overlying the continental shelf and underneath an ice shelf and an ice tongue in the Terra Nova Bay (TNB), East Antarctica. To explore the properties of water mass and circulation patterns in the TNB and the corresponding effects on sub ice shelf basal melting, we explicitly incorporate the dynamic-thermodynamic processes acting on the ice shelf in the Regional Ocean Modeling System. The simulated water mass formation and circulation in the TNB region agree well with previous studies. The model results show that the TNB circulation is dominated by the geostrophic currents driven by lateral density gradients induced by the releasing of brine or freshwater at the polynya of the TNB. Meanwhile, the circulation dynamics in the cavity under the Nansen Ice shelf (NIS) are different from those in the TNB. The gravity-driven bottom current induced by High Salinity Shelf Water (HSSW) formed at the TNB polynya flows towards the grounding line, and the buoyance-driven flow associated with glacial meltwater generated by the HSSW emerges from the cavity along the ice base. Both current systems compose the thermohaline overturning circulation in the NIS cavity. This study estimates the NIS basal melting rate to be 0.98 m/a, which is comparable to the previously observed melt rate. However, the melting rate shows a significant variation in space and time.

메탄올의 이산화탄소 흡수평형 추산에 대한 PC-SAFT모델식과 Two-model approach 모델식의 비교연구 (Comparative Study on the Estimation of CO2 absorption Equilibrium in Methanol using PC-SAFT equation of state and Two-model approach.)

  • 노재현;박회경;김동선;조정호
    • 한국산학기술학회논문지
    • /
    • 제18권10호
    • /
    • pp.136-152
    • /
    • 2017
  • 본 연구에서는 $CO_2$ 제거 용매로써 메탄올 수용액을 사용하는 $Rectisol^{(R)}$공정을 모델링하기 위한 열역학 모델식으로는 PC-SAFT(Perturebed-Chain Statistical Associating Fluid Theory) 상태방정식과 액체활동도계수 모델식을 기본으로 조합된 Two-model approach식{NRTL(Non Random Two Liquid) + Henry + Peng-Robinson}을 비교하였다. 또한 PC-SAFT 상태방정식의 이성분계 상호작용 매개변수와 Two-model approach식의 Henry 상수를 새롭게 결정하기 위해서 273.25K과 262.35K에서 $CO_2$와 메탄올 간의 흡수평형실험을 수행하고 회귀분석을 하였다. 그리고 새롭게 결정한 매개변수의 정확성은 실험 데이터의 추산결과를 통해 검증하였다. 이러한 모델식과 검증한 매개변수를 사용하여 $CO_2$ 제거공정을 모델링 하였다. 그 결과 Two-model approach식을 사용한 경우가 PC-SAFT EOS을 사용한 경우에 비해 $CO_2$ 99.00% 제거하기 위해 요구되는 메탄올 용매 유량이 약 43.72% 더 높게 추산되었으며, 증류탑에서의 냉각수 소모량은 39.22%정도, 스팀소모량은 43.09%정도 더 소요됨을 알 수 있었다. 결론적으로 고압에서 운전되는 $Rectisol^{(R)}$ 공정을 Henry관계식의 도움을 받는 액체활동도계수 모델식을 사용하여 모델링을 하는 경우 PC-SAFT 상태방정식을 사용한 경우에 비해서 크게 설계된 다는 것을 알 수 있었다. 이러한 이유는 액상에 대한 용해도가 낮은 가스성분이 일정한 온도에서 액상에 녹아드는 양은 기상의 분압에 비례하여 증가하는 것으로 계산되는 Henry 관계식의 특성 때문에 메탄올에 대해 용해도가 큰 $CO_2$의 경우 메탄올과 $CO_2$간의 흡수특성을 잘 예측하지 못하는 것을 알 수 있었다.