장비의 형상과 경계 조건에 의해 정의된 고유값 문제에서 기저 함수를 선택하십시오. 직사각형 파일럿 플랜트 영역의 경우, 일반적으로 데카르트 좌표계의 곱 형태로 사인 및 코사인 고유 함수를 사용합니다. 원통형 영역의 경우, 반경 방향에는 베셀 함수를, 각도 방향에는 푸리에 사인 또는 코사인 함수를, 축 방향에는 적절한 삼각 함수를 사용합니다. 물리적 경계 조건의 동차(homogeneous) 형태를 만족하는 함수가 올바른 함수이며, 이 함수들이 팽창 계수(expansion coefficients) 계산을 위한 직교 기저를 형성합니다.

형상은 고유 함수군을 결정하고, 경계 조건은 허용 가능한 고유값과 특정 사인, 코사인 또는 베셀 조합을 결정합니다. 형상과 경계 조건 모두에 일치하는 기저 함수는 독립적이고 직교하는 모드를 생성하여 확산 또는 전도 모델을 훨씬 쉽게 해결할 수 있게 합니다.

공간 고유값 문제로 시작하기

좌표계별로 지배 방정식 분리

변수 분리 후, 과도 전도 및 확산 모델은 일반적으로 다음과 같은 공간 고유값 문제로 축소됩니다.

[ \mathcal{L}_x\Phi=-\lambda^2\Phi, ]

여기서 (\mathcal{L}_x)는 열 또는 물질 전달 연산자의 공간 부분입니다.

적절한 기저 함수는 파일럿 플랜트의 경계 조건 하에서 이 연산자의 고유 함수입니다. 단순히 수학적으로 편리하다는 이유만으로 선택되는 것이 아닙니다.

동차 경계 조건에 기저 맞추기

경계 조건은 온도나 농도 값을 지정하거나, 플럭스를 지정하거나, 플럭스와 표면 값의 관계를 정의할 수 있습니다.

이는 대체로 다음과 같이 대응됩니다:

디리클레(Dirichlet): (\Phi=0), 정상 상태 성분이 제거된 후 고정된 값을 나타냄.
노이만(Neumann): (\Phi'=0), 단열 또는 제로 플럭스 경계를 나타냄.
로빈(Robin): (\Phi'+A\Phi=0), 대류 또는 계면 전달을 나타냄.

비동차 조건의 경우, 필드는 종종 정상 상태 또는 리프팅(lifting) 해와 과도 나머지(transient remainder)로 분해됩니다. 그런 다음 고유 함수는 나머지 부분의 동차 경계 조건을 만족하도록 선택됩니다.

직사각형 형상을 위한 함수 선택

데카르트 방향에서 사인 및 코사인 고유 함수 사용

1차원 직사각형 영역의 경우, 공간 방정식은 다음과 같은 형태를 가집니다.

[ \frac{d^2\Phi}{dx^2}=-\lambda^2\Phi, ]

일반 해는

[ \Phi(x)=a\sin(\lambda x)+b\cos(\lambda x). ]

경계 조건에 따라 어떤 조합이 허용되는지, 그리고 어떤 이산적인 (\lambda) 값이 양쪽 경계를 모두 만족하는지가 결정됩니다.

다차원 직사각형 장비를 위한 곱 사용

직사각형 슬래브, 채널, 촉매 섹션 또는 평판 구성 요소의 경우, 분리된 모드는 일반적으로 다음과 같은 형태를 가집니다.

[ \Phi_{mn}(x,y)=X_m(x)Y_n(y), ]

또는 3차원의 경우,

[ \Phi_{mnp}(x,y,z)=X_m(x)Y_n(y)Z_p(z). ]

각 요소는 해당 표면 쌍의 경계 조건으로부터 독립적으로 선택됩니다. 1차원 요소들이 직교하기 때문에 전체 모드도 직교합니다.

경계 조건을 적용하여 사인 또는 코사인 모드 선택

(x=0)에서의 제로 값 경계는 코사인 항을 제거하여 사인 모드를 생성합니다. (x=0)에서의 제로 플럭스 경계는 사인 항을 제거하여 코사인 모드를 생성합니다.

예를 들어, 구간 (0\leq x\leq L)에서:

(\Phi(0)=\Phi(L)=0)은 (\sin(n\pi x/L))과 같은 사인 모드로 이어집니다.
(\Phi'(0)=\Phi'(L)=0)은 상수 제로 고유값 모드를 포함하여 (\cos(n\pi x/L))과 같은 코사인 모드로 이어집니다.
한쪽 끝은 고정 조건, 다른 쪽 끝은 제로 플럭스 조건인 경우, 일반적으로 (\pi/L)의 반정수 배를 포함하는 이동된 고유값 시퀀스가 생성됩니다.
대류 경계는 가장 단순한 정수 배가 아닌 특성 방정식에 의해 결정되는 로빈 고유값을 생성합니다.

정확한 고유값은 장비의 치수와 경계 매개변수에 따라 달라집니다.

원통형 형상을 위한 함수 선택

반경 방향 변화에 베셀 함수 사용

원통 좌표계에서 확산 또는 전도 연산자의 반경 부분은 (1/r) 인자를 포함합니다. 따라서 반경 방향 고유값 방정식은 데카르트 방정식과 다르며 베셀 함수로 이어집니다.

축대칭 문제의 경우, 반경 방향 해는 일반적으로 다음을 기반으로 합니다.

[ J_0(\lambda r), ]

여기서 (J_0)는 0차 제1종 베셀 함수입니다. 영역과 경계 조건에 따라 다른 조합이나 차수가 필요할 수 있습니다.

축을 포함하는 고체 원통은 일반적으로 (r=0)에서 유한하지 않기 때문에 특이 베셀 해 (Y_\nu(\lambda r))를 제외합니다. 반면, 환형 영역은 도메인이 축을 포함하지 않으므로 (J_\nu)와 (Y_\nu)의 선형 조합이 필요할 수 있습니다.

비축대칭 공정일 때 각도 푸리에 함수 추가

온도나 농도가 각도 좌표에 따라 변하는 경우, 각도 고유 함수는 다음과 같습니다.

[ \cos(m\theta) \quad\text{및}\quad \sin(m\theta), ]

여기서 (m)은 정수 방위각 모드 번호입니다.

완전한 축대칭 관형 반응기, 쉘 또는 컬럼 모델의 경우 (m=0) 모드만 유지되므로 각도 의존성이 사라집니다. 비축대칭 가열, 벽면 조건 또는 공급 분포는 더 높은 각도 모드를 필요로 합니다.

원통 축을 따라 삼각 함수 사용

유한 원통 또는 관형 파일럿 유닛의 경우, 축 좌표는 일반적으로 입구, 출구, 끝 벽 또는 대칭 조건에 따라 선택된 사인 및 코사인 고유 함수를 생성합니다.

따라서 분리된 원통형 모드는 다음과 같은 구조를 가질 수 있습니다.

[ \Phi_{mnp}(r,\theta,z)

R_{mn}(r), \Theta_m(\theta), Z_p(z), ]

여기서 (R_{mn})은 베셀 기반 반경 함수, (\Theta_m)은 사인 또는 코사인 각도 함수, (Z_p)는 축 삼각 함수입니다.

경계 조건이 고유값을 결정하도록 하기

디리클레 조건은 고유 함수의 영점을 선택함

원통형 벽면에 지정된 온도나 농도가 있는 경우, 반경 방향 고유 함수는 균질화 후 해당 벽면에서 0이 되어야 합니다.

(\Phi(R)=0)인 고체 원통의 경우, 반경 방향 고유값은 관련 베셀 함수의 영점과 관련이 있습니다. 이는 두 개의 고정 값 경계가 있는 직사각형 영역에서 얻은 사인 고유값의 원통형 대응물입니다.

노이만 조건은 도함수의 영점을 선택함

단열 원통 벽면은 반경 방향 도함수가 0이어야 합니다:

[ \left.\frac{d\Phi}{dr}\right|_{r=R}=0. ]

허용 가능한 고유값은 함수 자체의 영점이 아니라 베셀 함수 도함수의 영점에 의해 결정됩니다.

동일한 원리가 직사각형 형상에도 적용되며, 여기서 제로 플럭스 조건은 코사인 유형 모드나 도함수 기반 특성 방정식을 선택합니다.

로빈 조건은 초월 특성 방정식을 생성함

대류 벽면 조건은 다음과 같은 형태를 가집니다.

[ -k\left.\frac{d\Phi}{dr}\right|_{r=R}

h\Phi(R), ]

또는 이에 상응하는 무차원 형태입니다. 원통에서 베셀 반경 방향 해를 대입하면 베셀 함수와 그 도함수를 모두 포함하는 특성 방정식이 생성됩니다.

이러한 고유값은 일반적으로 수치적으로 계산되어야 합니다. 이는 예상된 결과이며 해석적 확장을 무효화하지 않습니다.

직교성을 올바르게 유지하기

적절한 내적 사용

직사각형 고유 함수는 일반적으로 표준 데카르트 내적 하에서 직교합니다,

[ \langle f,g\rangle

\int f(x)g(x),dx. ]

다차원 직사각형 영역의 경우 내적은 전체 부피에 걸쳐 확장됩니다.

원통형 반경 함수는 반경 가중치 (r)와 함께 직교합니다:

[ \langle f,g\rangle_r

\int_0^R r,f(r)g(r),dr. ]

이 가중치는 원통형 부피 요소 (r,dr,d\theta,dz)에서 나오기 때문에 필수적입니다.

계수를 계산하기 전에 모드 정규화

고유 함수를 식별한 후, 다음과 같이 정규화하십시오.

[ \langle \Phi_i,\Phi_j\rangle=\delta_{ij}, ]

여기서 (\delta_{ij})는 크로네커 델타입니다. 초기 온도, 농도 또는 강제 분포 (u)에 대한 팽창 계수는 스칼라 투영에서 얻어집니다,

[ c_i=\langle u,\Phi_i\rangle. ]

원통형 모델의 경우, 이 투영은 올바른 반경 가중치와 해당되는 경우 각도 및 축 위치에 대한 적분을 포함해야 합니다.

물리학이 허용할 때 제로 모드 포함

순수 노이만 문제는 (\lambda=0)인 상수 고유 함수를 포함할 수 있습니다. 이 모드는 경계를 통해 순 열이나 물질 플럭스가 없을 때 공간 평균의 보존을 나타냅니다.

이를 생략하면 특히 단열 고체나 폐쇄 확산 시스템의 경우 잘못된 장기 해가 생성될 수 있습니다.

트레이드오프 이해하기

직사각형 기저는 구현하기 더 간단함

사인 및 코사인 기저는 표준 직사각형 경계 조건에 대해 간단한 고유값과 정규화 공식을 가집니다. 평판, 슬래브, 직사각형 덕트 및 데카르트 계산 영역에 대해 가장 효율적인 선택인 경우가 많습니다.

서로 다른 경계가 서로 다른 로빈 매개변수를 사용하거나 재료 특성이 공간적으로 변할 때 그 단순함은 감소합니다.

원통형 기저는 물리적으로 정확하지만 더 복잡함

베셀 함수는 원통형 장비의 반경 방향 확산과 변화하는 단면적을 정확하게 나타냅니다. 이를 데카르트 사인 함수로 대체하면 특히 축 근처나 반경 방향 기울기가 지배적일 때 반경 방향 해가 왜곡될 수 있습니다.

추가적인 복잡성은 베셀 특성 방정식을 풀고, 반경 방향 가중치를 적용하며, (r=0)에서의 정규성을 처리해야 한다는 점입니다.

일치하지 않는 함수는 경계 및 계수 오류를 생성함

필요한 동차 경계 조건을 만족하지 않는 기저는 부분 적분 후 원치 않는 경계 항을 남깁니다. 그러면 변환된 확산 또는 전도 방정식은 해결되지 않은 항이나 잘못된 강제 기여분을 포함하게 됩니다.

이는 단순히 우아함의 상실이 아닙니다. 표면 플럭스, 과도 감쇠율 및 장기 평형 예측을 변화시킬 수 있습니다.

이상적인 고유 함수 방법의 모델링 한계

표준 사인, 코사인 및 베셀 구성은 상수 또는 충분히 규칙적인 재료 특성을 포함하는 적절한 속성을 가진 선형 연산자를 가정합니다. 강하게 변하는 특성, 불규칙한 형상, 비선형 전달 계수 또는 복잡한 내부 소스는 유한 체적, 유한 요소 또는 수치 스펙트럼 기법과 같은 수치적 방법을 필요로 할 수 있습니다.

형상 기반 기저 선택은 그러한 경우에도 여전히 유용하지만, 완전한 정확한 해보다는 벤치마크, 축소 모델 또는 국소적으로 적절한 근사치 역할을 할 수 있습니다.

파일럿 플랜트에 적용하는 방법

경계 조건을 균질화하고 장비와 일치하는 좌표계로 공간 연산자를 작성하는 것부터 시작하십시오.

주요 초점이 직사각형 슬래브, 채널 또는 평판 유닛인 경우: 사인 및 코사인 고유 함수를 사용하고, 좌표 간의 곱을 취하며, 각 디리클레, 노이만 또는 로빈 경계를 만족하도록 각 요소를 선택하십시오.
주요 초점이 고체 원통형 또는 관형 유닛인 경우: 반경 방향 모드에는 베셀 함수를 사용하고, 축에서 특이 반경 방향 해를 버리고, 이를 축 사인 또는 코사인 모드와 결합하십시오.
주요 초점이 환형 또는 쉘 측 영역인 경우: 도메인이 축을 제외하는 곳에서 정규 및 특이 베셀 함수 분기를 모두 허용하고, 두 반경 방향 벽 조건을 모두 적용하여 조합과 고유값을 결정하십시오.
주요 초점이 축대칭 모델인 경우: 0 각도 모드만 유지하고 직교성을 적용하고 계수를 계산할 때 반경 가중치 (r)를 사용하십시오.
주요 초점이 대류 벽면 전달인 경우: 로빈 조건에서 특성 방정식을 형성하고 표준 푸리에 시퀀스를 강요하는 대신 수치적으로 고유값을 계산하십시오.
주요 초점이 단열 장비의 보존인 경우: 0 고유값을 명시적으로 확인하고 경계 조건이 허용할 때 상수 모드를 포함하십시오.

파일럿 플랜트 형상과 경계 조건 모두에 일치하는 기저를 선택하면 전달 모델이 유효하고 효율적인 모드 표현으로 바뀝니다.

요약 표:

파일럿 플랜트 형상	권장 기저 함수	주요 고려 사항
직사각형 슬래브, 채널 또는 평판	좌표 곱으로 결합된 사인 및 코사인 함수	디리클레, 노이만 또는 로빈 경계에 따라 각 요소를 선택하십시오.
고체 원통 또는 관형 유닛	축 사인 또는 코사인 함수와 결합된 반경 방향 베셀 함수	축에서 특이 베셀 해를 제외하고 올바른 반경 방향 고유값을 적용하십시오.
환형 또는 쉘 측 영역	제1종 및 제2종 베셀 함수의 조합	도메인이 축을 제외하므로 두 반경 방향 분기가 모두 필요할 수 있습니다.
축대칭 장비	0 각도 푸리에 모드	직교성 및 계수 계산에 반경 가중치 `r`을 사용하십시오.
대류 벽면	수치적으로 계산된 로빈 고유값을 가진 형상별 고유 함수	플럭스-값 관계를 적용하여 특성 방정식을 형성하십시오.
단열 장비	상수 모드를 포함한 코사인 또는 도함수 기반 모드	공간 평균 보존을 유지하기 위해 0 고유값을 유지하십시오.

파일럿 규모의 전달 또는 분리 시스템을 계획 중이신가요? 올바른 교육 및 직업 훈련 단위 조작 파일럿 플랜트를 선택하고 적용하는 데 대한 전문가 지원을 받으려면 LABPARK에 문의하십시오. LABPARK는 대학, 연구 기관 및 기업을 위한 화학 공학, 바이오 공정 및 생명 공학, 환경 및 수처리 파일럿 플랜트를 공급합니다. 당사의 솔루션은 실제 장비를 정확하게 모델링하고, 확산 및 전도 거동을 검증하며, 운영자와 학생을 교육하고, 형상에 적합한 설계, 신뢰할 수 있는 계측 및 실용적인 기술 지침을 통해 공정 연구를 가속화하도록 돕습니다.