열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 ·...

57
-1- 유체 연소 시스템 설계해석 S/W 개발 및 사업화 결과보고서 창업기업명 컴주 : ( ) : 한국산업기술평가원

Transcript of 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 ·...

Page 1: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 1 -

열 유체 연소 시스템 설계해석 S/Wㆍ ㆍ

개발 및 사업화 결과보고서

창업기업명 클 린 컴 주: ( )

대 표 김 용 모:

한국산업기술평가원

Page 2: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 2 -

결 과 보 고 서

년도 신기술창업보육사업에 의하여 완료된 열 유체 연소시스템 설계해석2001 “ ㆍ ㆍ

개발 및 사업화 의 결과 보고서를 별첨과 같이 제출합니다S/W ” .

첨 부 결과보고서 부: 10

년 월 일2002 7

창 업 기 업 명 클린컴 주: ( )

대 표 : 김 용 모 인( )

한국산업기술평가원장 귀하

Page 3: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 3 -

제 출 문

산업기술평가원장 귀하

본 보고서를 열 유체 연소시스템 설계 해석 소프트웨어 개발 및 사업화 의 결” ”ㆍ ㆍ

과 보고서로 제출합니다.

2002 . 7 . .

기 술 개 발 책 임 자 : 김 용 모

참 여 기 술 인 력 : 김 후 중

강 성 모

안 재 현

이 영 주

이 동 화

이 주 원

Page 4: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 4 -

목차 및 내용

기술개발 및 사업화 목적과 중요성I.

열 유체 연소 시스템 설계해석 개발의 경제 사회 기술적 중요성1. . . S/W ㆍ ㆍ

기술적 측면1.1

경제 산업적 측면1.2 ㆍ

사회문화적 측면1.3

앞으로의 전망2.

문제점3.

기술 개발 방안4.

기술개발 내용II.

제품의 개요1.

개발된 의 내용2. S/W

2.1 CAC

2.2 SCALE

성능비교3.

개발시 애로사항 및 해결 과정4.

양산사업화 추진내용III.

추진경과1.

매출실적과 향후 계획2.

사업화 추진상 애로사항 및 해결과정3.

기대효과IV.

첨부자료V.

제품에 관한 기사1.

국내 및 외국관련 업체와의 상호협조 양해각서2.

팜플렛3.

Page 5: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 5 -

기술개발 및 사업화 목적과 중요성I.

열 유체 연소 시스템 설계해석 개발의 경제 사회 기술적 중요성1. . . S/W ㆍ ㆍ

기술적 측면1.1

갈수록 심각해지는 환경오염과 에너지원 고갈로 인해 저공해 고효율 열 유체 연소/ / /

시스템에 대한 개발 수요가 급증하고 있으나 국내의 관련 기술에 대한 연구 개발

실적은 매우 빈약한 실정이다 열 유체 연소 시스템에 대한 설계해석기술은 미국. . . ,

유럽을 비롯한 선진국의 경우 오랜 경험에 의하여 축적된 기술과 전산유체 기반의

해석기술을 이용하여 상당한 수준에 도달해 있으나 국내의 경우 이러한 설계해석기

술을 체계적으로 개발한 적이 없다.

경제 산업적 측면1.2 ㆍ

열 유체 연소 시스템과 관련된 대부분의 원천기술 뿐 만 아니라 핵심 기술을 선ㆍ ㆍ

진국에 의존하고 있는 산업계의 현실을 생각해 볼 때 앞으로의 기업의 생존력을,

높이기 위해서는 이런 원천 기술 및 핵심 기술에 대한 연구 개발이 활발히 이루어

져야 한다 이를 위해 열 유체 연소 해석코드를 통한 최적설계로 열 유체 연소 시스. . . . .

템의 개발 기간을 대폭 단축할 수 있으며 그에 따른 개발 비용의 절감이 가능하며,

제조업등 전통 산업의 고부가가치 산업화를 위한 기술의 접목이 가능하고 막대한IT

기술 도입료의 절감 및 핵심 연소 기술의 배양을 통해 국가 경제 발전에 이바지할

수 있기 때문에 전문 기술의 상용화를 통하여 저공해 산업용 연소기 최적 설계 기

술을 보급하여 국내관련산업체의 기술경쟁력을 제고시킬 수 있다.

사회문화적 측면1.3

지금까지 외국의 설계해석 소프트웨어 개발 업체와 기술 컨설팅 업체에 의존해 왔

던 국내 시장의 상당 부분을 되찾을 수 있을 것으로 예측되며 이를 통해 국내 열유

체연소 시스템 산업에 상당한 파급 효과가 있을 것으로 보이며 단순히 코드 판매만

을 주력하는 외국 소프트웨어 개발사와 핵심 기술에 대한 전수를 기피하는 외국계

기술컨설팅 업체와 달리 본 사는 관련 핵심 기술의 전수와 심화된 기술 컨설팅을

통해 국내 관련업체의 기술 축적이 가능하므로 이에 따른 국내 관련 기술의 개발이

상당히 가속화 될 것으로 예상된다.

Page 6: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 6 -

앞으로의 전망2.

최근 컴퓨터와 전산유체역학 의 기술을 이용(Computational Fluid Dynamics, CFD)

한 유동해석 기술의 발달로 인해 해석 프로그램을 이용한 해석 및 설계가 점차 각

광을 받고 있다 이런 기반의 모의 해석을 이용하면 기존의 실험을 통한 성능. CFD

해석을 통해 소요되는 막대한 경비 인력 및 시간등의 손실을 대폭 절감할 수 있으,

며 측정이 곤란한 부분에 대한 정확한 예측을 통해 해석 및 설계의 정밀성과 효율

을 향상시킬 수 있다 그리고 이러한 전산 유체 역학기반의 설계해석 기법을 이용.

하면 각종 산업용 열 유체 연소 시스템의 개발 기간을 단축시킬 수 있으며 이에 따/ /

라 개발 비용을 절감할 수 있다 또한 막대한 기술 도입료을 절감할 수 있을 뿐만.

아니라 핵심 기술 배양을 통해 국가 경제 발전에 이바지할 수 있으며 이에 따른 경

제적 이익의 창출이 가능하다.

문제점3.

열 유체 연소 시스템의 특성을 해석하기 위한 상용코드의 개발은 수치 및 물리 모델. .

들에 대한 체계적이고 포괄적인 접근을 필요로 하기 때문에 전문인력에 의한 장기

간에 걸친 집중적인 연구개발이 필수적이다 그러므로 경쟁력이 기술을 지속적으로.

개발하기 위해서는 우수한 전문인력의 지속적인 수급 이러한 우수인력에 대한 적,

절한 보수 그리고 이를 위해서 안정적인 시장을 확보해야 하나 아직까지 국내 기,

술력으로 자체 개발된 코드에 대한 불신으로 기존의 외국의 범용 코드에 비해 보다

정확한 계산 결과와 수치적 효율성을 확보했음에도 외국의 범용 코드를 선호하므로

해서 사업초기에 다소 어려움이 예상된다.

기술개발 방안4.

열 유체 연소 시스템 설계해석기술의 개발에 있어서 국내의 한정된 전문인적자원등. .

의 문제점을 해결하기 위하여 핵심연구인력에 대한 교육 정기적인 세미나 국제학, ,

회 참가 및 발표 국제교류 그리고 산업체 대학 연구소와의 유기적인 협조체제의, , , ,

구축과 함께 정보의 공유와 공동 기술개발의 활성화를 통해 국제경쟁력을 가지는

기술을 계속해서 개발하고자 한다 또한 각 분야별로 최고의 권위를 자랑하는 전문.

가를 기술 자문단으로 영입하여 기술 컨설팅에 대한 자문뿐 만 아니라 과제에 직접

참여하여 기술 자문의 질적 향상을 추구할 것이다.

Page 7: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 7 -

기술개발 내용II.

제품의 개요1.

클린컴 주 에서 개발한 주요 는 열 유체 연소 시스템의 해석용 코드 와( ) S/W / / (CAC)

연소기 설계용 프로그램 로 되어 있다Scale-Up (SCALE) .

열 유체 연소 시스템의 해석용 코드인 는 일반적인 열유체해석과 함께 연소과/ / CAC

정 및 공해물질 생성특성 해석을 위해 개발된 전문화된 코드이다 코드는 혼합. CAC

기의 점화특성 화염전파 화염안정성 연소효율 분사노즐의 분무 특성 비평형 화, , , , ,

학반응 복사열전달 난류 유동 난류와 화학반응의 상세 작용 공해물질, , , , (NOx,

등 의 생성 메커니즘 액적의 미립화 및 분열 난류와 액적의 상호 작용 그리Soot ) , , ,

고 액적 증발 그리고 연소 불안정 현상 등의 상세 물리적 과정을 제대로 예측할,

수 있기 때문에 저공해 고성능 연소기 추진 및 동력시스템의 개발 및 운전에 유용,

하게 활용될 수 있으며 또한 비정렬 격자 유동해법을 이용하여 연소기내의 기하학

적 및 물리적으로 복잡한 난류 연소 유동장을 수치효율성을 유지하면서 정확하게

모델할 수 있기 때문에 실제적인 연소시스템의 설계해석을 위한 궁극적인 설계도구

로 활용될 수 있다 특히 버너 산업로 자동차엔진. , , (diesel, gasoline, GDI, HCCI),

가스터빈 엔진 로켓엔진, 그리고 소각로등의 연소과정 및 공해물질 배출특성을 해

석할 수 있고 유전자 알고리즘과 결합하여 최적설계에 활용될 수 있으며 연소과정

을 이용하는 에서 미소입자생성 과정material processing , 연소관련 폭발현상 그,

리고 건물내 화재시의 화염전파 소화 및 재점화 특성등의 해석에도 적용될 수 있,

다.

또한 연소기 설계용 프로그램인 코드는SCALE 난류 모델 액적 모델 열전달 모델, , ,

상세 화학 반응 모델등의 물리 수치 모델을 이용한 수치해석적 설계 기법과ㆍ

설계 기법을 이용하고computational scaling 난류화염구조를 기반으로 한

설계해석기술을 이용하여 소규모 연소실험용 연소기의 실험결과를Scale-Up 100

배 이상 큰 연소기의 연소특성 해석에 효율적으로 적용할 수 있다 산업용 연소 시.

스템을 통한 데이터 베이스의 구축에는 많은 어려움이 따르므로 설계에 필요한 대

부분의 데이터는 연구용 연소기의 연소실험 결과로부터 취득하기 때문에 이러한 화

염구조를 기반으로 한 설계해석기술은 유용하게 활용될 수 있으며 연료 공기의 유,

입 속도등의 기본적인 운전조건 단순한 기하학적 관계 등을 통한 간략한 실험식을,

사용한 기존의 이 해석에 큰 한계를 가지고 있는scaling model 비선형성이 큰 연

소시스템의 설계해석에 적절히 적용될 수 있다.

Page 8: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 8 -

개발된 의 내용2. S/W

열 유체 연소 특성 및 공해물질 생성 해석 전용 코드2.1 CAC ( )ㆍ ㆍ

코드의 열 및 물질전달을 수반한 각종 시스템의 열 유체 연소 특성 및 공해CAC , ,

물질의 생성 과정을 해석할 수 있는 전용 프로그램으로서 각 분야별 최신 모델을

적용하였기 때문에 기존의 범용 공학 해석 프로그램에서 제공할 수 없는 정확도와

수치효율이 보장되는 전용 코드로서 각종 시스템의 설계 최적 운전 조건 도출 성, ,

능 진단 및 개선점 도출에 적용할 수 있으며 코드의 주요 사양은 아래와 같다CAC .

주요 사양■다영역 유동해법을 적용한 유한체적법

정렬 및 비정렬 격자로 구성된 격자계Hybrid

모든 수에 적용 가능한 정상 및 비정상 층류 및 난류 유동장Mach , Solver

복사열전달 해석을 위한 구분종좌법 유한체적법 모델(DOM), (FVM)

자연 대류 및 강제 대류 열전달

이동 경계 수치 절차 및 자유 표면 유동해석을 위한 모델VOF

평형 및 비평형 상세 화학반응 모델

압력과 온도에 따른 가변 열화학 물성치

액체연료의 미립화 액적 분열 및 충돌 모델,

입자와 난류 유동장의 상호작용 해석 모델(Dispersion, Modulation)

아임계 및 초임계 영역에서 적용할 수 있는 유효전도도 액적증발 모델

액체 연료에 포함된 다성분 화학종을 고려할 수 있는 연속 물성치 모델

벽면에서의 연료액막의 영향을 고려한 연료분무의 벽면 충돌 모델

난류 분무 화염장 해석을 위한 모델RIF

를 포함한Thermal, Prompt, Fuel NOx Full NOx Chemistry

및 의 모델Linstedt Franklach Soot

병렬머신용 병렬 알고리즘PC Cluster

난류와 화학반응의 상호작용 해석을 위한 정상 및 비정상 화염편 모델

Page 9: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 9 -

기능■자동차 엔진 성능 해석(Diesel, Gasoline, GDI, HCCl)

버너 유동 및 화염 특성 해석 및 설계

산업 설비 시스템의 열전달 및 유동 해석

가스 터빈 및 로켓 엔진의 성능 및 안정성 해석

소각로 소각 및 공해물질 배출 특성 해석

산업용 연소로 용해로 가열로 열적 특성 해석, ,

연소과정을 이용한 미소입자 생성과정 해석(Nano-particle)

연소 관련 폭발 현상의 수치모사 및 사고 해석

건물내 화재시 화염전파 소화 및 재점화 특성 해석,

성능■

층류 및 난류 부상화염

부상화염은 많은 실용적인 연소기에서 흔히 나타나는 현상이나 기존의 상용코드는

이러한 부분예혼합 화염장을 제대로 해석할 수 없다 코드는 복잡한 형상의 연. CAC

소유동장에 대해 격자를 생성할 때 격자구조 즉 격자점을 서로 연결하는 기하학적

위상을 제거하여 임의의 형태로 격자를 만듦으로써 격자 생성이 유동장의 복잡성으

로부터 매우 자유롭고 국부적인 화염영역으로 격자 적응성이 용이한 비정렬 격자계

와 검사체적을 구성하기 위해 다면체 혼합격자의 사용 및 차원으로의 확장성과 경3

계조건 부여가 용이한 셀 중심법을 사용하여 물리적으로 복잡한 화염장 해석을 적-

용성 현저하게 개선하였다.

Page 10: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 10 -

중심부의 연료 제트 유동과 동축류 희박 연료 제트 유동의 속도 분포에 대해 완전

발달한 층류 파이프 유동 과 환형 유동 을 각각 가정하였(tube flow) (annular flow)

다 아래의 수치해석결과에서 보는 바와 같이 코드는. CAC 축방향 속도의 반경방향

분포 그리고 복잡한 화연구조를 가지는 삼지점부근 축방향 속도 및 온도분포을 잘,

예측하고 있음을 알 수 있다.

또한 아래 해석결과에서 보듯이 코드에 포함되어 있는 난류CAC strained flamelet

연소모델은 난류 부상 화염의 상세화염구조와 연소특성을 잘 예측하고 있음을 알

수 있다 이러한 난류부상화염은 공해물질 배출특성 및 화염안정화에 지대한 영향.

을 미치기 때문에 연소시스템 설계해석시 반드시 고려해야 할 필요가 있으나 기존

상용코드는 이러한 현상을 예측할 수 없다.

그림 층류 부상 화염1

Page 11: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 11 -

그림 난류 부상 화염2

비예혼합 난류제트화염내의 화염구조 몇 생성 예측CO/H2/N2 NO

난류연소 모델의 예측능력을 검증하기 위해 서로 다른 연료노즐지름, (4.58mm.

을 가지는 두가지 제트화염을 해석하였다 이때 두 화염의 레이놀즈수는7.72mm) . ,

으로 서로 같기 때문에 노즐지름으로 스케일링을 할 경우 온도 및 주요화학16,700 ,

종들은 거의 완벽하게 상사관계를 갖게 된다 그러나 전체유동변형율 및 체류시간. ,

이 다르게 되므로 난류혼합에 의한 비평형화학반응의 정도가 차이가 나게 되고 이,

에 민감한 라디칼 및 분포는 서로 상이한 결과를 나타내게 된다 두 화염간에NO .

존재하는 이러한 상반된 경향 즉 온도 및 주요화학종의 상사관계와 라디칼 및, NO

의 민감성은 난류연소모델의 예측능력을 검증하는 데 매우 유용하다. CO 40%, H2

그리고 로 구성된 연료가 의 속도를 가지는 공기제트로 분사30%, N2 30% 0.7 m/s

되며 이때 이론당량 혼합분율 은 이다 이론당량선을 따라 영역을 효과, , Zst 0.295 .

적으로 분할하기 위해서 축대칭 비정렬 삼각격자계를 사용하였다.

Page 12: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 12 -

연료노즐근방에서는 급격한 난류혼합에 의해서 스칼라소산율이 매우 크게 존재하며

하류로 가면서 난류에 의한 혼합이 완화됨에 스칼라소산율 역시 매우 급격하게 감

소하게 되는 일반적인 제트화염에서의 특징을 나타내고 있다 그러나 와. Flame A

비교할 때 는 낮은 속도와 큰 노즐지름을 갖고 있기 때문에 전체유동변형, Flame B

율이 훨씬 작게 된다 따라서 전체적으로 에 비해 훨씬 낮은 스칼라소산율. , Flame A

의 분포를 가지게 됨을 알 수 있다 반대로 속도가 낮기 때문에 체류시간은. Flame

가 훨씬 크게 나타남을 알 수 있다 아래결과에서 보는 바와 같이 코드에서B . CAC

사용된 비정상 화염편 난류연소 모델은 중심축을 따라 주요 화학종 라디칼 그리고, ,

의 농도분포등의 상세화염구조를 잘 예측함을 알 수 있으며 기존의 상용코드는NOx

이러한 예측능력을 전혀 보유하지 않은 상태에 있다 복사와 를. full NOx chemistry

고려한 비정상 층류화염편 모델은 실험에서 나타난 에 대한fluid-dynamic scaling

라디칼과 의 민감성을 잘 예측하고 있음을 알 수 있다OH NO .

Page 13: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 13 -

그림 3 Turbulent CO/H2/N2 Jet Flames (Barlow et.al.,2000)

Page 14: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 14 -

그림 4 Centerline Profiles : Mean Flames Structures Lagrangian Flamelet Model

(Red line :with radiation; Blue line :adiabatic)

그림 5 Sensitivity Of Minor Species To Fluid-Dynamic Scaling

Page 15: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 15 -

아세틸렌 공기 난류 확산화염에서의 분포 특성/ soot

아세틸렌 공기 비 예혼합 제트화염 실험에서 계측된 체적분율을 예측하기 위/ soot

하여 수치계산을 수행하였다 질소로 희석된 연료. (50% C2H2 + 50% N2 가 지름)

인 연료노즐을 통해 의 체적 유량을 가지고 유입되며 이 버너는 내부2mm 3.2 l/min

지름이 인 원통형 굴뚝에 의해 둘러싸여져 있으며 산화제125mm (80% N2 + 20%

O2 는 연료노즐과 굴뚝사이에서 의 체적유량을 가지고 분출된다 열역학적) 10 l/min .

상태량 및 각 화학종의 농도를 계산하기 위하여 을 이용하laminar flamelet model

였으며 주요 화학종 및 라디칼의 분포특성을 살펴보기 위하여 개의 화학종과101

개의 반응단계를 가지는 상세화학 반응모델을 이용하였다 의 분포특성을544 . soot

살펴보기 위하여 의 체적분율과 수 밀도에 대한 수송방정식을 계산하였으며 이soot

방정식의 생성항은 등에 의해 제안된 모델과 등에 의해 제안된Lindstedt Frenklach

모델을 이용하였다 상대적으로 간단한 모델에 비하여 코드에 포함. Lindstedt CAC

되어 있는 모델을 이용한 경우 의 분포와 체적분율의 분Frenklach soot precursor

포가 매우 홉사하기 때문에 의 분포 특성이 좀 더 상류쪽으로 분포함으로써pyrene

좀 더 개선된 분포특성을 가짐을 알 수 있다soot .

Page 16: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 16 -

그림 6 Temperature & Acetylen Mass

Fraction 그림 7 Pyrene Mass Fraction

그림 8 Lindstedt Model 그림 9 Franklach Model

Page 17: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 17 -

그림 10 Axial Profile Of Soot Volume Fraction

연소실내의 선회난류 유동장 및 공기 다단연소버너내의 화염구조 해석

과 의 선회난류 유동에 관한 예측능력을 검정하였다 선회기를 격Roback Johnson .

자생성 프로그램을 이용하여 모델링하고 개의 선회 베인중 한 개가 들어가도록8

각도로 계산 영역을 선정한 차원 격자계를 나타낸다 수치해석 결과에 나타45 3 .℃

낸 바와 같이 선회난류유동장의 예측된 속도장 난류강도 및 난류 응력 분포를 실,

험치와 비교해 본 결과 난류모델에 비하여 코드에서 채택하고standard k-e CAC

있는 난류모델이 실험치와 상대적으로 잘 일치함을 알 수 있었다extended k-e .

또한 코드는 대체적으로 실험치에 근접하게 예측하고 있으며 입구조건이 불확CAC

실한 경우 삼차원 모델을 효율적으로 적용할 수 있음을 의미한다.

또한 코드를 사용하여 선회기가 장착된 공기 다단 연소기의 연소특성를 하였CAC

다 실험에서 사용한 를 의 프로판과 의 메탄의 혼합가스로 가정하였으. LPG 95% 5%

며 생성메카니즘을 포함한 상세한 연소메카니즘을 고려하기 위해서 개의, NOx 55

화학종간의 개의 화학반응식을 사용하였다 난류에 의한 반응면내의 확산속도를266 .

지배하는 스칼라소산율에 의해 나타나는 화염구조의 영향을 잘 반영하고 있으며,

스칼라소산율이 xq 보다 커지는 경우 반응면으로부터의 열전달량과 반응면으=53/s

로 확산되어 들어가는 물질전달량이 커져서 더 이상 연소반응이 유지되지 못하고

국소소염이 나타나게 됨을 알 수 있다 수치해석 결과는 다단연소기의 선회화염구.

조를 정성적으로 잘 예측하고 있음을 알 수 있었다.

Page 18: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 18 -

그림 11 Swirler Solid Modeling 그림 난류 선회 연소기12

그림 난류 선회 연소기 유동 해석13

Page 19: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 19 -

그림 14 Air Staged Burner

차원 연소기 및 산업용 연소로내의 난류화염장 해석3 Syngas

를 통해 화염안정화를 시키는 차원 연소기 내에서 일어나는 연소과Flameholder 3 [6]

정을 정상 층류화염편 모델을 이용하여 해석하였으며 연료는 CO/H2/N2

의 혼합기를 사용하였다 기하학적 및 물리적으로 복잡한 삼차원 난류(40/40/20%) .

연소장을 상세화학반응을 고려하여 로 계산하는 데 걸리는 시간은 지나single CPU

치게 과도하기 때문에 병렬화된 코드를 이 문제를 해석하였다 아래 그림에서CAC .

보듯이 계산영역을 개로 나누어 병렬계산을 수행할 에 비해 계산시간8 single CPU

을 정도로 단축할 수 있었다 이는 병렬 코드를 이용할 경우 복잡한 실제의1/6 . CAC

차원 연소시스템의 최적설계해석이 가능함을 의미한다 안쪽으로 들어오는 연료는3 .

바깥공기제트와 혼합을 이루면서 연소반응이 일어나게 되고 에 의한 재flameholder

순환유동을 통해 화염안정화를 이루고 있다 라디칼의 분포로 알 수 있듯이 반. OH

응영역은 주로 후류의 재순환영역 바깥쪽 혼합층을 따라 이루어지고flameholder ,

부력의 영향에 의해 고온의 화염이 윗방향으로 크게 치우치는 경향이 뚜렷이 나타

나고 있음을 알 수 있으며 코드의 수치해석결과는 확산화염의 구조와CAC Syngas

의 생성특성에 있어서 실험치와 정량적으로도 대체적으로 잘 일치하고 있음을NO

알 수 있다 또한 매우 복잡한 산업용 연소로에도 비슷한 수치절차를 적용하여 효.

율적으로 계산할 수 있었다.

Page 20: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 20 -

그림 병렬 계산15 Syngas Combustor ( )

Page 21: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 21 -

그림 산업용 연소로16

그림 산업용 연소로 열유동 해석17

Page 22: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 22 -

가스터빈 연소기 및 내부 난류유동장Scroll

본 연구에서는 에서 와 가 수행한 물을 이용한 모Imperial College McGuirk Palma

델 연소기 실험을 풀어보았다 개의 주 연소 공기공과 희석공기공이 각각 한 개씩. 6

들어갈 수 있도록 로 잘랐고 입구측 에서부터 희석 공기공까지는 사면체격자60 ,˚

계를 이용하였고, 프리즘격자를 이용하여 계산을 수행하였다 이러한 비정형격자.

삼차원 모델은 가스터빈의 선회난류 유동특성을 대체로 잘 예측하고 있Can-type

음을 알 수 있다 또한 코드는 가스터빈 엔진에서 내부유동장 및 온도분. CAC Scroll

포의 기본특성을 잘 예측함을 알 수 있다.

그림 유동 특성 해석18 Gas Combustor

Page 23: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 23 -

그림 내부 유동장 해석19 Gas Turbine Engine Scroll

직접분사식 디젤엔진의 연소과정 및 공해물질 배출특성 해석

산업용 연소시스템에 있어서 액체연료의 분무연소특성은 액적의 미립화 액적과 유,

동장의 상호작용 액적의 운동과 증발과정 및 연료증기의 자발화 등에 의해 지대한,

영향을 받게 된다 최근들어 액체연료를 사용하는 고성능 엔진에서의 연소실 압력.

이 점점 높아져 연료의 임계압력에 가깝거나 또는 넘어서 운전 되고 있기 때문에,

저압상태의 증발과정에 비해 현저하게 다른 증발특성을 가지는 고압상태에서의 연

료 액적의 증발과정을 잘 모사할 수 있게 하기위해서 코드에 고압증발과정을CAC

모델링하였다 또한 기존의 사용모드는 을 사용하고 있지만 이. shell ignition model

모델은 문제의 조건에 따라 계수를 조정해야 하고 난류와 화학반응의 상호작용을

전혀 고려할 수 없는 결점을 가지고 있다 이러한 문제점을 해결하기 위하여 유동.

장과 화염장을 분리해서 계산하는 모델Representative Interactive Flamelet (RlF)

을 코드에 모델링하였으며 이 모델은 계수조정이 필요없을 뿐만아니라 상세CAC RlF

화학반응을 적용할 수 있는 이점을 가지고 있다.

디젤연소의 고압환경에서 분무 유동과 증발 연료분무의 자발화 및 연소를 수치적,

으로 분석하기위해서 액적의 가열 액적 내부의 유동 강제대류 실재, , , Stefan flow,

기체 효과 주위공기의 용해 등 고압에서 일어나는 효과들을 모두 고려한 고압증발,

모델을 이용하였다 또한 모델을 이용하여 유동장과 화염장을 분리해서 계산하. RlF

였고 스칼라 소산율의 공간적 비균일성을 고려하기 위하여 모델multiple RlF [11]

을 이용하였다 여기에 사용한 연료는 으로 자발화를 예측할 수 있는 저. n-heptane ,

온에서의 탄소고리의 끊김 고온에서의 연료의 산화 등을 고려한 화학종, 44 , 114

화학반응을 이용하였다 아래에 도시된 수치해석결과에서 보는 바와 같이 코드. CAC

는 직접분사식 디젤엔진의 분무연소과정 및 공해배출 특성을 대체로 잘 예측하였

다.

Page 24: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 24 -

그림 20 combustion processes & pollutant formation in DI Diesel Engine (a)

Mean Temperature (b) NO Concentration (c) OH Radical Mass Fraction (d)

Soot Volume Fraction

Page 25: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 25 -

기타 열 유체분야 적용사례 자동차 자동차 로켓등, ( , , )

클린컴 주 에서 개발한 코드의 전문성 과 범용성에 있어서 경쟁력을 갖추고 있( ) CAC

음을 설명하기 위하여 복잡한 기하학적 형상을 가지는 자동차 및 항공기의 외부 난

류유동장 그리고 물리적으로 복잡한 로켓내부의 연소유동장을 포함하였다 코, . CAC

드는 이러한 복잡한 문제들을 높은 수치효율성과 공학적인 정확도내에서 성공적으

로 해석하였다.

그림 자동차 외부 유동 해석21

Page 26: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 26 -

그림 자동차 외부 유동 해석22

그림 수치해석용 격자계 구성23 Boeing-747

Page 27: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 27 -

그림 외부 유동장 해석24 Boeng-747

그림 25 Vortex Hybrid Rocket Engine Mesh

Page 28: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 28 -

그림 26 Velocity Fields & Temperature Distribution

그림 27 Profiles Of Heat Flux At The Solid Propellant Surface

Page 29: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 29 -

연소 시스템 설계해석용 프로그램2.2 SCALE( )

코드는 수치해석적 설계 기법과 설계 기법을 접목하SCALE Computational Scaling

고 난류화염 구조에 기반을 둔 연소 시스템 설계 해석용 코드로서 대형 산업용 연,

소 시스템의 설계에 필요한 각종 데이터 및 향후 운전 성능에 대한 예측 결과를 제

공하여 시스템 설계 비용과 시간의 절약이 가능한 연소 시스템 설계 해석 전용 프

로그램이며 주요사양은 아래와 같다.

주요사양■기하학적 공기역학적 열적 비상사, , scaling model

기하학적 공기역학적 열적 상사, , scaling model

화염 구조 기반 기법scaling

연소 시스템 설계scale-up

연소시스템 설계해석용 프로그램

Page 30: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 30 -

기능■연소 시스템 설계scale-up

Air staging Burner scale-up

Fuel staging Burner scale-up

연소 시스템에서 발생하는 공해 배출물질 배출 특성 예측

성능■

코드의 예측 능력을 검증하기 위해 버너의 스케일을 에서 로SCALE 30KW 12MW

변화시켜가며 최대부하시의 계측치와 예측치를 비교하였다 용량에서EINOx . 30KW

약간의 정량적 불일치를 보이고 있지만 정성적으로 잘 일치하는 결과를 보이고 있

다 이는 에 필요한 각종 형상 및 파라메터들이. scaling constant velocity principle

에 의거하여 제시되었으나 유속이 작은 경우 유동의 스케일과 유동에 대한 상사성

의 이탈로 인해 정량적으로 일치하지 못한 결과를 보이고 있는 것으로 사료된다.

그러나 용량이 증가할수록 상사성이 향상되어 좀더 정량적으로 일치하는 결과를 보

이고 있다 이러한 해석결과는 실험실용 모델연소기의 상세계측결과를 이용하여.

코드로 해석할 경우 배 이상으로 큰 연소기의 설계해석이 가SCALE 100 Scale-Up

능함을 의미 하며 이러한 설계해석기술은 관련산업의 경쟁력과 부가가치를 한층 제

고시킬 것으로 사료된다.

Page 31: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 31 -

성능 비교3.

당사 코드와 상용 전산해석 코드의 전문성 범용성 수치효율성 정확성 결과유용, , , ,

성 설계신뢰성 가격대 성능비 사후지원성 사용편의성의 총 개 항목에 대한 평, , , , 9

가 비교를 수행하였으며 비교 결과는 그림 과 같다 비교 평가 결과 전문성 수치. ,

효율성 가격대비 성능비 등이 매우 우수하게 나타났으며 범용성에 있어서는 상용,

코드에 비해 낮게 나타났는데 이는 전문성을 강조한 전문 해석 코드이기 때문이다.

Page 32: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 32 -

성능 비교 :

비교대상 대상 코드중 각 부문별 최고 성능을 나타내는 코드를 으로 하여100

코드의 성능을 수치화 하였으며 타사 코드의 경우 평균치를 나타냄.

참고 일본 자동차 기술회 유체기술부문위원회 성능 비교 결과( : ,票山利彦

Comparison of CFD Commercial Codes , 1997)

비교 대상 코드 :

당사 코드, FLUENT, CFD-ACE, FIDAP, ICEPAK, STAR-CD, PHOENICS,

총 개 코드CFD-FASTRAN, FIRE 9

개발시 애로사항 및 해결과정4.

열 유체 연소 시스템의 특성을 해석하기 위한 상용코드의 개발은 수치 및 물리 모델. .

들에 대한 체계적이고 포괄적인 접근을 필요로 하기 때문에 전문인력에 의한 장기

간에 걸친 집중적인 연구개발을 필요로 하는 지난한 작업이다 본 클린컴 주 사는. ( )

이러한 애로점을 해결하기 위하여 핵심연구인력에 대한 교육 정기적인 세미나 국, ,

제학회 참가 및 발표 국제교류 그리고 참여 연구자들의 헌신적인 노력 및 강한 사, ,

명감에 의해 이러한 경쟁력있는 기술을 개발할 수 있었다 또한 산업체 대학 연구. , ,

소와의 유기적인 협조체제의 구축과 함께 정보의 공유와 공동 기술개발의 활성화를

통해 핵심 기술에 대한 개발을 성공적으로 진행하였다.

양산사업화 추진 내용III

추진 경과1.

관련 기술 개발 완료■코드 개발- Flamelet

설계 프로그램 개발- SCALE

이용 병렬처리 프로그램 개발- PC cluster

비정렬 하이브리드 격자 프로그램 개발-

개발- Graphic User Interface

코드와 타 코드와의 개발- Flamelet interface

이용 최적 설계 기술 개발- PC Cluster

개선된 수치 및 물리모델 개발-

유전자 알고리즘과 결합된 코드를 이용한 고성능 병렬계산을 통한 최적 설- CFD

계 기술 확립

예측능력 향상을 위한 수치 및 물리 모델 연구 적용- ,

Page 33: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 33 -

관련 첨단분야 에 적용할 수 있는 상세 모델 개발- (IT, BT, NT, ET, ST)

사용자 주문형 코드 개발-

최적 설계 및 운전 조건 도출을 위해 최첨단 기술인 가상현실 기법을 이용한- IT

개발Virtual Engine

코드 상용화■코드 상용화를 위해 코드의 핵심 서브루틴의 개발을 완료하였으며 최신 수치 및 물

리모델을 적용하고 이에 대한 벤치마킹을 완료하였으며 사용자의 편의성을 강조한

를 구축User Interface

기술 컨설팅 체제 구축■수년간 대학 연구실에서 열 유체 연소 관련 기술 개발을 진행해 왔으며 그동안ㆍ ㆍ

상당히 높은 수준의 기술 축적과 고급 기술 인력을 배양하였고 그동안 축적된 다양

한 분야에 대한 연구와 과제 수행 경험을 바탕으로 기술 컨설팅 체제를 구축

각계 전문가로 구성된 기술 자문단 구성-

각 분야별로 최고의 권위를 자랑하는 전문가를 기술 자문단으로 영입하여 기술 컨

설팅에 대한 자문 뿐 만 아니라 과제에 직접 참여하여 기술 자문의 질적 향상을 추

Page 34: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 34 -

성 명 소 속 기술 자문 분야

김세원 박사 한국생산기술연구원 소각로 산업용 연소로,

안국영 박사 한국기계연구원 화력 열병합 발전설비,

손정락 박사 서울대 항공 가스터빈,

최윤호 교수 아주대 전산유체역학(CFD)

최도철 교수 미국UTRC( ) 발전설비 가스터빈,

김종수 박사 KIST 화재 및 사고해석

이기형 교수 한양대 내연기관 자동차 관련,

안강호 교수 한양대 나노 입자 합성, MEMS

이성택 박사 주 한화( ) 로켓 미사일 추진체,

석 박사급 기술 컨설팅 인력 확보- ㆍ

전산유체역학에 대한 전문적 능력과 관련 과제 수행 경험과 기술적 지식이 풍부한

석 박사급 기술 컨설팅 인력을 확보하였으며 이를 바탕으로 다양한 분야에 대한,ㆍ

기술 컨설팅을 제공

기술 컨설팅 업무를 위한 초고속 병렬처리 컴퓨터 보유-

전산유체역학 기반의 기술 컨설팅과 클라이언트의 기술 개발에 필요한 연구개발 비

용을 최소화하고 보다 효율적인 기술 컨설팅을 위해 초고속 병렬처리 컴퓨터 확보

여러 대학 연구실 연구소 기업과 관련 핵심 기술 정보 공유 및 공동 연구 개발- , ,

폭넓은 분야에 대한 핵심 기술의 공유를 위해 관련 대학 연구실 연구소 기업과의, ,

연계를 진행하였으며 관련 기술 세미나를 개최하여 기술 및 주요 정보를 공유하고

공동 관심 분야에 대한 공동 연구개발을 통해 기술 축적도의 증진을 이룩

홍보 및 마케팅■

- Database Marketing

본 사업의 특성상 일반 고객을 대상으로 하는 것이 아닌 열 유체 연소 시스템 관ㆍ ㆍ

련 기업이나 관련 기술 인력이 주요 대상 고객이므로 대상 고객에 대한 고객

를 구축하고 이를 통한 전략을 수립하였고 구축된 를 통해Database Marketing DB

적극적 홍보를 진행해왔으며 상당수 기업 및 인력에 대한 홍보를 완료하였다 향후.

지속적인 홍보를 통해 기술성과 제품 우수성을 홍보.

Page 35: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 35 -

홈페이지 구축-

홈페이지를 통한 인터넷 홍보를 진행하고 있으며 홈페이지의 컨텐츠를 지속적으로

확충 보강 하여 열 유체 연소 시스템 관련 포털 사이트로 자리매김할 수 있도록, ㆍ ㆍ

다양한 정보와 기술 자료를 제공.

홈페이지: www.cleancombustion.co.kr

학술대회 홍보데스크 설치-

한국폐기물 학회 춘계 학술대회에 홍보 데스크를 설치하여 참석자들에게 관련 기술

및 제품 홍보를 진행하였으며 최근 기술에 대한 저변 및 관심의 확대로 인해, CFD

많은 수의 참석자가 관심을 보였으며 기술 컨설팅 의사를 타진.

해외 유명 학술대회 참가-

국내외 각종 학술대회에 참가하여 본 사의 코드를 이용한 해석 결과를 발표하였으

며 이를 통해 코드의 우수성을 홍보

Page 36: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 36 -

참여 학술대회:

미국자동차공학회 년 월 미국 디SAE(Society of Automotive Engineering, , 2002 2 ,

트로이트)

난류비예혼합화염 국제 워크샵 년 월 일본 삿포로TNF workshop( , 2002 7 )

국제연소학회 년 월 일본 삿포로Int. Symposium on Combustion( , 2002 7 )

학술대회 년 월 캐나다 밴쿠버ASME (2002 8 , )

프로그램 설명회 개최-

열 유체 연소 시스템 관련 기업의 임원 및 주요 인력을 초청하여 기술 세미나를ㆍ ㆍ

개최하여 본 사의 제품의 적용성과 정확성을 널리 알렸으며 참석자 대부분이, CFD

기반 설계해석 기술에 대한 가능성과 잠재력에 대해 인식의 폭을 확대시킬 수 있었

고 이를 통해 잠재 고객에 대한 홍보 효과를 얻음

참여업체 및 연구소 주 수국 귀뚜라미보일러 캠코엔지니어링 한국폐기물기술: ( ) , , , ,

벽산엔지니어링 한국기계연구원 두원공조 지역난방공사 대경기계 서울대, , KIST, , , , ,

한양대 건설 동부건설 등, LG ,

프로그램 설명회(2002. 5. 28)

Page 37: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 37 -

기술 세미나 개최-

각 분야의 전문가를 초청 기술 세미나를 개최하여 코드 개발 및 기술 컨설팅에 필,

요한 핵심기술의 공유와 각종 정보를 교환하고 열 유체 연소 시스템과 관련한 최ㆍ ㆍ

신 기술의 개발을 위한 체제 구축

Prof. Bilger(U. of Sidney)

: Combustion Modeling

신병록교수 도호쿠대학 기계공학과( , )

: CFD Modeling for Cavitation Flows

김종수박사(KIST)

화염안정화 메커니즘 해석:

한동희 박사 현대자동차 기술연구소( )

레이저 계측기술을 이용한 난류연소화염장 해석:

박외철 교수 부경대 산업안전공학과( )

화재사고 관련 연소 유동장 해석:

관련 내용 기사화-

월간지 보일러 에 열 유체 연소 시스템 설계해석기술 설명회 개최와 관련된 내용“ ” . .

과 당사에 대한 소개에 대한 기사가 게재 첨부자료 참조( )

관련 기업과의 전략적 제휴■관련 분야의 기업과의 상호업무협조 상호공동연구 협약 체결등을 통한 전략적 제,

휴를 통해 업무 영역의 확대와 원할한 연구 수행

세계적으로 기술력과 권위를 인정받고 있는 미국 사와 국제 공동 연구 협약 체ESI

결을 통해 국제무대로의 진출 가능성 확보

주 경원테크와 상호업무협조계약 체결- ( )

주 명진에어테크와 상호 업무 협조 계약 체결- ( )

사 미국 와 체결- ESI ( ) Agreement of International Research Cooperation

연구개발 인력 및 공간 확보■박사급 연구 개발 인력 확보-

전산유체역학에 대한 전문적 능력과 관련 과제 수행 경험과 기술적 지식이 풍부한

석 박사급 기술 컨설팅 인력을 확보하였으며 이를 바탕으로 다양한 분야에 대한,ㆍ

기술 컨설팅을 제공

Page 38: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 38 -

부설 연구소 공간 확보 한양종합기술원- ( )

독립적인 연구공간의 확보를 위해 한양종합기술원에 연구소 공간을 확보하여 원할

한 연구 수행

연구용 초고속 병렬처리 구비- PC cluster

국내에 극히 제한적으로 보급 활용되고 있는 초고속 병렬처리 를 구비, PC cluster

하여 코드 개발을 위해 사용하고 있으며 상당히 높은 수준의 병렬처리 기술을 이,

미 확보하였고 그동안 축적된 기술력으로 인해 국내 연구소 대학 등에서 기술 컨,

설팅 및 병렬계산 코드의 판매를 타진 해옴

매출 실적 및 향후 계획2.

코드 판매- CAC

본 사에서 개발한 열 유체 연소 특성 및 공해물질 생성 해석 전용 코드인 CACㆍ ㆍ

코드를 판매 서울대학교 기계공학부( )

기술자문 용역 수주-

자동차 엔진 실린더 내의 유동 및 연소특성의 해석과 관련한 기술 자문 용역 수주

Page 39: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 39 -

기술 자문 용역 수주 진행-

년 월 현재 현대자동차 한국전력 기계연구원 생산기술연구원의 개 업체2002 6 , , , 4

및 연구소에서 열 유체 연소 시스템과 관련한 기술 자문을 요청ㆍ ㆍ

연도구분

년2002월현재(6 )

년2002 년2003 년2004

매출액억원( )

만원2,360 억원3 억원12 억원35

주요제품

CAC™SCALE™

CAC™FLAME™SCALE™

CAC™FLAME™SCALE™

범용 CAC™CAC™FLAME™SCALE™

주요납품처

대학연구실서울대(

기계공학부)산업체(KWtech)

출연기관연구소 개(5 )

대학연구실 개(10 )산업체 개(5 )컨설팅 건이상(8 )

출연기관연구소 개(10 )대학연구실 개(15 )산업체 개(10 )컨설팅 건이상(15 )

출연기관연구소 개(10 )대학연구실 개(25 )산업체 개(15 )컨설팅 건이상(30 )

사업화 추진상 애로사항 및 해결과정3.

사업화 추진상 애로사항-

사업의 특성상 보유 기술력의 신뢰도가 사업화에 매우 큰 영향을 미침

그러나 아직까지 열 유체 연소 시스템과 관련한 국내 기술력에 대한 불신이 팽배ㆍ ㆍ

하여 사업 초기 이로 인한 어려움이 상당히 큼.

또한 국내 기술력으로 자체 개발된 코드에 대한 불신으로 기존의 외국의 범용 코드

에 비해 상당히 정확한 계산 결과와 수치적 효율성을 확보했음에도 외국의 범용 코

드를 선호

해결 과정-

이런 국내 기술력과 코드에 대한 불신의 한계를 극복하고 기술 신뢰성을 높이기 위

해 각종 학술대회 및 기술 세미나를 통해 본 사에서 확보하고 있는 기술력이 구미

선진국의 기술력에 비해 결코 뒤지지 않음을 입증하고 본 사에서 개발한 코드의 계

산 정확성과 수치적 효율성이 매우 월등함을 강조하여 국내 관련 업체 대학 연구, ,

소의 관련 핵심 기술 인력에게 상당히 높은 평가를 받음

또한 외국의 학술대회에 본 사의 코드를 이용한 계산 결과를 발표하여 상당수의 참

가자로부터 호평과 관심을 받음.

향후 지속적으로 기술력과 코드의 우수성에 대한 홍보를 통해 국내 기술력이 외국

의 기술력에 뒤지지 않음을 피력하고 이를 통해 코드의 판매와 기술 컨설팅의 증진

을 꾀함

Page 40: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 40 -

기대효과IV

현재까지의 사업화 추진을 통해 상당한 수준의 기술 인지도와 코드 우수성을 검-

증 홍보하였으며 이를 바탕으로 향후 코드 판매와 기술컨설팅 관련 사업을 주력,

사업분야로서 진행

지금까지 외국의 설계해석 소프트웨어 개발 업체와 기술 컨설팅 업체에 의존해-

왔던 국내 시장의 상당 부분을 되찾을 수 있을 것으로 예측되며 이를 통해 국내 열

유체연소 시스템 산업에 상당한 파급 효과가 있을 것으로 사료됨

단순히 코드 판매만을 주력하는 외국 소프트웨어 개발사와 핵심 기술에 대한 전-

수를 기피하는 외국계 기술컨설팅 업체와 달리 본 사는 관련 핵심 기술의 전수와

심화된 기술 컨설팅을 통해 국내 관련업체의 기술 축적이 가능하므로 이에 따른 국

내 관련 기술의 개발이 상당히 가속화 될 것으로 예상됨

그동안 국내의 코드 개발과 기술 컨설팅 업체가 없어 외국 업체 의존해 왔던 국-

내 시장의 진출로 인해 상당히 큰 수입 대체 효과가 기대되며 사 사와, ESI , CFDRC

같이 전세계 시장에서 상당한 영향력을 행사하는 업체와의 전략적 제휴를 통해 해

외 시장으로의 진출이 가능하므로 국내에서는 전무했던 수출 효과가 발생할 수 있

어 국가 경제에 크게 기여할 수 있을 것으로 사료됨.

본 사에서 개발한 코드와 이를 통한 기술 컨설팅을 통해 국내 관련 업체의 기술-

력이 향상되어 이를 바탕으로 한 생산성 향상과 품질 개선이 가능할 뿐 만 아니라

전통 산업의 고부가가치 산업으로의 전환이 가능하므로 국내 산업 발전에 상당한

기여가 가능할 것으로 사료됨.

국내 시장 규모 및 예상 매출 분석◆

Page 41: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 41 -

국내 해석 코드 시장 규모 년 현재 약 백억원으로 추정되었으며 년에: 2001 2 , 2002

는 성장 년에는 성장 년 이후는 성장할 것으로 추정12% , 2003 13% , 20004 12.2%

되며 열 유체 연소 관련 해석코드 시장 규모는 년 현재 총 해석코드 시장의/ / 2001

이며 향후 연평균 의 성장이 예상됨 한국소프트웨어진흥원 분류별34% 4.5% . ( S/W

시장 분석 자료 근거)

국내 기술 컨설팅 시장 규모 년 시장규모는 년 대비 성장한 약: 2001 2000 9.6%

억원으로 추정되며 향후 년 동안 연평균 의 성장세로 증가할 것으로 추680 , 5 12.7%

정되며 열 유체 연소관련 기술 컨설팅 시장 규모는 년 현재 총시장 규모의/ / 2000

로 추정 됨 년 동부증권 예측 자료 한국산업기술진흥협회 시장분석 자14% (2000 ,

료)

향후 예상 영업 이익 분석◆

매출원가는 매출액의 로 추정되며 판매비는 마케팅 홍보등에 소요된 금액으로20% ,

매출액의 로 가정함10%

인건비성 관리비는 관리부문의 인건비 복리후생비 교육훈련비 등 인력규모가, ,

가 되는 비용임 매년 임금 상승율 가정 과 인력 규모 증가 년Cost Driver . (10% ) (4

간 매년 씩 증가 를 반영하여 산정함20% ) .

기타 관리비는 세금과 공과 지급임차료 감가상각비 등으로 설비투자 등 기타 변수, ,

에 의해 결정되는 비용임 유사업체의 실제 발생된 비용을 기준으로 연간 규모를.

결정하고 연간 매출액 증가율 세금과공과 감가상각비 및 연간 증가율 지급임차, ( , ) (

료 을 반영하여 산정함) .

Page 42: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 42 -

첨부자료V.

제품에 대한 기사1.

월간 보일러 월호- 6

국내 및 외국관련 업체와의 상호협조 양해각서2.

경원테크 주- ( )

명진에어테크-

- Engineering Science lnc.

팜프렛3.

Page 43: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 43 -

월간 보일러

Page 44: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 44 -

Page 45: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 45 -

상호 업무 협조 양해 각서

Page 46: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 46 -

상호 업무 협조 양해 각서

Page 47: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 47 -

Agreement of International Research Cooperation

Page 48: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 48 -

Solution for Thermal Systems

Page 49: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 49 -

CFD Consulting Business

Page 50: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 50 -

CFD Consulting Business

Page 51: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 51 -

CFD code Sale Business

Page 52: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 52 -

열 유체 연소 특성 및 공해물질 생성 해석 전용 코드CAC ㆍ ㆍ

Page 53: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 53 -

열 유체 연소 특성 및 공해물질 생성 해석 전용 코드CAC ㆍ ㆍ

Page 54: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 54 -

연소 시스템 설계 해석용 프로그램SCALE

Page 55: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 55 -

요 약 서

항 목 현 황

일 반 현 황

설립일 자본금 원: 2001. 10. 12 : 50,000,000

총 종업원수 명: 7

경영 명: 1

사무 명: 2

연구 기술 명: 4ㆍ

사업장 현황

본사 소재지 서울특별시 성동구 행당동 한양대학교 공업센터: 17

M-307 ( 82 m2)

부설연구소 소재지 서울특별시 성동구 행당동 한양종합기술원: 17

호107-2 ( 35 m2)

개 발 현 황

열 유체 연소 시스템 설계해석용 의 개발- . . S/W(CAC, SCALE)

를 이용한 병렬프로그램 개발- PC cluster

코드와 타 코드와의 개발- Flamelet interface

수치 및 물리모델에 대한 벤치마킹의 완료-

양산체제등

설비구축현황

열 유체 분야의 기술 컨설팅 체제 구축- . .

연구개발 인력 및 공간 확보-

국 내외 관련 기업과의 전략적 제휴- .

홍보 및 마켓팅 기술설명회 개최- ,

사업화 현황

열 유체 연소 특성 및 공해물질 생성 해석 전용 코드인- CACㆍ ㆍ

코드를 판매

자동차 엔진 실린더 내의 유동 및 연소특성의 해석과 관련한-

기술 자문 용역 수주

년 월 현재 현대자동차 한국전력 기계연구원- 2002 6 , , ,

생산기술연구원의 개 업체 및 연구소에서 열 유체 연소4 ㆍ ㆍ

시스템과 관련한 기술 자문을 요청하여 용역 수주절차가 진행되고

있음

환경을 개선하여 상용코드 의 판매를 크게- GUI (CAC, SCALE)

증가시킬 계획임대외기관 인증

특허 외부자금/

활용현황

- 벤처기업 인증 건의 특 허 정책자금 활용 그리고 외부자본의, 1 , ,

유치등을 추진하거나 계획하고 있음

기 타 사 항

Page 56: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 56 -

Page 57: 열 유체 연소시스템설계해석S/Wㆍㆍ 개발및사업화결과보고서 · 이동경계수치절차및자유표면유동해석을위한 모델VOF 평형및비평형상세화학반응모델

- 57 -

이 보고서는 산업자원부 한국산업기술평가원에서

시행한 신기술창업보육사업의 기술개발 및 사업화

결과보고서입니다.