τ만 바꿨는데 정지한 기체의 밀도가 바뀌었다 — Shan–Chen 속도 이동 forcing의 (τ−½)² 응력
정지한 D1Q3 기체 기둥에 같은 힘을 걸면 Shan–Chen forcing만 τ에 따라 밀도가 달라진다. 평형 속도 이동이 (τ−½)²F²/ρ 응력을 남긴다.
v 항과 v² 항을 더했더니 양 끝은 맞고 가운데서 42% 빗나갔다 — Prony 두 항 저항법칙과 Re=60
Bossut은 저항이 v²에, Du Buat는 v에 비례한다고 했다. Prony의 합 a·v²+b·v는 Re=60에서 주도권이 바뀌고, 그 근처에서 종단속도를 42%까지 틀린다.
1 μm 알루미늄 입자가 연속체 공식보다 4.6배 늦게 데워졌다 — Knudsen 전이영역의 Sherman 조화 보간
1500 K, 1 atm에서 1 μm 알루미늄의 Kn은 0.67이다. 연속체 Nu=2는 가열 시간을 4.58배 짧게 본다.
σ₁₂에 1.00을 넣고 0.78을 돌려받았다 — 3상 Allen–Cahn에서 없는 상이 계면으로 새는 자리
3상 Allen–Cahn에서 Σc=1을 평균 투영으로 맞추면 없는 상이 계면으로 새어 표면장력이 21.9% 깎인다. 3중 벌점항이 그 자리를 막는다.
이동도를 10배 키웠는데 계면 두께는 그대로였다 — 보존형 Allen–Cahn의 반확산 항
반확산 항이 확산을 정확히 상쇄해 계면 두께를 이동도와 무관하게 고정하는 원리, 그리고 그 항이 LBM 평형분포로 들어가는 자리
항력계수를 9% 깎았더니 액적이 35% 더 느려졌다 — 변형과 기화 보정이 붙는 자리
구형·TAB 변형·Renksizbulut–Yuen 기화 보정 세 항력계수 모델을 같은 액적에 걸어 관통거리로 비교한다.
아무 일도 일어나면 안 되는 문제에서 압력이 21% 튀었다 — Abgrall 기준과 비보존항
균일 압력·속도 유동이 균일하게 유지되지 않는 원인은 플럭스가 아니라 비보존항 이산화다. Re–Abgrall 논문의 H_u 연산자를 같은 격자에서 직접 재봤다.
Δt를 5배로 키우니 1.9배 빨라졌고, 10배에서는 제자리였다 — 암시적 표면장력이 여는 창
모세관 제약을 깨도 상한은 남는다. Janodet(2025)이 밀도비 1000에서 잰 안정 한계와 속도 이득의 최적점을 따라간다.
격자를 절반으로 줄였더니 발산이 두 배 빨리 왔다 — 6방정식 two-fluid 모델의 복소 고유값
단일압력 6방정식 two-fluid 모델은 미끄럼 속도가 있으면 고유값이 복소수가 된다. 계면압력 항의 계수가 1을 넘는 순간 실수로 돌아온다.
[논문 리뷰] 이류와 압력을 완전히 갈라놓는다 — Baer–Nunziato flux의 TV 분리
BN flux를 이류와 압력으로 가르면 음속은 P-system에만 남고, Noh 정체 영역에서 압력이 음수로 떨어진다.
[논문 리뷰] 물방울은 음파를 거의 되돌려 보낸다 — 음향 임피던스와 sharp 계면 2상 솔버
공기-물 반사 에너지는 99.889%다. 기포 1%만 섞여도 Wood 음속은 1500에서 120.5 m/s로 떨어진다.
빨래는 왜 5도에서도 마르는가 — 확산 지배 증발과 열 지배 비등
5도 빨래는 포화 질량분율이 0이 아니라 마르고, 100도 비등은 열수지와 Lee의 r이 그 속도를 정한다.
번진 계면을 매 스텝 되돌리다 — 안티확산 계면 샤프닝(IST) 직접 구현
200스텝 뒤 헬륨 계면이 6칸에서 20칸으로 번지면, IST 후처리는 ε/a만으로 그 두께를 상수로 잠근다.
정지한 물방울이 스스로 흐르는 이유 — 기생 전류와 well-balanced 계면 장력
정지한 물방울에서 곡률을 3%만 틀리면 잔차 최대가 4.8e5로 뛰고, 맞으면 4.1e-10에 남는다. 에너지는 24자리가 벌어진다.
빗방울은 왜 눈물방울이 아닌가 — Weber 수와 Taylor 유비 분열(TAB)
지름 기준 Weber가 12를 넘으면 TAB 변형 y가 1을 찍어 방울이 갈라진다. 오버슈트가 임계를 낮춘다.
물방울 시뮬레이션이 Δt에 발목 잡히는 이유 — 표면장력을 음함수로 푼 논문
N=64에서 256으로 가면 Δt_σ는 8분의 1이 된다. ωΔt_σ는 늘 2.221이고, 명시적은 0.9배에서 터진다.
두 상에 각자의 방정식을 준다 — Baer–Nunziato 7-방정식 모델
인쇄된 324 bar는 강성 기체 두 큰 항이 1.2%만 남기고 상쇄된 결과다. 완화의 첫 스텝이 그 상쇄를 민다.
[논문리뷰] 단 한 스테이지로 2차 정확도 — Kapila 다상류를 위한 GRP 솔버
N=320에서 Godunov 오차는 2.31e-2, 단일 스테이지 GRP는 1.55e-5로 차수 2.03이다.
물방울은 왜 둥글고 컵 가장자리의 물은 왜 차오르는가 — 표면장력과 접촉각
코드 출력에서 5 µm 방울의 압력 점프는 28800 Pa, 0.1 mm 관은 30도에서 127.1 mm다.
물 위에서 외친 소리는 왜 물속에 닿지 않는가 — 분산 계면 음향 투과
공기–물 임피던스 비 약 3600에서 반사율이 99.9%에 붙는다. 임계각 13.4도와 약압축성 보간을 1D 관으로 재현한다.
[논문 리뷰] Poisson 없이 두 상이 산다 — 일반 압력 방정식과 MSTACS의 완전 명시적 결합
인공음속 c_s를 5에서 15로 키우면 Δt는 1/3이 되고, 부피오차 E_G는 6.57×10^{-8}이다.
SPH 입자법으로 유체 풀기 — 커널부터 인공 점성까지
밀도 합으로 압력을 만드는 Tait 식(γ=7)은 1%의 밀도 오차를 크게 증폭한다. 대칭 항을 빼면 그 압력이 운동량으로 샌다.
압축성 다상유동을 위한 Rhie-Chow 보간법 총정리
물-공기 면에서 밀도비 1000에서 산술 평균은 약 500.5, 조화 평균은 2.0이다. 필터는 ∇P에만 건다.
5-Equation Model 구현 가이드: Python으로 만드는 다상유동 솔버
α=0.5에서 혼합 음속은 약 24 m/s다. 같은 초기 상태의 공기는 약 374, 물은 약 1625 m/s다.
계면 포착 기법 비교: VOF vs Level Set vs Diffuse Interface
압축성 계면에서 α를 보존형으로 풀면 압력이 흔들린다. 5방정식은 α의 이류만 비보존형으로 남겨 평형을 지킨다.
압축성 다상유동 CFD 입문: 왜 어렵고, 왜 중요한가
계면에서 γ를 평균하면 압력이 흔들린다. α=0.5에서 혼합물의 음속은 단상값의 100분의 1까지 떨어진다.