콜브룩 방정식 단 한 번 로그 호출로 푸는 파데 근사법

초록

본 논문은 콜브룩 방정식의 마찰계수를 구할 때, 최초 한 번만 로그 함수를 호출하고 이후 반복에서는 파데(Padé) 다항식으로 대체하는 새로운 반복 알고리즘을 제시한다. 정확도와 수렴 횟수는 기존 뉴턴‑라프슨 방식과 동일하지만, 연산 비용이 크게 감소한다.

상세 분석

콜브룩 방정식은 로그와 비정수 거듭제곱을 포함해 직접 해를 구할 수 없는 암묵식이다. 전통적으로는 뉴턴‑라프슨이나 고정점 반복을 사용하지만, 매 반복마다 로그 연산이 필요해 CPU 사이클을 많이 소모한다. 저자들은 로그 함수를 파데 근사식으로 대체함으로써 연산량을 기본 사칙연산으로 축소한다. 파데 근사는 분자와 분모가 다항식인 유리함수 형태이며, 근사 차수가 높을수록 정확도가 향상된다. 논문에서는 (1,1)부터 (3,2)까지 다양한 차수의 파데식을 제시하고, 특히 (2,3) 차수가 0.61.6 구간에서 상대오차가 10⁻⁴% 이하임을 실험적으로 확인했다. 초기 추정값은 기존과 마찬가지로 도메인 내 임의값을 사용해도 수렴하지만, 고정값 혹은 간단한 다항식(예: 유전 알고리즘으로 도출된 40% 오차 식)으로 시작하면 평균 45회 반복으로 충분히 수렴한다. 중요한 점은 첫 번째 반복에서만 실제 로그를 계산하고, 이후에는 파데식으로 대체하므로 전체 로그 호출 횟수가 1회로 제한된다는 것이다. 이는 특히 대규모 배관망 시뮬레이션에서 메모리와 연산시간을 크게 절감한다. 또한 Lambert W 형태 변환 시 발생할 수 있는 오버플로우 문제를 회피할 수 있다. 실험 결과는 기존 Newton‑Raphson과 동일한 수렴 속도와 10⁻⁶ 이하의 상대오차를 보이며, CPU 벤치마크에서는 로그 호출을 제외한 경우 20~30% 정도의 속도 향상을 기록한다. 따라서 파데 기반 단일 로그 반복법은 정확도 손실 없이 계산 효율성을 높이는 실용적인 대안이다.