실리콘 주조 공정의 대형화가 산소 및 탄소 불순물 분포에 미치는 영향에 관한 수치 해석
Numerical Analysis of the Influence of Size Upgrading on Oxygen and Carbon Impurities in Casting Silicon
본 논문은 태양전지용 실리콘 잉곳 제조 공정인 방향성 응고법(Directional Solidification, DS)에서 생산 효율 증대와 원가 절감을 위해 필수적인 ‘장비 대형화(Size Upgrading)’가 품질의 핵심인 산소(O) 및 탄소(C) 불순물 분포에 미치는 영향을 수치 시뮬레이션을 통해 체계적으로 분석하였습니다. 연구진은 기존 G6 급에서 G7 급으로의 대형화 과정(잉곳 크기 확장 1124mm → 1320mm)에서 발생하는 고온 열장 및 유동장의 변화를 과도 상태(Transient phase) 글로벌 시뮬레이션으로 계산하였습니다. 아르곤 가스의 대류, 실리콘 용융액의 유동, 그리고 도가니 벽면 및 흑연 부품과의 화학 반응을 포함한 복합 모델을 구축하여 분석한 결과, 대형화된 G7 로(Furnace)에서는 G6 대비 실리콘 결정 내 평균 산소 및 탄소 농도가 각각 6.7% 및 7.3% 감소하는 긍정적인 효과가 나타남을 확인하였습니다. 본 연구는 대형화 과정에서 결정-용융액 계면의 불순물 편석 현상과 아르곤 가스 유동 구조의 최적화가 고품질 실리콘 잉곳 제조에 결정적인 역할을 함을 입증하는 정량적 지표를 제시합니다.
Metadata and Keywords
Paper Metadata
- Industry: 태양에너지 및 반도체 소재 산업 (Photovoltaics, Silicon Ingot Manufacturing)
- Material: 다결정/단결정 실리콘 (Grystalline Silicon)
- Process: 방향성 응고 공정 (Directional Solidification Process, DS), 수치 시뮬레이션
- System: G6 및 G7 급 방향성 응고로 (DS Furnace System)
- Objective: 장비 대형화에 따른 산소(O) 및 탄소(C) 불순물 제어 및 분포 최적화
Keywords
- 방향성 응고법 (Directional Solidification Method)
- 실리콘 (Silicon)
- 대형화 (Size Upgrade)
- 산소 및 탄소 불순물 (Oxygen and Carbon Impurities)
- 수치 시뮬레이션 (Numerical Simulation)
- 태양전지 (Solar Cells)
Executive Summary
Research Architecture
연구진은 G6(1124mm)와 G7(1320mm) 두 가지 크기에 대해 열전달, 유체 유동, 가스 대류 및 화학 반응이 결합된 2차원 과도 상태 글로벌 모델을 구축하였습니다. 시뮬레이션 격자수는 G6 34,248개, G7 42,994개로 구성하여 공간 해상도를 확보하였으며, 성장 단계별(결정화율 20%~100%) 불순물 거동을 정밀 추적하였습니다.
Method Snapshot
실리콘 용융액 내 산소의 용해(도가니 벽면), 아르곤 가스에 의한 SiO 증발 및 운송, 흑연 부품과의 반응을 통한 CO 형성, 그리고 CO의 용융액 재용해 과정을 수학적으로 모델링하였습니다. 특히 탄소의 낮은 편석 계수(0.07)와 산소의 편석 계수(1.25)를 적용하여 계면에서의 농도 변화를 계산하였습니다.
Key Findings
분석 결과, G7 로에서는 아르곤 가스 통로의 공간적 제한으로 인해 SiO 가스의 배출이 원활해지고 CO 가스의 재유입(Re-flux)이 억제됨을 확인하였습니다. 이로 인해 결정 내 평균 불순물 농도가 G6 대비 약 7% 내외로 감소하였습니다. 또한, 결정이 성장함에 따라 산소 농도는 점차 감소하는 반면, 탄소 농도는 용융액 내 축적 효과로 인해 점차 증가하는 상반된 경향성을 보였습니다.
Industrial Applications
본 연구 결과는 차세대 대면적 실리콘 잉곳 제조 장비 설계 시 핫존(Hot-zone) 구조 및 아르곤 가스 유동 가이드라인으로 활용될 수 있습니다. 특히 대형 도가니 사용 시 수평 방향의 농도 불균일성을 완화하기 위한 기계적 장치나 공정 변수 조절의 필요성을 시사합니다.
Limitations and Cautions
수치 모델은 2차원 축대칭(Axisymmetric) 가정을 기반으로 하므로, 실제 3차원적 유동의 복잡성을 완벽히 반영하지 못할 수 있습니다. 또한, 도가니 내벽 코팅의 상태나 원소재의 초기 순도 등 외부 환경 변수에 따라 실제 불순물 농도는 계산값과 차이가 발생할 수 있으므로 실험적 검증 작업이 병행되어야 합니다.
Paper Details
1. Overview
- Title: Numerical Analysis of the Influence of Size Upgrading on Oxygen and Carbon Impurities in Casting Silicon
- Author: Wenjia Su, Zhen Zhang, Jiulong Li, Zhicheng Guan, Jiaqi Li
- Year: 2023
- Journal: Silicon (Springer)
- DOI/Link: 10.1007/s12633-023-02323-8
2. Abstract
실리콘 잉곳의 방향성 응고 과정에서 장비 대형화는 생산 능력 증대와 비용 절감의 핵심 수단입니다. 본 연구에서는 G6 및 G7 방향성 응고로에서의 산소 및 탄소 불순물 분포를 연구하였습니다. 과도 상태 글로벌 시뮬레이션 기법을 사용하여 가스 및 용융액의 대류, 화학 반응, 결정-용융액 계면에서의 불순물 편석을 고려한 열 및 유동장을 계산하였습니다. 시뮬레이션 결과, G6 대비 G7로에서 실리콘 결정 내 평균 산소 및 탄소 농도가 각각 6.7%, 7.3% 감소하였습니다. 결정이 성장함에 따라 산소 농도는 점진적으로 감소하고 탄소 농도는 증가하는 경향을 확인하였습니다.
3. Methodology
3.1. 기하학적 구조 및 열전달 모델: G6(1124×1124×343 mm³)와 G7(1320×1320×385 mm³) 시스템의 흑연 가열기, 석영 도가니, 단열재 구조를 반영한 글로벌 모델링을 실시하였습니다.
3.2. 불순물 결합 수송 모델: O와 C의 용융액 및 결정 내 확산 계수를 각각 5.0×10⁻⁸ m²/s 및 5.0×10⁻¹¹ m²/s로 설정하고, 5단계 화학 반응(용해, 증발, 운송, 재용해, 편석)을 연동하였습니다.
3.3. 수치적 기법: 아르곤 가스 유동과 실리콘 용액 유동을 Navier-Stokes 방정식을 통해 풀이하고, 복사 열전달(View Factor 등)을 고려한 열 평형 계산을 수행하였습니다.
4. Results and Analysis
4.1. 50% 결정화 단계에서의 불순물 분포: 아르곤 가스 유동 구조 변화에 따라 SiO 가스가 도가니 상단으로 더 효율적으로 배출됨을 가시화하였습니다. G7에서는 좁아진 상부 공간 덕분에 아르곤 유동이 더 강력하게 유지되어 불순물 정체 구역이 감소하였습니다.
4.2. 결정 성장 단계별 농도 변화: 결정화율 20% 지점에서는 용융액 중심부의 탄소 농도가 높았으나, 80% 지점에서는 편석 현상으로 인해 계면 근처의 농도가 최고치에 도달하는 등 농도 분포의 위치 전이(Shift) 현상을 분석하였습니다.
5. Figure and Table List
- Fig. 1: G6 및 G7 방향성 응고로의 구조 및 시뮬레이션 격자도 비교.
- Table 1: 산소 및 탄소 수송 과정에서 발생하는 5가지 단계별 화학 반응식 요약.
- Fig. 3: 대형화 전후의 실리콘 내 산소 및 탄소 농도 분포 비교맵. 대형화 후 수평 방향 농도 편차가 증가하는 양상을 보여줍니다.
- Fig. 4: 용융액 표면 상부의 아르곤 가스 유동 구조. G7에서 와류(Vortex) 구조가 변화하여 불순물 제거 효율이 향상됨을 가시화.
- Fig. 8: 결정화 분율에 따른 평균 산소 및 탄소 농도 변화 그래프. O 농도 감소와 C 농도 증가의 정량적 추이를 표시합니다.
6. Reference List
- Wu et al. (2015) – Insulation partition design using transient global simulations.
- Fornies et al. (2021) – UMG silicon for solar PV: from defects engineering to PV module degradation.
- Liu et al. (2017) – Effect of crucible cover coating on global oxygen and carbon transport.
- Qi et al. (2021) – Effect of Argon flow on oxygen and carbon coupled transport.
- Kumar et al. (2020) – Comparative analysis of thermal stress in 7kg to 330kg silicon ingots.
Technical Q&A
Q: 장비가 대형화되었는데 왜 산소 농도가 오히려 감소하나요?
A: 대형화된 G7 로에서는 도가니 상단 덮개와 용융액 표면 사이의 거리가 좁아지면서 아르곤 가스의 유속이 더 잘 유지됩니다. 이 강력한 아르곤 유동이 용융액 표면에서 증발한 SiO(이산화규소 전구체 기체)를 더 신속하게 배출하고, 외부의 탄소 오염 기체(CO)가 용융액으로 재유입되는 것을 물리적으로 차단하기 때문입니다.
Q: 실리콘 결정 내에서 탄소 농도가 위로 갈수록 높아지는 이유는 무엇입니까?
A: 탄소는 실리콘 결정보다 용융액에 머물려는 성질(분배 계수 k=0.07로 매우 낮음)이 강합니다. 결정이 성장함에 따라 배출된 탄소 원자들이 남은 용융액에 계속 축적되고, 공정 후반부로 갈수록 고농도의 용융액으로부터 결정이 생성되기 때문에 잉곳 상단부의 탄소 농도가 하단부보다 월등히 높아집니다.
Q: 시뮬레이션 결과에서 나타난 ‘수평 방향 불균일성’은 현장에서 어떤 문제를 일으키나요?
A: 잉곳 중심부와 가장자리 사이의 불순물 농도 차이가 크면 태양전지 기판(Wafer)별로 효율 편차가 발생합니다. 특히 중심부의 산소 농도가 안정화되지 않으면 열처리 과정에서 산소 석출물(Oxygen Precipitates)이 불균일하게 형성되어 소자의 수명을 단축시킬 수 있습니다.
Q: G7 공정에서 탄소 함량을 더 줄이려면 어떤 조치가 추가로 필요할까요?
A: 본 논문의 결과를 응용한다면, 아르곤 가스의 유동 패턴을 더 정밀하게 제어하거나, 도가니 덮개(Cover)의 재질을 탄소 발생이 적은 내열 합금이나 코팅재로 교체하여 CO 가스의 발생원 자체를 억제하는 전략이 유효할 것입니다.
Conclusion
본 연구는 실리콘 잉곳 주조 장비의 대형화가 규모의 경제뿐만 아니라 불순물 제거 효율 면에서도 유리할 수 있음을 과학적으로 입증하였습니다. 특히 가스 정체 구역을 최소화하는 유동 설계가 대형 잉곳의 순도를 결정짓는 핵심 변수임을 밝혀냈습니다.
앞으로 실리콘 잉곳 산업은 G8, G10 급 이상의 초대형화로 나아갈 것이며, 이때 발생하는 더욱 극심한 농도 불균일성을 해결하기 위한 3차원 유동 제어 및 자력 가동(Magnetic Field) 지원 기술 등에 대한 추가 연구가 지속되어야 할 것입니다. 본 시뮬레이션 모델은 이러한 차세대 장비 최적화의 견고한 기초가 될 것입니다.
출처 정보 (Source Information)
Citation: Wenjia Su, et al. (2023). Numerical Analysis of the Influence of Size Upgrading on Oxygen and Carbon Impurities in Casting Silicon . Silicon.
DOI/Link: https://doi.org/10.21203/rs.3.rs-2096515/v1
Technical Review Resources for Engineers:
▶ 원문 논문 보기 (PDF)
▶ 기술 검토 및 적용 가능성 문의
This material is for informational purposes only. Unauthorized commercial use is prohibited.
Copyright © 2026 STI C&D. All rights reserved.