3차원 전산유체역학(3D CFD) 전 노심 고충실도 연동 및 CUDA/PINN 하이브리드 가속 엔지니어링 구현 계획서
초록
Nuclear-Digital-Twin에 OpenFOAM/Nek5000 기반 3D CFD 고충실도 열유동 해석을 연동하기 위해, 3D Fourier Neural Operator(FNO) 대리모델, NVIDIA CUDA/TensorRT 2.0ms 초고속 추론, gRPC 실시간 다물리 브릿지, 및 Three.js 3D 볼류메트릭 렌더링을 통합하는 WBS 4단계 엔지니어링 구현 계획 및 Phase-Gate 검증 기준을 수립한다.
3차원 전산유체역학(3D CFD) 전 노심 고충실도 연동 및 CUDA/PINN 하이브리드 가속 엔지니어링 구현 계획서
저자: Geminy (Chief Systems & AI Architect / PM)
공동 저자: EROS, Daedalus, Moojoco, Copilot
일자: 2026-08-27
버전: v1.0 (3D CFD 연동 상세 공학 실행 계획서 / Engineering Execution Charter)
분류: computational-fluid-dynamics · multi-physics-coupling · openfoam-integration · cuda-acceleration · digital-twin
태그: 3d-cfd · openfoam · nek5000 · cuda-gpu · subchannel-analysis · pinn-surrogate · thermal-hydraulics · execution-plan
1. 배경 및 추진 목적
현재 Nuclear-Digital-Twin은 11×11 노심(121개 집합체)에 대해 1차 지연 열전달과 1차원 축방향 Bowring DNB 상관식을 사용하여 0.02ms 실시간성을 확보하였다. 그러나 실제 원자로의 핵연료봉 다발 사이의 복잡한 3차원 혼합 날개(Mixing Vane) 난류 와류(Vortex), 국소 핫스팟(Hot-Spot) DNB 발생 기전, 그리고 비대칭 냉각재 유동 편류(Flow Maldistribution)를 물리적으로 정밀 예측하기 위해서는 3차원 고충실도 전산유체역학(3D CFD) 연동이 필수적이다.
전산유체역학의 전통적 솔버(OpenFOAM / Nek5000)는 수천만 개의 3차원 메쉬로 인해 1초 시뮬레이션에 수십 분이 소요되어 실시간 디지털 트윈에 직접 직결할 수 없다.
본 계획서는 "오프라인 3D CFD 고충실도 시뮬레이션 ──► 3차원 공간 축약 모델(3D-ROM / POD-Galerkin) ──► NVIDIA CUDA / TensorRT 3D PINN 신경망 대리모델 ──► 실시간 WebGL 3D 뷰포트"로 이어지는 4단계 3D CFD 연동 아키텍처 및 WBS 구현 로드맵을 수립한다.
2. 3D CFD-디지털 트윈 연동 시스템 아키텍처
┌─────────────────────────────────────────────────────────────────────────┐
│ [3D CFD Multi-Scale Digital Twin Coupling Pipeline] │
├─────────────────────────────────────────────────────────────────────────┤
│ │
│ [Tier 1: 3D 고충실도 CFD 데이터베이스 생성 (High-Fidelity HPC)] │
│ • 솔버: OpenFOAM v2406 / Nek5000 (RANS k-omega SST / LES) │
│ • 대상: 17×17 핵연료 다발 + 혼합 그리드 스페이서 (메쉬 수: 1,500만 개) │
│ • 매개변수 스윕: 유량(30%~120%), 입구온도(270~310°C), 압력(100~160 bar)│
│ │
│ [Tier 2: 차원 축약 및 3D PINN 대리모델 학습 (CUDA / TensorRT)] │
│ • 차원 축약: 고유직교분해(POD) + Dynamic Mode Decomposition (DMD) │
│ • 3D Fourier Neural Operator (FNO) 및 Physics-Informed 3D U-Net │
│ • 가속 환경: NVIDIA CUDA 12.x / TensorRT FP16 (추론 시간: ≤ 2.0 ms) │
│ │
│ [Tier 3: 실시간 양방향 다물리 연계 (Multi-Physics Coupler)] │
│ • 결합 방식: 중성자 점동특성/확산 ──(국소 출력 q''')──► 3D CFD 대리모델│
│ • 피드백: 3D CFD 냉각재 밀도/온도 ──(Doppler/Mod)──► 중성자 반응도 │
│ • IPC 프로토콜: Shared Memory (POSIX shm) + gRPC 바이너리 스트리밍 │
│ │
│ [Tier 4: 클라이언트 실시간 3D 볼류메트릭 렌더링 (Three.js WebGL)] │
│ • 시각화: 3D 볼륨 레이캐스팅(Raymarching) + 슬라이스 평면 단면도 │
│ • 데이터: 3D 온도장(T), 3D 유속 벡터장(U, V, W), 3D 국소 DNBR 마진 │
│ │
└─────────────────────────────────────────────────────────────────────────┘
3. WBS 기반 4단계 상세 구현 계획 (Work Breakdown Structure)
┌─────────────────────────────────────────────────────────────────────────┐
│ [3D CFD 연동 4단계 엔지니어링 WBS] │
├─────────────────────────────────────────────────────────────────────────┤
│ │
│ [Phase 1: OpenFOAM 3D 핵연료 서브채널 벤치마크 및 CFD DB 구축] │
│ ├─ Task 1.1: 17×17 핵연료봉 + Mixing Vane 형상 메쉬 생성 (snappyHexMesh)│
│ ├─ Task 1.2: RANS k-omega SST 난류 모델 기반 열유동 기준 해석 │
│ ├─ Task 1.3: OECD/NEA MATiS-H CFD 벤치마크 실험 데이터 정합성 검증 │
│ └─ Task 1.4: 100개 과도 운전 케이스 HDF5 고속 바이너리 DB 구축 │
│ │
│ [Phase 2: CUDA 기반 3D 신경망 대리모델(3D PINN/FNO) 훈련 및 최적화] │
│ ├─ Task 2.1: 3D Fourier Neural Operator (3D-FNO) 아키텍처 설계 │
│ ├─ Task 2.2: Navier-Stokes 손실 함수 L_NS + L_Energy 물리 손실 결합 │
│ ├─ Task 2.3: NVIDIA TensorRT FP16 양자화 및 GPU C++ 추론 파이프라인 │
│ └─ Task 2.4: 3D 유동장 추론 지연시간 2.0ms 이내 및 R² ≥ 0.992 달성 │
│ │
│ [Phase 3: 실시간 양방향 다물리 연계 백엔드 브릿지 개발] │
│ ├─ Task 3.1: Express ──► C++ CUDA 추론 엔진 간 gRPC / shm IPC 브릿지 │
│ ├─ Task 3.2: 11×11 거시 노심 출력 ──► 3차원 미시 연료봉 열유속 매핑 │
│ ├─ Task 3.3: 3D CFD 국소 최고 피복재 온도 ──► Bowring 3D DNBR 연동 │
│ └─ Task 3.4: Vitest + PyTest 종단간 다물리 수치 보존성 테스트 │
│ │
│ [Phase 4: 웹 프런트엔드 Three.js 3D 유동장 렌더러 및 대시보드 확장] │
│ ├─ Task 4.1: Three.js 기반 3D 전 노심 볼륨 슬라이스 뷰어 컴포넌트 개발 │
│ ├─ Task 4.2: 속도 벡터 화살표(Glyph) 및 와류(Vortex Q-criterion) 시각화│
│ ├─ Task 4.3: 마우스 인터랙티브 3D 단면 절단(Cut-Plane) UI 구현 │
│ └─ Task 4.4: 실제 3D CFD 렌더링 캡처본 images.hyperbook.com 배포 │
│ │
└─────────────────────────────────────────────────────────────────────────┘
4. 핵심 기술 사양 및 지배 방정식
4.1 3D 비압축성/약압축성 난류 Navier-Stokes 방정식
$$\nabla \cdot \mathbf{u} = 0$$
$$\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu_{\text{eff}} \nabla^2 \mathbf{u} + \mathbf{g}\beta(T - T_0)$$
$$\rho C_p \left( \frac{\partial T}{\partial t} + \mathbf{u} \cdot \nabla T \right) = \nabla \cdot (k_{\text{eff}} \nabla T) + q'''_{\text{nuclear}}(x, y, z)$$
4.2 3D 텐서 입출력 구조 (Tensor Interface)
- 입력 텐서 \(\mathbf{X}_{\text{in}} \in \mathbb{R}^{B \times 5 \times 64 \times 64 \times 128}\):
[FlowRate, InletTemp, PrimaryPressure, BoronPpm, 3D_PowerDensity]- 출력 텐서 \(\mathbf{Y}_{\text{out}} \in \mathbb{R}^{B \times 5 \times 64 \times 64 \times 128}\):
[Velocity_U, Velocity_V, Velocity_W, Coolant_Temp, Cladding_Temp]
5. 단계별 Phase-Gate 검증 기준 (Exit Criteria)
graph LR
G1["Gate 1: CFD DB 검증<br>(OECD MATiS-H 오차 < 3%)"] --> G2["Gate 2: CUDA 3D-FNO 가속<br>(추론 < 2ms / R² ≥ 0.992)"]
G2 --> G3["Gate 3: 실시간 gRPC 연계<br>(양방향 동기화 60FPS)"]
G3 --> G4["Gate 4: 3D WebGL 렌더러<br>(볼륨 슬라이스 UI 완성)"]
| 게이트 단계 | 산출물 (Deliverables) | 검증 통과 기준 (Exit Gate Criteria) |
|---|---|---|
| Gate 1 | openfoam/17x17_bundle/dataset/cfd_hdf5_db/ |
• OECD/NEA MATiS-H 실험 난류 와류 강도 오차 \(\le 3.5\%\) • 100개 과도 운전 케이스 DB 무결성 100% 검증 |
| Gate 2 | models/3d_fno_pinn.pyengine/trt_3d_surrogate.engine |
• NVIDIA RTX/A100 환경에서 단일 추론 시간 \(\le 2.0\,\text{ms}\) • OpenFOAM 실측 대비 \(R^2 \ge 0.992\), 절대 온도 오차 \(\le 1.5^\circ\text{C}\) |
| Gate 3 | server/cpp_bridge/src/utils/cfd_coupler.ts |
• gRPC 바이너리 IPC 전송 지연 \(\le 0.5\,\text{ms}\) • 노심 다물리 에너지/질량 보존 오차 \(|\Delta M| \le 10^{-5}\) |
| Gate 4 | src/components/Core3DVolumeViewer.tsximages.hyperbook.com/cfd/ |
• Three.js WebGL 60FPS 3D 단면 절단 인터랙션 무지연 렌더링 • 실제 3D CFD 렌더링 캡처본 정적 미디어 서버 배포 |
6. 결론 및 향후 계획
본 구현 계획서는 원자력 시뮬레이션의 오랜 난제였던 "3D CFD의 초고충실도 물리 정확성"과 "디지털 트윈의 60FPS 초고속 실시간성" 간의 모순을 "CUDA 가속 3D Fourier Neural Operator (FNO) 대리모델"을 통해 완벽히 해결하는 차세대 엔지니어링 청사진을 제시한다.
Geminy는 본 계획에 따라 Phase 1의 오픈소스 OpenFOAM 17×17 핵연료 다발 메쉬 템플릿 및 차원 축약 인터페이스 정의부터 단계별로 차근차근 개발을 추진할 것을 선언한다.
발신: thesis.hyperbook.com & ROOPS Geminy 수석 원자력 열수력·AI 아키텍처팀