다중척도법

수치해석 시뮬레이션 계산물리 마지막 수정: 2026-09-02 04:27:41

1. 개요[편집]

원자 하나하나로 교량을 계산할 수는 없다. 그렇다고 연속체 방정식이 균열 끝의 원자를 알 리도 없다. 그래서 둘을 붙인다.

다중척도법(multiscale method)은 하나의 문제 안에 크기가 몇 자릿수씩 차이 나는 시공간 척도가 공존할 때, 각 척도를 그에 맞는 모델로 따로 풀고 그 사이에 정보를 흘려 전체 답을 구성하는 방법론의 총칭이다. 미시 모델은 정확하지만 관심 영역 전체를 덮기엔 터무니없이 비싸고, 거시 모델은 값싸지만 미시가 만들어 내는 물성·구성식을 스스로 알지 못한다. 다중척도법은 이 둘 사이의 거래를 체계화한 것이다.

다중격자법과는 다른 물건이다. 이름이 비슷하고 둘 다 여러 해상도의 격자를 쓰지만, 다중격자법은 하나의 이산 방정식을 빨리 푸는 선형 해법 가속기이고 성긴 격자는 어디까지나 오차의 저주파 성분을 걷어내는 도구일 뿐이라 수렴한 답은 성긴 격자를 안 썼을 때와 정확히 같다. 다중척도법에서 성긴 척도와 미세 척도는 애초에 서로 다른 물리 모델이고, 성긴 쪽 답은 미세 쪽 정보 없이는 존재하지도 않는다. 우연히 알고리즘 구조가 닮은 경우는 있어도 목적이 다르다.¹

2. 척도 분리라는 전제[편집]

거의 모든 다중척도법은 작은 매개변수 $\varepsilon = \ell/L$ (미시 특성길이 / 거시 특성길이)이 존재하고 그것이 충분히 작다는 가정 위에 서 있다. 그러면 미지수를 빠른 변수와 느린 변수의 함수로 동시에 놓고 두 변수를 형식적으로 독립 취급할 수 있다.

이 기계 자체는 특이 섭동 계열의 고전 도구다. 시간축에 적용하면 빠른 진동 $T_0=t$ 와 느린 포락선 $T_1=\varepsilon t$ 를 분리해 세속항 제거 조건에서 진폭 방정식을 뽑는 다중척도 섭동전개가 되고, 그 계산 절차와 예제는 이중 평균법·더핑 방정식 문서에 이미 정리돼 있으니 여기서 반복하지 않는다. 이 문서는 같은 발상을 공간과 모델 계층에 적용하는 쪽, 즉 계산재료·계산역학에서 말하는 다중척도 모델링을 다룬다.

주의할 것은 척도 분리가 가정이지 사실이 아니라는 점이다. 난류처럼 스케일이 연속 스펙트럼으로 이어져 있으면 $\varepsilon$ 을 고를 자리가 없고, 그때는 분리 대신 필터링을 쓰게 된다(대와류 모사). 척도 사이의 관계를 계통적으로 다루는 가장 일반적인 문법은 재규격화군이고, 아래 기법들은 전부 그 특수한 실용판으로 읽을 수 있다.

3. 점근 균질화와 유효 물성[편집]

주기 미세구조를 가진 매질에서 $a(x/\varepsilon)$ 처럼 급격히 변하는 계수를 가진 방정식

-\nabla\cdot\left(a\!\left(\tfrac{x}{\varepsilon}\right)\nabla u^\varepsilon\right) = f

를 생각하자. 해를 두 척도 전개

u^\varepsilon(x) = u_0(x,y) + \varepsilon u_1(x,y) + \varepsilon^2 u_2(x,y)+\cdots, \qquad y = \frac{x}{\varepsilon}

로 놓고 $\nabla \to \nabla_x + \varepsilon^{-1}\nabla_y$ 를 대입한 뒤 차수별로 정리하면, $O(\varepsilon^{-1})$ 차수에서 단위 셀 문제(cell problem)가 떨어져 나온다. 그 셀 문제의 해로부터 계산되는 상수 텐서 $a^{*}$ 가 거시 유효 물성이고, $u_0$ 는 진동이 없는 매끄러운 방정식 $-\nabla\cdot(a^{*}\nabla u_0)=f$ 를 만족한다.

가장 유명한 결과 하나만 짚자. 1차원 층상 매질에서 유효 계수는 산술평균이 아니라 조화평균 $a^{*} = \langle 1/a\rangle^{-1}$ 이다. 전기회로의 직렬 저항과 정확히 같은 이유이며, “평균 물성을 넣으면 되겠지”라는 순진한 직관이 왜 틀리는지를 보여 주는 교과서 예제다. 고체역학에서 이 두 극단이 각각 등변형률 가정의 보이트 상한(산술평균)과 등응력 가정의 로이스 하한(조화평균)이며, 등방 이상 복합재에 대해 훨씬 조인 경계가 해신-슈트리크만 한계다.² 관련 재료 쪽 논의는 복합재 해석과 유효매질이론에 있다.

수치적으로는 단위 셀 문제를 유한요소법으로 풀되 셀 경계에 주기경계조건을 건다. 미세구조가 주기적이지 않으면 대표 체적 요소(RVE)를 잘라 쓰는데, 여기서 두 가지가 늘 말썽이다. 첫째, RVE가 충분히 커야 통계적으로 대표성이 생기고 그 크기는 미리 알 수 없어 수렴 연구를 돌려야 한다. 둘째, 유한한 RVE에서는 경계조건 선택이 답을 바꾼다 — 변위 규정 조건은 강성을 과대평가하고 응력 규정 조건은 과소평가하며, 주기 조건이 그 사이에서 가장 빨리 수렴한다. 재료 데이터시트에 물성 하나만 적혀 있는 것을 볼 때 이 사정을 기억해 두면 좋다.

4. 이질 다중척도법 (HMM)[편집]

균질화는 유효 계수를 미리 한 번 계산해 표로 만들어 두는 접근이다. 그런데 미시 거동이 거시 상태(변형률, 온도, 이력)에 비선형으로 의존하면 표를 만들 수가 없다. 그때 쓰는 것이 이질 다중척도법(heterogeneous multiscale method, HMM)이다.³

구조는 단순하다. 거시 솔버를 정상적으로 돌리되, 매 스텝 각 적분점에서 구성식이 필요한 순간에만 미시 솔버를 짧게 켠다.

거시 상태(변형률, 유동 구배 등)를 미시 상자의 구속조건으로 내려보낸다.
미시 모델(분자동역학, 격자 모형, 미세 FEM)을 작은 상자에서 짧은 시간만 돌린다.
필요한 거시량(응력, 플럭스)만 시공간 평균으로 추정해 올려 보낸다.
거시 스텝을 크게 밟는다.

핵심 이득은 미시 솔버가 거시 도메인 전체를 덮지 않는다는 것이다. 비용은 (적분점 수) × (미시 상자 부피 × 미시 시간)에 비례하고, 척도 분리가 성립하면 미시 상자는 $\ell$ 몇 배면 충분하다. 고체 쪽에서 각 적분점마다 RVE를 하나씩 매달아 도는 $\mathrm{FE}^2$ 방식이 사실상 같은 골격이며, 구성방정식을 손으로 못 쓰겠는 재료에 대한 표준 우회로다.

대가도 분명하다. 미시 솔버는 유한한 표본이므로 올라오는 응력에 통계 잡음이 섞이고, 그 잡음이 거시 뉴턴 반복의 수렴을 망친다. 잡음을 줄이려면 미시 시간을 늘려야 하는데 그게 곧 비용이라, HMM 구현의 절반은 “얼마나 대충 재도 되는가”의 오차 예산 문제다. 요즘은 미시 응답을 회귀 모델로 학습해 캐시하는 대리모델 결합이 흔하다.

5. 원자-연속체 결합과 준연속체법[편집]

균열 선단, 전위 핵, 접촉면처럼 국소적으로만 연속체 가정이 깨지는 문제가 있다. 도메인 전체를 원자로 채우는 것은 낭비고, 전부 연속체로 두면 정작 알고 싶은 것을 놓친다. 그래서 문제 영역을 나눈다.

준연속체법(quasicontinuum, QC)이 이 계열의 원형이다. 원자 전부를 자유도로 두는 대신 대표 원자(repatom)만 남기고 나머지는 유한요소 보간으로 위치를 결정한다. 변형이 완만한 영역에서는 셀 하나가 균일 변형한다고 보고 원자간 퍼텐셜로부터 에너지 밀도를 바로 계산하며(코시-보른 규칙), 결함 근처에서는 보간을 풀고 원자 하나하나를 살린다. 결과적으로 결함 주변만 원자 해상도, 나머지는 연속체 해상도인 적응 격자가 만들어진다.

여기서 이 분야의 가장 유명한 함정이 등장한다. 국소(연속체) 영역과 비국소(원자) 영역이 만나는 핸드셰이크 영역에서, 에너지를 두 방식으로 나눠 더한 뒤 미분하면 균일 변형 상태에서조차 힘이 0이 되지 않는다. 이것이 유령힘(ghost force)이다. 물리적 원인은 아무것도 없고 순전히 에너지 분할이 만든 인공물인데, 크기가 무시할 수준이 아니라 계면 근처에 가짜 응력장을 깔아 놓아 전위가 계면에서 걸리거나 튕기는 사고를 낸다. 대응책은 대체로 셋이다.

유령힘 보정. 균일 변형에서의 잔류력을 미리 계산해 상수 하중으로 빼 준다. 값싸지만 임시방편이고 에너지 보존이 깨진다.
준비국소 정식화. 계면 근처 원자에 국소/비국소가 섞인 특별한 에너지를 부여해 유령힘이 구조적으로 0이 되게 만든다. 무유령힘 결합의 정공법.
힘 기반 결합. 아예 에너지 범함수를 포기하고 힘만 블렌딩한다. 유령힘은 사라지지만 계가 비보존적이 되어 야코비안이 비대칭이 되고, 안정성 해석이 까다로워진다.

같은 문제는 원자-연속체뿐 아니라 밀도범함수이론 영역과 고전 힘장을 붙이는 QM/MM에서도 똑같이 나타난다. 경계에서 잘린 화학결합을 수소로 막고(link atom) 정전 상호작용을 어떻게 넘길지 정하는 문제가 QM/MM판 핸드셰이크다.

6. 조립과 시간 척도[편집]

공간만 문제가 아니다. 분자동역학의 시간 스텝은 결합 진동 주기에 묶여 펨토초 단위인데, 단백질 접힘이나 확산은 마이크로초~초에서 일어난다. 열 자릿수 이상의 간극이라 계산기를 더 사서 메울 수 있는 종류가 아니다. 두 갈래 대응이 있다.

조립(coarse-graining). 원자 여러 개를 비드 하나로 묶고, 원자 수준 분포를 재현하도록 비드 사이의 유효 퍼텐셜을 역설계한다(볼츠만 역전, 힘 정합, 상대엔트로피 최소화). 자유도가 줄어 계산이 싸지는 것은 물론이고, 없어진 자유도가 퍼텐셜 우물을 얕게 만들어 궤적이 실제 시간보다 빨리 흐른다. 이 가속은 공짜가 아니라 정보를 버린 대가다 — 지운 자유도가 하던 마찰과 요동이 사라졌으므로, 동역학을 제대로 되살리려면 랑주뱅 동역학처럼 마찰·잡음 항을 다시 넣어야 한다. 이 구조를 정확하게 쓴 것이 모리-츠방치히 형식론의 기억항이고, 같은 논리가 축소차수모델의 폐쇄항으로 그대로 반복된다.

희귀사건 가속. 상태 사이 전이가 드물다는 사실 자체를 이용한다. 안정 상태 안에서의 진동은 어차피 재미없으니 전이만 세는 동역학 몬테카를로(kMC)로 넘어가거나, 우물에 바이어스 퍼텐셜을 부어 탈출을 빠르게 하고 시간을 되돌려 보정하는 가속 MD 계열을 쓴다. 전이율 자체는 전이상태 이론이 준다. 결국 “미시 모델에서 속도상수를 뽑아 상위 모델에 넣는다”는 점에서 이것도 다중척도 결합이다.

7. 순차 결합과 동시 결합[편집]

지금까지 나온 것들을 두 부류로 정리할 수 있다.

	순차(hierarchical)	동시(concurrent)
정보 흐름	미시 → 거시, 오프라인 한 번	매 스텝 양방향
산출물	물성·구성식·속도상수	결합된 하나의 계
예	균질화, 조립 퍼텐셜, kMC 속도표	HMM, QC, QM/MM, $\mathrm{FE}^2$
비용	싸다(미리 계산)	비싸다(런타임 결합)
한계	거시 상태 의존성을 못 담음	계면 인공물·잡음·안정성

실무 판단 기준은 간단하다. 미시 응답이 거시 상태의 함수로 표로 만들어질 수 있으면 순차로 간다. 표가 안 만들어질 때(경로 의존, 강한 비선형, 국소 결함) 비로소 동시 결합의 값을 치른다. 처음부터 동시 결합으로 뛰어드는 것은 대개 과잉이고, 계면 인공물 디버깅에 인생을 쓰게 된다.

8. 현업에서의 현실[편집]

논문에는 “다중척도”라고 쓰여 있지만 실제로는 미시 시뮬레이션 결과에서 계수 두 개를 뽑아 거시 모델에 넣은 것인 경우가 압도적으로 많다. 그게 나쁜 것은 아니다. 대개 그걸로 충분하다.⁴
척도 분리가 실제로 성립하는지 검증하지 않고 기계를 돌리는 것이 가장 흔한 사고다. $\varepsilon$ 이 $0.3$ 쯤 되는 문제에 균질화를 쓰면 답이 그럴듯하게 틀린다.
계면 근처 결과는 일단 의심하는 것이 안전하다. 유령힘이든 반사되는 탄성파든, 다중척도 결합은 계면에서 대가를 치른다.
검증은 언제나 “미시 모델로 전부 푼 작은 문제”와의 비교다. 그 기준 계산이 불가능한 크기의 문제에 다중척도법을 쓰고 있다면, 검증되지 않은 결과를 들고 있는 셈이라는 것을 인정하고 시작해야 한다.

9. 관련 문서[편집]

헷갈리는 게 당연한 이유가 있다. 대수적 다중격자에서 성긴 공간을 만드는 보간 연산자를 “유효 방정식”으로 읽으면 정말 균질화처럼 보이고, 실제로 다중척도 유한요소법(MsFEM)은 셀 문제의 해를 기저함수로 쓴다는 점에서 두 세계의 중간에 있다. 그래도 선형 해법 가속기와 모델 축약은 다른 목적을 가진다는 구분은 지켜 두는 편이 안전하다. ↩
보이트-로이스 폭은 대비가 큰 복합재에서 웃음이 나올 만큼 넓다. 강성비 100배짜리 조합이면 상하한이 한 자릿수 차이가 나서 “이 구간 어딘가에 있습니다”라는 정보가 사실상 무정보에 가깝다. 그래서 해신-슈트리크만이 논문 한 편으로 유명해진 것이다. ↩
E Weinan과 Björn Engquist가 2003년에 이름을 붙였지만, 발상 자체는 훨씬 오래됐다. “필요한 값만 미시 계산으로 조달한다”는 구조는 이미 여러 분야에서 각자 재발명돼 있었고, HMM의 기여는 그것들을 하나의 추상 틀로 묶고 오차 해석을 붙인 데 있다. 좋은 프레임워크의 표본 같은 사례. ↩
“우리 코드는 원자부터 부품까지 이어 붙였습니다”라는 발표를 들으면 반드시 물어봐야 할 질문이 있다. 척도 사이를 몇 번 왕복했는지, 그리고 계면에서의 오차를 무엇으로 쟀는지. 답이 “한 번 내려보냈습니다”면 그건 순차 결합이고, 그렇게 부르면 아무 문제도 없다. ↩