최적화된 소스 함수 기반 차세대 풍파 모델

초록

본 논문은 기존 연구를 토대로 비선형 상호작용 근사(DIA)를 최적화하고, 대기 동적 경계층을 포함한 입력항, 그리고 파고 스펙트럼에 대해 2차 형태의 소산항을 결합한 새로운 소스 함수를 제시한다. 제안된 함수는 WAM·WAVEWATCH 대비 정확도와 연산 속도에서 우수함을 학술 테스트를 통해 검증했으며, 차세대(4세대) 파랑 예측 모델의 기반으로 활용 가능함을 논한다.

상세 요약

본 연구는 풍파 모델의 핵심인 소스 함수(source function)를 세 부분으로 재구성하였다. 첫 번째는 비선형 전이 메커니즘을 기술하는 이산 상호작용 근사(discrete interaction approximation, DIA)의 최적화이다. 기존 DIA는 파동 스펙트럼의 4차원 적분을 근사화하면서 계산량을 크게 줄였지만, 파라미터 선택에 따라 정확도가 크게 달라지는 단점이 있었다. 저자는 실험적 최적화 절차를 통해 상호작용 삼각형의 각도와 파수 비율을 재조정하고, 스펙트럼의 고주파 영역에서의 과소평가를 보정하는 가중치를 도입하였다. 이 과정에서 비선형 전이 항의 에너지 보존성을 유지하면서도 연산 비용을 기존 DIA 대비 30 % 이상 감소시켰다.

두 번째 요소는 입력항(input term)이다. 전통적인 입력항은 경험적 계수를 이용해 바람 속도와 파고 성장률을 연결했으나, 대기와 해양 사이의 동적 경계층(dynamic atmospheric boundary layer, DABL)을 무시한다는 한계가 있었다. 저자는 대기 경계층의 전단 프로파일과 난류 점성 계수를 파라미터화한 별도 블록을 추가함으로써, 바람의 유속 프로파일이 파면에 미치는 영향을 정량화하였다. 이 블록은 대기 안정도와 파면 거칠기(roughness)를 연계시켜, 동일한 10 m 풍속에서도 실제 파면에 전달되는 동압이 달라지는 현상을 재현한다. 결과적으로 입력항은 기존 경험식 대비 파고 성장률을 5 %~12 % 정도 정확하게 예측하였다.

세 번째는 소산항(dissipation term)이다. 대부분의 3세대 모델은 선형 혹은 1차 형태의 소산을 사용해 파고 감소를 기술했지만, 실제 파괴 과정은 파 스펙트럼의 에너지 밀도와 비선형 상호작용에 의존한다. 저자는 파고 스펙트럼 S(k,θ)와 직접 곱해지는 2차 형태 – D = β S² – 를 제안하였다. 이 항은 고주파 파동이 저주파 파동에 에너지를 전달하면서 발생하는 파괴 메커니즘을 수학적으로 표현한다. 이론적 유도는 파동-파동 상호작용 이론과 파면 파괴 실험 데이터를 기반으로 하며, β 계수는 파면 거칠기와 물리적 점성에 따라 변하도록 설계되었다. 실험 결과, 2차 소산항은 파고 스펙트럼의 고주파 감쇠를 기존 모델보다 15 %~20 % 더 정확히 재현하였다.

모델 검증을 위해 저자는 WAM·WAVEWATCH에서 사용되는 표준 학술 테스트(단일풍, 주기적 바람, 변동 바람 등)를 수행하였다. 최적화된 DIA와 DABL 입력항, 2차 소산항을 결합한 새로운 모델은 평균 제곱근 오차(RMSE)와 상관계수(R)에서 기존 모델 대비 각각 0.12와 0.06 정도 개선되었으며, 연산 시간은 0.78배로 단축되었다. 특히, 급격한 바람 변화 상황에서 파고와 주파수 스펙트럼의 전이 과정을 보다 현실적으로 포착함을 확인하였다.

이러한 결과는 소스 함수의 물리적 정밀도와 계산 효율성을 동시에 달성할 수 있음을 보여준다. 특히 대기 경계층 블록의 도입은 바람-파면 상호작용을 통합적으로 기술함으로써, 향후 기후 모델링이나 해양‑대기 연계 예보 시스템에 직접 적용 가능한 기반을 제공한다. 저자는 이러한 특성을 바탕으로 차세대(4세대) 풍파 모델은 “물리‑수치‑관측 통합 프레임워크”를 목표로 해야 한다고 제언한다.