Taylor–Green 와류와 점성 감쇠

에너지는 줄어도 와류의 모양은 남을 수 있다

부드럽게 회전하는 유체 셀이 격자처럼 반복된다고 상상해봅시다. 시간이 지나면 운동은 약해지지만 무늬는 알아볼 수 있게 남습니다. 이때 방정식의 어느 부분이 속력을 줄이고 어느 부분이 운동 방향을 꺾을까요? 점성이 모든 것을 맞춘다는 설명으로는 부족합니다. 이 예에서는 압력이 비선형 대류를 맞추고, 점성이 속도의 시간 감쇠를 맞춥니다.

이차원 Taylor–Green 가족은 명시한 비압축성 방정식의 매끄러운 정확 해입니다. 수치 솔버와 관찰할 장을 혼동하지 않고 속도, 압력, 미분, 물리적 수지와 입자 경로를 비교할 수 있습니다. 학습 모델의 일관성을 확인하는 유용한 예이지만 일반적인 삼차원 난류를 재현하는 모델은 아닙니다.

선수 지식, 가정과 단위

속도 기울기, 발산, 물질 가속도, 응력과 운동 에너지를 복습합니다. 밀도는 양의 상수이고 동점성계수는 음이 아닌 상수입니다. 주기 정사각형에는 벽이나 체적력이 없습니다. 명시한 깊이를 곱해 이차원 셀을 실제 체적으로 바꿉니다. 비압축성 모델에서는 압력의 상수 기준을 선택할 수 있으며 운동량에는 그 기울기가 들어갑니다.

기호	단위	의미
rho	kg/m³	일정한 질량 밀도
U	m/s	초기 속도 진폭
k	rad/m	양의 공간 파수
nu	m²/s	동점성계수
깊이	m	평면 밖 실제 두께

$$\partial_t\mathbf u+(\mathbf u\cdot\nabla)\mathbf u=-\frac{\nabla p}{\rho}+\nu\Delta\mathbf u,\qquad\nabla\cdot\mathbf u=0,\qquad\mu=\rho\nu.$$

실제 속도장 구성하기

진폭은 지수적으로 감소시키고 두 속도 성분에는 회전 대칭의 공간 구조를 줍니다. 부호가 중요합니다. 아래 관례에서 원점 부근의 와도는 양수입니다. U는 균일한 배경 병진 속도가 아니라 초기 무늬의 진폭입니다. 셀의 평균 속도는 영이고 개별 입자는 공간적으로 다른 속도장을 지나갑니다. U를 평균 유속으로 읽으면 수송과 에너지의 의미를 모두 틀리게 됩니다.

$$A(t)=Ue^{-2\nu k^2t},\qquad u=-A(t)\cos(kx)\sin(ky),\qquad v=A(t)\sin(kx)\cos(ky).$$

실제 정사각형 한 변은2pi/k이고 장은 경계를 넘어 반복됩니다. 주기적 동일시이지 미끄럼 없는 고체벽 조건이 아닙니다. 그림의 정사각형 테두리를 모든 변에 벽이 있다는 뜻으로 읽지 않습니다. 깊이를 포함한 체적과 질량도 따로 표시합니다.

$$L=\frac{2\pi}{k},\qquad V=L^2d,\qquad M=\rho V.$$

비압축성과 국소 회전 확인하기

두 속도 성분을 각자의 좌표로 미분합니다. 발산 기여가 셀 평균에서만이 아니라 모든 점에서 상쇄됩니다. 비압축성은 체적 변화율을 제한하는 조건이지 속도 기울기의 모든 성분이 영이라는 뜻이 아닙니다. 기울기의 반대칭 부분은 와도와 연결되고 대칭 부분은 변형률 속도를 나타냅니다.

$$u_x=Ak\sin(kx)\sin(ky),\quad v_y=-Ak\sin(kx)\sin(ky),\quad\omega_z=v_x-u_y=2Ak\cos(kx)\cos(ky).$$

x가 오른쪽이고 y가 위쪽인 축에서 원점의 양의 와도는 반시계 방향 국소 회전과 맞습니다. 화면의 세로축을 뒤집으면서 해석을 그대로 두면 회전 방향이 거꾸로 보일 수 있습니다. 와도만으로 변형이나 에너지 소산을 모두 설명하지 않습니다.

비선형 대류를 먼저 계산하기

국소 속도를 따라 각 속도 성분을 방향 미분합니다. 유동 무늬가 단지 약해지더라도 대류 항은 일반적으로 영이 아닙니다. 입자가 무늬 속에서 방향을 바꾸면서 가속하기 때문입니다. 따라서 시간에 따라 변하지 않는 Euler 관점의 장에서도 물질 가속도는 비영일 수 있습니다. 시각화에서는 시간 미분과 대류 가속도를 분리해 읽습니다.

$$(\mathbf u\cdot\nabla)\mathbf u=-\frac{A^2k}{2}\begin{pmatrix}\sin(2kx)\\\sin(2ky)\end{pmatrix}.$$

사인 내부의 두 배 파수는 속도와 속도 기울기를 곱해서 생깁니다. 압력을 속도와 같은 단일 공간 주파수로만 표현하는 학습 모델은 속도 그림을 그럴듯하게 만들더라도 압력 응답을 놓칠 수 있습니다.

압력이 운동 방향을 꺾는 가속도를 제공한다

압력 변화량은 속도보다 공간 주파수가 두 배이고 지수 감쇠율도 두 배입니다. 부호를 비교하기 전에 직접 미분합니다. 압력 가속도는 압력 기울기를 밀도로 나눈 값의 음수입니다. 대류를 양의 압력 기울기와 비교하면 힘의 부호를 뒤집습니다. 두 좌표가 모두 기여하므로 일차원 압력 곡선 하나만으로 전체 장을 설명할 수 없습니다.

$$p-p_0=-\frac{\rho A^2}{4}\big[\cos(2kx)+\cos(2ky)\big],\qquad-\frac{\nabla p}{\rho}=-\frac{A^2k}{2}\begin{pmatrix}\sin(2kx)\\\sin(2ky)\end{pmatrix}.$$

p-zero를 바꾸면 상수만 더해지고 운동은 변하지 않습니다. 기준에 대한 음의 압력 변화량을 음의 열역학적 절대 압력이라고 해석하면 안 됩니다. 이 비압축성 모델은 완전한 상태 방정식이나 공동현상 모델을 제공하지 않습니다.

점성이 시간 감쇠를 제공한다

각 속도 성분에는 두 좌표의 사인 또는 코사인 인자가 하나씩 있습니다. 따라서 라플라시안은 그 성분에 음의2k²를 곱합니다. 지수 진폭을 시간 미분하면 여기에 동점성계수가 곱해진 같은 벡터장이 나옵니다. 앞의 압력·대류 균형과 합치면 운동량 방정식이 점별로 성립합니다. 단지 전체 평균 잔차가 영인 것이 아닙니다.

$$\partial_t\mathbf u=-2\nu k^2\mathbf u,\qquad\Delta\mathbf u=-2k^2\mathbf u,\qquad\partial_t\mathbf u=\nu\Delta\mathbf u.$$

점성이 영이면 정상 비선형 Euler 와류가 됩니다. 압력과 대류는 일반적으로 여전히 비영입니다. 점성을 껐다고 입자가 정지하거나 압력 힘이 없어지는 것은 아닙니다.

한 점에서 하는 손계산

밀도와 파수를1, 초기 진폭을2, 동점성계수를0.1로 놓습니다. 시간 영, x와 y가 각각 pi/4인 점에서 속도는(−1,1)m/s입니다. 대류는 음의 x와 음의 y 쪽을 함께 가리킵니다. 시간 가속도는 속도의 반대 방향이며 이 점에서는 더 작습니다. 렌더링된 애니메이션에 기대지 않고 화살표와 라벨을 확인할 수 있는 수치 예입니다.

$$\rho=k=1,\ U=2,\ \nu=0.1,\ t=0,\ x=y=\pi/4:\quad\mathbf u=(-1,1),\quad(\mathbf u\cdot\nabla)\mathbf u=(-2,-2),\quad\mathbf u_t=(0.2,-0.2).$$

속도 단위는m/s이고 두 가속도 단위는m/s²입니다. 서로 다른 물리량의 화살표 길이를 비교하려면 각각의 축척을 밝혀야 합니다. 설명 없는 공통 벡터 축척 하나로 물리적 크기까지 비교할 수는 없습니다.

실제 주기 셀의 에너지 적분

운동 에너지 밀도는 밀도와 속력 제곱의 곱의 절반입니다. 두 삼각함수 기여를 깊이까지 포함한 한 주기 셀에서 적분합니다. 사인 제곱과 코사인 제곱은 각각 평균이 절반입니다. 따라서 평균 에너지 밀도의 계수는 절반이 아니라 사분의 일입니다. 압력과 대류는 국소적으로 에너지를 옮기지만 완전한 주기 셀의 경계를 통한 순수송은 상쇄됩니다.

$$\overline e_K=\frac{\rho A^2}{4},\qquad K=\frac{\rho U^2V}{4}e^{-4\nu k^2t},\qquad[K]=\mathrm J.$$

임의의 작은 관측 창에는 이 상쇄가 없습니다. 그 창의 에너지 변화를 모두 점성 소산이라고 해석하기 전에 경계 수송을 포함해야 합니다. 체적과 에너지 밀도를 곱한 실제 J와 면적당 값을 분리합니다.

소산은 회전만이 아니라 변형률 속도에서 나온다

속도 기울기의 대칭 부분 D를 만듭니다. 뉴턴 비압축성 유체에서 국소 운동 에너지의 내부 에너지 전환율은 동적 점도와 D의 텐서 노름 제곱을 곱한 값의 두 배입니다. 여기서 점도는 동점성계수 nu가 아니라 밀도를 곱한 mu입니다. 텐서 노름에는 두 비대각 성분을 모두 세어야 합니다. 이 특수한 장에서는 비대각 합이 영이고 두 대각 성분이 일반적으로 남습니다. 와도는 국소 회전의 관측이며 소산식의 국소 변형률 텐서를 그대로 대체할 수 없습니다.

$$D=\frac{\nabla\mathbf u+(\nabla\mathbf u)^T}{2},\qquad\Phi=\int_V2\rho\nu(D:D)\,dV=4\nu k^2K,\qquad\frac{dK}{dt}=-\Phi.$$

Phi의 단위는W입니다. 운동 에너지를 줄이는 것이지 전체 에너지를 파괴하는 것은 아닙니다. 이 교육 모델은 온도를 진화시키지 않으므로 전환된 내부 에너지의 공간 분포를 직접 계산하지는 않습니다.

유선, 입자 경로와 수치 오차 구분하기

유선은 한 시각의 속도장에 접하는 곡선이고 입자 경로는 시간 의존 초기값 문제의 해입니다. 이 가족은 고정된 공간 무늬에 시간 스칼라만 곱하므로 입자가 같은 기하학적 유선을 따라가며 느려집니다. 곡선의 기하가 일치하는 것은 특수한 성질입니다. 시간에 따른 매개화까지 같다는 뜻도 아니고 일반적인 비정상 흐름에 그대로 적용할 수도 없습니다.

$$\dot X=u(X,Y,t),\quad\dot Y=v(X,Y,t),\qquad\frac{d}{dt}\big[\cos(kX)\cos(kY)\big]=0.$$

RK4 추적기는 실제 초와 경계에서 감지 않은 좌표를 사용해 이 방정식을 근사합니다. 시간 간격을 절반으로 줄이고 독립적으로 계산한 끝점과 비교합니다. 기하학적 라벨의 보존도 확인하지만, 잘못된 시계로 움직이는 추적기조차 같은 라벨을 보존할 수 있습니다. 끝점 수렴과 시간 의존 벡터장의 사용은 별도로 확인해야 합니다.

압력을 제거하는 대조 실험

같은 실제 속도장을 두고 운동량 평가에서 압력 가속도만 제거합니다. 시간 항과 점성 항은 이미 맞으므로 남은 잔차는 대류 항입니다. 속도장의 발산이 영이라는 이유로 이 잔차가 없어지지는 않습니다. 다른 방정식을 만족하는 별개 흐름으로 바꾸지 않고, 빠진 물리 기여 하나를 드러내는 대조입니다.

$$R=\mathbf u_t+(\mathbf u\cdot\nabla)\mathbf u+\frac{\nabla p}{\rho}-\nu\Delta\mathbf u,\qquad R_{\mathrm{omit}\ p}=(\mathbf u\cdot\nabla)\mathbf u.$$

각 벡터를 독립적으로 정규화한 그림에서는 방향은 읽을 수 있어도 크기 비교는 불가능합니다. 방향 전용 도식 옆에는 실제 SI 성분 값을 함께 남깁니다.

작은 변화와 반올림된 같음

아주 작은 양의 점성이나 시간에서는 부동소수점 감쇠 인자가1로 반올림됩니다. 여기서1을 빼면 수학적으로는 비영인 변화가 영으로 나옵니다. expm1로 변화를 직접 계산하고 두 관측을 구분해 남깁니다. 에너지의 지수는 속도 감쇠 지수의 두 배입니다. 시간 영이나 점성 영의 정확한 영 변화와, 표현 범위 아래의 비영 곱을 구분합니다.

$$z=2\nu k^2t,\qquad\frac{A(t)-U}{U}=\operatorname{expm1}(-z),\qquad\frac{K(t)-K(0)}{K(0)}=\operatorname{expm1}(-2z).$$

이 비율의 해석은 U가 비영일 때 적용됩니다. 구현에는 절대 변화량도 남기므로 정지 상태의 변화는 정확히 영으로 읽을 수 있습니다. 양의 지수가 표시 범위 아래이면 상대 변화는 미정이어도 가중한 절대 변화는 표현 가능할 수 있습니다. 각 관측을 따로 계산해야 합니다. 속도가 표현 가능해도 진폭의 제곱은 언더플로할 수 있습니다. 작은 수치 잔차 하나만으로 이런 척도의 에너지와 운동까지 정확하다고 증명할 수는 없습니다.

에너지 밀도와 셀 전체의 에너지를 별도 관측으로 생각해 보세요. 작은 비영 밀도는 선형 표시 범위 아래여도 실제 셀 체적을 곱한 에너지는 표현 가능할 수 있습니다. 구현은 체적과 진폭 제곱을 먼저 결합하고 마지막에 반올림합니다. 작은 감쇠도 초기 진폭·에너지와 결합한 뒤 절대 변화를 표시합니다. binary64의 가장 작은 척도에서는 지수 변화의 이차 보정이 반올림 중점의 어느 쪽인지 결정하기도 합니다. 작은 변화를 항상 일차 항으로 바꾸면 마지막 비트를 잘못 고를 수 있습니다. 좌표와 삼각함수 인자의 일반적인 유한 정밀도는 유지됩니다.

실험 순서와 AI 연구 활용

시간 영의 속도, 압력, 와도를 먼저 봅니다. 대류와 압력 가속도를 따로 표시한 다음 시간 미분과 점성 가속도를 비교합니다. 점성을 영으로 놓아 정상 Euler 극한을 관찰합니다. 점성을 복원하고 속도 감쇠와 더 빠른 에너지 감쇠를 비교합니다. 압력 기준을 바꾸어 힘과 운동이 그대로인지 확인합니다. 잔차 대조에서만 압력 항을 제거합니다. 마지막으로 여러 시간 간격의 입자 끝점과 작은 점성의 안정적 변화·반올림 관측을 비교합니다.

학습 PDE 모델에서는 보지 않은 점성과 파수를 평가용으로 남깁니다. 발산, 단위를 가진 운동량 잔차, 압력 기준 처리, 적분 에너지와 실제 입자 운동을 따로 보고합니다. 진폭을 무조건 작게 만드는 모델도 낮은 잔차를 얻을 수 있으며, 보기 좋은 유선만으로 입자의 시간 진행을 검증할 수는 없습니다. 이 정확 해 가족은 통제된 벤치마크이며 벽, 충격파, 임의 초기 자료나 삼차원 난류 에너지 전달을 다룬다는 증거는 아닙니다.

해설이 있는 연습

연습 1: 실제 셀

파수가2rad/m이고 깊이가0.5m일 때 한 주기 셀의 변 길이와 체적을 구하세요.

해설: 변 길이는 pi m이고 체적은0.5pi²m³입니다. 밀도를 곱해 kg을 얻거나 에너지 밀도를 곱해 J를 얻기 전에 깊이를 포함해야 합니다.

연습 2: 점성이 영인 경우

점성을 영으로 놓으면 운동량 방정식의 모든 항이 영인가요?

해설: 아닙니다. 시간 항과 점성 항은 사라지지만 비선형 대류와 압력 가속도는 일반적으로 남아 서로 같습니다. Euler 관점의 장은 정상이지만 입자는 곡선을 따라 돌며 가속합니다.

연습 3: 두 감쇠율

속도 진폭이 처음의 절반이 되었습니다. 셀 운동 에너지는 처음의 얼마인가요?

해설: 사분의 일입니다. 밀도와 셀 기하가 고정된 상태에서 에너지는 속도 진폭의 제곱에 비례합니다. 소산율은m/s의 속도가 아니라J/s의 에너지 변화율과 비교합니다.

연습 4: 압력 기준

속도장을 유지하고 모든 점의 압력에3Pa를 더하면 운동량의 어느 항이 달라지나요?

해설: 어느 항도 달라지지 않습니다. 압력 기울기는 그대로입니다. 표시 압력만 이동하고 이 비압축성 모델의 압력 가속도, 궤적, 운동 에너지와 소산은 같습니다.

연습 5: 유선의 일치

두 수치 추적기가 올바른 기하학적 곡선에 머물지만 같은 시간에 다른 위치에 도달합니다. 둘 다 검증되었나요?

해설: 아닙니다. 기하학적 라벨의 보존은 곡선만 확인합니다. 시간 의존 ODE, 독립 끝점과 시간 간격 세분화를 비교해야 합니다. 잘못된 시계로 올바른 유선을 움직일 수도 있습니다.

연습 6: 반올림 감쇠 대조

양의 점성에서 표시된 감쇠 인자는1인데 expm1은 음의 변화를 줍니다. 어느 관측을 지워야 하나요?

해설: 둘 다 남깁니다. 반올림된 인자와 더 안정적으로 계산한 변화라는 정의가 다릅니다. 비영 변화 자체가 산술 범위 밖이면 해당 관측을 계산 불가로 남기고 정확히 점성이 영인 상황으로 바꾸어 해석하지 않습니다.

자료와 적용 범위

NVIDIA Warp의 Taylor–Green 예제는 단위 정사각형의 정확 속도·압력을 구현한1차 자료입니다. 여기서는 차원 척도, 실제 수지와 대조 실험을 별도로 유도했습니다. 주기적인 교육 셀은 명시적인 모델 선택이며 그 예제의 모든 수치 경계 조건이나 솔버 선택을 이어받는 것은 아닙니다. FEM 솔버 코드나 그 솔버 결과를 이 구현의 검증 근거로 가져오지 않습니다.