강직 시스템을 위한 내재 이차 안정성을 갖는 새로운 일반 선형 방법

초록

본 논문은 차수 p = q이며 단계 수 r = s = p + 1인 새로운 일반 선형 방법(GLM)을 제안한다. 내재 이차 안정성(IQS) 조건과 A‑안정성·L‑안정성 제약을 동시에 만족하도록 계수 행렬을 설계하고, 오류 상수를 최소화한다. 1차부터 4차까지의 암시적 GLM을 구성하고, Van der Pol 진동자, Burgers 방정식, Gray‑Scott 모델 등 세 가지 강직 문제에 적용해 기존 방법과 비교한다. 실험 결과는 제안된 스키마가 정확도와 안정성 면에서 경쟁력을 갖음을 보여준다.

상세 분석

이 연구는 일반 선형 방법(GLM)의 설계 자유도를 활용해, 특히 r = s = p + 1이라는 Nordsieck 입력 벡터 구조를 채택함으로써 단계와 외부 출력 차원을 동일하게 맞춘다. 이러한 구성은 가변 단계 크기와 가변 차수 구현에 유리하며, 기존의 s = p 구조에서 발생할 수 있는 큰 오류 상수 문제를 완화한다. 논문은 먼저 GLM의 기본 형태(2.1)를 제시하고, A, U, B, V 행렬을 각각 s × s, s × r, r × s, r × r 형태로 정의한다. 여기서 A와 V는 (2.4)와 같이 삼각형 형태를 갖으며, λ > 0인 암시적 파라미터와 V의 첫 행이 단위 행벡터인 특성을 이용해 안정성을 보장한다.

차수 조건은 Theorem 2.2에 따라 e c z = z A e c z + U W 등으로 정리되며, 이는 Nordsieck 벡터가 차수 p = r − 1을 만족하도록 하는 일련의 선형 방정식이다. 이후 내재 이차 안정성(IQS) 조건을 Definition 2.3의 행렬 등가 관계 BA ≡ XB, BU ≡ XV − VX 로 도입한다. X 행렬은 (2.16)과 같이 상위 2행만 비제로인 형태를 가정함으로써, M(ω) = V + ω B(I − ω A)⁻¹U 의 특성 다항식이 ηʳ⁻²·