페클레 수

편집 역사 토론
유체역학 수치해석 마지막 수정: 2026-07-12 04:09:14

1. 개요[편집]

페클레 수
Péclet Number
기호Pe
정의대류 수송 / 확산 수송
분야유체역학 · 열전달 해석 · 수치해석
관계식Pe = Re × Pr (열) = Re × Sc (물질)

강물에 잉크를 떨어뜨렸다. 잉크가 퍼져 나가는가, 흘러가 버리는가 — 그 답이 페클레 수다.

페클레 수(Péclet number, Pe)는 유동장에서 물리량(열 또는 물질)이 대류(convection)로 실려 가는 속도확산(diffusion)으로 퍼지는 속도의 비를 나타내는 무차원수다. 프랑스 물리학자 장 클로드 외젠 페클레의 이름을 땄으며, 정의는 다음과 같다.

Pe=ULαPe = \frac{U L}{\alpha}

여기서 UU는 대표 속도, LL은 대표 길이, α\alpha는 확산계수(열이면 열확산율, 물질이면 물질확산계수)다. Pe가 크면 대류가 지배해 물리량이 흐름을 타고 실려 가고, Pe가 작으면 확산이 지배해 물리량이 사방으로 고르게 퍼진다. Pe → 0의 극한은 순수 확산(라플라스 방정식류), Pe → ∞의 극한은 순수 대류(쌍곡형 이류방정식)에 해당한다.

2. 레이놀즈수·프란틀수와의 관계[편집]

페클레 수는 독립적으로 태어난 숫자가 아니라 두 무차원수의 곱으로 분해된다. 열 문제에서는

Pe=RePr=ULνναPe = Re \cdot Pr = \frac{UL}{\nu}\cdot\frac{\nu}{\alpha}

레이놀즈수(관성/점성)와 프란틀수(운동량확산/열확산)의 곱이다. 물질 전달 문제에서는 프란틀수 자리에 슈미트수(Sc)가 들어가 Pe = Re·Sc가 된다.1 이 관계는 왜 물과 공기가 열을 다르게 실어 나르는지를 설명한다. 물은 Pr ≈ 7, 공기는 Pr ≈ 0.7이라, 같은 레이놀즈수라도 물의 열 페클레 수가 열 배 크다. 열전달 해석에서 유동장(Re)과 열 수송(Pe)을 따로 봐야 하는 이유가 여기에 있다.

3. 셀 페클레 수: 수치해석의 지뢰밭[편집]

CFD와 열전달 해석에서 페클레 수가 진짜로 무서운 얼굴을 드러내는 곳은 격자 셀 단위다. 전체 유동의 Pe가 아니라, 격자 간격 Δx\Delta x를 대표 길이로 쓴 셀 페클레 수를 정의한다.

Pecell=UΔxαPe_{\text{cell}} = \frac{U \, \Delta x}{\alpha}

이 숫자가 차분 도식의 운명을 가른다. 1차원 정상 대류-확산 방정식을 중심차분(유한체적법)으로 이산화하면, 셀 페클레 수가 임계값을 넘는 순간 해가 물리적으로 말이 안 되는 진동(oscillation)을 일으킨다. 그 임계값이 바로

Pecell=UΔxα>2Pe_{\text{cell}} = \frac{U\,\Delta x}{\alpha} > 2

이다.2 중심차분은 상류·하류 정보를 대칭으로 취급하는데, 대류가 강해지면(셀 Pe > 2) 흐름은 명백히 상류 편향인데도 하류 값을 절반이나 반영하다 보니 계수 행렬이 대각지배성을 잃고 해가 톱니처럼 요동친다. 온도장에 물리적으로 존재할 수 없는 언더슛/오버슛이 나타나는 것.

4. 진동을 잡는 처방들[편집]

셀 페클레 수 문제의 해법은 결국 세 갈래다.

  • 격자 세분화Δx\Delta x를 줄여 셀 Pe를 2 아래로 내린다. 가장 정직하지만 격자 수가 폭증한다. 적응 격자 세분화가 이 부담을 덜어 준다.
  • 상류 도식(upwind) — 대류항을 흐름의 상류 쪽 값으로만 계산한다. 셀 Pe와 무관하게 진동이 사라지지만, 대신 수치확산(numerical diffusion)이라는 가짜 확산을 도입해 해를 뭉갠다. 즉 안정성을 정확도와 맞바꾼다.3
  • 고차·혼합 도식 — QUICK, TVD, 지수 도식(exponential scheme) 등은 상류의 안정성과 중심차분의 정확도 사이에서 절충한다. 특히 지수 도식은 1차원 대류-확산의 정확해를 셀 페클레 수의 함수로 그대로 담고 있어, 파탄카의 하이브리드/멱법칙(power-law) 도식의 이론적 뿌리가 된다.

실무 요령은 명확하다. 대류가 지배하는 영역(큰 셀 Pe)에서는 중심차분을 고집하지 말고, 격자를 조이거나 상류 계열 도식으로 갈아타라. “중심차분으로 돌렸더니 온도장에 이빨이 났다”면 십중팔구 셀 페클레 수가 범인이다.4

5. 어디에 쓰이나[편집]

페클레 수는 CFD를 넘어 폭넓게 등장한다. 히트싱크의 강제대류 냉각 설계, 원자로 냉각재 유동, 잉크젯 프린팅의 액적 확산, 미세유체 칩 안 화학종의 혼합 효율, 심지어 해양의 오염물 확산 예측까지 — 대류와 확산이 경쟁하는 모든 곳에 Pe가 있다. 특히 스토크스 유동 규모의 미세유체에서는 레이놀즈수가 작아도 페클레 수는 여전히 클 수 있어(확산이 워낙 느려서), “층류인데 안 섞인다”는 마이크로믹서 설계의 근본 난제로 이어진다.

6. 관련 문서[편집]

7. Footnotes[편집]

  1. 슈미트수 Sc = ν/D는 운동량확산과 물질확산의 비다. 열의 프란틀수와 완전히 대칭적인 구조라, 열-물질 유사성(heat-mass analogy)의 근거가 된다. 자연은 게을러서 같은 수학을 재활용한다.

  2. 정확히는, 중심차분 계수가 음수가 되지 않을 조건이 |Pe_cell| ≤ 2다. 이 값을 넘으면 이웃 셀 계수가 음수가 되어 최대원리(maximum principle)가 깨진다. 물리적으로 “온도가 두 이웃보다 더 뜨겁거나 차가운” 셀이 생긴다는 뜻.

  3. 수치확산의 크기는 대략 α_num ≈ U·Δx/2에 비례한다. 즉 상류 도식은 셀 페클레 수만큼의 가짜 확산을 공짜로 끼워 넣는 셈. 격자가 흐름에 비스듬할수록 이 오차는 더 커진다(교차확산).

  4. 반대로 확산이 지배하는 영역(작은 셀 Pe)에서 굳이 상류 도식을 쓰면 정확도만 손해다. 도식 선택은 로컬 셀 페클레 수를 보고 정하는 게 맞고, 그래서 하이브리드 도식이 존재한다.