이론적 근거는 확산 방정식과 그린 함수 중첩의 결합입니다. 연속적인 용질 방출은 무한히 작은 순간 점 방출의 연속으로 표현되며, 각각은 확산 방정식에 따라 확산됩니다. 이러한 순간적 해를 방출 이력에 걸쳐 적분하면 비정상 상태 농도장이 얻어지며, 장시간이 지나고 소스가 일정하며 주변 조건이 정상 유동을 지원할 때 정상 상태 분포가 형성됩니다.
연속 방출은 순간 확산 응답의 시간 적분입니다. 무한하고 정체된 3차원 매질에서 일정한 점 소스 방출의 경우, 시간이 충분히 흐르면 비정상 상태 해는 (C(r)=\dot{m}/(4\pi Dr))에 수렴합니다.
연속 방출 모델링 방법
지배 확산 방정식
점 소스가 있는 분자 확산의 경우, 농도는 픽의 제2법칙에 의해 지배됩니다:
[ \frac{\partial C}{\partial t}=D\nabla^2 C ]
여기서 (C)는 농도, (t)는 시간, (D)는 분자 확산 계수, (\nabla^2 C)는 농도장의 공간적 곡률을 나타냅니다.
원점에 위치한 소스의 경우, 구면 대칭성으로 인해 문제는 반경 거리 (r)로 축소됩니다. 따라서 농도는 용질이 소스를 둘러싼 구면 표면을 따라 어떻게 퍼지는지에 따라 결정됩니다.
순간 점 소스 해
질량 (m)의 순간 방출은 경계가 없는 3차원 매질에서 가우시안 그린 함수로 설명되는 확산 응답을 생성합니다:
[ C_{\text{inst}}(r,t)= \frac{m}{(4\pi Dt)^{3/2}} \exp\left(-\frac{r^2}{4Dt}\right) ]
이 해는 한순간에 방출된 유한한 질량을 나타냅니다. 질량이 점차 더 큰 부피를 차지함에 따라 피크 농도는 감소하면서 시간이 지남에 따라 외부로 퍼져 나갑니다.
방출 시간에 따른 중첩
확산은 주어진 가정하에 선형 전달 과정입니다. 따라서 여러 방출에 대한 응답은 개별 방출에 대한 응답의 합과 같습니다.
방출 속도가 (\dot{m}(\tau))일 때 (여기서 \tau는 각 미소 질량 증분이 방출되는 시간), 농도는 순간 해를 적분하여 얻습니다:
[ C(r,t)= \int_0^t \frac{\dot{m}(\tau)} {\left[4\pi D(t-\tau)\right]^{3/2}} \exp\left[ -\frac{r^2}{4D(t-\tau)} \right]d\tau ]
이 적분은 연속 방출로부터 비정상 상태 농도 분포를 계산하기 위한 이론적 토대입니다.
비정상 상태 농도 분포
일정한 속도의 방출
(t=0)에서 시작하여 일정한 속도 (\dot{m})로 질량을 연속적으로 방출하는 소스의 경우, 중첩 적분을 계산하면 다음과 같습니다:
[ C(r,t)= \frac{\dot{m}}{4\pi Dr} \operatorname{erfc} \left( \frac{r}{2\sqrt{Dt}} \right) ]
여기서 여오차 함수(\operatorname{erfc})는 시간 (t) 후에 거리 (r)에서 확산 전선이 얼마나 발달했는지를 설명합니다.
오차 함수 항의 물리적 의미
인자
[ \operatorname{erfc} \left( \frac{r}{2\sqrt{Dt}} \right) ]
는 정상 상태 반경 프로파일에 대한 시간 의존적 보정입니다. 이는 용질이 소스에서 관측 지점까지 확산되는 데 유한한 시간만이 경과했음을 고려합니다.
초기 시간대에는 특정 거리에서 (r/(2\sqrt{Dt}))가 크기 때문에 여오차 함수는 작습니다. 따라서 농도는 최종 정상 상태 값보다 훨씬 낮습니다.
확산 시간 척도
관련된 비교는 경과 시간 (t)과 특성 확산 시간 사이의 비교입니다:
[ t_D\sim\frac{r^2}{D} ]
(t)가 (r^2/D)에 비해 작을 때 측정값은 강한 비정상 상태를 나타냅니다. (t)가 (r^2/D)보다 클 때, 해당 거리의 농도는 정상 상태 값에 가깝습니다.
정상 상태 농도 분포
장시간 극한값 구하기
(t)가 무한대로 접근함에 따라,
[ \frac{r}{2\sqrt{Dt}}\rightarrow 0 ]
이고
[ \operatorname{erfc}(0)=1 ]
이므로 비정상 상태 식은 다음과 같이 됩니다:
[ C(r)=\frac{\dot{m}}{4\pi Dr} ]
이것이 무한 3차원 매질에서 일정한 점 방출에 대한 정상 상태 농도 분포입니다.
플럭스 해석
정상 상태 결과는 질량 보존 법칙에서 직접 도출됩니다. 반경 (r)에서 총 방출 속도는 면적 (4\pi r^2)인 구면을 통과해야 합니다.
따라서 반경 방향 확산 플럭스는 다음과 같습니다:
[ J_r=\frac{\dot{m}}{4\pi r^2} ]
픽의 제1법칙을 사용하면,
[ J_r=-D\frac{dC}{dr} ]
이고, 소스에서 멀어질수록 농도가 0에 수렴한다는 조건으로 적분하면 다음을 얻습니다:
[ C(r)=\frac{\dot{m}}{4\pi Dr} ]
(1/r) 의존성은 팽창하는 구면 표면에 걸친 기하학적 희석을 반영합니다.
정상 상태가 유한 질량 평형이 아닌 이유
정상 상태란 용질이 움직임을 멈췄다는 의미가 아닙니다. 외부 확산 플럭스가 일정하기 때문에 각 위치의 농도가 더 이상 시간에 따라 변하지 않는다는 의미입니다.
소스는 계속해서 질량을 추가하며, 이 질량은 모델링된 영역을 떠나거나 주변 매질로 무한히 확장되어야 합니다. 출구, 경계 싱크, 반응 또는 기타 제거 메커니즘이 없는 경우, 폐쇄된 파일럿 시스템에서 진정한 정상 상태는 물리적으로 도달하지 못할 수 있습니다.
파일럿 실험 측정값 해석
관측 시간과 거리 비교
측정값을 비정상 상태로 처리할지 정상 상태로 처리할지는 경과 시간과 방출 지점으로부터의 거리에 따라 달라집니다.
소스에 가까운 위치는 비교적 빠르게 정상 상태 프로파일에 도달할 수 있지만, 더 먼 위치는 동일한 실험 동안 비정상 상태로 남아 있을 수 있습니다.
관련 전달 영역 추정
무차원 양
[ \eta=\frac{r}{2\sqrt{Dt}} ]
은 직접적인 진단 도구를 제공합니다.
큰 (\eta) 값은 측정 위치에서 확산 전선이 완전히 발달하지 않았음을 나타냅니다. 작은 (\eta) 값은 비정상 상태 보정이 약하며 정상 상태 근사가 점점 더 적절함을 나타냅니다.
전체 비정상 상태 모델 사용
소스 가동 중에 농도 데이터를 수집할 때, 방출 지속 시간이 확산 시간과 비슷할 때, 또는 확산 전선이 정상 상태에 도달하지 못할 정도로 공간적으로 멀리 떨어져 있을 때는 비정상 상태 방정식을 사용해야 합니다.
또한 돌파 시간(breakthrough timing), 농도 축적, 또는 파일럿 시스템이 정상 상태 프로파일에 도달하는 속도를 추정하는 것이 목표일 때 적절한 모델입니다.
트레이드오프 이해
방정식 이면의 가정
제시된 해는 점 소스, 일정한 확산 계수, 구면 대칭, 무한하거나 충분히 큰 매질을 가정하며, 이류(advection), 반응, 흡착 또는 경계 효과가 없다고 가정합니다.
실제 파일럿 장비는 이러한 가정 중 일부를 위반할 수 있습니다. 유동, 용기 벽, 혼합, 유한한 기하학적 구조 및 불균일한 재료 특성은 농도 분포를 실질적으로 변화시킬 수 있습니다.
파일럿 유닛에서 "정상 상태"의 의미
식 (C(r)=\dot{m}/(4\pi Dr))은 이상적인 확산 제어 정상 상태 장을 설명합니다. 이를 유한한 반응기나 실험 용기의 최종 상태로 자동으로 해석해서는 안 됩니다.
경계가 용질을 반사하거나, 시스템이 닫혀 있거나, 방출이 중단되면 장시간 거동은 무한 매질 점 소스 해와 다를 수 있습니다.
연속 방출 대 펄스 방출
연속 방출은 이전 모든 방출의 효과를 누적하며, 이는 여오차 함수 해를 생성합니다. 반면 펄스 방출은 순간 가우시안 해를 따르며 일반적으로 정상 상태 반경 프로파일이 아닌 이동하는 농도 최대치를 가집니다.
따라서 연속 소스 방정식을 유한 펄스에 적용하면 농도 크기와 시간 의존성을 모두 잘못 나타낼 수 있습니다.
확산 대 이류
이 공식은 분자 확산만을 설명합니다. 벌크 유체 운동이 중요한 경우, 전달에는 이류 항도 포함되어야 하며 농도장은 순수 확산 방정식이 아닌 이류-확산 방정식에 의해 지배될 수 있습니다.
목표에 맞는 올바른 선택
물리적 질문과 파일럿 실험의 시간 척도에 맞는 모델을 사용하십시오.
- 비정상 상태 거동이 주된 관심사인 경우: (C(r,t)=\frac{\dot{m}}{4\pi Dr}\operatorname{erfc}\left(\frac{r}{2\sqrt{Dt}}\right))를 사용하고 경과 시간과 (r^2/D)를 비교하십시오.
- 장시간 반경 분포가 주된 관심사인 경우: 관련 위치가 거의 정상 상태에 도달했음을 확인한 후에만 (C(r)=\frac{\dot{m}}{4\pi Dr})를 사용하십시오.
- 모델 타당성이 주된 관심사인 경우: 이류, 경계, 반응, 유한 소스 크기 또는 불일정한 방출 속도가 더 완전한 전달 모델을 요구하는지 확인하십시오.
- 측정값 해석이 주된 관심사인 경우: 가동 또는 돌파 효과가 보일 때 비정상 상태 해를 먼저 피팅한 다음, 데이터가 정상 상태 근사를 지원하는지 평가하십시오.
구분은 전반적으로 동일한 원리에 의해 지배됩니다. 연속 방출 농도장은 시간 적분된 확산 응답이며, 정상 상태는 그들의 장시간, 일정 플럭스 극한값입니다.
요약 표:
| 모델 또는 개념 | 핵심 방정식 | 최적 사용처 |
|---|---|---|
| 확산 방정식 | ∂C/∂t = D∇²C | 분자 확산 설명 |
| 순간 방출 | 가우시안 그린 함수 | 펄스 또는 유한 질량 방출 모델링 |
| 연속 방출 | C(r,t) = [ṁ/(4πDr)] erfc[r/(2√Dt)] | 비정상 상태 축적 및 돌파 예측 |
| 정상 상태 방출 | C(r) = ṁ/(4πDr) | 장시간 반경 농도 프로파일 추정 |
| 확산 시간 척도 | tD ≈ r²/D | 측정값이 비정상 상태인지 정상 상태인지 판별 |
| 타당성 검사 | 이류, 경계, 반응, 흡착, 소스 기하학 고려 | 필요 시 더 완전한 전달 모델 선택 |
LABPARK으로 확산 공정 모델링
비정상 상태 농도 축적, 돌파 시간 또는 정상 상태 전달에 대한 신뢰할 수 있는 데이터가 필요하십니까? LABPARK에 문의하여 귀하의 응용 분야에 맞춘 교육용 또는 직업 훈련용 단위 조작 파일럿 플랜트에 대해 논의하십시오. LABPARK은 화학 공학, 바이오 공정 및 생명 공학, 환경 및 수처리 분야의 파일럿 플랜트로 대학, 연구 기관 및 기업을 지원합니다. 당사의 시스템은 실용적이고 제어된 실험을 사용하여 확산, 물질 전달, 이류-확산 거동, 경계 효과 및 공정 스케일업을 조사하는 데 도움을 줍니다.
관련 제품
- 입자 내 확산 유효인자 측정을 위한 교육용 단위조작 파일럿 플랜트
- 다기능 특수 증류 교육용 파일럿 플랜트
- 연속식 배치 추출 증류 교육용 파일럿 플랜트
- 실습용 단위 조작을 위한 이중 모드 정류 파일럿 플랜트
- 다단 교반 탱크 직렬 체류 시간 분포 및 혼합 성능 결정 교육용 파일럿 플랜트
