플라스틱 멀치 적용 토양 온도 전이 해석: 비선형 유한요소 모델과 최적 설계

초록

**
본 논문은 플라스틱 멀치를 덮은 토양의 1차원 비선형 열전달을 유한요소법(Galerkin)으로 모델링하고, 복사·대류 비선형 경계조건을 포함한 전이 방정식을 풀어 실험 데이터와 비교 검증하였다. 검증된 모델을 이용해 멀치의 광학적 투과·반사 특성을 Taguchi 방법으로 최적화하여, 현재 콜롬비아에서 널리 사용되는 LDPE 멀치보다 높은 토양 온도를 얻을 수 있는 최적 조합(τs > 0.7, ρs < 0.2, τl < 0.01, ρl > 0.7)을 제시한다. 또한 연간 기후 데이터를 적용해 깊이별 온도 변화를 시뮬레이션함으로써 1 m 이하에서는 일일 변동이 거의 사라짐을 확인하였다.

**

상세 분석

**
이 연구는 플라스틱 멀치가 토양 표면에 미치는 열효과를 정량화하기 위해 1‑차원 비선형 열전달 방정식을 설정하고, Galerkin 유한요소법을 이용해 공간을 50개의 2차 라그랑주 요소로 이산화하였다. 비선형성은 주로 멀치와 대기·토양 사이의 복사 교환에서 발생하는 T⁴ 항에 기인하며, 이를 고정점 반복법으로 해결한다. 경계조건은 (1) 멀치와 토양 사이의 대류·복사 열플럭스, (2) 멀치와 대기 사이의 복사·대류 열플럭스, (3) 토양 하부는 무한반체 가정으로 온도 구배가 0인 Neumann 조건을 적용한다. 복사 계수는 스윈뱅크식과 실험적으로 도출된 대류계수를 사용해 구체화하였다.

모델 검증을 위해 콜롬비아 보고타 사반나 지역에서 LDPE 멀치를 설치하고, 현장 기상 데이터(태양복사, 기온, 풍속)와 토양 표면 온도를 5분 간격으로 측정하였다. 광학 특성(τs = 0.733, ρs = 0.265, τl = 0.6, ρl = 0.398)은 FTIR 및 UV‑VIS 분광법으로 확보했으며, 토양 열전도도와 비열은 전기 저항 가열 실험으로 추정하였다. 시뮬레이션 결과는 평균 상대오차 5.76 %로 실험과 높은 일치를 보였으며, 급격한 구름 이동에 의한 순간적인 복사 변동은 측정 평균값의 한계로 설명된다.

Taguchi 직교표를 활용한 다변량 최적화에서는 τs, ρs, τl, ρl 네 가지 광학 파라미터를 0.01, 0.2, 0.7 등 3단계로 변동시켜 최대 토양 온도(Tmax)를 평가하였다. 결과는 τs > 0.7(고투과), ρs < 0.2(저반사), τl < 0.01(극저투과), ρl > 0.7(고반사) 조합이 최적임을 보여준다. 이는 현재 일반적으로 사용되는 LDPE(τl ≈ 0.6, ρl ≈ 0.4)보다 훨씬 높은 온도 상승을 기대할 수 있음을 의미한다.

연간 시뮬레이션에서는 깊이 1 m 이하에서 일일 온도 변동이 거의 사라지고, 전체 연간 변동은 2.5 m 이하에서 소멸한다는 결과를 얻어, 토양 온도 안정성을 위한 멀치 적용 깊이 기준을 제시한다.

본 연구는 비선형 FEM을 통한 토양‑멀치 열전달 모델링, 실험 검증, 그리고 설계 최적화를 일관되게 수행한 점에서 학술적·실용적 기여가 크다. 다만, 1‑차원 가정과 멀치의 무열용량 가정이 실제 다층 구조와 복합 재료 특성을 완전히 반영하지 못한다는 한계가 남는다. 향후 다차원 모델링과 멀치 재료의 동적 열용량 고려가 필요하다.

**