Figure 4.9 Flow analysis results using FLOW3D of the metal flow and solidification in the main cavity. (The velocity is in m/s.)

Numerical Analysis of Die-Casting Process in Thin Cavities Using Lubrication Approximation

Alexandre Reikher
A Dissertation Submitted in
Partial Fulfillment of the
Requirements for the Degree of
Doctor of Philosophy
In Engineering
at
The University of Wisconsin Milwaukee
December 2012

ABSTRACT

얇은 벽 부품의 주조는 오늘날 다이 캐스트 산업의 현실이 되었습니다. 전산 유체 역학 분석은 생산 개발 프로세스의 필수적인 부분입니다. 일반적으로 에너지 방정식과 결합 된 3 차원 Navier-Stokes 방정식은 유동 및 응고 패턴, 유동 선단의 위치, 함수로서 고체-액체 인터페이스의 위치를 ​​이해하기 위해 해결되어야 합니다.

캐비티 충전 및 응고 과정에서 시간. 얇은 벽 주조에 대한 지배 방정식의 일반적인 솔루션에는 많은 수의 계산 셀이 필요하므로 솔루션을 생성하는 데 비현실적으로 오랜 시간이 걸립니다.

Hele Shaw 유동 모델링 접근법을 사용하면 평면 외 유동을 무시함으로써 얇은 캐비티의 유동 문제 해결을 단순화 할 수 있습니다. 추가적인 이점으로, 문제는 3 차원 문제에서 2 차원 문제로 축소됩니다. 그러나 Hele-Shaw 근사는 흐름의 점성력이 관성력보다 훨씬 더 높아야하며,이 경우 Navier-Stokes 방정식은 Reynolds의 윤활 방정식으로 축소됩니다.

그러나 다이 캐스트 공정의 빠른 사출 속도로 인해 관성력을 무시할 수 없습니다. 따라서 윤활 방정식은 흐름의 관성 효과를 포함하도록 수정되어야 합니다.

이 박사 학위 논문에서는 얇은 공동에서 응고와 함께 액체 금속의 정상 상태 및 과도 흐름을 모델링하기 위한 빠른 수치 알고리즘이 개발되었습니다. 설명된 문제는 저온 챔버, 고압 다이 캐스트 공정, 특히 얇은 환기 채널에서 관찰되는 금속 흐름 현상과 밀접한 관련이 있습니다.

채널의 금속 흐름 속도가 고체-액체 계면 속도보다 훨씬 높다는 사실을 사용하여 두께에 따른 열 전달을 처리하면서 금속 흐름을 주어진 시간 단계에서 안정된 것으로 처리하여 새로운 수치 알고리즘을 개발했습니다.

일시적인 방향. 얇은 캐비티의 흐름은 채널 두께에 대한 운동량과 연속성 방정식을 통합 한 후 2 차원으로 처리되고 열 전달은 두께 방향의 1 차원 현상으로 모델링 됩니다. 엇갈린 격자 배열은 유동 지배 방정식을 이산화하는데 사용되며 결과적인 편미분 방정식 세트는 SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 알고리즘을 사용하여 해결됩니다.

상 변화를 수반하는 두께 방향 열 전달 문제는 제어 볼륨 공식을 사용하여 해결됩니다. 고체-액체 계면의 위치와 모양은 솔루션의 일부로 Stefan 조건을 사용하여 찾을 수 있습니다. 시뮬레이션 결과는 응고와 함께 전체 3 차원 흐름 및 열 전달 방정식을 해결하는 상용 소프트웨어 FLOW-3D®의 예측과 잘 비교되는 것으로 나타났습니다.

제안된 수치 알고리즘은 또한 얇은 채널에서 일시적인 금속 충전 및 응고 문제를 해결하기 위해 적용되었습니다. 움직이는 고체-액체 인터페이스의 존재는 이제 반복적으로 해결되는 일련의 흐름 방정식에 비선형 성을 도입합니다.

다시 한번, FLOW3D®의 예측과 잘 일치하는 것이 관찰되었습니다.

이 두 연구는 제안 된 관성 수정 레이놀즈의 윤활 방정식과 함께 두께를 통한 열 손실 및 응고 모델을 성공적으로 구현하여 다이 캐스트 공정 중에 얇은 채널에서 액체 금속의 유동 및 응고에 대한 빠른 분석을 제공 할 수 있음을 나타냅니다. CPU 시간을 대폭 절약하여 얻은 이러한 시뮬레이션 결과는 다이 캐스트 다이의 환기 채널을 설계하는 동안 빠른 초기 분석을 제공하는 데 사용할 수 있습니다.

Figure 1.3. Schematic representation of steps in the hot chamber die-cast process: a.  plunger pushes metal from the sleeve through the gating system into the cavity; b. after  solidification process is complete, the die opens; c. the part is ejected from the cavity.
Figure 1.3. Schematic representation of steps in the hot chamber die-cast process: a. plunger pushes metal from the sleeve through the gating system into the cavity; b. after solidification process is complete, the die opens; c. the part is ejected from the cavity.
Figure 1.5. Schematic representation of steps in the cold chamber die-cast process: a.  molten metal is ladled into the shot sleeve; b. hydraulic cylinder applies pressure on  plunger; c. plunger pushes metal from the sleeve through the gating system into the  cavity; d. high pressure is maintained during solidification; e. after solidification is  complete, the die opens; f. the part is ejected from the cavity.
Figure 1.5. Schematic representation of steps in the cold chamber die-cast process: a. molten metal is ladled into the shot sleeve; b. hydraulic cylinder applies pressure on plunger; c. plunger pushes metal from the sleeve through the gating system into the cavity; d. high pressure is maintained during solidification; e. after solidification is complete, the die opens; f. the part is ejected from the cavity.
Figure 4.6 A schematic of a die-cast die with shot sleeve and plunger: 1) Shot  sleeve, 2) Plunger, 3) Stationary half of the die-cast die, 4) Ejector half of the die-cast die,  5) Mold cavity, 6) Ventilation channel.
Figure 4.6 A schematic of a die-cast die with shot sleeve and plunger: 1) Shot sleeve, 2) Plunger, 3) Stationary half of the die-cast die, 4) Ejector half of the die-cast die, 5) Mold cavity, 6) Ventilation channel.
Figure 4.8 A picture (a ‘full shot’) of a part made using the die-cast process. The  overflows are created when the metal front, after filling the main cavity, fills up the  machined ‘overflow’ pockets in the die-cast mold. Ventilation channel is last to fill-up.
Figure 4.8 A picture (a ‘full shot’) of a part made using the die-cast process. The overflows are created when the metal front, after filling the main cavity, fills up the machined ‘overflow’ pockets in the die-cast mold. Ventilation channel is last to fill-up.
Figure 4.9 Flow analysis results using FLOW3D of the metal flow and solidification in the main cavity. (The velocity is in m/s.)
Figure 4.9 Flow analysis results using FLOW3D of the metal flow and solidification in the main cavity. (The velocity is in m/s.)
Figure 4.10 Temperature distribution in the considered cavity of the die-cast die, filled  with liquid metal at the end of the fill process. (The temperature is in 0C.)
Figure 4.10 Temperature distribution in the considered cavity of the die-cast die, filled with liquid metal at the end of the fill process. (The temperature is in 0C.)
Figure 4.16 Experimentally observed solidified metal in the ventilation channel; a)  Measured length of metal flow in the ventilation channel after solidification stops it; b)  Enlarged image of the solidified metal in the channel
Figure 4.16 Experimentally observed solidified metal in the ventilation channel; a) Measured length of metal flow in the ventilation channel after solidification stops it; b) Enlarged image of the solidified metal in the channel
Figure 3-3 E200U+ PMD sponge draining

CORRELATIONS BETWEEN NEUTRAL BUOYANCY TESTS AND CED

Abstract

NEUTRAL BUOYANCY 테스트는 표면 장력 장치의 기능 검증을 위해 Matra Marconi Space에서 잘 알려진 테스트 기능입니다. 새로운 Eurostar 3000 PMD를 검증하고 성능을 E2000 + PMD의 성능과 비교하기 위해 완전한 저수지가 중성 부력 테스트 벤치에서 테스트되었습니다. 또한, 일시적 및 고정 모세관 현상에 대한 우리의 지식을 연관시키고 개선하기 위해 수치 시뮬레이션이 수행되었습니다. 다양한 결과에 대한 토론과 상관 관계가 제시됩니다. 다양한 도구의 한계가 표시되고 E3000 개선의 간단한 예가 모델링됩니다.

Figure 2-2 Propellant position under acceleration
Figure 2-2 Propellant position under acceleration
Figure 2-3 E2000+ PMD
Figure 2-3 E2000+ PMD
Figure 3-1 Propellant position under acceleration
Figure 3-1 Propellant position under acceleration
Figure 3-2 Propellant position under acceleration
Figure 3-2 Propellant position under acceleration
Figure 3-3 E200U+ PMD sponge draining
Figure 3-3 E200U+ PMD sponge draining
Figure 3-5 Improved PMD sponge draining
Figure 3-5 Improved PMD sponge draining
Figure 4-1 E2000+ reservoir
Figure 4-1 E2000+ reservoir
Figure 1: Die configuration for a multi-attribute composite die for high die life and self-lubricating surface

Innovative Die Material and Lubrication Strategies for Clean and Energy Conserving Forging Technologies

청정 및 에너지 절약 단조 기술을 위한 혁신적인 다이 재료 및 윤활 전략

이 최종 기술 보고서에는 수상 번호 DE-FC07-01ID14206에 따라 미국 에너지 부에서 부분적으로 자금을 지원 한 “청정 및 에너지 절약 단조 기술을위한 혁신적인 다이 재료 및 윤활 전략”프로젝트에서 수행 된 작업이 포함되어 있습니다. 프로젝트 수행을위한 계약 시간은 2001 년 9 월 30 일부터 2005 년 9 월 29 일까지였습니다. 그러나 DOE / OIT는 2003 년과 2004 년 회계 연도 지난 2 년 동안 자금을 제공 할 수 없었고 프로젝트는 2003-04 회계 연도에 조기 종료되었습니다. 결과적으로 많은 주요 연구 과제가 특정 이정표를 달성하기 위해 수정되거나 완료되지 않고 종료되었습니다. Ohio State University의 산업, 용접 및 시스템 공학 교수 인 Rajiv Shivpuri 박사는이 프로젝트의 프로젝트 책임자이자 수석 조사자였습니다. 이상은 오하이오 주립 대학 연구 재단 (OSURF)에서 관리했습니다. OSURF는 모든 재정 및 행정 문제도 담당했습니다. 재정 보고서는 별도로 제출됩니다. 에너지 부서, 산업 기술 사무소의 프로그램 관리자는 Golden Office의 Mr. Ramesh Jain과 Mr. Dibyajyoti Aichbhowmik이었습니다.
이 프로젝트의 주요 성과는 다음과 같습니다.

• 단조 산업 및 해당 공급 업체와 함께 산업 응용 분야를위한 혁신적인 다이 재료 및 윤활 전략을 탐색하기위한 주요 협력 노력이 수립되었습니다. 여기에는 단조 산업과 협력하는 워크숍과 심포지엄이 포함되었습니다. 단조 산업 전체에 결과를 전파하기 위해 단조 산업 기술 컨퍼런스에서 발표되었습니다.

• 단조 산업 협회와 단조 산업 교육 연구 재단의 후원으로 단조 기술 우수 센터 설립. 이 센터의 일부로 산업, OSU, 오하이오 주 및 DOE 지원과 함께 2 개의 단조 셀이 설치되었습니다. 1300 톤 기계식 프레스 셀과 350 톤 유압 프레스 셀입니다. 이것은 단조 연구에 150 만 달러를 투입 한 것입니다.

• LENS (Laser Enhanced Net Shaping) 기반 니켈 알루미나 이드 코팅 오버레이 (자세한 내용은 부록 A 참조)를 포함하여 혁신적인 다이 코팅이 탐색되었습니다.

• 열간 단조 응용 분야를위한 금형 재료를 최적으로 선택하고 설계하기 위해 혁신적인 실험 설정 및 예측 열 연화 소프트웨어가 개발되었습니다 (부록 B, C 및 D).

• 윤활 전략 및 단일 액적 기반 윤활 모델은 확산 및 열 전달을위한 열간 단조 윤활제의 최적 증착을 위해 개발되었습니다 (부록 E 및 F).

• 윤활유 분해 및 바운스 용 모델이 개발되었습니다. 이 모델은 뜨거운 다이 표면의 흑연 윤활로 인한 공기 및 지하수 오염을 줄이는 데 사용할 수 있습니다.

(부록 G). 이 보고서는 Shivpuri 박사와 Yijun Zhu (연구원)가 작성했습니다. 여기에는 다른 외부 또는 내부 지원과 함께 프로젝트 종료 후 일부 연구 계획 및 프로젝트 기간 동안 완료된 작업에 대한 세부 정보가 포함되어 있습니다.

1.1 프로젝트 목표

이 프로젝트의 목표는 혁신적인 다이 재료 및 윤활 전략을 개발 및 구현하여 다이 수명을 8 배 늘리고, 에너지 투입량을 15 % 줄이며, 부품 당 에너지 비용을 50 % 줄이며, 윤활유에서 나오는 미립자 배출량을 90 % 줄이며, 다이 관련 가동 시간을 90 %까지 늘립니다.

단조 산업, 공급 업체 (철강 및 알루미늄 생산 업체 (IOF), 윤활유, 표면 기술 및 다이 소재 공급 업체) 및 고객 (OEM)에 미치는 최대의 광범위한 에너지 영향을 위해 전략이 선택되었습니다.

여기에는 최적의 윤활제 스프레이 기술, 고급 표면 엔지니어링에 의한 열간 단조의 흑연 제거, 경사 다이 재료 및 다이 엔지니어링, 열간 단조를위한 윤활 및 다이 활성화 등이 포함됩니다.

미국의 단조 산업은 1997 년에 약 120 억 달러였습니다 (DOD 국가 안보). 평가). 제품 총 판매 가치의 약 15 %가 에너지에 할당되며 연간 약 50 조 BTU입니다. 흑연 사용 (열간 단조) 및 냉간 단조 전환 코팅 사용으로 인한 환경 영향은 제품 비용에 20 % 이상 추가 될 것으로 예상됩니다.

Figure 1: Die configuration for a multi-attribute composite die for high die life and self-lubricating surface
Figure 1: Die configuration for a multi-attribute composite die for high die life and self-lubricating surface

BACKGROUND

실온 (저온) 및 고온 (온 및 고온)에서 수행되는 단조는 진화하는 야금, 공구 표면의 마찰 및 금속의 흐름 특성을 포함하는 잘 이해되지 않는 복잡한 현상입니다. 이 프로젝트에서 다루어 진 기술적 장애물은 다음과 같습니다.

• 냉간 및 열간 단조의 윤활 작용에 대한 지식 부족. 윤활유 및 윤활 기술의 선택은 윤활유 및 장비 공급 업체에 맡겨집니다. 이로 인해 윤활유의 과도하고 불량한 사용과 과도한 환경 오염이 발생합니다.

• 고급 단조 응용 분야를위한 새로운 표면 엔지니어링 및 다이 재료 기술의 성숙도가 부족합니다. 실제 생산에서이를 구현하는 데 따른 기술적 및 재정적 위험이 매우 높아 사용을 제한합니다. 이러한 기술의 시장 침투는 거의 존재하지 않습니다.

• 다이와 윤활 시스템의 설계 최적화를위한 계산 도구가 부족합니다.

윤활유 및 다이 소재 기술에서 다음과 같은 전략을 통해 프로젝트 목표를 실현할 계획이었습니다.

• 전략 # 1 : 오염을 제거하고, 윤활제 사용을 줄이며, 다이 냉각 감소로 인한 그물 성형을 가능하게하는 윤활제 스프레이 공정의 최적 설계를위한 시스템 개발. 또한 흑연 기반 윤활유의 필요성을 줄여줍니다.

• 전략 # 2 : 철 및 비철 부품의 온간 단조 (빌릿 가열이 1250F에서 900F로 감소)를위한 다이 수명과 공정을 개선하기 위한 윤활제 및 다이 코팅 가능 요소를 개발합니다. 단조 온도를 낮추면 공차가 개선되고 부품 당 에너지가 크게 절약됩니다.

• 전략 # 3 : 저 마찰 다이 표면 엔지니어링 (DLC (비철) 및 WC / C 코팅)을 사용하여 냉간 단조 빌릿에 인광 코팅을 사용하지 않습니다.

• 전략 # 4 : 열간 단조 금형을위한 고급 표면 클래딩 (렌즈 및 열 스프레이에 의한 단단한 표면) 및 이중 코팅 기술을 개발합니다. 기존의 코팅과 표면 공학 기술은 상당한 이득을 얻지 못했습니다.

• 전략 # 5 : 재료 및 공정 설계를 통해 냉간 및 열간 단조에서 공정 중 다이 고장을 제거하고 예측 다이 유지 보수를위한 소프트웨어를 개발합니다. 이는 스크랩 감소 및 다이 관련 다운 타임에 상당한 영향을 미칩니다.

개발중인 많은 기술은 수치 모델링, 윤활 및 냉각수 기술, 표면 기술, 재료의 신속한 프로토 타이핑, 레이저 기술 등과 같은 교차 절단 R & D 가능 요소를 다루고 있습니다. 이러한 기술은 지원 산업의 로드맵에서도 중요한 기술로 확인되었습니다.

미래의 산업으로. IOF를 위해 250 조 BTU의 에너지 절약과 3500 톤의 오염 물질이 예상됩니다. 프로젝트가 전액 지원을받지 못하고 프로젝트가 2004 년 9 월 30 일에 종료되었으므로 전략 # 1, # 4, # 5 만 추구했습니다. 연구 및 구현에 대한 세부 사항은 부록에 포함되어 있습니다.

Effect of lubricant heat

템퍼링, 마모 및 공구 열화에 대한 단조 윤활유의 효과를 평가하기 위해 다양한 열 전달 계수로 여러 시뮬레이션을 수행했습니다. 컴퓨터 시뮬레이션에 사용 된 열전달 계수의 값은 얻은 값과 일치하며 경우에 따라 Sridhar 등이 오하이오 주립 대학에서 수행 한 테스트에서 추정 한 값입니다. 사용 된 계면 열전달 계수의 값은 12 KW / m2 ° C, 24 KW / m2 ° C 및 33 KW / m2 ° C였으며, 이는 20 부, 30 부 및 100 부 물로 희석 된 수성 흑연 윤활제에 해당합니다 (희석 비율 1:20, 1:30 및 1 : 100). 이러한 각 희석 비율에 대해 3000 및 5000 샷 후 상부 다이의 경도 분포는 그림 C.3, C.4 및 C.5에 나와 있습니다. 희석 비 1:20에 대한 표면 경도 분포는 그림 C.6에 나와 있습니다.

Figure C. 2: stage gear blank forging sequence (Courtesy: Sypris Technologies
Figure C. 2: stage gear blank forging sequence (Courtesy: Sypris Technologies
Figure C. 3: Hardness distribution after 3000 and 5000 shots, heat transfer coefficient used = 12 KW/m2°C, press type: mechanical press
Figure C. 3: Hardness distribution after 3000 and 5000 shots, heat transfer coefficient used = 12 KW/m2°C, press type: mechanical press

F.5.3 Results of the Lubricant Properties

표 F.1은 윤활유의 측정 된 특성을 보여줍니다. DP는 107 및 CA 모세관 작용 방법에서 펜던트 드롭 방법을 나타냅니다. 테스트 된 액체에는 순수한 물이 포함됩니다. 다음과 같은 사실을 관찰 할 수 있습니다. a). 더 높은 표면 장력을 가진 더 높은 희석 비율 회사; 비). 희석 비율이 1 : 1보다 큰 액체의 경우 표면 장력이 물의 장력에 접근합니다. 드롭 펜던트 법으로 추정 한 모든 표면 장력은 동일한 경향을 공유하지만 약 10dynes / cm에 대해 모세관 작용법에 의한 것보다 작다는 것을 알 수 있습니다. 물의 표면 장력이 72.8dynes / cm라는 점을 감안할 때 모세관 작용법에서 얻은 결과가 실제 값에 더 가깝다고 생각합니다.

Figure F. 10: simulation results of lubricant 1:1 with 4mm diameter droplet at impact velocity 10cm/s.
Figure F. 10: simulation results of lubricant 1:1 with 4mm diameter droplet at impact velocity 10cm/s.
Figure F. 12: Experimental results of maxξ v.s. TD. We = 27.
Figure F. 12: Experimental results of maxξ v.s. TD. We = 27.
Figure G. 1: Dryoff process of a lubricant droplet at film boiling: (a)- (c) fluid dynamic process, (d). quasi-steady dryoff process.
Figure G. 1: Dryoff process of a lubricant droplet at film boiling: (a)- (c) fluid dynamic process, (d). quasi-steady dryoff process.
Figure 2.6 ESI apparatus for offline analysis with microscope imaging.

MODELING AND CHARACTERIZATION OF MICROFABRICATED EMITTERS: IN PURSUIT OF IMPROVED ESI-MS PERFORMANCE

미세 가공 방사체의 모델링 및 특성화 : 개선된 ESI-MS 성능 추구

by XINYUN WU

A thesis submitted to the Department of Chemistry in conformity with the requirements for the degree of Master of Science Queen’s University Kingston, Ontario, Canada December, 2011 Copyright © Xinyun Wu, 2011

Abstract

ESI (Electrospray ionization)는 특히 탁월한 감도, 견고성 및 단순성으로 대형 생체 분자를 분석하는 데있어 질량 분석 (MS)에 매우 귀중한 기술이었습니다. ESI 기술 개발에 많은 노력을 기울였습니다. 그 형태와 기하학적 구조가 전기 분무 성능과 추가 MS 감지에 중추적 인 것으로 입증 되었기 때문입니다.

막힘 및 낮은 처리량을 포함하여 전통적인 단일 홀 이미터의 본질적인 문제는 기술의 적용 가능성을 제한합니다. 이 문제를 해결하기 위해 현재 프로젝트는 향상된 ESI-MS 분석을위한 다중 전자 분무(MES) 방출기를 개발하는데 초점을 맞추고 있습니다.

이 논문에서는 스프레이 전류 측정을 위한 전기 분무와 오프라인 전기 분무 실험을 위한 전산 유체 역학 (CFD) 시뮬레이션의 공동 작업이 수행되었습니다. 전기 분무 성능에 대한 다양한 이미터 설계의 영향을 테스트하기 위해 수치 시뮬레이션이 사용되었으며 실험실 결과는 가이드 및 검증으로 사용되었습니다.

CFD 코드는 Taylor-Melcher 누설 유전체 모델(LDM)을 기반으로 하며 과도 전기 분무 공정이 성공적으로 시뮬레이션되었습니다.

이 방법은 750 μm 내경 (i.d.) 이미 터를 통해 먼저 검증되었으며 20 μm i.d.에 추가로 적용되었습니다. 모델. 전기 분무 공정의 여러 단계가 시각적으로 시연되었으며 다양한 적용 전기장 및 유속에서 분무 전류의 변화에 ​​대한 정량적 조사는 이전 시뮬레이션 및 측정과 잘 일치합니다.

단일 조리개 프로토 타입을 기반으로 2 홀 및 3 홀 이미터로 MES 시뮬레이션을 수행했습니다. 시뮬레이션 예측은 실험 결과와 유사하게 비교되었습니다. 이 작업의 증거는 CFD 시뮬레이션이 MES의 이미 터 설계를 테스트하는 효과적인 수치 도구로 사용될 수 있음을 입증했습니다.

이 작업에서 달성 된 마이크로 스케일 에미 터 전기 분무의 성공적인 시뮬레이션에 대한 벤치마킹 결과는 현재까지 발표 된 전기 분무에 대한 동적 시뮬레이션의 가장 작은 규모로 여겨집니다.

Co-Authorship

공동 저자: 이 논문에 대한 모든 연구는 Natalie M. Cann 박사와 Richard D. Oleschuk 박사의 지도하에 완료되었습니다. 다중 전자 분무에 관한 4 장에서 제시된 연구 작업의 일부는 Ramin Wright가 공동 저술했으며, 이 작업은 press에서 다음 논문에서 인용되었습니다.

ibson,G.T.T.; Wright, R.D.; Oleschuk, R.D. Multiple electrosprays generated from a single poly carbonate microstructured fibre. Journal of Mass Spectrometry, 2011, in press.

Chapter 1 Introduction

소프트 이온화 방법으로 ESI (electrospray ionization)의 도입은 질량 분석법 (MS)의 적용 가능성에 혁명을 일으켰습니다. 이 기술의 부드러운 특징은 상대적으로 높은 전하를 가진 이온을 생성하는 고유한 이점으로 인해 액상에서 직접 펩티드 및 단백질과 같은 큰 생체 분자를 분석 할 수 있게했습니다 [1].

지난 10 년 동안 ESI-MS는 놀라운 성장을 보였으며 현재는 단백질 체학, 대사 체학, 글리코 믹스, 합성 화학자를 위한 식별 도구 등 다양한 생화학 분야에서 광범위하게 채택되고 있습니다 [2-3].

ESI-MS는 겔 전기 영동과 같은 생물학적 분자에 대한 기존의 질량 측정 기술보다 훨씬 빠르고 민감하며 정확합니다. 또한, 액체상에서 직접 분석 할 수 있는 큰 비 휘발성 분자의 능력은 고성능 액체 크로마토 그래피 (HPLC) 및 모세관 전기 영동 (CE)과 같은 업스트림 분리 기술과의 결합을 가능하게합니다 [4].

일반적인 ESI 공정은 일반적으로 액적 형성, 액적 수축 및 기상 이온의 최종 형성을 포함합니다. 일렉트로 스프레이의 성능에 영향을 미치는 많은 요소 중에서 스프레이를 위한 이미터의 구조 (즉, 기하학, 모양 등)가 중요한 요소입니다.

전통적인 전기 분무 이미터는 일반적으로 풀링 또는 에칭 기술로 제작 된 단일 채널 테이퍼 형 또는 비 테이퍼 형입니다. 그러나 이러한 이미터는 종종 막힘, 부적절한 처리량 등과 같은 문제로 어려움을 겪습니다. [5]

향상된 감도 및 샘플 활용을 위해 다중 스프레이를 생성하는 새로운 이미터 설계 개발로 분명한 발전이 있었습니다. 새로운 ESI 이미터 설계에 대한 연구는 실험적으로나 이론적으로 큰 관심을 불러 일으켰습니다 [3]. 그러나 ESI의 복잡한 물리적 과정은 팁 형상 외에도 많은 다른 변수에 의존하기 때문에 연구간 직접 비교의 어려움은 장애물이 됩니다.

또한 새로운 나노 이미터 제조 및 테스트 비용이 상당히 높을 수 있습니다. 이 논문은 CFD 시뮬레이션 도구를 활용하여 가상 랩을 설정함으로써 이러한 문제를 해결합니다. 다른 매개 변수로 인해 상호 연결된 변경 없이 다양한 이미터 설계를 비교할 수 있도록 이상적으로 균일한 물리적 조건을 제공합니다.

맞춤 제작된 프로토 타입의 실험 측정 값도 수집되어 더 나은 계산 체계를 형성하는 데 도움이 되는 지침과 검증을 모두 제공합니다. 특히 이 분야의 주요 미래 플랫폼으로 여겨지는 다중 노즐 이미 터 설계에 중점을 둘 것입니다.

전기 분무 거동에 영향을 미치는 요인에 대한 추가 기본 연구는 다양한 기하학적 및 작동 매개 변수와 관련하여 수행됩니다. 이는 보다 효율적이고 견고한 이미터의 개발을 가능하게 할 뿐만 아니라 더 넓은 영역에서 ESI의 적용을 향상시킬 수 있습니다.

Figure 1.1Schematic setup for ESI-MS technique
Figure 1.1Schematic setup for ESI-MS technique
Figure 1.2 Schematic of major processes occurring in electrospray [5].
Figure 1.2 Schematic of major processes occurring in electrospray [5].
Figure 1.3 Illustration of detailed geometric parameters of a spraying Taylor cone wherera is the radius of curvature of the best fitting circle at the tip of the cone; re is the radius of the emission region for droplets at the tip of a Taylor cone;is the liquid cone angle.
Figure 1.3 Illustration of detailed geometric parameters of a spraying Taylor cone wherera is the radius of curvature of the best fitting circle at the tip of the cone; re is the radius of the emission region for droplets at the tip of a Taylor cone;is the liquid cone angle.
Figure 1.4 (A)Externally tapered emitter  (B) Optical image of a clogged tapered emitter with normal use [46].
Figure 1.4 (A)Externally tapered emitter (B) Optical image of a clogged tapered emitter with normal use [46].
Figure 1.5 (A)Three by three configuration of an emitter array made with polycarbonate using laser ablation; (B) Photomicrograph of nine stable electrosprays generated from the nine-emitter array [52]
Figure 1.5 (A)Three by three configuration of an emitter array made with polycarbonate using laser ablation; (B) Photomicrograph of nine stable electrosprays generated from the nine-emitter array [52]
Figure 1.6 SEM images of the distal ends of four multichannel nanoelectrospray emitters and a tapered emitter: (A) 30 orifice emitter; (B) 54 orifice emitter; (C) 84 orifice emitter; (D) 168 orifice emitter; Scale bars in A, B, and C represent 50 μm, and 100 μm in D[54]
Figure 1.6 SEM images of the distal ends of four multichannel nanoelectrospray emitters and a tapered emitter: (A) 30 orifice emitter; (B) 54 orifice emitter; (C) 84 orifice emitter; (D) 168 orifice emitter; Scale bars in A, B, and C represent 50 μm, and 100 μm in D[54]
Figure 1.7 Photomicrographs of electrospray from of a 168-hole MCN emitter at different flow rates. (A) A traditional integrated Taylor cone observed from offline electrospray of water with 0.1% formic acid at 300 nL/min; (B) A mist of coalesced Taylor cones observed from offline electrospray at 25 nL/min[54]
Figure 1.7 Photomicrographs of electrospray from of a 168-hole MCN emitter at different flow rates. (A) A traditional integrated Taylor cone observed from offline electrospray of water with 0.1% formic acid at 300 nL/min; (B) A mist of coalesced Taylor cones observed from offline electrospray at 25 nL/min[54]
Figure 1.8 Circular arrays of etched emitters for better electric field homogeneity [53].
Figure 1.8 Circular arrays of etched emitters for better electric field homogeneity [53].
Figure 2.6 ESI apparatus for offline analysis with microscope imaging.
Figure 2.6 ESI apparatus for offline analysis with microscope imaging.
Figure 3.9 Typical panel for displaying instant simulation result during simulation process.
Figure 3.9 Typical panel for displaying instant simulation result during simulation process.
Figure 5.3 Generation of a Taylor cone-jet mode (simulation) plotted with iso-potential lines at times    (Top to bottom panels correspond to 0.002 s, 0.012 s, 0.018 s, 0.08 s respectively).
Figure 5.3 Generation of a Taylor cone-jet mode (simulation) plotted with iso-potential lines at times (Top to bottom panels correspond to 0.002 s, 0.012 s, 0.018 s, 0.08 s respectively).
Figure 5.8 (A) Taylor cone-jet profiles with different contact angle of 30 degrees and 20 degrees (B) under the same physical conditions of 6 kV and 0.04 m/s. (C) Cone-jet profile generated from a tapered tip with a 20 degree contact angle at 6 kV and 0.04 m/s (as a comparison with (B)).
Figure 5.8 (A) Taylor cone-jet profiles with different contact angle of 30 degrees and 20 degrees (B) under the same physical conditions of 6 kV and 0.04 m/s. (C) Cone-jet profile generated from a tapered tip with a 20 degree contact angle at 6 kV and 0.04 m/s (as a comparison with (B)).

Omit below: Please refer to the original text for the full content.

Bibliography

1. Mclafferty, F.W., Tandem Fourier-Transform Mass-Spectrometry of Large Molecules.Abstracts of Papers of the American Chemical Society, 1986. 192: p. 21-Anyl. 2. Griffiths, W.J. and Y.Q. Wang, Mass spectrometry: from proteomics to metabolomics and lipidomics. Chemical Society Reviews, 2009. 38(7): p. 1882-1896. 3. Gibson, G.T.T., S.M. Mugo, and R.D. Oleschuk, Nanoelectrospray Emitters: Trends and Perspective. Mass Spectrometry Reviews, 2009. 28(6): p. 918-936. 4. Cech, N.B. and C.G. Enke, Practical implications of some recent studies in electrospray ionization fundamentals. Mass Spectrometry Reviews, 2001. 20(6): p. 362-387. 5. Su, S., Development and Application of Non-tapered Electrospray Emitters for Nano-ESI Mass Spectrometry, in Chemistry. 2008, Queen’s University: Kingston. p. 185. 6. Zeleny, J., The electrical discharge from liquid points, and a hydrostatic method of measuring the electric intensity at their surfaces. Physical Review, 1914. 3(2): p. 69-91. 7. Dole, M., L.L. Mack, and R.L. Hines, Molecular Beams of Macroions. Journal of Chemical Physics, 1968. 49(5): p. 2240-&. 8. Yamashita, M. and J.B. Fenn, Negative-Ion Production with the Electrospray Ion-Source.Journal of Physical Chemistry, 1984. 88(20): p. 4671-4675. 9. Kebarle, P. and U.H. Verkerk, Electrospray: From Ions in Solution to Ions in the Gas Phase, What We Know Now. Mass Spectrometry Reviews, 2009. 28(6): p. 898-917. 10. Taylor, G., Disintegration of Water Drops in Electric Field. Proceedings of the Royal Society of London Series a-Mathematical and Physical Sciences, 1964. 280(138): p. 383. 11. Cole, R.B., Some tenets pertaining to electrospray ionization mass spectrometry. Journal of Mass Spectrometry, 2000. 35(7): p. 763-772. 12. Rayleigh, L., On the equilibrium of liquid conducting masses charged with electricity.Philos. Mag., 1882. 14: p. 184-186. 13. Mack, L.L., et al., Molecular Beams of Macroions .2. Journal of Chemical Physics, 1970. 52(10): p. 4977-&. 14. Gamero-Castano, M. and J.F. de la Mora, Kinetics of small ion evaporation from the charge and mass distribution of multiply charged clusters in electrosprays. Journal of Mass Spectrometry, 2000. 35(7): p. 790-803. 15. Gamero-Castano, M. and J.F. de la Mora, Modulations in the abundance of salt clusters in electrosprays. Analytical Chemistry, 2000. 72(7): p. 1426-1429. 16. Loscertales, I.G. and J.F. Delamora, Experiments on the Kinetics of Field Evaporation of Small Ions from Droplets. Journal of Chemical Physics, 1995. 103(12): p. 5041-5060. 17. Rohner, T.C., N. Lion, and H.H. Girault, Electrochemical and theoretical aspects of electrospray ionisation. Physical Chemistry Chemical Physics, 2004. 6(12): p. 3056-3068.

18. Iribarne, J.V. and B.A. Thomson, Evaporation of Small Ions from Charged Droplets.Journal of Chemical Physics, 1976. 64(6): p. 2287-2294. 19. Meng, C.K. and J.B. Fenn, Formation of Charged Clusters during Electrospray Ionization of Organic Solute Species. Organic Mass Spectrometry, 1991. 26(6): p. 542-549. 20. Nohmi, T. and J.B. Fenn, Electrospray Mass-Spectrometry of Poly(Ethylene Glycols) with Molecular-Weights up to 5 Million. Journal of the American Chemical Society, 1992. 114(9): p. 3241-3246. 21. de la Mora, J.F., Electrospray ionization of large multiply charged species proceeds via Dole’s charged residue mechanism. Analytica Chimica Acta, 2000. 406(1): p. 93-104. 22. Iavarone, A.T., J.C. Jurchen, and E.R. Williams, Supercharged protein and peptide lone formed by electrospray ionization. Analytical Chemistry, 2001. 73(7): p. 1455-1460. 23. Hogan, C.J., et al., Charge carrier field emission determines the number of charges on native state proteins in electrospray ionization. Journal of the American Chemical Society, 2008. 130(22): p. 6926-+. 24. Nguyen, S. and J.B. Fenn, Gas-phase ions of solute species from charged droplets of solutions. Proceedings of the National Academy of Sciences of the United States of America, 2007. 104(4): p. 1111-1117. 25. Luedtke, W.D., et al., Nanojets, electrospray, and ion field evaporation: Molecular dynamics simulations and laboratory experiments. Journal of Physical Chemistry A, 2008. 112(40): p. 9628-9649. 26. Enke, C.G., A predictive model for matrix and analyte effects in electrospray ionization of singly-charged ionic analytes. Analytical Chemistry, 1997. 69(23): p. 4885-4893. 27. Maze, J.T., T.C. Jones, and M.F. Jarrold, Negative droplets from positive electrospray.Journal of Physical Chemistry A, 2006. 110(46): p. 12607-12612. 28. Kebarle, P. and M. Peschke, On the mechanisms by which the charged droplets produced by electrospray lead to gas phase ions. Analytica Chimica Acta, 2000. 406(1): p. 11-35. 29. Loeb, L.B., A.F. Kip, and G.G. Hudson, Pulses in negative point-to-plane corona.Physical Review, 1941. 60(10): p. 714-722. 30. Cole, R.B., Electrospray ionization mass spectrometry : fundamentals, instrumentation, and applications. 1997, New York: Wiley. xix, 577 p. 31. Smith, D.P.H., The Electrohydrodynamic Atomization of Liquids. Ieee Transactions on Industry Applications, 1986. 22(3): p. 527-535. 32. Taylor, G.I. and A.D. Mcewan, Stability of a Horizontal Fluid Interface in a Vertical Electric Field. Journal of Fluid Mechanics, 1965. 22: p. 1-&. 33. Ikonomou, M.G., A.T. Blades, and P. Kebarle, Electrospray Mass-Spectrometry of Methanol and Water Solutions Suppression of Electric-Discharge with Sf6 Gas. Journal of the American Society for Mass Spectrometry, 1991. 2(6): p. 497-505.

34. Wampler, F.M., A.T. Blades, and P. Kebarle, Negative-Ion Electrospray Mass-Spectrometry of Nucleotides – Ionization from Water Solution with Sf6 Discharge Suppression. Journal of the American Society for Mass Spectrometry, 1993. 4(4): p. 289-295. 35. Marginean, I., P. Nemes, and A. Vertes, Order-chaos-order transitions in electrosprays: The electrified dripping faucet. Physical Review Letters, 2006. 97(6): p. -. 36. Marginean, I., P. Nemes, and A. Vertes, Astable regime in electrosprays. Physical Review E, 2007. 76(2): p. -. 37. Nemes, P., I. Marginean, and A. Vertes, Spraying mode effect on droplet formation and ion chemistry in electrosprays. Analytical Chemistry, 2007. 79(8): p. 3105-3116. 38. Marginean, I., et al., Electrospray characteristic curves: In pursuit of improved performance in the nanoflow regime. Analytical Chemistry, 2007. 79(21): p. 8030-8036. 39. Page, J.S., et al., Subambient pressure ionization with nanoelectrospray source and interface for improved sensitivity in mass spectrometry. Analytical Chemistry, 2008. 80(5): p. 1800-1805. 40. Delamora, J.F. and I.G. Loscertales, The Current Emitted by Highly Conducting Taylor Cones. Journal of Fluid Mechanics, 1994. 260: p. 155-184. 41. Ganan-Calvo, A.M., On the general scaling theory for electrospraying. Journal of Fluid Mechanics, 2004. 507: p. 203-212. 42. Smith, D.R., G. Sagerman, and T.D. Wood, Design and development of an interchangeable nanomicroelectrospray source for a quadrupole mass spectrometer.Review of Scientific Instruments, 2003. 74(10): p. 4474-4477. 43. Barnidge, D.R., S. Nilsson, and K.E. Markides, A design for low-flow sheathless electrospray emitters. Analytical Chemistry, 1999. 71(19): p. 4115-4118. 44. Guzzetta, A.W., R.A. Thakur, and I.C. Mylchreest, A robust micro-electrospray ionization technique for high-throughput liquid chromatography/mass spectrometry proteomics using a sanded metal needle as an emitter. Rapid Communications in Mass Spectrometry, 2002. 16(21): p. 2067-2072. 45. Wilm, M. and M. Mann, Analytical properties of the nanoelectrospray ion source.Analytical Chemistry, 1996. 68(1): p. 1-8. 46. Covey, T.R. and D. Pinto, Practical Spectroscopy. Vol. 32. 2002. 47. Kelly, R.T., et al., Nanoelectrospray emitter arrays providing interemitter electric field uniformity. Analytical Chemistry, 2008. 80(14): p. 5660-5665. 48. Choi, Y.S. and T.D. Wood, Polyaniline-coated nanoelectrospray emitters treated with hydrophobic polymers at the tip. Rapid Communications in Mass Spectrometry, 2007. 21(13): p. 2101-2108. 49. Tojo, H., Properties of an electrospray emitter coated with material of low surface energy. Journal of Chromatography A, 2004. 1056(1-2): p. 223-228.

50. Liu, J., et al., Electrospray ionization with a pointed carbon fiber emitter. Analytical Chemistry, 2004. 76(13): p. 3599-3606. 51. Sen, A.K., et al., Modeling and characterization of a carbon fiber emitter for electrospray ionization. Journal of Micromechanics and Microengineering, 2006. 16(3): p. 620-630. 52. Tang, K.Q., et al., Generation of multiple electrosprays using microfabricated emitter arrays for improved mass spectrometric sensitivity. Analytical Chemistry, 2001. 73(8): p. 1658-1663. 53. Deng, W. and A. Gomez, Influence of space charge on the scale-up of multiplexed electrosprays. Journal of Aerosol Science, 2007. 38(10): p. 1062-1078. 54. Su, S.Q., et al., Microstructured Photonic Fibers as Multichannel Electrospray Emitters.Analytical Chemistry, 2009. 81(17): p. 7281-7287. 55. Sen, A.K., J. Darabi, and D.R. Knapp, Simulation and parametric study of a novel multi-spray emitter for ESI-MS applications. Microfluidics and Nanofluidics, 2007. 3(3): p. 283-298. 56. Hayati, I., A. Bailey, and T.F. Tadros, Investigations into the Mechanism of Electrohydrodynamic Spraying of Liquids .2. Mechanism of Stable Jet Formation and Electrical Forces Acting on a Liquid Cone. Journal of Colloid and Interface Science, 1987. 117(1): p. 222-230. 57. Glonti, G.A., On the Theory of the Stability of Liquid Jets in an Electric Field. Soviet Physics Jetp-Ussr, 1958. 7(5): p. 917-918. 58. Nayyar, N.K. and G.S. Murty, The Stability of a Dielectric Liquid Jet in the Presence of a Longitudinal Electric Field. Proceedings of the Physical Society of London, 1960. 75(483): p. 369-373. 59. Allan, R.S. and S.G. Mason, Particle Behaviour in Shear and Electric Fields .1. Deformation and Burst of Fluid Drops. Proceedings of the Royal Society of London Series a-Mathematical and Physical Sciences, 1962. 267(1328): p. 45-&. 60. Melcher, J.R. and G.I. Taylor, Electrohydrodynamics – a Review of Role of Interfacial Shear Stresses. Annual Review of Fluid Mechanics, 1969. 1: p. 111-&. 61. Saville, D.A., Electrohydrodynamics: The Taylor-Melcher leaky dielectric model. Annual Review of Fluid Mechanics, 1997. 29: p. 27-64. 62. Carretero Benignos, J.A. and Massachusetts Institute of Technology. Dept. of Mechanical Engineering., Numerical simulation of a single emitter colloid thruster in pure droplet cone-jet mode. 2005. p. 117 leaves. 63. Hartman, R.P.A., et al., The evolution of electrohydrodynamic sprays produced in the cone-jet mode, a physical model. Journal of Electrostatics, 1999. 47(3): p. 143-170. 64. Hartman, R.P.A., et al., Electrohydrodynamic atomization in the cone-jet mode physical modeling of the liquid cone and jet. Journal of Aerosol Science, 1999. 30(7): p. 823-849.

65. Yoon, S.S., et al., Modeling multi-jet mode electrostatic atomization using boundary element methods. Journal of Electrostatics, 2001. 50(2): p. 91-108. 66. Zeng, J., D. Sobek, and T. Korsmeyer, Electro-hydrodynamic modeling of electrospray ionization: Cad for a mu fluidic device – Mass spectrometer interface. Boston Transducers’03: Digest of Technical Papers, Vols 1 and 2, 2003: p. 1275-1278, 1938. 67. Lastow, O. and W. Balachandran, Numerical simulation of electrohydrodynamic (EHD) atomization. Journal of Electrostatics, 2006. 64(12): p. 850-859. 68. http://www.flow3d.com. 69. Valaskovic, G.A., et al., Attomole-Sensitivity Electrospray Source for Large-Molecule Mass-Spectrometry. Analytical Chemistry, 1995. 67(20): p. 3802-3805. 70. Kriger, M.S., K.D. Cook, and R.S. Ramsey, Durable Gold-Coated Fused-Silica Capillaries for Use in Electrospray Mass-Spectrometry. Analytical Chemistry, 1995. 67(2): p. 385-389. 71. Fang, L.L., et al., Online Time-of-Flight Mass-Spectrometric Analysis of Peptides Separated by Capillary Electrophoresis. Analytical Chemistry, 1994. 66(21): p. 3696-3701. 72. Cao, P. and M. Moini, A novel sheathless interface for capillary electrophoresis/electrospray ionization mass spectrometry using an in-capillary electrode. Journal of the American Society for Mass Spectrometry, 1997. 8(5): p. 561-564. 73. Fong, K.W.Y. and T.W.D. Chan, A novel nonmetallized tip for electrospray mass spectrometry at nanoliter flow rate. Journal of the American Society for Mass Spectrometry, 1999. 10(1): p. 72-75. 74. Emmett, M.R. and R.M. Caprioli, Micro-Electrospray Mass-Spectrometry – Ultra-High-Sensitivity Analysis of Peptides and Proteins. Journal of the American Society for Mass Spectrometry, 1994. 5(7): p. 605-613. 75. Gatlin, C.L., et al., Protein identification at the low femtomole level from silver-stained gels using a new fritless electrospray interface for liquid chromatography microspray and nanospray mass spectrometry. Analytical Biochemistry, 1998. 263(1): p. 93-101. 76. Aturki, Z., et al., On-line CE-MS using pressurized liquid junction nanoflow electrospray interface and surface-coated capillaries. Electrophoresis, 2006. 27(23): p. 4666-4673. 77. Edwards, J.L., et al., Negative mode sheathless capillary electrophoresis electrospray ionization-mass spectrometry for metabolite analysis of prokaryotes. Journal of Chromatography A, 2006. 1106(1-2): p. 80-88. 78. http://www.kiriama.com/kiriama%20single-mode%20polymer%20fibers_009.htm. 79. Wilm, M.S. and M. Mann, Electrospray and Taylor-Cone Theory, Doles Beam of Macromolecules at Last. International Journal of Mass Spectrometry, 1994. 136(2-3): p. 167-180.

80. Hirt, C.W. and B.D. Nichols, Volume of Fluid (Vof) Method for the Dynamics of Free Boundaries. Journal of Computational Physics, 1981. 39(1): p. 201-225. 81. Melcher, J.R., Continuum electromechanics. 1981, Cambridge, Mass.: MIT Press. 1 v. (various pagings). 82. http://www.flow3d.com/cfd-101/cfd-101-FAVOR.html. 83. http://www.flow3d.com/cfd-101/cfd-101-FAVOR-no-loss.html. 84. Savage, B.M. and M.C. Johnson, Flow over ogee spillway: Physical and numerical model case study. Journal of Hydraulic Engineering-Asce, 2001. 127(8): p. 640-649. 85. http://www.flow3d.com/cfd-101/cfd-101-free-surface-fluid-flow.html. 86. Graham T. T. Gibson, R.D.W.a.R.D.O., Multiple electrosprays generated from a single poly carbonate microstructured fibre. Mass Spectrometry, 2011. 87. Smith, R.D., et al., Analytical characterization of the electrospray ion source in the nanoflow regime. Analytical Chemistry, 2008. 80(17): p. 6573-6579. 88. Hirt, C.W., Electro-hydrodynamics of semi-conductive fluids: with application to electro-spraying. Flow Science Technical Note, 2004. 70(FSI–04–TN70): p. 1-7. 89. de la Mora, J.F., The fluid dynamics of Taylor cones. Annual Review of Fluid Mechanics, 2007. 39: p. 217-243. 90. Cloupeau, M. and B. Prunetfoch, Electrostatic Spraying of Liquids in Cone-Jet Mode.Journal of Electrostatics, 1989. 22(2): p. 135-159. 91. Hayati, I., A.I. Bailey, and T.F. Tadros, Investigations into the Mechanisms of Electrohydrodynamic Spraying of Liquids .1. Effect of Electric-Field and the Environment on Pendant Drops and Factors Affecting the Formation of Stable Jets and Atomization. Journal of Colloid and Interface Science, 1987. 117(1): p. 205-221. 92. FLOW-3D User Manual, Ver. 9.4. 93. Sen, A.K., J. Darabi, and D.R. Knapp, Analysis of Droplet Generation in Electrospray Using a Carbon Fiber Based Microfluidic Emitter. Journal of Fluids Engineering-Transactions of the Asme, 2011. 133(7).

Figure 20. Top: image of electrospray, bottom: cone-jet profile using the CF emitter. Distance between the carbon fiber tip and the counter electrode is 4.0 mm, potential difference is 3500 V, flow rate is 300 nL min−1 .

Modeling and characterization of a carbon fiber emitter for electrospray ionization

A K Sen1, J Darabi1, D R Knapp2 and J Liu2
1 MEMS and Microsystems Laboratory, Department of Mechanical Engineering,
University of South Carolina, 300 Main Street, Columbia, SC 29208, USA
2 Department of Pharmacology, Medical University of South Carolina, 173 Ashley Avenue,
Charleston, SC 29425, USA
E-mail: darabi@engr.sc.edu

뾰족한 탄소 섬유(CF)를 사용하는 새로운 마이크로 스케일 이미터는 질량 분석 (MS) 분석에서 전기 분무에 사용할 수 있습니다. 탄소 섬유는 360 µm OD 및 75 µm ID의 용융 실리카 모세관과 동축에 위치하며 날카로운 팁은 튜브 말단에서 30 µm 연장됩니다.

Abstract

전기 분무 이온화 (ESI) 프로세스는 전기 유체 역학을 해결하기 위한 Taylor–Melcher 누설 유전체 유체 모델 및 액체-가스 인터페이스 추적을 위한 유체 부피 (VOF) 접근 방식을 기반으로 하는 전산 유체 역학 (CFD) 코드를 사용하여 시뮬레이션 됩니다. CFD 코드는 먼저 기존 지오메트리에 대해 검증한 다음 CF 이미터 기반 ESI 모델을 시뮬레이션하는데 사용됩니다.

시뮬레이션된 전류 흐름 및 전류 전압 결과는 CF 이미터의 실험 결과와 잘 일치합니다. 이미터 형상, 전위차, 유속 및 액체의 물리적 특성이 CF 이미터의 전기 분무 거동에 미치는 영향을 철저히 조사합니다.

스프레이 전류와 제트 직경은 액체의 유속, 전위차 및 물리적 특성과 상관 관계가 있으며 상관 결과는 문헌에 보고된 결과와 정량적으로 비교됩니다. (이 기사의 일부 그림은 전자 버전에서만 색상입니다)

Introduction

1980 년대 후반부터 매트릭스 보조 레이저 탈착 이온화 (MALDI)와 전기 분무 이온화 (ESI)의 두 가지 이온화 기술을 구현하여 감도, 속도 및 구조 정보 수준 측면에서 MS 분석이 엄청나게 성장했습니다. 1980 년대 초까지 전자 충격 (EI) 또는 화학 이온화 (CI) 방법은 가스 크로마토 그래피에 적합한 작은 생체 분자를 이온화 하는 데 사용되었습니다.

그러나 크고 열에 민감한 비 휘발성 샘플은 적절한 사전 처리 없이 EI 또는 CI-MS 기술로 분석 할 수 없습니다 [1]. ESI 기술을 사용하면 액체상에서 직접 이러한 큰 분자를 분석 할 수 있습니다 [2]. Zeleny [3, 4]는 출구에 높은 전위를 적용하여 모세관에서 액체 용액을 분사 할 수 있음을 보여주었습니다.

Dole [5, 6] 및 Fenn [7]의 선구적인 연구는 ESI를 고분자 및 생체 분자와 같은 대형 화합물의 이온화 방법으로 표시했습니다. 이에 이어이 기술에 의한 기상 이온 발생에 관련된 과정과 메커니즘이 널리 조사되고 있습니다.

ESI 방법에서 기체 이온화 된 분자는 강한 전계가 있는 상태에서 미세한 물방울을 생성하여 액체 용액에서 생성됩니다. ESI 프로세스의 이러한 능력은 단백질 및 기타 생체 분자 연구에 자연적으로 적용됨을 발견했습니다. ESI 방법과 관련된 다양한 프로세스가 그림 1에 나와 있습니다.

Figure 1. Schematic of an ESI process.
Figure 1. Schematic of an ESI process.

ESI 전위는 일반적으로 전도성 물질로 코팅 된 이미 터 튜브를 통해 외부에서 샘플 액체에 적용되지만 액체 샘플 내부에 적용될 수도 있습니다. Herring과 Qin [8]은 이미 터 팁에 삽입된 팔라듐 와이어를 통해 전기 분무 전위가 적용되는 모세관 전기 영동 (CE)을위한 ESI 인터페이스를 보여주었습니다.

Chiou의 설계 [9]에서는 작은 PDMS 칩에 있는 샘플 저장소, 마이크로 채널 및 실리카 모세관 노즐과 통합 된 내장 전극을 통해 전기 분무를 위한 고전압이 적용되었습니다.

Cao and Moini [10]는 ESI 전압이 모세관 내부에 위치한 전극을 통해인가되고 전기적 접촉이 출구 근처 모세관 벽의 작은 구멍을 통해 유지되는 전기 분무 방출기를 설계했습니다. 작은 모세관 직경 (~ 10 µm)을 가진 이미 터를 사용하여 낮은 전압에서 전기 분무가 가능하지만, 더 작은 구멍은 과도한 배압으로 인해 쉽게 막힐 수 있습니다.

직경이 더 큰 (> 50µm) 이미 터를 처리하는 것이 더 쉽습니다. 그러나 그들은 더 작은 직경의 이미 터만큼 효율적이지 않습니다 [11]. 일반적으로 ESI 전압을 적용하기 위해 유리 또는 용융 실리카와 같은 절연 재료로 제작 된 저 유량 이미 터의 외주에 전도성 코팅이 적용됩니다.

용융 실리카 모세관의 끝 부분에있는 스퍼터 코팅 된 귀금속 층은 내구성에 빠르게 영향을 미치는 것으로 관찰되었습니다. 코팅의 빠른 열화는 방전, 전기 화학적 반응 및 층과 용융 실리카 표면 사이의 불량한 기계적 결합으로 인해 발생할 수 있습니다.

이러한 에미 터의 수명은 스퍼터 코팅 후에 금을 전기 도금하거나 [12] 스퍼터 코팅 된 금 위에 SiOx를 코팅하여 증가시킬 수 있습니다 [13]. 크롬 또는 니켈 합금의 접착층 위에 금으로 코팅 된 이미 터는 우수한 결합력을 제공 할 수 있으며 음극으로 작동 할 때 내구성이 있습니다.

그러나 양극으로 작동하는 동안 접착층은 금 막을 통해 화학적으로 용해됩니다. 이미 터의 안정성과 내구성을 향상시키기 위해 대체 전도성 코팅이 평가되었습니다.

안정적인 ESI 작동을 위해 콜로이드 흑연 코팅 이미 터가 사용되었으며 수명이 길었습니다 [14]. 폴리아닐린 (PANI) 코팅 이미 터는 두꺼운 코팅으로 인해 높은 내구성을 보여주고 방전에 강합니다. PANIcoated와 gold-coated nanospray emitter의 electrospray ionization 거동을 비교 한 결과 PANIcoated emitter는 goldcoated emitter와 비슷한 향상된 감도를 제공합니다 [15].

그라파이트-폴리이 미드 혼합물은 또한 무 접착 전기 분무 방출기의 경우 전도성 코팅으로 사용되었습니다. 전도성 코팅의 안정성은 산화 스트레스 동안 좋은 성능을 나타내는 전기 화학적 방법에 의해 조사되었습니다 [16].

탄소 코팅 이미 터의 기능은 마이크로 스프레이 및 시스리스 CE 및 ESI 응용 분야에서 입증되었습니다. 이 이미 터는 견고하지는 않지만 방수가 되지 않는 CE 또는 ESI 애플리케이션에 충분히 내구성이있었습니다 [17].

우리는 막힘 문제를 제거하고 시료 액체와 금층 사이의 접촉 문제를 피할 수있는 뾰족한 탄소 섬유 기반의 새로운 ESI 방출기를 도입하여 ESI 시스템의 적용 성, 신뢰성 및 내구성을 향상 시켰습니다 [18]. 이 작업에서 탄소 섬유 기반 ESI 이미 터는 전산 유체 역학 (CFD) 소프트웨어 패키지 FLOW-3D [19]를 사용하여 시뮬레이션됩니다.

실험은 새로운 CF 이미 터를 사용하여 수행됩니다. 모델 예측은 실험 결과와 비교됩니다. 새로운 이미 터의 ESI 성능은 이미 터의 기하학적 구조, 유속, 액체의 물리적 특성과 같은 다양한 매개 변수에 대한 반응을 연구하여 평가됩니다.

스프레이 전류 및 제트 직경은 유량 및 액체의 특성과 상관 관계가 있으며 상관 결과는 문헌에보고 된 결과와 정량적으로 비교됩니다. 다음 섹션에서 ESI 공정을 지배하는 전기 유체 역학 이론은 Taylor–Melcher 누설 유전체 모델 [20]을 참조하여 설명됩니다.

그런 다음 Hartman 등이 사용하는 ESI 구성을 고려하여 CFD 코드의 유효성을 확인합니다 [21]. 또한 CF 기반 ESI 모델에 대한 시뮬레이션 및 실험 결과가 제시되고 논의됩니다. 마지막으로 모수 연구 결과와 상관 관계를 제시하고 논의합니다.

Figure 2. Forces in the liquid cone.
Figure 2. Forces in the liquid cone.
Figure 3. Schematic of the ESI model studied by Hartman et al [21].
Figure 3. Schematic of the ESI model studied by Hartman et al [21].
Figure 6. Cone-Jet profile and the electric potential contours at 19 kV; cone length is 4.3 mm.
Figure 6. Cone-Jet profile and the electric potential contours at 19 kV; cone length is 4.3 mm.
Figure 7. A photograph of the experimental cone shape; cone length is 4.2 ± 0.2 mm [21].
Figure 7. A photograph of the experimental cone shape; cone length is 4.2 ± 0.2 mm [21].
Figure 15. Electric field contours at various time steps
Figure 15. Electric field contours at various time steps
Figure 20. Top: image of electrospray, bottom: cone-jet profile using the CF emitter. Distance between the carbon fiber tip and the counter electrode is 4.0 mm, potential difference is 3500 V, flow rate is 300 nL min−1 .
Figure 20. Top: image of electrospray, bottom: cone-jet profile using the CF emitter. Distance between the carbon fiber tip and the counter electrode is 4.0 mm, potential difference is 3500 V, flow rate is 300 nL min−1 .

References

[1] Siuzdak M 1996 Mass Spectrometry for Biotechnology (New York: Academic)
[2] Cole R B (ed) 1997 Electrospray Ionization Mass Spectrometry (New York: Wiley-Interscience)
[3] Zeleny J 1914 Phys. Rev. 3 69–91
[4] Zeleny J 1917 Phys. Rev. 10 1–6
[5] Dole M, Mack L L, Hines R L, Mobley R C, Ferguson L D and Alice M B 1968 Molecular beams of macroions
J. Chem. Phys. 49 2240–9
[6] Clegg G A and Dole M 1971 Molecular beams of macroions: III. Zein and polyvinylpyrrolidone Biopolymers
10 821–6
[7] Fenn J B, Mann M, Meng C K, Wong S F and Whitehouse C M 1989 Electrospray ionization for mass
spectrometry of large biomolecules Science 246 64–71
[8] Herring C J and Qin J 1999 An on-line preconcentrator and the evaluation of electrospray interfaces for the capillary
electrophoresis/mass spectrometry of peptides Rapid Commun. Mass Spectr. 13 1–7
[9] Chiou C H, Lee G B, Hsu H T, Chen P W and Liao P C B 2002 Microscale Tools for Sample Preparation, Separation
and Detection of Neuropeptides Sensors Actuators B 86 280–6
[10] Cao P and Moini M 1997 A novel sheathless interface for capillary electrophoresis/electrospray ionization mass
spectrometry using an in-capillary electrode J. Am. Soc. Mass Spectrom 8 561–4
[11] Janini G M, Conards T P, Wilkens K L, Issaq H J and Veenstra T D 2003 A sheathless nanoflow electrospray
interface for on-line capillary electrophoresis mass spectrometry Anal. Chem 75 1615–9
[12] Barroso M B de Jong and Ad P 1999 Sheathless preconcentration-capillary zone electrophoresis-mass
spectrometry applied to peptide analysis J. Am. Soc. Mass Spectrom 10 1271–8
[13] Valaskovic G A and McLafferty F W 1996 Long-lived metallized tips for nanoliter electrospray mass spectrometry
J. Am. Soc. Mass Spectrom. 7 1270–2
[14] Zhu X, Thiam S, Valle B C and Warner I M 2002 A colloidal graphite coated emitter for seathless capillary
electrophoresis/nanoelectrospray ionization mass spectrometry Anal. Chem 74 5405–9
[15] Maziarz E P I II, Lorenz S A, White T P and Wood T D 2000 Polyaniline: a conductive polymer coating for durable
nanospray emitters J. Am. Soc. Mass. Spectrom 11 659–63
[16] Nilsson S, Wetterhall M, Bergquist J, Nyholm L and Markides K E 2001 A simple and robust conductive
graphite coating for sheathless electrospray emitters used in capillary electrophoresis/mass spectrometry Rapid
Commun. Mass Spectr. 15 1997–2000
[17] Chang Y Z and Her G R 2000 Sheathless capillary electrophoresis/electospray mass spectrometry using a
carbon-coated tapered fused silica capillary with a beveled edge Anal. Chem. 72 626–30
[18] Liu J, Ro K W, Busman M and Knapp D R 2004 Electrospray ionization with a pointed carbon fiber emitter Anal. Chem. 76 3599–606
[19] Hirt C W 2004 Electro-hydrodynamics of semi–conductive fluids: with application to electro–spraying Flow Science
Technical Note 70 FSI–04–TN70 1–7
[20] Saville D A 1997 Electrohydrodynamcis: the Taylor–Melcher leaky dielectric model Annu. Rev. Fluid Mech. 29 27–64
[21] Hartman R P A, Brunner D J, Camelot D M A, Marijnissen J C M and Scarlett B 1999
Electrohydrodynamic atomization in the cone-jet mode physical modeling of the liquid cone and jet J. Aerosol Sci.
30 823–49
[22] Castellanos A 1998 Basic Concepts and Equations in Electrohydrodynamics Electrohydrodynamics
ed A Castellanos (Berlin: Springer)
[23] Melcher J R 1981 Continuum Electromechanics (Cambridge, MA: MIT Press)
[24] Hirt C W and Nichols B D 1981 Volume of fluid (VOF) method for the dynamics of free boundaries J. Comp. Phys.
39 201–25
[25] De la Mora F J and Loscertales I G 1994 The current emitted by highly conducting Taylor cones J. Fluid Mech. 260
155–84
[26] Ganan-Calvo A M 1997 Cone–jet analytical extension of Taylor’s electrostatic solution and the asymptotic universal
scaling laws in electrospraying Phys. Rev. Lett. 79 217–20
[27] Higuera F J 2004 Current/flow–rate characteristic of an electrospray with a small meniscus J. Fluid Mech.
513 239–46
[28] Zeng J, Sobek D and Korsmeyer T Electro-hydrodynamic modeling of electrospray ionization: cad for a microfluidic
device-mass spectrometer interface Transducers ’03: 12th Int. Conf. on Solid State Sensors, Actuators and
Microsystems 2 1275–8
[29] Ganan–Calvo A M, Davila J and Barrero A 1997 Current and droplet size in the electrospraying of liquids. Scaling laws J. Aerosol Sci. 28 249–75
[30] Cloupeau M and Prunet-Foch B 1989 Electrostatic spraying of liquids in cone–jet mode J. Electrost. 22 135–59

igure 1:Essential componentsof the MHD printhead (a) cross-sectional view of printhead showing flow of liquid metal.(b) simulation model showing the magneticfield generated by a pulsed magnetic coil as well as an ejecteddroplet of liquid aluminum.

Timeline of molten metal droplet ejection

용융 금속 액적 분출 타임 라인

Keywords: Magnetohydrodynamicdroplet ejection, droplet on demandprinting, 3D printing of molten metal, additive manufacturing, thermo-fluidic analysis, molten aluminum.

우리는 액체 금속 방울을 사용하여 3D 고체 금속 구조의 DOD (drop-on-demand) 프린팅을 위한 새로운 방법을 제시합니다. 이 방법은 MHD (Magnetohydrodynamic) 기반 방울 생성에 의존합니다. 특히, 외부 코일에 의해 공급되는 맥동 자기장은 액체 금속으로 채워진 분사 챔버 내에서 MHD 기반 힘 밀도를 유도하여 물방울이 노즐을 통해 분사되도록 합니다.

임의의 모양의 3 차원 (3D) 고체 금속 구조는 드롭 방식의 유착 및 응고와 함께 방울의 층별 패턴 증착을 통해 인쇄 할 수 있습니다. 샘플 인쇄 구조와 함께 이 프로토 타입 MHD 인쇄 시스템을 소개합니다. 또한 드롭 생성을 제어하는 기본 물리학에 대해 논의하고 장치 성능을 예측하기 위한 계산 모델을 소개합니다.

주문형 잉크젯 인쇄는 상업용 및 소비자용 이미지 재생을 위한 잘 정립된 방법입니다. 이 기술을 주도하는 동일한 원리는 기능 인쇄 및 적층 제조 분야에도 적용될 수 있습니다. 기존의 잉크젯 기술은 폴리머에서 살아있는 세포에 이르기까지 다양한 재료를 증착하고 패턴 화하여 다양한 기능성 매체, 조직 및 장치를 인쇄하는 데 사용되어 왔습니다 [1, 2]. 이 작업의 초점은 잉크젯 기반 기술을 3D 솔리드 금속 구조 프린팅으로 확장하는 데 있습니다 [3, 4]. 현재 대부분의 3D 금속 프린팅 응용 분야에는 레이저 (예 : 선택적 레이저 소결 [5] 및 직접 레이저 금속 소결 [6]) 또는 전자 빔 (예 : 레이저 소결 [6])과 같은 외부 지향 에너지 원의 영향으로 증착 된 금속 분말 소결 또는 용융이 포함됩니다. 전자빔 용융 [7])을 사용하여 고체 물체를 형성합니다. 그러나 이러한 방법은 비용 및 복잡성 측면에서 단점이 있습니다. 3D 프린팅 프로세스에 앞서 금속을 밀링해야 합니다.

igure  1:Essential componentsof the MHD printhead (a) cross-sectional  view  of  printhead  showing  flow  of liquid metal.(b) simulation model showing the magneticfield  generated  by  a  pulsed  magnetic  coil  as  well  as an ejecteddroplet of liquid aluminum.
igure 1:Essential componentsof the MHD printhead (a) cross-sectional view of printhead showing flow of liquid metal.(b) simulation model showing the magneticfield generated by a pulsed magnetic coil as well as an ejecteddroplet of liquid aluminum.

이 작업에서 우리는 자기 유체 역학의 원리에 기반한 금속 구조물의 적층 제조에 대한 새로운 접근 방식을 소개합니다. 이 방법에서는 감긴 고체 금속 와이어가 MHD 프린트 헤드의 아세라 미치 팅 챔버에 연속적으로 공급되고 용융되어 그림 1에 표시된 것처럼 모세관 힘을 통해 배출 챔버에 공급되는 액체 금속 저장소를 형성합니다. 코일이 배출 챔버를 둘러싸고 있습니다. 액체 금속 내에서 과도 전기장을 유도하는 과도 자기장을 생성하도록 전기적으로 펄스됩니다. 전기장은 유도 된 순환 전류 밀도를 생성하며, 이는 적용된 자기장과 결합하여 챔버 내에서 오리피스의 액체 금속 방울을 방출하는 역할을하는 로렌츠 힘 밀도 (fMHD)를 생성합니다. 분출 된 액 적은 기질로 이동하여 결합 및 응고되어 확장 된 고체 구조를 형성합니다. 임의의 형태의 3 차원 구조는 입사 액 적의 정확한 패턴 증착을 가능하게하는 움직이는 기판을 사용하여 층별로 인쇄 될 수 있습니다. 이 기술은 Vader Systems (www.vadersystems.com)에서 MagnetoJet이라는 상표명으로 개척하고 상용화했습니다. MagnetoJet 인쇄 공정의 장점은 상대적으로 높은 증착 속도와 낮은 재료 비용으로 임의의 모양의 3D 금속 구조를 인쇄하는 것입니다. 이 작업에서는 MagnetoJet 프로토 타입 프린팅 프로세스에 대해 논의하고 샘플 3D 프린팅 구조를 시연하며 합리적인 설계 및 장치 성능 예측을 가능하게하는 계산 모델을 소개합니다.

Figure 2:Printed     3D structures: (a) ring showing as printed base and processed    upper    portion, and (b) cat
Figure 2:Printed 3D structures: (a) ring showing as printed base and processed upper portion, and (b) cat
Figure 2: Computational domain and boundary conditions for the two-dimensional flow problem

3-D transient simulation of viscoelastic coating flows

점탄성 코팅 흐름의 3-D 과도 시뮬레이션

James M. Brethour
Flow Science, Inc.
Santa Fe, New Mexico USA 87505
Presented at the 13th International Coating Science and Technology Symposium, September 10-
13, 2006, Denver, Colorado1

일시적인 프로세스의 3 차원 시뮬레이션은 자유 표면 이동 중에 왜곡을 방지하기 위해 시뮬레이션 중에 업데이트 해야 하는 복잡한 메시를 생성하기 때문에 일반적으로 사용자와 컴퓨터 모두에게 매우 어렵고 지루합니다.

고정된 규칙적인 메시를 통해 유체 운동을 추적하는 Eulerian 기술을 사용하면 이러한 어려움이 제거됩니다. 이러한 방식으로, 큰 유체 변형과 심지어 분열을 계산할 수 있습니다.

이 작업에 사용된 계산 소프트웨어인 FLOW-3D® [1]는 지속적으로 변화하는 유체 영역의 자유 표면을 추적하기 위해 Volume-of-Fluid 기반 기술의 독창적이고 진정한 형태 인 TruVOF®를 사용합니다.

이 모델에 추가 된 것은 점탄성 흐름의 시뮬레이션을 가능하게 하는 사용자 정의입니다. 점탄성 모델은 형태 텐서 [2]를 사용하여 각 유체 요소의 변형 및 회전 이력을 추적합니다. 이러한 계산은 이미 흐름 모델에 존재하는 질량 보존 및 운동량 방정식과 함께 해결됩니다. 필요한 추가 매개 변수는 탄성 계수와 이완 시간입니다.

계산 결과는 슬롯 코팅 [3]에서 하류 접촉 라인이 불안정해질 때까지 코팅액의 공급이 점차 감소하는 저 유량 한계의 실험 결과와 비교됩니다. 계산 결과는 모세관 수의 변화와 유체의 탄성 모두에 대한 실험과 잘 연관되어 있습니다.

Figure 1: Two-dimensional slice of slot coating process; in the experiments, the coating gap was maintained at 100 μm, the slot gap was 125 μm, and the vacuum pressure and web speed were continously varied.
Figure 1: Two-dimensional slice of slot coating process; in the experiments, the coating gap was maintained at 100 μm, the slot gap was 125 μm, and the vacuum pressure and web speed were continously varied.
Figure 2: Computational domain and boundary conditions for the two-dimensional flow problem
Figure 2: Computational domain and boundary conditions for the two-dimensional flow problem
Figure 3: Plot of low flow limits in slot coating as a function of capillary number and fluid elasticity. The solid markers indicate simulation results while the open markers indicate experimental results [3]. The lines represent best-fit power-law curves.
Figure 3: Plot of low flow limits in slot coating as a function of capillary number and fluid elasticity. The solid markers indicate simulation results while the open markers indicate experimental results [3]. The lines represent best-fit power-law curves.
Damascene templates

High-Rate Nanoscale Offset Printing Process Using Directed Assembly and Transfer of Nanomaterials

지난 10 년 동안 나노 크기의 재료와 공정을 제품에 통합하는 데 제한적인 성공을 거두면서 나노 기술에 상당한 투자와 발전이 있었습니다.

잉크젯, 그라비아, 스크린 프린팅과 같은 접근 방식은 나노 물질을 사용하여 구조와 장치를 만드는 데 사용됩니다. [1–7] 그러나 상당히 느리고 µm 스케일 분해능 만 제공 할 수 있습니다. 다양한 모양과 크기의 100nm 미만의 특징을 달성하기 위해 딥펜 리소그래피 (DPN) [8-11] 및 소프트 리소그래피 [12-16]와 같은 다양한 기술이 개발되고 광범위하게 연구되었습니다.

DPN은 직접 쓰기 기술로, atomic force microscopy 현미경 팁을 사용하여 다양한 기판에 여러 패턴을 생성합니다. DPN을 사용한 확장 성을 해결하기 위해 단일 AFM 팁 대신 2D 형식으로 배포 된 AFM (Atomic Force Microscopy) 팁 [17,18]이 사용되었습니다. 소프트 리소그래피에서는 나노 물질을 포함하는 잉크로 적셔진 원하는 릴리프 패턴을 가진 경화된 엘라스토머가 기판과 컨 포멀 접촉하게 되며, 여기서 패턴 화 된 나노 물질이 전달되어 기판에서 원하는 특징을 달성합니다.

이 논문에서는 작거나 큰 영역에서 몇 분 만에 나노, 마이크로 또는 거시적 구조를 인쇄 할 수 있는 다중 스케일 오프셋 인쇄 접근 방식을 제시합니다. 이 프로세스는 나노 입자 (NP), 탄소 나노 튜브 (CNT) 또는 용해 된 폴리머를 포함하는 서스펜션 (잉크)에서 나노 물질의 전기 영동 방향 조립을 사용하여 특별히 제작 된 재사용 가능한 Damascene 템플릿에 패턴을 “inking” 하는 것으로 시작됩니다. 이 잉크 프로세스는 실온과 압력에서 수행됩니다.

두 번째 단계는 템플릿에 조립된 나노 물질이 다른 기판으로 전송되는 “printing”로 구성됩니다. 전송 프로세스가 끝나면 템플릿은 다음 조립 및 전송주기에서 즉시 재사용 할 수 있습니다. 이 오프셋 인쇄 프로세스를 통해 NP (폴리스티렌 라텍스 (PSL), 실리카,은) 및 CNT (다중 벽 및 단일 벽)를 100μm에서 500nm까지의 크기 범위를 가진 패턴에 조립하고 유동성 기판에 성공적으로 옮깁니다.

다양한 나노 물질을 다양한 아키텍처로 조립하기 위해 템플릿 유도 유동, 대류, 유전 영동 (DEP) 및 전기 영동 조립과 같은 몇 가지 직접 조립 프로세스가 조사되었습니다. 모세관력이 지배적인 조립 메커니즘인 유체 조립 공정은 다양한 나노 물질에 적용 할 수 있습니다.

대류 조립 공정은 현탁 메니 스커 스와 증발을 활용하여 단일 나노 입자 분해능으로 정밀 조립을 가능하게 합니다. 이러한 조립 공정 중 많은 부분이 트렌치와 같은 마이크로 및 나노 스케일 기능으로 고해상도의 직접 조립을 보여 주었지만, 확장성 부족, 느린 공정 속도 및 반복성과 같은 많은 단점이 있습니다.

DEP 어셈블리는 NP와 전극 사이에 고배향 탄소 나노 튜브 어셈블리를 사용하여 나노 와이어 및 구조를 만드는 데 사용되었습니다. 조립 효율은 전기장과 전기장 구배에 상당한 영향을 미치는 전극의 기하학적 구조와 간격에 크게 좌우됩니다. 전기 영동 기반 조립 공정은 유체 조립에 비해 훨씬 짧은 시간에 전도성 표면에 표면 전하를 가진 나노 물질을 조립하는 것을 포함합니다. [34–37]

그러나 전기 영동 조립은 조립이 전도성 표면에 발생해야 하므로 다양한 장치를 만드는 데 실용적이지 않습니다. 한 가지 해결책은 원하는 나노 스케일 구조를 기반으로 전도성 패턴이 있는 템플릿을 만들고, 전기 영동 공정을 사용하여 패턴 위에 나노 물질을 조립 한 다음 조립 된 구조를 수용 기판에 옮기는 것입니다.

그림 1a와 같이 절연 필름에 전도성 와이어와 같은 패턴 구조가있는 기존 템플릿을 사용하면 나노 스케일 와이어의 잠재적 인 큰 강하로 인해 어셈블리가 불균일 해지며 대부분의 입자는 그림 1에 표시된 마이크로 와이어 b. 또한 NP는 3D 와이어의 측벽에도 조립되므로 바람직하지 않습니다. 또한 나노 스케일 와이어와 템플릿 사이의 작은 접촉 면적으로 인해 나노 스케일 와이어는 이송 과정에서 쉽게 벗겨집니다.

Damascene templates
Figure 1. Damascene templates: a) A schematic of a conventional wire template used for electrophoretic assembly. In these templates nanowire are connected to a micrometer scale electrodes, which are in turn connected, to a large metal pad through which the potential is applied. b) SEM images of a typical nanoparticle assembly result obtained for confi guration shown in (a). c) A schematic of a Damascene template where all of the wires (nano- or micrometer scale) and the metal pad are connected to a conductive fi lm underneath the insulating fi lm. d) A schematic of Damascene template fabrication. Inset is artifi cially colored cross-sectional SEM image showing the metal nanowires to be at the same height as that of the SiO 2 and showing the conductive fi lm underneath the insulator. e) An optical image of a 3 inch Damascene template.
Offset printing
Figure 2. Offset printing: a) A schematic of the nanoscale offset printing approach. The insulating (SiO 2 ) surface of the Damascene template is selectively coated with a hydrophobic SAM (OTS). Using electrophoresis, nanomaterials are assembled on the conductive patterns of the Damascene template (“inking”), which are then transferred to a recipient substrate (“printing”). After the transfer, the template is ready for the next assembly and transfer cycle. b) SEM image of 50 nm PSL particles assembly with high density on 1 µm wide electrodes. c) Silica particles (20 nm) assembly on crossbar 2D patterns demonstrating the versatility of the Damascene template. Inset fi gure is a high-resolution image of assembled silica particles. d) SEM image of assembled SWNTs on micrometer scale patterns. e) MWNTs assembled on 100 µm features. f) Cellulose assembled on 2 µm electrodes. g) SWNTs assembled in cross bar architecture patterns. h) Flexible devices with array of transferred SWNTs and metal electrodes (printed on PEN). Inset is the microscopy image of two electropads and transferred SWNTs on PEN fi lm.
Analysis of nanomaterial assembly on electrodes
Figure 3. Analysis of nanomaterial assembly on electrodes

이것은 또한 그림 3b에 표시된대로 유한 체적 모델링 (Flow 3D)을 사용하는 전기장 윤곽 시뮬레이션 결과에 의해 확인됩니다. 전기장 강도의 윤곽은 전도성 패턴의 가장자리에있는 전기장이 중앙에있는 것보다 더 강하다는 것을 나타냅니다. 그러나 적용된 전위가 2.5V로 증가하면 그림 3c에 표시된대로 100nm 실리카 입자가 Damascene 템플릿을 가로 질러 전도성 패턴의 표면에 완전히 조립되어 조립을위한 임계 전기장 강도에 도달했음을 나타냅니다. 정렬 된 SWNT는 여과 전달 경로를 피하고 나노 튜브 사이의 접합 저항을 최소화하여 소자 성능의 최소 변화를 가져 오기 때문에 많은 응용 분야에서 고도로 조직화 된 SWNT가 필요합니다.

References

[1] M.Abulikemu, E.H.Da’as, H.Haverinen, D.Cha, M.A.Malik, G.E.Jabbour, Angew.Chem.Int.Ed.2014, 53, 599.
[2] a) Z.Lu, M.Layani, X.Zhao, L.P.Tan, T.Sun, S.Fan, Q.Yan, S.Magdassi, H.H.Hng, Small 2014, 10, 3551; b) H.Ko, J.Lee, Y.Kim, B.Lee, C.H.Jung, J.H.Choi, O.S.Kwon, K.Shin, Adv.Mater.2014, 26, 2286.
[3] C.J.Hansen, R.Saksena, D.B.Kolesky, J.J.Vericella, S.J.Kranz, G.P.Muldowney, K.T.Christensen, J.A.Lewis, Adv.Mater.2013, 25, 2.
[4] F.C.Krebs, N.Espinosa, M.Hösel, R.R.Søndergaard, M.Jørgensen, Adv.Mater.2014, 26, 29.
[5] W.Honda, S.Harada, T.Arie, S.Akita, K.Takei, Adv.Funct.Mater. 2014, 24, 3298.
[6] R.Guo, Y.Yu, Z.Xie, X.Liu, X.Zhou, Y.Gao, Z.Liu, F.Zhou, Y.Yang, Z.Zheng, Adv.Mater.2013, 25, 3343.
[7] A.Dzwilewski, T.Wågberg, L.Edman, J.Am.Chem.Soc.2009, 131, 4006.
[8] R.D.Piner, J.Zhu, F.Xu, S.Hong, C.A.Mirkin, Science 1999, 283, 661.
[9] J.-H.Lim, C.A.Mirkin, Adv.Mater.2002, 14, 1474.
[10] X.Liu, L.Fu, S.Hong, V.P.Dravid, C.A.Mirkin, Adv.Mater.2002,14, 231.
[11] D.A.Weinberger, S.Hong, C.A.Mirkin, B.W.Wessels, T.B.Higgins, Adv.Mater.2000, 12, 1600.
[12] J.P.Rolland, E.C.Hagberg, G.M.Denison, K.R.Carter, J.M.DeSimone, Angew.Chem.2004, 116, 5920.
[13] T.Granlund, T.Nyberg, L.S.Roman, M.Svensson, O.Inganäs, Adv.Mater.2000, 12, 269.
[14] Y.Xia, G.M.Whitesides, Annu.Rev.Mater.Sci.1998, 28, 153.
[15] W.S.Beh, I.T.Kim, D.Qin, Y.Xia, G.M.Whitesides.Adv.Mater. 1999, 11, 1038.
[16] Y.Yin, B.Gates, Y.Xia.Adv.Mater.2000, 12, 1426.
[17] K.Salaita, Y.Wang, J.Fragala, R.A.Vega, C.Liu, C.A.Mirkin,Angew.Chem.2006, 118, 7378.
[18] D.Bullen, S.-W.Chung, X.Wang, J.Zou, C.A.Mirkin, C.Liu, Appl.Phys.Lett.2004, 84, 789.
[19] Y.L.Kim, H.Y.Jung, S.Park, B.Li, F.Liu, J.Hao, Y.-K.Kwon, Y.J.Jung, S.Kar, Nat.Photonics 2014, 8, 239.
[20] X.Xiong, L.Jaberansari, M.G.Hahm, A.Busnaina, Y.J.Jung, Small 2007, 3, 2006.
[21] A.B.Marciel, M.Tanyeri, B.D.Wall, J.D.Tovar, C.M.Schroeder, W.L.Wilson, Adv.Mater.2013, 25, 6398.
[22] J.T.Wang, J.Wang, J.J.Han, Small 2011, 7, 1728.
[23] S.Y.Lee, S.H.Kim, H.Hwang, J.Y.Sim, S.M.Yang, Adv.Mater. 2014, 26, 2391.
[24] J.Y.Oh, J.T.Park, H.J.Jang, W.J.Cho, M.S.Islam, Adv.Mater. 2014, 26, 1929.
[25] K.W.Song, R.Costi, V.Bulovi, Adv.Mater.2013, 25, 1420.
[26] P.Maury, M.Escalante, D.N.Reinhoudt, J.Huskens, Adv.Mater. 2005, 17, 2718.
[27] Y.Xia, Y.Yin, Y.Lu, J.McLellan, Adv.Funct.Mater.2003, 13, 907.
[28] L.Jaber-Ansari, M.G.Hahm, S.Somu, Y.E.Sanz, A.Busnaina, Y.J.Jung, J.Am.Chem.Soc.2008, 131, 804.
[29] T.Kraus, L.Malaquin, H.Schmid, W.Riess, N.D.Spencer, H.Wolf,Nat.Nanotechnol.2007, 2, 570.
[30] K.D.Hermanson, S.O.Lumsdon, J.P.Williams, E.W.Kaler, O.D.Velev, Science 2001, 294, 1082.
[31] H.-W.Seo, C.-S.Han, D.-G.Choi, K.-S.Kim, Y.-H.Lee, Microelectron.Eng.2005, 81, 83.
[32] E.M.Freer, O.Grachev, X.Duan, S.Martin, D.P.Stumbo, Nat.Nanotechnol.2010, 5, 525.
[33] D.Xu, A.Subramanian, L.Dong, B.J.Nelson, IEEE Trans.Nanotechnol.2009, 8, 449.
[34] X.Xiong, P.Makaram, A.Busnaina, K.Bakhtari, S.Somu, N.McGruer, J.Park, Appl.Phys.Lett.2006, 89, 193108.
[35] R.C.Bailey, K.J.Stevenson, J.T.Hupp, Adv.Mater.2000, 12, 1930.
[36] Q.Zhang, T.Xu, D.Butterfi eld, M.J.Misner, D.Y.Ryu, T.Emrick, T.P.Russell, Nano Lett.2005, 5, 357.
[37] E.Kumacheva, R.K.Golding, M.Allard, E.H.Sargent, Adv.Mater. 2002, 14, 221.
[38] M.Wei, Z.Tao, X.Xiong, M.Kim, J.Lee, S.Somu, S.Sengupta, A.Busnaina, C.Barry, J.Mead, Macromol.Rapid Commun.2006, 27, 1826.
[39] a) D.Schwartz, S.Steinberg, J.Israelachvili, J.Zasadzinski, Phys.Rev.Lett.1992, 69, 3354; b) W.Yang, P.Thordarson, J.J.Gooding, S.P.Ringer, F.Braet, Nanotechnology 2007, 18, 412001.
[40] S.Siavoshi, C.Yilmaz, S.Somu, T.Musacchio, J.R.Upponi, V.P.Torchilin, A.Busnaina, Langmuir 2011, 27, 7301.
[41] E.Artukovic, M.Kaempgen, D.Hecht, S.Roth, G.Grüner, NanoLett.2005, 5, 757.
[42] L.Hu, D.Hecht, G.Grüner, Nano Lett.2004, 4, 2513.
[43] M.Fuhrer, J.Nygård, L.Shih, M.Forero, Y.G.Yoon, H.J.Choi, J.Ihm, S.G.Louie, A.Zettl, P.L.McEuen, Science 2000, 288,
494.
[44] J.J.Gooding, A.Chou, J.Liu, D.Losic, J.G.Shapter, D.B.Hibbert,Electrochem.Commun.2007, 9, 1677.
[45] A.Chou, T.Böcking, N.K.Singh, J.J.Gooding, Chem.Commun. 2005, 7, 842.
[46] D.Hines, V.Ballarotto, E.Williams, Y.Shao, S.Solin, J.Appl.Phys. 2007, 101, 024503.
[47] H.Park, A.Afzali, S.-J.Han, G.S.Tulevski, A.D.Franklin, J.Tersoff, J.B.Hannon, W.Haensch, Nat.Nanotechnol.2012, 7, 787.
[48] S.Somu, H.Wang, Y.Kim, L.Jaberansari, M.G.Hahm, B.Li, T.Kim, X.Xiong, Y.J.Jung, M.Upmanyu, A.Busnaina, ACS Nano 2010, 4, 4142.
[49] L.Jaber-Ansari, M.G.Hahm, T.H.Kim, S.Somu, A.Busnaina, Y.J.Jung, Appl.Phys.A 2009, 96, 373.
[50] B.Li, M.G.Hahm, Y.L.Kim, H.Y.Jung, S.Kar, Y.J.Jung, ACS Nano 2011, 5, 4826.
[51] B.Li, H.Y.Jung, H.Wang, Y.L.Kim, T.Kim, M.G.Hahm, A.Busnaina, M.Upmanyu, Y.J.Jung, Adv.Funct.Mater.2011, 21, 1810.
[52] M.A.Meitl, Z.T.Zhu, V.Kumar, K.J.Lee, X.Feng, Y.Y.Huang, I.Adesida, R.G.Nuzzo, J.A.Rogers, Nat.Mater.2005, 5, 33.
[53] F.N.Ishikawa, H.Chang, K.Ryu, P.Chen, A.Badmaev, L.GomezDe Arco, G.Shen, C.Zhou, ACS Nano 2008, 3, 73.
[54] N.Inagaki, Plasma Surface Modifi cation and Plasma Polymerization, CRC, Boca Raton, FL, USA 1996.
[55] E.Liston, L.Martinu, M.Wertheimer, J.Adhes.Sci.Technol.1993, 7, 1091.
[56] T.Tsai, C.Lee, N.Tai, W.Tuan, Appl.Phys.Lett.2009, 95, 013107.
[57] J.G.Bai, Z.Z.Zhang, J.N.Calata, G.-Q.Lu, IEEE Trans.Compon.Packag.Technol.2006, 29, 589.
[58] J.G.Toffaletti, Crit.Rev.Clin.Lab.Sci.1991, 28, 253.
[59] J.-L.Vincent, P.Dufaye, J.Berré, M.Leeman, J.-P.Degaute, R.J.Kahn, Crit.Care Med.1983, 11, 449.
[60] R.Henning, M.Weil, F.Weiner, Circ.Shock 1982, 9, 307.

FIG. 2. Sequence of images showing capillary-driven neck evolution and droplet formation for low-viscosity fluids

Computational analysis of self-similar capillary-driven thinning and pinch-off dynamics during dripping using the volume-of-fluid method

낙하 형성 및 분리는 표면 장력 구동 흐름으로 인해 가늘어지는 유체 목의 형성을 포함하여 큰 위상 변화를 수반하며, 목의 pinch-off에서 Laplace pressure와 같은 속성은 유한한 시간 특이성을 나타냅니다. 드롭 형성 중에 발생하는 큰 위상 변형과 비선형성을 정확하게 시뮬레이션하는 것은 일반적으로 pinch-off 순간에 가까운 작은 특징을 해결하기 위해서는 고해상도 및 정확도가 필요하기 때문에 수치 시뮬레이션이 계산적으로 요구됩니다.

필요한 질량 및 계산 시간을 보존하고 인터페이스를 추적하는 데 내재된 이점에도 불구하고, 초기 실무자들이 물 점도가 10배 이상인 유체에 대한 수렴 문제를 보고했기 때문에 낙하 형성 연구에 VOF(Volume-of-fluid) 방법을 활용하는 연구는 거의 없습니다.

이 기여에서, 우리는 FLOW-3D에 구현된 VOF 방법을 사용하여 물 점도보다 4배 더 높은 점도 값을 포함하여 뉴턴 유체에 대한 드리프트의 원형 자유 표면 흐름을 시뮬레이션합니다. 우리는 이 연구의 일부로 수행된 실험에 대해 시뮬레이션된 목 모양, 목 진화 속도 및 헤어짐 길이를 벤치마킹합니다.

핀치오프 역학은 관성, 점성 및 모세관 응력의 복잡한 상호 작용에 의해 결정되며, 여기서 실험과 시뮬레이션 모두에서 대조되는 자기 유사 스케일링 법칙은 종종 역학에 대해 설명합니다. 우리는 시뮬레이션된 반지름 진화 프로파일이 축 대칭 흐름에 대한 뉴턴 유체에 대해 실험적으로 관찰되고 이론적으로 예측되는 핀치오프 역학과 일치함을 보여준다. 또한, 우리는 가는 목 안에서 법칙, 속도 및 변형 필드의 스케일링에 대한 사전 요인을 결정하고, 우리는 실험과 비교할 수 있는 중단 시간과 길이뿐만 아니라 사전 요인을 VOF 방법을 사용하여 시뮬레이션할 수 있음을 보여줍니다.

experimental setup, as shown schematically in Fig. 1(a), includes a dispensing system
experimental setup, as shown schematically in Fig. 1(a), includes a dispensing system
 A numerical simulation of drop formation from a cylindrical nozzle at a constant flow rate is performed. (c) Graphical representation of the VOF approach
A numerical simulation of drop formation from a cylindrical nozzle at a constant flow rate is performed. (c) Graphical representation of the VOF approach
FIG. 2. Sequence of images showing capillary-driven neck evolution and droplet formation for low-viscosity fluids
FIG. 2. Sequence of images showing capillary-driven neck evolution and droplet formation for low-viscosity fluids. (a) A sequence of simulated images of water (0 wt. % glycerol) shows neck formation and subsequent thinning and pinch-off dynamics including the formation of the satellite drop. (b) A sequence of images shows neck radius evolution and drop detachment for the low viscosity fluid composed of 50 wt. % glycerol in water. The time step between images is 500 µs, and the scale bar represents a length of 1 mm for the two cases shown. The color bar shows the velocity field in units of cm/s. The addition of glycerol seems to exercise a relatively minor influence on pinch-off dynamics despite a five-fold increase in viscosity.
FIG. 3. Computed evolution of the minimum radius of the water neck during the drop formation and detachment process
FIG. 3. Computed evolution of the minimum radius of the water neck during the drop formation and detachment process. The instantaneous neck radius of water and the inertio-capillary fit are shown. The inset shows a self-similar nature of neck thinning dynamics close to a pinch-off moment. The characteristic cone angle of 18.1◦ as predicted by Day et al.50 and visualized in experiments52 is captured well using the VOF method.
FIG. 5. Glycerol thinning image sequence and break-up length visualization for three cases
FIG. 5. Glycerol thinning image sequence and break-up length visualization for three cases. (a) Glycerol thinning is shown through a sequence of snapshots with a time step ∆t = 5 ms and reveals quite different dynamics compared to previously seen for low viscosity fluids. The length of a filament changes significantly when the glycerol content increases above 70 wt. %. (b) Final lengths of the simulated liquid filaments before pinch-off for three cases of glycerol + water mixtures (0 wt. %, 70 wt. %, and 100 wt. %).
FIG. 8. Comparison of experiments and simulations for the case of a drop formation for 80 wt. % glycerol and water mixture
FIG. 8. Comparison of experiments and simulations for the case of a drop formation for 80 wt. % glycerol and water mixture. (a) A set of images obtained from experiments (upper row) and simulations (bottom row) with a time step of 1 ms show good agreement. The simulated drop profiles shown in the bottom row are colored by the velocity magnitude [ranging from 0 (dark blue) to 100 cm/s (red) and colored online], and velocity vectors are shown in the images. (b) Radius evolution with time of liquid filament formed during the drop formation process is shown on a log-log plot for the two cases.

FLOW-3D AM

flow3d AM-product
flow3d AM-product

FLOW-3D AM 은 레이저 파우더 베드 융합 (L-PBF), 바인더 제트 및 DED (Directed Energy Deposition)와 같은 적층 제조 공정 ( additive manufacturing )을 시뮬레이션하고 분석하는 CFD 소프트웨어입니다. FLOW-3D AM 의 다중 물리 기능은 공정 매개 변수의 분석 및 최적화를 위해 분말 확산 및 압축, 용융 풀 역학, L-PBF 및 DED에 대한 다공성 형성, 바인더 분사 공정을위한 수지 침투 및 확산에 대한 매우 정확한 시뮬레이션을 제공합니다.

3D 프린팅이라고도하는 적층 제조(additive manufacturing)는 일반적으로 층별 접근 방식을 사용하여, 분말 또는 와이어로 부품을 제조하는 방법입니다. 금속 기반 적층 제조 공정에 대한 관심은 지난 몇 년 동안 시작되었습니다. 오늘날 사용되는 3 대 금속 적층 제조 공정은 PBF (Powder Bed Fusion), DED (Directed Energy Deposition) 및 바인더 제트 ( Binder jetting ) 공정입니다.  FLOW-3D  AM  은 이러한 각 프로세스에 대한 고유 한 시뮬레이션 통찰력을 제공합니다.

파우더 베드 융합 및 직접 에너지 증착 공정에서 레이저 또는 전자 빔을 열원으로 사용할 수 있습니다. 두 경우 모두 PBF용 분말 형태와 DED 공정용 분말 또는 와이어 형태의 금속을 완전히 녹여 융합하여 층별로 부품을 형성합니다. 그러나 바인더 젯팅(Binder jetting)에서는 결합제 역할을 하는 수지가 금속 분말에 선택적으로 증착되어 층별로 부품을 형성합니다. 이러한 부품은 더 나은 치밀화를 달성하기 위해 소결됩니다.

FLOW-3D AM 의 자유 표면 추적 알고리즘과 다중 물리 모델은 이러한 각 프로세스를 높은 정확도로 시뮬레이션 할 수 있습니다. 레이저 파우더 베드 융합 (L-PBF) 공정 모델링 단계는 여기에서 자세히 설명합니다. DED 및 바인더 분사 공정에 대한 몇 가지 개념 증명 시뮬레이션도 표시됩니다.

레이저 파우더 베드 퓨전 (L-PBF)

L-PBF 공정에는 유체 흐름, 열 전달, 표면 장력, 상 변화 및 응고와 같은 복잡한 다중 물리학 현상이 포함되어 공정 및 궁극적으로 빌드 품질에 상당한 영향을 미칩니다. FLOW-3D AM 의 물리적 모델은 질량, 운동량 및 에너지 보존 방정식을 동시에 해결하는 동시에 입자 크기 분포 및 패킹 비율을 고려하여 중규모에서 용융 풀 현상을 시뮬레이션합니다.

FLOW-3D DEM FLOW-3D WELD 는 전체 파우더 베드 융합 공정을 시뮬레이션하는 데 사용됩니다. L-PBF 공정의 다양한 단계는 분말 베드 놓기, 분말 용융 및 응고,이어서 이전에 응고 된 층에 신선한 분말을 놓는 것, 그리고 다시 한번 새 층을 이전 층에 녹이고 융합시키는 것입니다. FLOW-3D AM  은 이러한 각 단계를 시뮬레이션하는 데 사용할 수 있습니다.

파우더 베드 부설 공정

FLOW-3D DEM을 사용하면 아래 동영상처럼 입자의 분포를 무작위로 떨어뜨려 파우더 베드 배치 프로세스를 시뮬레이션할 수 있습니다.

다양한 파우더 베드 압축을 달성하는 한 가지 방법은 베드를 놓는 동안 다양한 입자 크기 분포를 선택하는 것입니다. 아래에서 볼 수 있듯이 세 가지 크기의 입자 크기 분포가 있으며, 이는 가장 높은 압축을 제공하는 Case 2와 함께 다양한 분말 베드 압축을 초래합니다.

파우더 베드 분포 다양한 입자 크기 분포
세 가지 다른 입자 크기 분포를 사용하여 파우더 베드 배치
파우더 베드 압축 결과
세 가지 다른 입자 크기 분포를 사용한 분말 베드 압축

입자-입자 상호 작용, 유체-입자 결합 및 입자 이동 물체 상호 작용은 FLOW-3D DEM을 사용하여 자세히 분석 할 수도 있습니다 . 또한 입자간 힘을 지정하여 분말 살포 응용 분야를 보다 정확하게 연구 할 수도 있습니다.

FLOW-3D AM  시뮬레이션은 이산 요소 방법 (DEM)을 사용하여 역 회전하는 원통형 롤러로 인한 분말 확산을 연구합니다. 비디오 시작 부분에서 빌드 플랫폼이 위로 이동하는 동안 분말 저장소가 아래로 이동합니다. 그 직후, 롤러는 분말 입자 (초기 위치에 따라 색상이 지정됨)를 다음 층이 녹고 구축 될 준비를 위해 구축 플랫폼으로 펼칩니다. 이러한 시뮬레이션은 저장소에서 빌드 플랫폼으로 전송되는 분말 입자의 선호 크기에 대한 추가 통찰력을 제공 할 수 있습니다.

파우더 베드 용해

파우더 베드를 놓은 후 FLOW-3D  WELD 에서 레이저 빔 공정 매개 변수를 지정 하여 고 충실도 용융 풀 시뮬레이션을 수행 할 수 있습니다  . 온도, 속도, 고체 분율, 온도 구배 및 고체 속도의 플롯을 자세히 분석 할 수 있습니다.

레이저 출력 200W, 스캔 속도 3.0m / s, 스폿 반경 100μm에서 파우더 베드의 용융 풀 분석.

용융 풀이 응고되면 FLOW-3D AM  압력 및 온도 데이터를 Abaqus 또는 MSC Nastran과 같은 FEA 도구로 가져와 응력 윤곽 및 변위 프로파일을 분석 할 수도 있습니다.

다층 적층 제조

첫 번째 용융 층이 응고되면 입자의 두 번째 층이 응고 층에 증착됩니다. 새로운 분말 입자 층에 레이저 공정 매개 변수를 지정하여 용융 풀 시뮬레이션을 다시 수행합니다. 이 프로세스를 여러 번 반복하여 연속적으로 응고 된 층 간의 융합, 빌드 내 온도 구배를 평가하는 동시에 다공성 또는 기타 결함의 형성을 모니터링 할 수 있습니다.

다층 적층 적층 제조 시뮬레이션

바인더 분사 (Binder jetting)

Binder jetting 시뮬레이션은 모세관 힘의 영향을받는 파우더 베드에서 바인더의 확산 및 침투에 대한 통찰력을 제공합니다. 공정 매개 변수와 재료 특성은 증착 및 확산 공정에 직접적인 영향을 미칩니다.

방향성 에너지 증착

FLOW-3D AM 의 내장 입자 모델 을 사용하여 직접 에너지 증착 프로세스를 시뮬레이션 할 수 있습니다. 분말 주입 속도와 고체 기질에 입사되는 열유속을 지정함으로써 고체 입자는 용융 풀에 질량, 운동량 및 에너지를 추가 할 수 있습니다. 다음 비디오에서 고체 금속 입자가 용융 풀에 주입되고 기판에서 용융 풀의 후속 응고가 관찰됩니다.

dem9
dem10

FLOW DEM

FLOW DEM 은 FLOW-3D 의 기체 및 액체 유동 해석에 DEM(Discrete Element Method : 개별 요소법)공법인 입자의 거동을 분석해주는 모듈입니다.

주요 기능 :고체 요소의 충돌, 스프링(Spring) / 대시 포트(Dash Pot) 모델 적용 Void, 1 fluid, 2 fluid(자유 계면 포함) 각각의 모드에 대응 가변 밀도 / 가변 직경 입자 크기조절로 입자 특성을 유지하면서 입자 수를 감소 독립적인 DEM의 Sub Time Step 이용

Discrete Element Method : 개별 요소법

다수의 고체 요소의 충돌 운동을 분석하는 데 유용합니다. 유동 해석과 함께 사용하면 광범위한 용도에 응용을 할 수 있습니다.

dem1

입자 간의 충돌

Voigt model은 스프링(Spring) 및 대시 포트(Dash pot)의 조합에 의해 입자 충돌 시의 힘을 평가합니다. 탄성력 부분은 스프링 모델에서,
비탄성 충돌의 에너지 소산부분은 대시 포트 모델에서 시뮬레이션되고 있으며, 중량 및 항력은 작용하는 외력으로 고려 될 수 있습니다.

분석 모드

기본적으로 이용하는 운동 방정식은 FLOW-3D 에 사용되는 질량 입자의 운동 방정식과 같은 것이지만, 여기에 DEM으로
평가되는 항목이 추가되기 형태로되어 있으며, 실제 시뮬레이션으로는 ‘void + DEM’, ‘1 Fluid + DEM’ , ‘ 1 Fluid 자유계면 + DEM ‘을 기본 유동 모드로 취급이 가능합니다.

dem4

입자 유형

입자 타입도 표준 기능의 질량 입자 모델처럼 입자 크기 (반경)와 밀도가 동일한 것 외, 크기는 같지만 밀도가 다른 것이나 밀도는 같지만 크기가 다른 것 등도 취급 가능합니다. 이로 인해 표준 질량 입자 모델에서는 입자 간의 상호 작용이 고려되어 있지 않기 때문에 모든 아래에 가라 앉아 버리고 있었지만, FLOW DEM을 이용하여 기하학적 관계를 평가하는 것이 가능합니다.

dem7

응용 분야

1. Mechanical Engineering 분야

수지 충전, 스쿠류 이송, 분말 이송 / Resin filling, screw conveyance, powder conveyance

2. Civil Engineering분야

3. Civil Engineering 분야

파편, 자갈, 낙 성/ Debris flow, gravel, falling rock

dem11

3. Chemical Engineering, Pharmaceutics 분야

유동층, 사이클론, 교반기 / Fluidized bed, cyclone, stirrer

dem12

4. MEMS, Electrical Engineering 분야

하전 입자를 포함한 전기장 해석 등

dem15

입자 그룹 가시화

그룹 가시화

DEM은 일반적으로 다수의 입자를 필요로하는 분석을 상정하고 있습니다. 
다만 이 경우, 계산 부하가 높아 지므로 현실적인 계산자원을 고려하면, 입자 수가 너무 많아 현실적으로 취급 할 수 없는 경우 입자의 특성은 유지하고 숫자를 줄여 가시화할 필요가 있습니다 .
일반적인 유동해석 계산의 메쉬 해상도에 해당합니다.
메쉬 수 많음 (계산 부하 큼) → 소 (계산 부하 적음)
입자 수 다 (계산 부하 큼) → 소 (계산 부하 적음)

원래 입자수
입자 사이즈를 키운경우
그룹 가시화
  • 입자 수를 줄이기 위해 그대로 입경을 크게했을 경우와 그룹 가시화 한 경우의 비교.
  • 입자 크기를 크게하면 개별 입자 특성이 달라지기 때문에 거동이 달라진다. (본 사례에서는 부력이 커진다.)
  • 그룹 가시화의 경우 개별 특성은 동일 원래의 거동과 대체로 일치한다.

주조 시뮬레이션에 DEM 적용

그룹 가시화 비교 예

그룹 가시화한 경우와 입경을 크게하여 수를 줄인 경우, 입경을 크게하면
개별 입자 특성이 변화하여 거동이 바뀌어 버리기 때문에 실제 계산으로는 사용할 수 어렵습니다.

중자 모래 분사 분석

DEM에서의 계산부하를 생각할 때는 입자모델에 의한 안정제한을 고려해야 하지만 서브타임스텝이라는 개념을 도입함으로써 입자의 경우와 유체의 경우의 타임스텝을 바꾸고 필요이상으로 계산시간을 들이지 않고 효율적으로 계산하는 것을 가능하게 하고 있습니다.

이를 통해 예를 들어 중자사 분사 시뮬레이션 실험에서는 이러한 문제로 자주 이용되는 빙엄 유체에서는 실험과의 정합성이 별로 좋지 않기 때문에 당사에서는 이전부터 입상류 모델이라는 모델을 개발하고 연속체로부터의 접근에서도 실험과의 높은 정합성을 실현할 수 있는 모델화를 해왔는데, 이번에 DEM을 사용해도 그것과 거의 같은 결과를 얻습니다. 할 수 있음을 확인할 수 있었다.

Reference :

  • Lefebvre D., Mackenbrock A., Vidal V., Pavan V. and Haigh PM, 2004,
  • Development and use of simulation in the Design of Blown Cores and Moulds
그림 2 : FLOW-3D를 사용한 흐름 및 형태 시뮬레이션. 파스칼 단위의 압력 및 mm 단위의 거리.

Microscopic Bubbles Switch Fiber-Optic Circuits

Figure 1: The Agilent Photonic Switching Platform
Figure 1: The Agilent Photonic Switching Platform

컴퓨터 시뮬레이션은 광섬유 회로에서 광 신호를 전환하는데 사용되는 혁신적인 스위치에서 미세 기포 문제를 이해하고 해결하는 데 중요한 역할을 했습니다. Agilent Photonic Switching Platform은 평면 광파 회로에서 잘린 작은 트렌치의 올바른 지점에 거품을 불어서 작동합니다. 버블은 광섬유 네트워크를 재구성하기 위해 광선을 다른 경로로 리디렉션 합니다. 초기 프로토타입은 기포 반사로 인해 무언가 불안정하다는 것을 나타내는 성능 문제를 보여주었습니다. 그러나 거품의 크기가 작기 때문에 문제를 진단하고 해결하는데 필요한 포괄적인 물리적 측정을 수행할 수 없었습니다.

애질런트의 선임 과학자인 John Uebbing은 전산 유체 역학 (CFD) 소프트웨어를 사용하여 거품을 시뮬레이션했습니다. 기포는 실리콘 기판에 위치한 전기 히터에 의해 유도된 증발에 의해 유지됩니다. 애질런트 팀은 트렌치 벽의 응결로 인해 유체가 축적된다는 사실을 발견했습니다. 스위치 동작의 대부분을 결정하는 것은 이러한 축적입니다. 추가 시뮬레이션을 통해 연구원들은 안정적인 신호를 제공하기 위해 장치를 변경하는 두 가지 다른 방법을 검증 할 수 있었습니다.  “처음에 우리 팀원 중 일부는 이러한 결과를 믿지 않았지만 계속된 물리적 테스트를 통해 사실이 입증되었습니다.”라고 Uebbing은 말했습니다. “CFD가 없었다면 이 문제의 해결책에 도달하지 못했을 것입니다.”

신기술 개발

광섬유 케이블은 데이터 통신 처리량을 크게 증가 시켰으며, 광 신호 전환을 위한 전기 신호로 전환한 다음 다시 광 신호로 전환하지 않고도 대량의 광섬유 데이터를 전환 할 수 있기를 원했습니다. 1990 년대 중반 Agilent Laboratories (Hewlett-Packard Labs 소속)는 전광 회로 스위치의 중요성을 인식하고 이러한 기술을 개발하기 위한 연구 프로그램을 시작했습니다. 현재 Agilent Labs의 CORL (Communications and Optical Research Laboratory) 내에 엔지니어와 과학자 팀이 구성되어 컴팩트하고 확장 가능하며 광 신호에 최소한의 영향을 미치는 이 고유한 스위치 패브릭을 개발했습니다.

 시뮬레이션은 딤플의 원인을 정확히 파악하는데 도움이 되었으며 여러 대안 솔루션을 식별하고 평가하는 데 도움이되었습니다. 버블 스위치 엔지니어링의 이러한 발전은 FLOW-3D  소프트웨어 에서 사용할 수 있는 고급 모델링 기능 없이는 불가능했을 것  입니다. 우리에게 중요한 것은 프로젝트 시작부터 Flow Science 팀이 입증한 지식과 무결성이었습니다. 우리가 이야기 한 다른 소프트웨어 회사에는 관련된 문제에 대한 표면적 이해만 있는 영업 담당자가 있었지만 Flow Science는 전문 지식을 갖춘 기술 직원을 고용하여 우리가 달성하고자 하는 것을 정확히 이해했습니다. 프로세스의 여러 단계에서 중요한 장애물을 극복 할 수 있는 중요한 도움을 제공했습니다.
– John Uebbing, 애질런트 선임 과학자

작동을 위해 Agilent Photonic Switching Platform은 두 개의 광섬유 네트워크의 교차점에 배치됩니다 (그림 1). 광 신호가 광섬유를 통해 들어 오면 직선 도파관을 통해 방해받지 않고 평면 광파 회로를 통과 할 수 있습니다. 그러나 신호가 다른 광섬유로 리디렉션되어야하는 경우 잉크젯 기술은 두 도파관 경로의 교차점에 거품을 삽입하여 광학 특성을 변경하고 신호를 출력 광섬유로의 경로 아래로 반사합니다. 기포는 거울이나 기계적으로 움직이는 부품을 사용하지 않고도 5 밀리 초 이내에 형성 및 제거 할 수 있습니다. 이 스위치는 교차된 광 도파관 배열과 인덱스 매칭되는 유체에 거품을 불어서 작동합니다. 기포는 소자 기판의 전기 히터에 의해 유도 된 증발에 의해 형성됩니다. 유체는 도파관의 교차점에 위치한 일련의 마이크로 트렌치를 채웁니다. 기포 벽으로부터의 내부 전반사로 인해 빛이 한 도파관에서 다른 도파관으로 전환됩니다. 문제는 광 도파관의 수용 각 또는 개구 수가 상당히 낮다는 것입니다. 기포의 수직 반사벽이 도파관의 축에 수직이 아니면 빛이 출력 도파관으로 제대로 반사되지 않고 신호 손실이 발생합니다.

프로토 타입의 딤플 충격 성능

초기 프로토 타입에서 광범위한 실험 테스트를 수행하여 히터 전력 및 주변 압력이 광학 반사 특성과 기포 모양 및 크기에 미치는 영향을 보여주었습니다. 이 테스트는 반사된 광 신호 대 히터 전력 곡선이 효과적인 광 스위칭에 필요한 엄격한 요구 사항을 충족하지 못하고 반사된 광 신호에 불안정성이 있음을 보여주었습니다.

그림 2 : FLOW-3D를 사용한 흐름 및 형태 시뮬레이션. 파스칼 단위의 압력 및 mm 단위의 거리.
그림 2 : FLOW-3D를 사용한 흐름 및 형태 시뮬레이션. 파스칼 단위의 압력 및 mm 단위의 거리.

컴퓨터 시뮬레이션에서 그림 2와 같이 버블의 각 면에 딤플이 형성되어 있음을 보여 주었을 때, 딤플이 전력 곡선의 혹의 원인이 되었고 반사된 신호가 그렇게 불안정한 이유 일 수 있다는 사실이 애질런트 연구팀에 나타났습니다. 센서로 물리적 측정을 수행하는 팀의 능력은 MEMS 장치의 규모까지 확장되지 않았습니다. 그들이 할 수 있는 최선은 특수 광학 장치를 사용하여 현미경 사진을 찍는 것입니다. 이 사진은 딤플이 파장 스케일에서 매우 얇기 때문에 딤플을 직접 보여줄 수 없습니다.

거품 시뮬레이션

처음에는 버블의 작동을 시뮬레이션하기 위한 여러 가지 대안이 고려되었습니다. 팀은 다양한 분석 모델을 사용하여 기포 형성을 조사했지만 이 모델은 현재 프로토 타입이 좋은 기포를 생성해야 한다고 예측했기 때문에 문제를 포착하기에는 너무 단순했습니다. 맞춤 소프트웨어를 작성하기 위해 대학 교수를 고용했지만 이 프로젝트를 완료하는 데 상당한 시간이 소요되었습니다. 그 동안 Uebbing은 문제의 복잡한 물리학을 처리 할 수 있는 상용 소프트웨어 패키지를 찾기 시작했습니다. “저는 여러 CFD 소프트웨어 개발자들과 이야기를 나눴지만 그들 중 누구도 광범위한 수정 없이 문제를 해결할 수 있는 버블 모델을 가지고 있지 않다고 판단했습니다.”라고 Uebbing은 말했습니다. “반면에 Flow Science는

Flow Science의 새로운 균질 기포 모델은 균일한 기포 압력과 온도를 가정합니다. 이것은 현실에 대한 좋은 근사치입니다. 주요 문제 중 하나는 액체, 증기 및 고체가 모두 결합되는 접점 라인의 모델링입니다. 동질 버블 모델은 이 시점에서 계산 셀의 힘과 플럭스의 균형을 맞춥니 다. Uebbing은 이전 버전의 소프트웨어를 사용하기 시작했지만 새 모델이 출시 되자마자 Uebbing은 문제를 해결해 보았습니다. “시뮬레이션 결과는 결국 실험을 설명하는 데 매우 중요한 dimple 을 보여주었습니다.”라고 Uebbing은 말했습니다. 흥미롭게도 시뮬레이션 결과 버블이 35kHz에서 진동하는 것으로 나타났습니다. 우리는 그것이 실제로 그 주파수에서 진동한다는 것을 보여주는 실험 데이터를 가져 왔지만 우리는 이유를 몰랐습니다.

현실과의 다소 예상치 못한 상관 관계는 팀에게 시뮬레이션 결과에 대한 확신을주었습니다. 시뮬레이션 결과는 문제 영역의 모든 지점에서 유속, 압력 및 온도를 보여줌으로써 테스트에서 측정 할 수 있었던 것 이상이었습니다. 이 결과로 우리는 무슨 일이 일어나고 있는지 파악할 수 있었습니다. dimple은 모세관 현상으로 인해 발생합니다. 응축액이 거품 벽에 쌓입니다. 트렌치 벽에 있는 액체의 얇은 막을 통해 빠져 나 가려고 합니다. 이러한 얇은 층을 통해 액체를 밀어 넣으려면 상당한 압력 차이가 필요합니다. 기포 벽 중앙의 높은 압력으로 인해 기포가 dimple을 형성합니다.”

문제 해결

딤플이 어떻게 형성되었는지 이해하면 안정적인 신호를 제공하기 위해 거품 모양을 수정하는 두 가지 방법이 제안되었습니다. 첫 번째는 트렌치의 유리 측벽 아래로 버블 히터를 확장하는 것입니다. 그런 다음 열이 마이크로 트렌치의 벽 위로 흘러 표면을 건조시킵니다. FLOW-3D를 사용한 시뮬레이션   은 건식 벽 거품이 매우 안정적인 스위치 신호를 제공함을 보여줍니다. 기본 물리학에 따르면 기포 온도가 벽 온도보다 낮 으면 벽이 건조해질 것입니다. 이러한 기대는 FLOW-3D  시뮬레이션 으로 확인되었습니다  .

FLOW-3D로 확인 된 두 번째 방법은 마이크로 트렌치에 소위 정적 버블을 만드는 것입니다. 장치 온도가 압력 설정 저장소 온도보다 약간 더 높으면 정적 거품이 존재합니다. 이 장치 온도는 기포를 트렌치의 모서리로 밀어 넣을 수있는 충분한 압력을 생성하지만 기포가 도파관 어레이와 히터 기판 사이의 틈을 통해 불어 나기에는 충분하지 않습니다. 이러한 정적 기포는 근처의 “crusher”기포를 사용하여 끌 수 있습니다. 이 기포는 일시적으로 충분한 과압을 생성하여 정적 기포가 붕괴되도록합니다. 분쇄기 거품 자체는 더 작은 트렌치에 있으므로 표면 장력이 작업을 완료 한 후 붕괴 될 수 있습니다. FLOW-3D 시뮬레이션은 이 모드에서 스위치 작동을 보여주기 위해 사용되었습니다.

FLOW-3D를 사용 하여 미세 유체 애플리케이션 모델링  의 성능과 다양성에 대해 자세히 알아보십시오. 

collapsed-raised-fluid-column-figure-1-1

Steady-State Accelerator for Free-Surface Flows

자유 표면 흐름을 위한 정상 상태 가속기

이 기사에서 Tony Hirt 박사는 다가오는 FLOW-3D  v12.0 릴리스에서 사용할 수있는 새로운 Steady-State Accelerator에 대해 설명합니다  .

일시적인 흐름의 점근 적 상태를 계산하는 것보다 안정된 자유 표면 흐름을 생성하는 더 빠른 방법이 자주 필요합니다. 상황은 압축성 흐름 솔버를 사용하여 비압축성 흐름을 해결하는 것과 유사합니다. 후자의 경우 압축 파는 붕괴하는 데 오랜 시간이 걸리고 결과적으로 비압축성 흐름을 남길 수 있습니다. 이에 따라 자유 표면 흐름에서 유체는 비압축성이지만 표면 파동은 안정된 자유 표면 구성을 생성하는 데 오랜 시간이 걸릴 수 있습니다.

비압축성 흐름의 경우, 압축 파를 심각하게 감쇠시키는 반복적 인 프로세스 (즉, 압력-속도 반복)를 사용합니다. 물리적으로 반복은 압력과 같은 파동이 국부적 인 영역에 영향을 미치는 짧은 거리를 이동하도록 허용하지만 압력 장에 상당한 노이즈를 유발할 수있는 장거리 전파 및 반사를 피할 수있을만큼 빠르게 감쇠됩니다.

이 노트에서 자유 표면 셀에 적용된 간단한 압력 조정은 표면 교란에 대한 감쇠력으로 작용합니다. 이 댐핑은 안정적인 자유 표면 구성에 대한 접근을 가속화합니다.

Steady-State Accelerator Idea

유체 인터페이스 또는 자유 표면은  VOF (Volume-of-Fluid) 기술을 사용하여 FLOW-3D 에서 추적됩니다 . 유체 변수 F의 비율은 유체가 차지하는 영역을 찾습니다. 유체에 고정 된 자유 표면이있는 경우 유체를 정의하는 F 값도 안정된 값을 유지해야합니다. F가 일정하려면 표면에 수직 인 유체 속도가 0이어야합니다. 물론 표면에서의 접선 유체 속도는 0 일 필요는 없습니다. 예를 들어, 위어 위의 흐름에는 일정한 흐름이 있지만 계단에서 나오는 흐름의 위치와 모양은 변하지 않습니다.

자유 표면 흐름에 대한 정상 상태 솔버를 사용하려면 흐름의 비압축성을 유지하면서 정상 표면 속도를 0으로 유도하는 방법을 찾아야합니다.

이를 수행하는 한 가지 방법은 정상 속도를 0으로 유도하는 방식으로 표면 압력을 조정하는 것입니다. 특히 정상 속도에 비례하는 총 표면 압력에 “댐핑”압력 기여를 추가하는 것입니다. 속도는 표면 밖으로 향하고 그렇지 않으면 음수입니다.

정상 속도가 0에 가까워지면 수정 압력도 0이되어야 표면이 고정 위치를 초과하지 않게됩니다. 물론 보정이 너무 크면 오버 슈트가 발생할 수 있습니다. 따라서 안정적인 보정 적용을 위해서는 몇 가지 제한 요소가 있어야합니다.

계수 약어 ssacc 을 나타내며, S는 teady- S 테이트 액세서리 elerator이 새로운 옵션을 활성화하는 프로그램 입력에 추가되었다. ssacc 의 값 은 편리한 상한 인 1.0보다 작거나 같아야합니다. 프로그램 내에서 댐핑 압력에 자동으로 적용되는 여러 제한 기가 불안정 해 지거나 일시적인 현상에 악영향을 미치는 것을 방지합니다.

안정성 및 댐핑 리미터에 대한 이전 문제는 강조되어야합니다. 정상 상태 가속기를 사용하면 자유 표면 흐름의 모든 과도 현상이 더 이상 완전히 사실적인 것으로 볼 수 없습니다. 댐핑 압력은 물리적 인 힘이 아니라 파동 전파와 반사를 줄이는 메커니즘입니다. 댐퍼는 큰 과도 현상의 발생을 방해하지 않도록 고안되었으며 흐름이 안정됨에 따라 안정된 결과를보다 빠르게 얻는 데에만 기여해야합니다. 그러나 사용자는 리미터가 예상하지 못한 초과 댐핑에 대해 주의를 기울여야 합니다. 이는 댐핑 계수 ssacc 의 입력 값을 줄임으로써 제거 할 수 있습니다 .

두 가지 예는 정상 상태 가속기의 댐핑 메커니즘이 어떻게 작동하는지 보여줍니다.

Steady-State Accelerator Examples

Collapse of Raised Fluid Column

첫 번째 예는 길이 100cm, 깊이 5cm의 2 차원 물 웅덩이로 구성됩니다. 물을 담은 탱크의 모든 경계는 대칭 경계입니다. 수영장 중앙에는 폭 10cm, 높이 3cm의 수영장 위에 물 블록이 있습니다. 이 블록은 중력으로 인해 물에 떨어지고 충돌 지점에서 멀리 이동 한 다음 탱크 끝에서 반사되는 파도를 생성합니다. 100 초 후에도 반복되는 반사 때문에 여전히 상당한 파동 작용이 있습니다 (그림 1).

새로운 정상 상태 가속기를 계수 ssacc = 1.0 과 함께 사용하면 모든 파동이 빠르게 감쇠되어 거의 평평한 표면이됩니다. 일부 잔류 흐름은 표면 아래에 남아 있지만 점도의 작용으로 서서히 감쇠됩니다 (그림 2). 이 예에서 추가 된 댐핑은 특히 인상적입니다.

Figure 1. Column collapse without damping. Times of flow plots are 0.0, 10.0, and 100.0s. Bottom figure is the mean kinetic energy vs. time.
Figure 2. Column collapse with damping coefficient ssacc=1.0 at times of 0.0, 10.0 and 100.0s. Bottom figure is the mean kinetic energy vs. time.

 

사각형 격자에서 45 °의 정사각형 채널에서 모세관 상승

수직 채널에서 유체의 모세관 상승은 간단한 분석할 수 있으며 솔루션이 있는 양호한 정상 상태 문제입니다. 중력에 대해 상승 된 유체의 양은 벽의 접착력, 즉 접촉각의 코사인에 표면 장력 곱하기 접촉 선 길이에 의해 결정됩니다. 이 예에서 유체는 물이며 표면 장력은 70 dynes / cm이고 접촉각은 30 °입니다. 채널은 단면이 정사각형이며 가장자리 길이가 0.707cm이고 직사각형 격자에서 45 ° 회전합니다. 문제가 x 및 y 방향으로 대칭을 이루기 때문에 그리드의 사분면 만 모델링됩니다. 그리드의 바닥에는 제로 게이지 압력의 물이 있으며 그리드의 가장자리 길이는 0.0125cm (41x41x80 셀)입니다. 상승시켜야하는 이론적 유체 량은 0.04373cc입니다. 그림 3a는 정상 상태 결과를 보여줍니다. 이는 감쇠 사용 여부와 비슷합니다. 댐핑없이 계산된 유체의 양은 이론 값보다 1.74 % 높습니다. 그림 3b와 같이 댐핑이 있는 경우에는 2.24 %가 너무 높습니다. 가속기를 사용하면 정상 상태는 약 0.15 초에 도달하는 반면 표준 솔버는 0.8 초 후에 만 ​​정상 상태 솔루션을 생성하므로 5 배 이상 더 오래 걸립니다.

Figure 3a. Capillary rise in square channel without damping pressures.
Figure 3b. Histories of fluid volume in the two simulations (blue is with damping).

ssacc가 1.0보다 작으면 댐핑이 적어 수렴에 더 빨리 도달합니다. 1.0을 포함한 모든 ssacc 값은 댐핑되지 않은 ssacc = 0.0 경우와 비교하여 이론과 밀접하게 일치하고 후면 벽에 적은 양의 유체를 나타내는 수렴된 솔루션을 만듭니다.

뒤쪽 벽에있는 작은 유체 조각은 평형 위치를 초과하는 유체의 오버 슈트에서 발생하며, 이는 점성력으로 인해 정착하는 데 오랜 시간이 필요한 소량의 유체를 벽에 남기고 뒤로 떨어집니다. 이 오버 슈트는 ssacc 가 0이 아닐 때 제거됩니다 .

Rivulet Formation in Slide Coating

Simulation of Transient and Three-Dimensional Coating Flows Using a Volume-of-Fluid Technique

Volume-of-Fluid 기법을 사용한 과도 및 3 차원 코팅 흐름 시뮬레이션

슬라이드 코팅 흐름은 정밀 필름 코팅 제품의 제조에 널리 사용됩니다. 코팅 속도를 높이고 코팅 필름의 성능을 향상시키기 위해 슬라이드 코팅 공정을 더 잘 이해하기 위해 상당한 노력을 기울이고 있습니다. 예를 들어 Chen1과 같이 잘 정의 된 한계 이상으로 코팅 속도를 높이면 코팅 비드가 완전히 파손될 수 있음이 입증되었습니다.

이 논문에서는 유체 표면의 임의, 3 차원 및 시간에 따른 변형을 설명 할 수있는 계산 방법에서 얻은 슬라이드 코팅 흐름의 시뮬레이션 결과를 제시합니다. 상용 프로그램에서 사용할 수있는이 방법은 VOF (Volume-of-Fluid) 기술 3,4로 유체를 추적하는 고정 그리드를 사용합니다. 표면 장력, 벽 접착력, 유체 운동량 및 점성 응력은 분석에서 완전히 설명됩니다.

기본 방법은 딥 코팅 데이터와의 비교를 통해 설명됩니다 5. 그런 다음 접촉 선과 동적 접촉각이 우리의 방법에서 암시 적으로 처리되는 방법에 대한 논의를 제시합니다. VOF 기술을 사용하기 때문에 유체를 포함하는 각 제어 볼륨에 작용하는 힘의 합계 만 필요합니다. 그러면 접촉 선의 위치와 동적 접촉각이 계산 된 힘 균형에서 자동으로 발생합니다. 우리의 기술은 코팅 흐름에서 시작 및 비드 분해 현상의 예와 함께 설명됩니다.

그림에서 볼 수 있듯이 신속한 공정의 경우 당사의 접근 방식은 기존 분석 방법으로는 달성하기 어려운 코팅 공정 설계 및 최적화 시뮬레이션을위한 효율성과 견고성을 제공합니다.

Introduction

모든 코팅 공정에는 일정한 조건을 달성하기 전에 코팅 재료가 큰 변형을 겪는 일종의 시작 기간이 포함됩니다. 시작 프로세스의 우수한 특성화는 낭비를 줄이고 프로세스가 원하는 한계 내에서 작동하는지 확인하는 데 종종 중요합니다.

다양한 섭동에 대한 코팅 흐름의 과도 ​​응답에 대한 유사한 이해가 또한 바람직하여 코팅 비드의 파손 및 코팅의 불균일성을 피할 수 있습니다. 코팅 흐름의 역학은 일반적으로 비선형이고 다양한 경쟁 물리적 프로세스의 결합 된 상호 작용을 포함하기 때문에 이론적 조사를 수행하기 위해 특수한 계산 도구에 의존해야합니다.

이 작업을 위해 선택한 모델링 도구의 장점은 고정 그리드를 통해 임의의 유체 변형을 추적 할 수있는 강력한 수치 기법 인 VOF (Volume-of-Fluid) 방법을 사용한다는 것입니다. 코팅 흐름 분석에 중요한 프로그램의 다른 기능과 함께 이것이 수행되는 방식은 다음 섹션에서 설명합니다.

Overview of Numerical Method

여기에 사용 된 수치 프로그램 FLOW-3D®는 1960 년대 중반 Los Alamos National Laboratory에서 개발 된 Marker-and-Cell (MAC) 방법 6에서 유래되었습니다. 원래 MAC 방법에 대한 많은 개선이 수년에 걸쳐 이루어졌습니다.

본 출원에서 가장 흥미로운 것은 유체 영역을 찾기 위해 연속적인 유체 부피 함수에 의해 개별 마커 입자를 대체하는 것입니다. VOF 방법에서는 관심있는 계산 영역을 포함하는 사각형 제어 볼륨의 고정 그리드가 구성됩니다. 각 제어 볼륨에 대해 숫자 F는 액체가 차지하는 볼륨의 비율을 표시하기 위해 유지됩니다.

F 함수를 사용하는 것 외에도 VOF 방법은 날카로운 액체-가스 인터페이스를 유지하는 방식으로 직사각형 셀의 고정 그리드를 통해 F 함수를 전진시키기 위해 특수 수치 기법을 사용합니다. 마지막으로 VOF 방법은 경계면에서 적절한 법선 및 접선 응력 조건을 충족하기 위해 신중하게 구현 된 자유 표면 경계 조건 세트를 사용합니다. 접근 방식의 또 다른 특징은 복잡한 기하학적 영역을 정의하는 방식입니다.

장애물은 제어 볼륨의 일부를 차단할 수 있도록하여 고정 그리드에 포함됩니다. 각 제어 볼륨에서 흐름을 위해 열린 분수 영역 및 볼륨은 지오메트리 표현으로 저장됩니다. FAVOR 방법 7이라고하는이 방법은 형상을 질량, 운동량 및 에너지에 대한 이산화 된 방정식에 자동으로 통합합니다. VOF 및 FAVOR 방법을 사용하면 코팅 문제에 대한 지오메트리 및 초기 유체 구성을 정의하는 데 필요한 복잡한 그리드 생성 프로세스가 없기 때문에 시간과 노력이 절약됩니다.

다음 섹션에서는 플랫 시트에 코팅을 담그는 응용 프로그램과 함께 기본적인 수치 방법의 유용성을 설명합니다.

Dip Coating – A Validation Test

Lee와 Tallmadge는 액체 수조에서 수직으로 인출 된 평판에 딥 코팅하는 과정에 대해 광범위한 조사를 수행했습니다.

이 프로세스는 다양한 상업용 응용 프로그램에서 널리 사용됩니다. 그들의 연구는 2 차원 흐름 (즉, 가장자리 효과 없음)에 초점을 맞추고 실험 데이터에 맞는 경험적 매개 변수를 포함하는 분석 표면 프로파일로 구성되었습니다. 0.085에서 23.9 사이의 모세관 수에 대한 실험 데이터가 수집되었으며, 레이놀즈 수는 0.044에서 12.7 사이입니다. 필름 두께에 대한 실험 데이터는 약 10 % 이하로 추정되는 오류를 가졌습니다.

이 실험에 대한 계산 모델은 코팅 할 시트의 수직 (접선) 속도와 동일한 수직 (접선) 속도가 주어진 직사각형 욕조로 구성되어 매우 간단합니다. 처음에 코팅액은 수평면을 가지며 시트는 충동 적으로 시작됩니다 (그림 1c 참조). 다양한 모세관 수 사례가 시뮬레이션되었으며 모든 경우에 예측 된 필름 두께는 실험 오차 범위 내에있었습니다. 예를 들어 모세관 번호 1.17에 해당하는 경우를 고려하십시오. 시트를 3.31cm / s에서 수조 (밀도 0.885gm / cc, 표면 장력 32.7dynes / cm 및 점도 1159.4cp를 갖는 점성 윤활유)에서 꺼냈다. 우리는 2.5cm의 욕조 너비와 2.0cm의 깊이 (35 x 25 그리드 셀)를 사용했습니다.

필름 흐름을 캡처하기 위해 욕조 위의 2.0cm 영역이 모델에 포함되었습니다 (수직으로 추가 25 개 셀 필요). 수조의 오른쪽은 유체 높이가 일정하게 유지되고 압력이 수압이고 흐름이 계산 영역으로 들어갈 수있는 열린 경계 였지만 휴식에서 시작해야했습니다. 이른바 “정체”경계 조건은 움직이는 시트의 오른쪽으로 충분히 멀리 떨어져있는 경우 수평 무한 욕조에 대한 좋은 근사치입니다. 모델링이 필요한 수조의 폭을 설정하기 위해 여러 가지 계산이 수행되었으며, 필름 두께가이 폭에 크게 민감하지 않다는 것이 밝혀졌으며 그 결과는 실험에서도 발견되었습니다.

그림 1a는 초기 조건, 그림 1b는 계산 된 과도 상태의 스냅 샷, 그림 1c는 최종 정상 상태 결과를 보여줍니다. 처음에 시트에 의해 그려지는 액체 팁의 모양은 정적 접촉각 (즉, 시트와 액체 사이의 접착력)에 따라 달라지며 임의로 10 도로 취해졌습니다. 액체가 끌어 올려짐에 따라, 배출되는 액체 필름을 대체하기 위해 시트쪽으로 흐름이 시작되어야한다는 신호로서 함몰 파가 나머지 수조에 대한 신호로 오른쪽으로 이동합니다. 약 5.0 초만에 정상 상태에 도달합니다. 필름 두께는 0.145cm로 계산되었으며, 이는 0.142cm의 측정 값과 매우 일치합니다.

Rivulet Formation in Slide Coating
Rivulet Formation in Slide Coating

자세한 내용은 본문을 참고하시기 바랍니다.

Capillary Flows/Capillary Filling/Thermocapillary Switch/Capillary Absorption/Marangoni flow

Capillary Flows

모세관 흐름은 일반적으로 미세 유체 장치에서 발생합니다. 예를 들어, 바이오 칩 설계에서는 한 곳에서 다른 곳으로 액체 용액을 전달하기 위해 긴 마이크로 채널이 자주 사용됩니다. 입구 채널은 액체 저장소에 연결되고 표면 장력은 액체를 마이크로 채널로 끌어 당깁니다 (액체가 칩 표면에 “젖은”경우). 이 페이지에서는 충진, 흡수 및 전환과 같은 모세관 흐름 분석에서 FLOW-3D의 특정 응용 분야를 다룹니다.

Marangoni flow in a dish of water that is heated at its center.

Marangoni flow는 중앙에서 데워진 물이 담긴 접시에 흐릅니다. 불균일한 표면 장력에 의해 생성 된 흐름은 20ºC의 초기 온도에서 0.75cm 깊이의 얕은 8.0cm 직경의 물 접시에 의해 입증됩니다. 원형 접시의 중앙에는 직경 0.5cm의 원통형 막대가 있습니다. 80 Cº의 온도로 가열하고 0.05 cm 깊이까지 수면에 담근다. 핫로드 근처의 물이 가열됨에 따라 표면 장력이 0.1678 dyne / cm / ºC만큼 감소하여 표면이 접시의 바깥 쪽 테두리쪽으로 후퇴합니다. Retraction는 처음에 표면에 뿌려진 질량없는 마커 입자로 표시됩니다.

Capillary Absorption

고체 물질의 기공에 모세관 흡수 때문에 액체와 고체 사이의 접착 발생합니다. 이 같은 흡수의 간단하면서도 유용한 시험은 핀란드 ABO Akademi 대학의 마르티 Toivakka에 의해 제안되었습니다. 테스트 기공은 ± 1.0 μm의 측면 벽 1.0 μm의 반경 원호입니다. 팽창 목에 연결된 넓은 2차원 채널로 구성되어 있습니다. 체적력의 부재 하에서, 표면장력 과 wall adhesion pull liquid 는 액체와 고체 사이의 static contact angle에 의해 결정됩니다. 첨부된 그림은 FLOW-3D가 올바르게 특정 접촉 각도 (유체는 적색표현) 충전 레벨을 계산하는 것을 나타냅니다.

Thermocapillary Switch

액체의 작은 덩어리나 가벼운 빔의 경로에서 움직이는 굴절, 혹은 반사로 다른 길로 리디렉션 할 수 있습니다. 이 개념은 특히나 한번 빔 내부 반사로 인해 갇혀 있는 섬유에 들어가 광학 섬유로 연결에서 매력적입니다. 어떠한 복잡성의 광 회로를 만들려면, 하나의 광섬유에서 다른 가벼운 방향을 바꿀 수 있는“스위치”를 둘 필요가 있습니다.

The animation above shows a FLOW-3D simulation of a drop of water in a 14mm-wide channel that is being heated at the bottom.

Capillary Filling

모세관 충전 과정을 이해하는 것은 칩 설계에 중요합니다.. 액체 흐름 통로의 다른 형상 포획 기포의 가능성 등의 충전 공정의 기술은, 같은 챔버와 칩의 내부 구조를 배치 기둥 분할하고, 밸브 결합에 설계자 안내 등 다양한 모세관 충전 동작이 발생할 수 있습니다.

시뮬레이션은 아래의 모세관 작용의 분석 예측의 유효성을 검사합니다. 모세관 채우기는 정확하게 표면 장력과 중력에 의해 균형을 잡습니다.이것은 FLOW-3D에 의해서 정확하게 예측되는 기본적인 과정입니다.

Porous Materials

Porous Materials

다공성 매체에서의 유체 수송을 수치적으로 모델링하는 것은 어려운 일이 될 수 있지만, FLOW-3D에는 다공성 물질과 관련된 문제를 해결하기 위한 많은 중요한 기능이 포함되어 있습니다. FAVOR™  기술에는 사용자가 연속적인 다공성 미디어를 나타낼 수 있도록 필요한 다공성 변수가 포함되어 있습니다. FLOW-3D를 사용하면 포화 및 불포화 유량 조건을 모두 시뮬레이션할 수 있습니다. 전력 법칙 관계를 통해 사용자는 불포화 유량 조건에서 모세관 압력(capillary pressure)과 포화도 사이의 비선형 관계를 모형화할 수 있습니다.

별도의 주입 및 배출 곡선을 사용하여 이력 현상을 모델링할 수 있습니다. 서로 직접 접촉하는 경우에도 서로 다른 다공성, 투과성 및 습윤성 특성을 서로 다른 장애물에 할당할 수 있습니다. 투과성은 흐름 방향에 따라 지정할 수 있으므로 다공성 매체의 비등방성 동작을 모델링할 수 있습니다. 유체와 다공성 매체 사이의 열 전달이 고려될 수 있습니다.

Precision Droplet Creation

Precision Droplet Creation

신속하고 반복적이며 정밀한 액적 생성 및 증착, 작은 피처의 인쇄 또는 패턴 화 (예 : l = 10-3-1 mm), 분무를 통해 제어되고 균일 한 두께로 박막을 형성하는 것이 매우 중요합니다.

다양한 이전 및 새로운 산업 응용 프로그램 (1-5). 액체 이동 및 액적 형성 / 증착 공정에는 복잡한 자유 표면 흐름과 자발적인 모세관 구동 불안정, 얇아 짐 및 핀치 오프 (1-5)를 겪는 원주 형 목의 형성이 포함됩니다. 단순한 뉴턴 유체 및 비탄성 유체에 대한 액적 형성 및 액체 전달을 분석하기위한 실험적, 이론적 및 1 차원 시뮬레이션 연구를 사용하여 이루어진 진전에도 불구하고 인쇄 및 스프레이에 대한 기계적 이해는 여전히 어려운 과제입니다.

현재 계산 노력의 주요 동기는 FLOW-3D에 포함 된 VOF (체적-유체) 접근 방식을 사용하여 뉴턴 유체의 핀치 오프 역학에 대한 기계 론적 이해를 얻을 수있는 가능성을 조사하는 것입니다. 우리는 계산 분석이 자기 유사 모세관 얇아 짐 및 핀치 오프 역학을 결정하는 모세관, 관성 및 점성 응력의 복잡한 상호 작용을 포착 함을 보여줍니다. 뉴턴 유체의 낙하 형성 및 분리를 위해, 우리는 이론 및 1D 시뮬레이션 (1-7)뿐만 아니라 실험 (1, 1)에서 예상되는 범용 스케일링 법칙을 사용하여 계산 분석에서 얻은 자기 유사 목 진화를 설명 할 수 있음을 보여줍니다.

2, 8-12). 이러한 프로토 타입 흐름을 시뮬레이션하는 데있어 우리의 성공은 FLOW-3D를 사용하여 유한 시간 특이점, 위성 낙하 형성 및 더 복잡한 형상의 인쇄 가능성에 대한 신중한 계산 분석을 위해 필요한 단계이며 설명이나 연구가 훨씬 더 어렵습니다. 1D 모델과 실험을 사용합니다.

Courtesy of University of Illinois at Chicago

Read the Computational Analysis of Drop Formation and Detachment case study.

Digital Microfluidics

Electrowetting은 전기장을 사용하여 표면 습윤 특성을 변경하는 과정입니다. Digital microfluidics는 전기 습식이 개별 유체 방울을 제어하고 조작하는데 사용되는 미세 유체 분야입니다. 이 아이디어는 디지털 마이크로 일렉트로닉스에서 영감을 얻었지만 전류 대신 이산 (또는 디지털화 된)액적을 사용하여 특정 시간 내에 특정 거리에 포함된 특정 양의 유체 또는 반응물을 이동합니다. 디지털 마이크로 플루이딕스는 높은 재구성 가능성과 대규모 병렬화를 통해 프로세스 속도를 높일 수있는 능력 때문에 다양한 바이오칩 설계에서 응용 분야를 찾습니다.

가장 중요한 표면 습윤 특성은 유체와 표면 사이의 접촉각입니다. FLOW-3D의 강력한 표면장력 모델은 전기 운동 모델과 함께 유전 영동, 열 모세관 작동 (온도에 따른 표면 장력을 통한 작동) 및 전기 습윤 자체와 같은 디지털 미세 유체 공정에서 습윤 역학을 포착하는 데 사용됩니다.

Continuous Flow Microfluidics

Continuous Flow Microfluidics

연속 흐름 미세 유체는 연속성을 깨지 않고 제작 된 마이크로 채널을 통해 액체 흐름을 조작하는 것입니다. 유체 흐름은 마이크로 펌프 (예 : 연동 펌프 또는 주사기 펌프)와 같은 외부 소스 또는 전기, 자기 또는 모세관 힘과 같은 내부 메커니즘에 의해 설정됩니다. 연속 유동 미세 유체 학은 미세 및 나노 입자 분리기, 입자 집속, 화학적 분리는 물론 단순한 생화학 적 응용을 포함한 다양한 응용 분야에서 응용 분야를 찾아 내지 만 높은 수준의 제어가 필요한 경우에는 선택 방법이 아닐 수 있습니다.

이 범주에 속하며 FLOW-3D를 사용하여 성공적으로 시뮬레이션한 프로세스 또는 장치로는 Joule 가열, 액체 게이트, 마이크로 유체 회로, 전기-오토믹 밸브, 입자 집중, 분류 및 분리, POC(Point-of-Care) 모세관 유량 장치 및 패턴 있는 표면 장치가 있습니다.

Sketch of cross section of the device
Capillary Flows
Electro osmosis
Electro-osmosis
Simulating joule heating
Joule Heating
Patterned surfaces in micro channels
Lab-on-a-chip
Magnetic fields
Magnetic Fields
Pneumatic valve
Microfluidic Circuits
Hong chamber simulations
Mixing Dynamics
Buoyancy dominant sorting
Particle Sorting

Capillary Flows

Capillary Flows

모세관 흐름은 일반적으로 미세 유체 장치에서 발생합니다. 예를 들어, 바이오 칩 설계에서 긴 마이크로 채널은 종종 액체 용액을 한 장소에서 다른 장소로 전달하는 데 사용됩니다. 입구 채널은 액체 저장소에 연결되고 표면 장력이 액체를 마이크로 채널로 당깁니다(액체가 칩 표면에 “습기”되는 경우). 이 페이지에서는 충전, 흡수 및 전환과 같은 모세관 흐름 분석에서 FLOW-3D에 대한 몇 가지 특정 용도에 대해 다룹니다.

Marangoni Flows

마랑고니는 그 중심에 가열된 물 접시에 흐릅니다. 균일하지 않은 표면 장력에 의해 발생하는 흐름은 20ºC의 초기 온도에서 깊이 0.75cm의 얕은 8.0cm의 물 접시에 의해 입증됩니다. 원형 접시 중앙에 놓인 원통형 막대는 직경 0.5cm로 80Cº의 온도로 가열되고 0.05cm의 깊이까지 수면에 잠깁니다. 핫 로드 주변의 물이 가열되면 표면 장력이 0.1678dyne/cm/ºC만큼 감소하여 표면이 접시의 바깥쪽 림 쪽으로 수축됩니다. 수축은 처음에 표면에 뿌려진 질량이 없는 마커 입자에 의해 나타납니다.

Capillary Filling

모세관 충전 공정을 이해하는 것은 칩 설계에 중요합니다. 액체 흐름 경로의 기하학적 구조가 다르면 기포를 고정할 수 있는 등의 모세관 충진 동작이 달라질 수 있습니다. 충전 프로세스에 대한 지식은 설계자가 챔버, 결합 기둥, 분할 및 밸브와 같은 칩의 내부 구조를 정렬하는 데 도움이 됩니다. 오른쪽의 시뮬레이션은 모세관 작용의 분석적 예측을 검증합니다. 모세관 충전은 표면 장력과 중력에 의해 균형을 이루며, 이는 FLOW-3D로 정확하게 예측되는 기본 공정입니다.

Thermocapillary Switch

910/5000광선의 경로 안팎으로 이동하는 소량의 액체는 굴절이나 반사를 통해 다른 경로로 방향을 바꿀 수 있습니다. 이 개념은 광선이 광섬유에 들어가면 내부 반사에 의해 포착되는 광섬유와 관련하여 특히 매력적입니다. 복잡한 광학 회로를 만들려면 한 광섬유에서 다른 광섬유로 빛을 리디렉션 할 수있는 “스위치”가 필요합니다.

제안 된 한 가지 개념은 열 모세관을 기반으로합니다. 광섬유 광선을 교차하는 마이크로 채널에 액체의 작은 방울을 놓습니다. 방울이 채널을 따라 빔이 통과해야하는 곳으로 이동하면 빔이 다른 섬유로 반사됩니다. 방울은 양면을 다르게 가열하여 이동합니다. 이것은 방울이 채널의 더 차가운 끝쪽으로 당겨 지도록 방울의 양쪽에있는 반월판의 표면 장력의 변화를 일으 킵니다.

Whole Blood Spontaneous Capillary Flow

Sketch of the cross section of the device (w=150 µm, h1=300 µm, h2=1200 µm, α=14.5o)

모세관 기반 마이크로 시스템은 추가 작동 메커니즘이 필요하지 않기 때문에 저렴하고 제작하기 쉽습니다. 마이크로펌프나 주사기와 같은 일반적인 마이크로 시스템은 부피가 크고 휴대할 수 없는 흐름 작동을 필요로 합니다.

버팔로 대학의 최근 연구는 모세관 유동 작용을 사용하여 미세 기기에서 액체를 이동시키는 간단한 해결책을 연구했습니다. 이 작업은 FLOW-3D를 사용하여 수정된 V-그루브 채널에서 자발적 모세관 흐름을 시뮬레이션합니다. 좁은 V-그루브 기하학(왼쪽)은 전혈과 같은 높은 점도의 유체도 이 유체를 통해 이동할 수 있기 때문에 좋은 솔루션을 제공합니다. 홈의 끝부분은 자발적인 모세관 흐름을 촉진하고 평행판은 충분한 혈액수송을 보장합니다.

본 연구에서는 FLOW-3D를 사용하여 채널 내 유체 헤드의 유속과 액체 전방의 진행을 추정합니다.

결과는 실험 및 분석(간단한) 결과와 비교됩니다. 아래 그림은 수치, 실험 및 분석 결과의 비교를 보여줍니다. FLOW-3D 결과는 실험 결과와 매우 일치합니다.

FLOW-3D Results

Analysis A: FLOW-3D results in red circles at the mid flow height, experimental results in green dots recorded at the medium fluid height, analytical results in green dashes
Analysis B: FLOW-3D results in red circles at the mid flow height, experimental results in green dots recorded at the medium fluid height, analytical results in green dashes

Animation of the results post-processed in FlowSight.

References

J. Berthiera, K.A. Brakke, E.P. Furlani, I.H. Karampelas, V. Pohera, D. Gosselin, M. Cubizolles, P. Pouteau, Whole blood spontaneous capillary flow in narrow V-groove microchannels, Sensors and Actuators B: Chemical, 2014

다공성 매체 / Porous Media

다공성 매체 / Porous Media

FLOW-3D는 다공성 매체 내의 포화 및 불포화 흐름을 모두 시뮬레이션할 수 있습니다. 포화된 다공성 미디어 흐름은 포화 구역과 불포화 구역 사이에 예리한(또는 거의 날카로운) 계면이 있고 계면에 특정 모세관 압력이 존재하는 상황에 적용됩니다. 이러한 상황은 지하수 흐름에서 발생합니다. 불포화 다공성 미디어 흐름은 포화 구역에서 불포화 구역으로 점진적으로 전환되는 상황에 적용됩니다. 이러한 상황에서는 설정된 모세관 압력이 없습니다. 모세관 압력은 현재 포화 수준과 다공성 물질 내 포화 이력의 함수입니다.

두 경우 모두 각 성분에 대해 서로 다른 다공성, 투과성 및 습윤성(모세관 압력 또는 모세관 압력 대 포화도)을 독립적으로 지정할 수 있으며 투과성은 등방성(모든 방향에서 동일) 또는 비등방성(흐름 방향에 따라 달라짐)일 수 있습니다.

아래 동영상은 종이와 같은 다공성 물질로 스며드는 물방울의 경우를 보여줍니다. 이 경우 다공성 물질은 불포화 상태로 모델링되므로 습윤성은 국부 포화에 따라 달라집니다. 이미 젖은 영역은 모세관 압력이 더 강한 반면, 낙하 가장자리에 있는 영역은 모세관 압력이 더 낮습니다. 이 작업은 별도의 주입 및 배출 곡선을 사용하여 수행됩니다. 따라서 방울이 재료에 균일하게 퍼지지 않습니다. 이러한 행동은 젖은 종이 타월을 짜는 것으로 볼 수 있습니다; 모든 물을 짜내는 것보다 종이를 적시는 것이 훨씬 쉽습니다.


다공성 매질에 흡수 된 물방울 시뮬레이션

다공성 매체에서의 불포화 흐름은 포화 흐름 조건에서는 존재하지 않는 많은 복잡한 현상을 수반합니다. 예를 들어 구성을 알 수 없는 자유 경계와 모세관 힘이 존재하여 액체를 포화도가 낮은 영역으로 끌어들이는 큰 음압을 발생시킬 수 있습니다. 또한 모세관 압력은 실험적인 판단과 모델링을 더욱 어렵게 만드는 이력(hysteresis) 동작을 보일 수 있습니다. 불포화 흐름과 관련된 합병증은 가장 간단한 경우를 제외한 모든 상황에서 수치적 해결 절차의 필요성을 나타냅니다. 이러한 유형의 흐름의 자유 표면적인 측면 때문에, FLOW-3D® 프로그램을 불포화 흐름의 일반적인 사례로 확장할 것을 생각하는 것은 당연합니다. 이 확장 작업을 수행하는 데 필요한 수정 사항은 아래 보고서에 설명되어 있습니다. 이후의 섹션에서 더 자세히 설명했듯이, 물질을 통과하는 흐름을 정확하게 모형화하기 위해서는 다공성 물질과 이를 관통하는 액체에 대한 충분한 경험적 데이터가 필요합니다. 모세관 압력과 투과성을 위해 여기에 보고된 모델에는 일부 재료에 대한 수정이 필요할 수 있는 일반적인 특성이 있습니다.

보고서 원문  : UNSATURATED FLOW IN POROUS MEDIA

FLOW-3D What’s New Ver.12.0

FLOW-3D v12.0은 그래픽 사용자 인터페이스 (GUI)의 설계 및 기능에서 매우 큰 변화를 이룬 제품으로 모델 설정을 단순화하고 사용자 워크 플로를 향상시킵니다. 최첨단 Immersed Boundary Method(침수경계 방법)은 FLOW-3D v12.0 솔루션의 정확성을 높여줍니다. 다른 주요 기능으로는 슬러지 침강 모델, 2-Fluid 2-Temperature 모델 및 Steady State Accelerator가 있으며,이를 통해 사용자는 자유 표면 흐름을 더욱 빠르게 모델링 할 수 있습니다.

Physical and Numerical Model

Immersed boundary method

힘과 에너지 손실에 대한 정확한 예측은 고체 주위의 흐름과 관련된 많은 엔지니어링 문제를 모델링하는 데 중요합니다. 새 릴리스 FLOW-3 Dv1.2.0에는 이러한 문제점 해결을 위해 설계된 새로운 고스트 셀 기반 Immersed Boundary Method (IBM)가 있습니다. IBM은 내 외부 흐름 해석을 위해, 벽 근처에서 보다 정확한 해를 제공하여 드래그 앤 리프트 힘의 계산을 향상시킵니다.힘과 에너지 손실의 정확한 예측은 고체 주위의 흐름을 포함하는 많은 공학적 문제를 모델링 하는데 중요합니다.

Two-field temperature for the two-fluid model

2 유체 열전달 모델은 각 유체에 대한 에너지 전달 방정식을 분리하기 위해 확장되었습니다. 각 유체는 이제 자체 온도 변수를 가지므로 인터페이스 근처의 열 및 물질 전달 솔루션의 정확도가 향상됩니다. 인터페이스에서의 열전달은 이제 시간의 표 함수가 될 수 있는 사용자 정의 열전달 계수에 의해 제어됩니다.

블로그 보기

Sludge settling model

새로운 슬러지 정착 모델은 수처리 애플리케이션에 부가되어 사용자들이 수 처리 탱크와 클래리퍼의 고형 폐기물 역학을 모델링 할 수 있게 해 줍니다. 침전 속도가 분산상의 액적 크기의 함수 인 드리프트-플럭스 모델과 달리, 침전 속도는 슬러지 농도의 함수이며 기능 및 표 형식으로 입력 할 수 있습니다.

개발노트 읽기

Steady-state accelerator for free surface flows

이름에서 알 수 있듯이 정상 상태 가속기는 정상 상태 솔루션에 대한 접근을 빠르게합니다.
이것은 작은 진폭 중력과 모세관 표면파를 감쇠시킴으로써 달성되며 자유 표면 흐름에만 적용 할 수 있습니다.

개발노트 읽기

Void particles

Void particles 가 기포 및 상 변화 모델에 추가되었습니다. Void particles는 붕괴 된 Void 영역을 나타내며, 항력 및 압력을 통해 유체와 상호 작용하는 작은 기포로 작용합니다. 주변 유체 압력에 따라 크기가 변하고 시뮬레이션이 끝날 때의 최종 위치는 공기 유입 가능성을 나타냅니다.

Sediment scour model

퇴적물 수송 및 침식 모델은 정확성과 안정성을 향상시키기 위해 정비되었습니다. 특히 퇴적물 종의 질량 보존이 크게 개선되었습니다.

개발 노트 읽기>

Outflow pressure boundary condition

고정 압력 경계 조건에는 압력 및 유체 분율을 제외한 모든 유량이 해당 경계의 상류의 유량 조건을 반영하는 ‘유출’옵션이 포함됩니다. 유출 압력 경계 조건은 고정 압력 및 연속 경계 조건의 하이브리드입니다.

Moving particle sources

시뮬레이션 중에 입자 소스를 이동할 수 있습니다. 시간에 따른 병진 및 회전 속도는 표 형식으로 정의됩니다. 입자 소스의 운동은 소스에서 방출 된 입자의 초기 속도에 추가됩니다.

Variable center of gravity

기변 무게중심은 중력 및 비관 성 기준 프레임 모델에서, 시간의 함수로서 무게 중심의 위치는 외부 파일에서 테이블로서 정의 될 수있다. 이 기능은 연료를 소비하고 분리 단계를 수행하는 로켓과 같은 모형을 모델링 할 때 유용합니다.

공기 유입 모델

가장 간단한 부피 기반 공기 유입 모델 옵션이 기존 질량 기반 모델로 대체되었습니다. 질량 기반 모델은 부피와 달리 주변 유체 압력에 따라 부피가 변화하는 동안 흡입된 공기량이 보존되기 때문에 물리학적 모델입니다.

Tracer diffusion

유동 표면에서 생성된 추적 물질은 분자 및 난류 확산 과정에 의해 확산될 수 있으며, 예를 들어 실제 오염 물질의 동작을 모방한다.

Model Setup

Simulation units

온도를 포함하여 단위 시스템은 완전히 정의해야하는데 표준 단위 시스템이 제공됩니다. 또한 사용자는 다양한 옵션 중에서 질량, 시간 및 길이 단위를 정의 할 수 있으므로 사용자 정의가 가능한 편리한 단위를 사용할 수 있습니다. 사용자는 압력이 게이지 또는 절대 단위로 정의되는지 여부도 지정해야합니다. 기본 시뮬레이션 단위는 기본 설정에서 설정할 수 있습니다. 단위를 완전히 정의하면 FLOW-3D 가 물리량의 기본값을 정의하고 범용 상수를 설정하여 사용자가 요구하는 작업량을 최소화 할 수 있습니다.

Shallow water model

Manning’s roughness in shallow water model

Manning의 거칠기 계수는 지형 표면의 전단 응력 평가를 위해 얕은 물 모델에서 구현되었습니다. 표면 결함의 크기를 기반으로 기존 거칠기 모델을 보완하며 이 모델과 함께 사용할 수 있습니다. 표준 거칠기와 마찬가지로 매닝 계수는 구성 요소 또는 하위 구성 요소의 속성이거나 지형 래스터 데이터 세트에서 가져올 수 있습니다.

Mesh generation

하단 및 상단 경계 좌표의 정의만으로 수직 방향의 메시 설정이 단순화되었습니다.

Component transformations

사용자는 이제 여러 하위 구성 요소로 구성된 구성 요소에 회전, 변환 및 스케일링 변환을 적용하여 복잡한 형상 어셈블리 설정 프로세스를 단순화 할 수 있습니다. GMO (General Moving Object) 구성 요소의 경우, 이러한 변환을 구성 요소의 대칭 축과 정렬되도록 신체에 맞는 좌표계에 적용 할 수 있습니다.

Changing the number of threads at runtime

시뮬레이션 중에 솔버가 사용하는 스레드 수를 변경하는 기능이 런타임 옵션 대화 상자에 추가되어 사용 가능한 스레드를 추가하거나 다른 태스크에 자원이 필요한 경우 스레드 수를 줄일 수 있습니다.

Probe-controlled heat sources

활성 시뮬레이션 제어가 형상 구성 요소와 관련된 heat sources로 확장되었습니다. 히스토리 프로브로 열 방출을 제어 할 수 있습니다.

Time-dependent temperature at sources     

질량 및 질량 / 운동량 소스의 유체 온도는 이제 테이블 입력을 사용하여 시간의 함수로 정의 할 수 있습니다.

Emissivity coefficients

공극으로의 복사 열 전달을위한 방사율 계수는 이제 사용자가 방사율과 스테판-볼츠만 상수를 지정하도록 요구하지 않고 직접 정의됩니다. 후자는 이제 단위 시스템을 기반으로 솔버에 의해 자동으로 설정됩니다.

Output

  • 등속 필드 솔버 옵션을 사용할 때 유량 속도를 선택한 데이터 로 출력 할 수 있습니다 .
  • 벽 접착력으로 인한 지오메트리 구성 요소의 토크 는 기존 벽 접착력의 출력 외에도 일반 이력 데이터에 별도의 수량으로 출력됩니다.
  • 난류 모델 출력이 요청 될 때 난류 에너지 및 소산과 함께 전단 속도 및 y +가 선택된 데이터로 자동 출력됩니다 .
  • 공기 유입 모델 출력에 몇 가지 수량이 추가되었습니다. 자유 표면을 포함하는 모든 셀에서 혼입 된 공기 및 빠져 나가는 공기의 체적 플럭스가 재시작 및 선택된 데이터로 출력되어 사용자에게 공기가 혼입 및 탈선되는 위치 및 시간에 대한 자세한 정보를 제공합니다. 전체 계산 영역 및 각 샘플링 볼륨 에 대해이 두 수량의 시간 및 공간 통합 등가물 이 일반 히스토리 로 출력됩니다.
  • 솔버의 출력 파일 flsgrf 의 최종 크기 는 시뮬레이션이 끝날 때보 고됩니다.
  • 2 유체 시뮬레이션의 경우, 기존의 출력 수량 유체 체류 시간 및 유체 가 이동 한 거리는 이제 유체 # 1 및 # 2와 유체의 혼합물에 대해 별도로 계산됩니다.
  • 질량 입자의 경우 각 종의 총 부피와 질량이 계산되어 전체 계산 영역, 샘플링 볼륨 및 플럭스 표면에 대한 일반 히스토리 로 출력되어 입자 종 수에 대한 현재 출력을 보완합니다.
  • 예를 들어 사용자가 가스 미순환을 식별하고 연료 탱크의 환기 시스템을 설계하는 데 도움이 되도록 마지막 국부적 가스 압력이 옵션 출력량으로 추가되었습니다. 이 양은 유체가 채워지기 전에 셀의 마지막 간극 압력을 기록하며, 단열 버블 모델과 함께 사용됩니다.

New Customizable Source Routines

사용자 정의 가능한 새로운 소스 루틴이 추가되었으며 사용자의 개발 환경에서 액세스 할 수 있습니다.

소스 루틴 이름설명
cav_prod_cal캐비 테이션 생산 및 확산 속도
sldg_uset슬러지 정착 속도
phchg_mass_flux증발 및 응축에 의한 질량 흐름
flhtccl유체#1과#2사이의 열 전달 계수
dsize_cal2상 유동에서의 동적 낙하 크기 모델의 충돌 및 이탈율
elstc_custom.점탄성 유체에 대한 응력 방정식의 소스 용어

Brand New User Interface

FLOW-3D의 사용자 인터페이스가 완전히 재설계되어 사용자의 작업 흐름을 획기적으로 간소화하는 최신의 타일 구조를 제공합니다.

Dock widgets 설정

Physics, Fluids, Mesh 및 FAVOR ™를 포함한 모든 설정 작업이 형상 창 주위의 dock widgets으로 변환되어 모델 설정을 단일 탭으로 압축 할 수 있습니다. 이 전환을 통해 이전 버전의 복잡한 트리가 훨씬 깔끔하고 효율적인 메뉴 표시로 바뀌어 모델 설정 탭을 떠나지 않고도 모든 매개 변수에 쉽게 액세스 할 수 있습니다.

New Model Setup icons
With our new Model Setup design comes new icons, representing each step of the setup process.
New Physics icons
Our Physics icons are designed to be easily differentiated from one another at a glance, while providing clear visual representation of each model’s purpose and use.

RSS feed

새 RSS 피드부터 FLOW-3D v12.0 의 시뮬레이션 관리자 탭이 개선되었습니다 . FLOW-3D 를 시작하면 사용자에게 Flow Science의 최신 뉴스, 이벤트 및 블로그 게시물이 표시됩니다.

Configurable simulation monitor

시뮬레이션을 실행할 때 중요한 작업은 모니터링입니다. FLOW-3Dv12.0에서는 사용자가 시뮬레이션을 더 잘 모니터링할 수 있도록 Simulation Manager의 플로팅 기능이 향상되었습니다. 사용자는 시뮬레이션 런타임 그래프를 통해 모니터링할 사용 가능한 모든 일반 기록 데이터 변수를 선택하고 각 그래프에 여러 변수를 추가할 수 있습니다. 이제 런타임에서 사용할 수 있는 일반 기록 데이터는 다음과 같습니다.

  • 최소/최대 유체 온도
  • 프로브 위치의 온도
  • 유동 표면 위치에서의 유량
  • 시뮬레이션 진단(예:시간 단계, 안정성 한계)
Runtime plots of the flow rate at the gates of the large dam / Large dam with flux surfaces at the gates

Conforming mesh visualization

사용자는 이제 새로운 FAVOR ™ 독 위젯을 통해 적합한 메쉬 블록을 시각화 할 수 있습니다 .

Large raster and STL data

데이터를 처리하는 데 걸리는 시간으로 인해 큰 형상 데이터를 처리하는 것은 어려울 수 있습니다. 대형 지오메트리 데이터를 처리하는 데 여전히 상당한 시간이 소요될 수 있지만 FLOW-3D는 이제 이러한 대형 데이터 세트를 백그라운드 작업으로로드하여 사용자가 데이터를 처리하는 동안 완벽하게 응답하고 중단없는 인터페이스에서 계속 작업 할 수 있습니다.

Spontaneous capillary flow (자발적인 모세관 유동)

  • 미세 홈은 모세관 효과를 사용하여 유동을 유도할 수 있음
  • 전혈과 같은 고점도 유체도 가능
  • 자발적인 모세관 유동의 시작 조건?
    – sinα=cosθ (여기서, α는 V-홈 채널의 반각이고 θ는 접촉각)
    – For α=14.5도, θ<75.5도
    – 혈액의 경우, θ<75.5도 에서 장치가 작동해야함


  • FLOW-3D의 강력한 표면 장력 모델을 사용하여 모세관 효과를 캡처함
  • 수치적 예측과 실험간의 탁월한 일치

Micro/Biofluidics with FLOW-3D – Liquid handling (액체 취급)

나노리터 물방울의 정밀 분배

  • 섬세하고 정확한 분석
  • 원액의 소비를 정확하게 제어할 수 있음
  • 유체 특성/동역학에 기반한 공정 파라미터
    – 자유 표면 흐름의 복잡성을 고려
    – 자연스러운 모세관 중심 불안정을 고려
    – 씨닝 및 핀치 오프를 고려

방울의 형성 및 분리

  • 모세관, 관성, 점성 및 중력의 복잡한 상호 작용
  • 표면 장력과 점성력이 “핀치 오프”를 넘어가면 분리가 발생
  • FLOW-3D는 예측할 수 있음
    – 근본적인 응력
    – 확장된 유동장
    – 희석된 액체 필라멘트 내의 유동장을 시각화

미세 방울의 병합을 위한 유전영동

  • 유전영동력은 불균일한 전기장(일반적으로 AC전기장)에서 움직임을 유발함
  • 나노리터 유체 또는 나노 규모 입자의 특성을 다루고 처리하는데 사용

유동 집중

  • 다유체 계면 장력 파악
  • 방울 형성의 세부 사항 확인
  • 미세 방울의 진화 파악 (형태/크기)

FLOW-3D 교육 안내

education_banner

HIGH-END TOP CLASS
FLOW-3D CFD EDUCATION


(주)에스티아이씨앤디에서는 FLOW-3D 제품군의 사용자 교육을 지원하고 있습니다. 홈페이지에 안내되어 있는 교육 일정과 교육신청 절차를 참고하시어 교육을 받으실 수 있습니다.

  • 교육 과정명 : 수리 분야

댐, 하천의 여수로, 수문 등 구조물 설계 및 방류, 월류 등 흐름 검토를 하기 위한 유동 해석 방법을 소개하는 교육 과정입니다. 유입 조건(수위, 유량 등)과 유출 조건에 따른 방류량 및 유속, 압력 분포 등 유체의 흐름을 검토를 할 수 있도록 관련 예제를 통해 적절한 기능을 습득하실 수 있습니다.

  • 교육 과정명 : 수처리 분야

정수처리 및 하수처리 공정에서 각 시설물들의 특성에 맞는 최적 운영조건 검토 및 설계 검토을 위한 유동해석 방법을 소개하는 교육 과정입니다. 취수부터 시작하여 혼화지, 분배수로, 응집지, 침전지, 여과지, 정수지, 협기조, 호기조, 소독조 등 각 공정별 유동 특성을 검토하기 위한 해석 모델을 설정하는 방법에 대해 알려드립니다.

  • 교육 과정명 : 주조 분야

주조 분야 사용자들이 쉽게 접근할 수 있도록 각 공정별로 해석 절차 및 해석 방법을 소개하는 교육 과정입니다. 고압다이캐스팅, 저압다이캐스팅, 경동주조, 중력주조, 원심주조, 정밀주조 등 주조 공법 별 관련 예제를 통해 적절한 기능들을 습득할 수 있도록 도와 드립니다.

  • 교육 과정명 : Micro/Bio/Nano Fluidics 분야

점성력 및 모세관력 같은 유체 표면에 작용하는 힘이 지배적인 미세 유동의 특성을 정확하게 표현할 수 있는 해석 방법에 대해 소개하는 교육 과정입니다. 열적, 전기적 물리 현상을 구현할 수 있도록 관련 예제와 함께 해석 방법을 알려드립니다.

  • 교육 과정명 : 코팅 분야 과정

코팅 공정에 따른 코팅액의 두께, 균일도, 유동 특성 분석을 위한 해석 방법을 소개하는 교육 과정입니다. Slide coating, Dip coating, Spin coating, Curtain coating, Slot coating, Roll coating, Gravure coating 등 각 공정별 예제와 함께 적절한 기능을 습득하실 수 있도록 도와 드립니다.

  • 교육 과정명 : 레이저 용접 분야

레이저 용접 해석을 하기 위한 물리 모델과 용접 조건들을 설정하는 방법에 대해 소개하는 교육 과정입니다. 해석을 통해 용접 공정을 최적화할 수 있도록 관련 예제와 함께 적절한 기능들을 습득할 수 있도록 도와 드립니다.

  • 교육 과정명 : 3D프린팅 분야 과정

Powder Bed Fusion(PBF)와 Directed Energy Deposition(DED) 공정에 대한 해석 방법을 소개하는 교육 과정입니다. 파우더 적층 및 레이저 빔을 조사하면서 동시에 금속 파우더 용융지가 적층되는 공정을 해석하는 방법을 관련 예제와 함께 습득하실 수 있습니다.

  • 교육 과정명 : 해안/해양 분야

해안, 항만, 해양 구조물에 대한 파랑의 영향 및 유체의 수위, 유속, 압력의 영향을 예측할 수 있는 해석 방법을 소개하는 과정입니다. 항주파, 슬로싱, 계류 등 해안, 해양, 에너지, 플랜트 분야 구조물 설계 및 검토에 필요한 유동해석을 하실 수 있는 방법을 알려드립니다. 각 현상에 대한 적절한 예제를 통해 기능을 습득하실 수 있습니다.

  • 교육 과정명 : 우주/항공 분야

항공기 및 우주선의 연료 탱크와 추진체 관리장치의 내부 유동, 엔진 및 터빈 노즐 내부의 유동해석을 하실 수 있도록 관련 메뉴에 대한 설명, 설정 방법을 소개하는 과정입니다. 경계조건 설정, Mesh 방법 등 유동해석을 위한 기본적인 내용과 함께 관련 예제를 통해 기능들을 습득하실 수 있습니다.

기타 고객 맞춤형 과정

상기 과정 이외의 경우 고객의 사업 업무 환경에 적합한 사례를 중심으로 맞춤형 교육을 실시합니다. 필요하신 부분이 있으시면 언제든지 교육 담당자에게 연락하여 협의해 주시기 바랍니다.

고객센터 및 교육 담당자

  • 전화 : 02)2026-0455, 02)2026-0450
  • 이메일 : flow3d@stikorea.co.kr

교육은 정해진 일정에 시행되는 정기 교육과 고객의 요청에 의해 시행되는 특별 교육이 있습니다. 특별 교육이 실시될 경우 홈페이지를 통해 사전 공지를 합니다.

1. 연간교육 일정
FLOW-3D 연간교육일정

2. 교육 내용 : FLOW-3D Basic
  1. FLOW-3D 소개 및 이론
    • FLOW-3D 소개  – 연혁, 특징 등
    • FLOW-3D 기본 개념
      • VOF
      • FAVOR
    • 해석사례 리뷰
  2. GUI 소개 및 사용법
    • 해석 모델 작성법  – 물리 모델 설정
      • 모델 형상 정의
      • 격자 분할
      • 초기 유체 지정
      • 경계 조건 설정
    • 해석 결과 분석 방법  – 해석 모델 설명
  3. 해석 모델 작성 실습
    • 해석 모델 작성 실습  – 격자 분할
      • 물리 모델 설정
      • 모델 형상 및 초기 조건 정의
      • 경계 조건 설정
      • 해석 과정 모니터링
      • 해석 결과 분석
    • 질의 응답 및 토의

3. 교육 과정 : FLOW-3D Advanced
  1. Physics Ⅰ
    • Density evaluation
    • Drift flux
    • Scalars
    • Sediment scour
    • Shallow water
  2. Physics Ⅱ
    • Gravity and non-inertial reference frame
    • Heat transfer
    • Moving objects
    • Solidification
  3. FLOW-3D POST (Post-processor)
    • FLOW-3D POST 소개
    • Interface Basics
    • 예제 실습
Education Banner
  • 교육 신청은 홈페이지의 교육 신청 창에서 최소 3일 전에 신청합니다.
  • 모든 교육과정은 신청 인원이 2인 이상일때 개설되며, 선착순 마감입니다.
  • 교육 신청을 완료하시면, 신청시 입력하신 메일주소로 교육 담당자가 확인 메일을 보내드립니다.
  • 교육 시간은 Basic : 오전10시~오후5시, Advanced : 오후1시30분~오후5시30분까지입니다.
  • 교육비 안내
    • FLOW-3D Basic (2일) : 기업 66만원, 학생 55만원
    • FLOW-3D Basic 레이저용접, 3D 프린팅(2일) : 기업 88만원, 학생 66만원
    • FLOW-3D Advanced (1일) : 기업 33만원, 학생 25만원
    • 상기 가격은 부가세 포함 가격입니다.
  • 교육비는 현금(계좌이체)로 납부 가능하며, 교재 및 중식이 제공됩니다.
  • 세금계산서 발급을 위해 사업자등록증 또는 신분증 사본을 함께 첨부하여 신청해 주시기 바랍니다.
  • 교육 종료 후 이메일로 수료증이 발급됩니다.
고객센터 및 교육 담당자
  • 전화 : 02)2026-0455, 02)2026-0450
  • 이메일 : flow3d@stikorea.co.kr
교육 장소 안내
  • 지하철 1호선/가산디지털단지역 (8번출구), 지하철 7호선/가산디지털단지역 (5번출구)
  • 우림라이온스밸리 B동 302호 또는 교육장
  • 당사 건물에 주차할 경우 무료 주차 1시간만 지원되오니, 가능하면 대중교통을 이용해 주시기 바랍니다.
오시는 길

금속 3D 프린팅 수치해석

금속 3D 프린팅 / 적층 제조 수치해석(FLOW-3D Weld/DEM)

Binder jetting

바인더 분사 시뮬레이션은 모세관 힘의 영향을 받는 파우더 베드의 바인더 확산 및 침투에 대한 통찰력을 제공합니다. 공정 매개 변수와 재료 특성은 증착 및 확산 공정에 직접적인 영향을 미친다.

Direct energy deposition

FLOW-3D의 Particle 모델을 사용하여 직접 에너지 축적 프로세스를 시뮬레이션 할 수도 있습니다. 고체 기판에 분말 주입 속도와 열유속 입사를 지정함으로써, 고체 입자는 용융풀을 통해 질량, 추진력 및 에너지를 추가할 수 있습니다. 다음 동영상에서는 용융풀을 통해 고체 금속 입자가 주입되고 이어서 기판에 용융풀응 응고시키는 과정이 관찰됩니다.


관련 기술자료

Effect of Y2O3 on microstructure

Hierarchical grain refinement during the laser additive manufacturing of Ti-6Al-4V alloys by the addition of micron-sized refractory particles

미크론 크기의 내화물 입자를 추가하여 Ti-6Al-4V 합금의 레이저 적층 제조중 계층적 입자 미세 조정 Xiang Wang, Lin-Jie Zhang, Jie Ning, ...
더 보기
Review on the evolution and technology of State-of-the-Art metal additive manufacturing processes

Review on the evolution and technology of State-of-the-Art metal additive manufacturing processes

최첨단 금속 적층 제조 공정의 진화 및 기술 검토 S.Pratheesh KumarS.ElangovanR.MohanrajJ.R.Ramakrishna Abstract Nowadays, the requirements of customers undergo dynamic changes ...
더 보기
Numerical simulation of slag movement from Marangoni flow for GMAW with computational fluid dynamics Figures

Numerical simulation of slag movement from Marangoni flow for GMAW with computational fluid dynamics

전산 유체 역학을 사용하여 GMAW에 대한 Marangoni 흐름에서 슬래그 이동의 수치 시뮬레이션 Dae-WonChoaYeong-DoParkbMuralimohanCheepucaBusan Machinery Research Center, Korea Institute of Machinery ...
더 보기
Weld bead surface images showing the slag formation location for (a) wire 1 and (b) wire 2.

The effect of alloying elements of gas metal arc welding (GMAW) wire on weld pool flow and slag formation location in cold metal transfer (CMT)

가스 금속 아크 용접 (GMAW) 와이어의 합금 원소가 CMT (Cold Metal Transfer)에서 용접 풀 흐름 및 슬래그 형성 위치에 미치는 ...
더 보기
Figure 2. Experimental setups for the (a) Al/Cu overlap joint and (b) laser welding process.

Investigation on Laser Welding of Al Ribbon to Cu Sheet: Weldability, Microstructure, and Mechanical and Electrical Properties

알루미늄 리본과 구리 시트의 레이저 용접에 대한 조사 : 용접성, 미세 구조, 기계적 및 전기적 특성 Won‐Sang Shin 1,†, Dae‐Won ...
더 보기
Fig. 9 (a) Velocity field, keyhole profile, and breakage of the keyhole to form bubble and (b) 2D temperature and velocity field along the longitudinal section

A Numerical Study on the Keyhole Formation During Laser Powder Bed Fusion Process

Keyhole에 대한 수치적 연구 : 레이저 분말 중 형성 베드 퓨전 공정 Subin Shrestha1J.B. Speed School of Engineering,University of Louisville,Louisville, ...
더 보기
Figure 5.6 Experimental set-up equipped with high-speed camera system

COMPUTATIONAL FLUID DYNAMIC MODELLING OF LASER ADDITIVE MANUFACTURING PROCESS AND EFFECT OF GRAVITY

전산 유체 역학 레이저 첨가제 모델링 제조 공정 및 중력의 영향 A thesis submitted to The University of ManchesterFor the ...
더 보기
마란 고니 효과 있음

Ti-6Al-4V 금속 분말에 의한 선택적 레이저 용융법 수치 해석

Ti-6Al-4V 금속 분말에 의한 선택적 레이저 용융법 수치 해석 선택적 레이저 용융법(SLM: Selective Laser Melting)은 3D 프린팅 기술의 하나로 최근 ...
더 보기
CASE2-실험 결과와 FLOW-3D WELD에 의한 해석 결과와의 비교(단면 형상)

FLOW-3D WELD 용접 사례

FLOW-3D WELD를 이용한 용접 해석 사례를 소개합니다. 열전도 형 용접 (레이저)  두께가 다른 모재 맞대기하이브리드  레이저 / 아크 하이브리드깊이 용해 형 ...
더 보기
분석 모델 (위)과 실측 결과 (아래)

납땜(Soldering) 영역 제어

납땜 후 납땜 형상은 기본적으로 용융 상태에서 형성된 형상이 유지됩니다.특히 납땜의 미세화에 따라 용융 상태의 납땜 형상을 결정하는 요인으로서 표면 ...
더 보기

Home

FLOW-3D 는 세계에서 가장 어려운 CFD문제를 해결하는 소프트웨어로, 3차원 자유표면 해석 분야에서 널리 사용되는 최적의 수치해석 소프트웨어 입니다. 특히 자유표면(자유수면)을 가진 유동흐름을 정확하게 예측하는 분야에서는 타의 추종을 불허하는 정확성을 자랑합니다.

FLOW-3D 는 핵폭탄 개발 프로젝트로 유명한 미국 국립 연구소 LANL(LosAlamos National Laboratory)의 허트(C. W. Hirt) 박사가 새로운 자유표면 추적기법(free surface tracking method)인 VOF(Volume ofFluid) 방법을 연구 개발한 후, 수 많은 유동현상에 대한 물리 모델을 추가하고 성능을 개선하여, 설계 및 운영단계에서 사용되면서 엔지니어에게 귀중한 통찰력을 제공하는 세계적인 CFD 소프트웨어 입니다.

FLOW-3D 는 정확한 자유표면 추적, 압축성/비압축성 유동, 층류/난류, 열전달(전도, 대류, 복사), 점성발열, 상변화(응고,증발)/공동현상, 표면장력, 다상유동, 물질확산, 자연대류/밀도류, 뉴턴/비뉴턴유체, 틱소트로피, 다공성매질, 가속도계/관성계, 입자추적, 전기섭동/전기삼투압/주울발열, 열모세관현상 등 수많은 물리 모델을 제공합니다.

수치해석과 관련하여 궁금하신 사항은 언제든지 부담없이 문의 해주십시오.
감사합니다.

 
FLOW-3D Product
FLOW-3D 는 당사의 주력 제품으로 강력하고 매우 정확한 다중 물리 전산 유체역학(CFD) 패키지 프로그램…
Read More >>
FLOW-3D  WELDAM은 레이저 용접 공정 제어를 통해 다공성, 열 영향 영역을 최소화하고, 미세 구조 변화를 제어 할 수 있습니다. …
Read More >>
FLOW-3D Cast 는 다양한 주조 공정의 충전 및 응고, 결함 분포 예측이 가능한 3차원 유동해석 프로그램…
Read More >>
FLOW-3D 제품에 대한 기술자료와 이론 및 논문 등 다양한 기술자료를 제공합…

 

Read More >>

신규소식 기술자료
图 6 各流量监测断面位置

Study on the downstream impact of the numerical simulation of tailings library based on FLOW-3D

Jiahao Hu1, Chengwei Na1 and Yi Wang1 Published under licence by IOP Publishing LtdIOP Conference Series: Earth and Environmental Science, Volume 643, 2020 6th International Conference on Hydraulic and Civil Engineering 11-13 ...
자세한 내용 보기
Interaction between oblique waves and arc-shaped breakwater

Interaction between oblique waves and arc-shaped breakwater: Wave action on the breakwater and wave transformation behind it

XinyuHanaShengDongaYizhiWangbaCollege of Engineering, Ocean University of China, Qingdao, 266100, ChinabShandong Harbour Engineering Group Co., Ltd., Rizhao, 276826, China Highlights Interaction of oblique waves and the arc-shaped breakwater was simulated. Wave ...
자세한 내용 보기
Gambar 1-6 : Banjir Di Jalan Pagarsih Akibat Banjir Sungai Citepus

연구 제안-혁신적 생산적 연구 자금 (RISPRO) INVITATION

홍수 재해 완화 기술 정책 개발 반둥 분지 지역 이 제안된 연구의 제목, 즉 반둥 분지 지역의 홍수 재해 완화를 위한 기술 정책의 준비에서 알 수 있듯이 연구 대상 지역은 ...
자세한 내용 보기
Fig. 2. Semi-Lagrangian cellwise advection. (a) Forward advection scheme, (b) Backward advection scheme.

Three-dimensional cellwise conservative unsplit geometric VOF schemes

3차원 셀별 보수 미분할 기하학적 VOF 체계 Raphaël Comminal, JonSpangenberg Abstract This work presents two unsplit geometric VOF schemes that extend the two-dimensional cellwise conservative unsplit (CCU) scheme [Comminal et ...
자세한 내용 보기
Effect of Y2O3 on microstructure

Hierarchical grain refinement during the laser additive manufacturing of Ti-6Al-4V alloys by the addition of micron-sized refractory particles

미크론 크기의 내화물 입자를 추가하여 Ti-6Al-4V 합금의 레이저 적층 제조중 계층적 입자 미세 조정 Xiang Wang, Lin-Jie Zhang, Jie Ning, Sen Li, Liang-Liang Zhang, Jian LongState Key Laboratory for Mechanical ...
자세한 내용 보기
Figure 47: The course of the level on the physical model [22]

NUMERICAL MODELLING OF FLOW IN SPILLWAY

Author Svoboda, JiříContributors Jandora, Jan (advisor); Holomek, Petr (referee) Abstract 이 학위 논문의 주제는 Boskovice 상수도의 안전 배수로에서 유량 수치 모델링 솔루션입니다. 디플로마 논문의 소개에서는 기본 오버플로를 일반적으로 설명하고 모양과 ...
자세한 내용 보기
Figure 1. Steady-state shear stress a as a function of shear rate y in Sn-Pb alloy [10).

Numerical Modelling of Semi-Solid Flow under Processing Conditions

처리조건에서의 반고체유동의 수치모델링 David H. Kirkwood and Philip J. WardDepartment of Engineering Materials, University of Sheffield, Sheffield I UK Keywords: semi-solid alloys, thixotropy, flow modelling. Abstract During the industrial process ...
자세한 내용 보기
Review on the evolution and technology of State-of-the-Art metal additive manufacturing processes

Review on the evolution and technology of State-of-the-Art metal additive manufacturing processes

최첨단 금속 적층 제조 공정의 진화 및 기술 검토 S.Pratheesh KumarS.ElangovanR.MohanrajJ.R.Ramakrishna Abstract Nowadays, the requirements of customers undergo dynamic changes and industries are heading towards the manufacturing of customized end-user ...
자세한 내용 보기
The Simulation of Droplet Impact on the Super-Hydrophobic Surface with Micro-Pillar Arrays Fabricated by Laser Irradiation and Silanization Processes

The simulation of droplet impact on the super-hydrophobic surface with micro-pillar arrays fabricated by laser irradiation and silanization processes

레이저 조사 및 silanization 공정으로 제작된 micro-pillar arrays를 사용하여 초 소수성 표면에 대한 액적 영향 시뮬레이션 ZhenyanXiaaYangZhaoaZhenYangabcChengjuanYangabLinanLiaShibinWangaMengWangabaSchool of Mechanical Engineering, Tianjin University, Tianjin, 300054, ChinabKey Laboratory of Mechanism Theory and ...
자세한 내용 보기
Modeling of contactless bubble–bubble interactions in microchannels with integrated inertial pumps

Modeling of contactless bubble–bubble interactions in microchannels with integrated inertial pumps

통합 관성 펌프를 사용하여 마이크로 채널에서 비접촉식 기포-기포 상호 작용 모델링 Physics of Fluids 33, 042002 (2021); https://doi.org/10.1063/5.0041924 B. Hayesa), G. L. Whitingb), and R. MacCurdyc) ABSTRACT In this study, ...
자세한 내용 보기
Figure 4. Calculate and simulate the injection of water in a single-channel injection chamber with a nozzle diameter of 60 μm and a thickness of 50 μm, at an operating frequency of 5 KHz, in the X-Y two-dimensional cross-sectional view, at 10, 20, 30, 40 and 200 μs.

DNA Printing Integrated Multiplexer Driver Microelectronic Mechanical System Head (IDMH) and Microfluidic Flow Estimation

DNA 프린팅 통합 멀티플렉서 드라이버 Microelectronic Mechanical System Head (IDMH) 및 Microfluidic Flow Estimation by Jian-Chiun Liou 1,*,Chih-Wei Peng 1,Philippe Basset 2 andZhen-Xi Chen 11School of Biomedical Engineering, Taipei Medical ...
자세한 내용 보기
Figure 6. Maximum inundation field in simulations with (a) no barrier on the seawall (red line), (b) a 1 m barrier across the entire sea wall, and (c) a 1.7 m barrier partially installed on the seawall.

Storm surge inundation simulations comparing three-dimensional with two-dimensional models based on Typhoon Maemi over Masan Bay of South Korea

Jae-Seol Shim†, Jinah Kim†, Dong-Chul Kim‡, Kiyoung Heo†, Kideok Do†, Sun-Jung Park ‡† Coastal Disaster Research Center,Korea Institute of Ocean Science &Technology, 426-744, Ansan, Gyeonggi,Koreajsshim@kiost.acjakim@kiost.ackyheo21@kiost.ackddo@kiost.ac‡ Technology R&D InstituteHyein E&C Co., ...
자세한 내용 보기
Fig. 2 Temperature distributions of oil pans (Cycling)

내열마그네슘 합금을 이용한 자동차용 오일팬의 다이캐스팅 공정 연구

A Study on Die Casting Process of the Automobile Oil Pan Using the Heat Resistant Magnesium Alloy 한국자동차공학회논문집 = Transactions of the Korean Society of Automotive Engineersv.17 no.3 = no.99 ...
자세한 내용 보기
Fig. 1.Schematic of wire feeding in a melting line.

Evaluation on the Efficiency of Cored Wire Feeding in Addition of Alloying Elements into Cu Melt

Bok-Hyun Kang*, Ki-Young KimKorea University of Technology and Education 코어드 와이어 피딩에 의한 Cu 용탕에의 합금 첨가 시 효율 평가 Abstract To add alloying elements into a pure copper melt, ...
자세한 내용 보기
Fig. 3. Mesh and depth map for the storm surge model of ADCSWAN model.

ADCSWAN과 FLOW-3D 모델을 이용한 태풍 차바 내습 시 부산 마린시티의 침수범람 재현

최흥배․엄호식†․박종집․강태욱*, *** ㈜지오시스템리서치 선임, ** ㈜지오시스템리서치 책임, **** 부경대학교 박사 Reproduction of Flood Inundation in Marine City, Busan During the Typhoon Chaba Invasion Using ADCSWAN and FLOW-3D Models 요 약 ...
자세한 내용 보기
Figure 14. Results of 3D flow simulation for V = 0.82 m/s: (a) perspective view of velocity field on the free surface, (b) top view of velocity field on the free surface, (c) velocity field in the horizontal plane at half-length section of the rotor, and (d) velocity field in the rotor symmetry plane.

Experimental Method for the Measurements and Numerical Investigations of Force Generated on the Rotating Cylinder under Water Flow

by Teresa Abramowicz-Gerigk 1,*,Zbigniew Burciu 1,Jacek Jachowski 1,Oskar Kreft 2,Dawid Majewski 3,Barbara Stachurska 3,Wojciech Sulisz 3 andPiotr Szmytkiewicz 3 1Faculty of Navigation, Gdynia Maritime University, 81-225 Gdynia, Poland2AREX Ltd., 81-212 ...
자세한 내용 보기
A photo of HeMOSU-1.

FLOW-3D를 이용한 해상 자켓구조물 주변의 세굴 수치모의 실험

Numerical Simulation Test of Scour around Offshore Jacket Structure using FLOW-3D J Korean Soc Coast Ocean Eng. 2015;27(6):373-381Publication date (electronic) : 2015 December 31doi : https://doi.org/10.9765/KSCOE.2015.27.6.373Dong Hui Ko*, Shin Taek ...
자세한 내용 보기
Numerical simulation of slag movement from Marangoni flow for GMAW with computational fluid dynamics Figures

Numerical simulation of slag movement from Marangoni flow for GMAW with computational fluid dynamics

전산 유체 역학을 사용하여 GMAW에 대한 Marangoni 흐름에서 슬래그 이동의 수치 시뮬레이션 Dae-WonChoaYeong-DoParkbMuralimohanCheepucaBusan Machinery Research Center, Korea Institute of Machinery and Materials, 48, Mieumsandan 5-ro 41beon-gil, Gangseo-gu, Busan 46744, Republic ...
자세한 내용 보기
Weld bead surface images showing the slag formation location for (a) wire 1 and (b) wire 2.

The effect of alloying elements of gas metal arc welding (GMAW) wire on weld pool flow and slag formation location in cold metal transfer (CMT)

가스 금속 아크 용접 (GMAW) 와이어의 합금 원소가 CMT (Cold Metal Transfer)에서 용접 풀 흐름 및 슬래그 형성 위치에 미치는 영향 Md. R. U. Ahsan1,3, Muralimohan. Cheepu2, Yeong-Do Park* 2,31Department ...
자세한 내용 보기
Figure 2. Experimental setups for the (a) Al/Cu overlap joint and (b) laser welding process.

Investigation on Laser Welding of Al Ribbon to Cu Sheet: Weldability, Microstructure, and Mechanical and Electrical Properties

알루미늄 리본과 구리 시트의 레이저 용접에 대한 조사 : 용접성, 미세 구조, 기계적 및 전기적 특성 Won‐Sang Shin 1,†, Dae‐Won Cho 2,†, Donghyuck Jung 1, Heeshin Kang 3, Jeng O ...
자세한 내용 보기
Fig. 6: Proposed Pattern Layout

Casting Defect Analysis on Caliper Bracket using Mold flow Simulation

금형 흐름 시뮬레이션을 사용한 캘리퍼 브래킷의 주조 결함 분석 Abstract 이 작업에서는 컴퓨터 보조 주조 시뮬레이션 기술을 사용하여 Green sand 주조의 모래, 기계 및 설계 관련 결함을 분석합니다. 자동차 브레이크 ...
자세한 내용 보기
Figure 4. Structure of artificial neural network [37]

Turbulent Flow Modeling at Tunnel Spillway Concave Bends and Prediction of Pressure using Artificial Neural Network

터널 배수로 오목 굴곡에서 난류 유동 모델링 인공 신경망을 이용한 압력 예측 및 예측 Zeinab Bashari Moghaddam 1Hossein Mohammad Vali Samani2Seyed Habib Mousavi Jahromi 3 Abstract 터널 배수로는 높은 자유 ...
자세한 내용 보기
Fig. 1 Fixed staff gauge for head measurement at the upstream side of the Yuanshanzi Flood Diversion Work in the Keelung River, Taiwan

Velocity distribution and discharge calculation at a sharp-crested weir

Shun-Chung Tsung • Jihn-Sung Lai •Der-Liang Young sharp-crested weir에서 속도 분포 및 배출 계산 개방 수로의 harp-crested 위어는 수두-방류 관계를 통해 방류를 계산하는데 유용한 장치입니다. 그러나 수위 측정 사이트와 배출 ...
자세한 내용 보기
planar representation (cross-section at tank centre).

Analysis of cryogenic propellant behaviour in microgravity and low thrust environments*

미세 중력 및 저 추력 환경에서 극저온 추진체 거동 분석 M.F. Fisher, G.R. Schmidt and J.J. MartinNASA Marshall Space Flight Center, Huntsville, AL 35824, USA Abstract 우주선 비행 작업 (예 ...
자세한 내용 보기
Figure 2.1. Test Setup.The test setup consists of a clear plastic scale model tank attached to a rigid aluminum frame by three multi-axis load cells driven by a position-controlled servo hydraulic system.(Data acquisition cabling removed for clarity).

Coupled Simulation of Vehicle Dynamics and Tank Slosh. Phase 1 Report. Testing and Validation of Tank Slosh Analysis

Prepared byGlenn R. WendelSteven T. GreenRussell C. Burkey Abstract: 차량 동력학의 컴퓨터 시뮬레이션은 차량 설계에서 귀중한 도구가 되었다. 그러나 그들은 차량의 탱크에서 유체 슬로싱의 복잡한 역학을 정확하게 시뮬레이션할 수 없다 ...
자세한 내용 보기
Fig. 4 Current lines in the horizontal level in: a 0.70 and b 14 cm from the streambed in tandem pies

3D numerical simulation of flow field around twin piles

트윈 말뚝 주위의 유동장 3D 수치 시뮬레이션 Amini, A; Parto, AAAmini, A (reprint author), AREEO, Kurdistan Agr & Nat Resources Res & Educ Ctr, Sanandaj, Iran., 2017; 65 (6): 1243 ...
자세한 내용 보기
Figure 1. Sketch map of the port Laozi on Lake Hongze

FLOW-3D software for substantiation the layout of the port water area

항구 수역의 레이아웃을 입증하기 위한 FLOW-3D 소프트웨어 B Pan1 and N Belyaev21 Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russiapan3.b@edu.spbstu.ru Abstract 방파제 설계시 항만 내 수역 및 접안 ...
자세한 내용 보기
Capsule-type Vane Tank

Numerical Simulation Analysis of Liquid Transportation in Capsule-type Vane Tank under Microgravity

Microgravity 하에서 캡슐형 베인 탱크의 액체 수송에 대한 수치 시뮬레이션 분석 Abstract Li Yong-Qiang1,2 & Dong Jun-Yan1 & Rui Wei1Received: 17 June 2019 /Accepted: 4 December 2019#Springer Nature B.V. 2020 ...
자세한 내용 보기
Figure Top view of velocity distribution of tailings mortar

Study on discharge velocity of tailings mortar in dam break based on FLOW-3D

Jiahao Hu1, Chengwei Na1, Yi Wang1*College of Water Conservancy,Shenyang Agricultural University ,Shenyang,Liaoning, 110866, China*Corresponding author’s e-mail: yiwang@syau.edu.cn Abstract Tailings pond is used to store the tailings discharged from the mine ...
자세한 내용 보기
 

전체 기술자료로 바로가기

컨설팅 절차

컨설팅 절차

  • 해석 컨설팅을 저희에게 의뢰하시면, 상세한 상담 후 견적을 작성하여 보내 드립니다. 상담은 전화, 이메일, 방문 등의 방법으로 진행됩니다.
  • 계약이 체결된 후 수치해석을 위한 자료 및 데이터를 받아, 협의된 안으로 수치해석을 수행합니다.
  • 컨설팅 진행 과정 중에 수시로 해석 결과 및 진행 상황에 대해 연락 드리며, 변경, 수정 사항을 협의하여 반영할 수 있습니다.
  • 수치해석이 완료되면 최종 보고서를 작성하여 제출하며, 필요시 방문하여 결과를 상세히 설명 드립니다.
  • 수치해석 기술 전수가 포함된 계약일 경우, 최종 보고서 제출 이후에 기술 전수 교육을 진행합니다.
  • 모든 기술 자료는 대외비로 취급되며, 철저하게 보안을 유지해드립니다.

컨설팅 분야

수자원 분야

  • 댐체, 수문, 제반 구조물 안정성 검토
  • 댐, 여수로 유동 해석
  • 여수로 수위별 방류량 해석
  • 여수로 월류 및 수위 검토 해석
  • 발전소 취수로 유동 해석
  • 배수터널 방류향 해석
  • 취수탑 유입 유량 해석
  • 교각주위 세굴 해석
  • 수문 수차 유량 해석
  • 저수지 수위별 유동해석
  • 배수암거 부정류 해석
  • 저수지 연결 터널 유동 해석
  • 교각 유동 작용 힘 검토
  • 도수터널 통수 능력 해석
  • 부유사 확산 검토
  • 냉각수 취수로 유량 해석
  • 수문 유동 양상 분석
  • 배수터널 방류량 해석
  • 월류 수위별 유량 유속 해석

수처리 분야

  • 정수지 유동해석
  • 분배수로 유량분배 해석
  • 침전지 유동 및 유속 분포 해석
  • 반응조 농도 및 반응시간 해석
  • 응집지 유동해석
  • 하수처리시설 슬러지 농도 해석
  • DAF 응집제 농도 해석
  • 수조 최적 교반 해석
  • 여과지 유동해석
  • 혼화지 유동해석
  • 호기조 담체 거동해석
  • 수처리 구조물 유동 양상 분석
  • 하수처리시설 유동해석
  • 분말활성탄 접촉조 해석
  • PSBR 반응조 해석
  • 지하수 ICE RING 형성 해석
  • 절리면 모세관 열유동 해석
  • DAF 실증시설 부상조 해석
  • 착수정 유량 분배 해석

우주 항공분야

  • 발사체 탱크 슬로싱 댐핑 평가 해석
  • 항공기 비행 및 급유 시 연료 탱크 내부 유동 해석
  • 항공기 날개 연료 탱크 내부 유동 해석
  • 항공기 연료 탱크 내부 유동 해석
  • 추진체 관리 장치 내부 유동 해석
  • 엔진 및 터빈 노즐 내부 유동 및 캐비테이션 해석

자동차 분야

FLOW-3D POST Gears
  • 자동차 연료 탱크에 연료 주입 시 탱크 내부 유동 해석
  • 피스톤 쿨링젯 시스템 해석
  • 전착 도장 해석
  • 자동차 연료 주입구의 주입 유량별 유동 특성 분석
  • 기어 펌프의 로터 회전에 따른 오일 유동 양상 분석
  • 엔진 실린더 내 피스톤 운동과 배기가스 유동 패턴 해석
  • 베어링 내 윤활을 위한 오일의 유동 양상 해석

해양분야

  • 해양 컨테이너 연료 탱크 슬로싱 해석
  • 방파제 구조물 주변 유동 해석
  • 선박 운항에 따른 항주파 및 유동 특성 분석
  • 사석 방파제 등 구조물 주변 유동 해석
  • 진동수주형 파력 발전 구조물 최적화 모델 해석
  • 선박 및 부유체 계류 시 계류 안정성 및 계류력 해석
  • 발전소 부근 해역 온배수 영향 예측
  • 지진 해일에 의한 영향 해석

주조 분야

  • 고압다이캐스팅  충진 거동 및 응고 해석
  • 저압주조 충진 거동 및 응고 해석
  • 경동주조 충진 거동 및 응고 해석
  • 중력주조 충진 거동 및 응고 해석
  • 원심주조 충진 거동 및 응고 해석
  • 금형온도 분포 해석
  • 제품 및 금형 열응력, 변형 해석
  • 주조 공법 별 온도 분포, 산화물 분포 및 결함 분석
  • 금형 및 몰드 냉각방안 최적화 검토

Micro/Bio/Nano Fluidics 분야

  • Slit 및 Slot 코팅 해석
  • Roll 코팅 해석
  • Gravure / Gravure-offset 프린팅 해석
  • Curtain 코팅 해석
  • Multi-layer Slide 코팅 해석
  • 전기 삼투를 이용한 마이크로 펌프 전위 및 유동해석
  • 마이크로 채널 액적 생성 연속성 및 혼합 해석
  • 잉크젯 헤드 조건에 따른 잉크 분사 성능 해석
  • 열모데관 유동해석과 모세관 충진 해석
  • 유전 영동 현상을 이용한 액적 융합 해석

레이저 용접 분야

  • 이종재 레이저 용접 해석
  • 용접속도와 경사도에 따른 키홀 내부의 기공 거동 해석
  • 이종재의 레이저 용접 시 wobbling 해석
  • 레이저 용접 Melt Pool 거동 해석
  • 레이저 파워, 속도에 따른 balling 결함 영향 해석

HVAC System Designs

HVAC(난방, 냉방 및 환기)시스템 엔지니어가 고려해야 하는 최적 설계 배치에 대한 검토를 수행

발전소의 경우 대형(길이 90m, 너비 33m, 높이 26m)건물로 변압기, 전력선, 조명 등 열 발생 장비를 갖추고 있어서 여러가지 시설물의 상황을 고려할 수 있음

건물 내 공기를 올바르게 분배하고 적절한 쾌적한 온도를 확보하기 위해 건물 구조와 흡입그 크기 등의 검토 가능

수치해석 기술 컨설팅 안내

FLOW-3D Case Studies

수치해석 기술 컨설팅 안내

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD 컨설팅 서비스를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.
상담에는 비용은 전혀 들지 않습니다.

CFD는 엔지니어가 공기, 물 또는 모든 유체와의 상호 작용을 이해할 수 있게 하는 매우 효과적인 기술로 대부분의 유동현상에 해답을 제시 할 수있는 막대한 잠재력을 가지고 있습니다.
다양한 유체 흐름 현상이나 온도 및 열전달 분석 등 필요한 시나리오에 대한 맞춤 솔루션을 제공합니다.

당사에는 20년 이상 수치해석 연구에 전념하고 있는 전문 연구인력과 다양한 기술적 경험과 전문 시뮬레이션 기술을 제공하는 숙련된 기술컨설팅팀이 준비되어 있습니다.
귀하의 프로젝트 성공 가능성을 기술시연을 통해 제공 할 수 있습니다.
프로그램 소개나 자문이 필요하신 분들은 언제든지 아래 연락처로 문의하시기 바랍니다.

  • 전화 :   02-2026-0455
  • Email : flow3d@stikorea.co.kr

컨설팅 형태

수치해석 의뢰

  • 고객이 당면한 문제를 분석 /검토/협의 후, 가장 적절한 수치해석 방법을 수립합니다.
  • 주로 상호 협의된 설계안 및 해석 조건에 대해 수치해석을 수행하여 결과를 도출 분석, 검토합니다.
  • 설계 변경 인자 및 해석 횟수는 고객과 협의하여 진행합니다. 수치해석 결과를 분석 검토하여 설계에 반영하기 위한 의견을 제시하여 드립니다.

해석 대행 의뢰

  • 고객사에 해석 프로세스가 정립되어 있는 경우에 대해, 계산 장비와 수치해석 인력을 이용하여 해석 대행 및 해석 결과물을 제출합니다.

컨설팅 절차

  • 해석 컨설팅을 저희에게 의뢰하시면, 상세한 상담 후 견적을 작성하여 보내 드립니다. 상담은 전화, 이메일, 방문 등의 방법으로 진행됩니다.
  • 계약이 체결된 후 수치해석을 위한 자료 및 데이터를 받아, 협의된 안으로 수치해석을 수행합니다.
  • 컨설팅 진행 과정 중에 수시로 해석 결과 및 진행 상황에 대해 연락 드리며, 변경, 수정 사항을 협의하여 반영할 수 있습니다.
  • 수치해석이 완료되면 최종 보고서를 작성하여 제출하며, 필요시 방문하여 결과를 상세히 설명 드립니다.
  • 수치해석 기술 전수가 포함된 계약일 경우, 최종 보고서 제출 이후에 기술 전수 교육을 진행합니다.
  • 모든 기술 자료는 대외비로 취급되며, 철저하게 보안을 유지해드립니다.

주요 컨설팅 의뢰 분야

수자원 분야

  • 댐체, 수문, 제반 구조물 안정성 검토
  • 댐, 여수로 유동 해석
  • 여수로 수위별 방류량 해석
  • 여수로 월류 및 수위 검토 해석
  • 발전소 취수로 유동 해석
  • 배수터널 방류향 해석
  • 취수탑 유입 유량 해석
  • 교각주위 세굴 해석
  • 수문 수차 유량 해석
  • 저수지 수위별 유동해석
  • 배수암거 부정류 해석
  • 저수지 연결 터널 유동 해석
  • 교각 유동 작용 힘 검토
  • 도수터널 통수 능력 해석
  • 부유사 확산 검토
  • 냉각수 취수로 유량 해석
  • 수문 유동 양상 분석
  • 배수터널 방류량 해석
  • 월류 수위별 유량 유속 해석

수처리 분야

Wastewater Treatment Plant
Wastewater Treatment Plant
  • 정수지 유동해석
  • 분배수로 유량분배 해석
  • 침전지 유동 및 유속 분포 해석
  • 반응조 농도 및 반응시간 해석
  • 응집지 유동해석
  • 하수처리시설 슬러지 농도 해석
  • DAF 응집제 농도 해석
  • 수조 최적 교반 해석
  • 여과지 유동해석
  • 혼화지 유동해석
  • 호기조 담체 거동해석
  • 수처리 구조물 유동 양상 분석
  • 하수처리시설 유동해석
  • 분말활성탄 접촉조 해석
  • PSBR 반응조 해석
  • 지하수 ICE RING 형성 해석
  • 절리면 모세관 열유동 해석
  • DAF 실증시설 부상조 해석
  • 착수정 유량 분배 해석

우주 항공분야

  • 발사체 탱크 슬로싱 댐핑 평가 해석
  • 항공기 비행 및 급유 시 연료 탱크 내부 유동 해석
  • 항공기 날개 연료 탱크 내부 유동 해석
  • 항공기 연료 탱크 내부 유동 해석
  • 추진체 관리 장치 내부 유동 해석
  • 엔진 및 터빈 노즐 내부 유동 및 캐비테이션 해석

자동차 분야

FLOW-3D POST Gears
  • 자동차 연료 탱크에 연료 주입 시 탱크 내부 유동 해석
  • 피스톤 쿨링젯 시스템 해석
  • 전착 도장 해석
  • 자동차 연료 주입구의 주입 유량별 유동 특성 분석
  • 기어 펌프의 로터 회전에 따른 오일 유동 양상 분석
  • 엔진 실린더 내 피스톤 운동과 배기가스 유동 패턴 해석
  • 베어링 내 윤활을 위한 오일의 유동 양상 해석

해양분야

  • 해양 컨테이너 연료 탱크 슬로싱 해석
  • 방파제 구조물 주변 유동 해석
  • 선박 운항에 따른 항주파 및 유동 특성 분석
  • 사석 방파제 등 구조물 주변 유동 해석
  • 진동수주형 파력 발전 구조물 최적화 모델 해석
  • 선박 및 부유체 계류 시 계류 안정성 및 계류력 해석
  • 발전소 부근 해역 온배수 영향 예측
  • 지진 해일에 의한 영향 해석

주조 해석 분야

  • 고압다이캐스팅  충진 거동 및 응고 해석
  • 저압주조 충진 거동 및 응고 해석
  • 경동주조 충진 거동 및 응고 해석
  • 중력주조 충진 거동 및 응고 해석
  • 원심주조 충진 거동 및 응고 해석
  • 금형온도 분포 해석
  • 제품 및 금형 열응력, 변형 해석
  • 주조 공법 별 온도 분포, 산화물 분포 및 결함 분석
  • 금형 및 몰드 냉각방안 최적화 검토

Micro/Bio/Nano Fluidics 분야

  • Slit 및 Slot 코팅 해석
  • Roll 코팅 해석
  • Gravure / Gravure-offset 프린팅 해석
  • Curtain 코팅 해석
  • Multi-layer Slide 코팅 해석
  • 전기 삼투를 이용한 마이크로 펌프 전위 및 유동해석
  • 마이크로 채널 액적 생성 연속성 및 혼합 해석
  • 잉크젯 헤드 조건에 따른 잉크 분사 성능 해석
  • 열모데관 유동해석과 모세관 충진 해석
  • 유전 영동 현상을 이용한 액적 융합 해석

레이저 용접 분야

  • 이종재 레이저 용접 해석
  • 용접속도와 경사도에 따른 키홀 내부의 기공 거동 해석
  • 이종재의 레이저 용접 시 wobbling 해석
  • 레이저 용접 Melt Pool 거동 해석
  • 레이저 파워, 속도에 따른 balling 결함 영향 해석

공기/열 흐름 분야 (HVAC System Designs)

HVAC(난방, 냉방 및 환기)시스템 엔지니어가 고려해야 하는 최적 설계 배치에 대한 검토를 수행

발전소의 경우 대형(길이 90m, 너비 33m, 높이 26m)건물로 변압기, 전력선, 조명 등 열 발생 장비를 갖추고 있어서 여러가지 시설물의 상황을 고려할 수 있음

건물 내 공기를 올바르게 분배하고 적절한 쾌적한 온도를 확보하기 위해 건물 구조와 흡입그 크기 등의 검토 가능

Computational Analysis of Drop Formation and Detachment

Computational Analysis of Drop Formation and Detachment

Introduction and Problem Statement

신속, 반복, 작은 물방울의 생성 및 증착, 작은 형상의 프린팅 또는 패터닝 (예 : l = 10-3-1 mm), 스프레이로  균일한 두께의 박막 형성은 다양한 산업에 매우 중요합니다(1-5). 액체 이동과 액적 형성 / 증착 공정은 복잡한 자유 표면 흐름, 자연적인 모세관운동 형성, thinning, pinch-off를 수반한다 (1-5). 단순한 뉴턴 및 비탄성 유체에 대해 액적 생성 및 액적 이동을 분석하기위한 실험적, 이론적 및 1 차원 시뮬레이션 연구가 진행되었지만 프린팅 또는 패터닝에 대한 기계론적인 이해는 여전히 과제로 남아 있습니다. 현재의 계산에 대한 주된 목표는 뉴턴 유체의 pinch-off에 대한 기계론적 이해를 얻기 위해 FLOW-3D에 내장된 VOF(volume-of-fluid) 접근법으로 시험하는 것입니다. 전산해석은 모세관, 관성, 점성 응력의 복잡한 상호 작용을 포착하여 자기유사 모세관의 thinning and pinch-off를 결정합니다. 뉴턴 유체의 물방울 형성 ​​및 분리현상은  전산해석으로부터 얻어진 자기유사 모세관현상 이론, 보편적인 축소화 기법인 1D 시뮬레이션 (1-7)과 실험 (1, 2, 8-12)을 이용하여 설명될 수 있음을 보여준다. 이러한 우리가 진행한 원형흐름 시뮬레이션은 유한한 시간의 비선형 역학, 위성 낙하현상, 복잡한 형상의 프린팅과 같이 어려운 전산해석의 기반이 될 것 입니다.

방울 형성의 전산 분석
그림 1 : FLOW-3D를 사용하여 시뮬레이션 한 저점도 유체의 드롭 형성 및 분리에 대한 전산해석 : (a) 5개의 저점도 유체에 대한 물방울의 necking에 대한 반경이 시간변화에 따라 표시됩니다. 물방울 necking의 반지름이 오른쪽에서 왼쪽으로 시간에 따른 전개를 보여줍니다. 마찬가지로 스냅 샷은 necking의 반경이 오른쪽에서 왼쪽으로 줄어듭니다. 속도의 크기 (단위 : cm/s) 와 화살표의 방향에 대한 컬러 맵을 사용하면 변형장을 결정할 수 있으며 Fluid 5 (표 1 참조)의 경우에는 순식간에 신장이됩니다. 이미지 II에 캡처 된 pinch-off 하기 전에 형성된 원추형 necking은 실험을 통해 얻은 necking 모양과 유사합니다.

Modeling Approach and Parameter Space

표면 장력 및 중력 모델을 적용한 FLOW-3D 에서 균일한 메쉬 크기를 사용하여 노즐에서 드롭 형성 및 분리에 대한 시뮬레이션을 수행하였습니다. 유한 체적의 유체를 떨어뜨리거나 분리하는 일은 물방울의 성장과 드롭, 노즐에 연결되는 모세관 현상, 관성, 점도 및 중력에 대한 상호 작용을 수반합니다. 시뮬레이션에서 스테인레스 강 노즐 ( {{D} _ {0}} = 2 {{R} _ {0}} = 1.7 \, \ text {mm}) 에서 유한 체적의 뉴턴 유체가 발생합니다. 표면 장력이 중력을 겪으면 새로 형성된 액적 분리가 발생합니다 (mg> 2 \ pi \ sigma {{R} _ {0}}). 시뮬레이션은 유체점도의 영향을 설명하기 위해 두 그룹으로 나누어져 있습니다: 저점도 유체 (글리세롤 함량이 40 % 미만인 물과 글리세롤/물 혼합물) 및 점도가 높은 유체 (예 : 글리세롤과 글리세롤/물 혼합물 점도 > 100x 물 점도). 두 그룹의 유체 특성은 각각 표 1과 2에 나와 있습니다.

계산 분석 드롭 형성 저점도

그림 2 : FLOW-3D를 사용하여 시뮬레이션 한 저점도 유체의 드롭형성 및 분리에 대한 전산 해석 : 반경 플롯에서 4개의 고점도 뉴톤유체에 대해 necking 반경을 시간변화에 따라 표시합니다. 낙하 분리 중 모세관 현상이 스냅 샷으로 표시됩니다. 컬러 맵은 Fluid 8의 속도 크기 (단위 : cm/s)의 변화를 포착합니다 (표2 참조). 화살표는 성장하는 물방울과 얇아지는 물방울내에서 흐름방향을 나타냅니다. FLOW-3D 시뮬레이션으로 얻은 necking 모양은 고점도의 뉴턴유체에 대한 특징인 원통형 유체요소로 이어집니다.

 

<표 1 : FLOW-3D를 사용하여 시뮬레이션 된 저점도 유체의 특성>
Fluid Property Fluid 1 Fluid 2 Fluid 3 Fluid 4 Fluid 5
Viscosity [Pa · s] 0.05 0.02 0.01 0.0075 0.005
Surface Tension  [mN / m] 68 68 68 68 68
Density [g / cm 3 ] 1 1 1 1 1
Ohnesorge Number 0.21 0.08 0.04 0.03 0.021
 저점도 유체 (표 1의 유체 2) 가 노즐에서 떨어지는 것을 시뮬레이션 합니다. 색상변수는 속도크기 (단위 : cm / s)이며 속도벡터가 표시됩니다.

 

<표 2 : FLOW-3D를 사용하여 시뮬레이션 된 고점도 유체의 특성>
Fluid Property Fluid 6 Fluid 7 Fluid 8 Fluid 9
Viscosity [Pa · s] 1.5 0.8 0.5 0.25
Surface Tension  [mN / m ] 68 68 68 68
Density [g / cm 3 ] 1 1 1 1
Ohnesorge Number 6.24 3.33 2.08 1.04

고점도 유체 (표 2의 유체 8) 가 노즐에서 떨어지는 것을 시뮬레이션 합니다. 색상변수는 속도크기 (단위 : cm / s) 이며 속도 벡터가 표시됩니다.

Discussion of the Simulation Results

드롭 형성 및 분리는 표1과 표2에 열거 된 유체에 대해 FLOW-3D 를 사용하여 시뮬레이션 하였고, 시간 경과에 따른 necking 모양, 반경을 분석하였습니다. 물방울의 necking 모양과 저점도에서의 necking에 대한 역학(그림 1 참조)은 실험, 흐름 이론, 1D 시뮬레이션, 자기유사 관성에 대한 모세현상의 특성을 나타냅니다 (1, 2, 6, 7, 13) :

(1)  \ displaystyle \ frac {{R (t)}} {{{{R} _ {0}}}} \ approx 0.8 R {{{{왼쪽} {R} {0} 3}}} 오른쪽}) ^ {{{{frac {1} {3}}} {{왼쪽 {{{{왼쪽}}} {2} {3}}}}

여기서 R (t)가  necking의 순간 반경이고, R0는 노즐의 외부반경이며,  \ displaystyle \ sigma 는 표면 장력,  \ displaystyle \ rho 는 유체의 밀도 tC 는 pinch-off 시간이다. 마찬가지로, 이러한 더 높은 점도의 뉴턴유체에 대한 반경 변화데이터는 시간에 따른 반경의 감소를 나타내는 것이며,  Papageorgiou’s visco-capillary scaling (8, 9)은 아래의 식으로 표현된다.

(2)  \ {0 \} {} {} {} {} {} {} {} {} {} {} {} {} { } ({{t} _ {p}} - t)

모세관 속도(표면 장력과 점도의 비)의 측정 값은 McKinley와 Tripathi (8)에 의해 Capillary Break-Up Extensional Rheometer (CaBER)라고 불리는 상업적으로 이용 가능한 장비를 사용하여 얻은 값과 모세관 속도는 공칭 표면 장력과 점도를 사용하여 계산됩니다.

FLOW-3D 는 물방울의 necking부분을 속도 벡터로 시각화하여 유체의 흐름을 나타낼 수 있습니다. 또한, 이는 그림 1과 같이 전단, 확장을 겪은 후 얇아지는 물방울이 흐르는 과정의 순간을 결정할 수 있는 가능성을 줍니다. 추가로, 낮은 점도의 뉴턴유체는 높은 점도의 뉴턴 유체에 비해 질적으로 다른 거동을 보여준다(그림 2참조). 낮은 점도의 뉴턴 유체에 대한 necking 프로파일은 이론(6,13)에 따라 자기 유사성이 됩니다.

Conclusions, Outlook and Ongoing work

우리의 예비결과는 FLOW-3D 기반의 전산해석이 액적 형성과 탈착의 기초가 되는 프로토타입의 자유 표면흐름을 시뮬레이션하는데 사용될 수 있음을 보여줍니다 . 시뮬레이션된 반경변화 프로파일이 실험적으로 관찰된 높은 유체 및 이론적으로 예측된 유체인 스케일링 법칙 및 pinch-off dynamics과 일치하는 것을 발견하였습니다.

자주 사용되는 1D 또는 2D 모델과 달리 FLOW-3D 는 기본 응력 및 확장 유동장 (균일도 및 크기)의 강도와 얇은 액체 필라멘트 내 흐름에 대한 시각화를 나타낼 수 있습니다(그림1과 2 참조). 확장 유동장과 연관된 흐름 방향 속도 구배는 모세관현상이 나타나는 물방울의 얇은 부분 내에서 발생합니다. 유동학적으로 복잡한 유체에서 non Newtonian shear 및 신장, 점도뿐만 아니라 그외의 탄성 응력이 nonlinear pinch-off dynamics을 급격하게 변화시킵니다(2, 10-12). 우리는 현재 점탄성과 non-Newtonian 유동학을 사용하여 FLow-3D에 복합 유체의 처리 성능평가를 위한 강력한 연산 프로토콜을 개발하고 있습니다.

References

  1. J. Eggers, Nonlinear dynamics and breakup of free-surface flows. Rev. Mod. Phys. 69, 865-929 (1997).
  2. G. H. McKinley, Visco-elasto-capillary thinning and break-up of complex fluids. Rheology Reviews, 1-48 (2005).
  3. B. Derby, Inkjet Printing of Functional and Structural Materials: Fluid Property Requirements, Feature Stability, and Resolution. Annual Review of Materials Research 40, 395-414 (2010).
  4. O. A. Basaran, H. Gao, P. P. Bhat, Nonstandard Inkjets. Annual Review of Fluid Mechanics 45, 85-113 (2013).
  5. S. Kumar, Liquid Transfer in Printing Processes: Liquid Bridges with Moving Contact Lines. Annual Review of Fluid Mechanics 47, 67-94 (2014).
  6. R. F. Day, E. J. Hinch, J. R. Lister, Self-similar capillary pinchoff of an inviscid fluid. Phys. Rev. Lett. 80, 704-707 (1998).
  7. J. Eggers, M. A. Fontelos, Singularities: Formation, Structure, and Propagation. (Cambridge University Press, Cambridge, UK, 2015), vol. 53.
  8. G. H. McKinley, A. Tripathi, How to extract the Newtonian viscosity from capillary breakup measurements in a filament rheometer. J. Rheol. 44, 653-670 (2000).
  9. D. T. Papageorgiou, On the breakup of viscous liquid threads. Phys. Fluids 7, 1529-1544 (1995).
  10. J. Dinic, L. N. Jimenez, V. Sharma, Pinch-off dynamics and dripping-onto-substrate (DoS) rheometry of complex fluids. Lab on a Chip 17, 460-473 (2017).
  11. J. Dinic, Y. Zhang, L. N. Jimenez, V. Sharma, Extensional Relaxation Times of Dilute, Aqueous Polymer Solutions. ACS Macro Letters 4, 804-808 (2015).
  12. V. Sharma et al., The rheology of aqueous solutions of Ethyl Hydroxy-Ethyl Cellulose (EHEC) and its hydrophobically modified Analogue (hmEHEC): Extensional flow response in capillary break-up, jetting (ROJER) and in a cross-slot extensional rheometer. Soft Matter 11, 3251-3270 (2015).
  13. J. R. Castrejón-Pita et al., Plethora of transitions during breakup of liquid filaments. Proc. Natl. Acad. Sci. U.S.A. 112, 4582-4587 (2015).

FLOW-3D 제품소개

About FLOW-3D


HPC-enabled FLOW-3D v12.0

FLOW-3D 개발 회사

Flow Science Inc Logo Green.svg
IndustryComputational Fluid Dynamics Software
Founded1980
FounderDr. C.W. “Tony” Hirt
Headquarters
Santa Fe, New Mexico, USA
United States
Key people
Dr. Amir Isfahani, President & CEO
ProductsFLOW-3D, FLOW-3D CAST, FLOW-3D AM, FLOW-3D CLOUD, FlowSight
ServicesCFD consultation and services

FLOW-3D 개요

FLOW-3D는 미국 뉴멕시코주(New Mexico) 로스알라모스(Los Alamos)에 있는 Flow Scicence, Inc에서 개발한 범용 전산유체역학(Computational Fluid Dynamics) 프로그램입니다. 로스알라모스 국립연구소의 수치유체역학 연구실에서 F.Harlow, B. Nichols 및 T.Hirt 등에 의해 개발된 MAC(Marker and Cell) 방법과 SOLA-VOF 방식을 기초로 하여, Hirt 박사가 1980년에 Flow Science, Inc사를 설립하여 계속 프로그램을 발전시켰으며 1985년부터 FLOW-3D를 전세계에 배포하였습니다.

유체의 3차원 거동 해석을 수행하는데 사용되는 CFD모형은 몇몇 있으나, 유동해석에 적용할 물리모델 선정은 해석의 정밀도와 밀접한 관계가 있으므로, 해석하고자 하는 대상의 유동 특성을 분석하여 신중하게 결정하여야 합니다.

FLOW-3D는 자유표면(Free Surface) 해석에 있어서 매우 정확한 해석 결과를 제공합니다. 해석방법은 자유표면을 포함한 비정상 유동 상태를 기본으로 하며, 연속방정식, 3차원 운동량 보전방정식(Navier-Stokes eq.) 및 에너지 보존방정식 등을 적용할 수 있습니다.

FLOW-3D는 유한차분법을 사용하고 있으며, 유한요소법(FEM, Finite Element Method), 경계요소법(Boundary Element Method)등을 포함하여 자유표면을 포함하는 유동장 해석(Fluid Flow Analysis)에서 공기와 액체의 경계면을 정밀하게 표현 가능합니다.

유체의 난류 해석에 대해서는 혼합길이 모형, 난류 에너지 모형, RNG(Renormalized Group Theory)  k-ε 모형, k-ω 모형, LES 모형 등 6개 모형을 적용할 수 있으며, 자유표면 해석을 위하여 VOF(Volume of Fluid) 방정식을 사용하고, 격자 생성시 사용자가 가장 쉽게 만들 수 있는 직각형상격자는 형상을 더욱 정확하게 표현하기 위해 FAVOR(Fractional Area Volume Obstacle Representation) 기법을 각 방정식에 적용하고 있습니다.

FLOW-3D는 비압축성(Incompressible Fluid Flow), 압축성 유체(Compressible Fluid Flow)의 유동현상 뿐만 아니라 고체와의 열전달 현상을 해석할 수 있으며, 비정상 상태의 해석을 기본으로 합니다.

FLOW-3D v12.0은 모델 설정을 간소화하고 사용자 워크 플로우를 개선하는 GUI(그래픽 사용자 인터페이스)의 설계 및 기능에 있어 중요한 변화를 가져왔습니다. 최첨단 Immersed Boundary Method는 FLOW-3Dv12.0솔루션의 정확도를 높여 줍니다. 다른 특징적인 주요 개발에는 슬러지 안착 모델, 2-유체 2-온도 모델, 사용자가 자유 표면 흐름을 훨씬 더 빠르게 모델링 할 수 있는 Steady State Accelerator등이 있습니다.

물리 및 수치 모델

Immersed Boundary Method

힘과 에너지 손실에 대한 정확한 예측은 솔리드 바디 주변의 흐름과 관련된 많은 엔지니어링 문제를 모델링하는 데 중요합니다. FLOW-3D v12.0의 릴리스에는 이러한 문제 해결을 위해 설계된 새로운 고스트 셀 기반 Immersed Boundary Method (IBM)가 포함되어 있습니다. IBM은 내부 및 외부 흐름을 위해 벽 근처 해석을 위해 보다 정확한 솔루션을 제공하여 드래그 앤 리프트 힘의 계산을 개선합니다.

Two-field temperature for the two-fluid model

2유체 열 전달 모델은 각 유체에 대한 에너지 전달 공식을 분리하도록 확장되었습니다. 이제 각 유체에는 고유한 온도 변수가 있어 인터페이스 근처의 열 및 물질 전달 솔루션의 정확도를 향상시킵니다. 인터페이스에서의 열 전달은 시간의 표 함수가 될 수 있는 사용자 정의 열 전달 계수에 의해 제어됩니다.

슬러지 침전 모델 / Sludge settling model

중요 추가 기능인 새로운 슬러지 침전 모델은 도시 수처리 시설물 응용 분야에 사용하면 수처리 탱크 및 정화기의 고형 폐기물 역학을 모델링 할 수 있습니다. 침전 속도가 확산된 위상의 방울 크기에 대한 함수인 드리프트-플럭스 모델과 달리, 침전 속도는 슬러지 농도의 함수이며 기능적인 형태와 표 형태로 모두 입력 할 수 있습니다.

Steady-state accelerator for free surface flows

이름이 암시하듯이, 정상 상태 가속기는 안정된 상태의 솔루션에 대한 접근을 가속화합니다. 이는 작은 진폭의 중력과 모세관 현상을 감쇠하여 이루어지며 자유 표면 흐름에만 적용됩니다.

꾸준한 상태 가속기

Void particles

보이드 입자가 버블 및 위상 변경 모델에 추가되었습니다. 보이드 입자는 항력과 압력 힘을 통해 유체와 상호 작용하는 작은 기포의 역할을 하는 붕괴된 보이드 영역을 나타냅니다. 주변 유체 압력에 따라 크기가 변경되고 시뮬레이션이 끝난 후 최종 위치는 공기 침투 가능성을 나타냅니다.

Sediment scour model

침전물의 정확성과 안정성을 향상시키기 위해 침전물의 운반과 침식 모델을 정밀 조사하였다. 특히, 침전물 종에 대한 질량 보존이 크게 개선되었습니다.

Outflow pressure boundary condition

고정 압력 경계 조건에는 이제 압력 및 유체 비율을 제외한 모든 유량이 해당 경계의 상류에 있는 흐름 조건을 반영하는 ‘유출’ 옵션이 포함됩니다. 유출 압력 경계 조건은 고정 압력 및 연속성 경계 조건의 혼합입니다.

Moving particle sources

시뮬레이션 중에 입자 소스는 이동할 수 있습니다. 시간에 따른 변환 및 회전 속도는 표 형식으로 정의됩니다. 입자 소스의 운동은 소스에서 방출 된 입자의 초기 속도에 추가됩니다.

Variable center of gravity

중력 및 비 관성 기준 프레임 모델에서 시간 함수로서의 무게 중심의 위치는 외부 파일의 표로 정의할 수 있습니다. 이 기능은 연료를 소모하는 로켓을 모델링하고 단계를 분리할 때 유용합니다.

공기 유입 모델

가장 간단한 부피 기반 공기 유입 모델 옵션이 기존 질량 기반 모델로 대체되었습니다.  질량 기반 모델은 부피와 달리 주변 유체 압력에 따라 부피가 변화하는 동안 흡입된 공기량이 보존되기 때문에 물리학적 모델입니다.

Air entrainment model in FLOW-3D v12.0

Tracer diffusion / 트레이서 확산

유동 표면에서 생성된 추적 물질은 분자 및 난류 확산 과정에 의해 확산될 수 있으며, 예를 들어 실제 오염 물질의 거동을 모방합니다.

모델 설정

시뮬레이션 단위

이제 온도를 포함하여 단위계 시스템을 완전히 정의해야 합니다. 표준 단위 시스템이 제공됩니다. 또한 사용자는 선택한 옵션에서 질량, 시간 및 길이 단위를 정의하여 편리하며, 사용자 정의된 단위를 사용할 수 있습니다. 사용자는 또한 압력이 게이지 단위로 정의되는지 절대 단위로 정의되는지 여부를 지정해야 합니다. 기본 시뮬레이션 단위는 Preferences(기본 설정)에서 설정할 수 있습니다. 단위를 완벽하게 정의하면 FLOW-3D는 물리적 수량에 대한 기본 값을 정의하고 범용 상수를 설정할 수 있으므로 사용자가 필요로 하는 작업량을 최소화할 수 있습니다.

Shallow water model

얕은 물 모델에서 매닝의 거칠기

Manning의 거칠기 계수는 지형 표면의 전단 응력 평가를 위해 얕은 물 모델에서 구현되었습니다. 표면 결함의 크기를 기반으로 기존 거칠기 모델을 보완하며이 모델과 함께 사용할 수 있습니다. 표준 거칠기와 마찬가지로 매닝 계수는 구성 요소 또는 하위 구성 요소의 속성이거나 지형 래스터 데이터 세트에서 가져올 수 있습니다.

메시 생성

하단 및 상단 경계 좌표의 정의만으로 수직 방향의 메시 설정이 단순화되었습니다.

구성 요소 변환

사용자는 이제 여러 하위 구성 요소로 구성된 구성 요소에 회전, 변환 및 스케일링 변환을 적용하여 복잡한 형상 어셈블리 설정 프로세스를 단순화 할 수 있습니다. GMO (General Moving Object) 구성 요소의 경우, 이러한 변환을 구성 요소의 대칭 축과 정렬되도록 신체에 맞는 좌표계에 적용 할 수 있습니다.

런타임시 스레드 수 변경

시뮬레이션 중에 솔버가 사용하는 스레드 수를 변경하는 기능이 런타임 옵션 대화 상자에 추가되어 사용 가능한 스레드를 추가하거나 다른 태스크에 자원이 필요한 경우 스레드 수를 줄일 수 있습니다.

프로브 제어 열원

활성 시뮬레이션 제어가 형상 구성 요소와 관련된 heat sources로 확장되었습니다.  history probes로 열 방출을 제어 할 수 있습니다.

소스에서 시간에 따른 온도

질량 및 질량/모멘트 소스의 유체 온도는 이제 테이블 입력을 사용하여 시간의 함수로 정의 할 수 있습니다.

방사율 계수

공극으로의 복사 열 전달을위한 방사율 계수는 이제 사용자가 방사율과 스테판-볼츠만 상수를 지정하도록 요구하지 않고 직접 정의됩니다. 후자는 이제 단위 시스템을 기반으로 솔버에 의해 자동으로 설정됩니다.

Output

  • 등속 필드 솔버 옵션을 사용할 때 유량 속도를 선택한 데이터로 출력 할 수 있습니다.
  • 벽 접착력으로 인한 지오메트리 구성 요소의 토크는 기존 벽 접착력 출력과 함께 별도의 수량으로 일반 이력 데이터에 출력됩니다.
  • 난류 모델 출력이 요청 될 때 난류 에너지 및 소산과 함께 전단 속도 및 y +가 선택된 데이터로 자동 출력됩니다.
  • 공기 유입 모델 출력에 몇 가지 수량이 추가되었습니다. 자유 표면을 포함하는 모든 셀에서 혼입 된 공기 및 빠져 나가는 공기의 체적 플럭스가 재시작 및 선택된 데이터로 출력되어 사용자에게 공기가 혼입 및 탈선되는 위치 및 시간에 대한 자세한 정보를 제공합니다. 전체 계산 영역 및 각 샘플링 볼륨 에 대해이 두 수량의 시간 및 공간 통합 등가물이 일반 히스토리 로 출력됩니다.
  • 솔버의 출력 파일 flsgrf 의 최종 크기는 시뮬레이션이 끝날 때 보고됩니다.
  • 2 유체 시뮬레이션의 경우, 기존의 출력 수량 유체 체류 시간 및 유체 가 이동 한 거리는 이제 유체 # 1 및 # 2와 유체의 혼합물에 대해 별도로 계산됩니다.
  • 질량 입자의 경우, 각 종의 총 부피 및 질량이 계산되어 전체 계산 영역, 샘플링 볼륨 및 플럭스 표면에 대한 일반 히스토리 로 출력되어 입자 종 수에 대한 현재 출력을 보완합니다.
  • 최종 로컬 가스 압력 은 사용자가 가스 포획을 식별하고 연료 탱크의 배기 시스템 설계를 지원하는 데 도움이되는 선택적 출력량으로 추가되었습니다. 이 양은 유체로 채워지기 전에 셀의 마지막 공극 압력을 기록하며 단열 버블 모델과 함께 사용됩니다.

새로운 맞춤형 소스 루틴

새로운 사용자 정의 가능 소스 루틴이 추가되었으며 사용자의 개발 환경에서 액세스 할 수 있습니다.

소스 루틴 이름기술
cav_prod_calCavitation 생성과 소산 비율
sldg_uset슬러지 침전 속도
phchg_mass_flux증발 및 응축으로 인한 질량 플럭스
flhtccl유체 # 1과 # 2 사이의 열전달 계수
dsize_cal2 상 흐름에서 동적 액적 크기 모델의 응집 및 분해 속도
elstc_custom점탄성 유체에 대한 응력 방정식의 Source Terms

새로운 사용자 인터페이스

FLOW-3D 사용자 인터페이스는 완전히 새롭게 디자인되어 현대적이고 평평한 구조로 사용자의 작업 흐름을 획기적으로 간소화합니다.

Setup dock widgets

Physics, Fluids, Mesh 및 FAVOR ™를 포함한 모든 설정 작업이 지오 메트리 윈도우 주변에서 독 위젯으로 변환되어 모델 설정을 단일 탭으로 요약할 수 있습니다. 이러한 전환으로 인해 이전 버전의 복잡한 접이식 트리가 훨씬 깨끗하고 효율적인 메뉴 프레젠테이션으로 대체되어 사용자는 ModelSetup탭을 떠나지 않고도 모든 매개 변수에 쉽게 액세스 할 수 있습니다.

New Model Setup icons

새로운 모델 설정 디자인에는 설정 프로세스의 각 단계를 나타내는 새로운 아이콘이 있습니다.

Model setup icons - FLOW-3D v12.0

New Physics icons

RSS feed

새 RSS 피드부터 FLOW-3D v12.0의 시뮬레이션 관리자 탭이 개선되었습니다. FLOW-3D 를 시작하면 사용자에게 Flow Science의 최신 뉴스, 이벤트 및 블로그 게시물이 표시됩니다.

RSS feed - FLOW-3D

Configurable simulation monitor

시뮬레이션을 실행할 때 중요한 작업은 모니터링입니다. FLOW-3Dv1.0에서는 사용자가 시뮬레이션을 더 잘 모니터링할 수 있도록 SimulationManager의 플로팅 기능이 향상되었습니다. 사용자는 시뮬레이션 런타임 그래프를 통해 모니터링할 사용 가능한 모든 일반 기록 데이터 변수를 선택하고 각 그래프에 여러 변수를 추가할 수 있습니다. 이제 런타임에서 사용할 수 있는 일반 기록 데이터는 다음과 같습니다.

  • 최소/최대 유체 온도
  • 프로브 위치의 온도
  • 유동 표면 위치에서의 유량
  • 시뮬레이션 진단(예:시간 단계, 안정성 한계)
출입문에 유동 표면이 있는 대형 댐
Runtime plots of the flow rate at the gates of the large dam

Conforming 메쉬 시각화

사용자는 이제 새로운 FAVOR ™ 독 위젯을 통해 적합한 메쉬 블록을 시각화 할 수 있습니다.Visualize conforming mesh blocks

Large raster and STL data

데이터를 처리하는 데 걸리는 시간 때문에 큰 지오 메트리 데이터를 처리하는 것은 수고스러울 수 있습니다. 대형 지오 메트리 데이터를 처리하는 데는 여전히 상당한 시간이 걸릴 수 있지만, FLOW-3D는 이제 이러한 대규모 데이터 세트를 백그라운드 작업으로 로드하여 사용자가 데이터를 처리하는 동안 완전히 응답하고 중단 없는 인터페이스에서 작업을 계속할 수 있습니다

표면 장력 / Surface Tension

표면 장력 / Surface Tension

FLOW-3D에 추가 된 최초의 물리 모델 중 하나는 표면 장력이었습니다.

이 모델은 잉크젯, 무중력 환경에서의 액체 연료 거동 및 다양한 MEMS (마이크로 전자 기계 시스템) 장치와 같이 다양한 종류의 응용 분야에서 수년 동안 널리 사용되어 왔습니다. 이 후에 모델의 개선 및 확장에 대한 많은 사용자 요청이 처리되었습니다.
표면 장력에 대해 보다 나은 성능개선을 위해 FLOW-3D 버전 11에 대한 새로운 모델이 개발되었습니다. 이 모델은 계산된 모든 표면 장력의 정확성과 임의 형상의 솔리드 표면을 잡아 당기는 접착력의 정확성을 향상시킵니다. 또한 이 새로운 모델은 다공성 물질의 모세관 압력과 비 균일한 표면 장력으로 인한 접선 표면 장력을 가지고 있습니다.

새로운 모델의 예는 무중력에 포함된 원형 벽을 적시는 단순한 문제입니다.

그림 1은 실린더와 접촉각이 0 도인 물로 채워진 0.25m 직경의 실린더 75 %의 경우를 보여줍니다. 버블은 10 초 전에 벽에서 깨끗하게 분리되어 탱크를 가로 질러 움직입니다. 비 구형은 기포 표면에서 모세관 파가 전파되기 때문입니다.

그림 1. 0.0, 2.5, 5.0 및 10.0 초에 무중력에서 접촉 각이 0 인 실린더 표면의 유체 (적색) 습윤 표면.

다른 예가 그림5에 도시되어 있습니다. 2에서 서로 다른 밀도의 2 개의 초기 구형 방울이 (플롯의 색으로 표시됨) 단단한 벽을 향해 아래로 이동합니다. 플롯의 시간은 0.0, 0.01, 0.02 및 0.03 초입니다. 방울은 직경이 0.0017m, 밀도가 다르지만 표면 장력 계수는 1.872 뉴턴 / m입니다.

그림 2. 접시쪽으로 움직이는 구형의 물방울. 새로운 표면 장력 모델로 시뮬레이션. 색상은 밀도를 나타냅니다.

표면 장력 모델에 대해 자세히 알아보십시오.

Download the Flow Science Report on Surface Tension

Download Surface Tension Validation – Simple Test Problems