수학적 기저 함수는 각도 방향의 구면 조화 함수와 반경 방향의 구면 베셀 함수이며, 장비의 경계 조건을 만족하도록 선택됩니다. 축대칭 구형 모듈의 경우 각도 의존성은 르장드르 다항식(P_\ell(\cos\theta))으로 전개하고, 반경 의존성에는 구면 베셀 함수 또는 이에 상응하는 사인-반경 나눗셈 형태를 사용합니다. 허용되는 고유값은 구의 표면 값이 고정되어 있는지, 플럭스가 0인지, 또는 대류 열전달이나 물질전달이 발생하는지에 따라 결정됩니다.

지배적인 원리는 구면 미분 연산자와 물리적 경계 조건을 모두 만족하는 직교 고유함수 집합을 선택하는 것입니다. 디리클레, 노이만 및 로빈 조건은 서로 다른 특성 방정식, 고유값 및 전개 계수를 생성합니다.

구형 수송 문제의 정식화 방법

열확산과 물질확산은 동일한 수학적 구조를 공유합니다

물성이 일정하고 내부 발생이 없는 경우, 온도 또는 농도는 일반적인 장(\Theta(r,\theta,\phi,t))으로 나타낼 수 있습니다. 과도 방정식은 다음과 같은 확산 형태를 가집니다.

[ \frac{\partial \Theta}{\partial t}

D\nabla^2\Theta, ]

여기서 (D)는 해당 확산계수입니다.

열전도의 경우,

[ D=\alpha=\frac{k_T}{\rho \hat C_p}, ]

이며, 물질확산의 경우 (D)는 종의 확산계수입니다.

구면 좌표계에서는 자연스럽게 변수분리가 이루어집니다

구면 좌표계의 라플라스 연산자는 반경 방향과 각도 방향으로 분리됩니다. 변수분리해는 개략적으로 다음과 같이 쓸 수 있습니다.

[ \Theta(r,\theta,\phi,t)

R(r)Y(\theta,\phi)G(t). ]

시간 의존 인자는 일반적으로 다음과 같은 형태를 가집니다.

[ G(t)=\exp(-D\lambda^2t), ]

따라서 공간 고차 모드는 더 빠르게 감쇠합니다.

축대칭은 각도 기저를 단순화합니다

장비와 부하가 방위각(\phi)에 무관하다면 각도 문제는 다음 변수에만 의존합니다.

[ \eta=\cos\theta. ]

이때 각도 방정식은 르장드르 방정식이 되며, 정칙해는 다음의 르장드르 다항식입니다.

[ P_\ell(\eta)=P_\ell(\cos\theta). ]

이 다항식은 (-1\leq\eta\leq1)에서 직교하므로, 장은 다음과 같이 전개할 수 있습니다.

[ \Theta(r,\theta,t)

\sum_{\ell=0}^{\infty} \sum_{n} A_{\ell n} R_{\ell n}(r) P_\ell(\cos\theta) e^{-D\lambda_{\ell n}^2t}. ]

완전한 3차원 부하의 경우 르장드르 다항식은 구면 조화 함수(Y_\ell^m(\theta,\phi))로 대체됩니다.

어떤 반경 방향 기저 함수를 사용하는가?

구면 베셀 함수는 일반적인 반경 방향 기저입니다

반경 방향 고유값 방정식에는 각도 분리 지수(\ell)가 포함됩니다. 구의 중심에서 정칙인 해는 다음의 구면 베셀 함수입니다.

[ R_\ell(r)=j_\ell(\lambda r). ]

두 번째 독립해인 구면 노이만 함수(y_\ell(\lambda r))는 (r=0)에서 특이하므로 원점을 포함하는 고체 구에서는 제외됩니다.

구면 대칭 모드(\ell=0)의 경우,

[ j_0(\lambda r)=\frac{\sin(\lambda r)}{\lambda r}. ]

따라서 일반적으로 사용되는 형태인

[ R_0(r)\propto\frac{\sin(\lambda r)}{r} ]

는 특히 (\ell=0) 반경 방향 해이거나 다음 변환의 결과입니다.

[ \Theta(r,t)=\frac{\Psi(r,t)}{r}. ]

(\Psi/r) 변환은 반경 방향 방정식을 단순화합니다

구면 대칭 문제에서 반경 방향 확산 방정식은 다음과 같습니다.

[ \frac{\partial \Theta}{\partial t}

D\left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial\Theta}{\partial r} \right) \right]. ]

다음을 사용하면

[ \Theta(r,t)=\frac{\Psi(r,t)}{r} ]

구면 반경 연산자는 1차원 확산 방정식으로 변환됩니다.

[ \frac{\partial\Psi}{\partial t}

D\frac{\partial^2\Psi}{\partial r^2}. ]

그 결과 얻어지는 공간 고유함수는 (\sin(\lambda r))와 (\cos(\lambda r))의 조합입니다. 일반적인 고체 구 정식화에서는 중심에서의 정칙성으로 인해 코사인 항이 제거됩니다.

직교성은 급수해를 가능하게 합니다

반경 방향 및 각도 방향 고유함수는 적절한 구면 가중치 아래에서 직교 집합을 이룹니다. 반경 방향 모드의 경우 체적 요소에 의해 (r^2) 가중치가 도입됩니다.

[ \int_0^l r^2 R_{\ell n}(r)R_{\ell m}(r),dr=0, \qquad n\ne m. ]

이 직교성을 이용하면 초기 온도 또는 농도 분포를 고유함수에 투영하여 계수 (A_{\ell n})을 구할 수 있습니다.

경계 조건이 고유값을 결정하는 방법

중심에는 정칙성 또는 대칭 조건이 필요합니다

고체 구의 중심에서 장은 유한하게 유지되어야 합니다. 물리적인 반경 방향 플럭스 또한 0이어야 합니다.

[ \left.\frac{\partial\Theta}{\partial r}\right|_{r=0}=0. ]

이는 구면 대칭 조건입니다. 이는 선택적인 수치 편의가 아니라 중심을 통한 연속성을 보장하고 특이한 반경 방향 해를 제거하는 조건입니다.

디리클레 조건은 고정된 표면 값을 부과합니다

고정된 표면 온도 또는 농도는 다음과 같이 나타냅니다.

[ \Theta(l,t)=\Theta_s. ]

정상 상태의 경계값을 뺀 후, 균질 고유함수 조건은 다음과 같이 됩니다.

[ R_{\ell n}(l)=0. ]

구면 대칭 모드의 경우 이는 다음을 의미합니다.

[ \sin(\lambda_n l)=0, ]

따라서

[ \lambda_n l=n\pi, ]

가 되며, 0 모드를 포함하는지 여부에 따라 인덱싱을 선택합니다.

이 조건은 제어가 잘되는 욕조, 재킷 또는 외부 농도 저장조에 의해 표면이 일정하게 유지되는 구를 모델링합니다.

노이만 조건은 지정된 플럭스를 부과합니다

단열 또는 불투과 표면은 법선 방향 플럭스가 0입니다.

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}=0 ]

이는 열전달의 경우에 해당하며, 물질확산의 경우에는 다음과 같습니다.

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}=0 ]

변환된 구면 대칭 해에 대해 이 조건은 다음 특성 방정식을 생성합니다.

[ \lambda_n l=\tan(\lambda_n l). ]

보다 일반적으로는 해당 구면 베셀 함수의 도함수가 0이어야 합니다.

[ \left.\frac{d}{dr}j_\ell(\lambda_{\ell n}r)\right|_{r=l}=0. ]

단열 또는 불투과 경계는 내부 발생원이나 소멸원이 있는 경우를 제외하면 구 내부의 총 열에너지 또는 종 질량을 보존합니다.

로빈 조건은 대류 전달을 나타냅니다

대류 표면은 혼합 조건으로 설명됩니다. 열전달의 경우,

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}

h_T\left[\Theta(l,t)-\Theta_\infty(t)\right]. ]

물질전달의 경우 이에 상응하는 조건은 다음과 같습니다.

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}

k_m\left[C(l,t)-C_\infty(t)\right]. ]

무차원화 후 전달계수는 비오형 매개변수를 통해 나타납니다. 변환된 구면 대칭 해에 대한 특성 방정식은 다음과 같은 형태로 쓸 수 있습니다.

[ \lambda_n l

\left(1-A l^2\right)\tan(\lambda_n l), ]

(A)의 정의, 부호 규약 및 변환된 종속 변수에 따라 달라집니다.

중요한 물리적 의미는 로빈 조건이 고체 외부의 저항을 고려한다는 점입니다. 여기에는 유체 경계층, 재킷 측 필름 또는 외부 물질전달 필름이 포함됩니다.

비균질 조건에는 정상 상태 조정이 필요합니다

고유함수 전개는 경계 조건이 균질일 때 가장 간단합니다. 표면 온도 또는 농도가 시간에 따라 변하는 경우, 장은 일반적으로 특수해 또는 준정상 부분과 균질 조건을 만족하는 과도 부분으로 분해합니다.

이 단계는 확산 방정식을 변환할 때 경계항이 미해결 상태로 남는 것을 방지합니다.

기저 함수가 경계 조건과 일치해야 하는 이유

경계 조건이 고유값 문제를 정의합니다

미분 연산자만으로는 기저가 결정되지 않습니다. 연산자와 중심 조건, 표면 조건 및 모든 대칭 제약 조건이 함께 Sturm-Liouville 고유값 문제를 정의합니다.

따라서 동일한 물성을 가진 두 구라도 하나는 고정 표면 온도를 가지고 다른 하나는 대류 가열을 받는다면 서로 다른 고유값이 필요할 수 있습니다.

일치하는 기저 함수는 불필요한 경계항을 제거합니다

유한 푸리에 변환 또는 고유함수 변환을 부분적분으로 유도할 때 경계항이 나타납니다. 동일한 균질 경계 조건을 만족하는 기저를 사용하면 이러한 항이 사라지거나 알려진 상수로 변환됩니다.

예를 들어 대류 경계에 사인 기저를 사용하면 실제 표면 물리를 올바르게 나타내지 못하는 경우가 일반적입니다. 그 결과 변환된 문제의 고유값, 계수 또는 둘 다 잘못될 수 있습니다.

경계 조건은 실제 모듈을 반영해야 합니다

파일럿 플랜트 장비에서 수학적 선택은 물리적 계면을 따라야 합니다.

제어된 욕조 또는 지정된 표면 조성은 디리클레 거동을 시사합니다.
단열 벽, 대칭면 또는 불투과 표면은 노이만 거동을 시사합니다.
재킷, 유동 기체, 액체 필름 또는 외부 벌크상은 일반적으로 로빈 거동을 시사합니다.
입구, 출구 또는 반응성 계면에는 구면 표면 조건이 아닌 별도의 지정값 또는 플럭스 조건이 필요할 수 있습니다.

상충 관계 이해하기

폐형식 고유값에는 한계가 있습니다

디리클레 문제와 일부 노이만 문제는 간단한 특성 방정식과 명시적인 고유값 수열을 생성합니다. 로빈 문제는 일반적으로 수치적 근 찾기가 필요한 초월 방정식으로 이어집니다.

외부 전달 저항이 공정에 실질적인 영향을 미치는 경우에는 추가적인 수치 작업이 정당화됩니다.

사인-반경 나눗셈 형태는 적용 범위가 제한적입니다

반경 방향 해를 다음과 같이 쓰는 것은

[ \frac{\sin(\lambda r)}{r} ]

구면 대칭 (\ell=0) 모드에 적합합니다. 이는 임의의 각도 변화를 위한 완전한 기저가 아니며, 이 경우에는 (j_\ell(\lambda r))를 사용해야 합니다.

비축대칭 또는 매우 불균일한 부하에 단순한 형태를 사용하면 물리적으로 중요한 각도 모드를 누락할 수 있습니다.

이상적인 가정은 파일럿 장비에서 성립하지 않을 수 있습니다

표준 해석적 기저는 일정한 물성, 구형 기하 및 선형 확산 또는 전도 방정식과 같은 조건을 가정합니다. 실제 모듈에는 가변 확산계수, 내부 반응, 유체 운동, 접촉 저항, 불균일한 표면 전달 또는 완전하지 않은 구형 형상이 포함될 수 있습니다.

이러한 영향에는 유한체적법, 유한요소법 또는 스펙트럼 계산과 같은 수치 방법이 필요할 수 있습니다.

실험적 경계 조건은 거의 완벽하지 않습니다

고정된 욕조 온도라고 가정한 값도 변할 수 있으며, 추정된 필름 계수는 유량, 점도 또는 농도에 따라 달라질 수 있습니다. 따라서 모델 검증에서는 가정한 경계 조건을 정확한 것으로 취급하기보다 측정된 표면 및 내부 데이터를 비교해야 합니다.

경계 조건은 추정된 확산계수 또는 전달계수에서 가장 큰 불확실성의 원인 중 하나인 경우가 많습니다.

목표에 맞는 올바른 선택

기저는 파일럿 모듈의 기하, 대칭성 및 계면 물리로부터 선택해야 합니다.

주요 관심사가 구면 대칭 과도 거동인 경우: (\ell=0) 구면 베셀 기저 또는 이에 상응하는 (\sin(\lambda r)/r) 표현을 사용하고, 중심 대칭 조건과 적절한 표면 특성 방정식을 적용하십시오.
주요 관심사가 각도 방향으로 불균일한 가열 또는 농도인 경우: 축대칭 시스템에서는 르장드르 다항식으로, 완전한 3차원 시스템에서는 구면 조화 함수로 해를 전개하고, 반경 방향에는 (j_\ell(\lambda r))를 사용하십시오.
주요 관심사가 제어된 표면 온도 또는 농도인 경우: 디리클레 조건을 적용하고, 지정된 경계값을 뺀 후 표면에서 0이 되는 고유함수를 사용하십시오.
주요 관심사가 단열, 불투과성 또는 대칭인 경우: 법선 플럭스가 0인 노이만 조건을 적용하고, 도함수와 일치하는 고유함수를 사용하십시오.
주요 관심사가 재킷 측 또는 필름 지배 전달인 경우: 측정되거나 추정된 열전달계수 또는 물질전달계수를 사용하여 로빈 조건을 적용한 다음, 고유값을 수치적으로 결정하십시오.
주요 관심사가 파일럿 플랜트 측정값과의 일치인 경우: 수송 매개변수를 추정하기 전에 가정한 경계 조건이 실제 입구, 출구, 벽, 대칭 및 외부 유체 조건을 반영하는지 확인하십시오.

구면 베셀 함수 또는 조화 함수 기저를 물리적 경계 조건에 올바르게 맞추면 과도 모델은 단순한 형식적 수학 해가 아니라 파일럿 플랜트를 타당하게 나타내는 모델이 됩니다.

요약 표:

항목	핵심 선택	적용
각도 기저	르장드르 다항식 또는 구면 조화 함수	축대칭 또는 3차원 부하
반경 방향 기저	구면 베셀 함수(j_\ell)	고체 구 내부의 정칙해
구면 대칭	(\sin(\lambda r)/r) 형태	(\ell=0) 과도 거동
디리클레 조건	고정된 표면 온도 또는 농도	제어된 욕조 또는 저장조
노이만 조건	0 또는 지정된 플럭스	단열 또는 불투과 표면
로빈 조건	대류 열전달 또는 물질전달	재킷, 필름 및 유동상
고유값	특성 방정식의 근	경계 조건에 의해 결정
모델 검증	예측 데이터와 측정 데이터 비교	파일럿 플랜트 매개변수 추정

파일럿 장비의 과도 열전달 또는 물질전달을 위한 신뢰할 수 있는 모델이 필요하십니까? LABPARK에 문의하여 귀하의 적용 분야를 상담해 보십시오. LABPARK는 화학공학, 바이오공정 및 생명공학, 환경 및 수처리를 위한 교육·직업훈련용 단위조작 파일럿 플랜트를 제공합니다. 당사의 시스템은 대학, 연구기관 및 기업이 구형 수송 거동을 연구하고, 열전달 및 물질전달 성능을 평가하며, 수학적 모델을 실제 실험과 연결할 수 있도록 지원합니다. 귀하의 공정 목표에 맞는 적절한 모듈을 선택하고, 현실적인 경계 조건을 정의하며, 교육 또는 연구용 장치를 구축할 수 있도록 도와드립니다.