FLOW-3D HYDRO

응용 프로그램 개요

FLOW-3D HYDRO

FLOW-3D HYDRO 는 더 높은 수준의 정확도와 모델 해상도를 제공하기 위해 3D 비 유압 모델링 기능이 필요한 경우 고급 모델링 도구로 사용할 수 있습니다. 일반적인 모델링 응용 분야는 소형 댐 / 인프라, 운송 수력학, 복잡한 3D 하천 수력학, 열 부력 연기, 배수구 및 오염 물질 수송과 관련됩니다. 

FLOW-3D HYDRO의 핵심 기능은 전체 3D 모델과 동적으로 연결될 수있는 얕은 물 모델입니다. 

이 기능을 통해 사용자는 멀티 스케일 모델링 애플리케이션을위한 모델 도메인을 확장하여 필요한 모델 해상도로 계산 효율성을 극대화 할 수 있습니다. FLOW-3D HYDRO  또한 강 및 환경 응용 분야에 특화된 추가 기능과 고급 물리학을 포함합니다.

시뮬레이션 템플릿

FLOW-3D HYDRO 의 작업 공간 템플릿으로 시간을 절약하고 실수를 방지하며 일관된 모델을 실행하십시오 . 작업 공간 템플릿은 일반적인 응용 분야에 대한 유체 속성, 물리적 모델, 수치 설정 및 시뮬레이션 출력을 미리로드합니다.

작업 공간 템플릿은 7 가지 모델 클래스에 사용할 수 있습니다.

  • 자유 표면 – TruVOF (기본값)
  • 공기 유입
  • 열 기둥
  • 퇴적물 수송
  • 얕은 물
  • 자유 표면 – 2 유체 VOF
  • 자유 표면 없음

사전로드 된 예제 시뮬레이션

FLOW-3D HYDRO 의 40 개 이상의 사전로드 된 물 중심 예제 시뮬레이션 라이브러리는 애플리케이션 모델링을위한 훌륭한 시작점을 제공합니다. 사전로드 된 예제 시뮬레이션은 모델러에게 모델 설정 및 모범 사례의 로드맵뿐만 아니라 대부분의 애플리케이션에 대한 자세한 시작점을 제공합니다.이전다음

비디오 튜토리얼

비디오 자습서는 새로운 사용자가 다양한 응용 프로그램을 모델링하는 방법을 빠르게 배울 수있는 훌륭한 경로를 제공합니다. FLOW-3D HYDRO 비디오 튜토리얼 기능 :

  • 광범위한 응용 및 물리학을위한 AZ 단계별 기록
  • “사용 방법”정보
  • 모범 사례를위한 팁
  • CAD / GIS 데이터, 시뮬레이션 파일 및 후 처리 파일

고급 솔버 개발

Tailings Model

새로운 Tailings Model은 tailings dam failure로 인한 tailings runout을 시뮬레이션하기위한 고급 기능을 제공합니다. tailings정의에 대한 다층 접근 방식과 함께 미세하고 거친 입자 구성을 나타내는 이중 모드 점도 모델은 모든 방법으로 건설 된 tailings 댐의 모델링을 허용합니다. 

얕은 물, 3D 및 하이브리드 3D / 얕은 물 메싱을 포함한 유연한 메싱을 통해 얕은 지역에서 빠른 솔루션을 제공하면서 다층 tailings의 복잡성을 정확하게 모델링 할 수 있습니다. 점성 경계층의 정확한 표현을 위해 얕은 물 메시에 2 층 Herschel-Bulkley 점도 모델을 사용할 수 있습니다.

모델 하이라이트

  • 미세 입자 및 거친 입자 광미 조성물을위한 이중 모드 점도 모델
  • 침전, 패킹 및 입자 종의 난류 확산을 포함한 Tailings  수송
  • 얕은 물 메시를위한 2 층 Herschel-Bulkley 점도 모델
  • 3D, 얕은 물, 3D / 얕은 물 하이브리드 메시를 포함한 유연한 메시 접근 방식
  • Multi-layer, variable composition tailings for general definition of tailings dam construction

Shallow Water

FLOW-3D HYDRO 의 얕은 물 모델링 기능은 3D 메시를 얕은 물 메시와 결합하여 탁월한 모델링 다양성을 제공하는 고유 한 하이브리드 메시를 사용합니다. 압력 솔버의 수치 개선으로 더 안정적이고 빠른 시뮬레이션이 가능합니다. 하이브리드 메쉬의 하단 전단 응력 계산이 크게 향상되어 정확도가 더욱 향상되었습니다. 지형에 거칠기를 적용하는 새로운 방법에는 Strickler, Chezy, Nikuradse, Colebrook-White, Haaland 및 Ramette 방정식이 포함됩니다.

Two-Fluid VOF Model

sharp 인터페이스가 있거나 없는 압축 가능 또는 비압축성 2 유체 모델은 항상 1 유체 자유 표면 모델과 함께 FLOW-3D 에서 사용할 수 있습니다 . 사실, sharp 인터페이스 처리는 TruVOF 기술을 자유 표면 모델과 공유하며 상용 CFD 소프트웨어에서 고유합니다. 최근 개발에는 2- 필드 온도 및 인터페이스 슬립 모델이 포함되었습니다. 이 모델은 오일 / 물, 액체 / 증기, 물 / 공기 및 기타 2 상 시스템에 성공적으로 적용되었습니다.

FLOW-3D HYDRO 는 2- 유체 솔루션의 정확성과 안정성에서 두 가지 중요한 발전을보고 있습니다. 운동량과 질량 보존 방정식의 강화 된 결합은 특히 액체 / 기체 흐름에서 계면에서 운동량 보존을 향상시킵니다. 연속성 방정식에서 제한된 압축성 항의 확장 된 근사값은 더 빠르고 안정적인 2 유체 압력 솔버를 만듭니다.

예를 들어, 터널 및 드롭 샤프트 설계와 같은 유압 응용 분야에서 공기가 종종 중요한 역할을 하기 때문에 두 개발 모두 FLOW-3D HYDRO 릴리스에 적시에 적용됩니다. 일반적으로 낮은 마하 수로 인해 이러한 경우 물과 공기에 제한된 압축성이 사용됩니다.

고성능 컴퓨팅 및 클라우드

고성능 컴퓨팅 FLOW-3D HYDRO

일반 워크스테이션 또는 랩톱으로 많은 작업을 수행 할 수 있지만, 대형 시뮬레이션과 고화질 시뮬레이션은 더 많은 CPU 코어를 활용함으로써 엄청난 이점을 얻을 수 있습니다. FLOW-3D CLOUD 및 고성능 컴퓨팅은 더 빠르고 정확한 모델을 실행할 수있는 더 빠른 런타임과 더 많은 선택권을 제공합니다.

하천 및 환경 중심 애플리케이션

TRANSPORTATION HYDRAULICS
SMALL DAMS AND DIVERSIONS
RIVER HYDRAULICS
SEDIMENT TRANSPORT AND DEPOSITION
OUTFALLS EFFLUENTS
THERMAL PLUMES BUOYANT FLOWS

Case Studies

FLOW-3D Glossary

FLOW-3D 용어 사전 / 용어 설명

FLOW-3D 용어 사전 / 용어 설명

Drift Flux

드리프트 모델은 밀도가 서로 다른 두 혼합 유체 구성 요소의 상대적 흐름을 설명합니다. 구성 요소는 상이 다를 수도 있고, 상이 같지만(불가침) 유체가 다를 수도 있습니다. 분산된 위상 입자 크기가 클 경우 드리프트 모델의 적용성에 대한 제한이 존재할 수 있습니다. 이러한 제한은 일반적으로 메쉬 셀 크기의 10% 미만으로 분산된 위상 입자 크기를 유지함으로써 방지할 수 있습니다.

배플

얇은 기하학적 조각을 나타내는데 사용되는 2 차원 객체입니다. 이들은 전처리기에 의해 셀면으로 이동하고 유체의 흐름을 부분적으로 또는 완전히 차단하도록 작용합니다. 배플은 지정된 열 전달 계수를 지정할 수 있으며, 배플을 통과하는 유량(유속 표면)을 측정하는 데 사용할 수 있습니다.

경계 조건

도메인의 범위에서 솔루션을 정의합니다. 경계 위치에서 흐름의 실제 조건을 나타내는 경계 조건을 선택하는 것이 중요합니다.

CFD

CFD (Computational Fluid Dynamics)는 수치적 솔루션을 통해 컴퓨터의 유체 흐름을 시뮬레이션 하는 유체 역학의 학문적 분기입니다.

Complements

The inverse of a shape defines the complement. For example, the complement of a solid sphere is a spherical hole surrounded by solid material.

형상의 역은 complement를 정의합니다. 예를 들어, 솔리드 구의 complement은 솔리드 재료로 둘러싸인 구형 구멍입니다.

Client

클라이언트 컴퓨터는 FLOW-3D를 실행하지만 다른 컴퓨터 (서버 컴퓨터)에서 소프트웨어 라이센스를 획득하는 컴퓨터입니다.

Components

Components는 공간의 개체를 정의하며 하위 구성 요소로 구성됩니다. 구성 요소는 열 전도성, 비열 및 표면 거칠기와 같은 재료 특성을 가질 수 있습니다.

Custom result

시뮬레이션 중 또는 완료 후 사용자가 생성한 데이터를 그래픽으로 표시합니다. 생성하려면 사용자가 flsgrf결과 파일을 연 다음 플로팅 매개 변수(예 : 플로팅 할 도메인 부분, 플로팅 할 수량 등)를 선택해야합니다.

Domain

지배 방정식이 solved되는 영역. 이것은 메쉬의 범위에 의해 정의됩니다.

Diagnostics

전 처리기 및 솔버의 진행 상황과 오류 및 경고에 대한 정보가 포함된 파일 세트입니다.

EPSI

압력/연속 반복이 어느 지점에서 수렴되는지를 결정하는데 사용된 수렴 기준입니다. 기본 숫자 설정을 사용하면 이 값은 FLOW-3D에 의해 자동으로 계산되며 시간 단계가 증가함에 따라 작아집니다.

Existing result

prpplt.* 또는 flsplt.* 파일은 전처리 종료 솔버 실행 종료시 또는 자동으로 생성되는 플롯 파일입니다.

F3D_HOME

FLOW-3D 프로그램 파일이 있는 디렉토리를 정의하는 환경 변수.

Floating license

FLOW-3D는 서버 시스템에 라이센스를 액세스하는 각 클라이언트 컴퓨터와 컴퓨터 네트워크에서 실행합니다. 허용하는 라이센스 최대 동시 시뮬레이션 수는 구매한 솔버 토큰 수에 의해 제한됩니다.

Flsgrf file

솔버가 생성한 결과 파일. 이 파일은 사전 정의된 시간 간격으로 생성된 정보를 포함하며 그래픽 디스플레이를 생성하는 데 사용됩니다. 사용자 지정 플로팅 중에 포스트 프로세서에서 사용합니다.

Flsplt file

솔버가 자동으로 생성한 플롯 파일입니다. 이 파일에는 시뮬레이션의 히스토리 데이터, 메시 등에 대한 기본 정보와의 $GRAFIC 이름 목록에 사전 정의된 그래픽 요청이 포함되어 prepin.* 파일 안에 있습니다.

Fluid #1 surface area

선택한 길이 단위의 자유 표면 영역이 제곱 됩니다. 인터페이스가 예리한 문제에만 해당됩니다.

Fluid thermal energy

영역에 존재하는 모든 유체에 포함된 총 열 에너지 (에너지 전송이 켜져 있는 시뮬레이션에만 해당).

Free surface

유체와 유체 사이의 인터페이스. FLOW-3D에서 이 인터페이스는 전단이 없는 것으로 가정되며, 이는 빈 공간에 있는 가스가 유체에 무시할 수 있는 트랙션을 발휘함을 의미한다.

GUI

” Graphical User Interface”.  GUI는 사용자가 FLOW-3D를 제어할 수 있는 그래픽 패널, 대화 상자 및 창을 제공합니다.

Iteration count

각 시간 단계에서 필요한 압력/연속 반복 횟수입니다. 유체량을 유지하고 유체 전체의 정확한 압력을 계산하려면 압력/연속 반복이 필요합니다.

License file

사용자가 FLOW-3D 를 실행할 수 있도록 암호화된 정보가 포함된 Flow Science에서 제공 한 전자 파일 입니다.

License server

플로팅 라이센스 시스템의 작동을 활성화하기 위해 FLEXlm 라이센스 소프트웨어가 설치된 시스템. FLOW-3D는 License Server에 설치할 필요가 없습니다.

Licensing

FLOW-3D 실행을 제어하는 ​​FLEXlm 소프트웨어.

Max. residual

압력/연속 반복의 최종 반복에 대한 연속성 방정식의 실제 차이. 이 값은 일반적으로 xxxx 반복에서 압력 반복이 수렴되지 않는 한 epsi보다 작다.

Mean kinetic energy

도메인에 존재하는 유체의 총 질량으로 나눈 모든 계산 셀의 운동 에너지의 합계. 시간이 지남에 따라 이 양이 변동이 멈추면 정상 상태에 도달했다는 좋은 지표다.

Node-locked license

특정 컴퓨터에 고정된 라이센스. 노드 잠금 라이센스는 네트워크를 통해 액세스 할 수 없으므로 일반적으로 모든 작업을 한 컴퓨터에서 수행해야하는 경우에만 사용됩니다.

Non-inertial reference frame

가속하는 기준 프레임. 비 관성 기준 프레임은 이동 컨테이너를 모방하는데 사용될 수 있습니다.

Pltfsi

1D 및 2D 플롯을 생성하는 FLOW-3D에 포함된 그래픽 디스플레이 프로그램.

Postprocessor

내 프로그램 FLOW-3D 읽을 수 있는 생산 데이터 파일에 해결사 출력 데이터 처리 FLOW-3D 의 또는 타사의 시각화 프로그램 또는 생산 텍스트 데이터는 타사 소프트웨어 프로그램에서 읽을 수 있습니다.

Prepin file

FLOW-3D 시뮬레이션을 생성하는 데 필요한 모든 정보가 포함된 텍스트 파일. GUI를 사용하거나 텍스트 편집기를 사용하여 수동으로 작성할 수 있습니다.

Preprocessor

솔버의 실행을 준비하기 위해 입력 파일을 기반으로 메쉬 및 초기 조건을 생성하는 FLOW-3D 내의 프로그램.

Prpgrf file

전처리기에 의해 생성된 결과 파일. 전 처리기의 정보를 포함하며 후 처리기에서 사용자 플롯을 생성하는 데 사용할 수 있습니다. 이 파일은 미리보기 버튼을 선택하거나 시뮬레이션에서 사전 프로세서(runpre 사용)를 실행하는 경우에만 실행됩니다.

Prpplt file

전처리기에 의해 자동으로 생성된 파일을 플롯 합니다. 메시, 구성 요소, 초기 조건 및 재료 특성에 대한 정보가 포함되어 있습니다.

Restart simulation

이전 시뮬레이션에서 계속 이어지는 시뮬레이션. 이전 시뮬레이션의 결과는 재시작을 위한 초기 조건 및 (선택적으로) 경계 조건을 생성하는 데 사용된다.

Server

라이센스 서버를 호스팅하는 시스템

Stability limit

각 시간 단계에서 사용할 수 있는 최대 시간 단계. 더 큰 시간 단계는 수치적 불안정성과 비물리적 결과로 이어질 것이다.

STL (Stereolithography) File

.STL 파일 형식은 일련의 삼각형이 있는 솔리드 모델의 표면에 근접한 표준 데이터 전송 형식이다. 삼각형은 가장자리에서 결합해야 하며 일관된 방향을 가리키는 정규식이 있어야 한다.

Solid fraction

응고된 영역의 유체 분율 (응고 모델이 켜져 있는 시뮬레이션에만 해당).

Solver

입력 파일에 정의된 흐름 문제를 시뮬레이션 하는 방정식 시스템을 해결하는 FLOW-3D 내의 프로그램입니다.

STL Viewer

스테레오리소그래피(STL) 파일을 표시하는 특수 유틸리티입니다. STL 파일은 CAD 소프트웨어로 제작되며 3 차원 객체의 표면을 형성하는 많은 삼각형으로 구성됩니다. 의 STL 뷰어 FLOW-3D는 메인 메뉴에서 유틸리티/STL 뷰어를 클릭하여 GUI를 통해 액세스 할 수 있습니다. 그러면 뷰어가 별도의 창에서 열립니다. 메쉬 및 형상 탭에서 STL 파일을 열고 볼 수도 있습니다.

Subcomponents

하위 구성 요소는 구성 요소라고하는 더 큰 모양을 형성하기 위해 결합할 수 있는 기하학적 모양입니다. 하위 구성 요소는 재료를 추가하거나 (고체로) 다른 하위 구성 요소에서 재료를 제거하거나 (구멍으로) 또는 모양 외부에 재료를 추가하도록 정의할 수 있습니다.

Time-step size

계산에 사용된 실제 시간 단계. 이 값은 안정성 한계와 같거나 작을 수 있습니다.

단위

단위는 물리적 특성에 설정된 값을 기반으로 합니다. 메쉬 블록 범위 및 셀 길이와 같은 항목은 이러한 물리적 속성을 설정하는 데 사용되는 단위를 자동으로 따릅니다.

Volume error (%)

유체 부피의 백분율은 주어진 시간에 도메인에 존재하는 총 유체의 백분율로 설명되지 않습니다. 따라서, 존재하는 총 부피가 작기 때문에 유체가 시스템에서 배출되는 시뮬레이션의 경우 부피 백분율 오차가 발생할 수 있습니다.

Volume of fluid #1

선택한 길이 단위로 입방체에 존재하는 유체 # 1의 총 부피입니다. 2 유체 문제의 경우, 유체 # 2의 부피는 항상 도메인 부피에서 유체 # 1의 부피를 뺀 값입니다.

Wall shear stress

FLOW-3D의 사용자가 또는 벽과 객체 인터페이스에서의 전단 응력의 계산 끌 수 있습니다 옵션을 선택합니다. “no-slip” 인터페이스의 효과를 모델링하려면 벽 전단 응력을 설정해야 합니다.

Workspace

작업 공간은 시뮬레이션 프로젝트를 위한 파일 컨테이너입니다. 작업 공간은 사용자가 FLOW-3D 뿐만 아니라 하드 드라이브에서도 작업을 구성하는 데 도움아 됩니다.F

FLOW-3D 및TruVOF는 미국 및 기타 국가에서 등록 상표입니다.

레이저 용접 수치해석 (FLOW-3D WELD)

FLOW-3D WELD Products

레이저 용접 수치해석 (FLOW-3D WELD)

FLOW-3D@ WELD는 레이저 용접 공정에 대한 정확한 시뮬레이션 기능을 제공하여 최적화된 공정을 개발하게 합니다. 더 나은 공정 제어를 통해 기공, 열 영향 영역을 최소화하고 미세 구조 변화를 제어할 수 있습니다.

레이저 용접 프로세스를 정확하게 시뮬레이션하기 위해 FLOW-3D@ WELD는 레이저 열원, 레이저-재료 상호 작용, 유체 흐름, 열 전달, 표면 장력, 응고, 다중 레이저 반사 및 위상 변화와 같은 모든 관련 물리 모델을 제공합니다.

Laser Welding

최근에는 뛰어난 생산성과 속도, 낮은 열 입력이 결합되어 기존의 용접 프로세스를 대체하는 레이저 용접 프로세스가 주목 받고 있습니다. 레이저 용접이 제공하는 장점은 용접강도가 좋고, 열 영향 부위가 작으며, 정밀도가 낮고 변형이 적으며, 강철, 알루미늄, 티타늄 및 이종 금속을 포함한 광범위한 금속 및 합금을 용접 할 수 있는 기능이 있습니다.

FLOW-3D@는 레이저 용접 공정에 대한 강력한 통찰력을 제공하고 궁극적으로 프로세스 최적화를 달성하는 데 도움이 됩니다.

보다 나은 프로세스 제어를 통해 기공을 최소화할 수 있습니다. 열 영향부위 및 미세조직을 제어가 가능합니다. FLOW-3D는 자유표면 추적 알고리즘을 통해 매우 복잡한 용접 POOL 시뮬레이션을 해석하는데 매우 적합합니다.

용접 모듈은 레이저 소스에 의해 생성된 Heat flux, 용융 금속에 대한 증발압력, shield gas 효과, 용융 풀의 반동압력 및 다중 레이저 반사와 같은 물리적 모델을 FLOW-3D에 적용하기 위해 개발되었습니다. 키홀 용접과 같은 현실적인 프로세스 시뮬레이션을 위해서는 모든 관련 물리적 현상을 적용하는 것이 중요합니다.

FLOW-3D는 레이저 용접의 conduction and keyhole 방식을 시뮬레이션 할 수 있습니다. 전 세계의 연구원들은 FLOW-3D를 사용하여 용접역학을 분석하고, 공정 매개 변수를 최적화하여 기공을 최소화하며, 레이저 용접공정에서의 dendrite 결정 성장 양상을 예측합니다.

Shallow penetration weld (top left); deep penetration weld with shield gas effects (top right); deep penetration weld with shield gas and evaporation pressure (bottom left); and deep penetration weld with shield gas, evaporation pressure and multiple laser reflections effects (bottom right).

Full Penetration Laser Welding Experiments

한국 카이스트와 독일 BAM은 16K kW레이저를 사용하여 10mm강판에 완전 침투 레이저 용접 실험을 수행하였습니다. CCD카메라의 도움을 받아 완전 용입 레이저 용접으로 형성된 상단 및 하단 용융풀 거동을 확인할 수 있었습니다. 그들은 또한 FLOW-3D 로 용접 공정 해석으로 해석과 실험결과의 경향이 일치하는 것을 알 수 있었습니다.

Experimental setup with CCD cameras observing the top and bottom molten pools
Schematic of computation domain in FLOW-3D

 

Simulation results at the top show melt pool lengths of 8mm and 15mm, whereas experiments indicated melt pool lengths of 7mm and 13mm

Laser Welding Porosity Case Study

General Motors, Michigan, 중국의 상하이 대학교는 용접 공정 변수, 즉 keyhole 용접에서 기공의 발생에 대해 용접 속도 및 용접 각도와 같은 공정 매개 변수가 미치는 영향을 알아보기 위해 협력하여 연구를 진행하였습니다.

레이저 용접된 Al 접합부 단면의 기공을 분석합니다. Keyhole이 유도 된 기공들은 유동 역학으로 인해 발생되고 균열을 일으킬 수 있습니다. 최적화 공정의 매개변수는 이러한 종류의 기공을 완화할 수 있습니다. FLOW-3D를 사용하여 연구원들은 증발 및 반동 압력, 용융풀, 온도에 따른 표면장력 및 Keyhole내의 다중 레이저 반사, 프레넬 흡수를 포함한 모든 중요한 물리적 현상을 설명했습니다.

연구진은 시뮬레이션 모델을 기반으로 Keyhole 용접에서 생성된 기공들의 주요 원인으로 불안정한 Keyhole을 규정하였습니다. 아래 이미지에서 볼 수 있듯이 뒤쪽 용융 풀의 과도한 재순환은 뒤쪽 용융 풀이 앞쪽 용융 풀 경계를 무너뜨리며 기공들을 생성시킵니다. 갇힌 공간이 증가하는 응고 전면에 의해 갇혔을때 기공들이 발생되었습니다.

Distribution of porosity in longitudinal welding sections as seen in simulations (top) and experiments (bottom)

용접 속도가 빠를수록 더 큰 keyhole이 생성되며 이로 인해, 보다 안정적인 keyhole이 생성됩니다. 연구진은 FLOW-3D를 사용하여 용접 속도와 용접 경사각으로 기공들의 생성을 완화시킬 수 있었습니다.

컴팩트 디스크 ELISA 칩 [2]

컴팩트 디스크 미세 유체 장치: Optimizing Real Estate

Compact Disc Microfluidic Devices: Optimizing Real Estate

미세 유체 장치 사용자의 증가하는 기대를 충족하려면 작은 미세 유체 장치에서 제한된 공간을 최적화하는 것이 중요합니다. 사용자는 단일 미세 유체 장치에서 최대의 기능과 여러 병렬 작업을 기대합니다. 제한된 공간을 최적화하는 문제는 이러한 장치의 많은 물리적 이점에도 불구하고 회전하는 미세 유체 장치로 확장됩니다. 회전 에너지를 이용하여 미세 유체 작업을 수행하는 회전 장치를 컴팩트 디스크 (CD) 미세 유체 장치라고합니다.

컴팩트 디스크 ELISA 칩 [1]
컴팩트 디스크 ELISA 칩 [2]
컴팩트 디스크 ELISA 칩 [2]

10 년 넘게 CD는 혈액 진단을위한 신속한 면역 분석 및 임상 생화학에서 지속적으로 장점을 보여 왔습니다. 마이크로 토탈 분석 시스템 (μTAS)으로 사용되며, 여러 개별 분석이 내장되어 단일 칩에서 동시에 실행됩니다. 핸즈프리 제어를 위해 프로그래밍 된 간단하고 저렴한 모터에서 작동하며 자석이나 표면 처리와 같은 외부 액추에이터가 필요하지 않습니다. 기본적으로 CD는 훌륭합니다! 그러나 공짜 점심 같은 것은 없습니다. 단방향 (방사형) 원심력으로 인해 CD는 회전하지 않는 미세 유체 장치보다 빠르게 공간이 부족합니다. 유체는 방사형으로 바깥쪽으로 만 이동하므로 CD가 수행 할 수있는 분석 단계의 수가 제한됩니다.

그림 3. CD 채널 내부의 방사형 물 기둥에 적용되는 다양한 신체 힘을 강조하는 회로도. 방사상으로 바깥쪽으로 작용하는 원심력을 확인합니다.
그림 3. CD 채널 내부의 방사형 물 기둥에 적용되는 다양한 신체 힘을 강조하는 회로도. 방사상으로 바깥쪽으로 작용하는 원심력을 확인합니다.

CD의 단 방향성 극복

Gorkin    [3]에서는 CD의 단 방향성 제약을 극복하기 위해 공압 펌핑이 제안되었습니다. 아이디어는 원심 에너지를 압축 에너지로 저장하고 다시 풀어서 유체를 중심으로 발사하는 것입니다. 아래 이미지는 로딩 챔버, 흡입 하위 구획 및 압축 하위 구획의 세 개의 챔버가있는 비교적 간단한 미세 유체 칩을 보여줍니다.

그림 4. CD 사진
그림 4. CD 사진
그림 5. FLOW-3D에서 모방 된 CD 디자인
그림 5. FLOW-3D에서 모방 된 CD 디자인

공압 펌핑 프로세스

유체가 로딩 챔버로 들어간 다음, 흡입 하위 구획을 통해 공기가 갇힌 압축 하위 구획으로 이동합니다. 공기가 갇 히면 CD가 특정 각속도로 회전하여 갇힌 공기가 압축됩니다. 공기가 더 이상 압축 할 수없는 경우 (안정 상태에 도달했기 때문에), 회전 속도가 감소하거나 완전히 꺼져 (누군가이 작업을 수행하고 있습니까? 아니면 장치가 수행하고 있습니까?) 유체가 로딩 챔버로 다시 펌핑됩니다. 이 마지막 단계는 이완 단계입니다. 공압 펌핑 공정의 5 단계는 다음과 같습니다.

그림 6. CD의 5 단계 공압 펌핑 [3]
그림 6. CD의 5 단계 공압 펌핑 [3]

회전 속도의 영향

회전 속도가 다르면 압축 하위 구획에서 공기의 압축 수준이 다릅니다. 회전 속도가 높을수록 유체가 공기에 더 세게 밀려 공기가 더 많이 압축됩니다. 그러나 공기가 압축 될 수있는 양에는 한계가 있습니다. 사실, 공기의 압축은 특정 회전 속도 이상으로 점진적으로 증가합니다. 압축 하위 구획의 부피는 회전 속도가 증가함에 따라 감소합니다. 흡입구의 액체 위치는 디스크 중앙에서 흡입 하위 구획의 유체 수준까지의 거리입니다. 이 거리는 증가합니다. 즉, 회전 속도가 증가함에 따라 유체가 중심에서 멀어집니다.

그림 7. 회전 주파수가 증가하면 압축이 증가합니다. [3]
그림 7. 회전 주파수가 증가하면 압축이 증가합니다. [3]

CD 미세 유체 장치 모델링

실험은 미세 유체 장치 설계의 핵심입니다. 그러나 충분한 실험을 수행하고 각 실험에 대한 완벽한 제어 환경을 유지하는 것은 불가능할 수 있습니다. 복잡한 설계에는 복잡한 실험 설정 및 분석이 필요합니다. FLOW-3D 의 정확하고 포괄적 인 다중 물리  모델링 기능 은 미세 유체 설계에 대한 통찰력과이를 최적화하는 방법을 제공합니다. FLOW-3D가  위에서 논의한 CD 미세 유체 장치에 대한 실험적 및 이론적 결과와 어떻게 비교되는지 보여 드리겠습니다  .

CD 미세 유체 장치에 대한 실험적 및 이론적 결과와 비교
CD 미세 유체 장치에 대한 실험적 및 이론적 결과와 비교

이미지 시퀀스는 실험 및 FLOW-3D  시뮬레이션 결과 의 시각적 비교를 제공합니다  . 두 유체 (공기 및 물) 압축 가능 모델을 사용하여 서로 다른 회전 속도에 대해 챔버 내부의 유체 역학을 시뮬레이션했습니다. 회귀 분석을 사용하여 아래 플롯에서 이러한 시각적 비교를 정량화하면 FLOW-3D  와 실험 결과,  FLOW-3D  및 분석 결과 간에 탁월한 상관 관계 (R 2 > 0.99)가 제공  됩니다.

그림 9. FLOW-3D 데이터와 실험 데이터의 비교. (Poly는 다항식 곡선 맞춤을 의미합니다.)
그림 9. FLOW-3D 데이터와 실험 데이터의 비교. (Poly는 다항식 곡선 맞춤을 의미합니다.)

시뮬레이션은 또한 다양한 회전 속도에 대한 정상 상태에 대한 접근 방식을 보여줍니다. 아래의 애니메이션은 CD의 운동 에너지 변동을 1000rpm nd 7000rpm에서 보여줍니다. 더 빠른 속도는 더 빠른 정상 상태를 강제하지만 정상 상태에 도달할 때까지 수위를 빠르게 변동시킵니다. 저속 시뮬레이션의 경우 그 반대입니다.

Mean kinetic energy fluctuations for the CD rotating at 1000 rpm
Mean kinetic energy fluctuations for the CD rotating at 1000 rpm
Mean kinetic energy fluctuations for the CD rotating at 7000 rpm
Mean kinetic energy fluctuations for the CD rotating at 7000 rpm

전반적으로  FLOW-3D  는 실험 결과를 정확하게 검증합니다. 사소한 오류는 부정확 한 지오메트리 (CAD) 생성 및 / 또는 물과 공기 사이의 인터페이스를 엄격하게 정의하기 때문일 수 있습니다. 이 사례 연구는 FLOW-3D  가 실험 결과를 검증하고 컴팩트 디스크 설계의 신뢰도를 높이는 데 효과적으로 사용될 수 있음을 보여줍니다  .

References

[1] He, Hongyan et al. “Design and Testing of a Microfluidic Biochip for Cytokine Enzyme-Linked Immunosorbent Assay”. Biomicrofluidics 3(2):22401 February 2009

[2] Roy, Emmanuel, et al. “From Cellular Lysis to Microarray Detection, an Integrated Thermoplastic Elastomer (TPE) Point of Care Lab on a Disc.” Lab on a Chip, vol. 15, no. 2

[3] Gorkin III, Robert et al. “Pneumatic pumping in centrifugal microfluidic platforms”. February 2010 Springerlink.com

벨기에 Zele에서 나온 WWTP의 개략도

활성화 된 슬러지 모델링

Activated Sludge Model

폐수 처리 플랜트 (WWTP) 내부의 생화학 적 반응 및 유체 역학에 대한 자세한 이해는 설계자와 엔지니어가 새로운 플랜트 설계를 평가하고, 관리 결정을 정량화하고, 새로운 제어 계획을 개발하고, 안전한 작업자 교육을 제공하는 데 도움이 될 수 있습니다. 이 블로그에서는 독자들에게 대규모 생화학 반응 시스템을 동적으로 해결 하는 FLOW-3D 의 새로운 ASM (Activated Sludge Model)을 소개합니다.

폭기조

폭기조는 대부분의 생화학 반응이 WWTP의 2 차 처리 부분에서 발생하는 곳입니다. 일반적으로 폭기 탱크는 대부분의 생화학 반응이 완료되는 데 걸리는 시간을 허용하는 긴 경로를 가지고 있습니다. 종이 폭기조의 전체 길이를 횡단하는 데 걸리는 시간을 잔류 시간이라고합니다. 폭기조에 산소가 주입되어 폐수에서 박테리아가 증식합니다. 박테리아는 산소를 사용하여 물에있는 폐기물을 분해하고 그렇게하면서 플록 또는 슬러지 블랭킷이라고하는 응집체를 형성합니다. 활성화 된 슬러지의 일부는 폐수의 생화학 적 처리를 더욱 촉진하기 위해 폭기조로 다시 재활용됩니다.

벨기에 Zele에서 나온 WWTP의 개략도
벨기에 Zele에서 나온 WWTP의 개략도

생화학 반응의 표준 시스템

국제 물 협회 (IWA)는 지난 40 년간 생화학 적 반응을 설명하는 세 가지 주요 수학적 시스템을 제안했다. 이러한 각 시스템 인 ASM-1, ASM-2 및 ASM-3은 폭기조 내부의 다양한 종의 성장 및 붕괴 역학을 다양한 세부 수준으로 포착합니다. ASM-3이 가장 포괄적입니다. 첫 번째 시스템 인 ASM-1은 아래 표 형식과 확장 형식으로 표시됩니다.

결합 편미분 방정식의 확장 시스템으로서의 생화학 반응의 ASM-1 시스템
결합 편미분 방정식의 확장 시스템으로서의 생화학 반응의 ASM-1 시스템

ASM 솔버 기능

대부분의 생화학 반응은 Monod 모델 또는 유사한 모델을 기반으로합니다. Monod 모델은 미생물의 성장 및 붕괴 속도를 예측하는 수학적 모델이며 간단한 방정식으로 설명됩니다.

여기서 a 와 k 는 최대 비 성장률 상수이고 기질 농도는 최대 비 성장률의 절반에 해당합니다. C 는 시간에 따라 변화하는 미생물 종의 농도 t 입니다. Monod 모델은 종의 농도에 따라 반응의 순서를 동적으로 변경하는 특성이 있습니다.

For C   >> A는 , 변화율 C는  0 차에 접근한다.

For C   << a는 , 변화율 C는 일차 접근한다.

이 모든 것은 미생물 종의 농도가 높으면 썩고 자라는 속도가 빨라지고, 종의 양이 적으면 썩거나 자라는 속도가 느리다는 것입니다. Monod 방정식의 해는 다음과 같이 Lambert 함수에 의해 제공됩니다.

간단한 Monod 방정식에 대한 분석 솔루션과 FLOW-3D 솔루션의 비교
간단한 Monod 방정식에 대한 분석 솔루션과 FLOW-3D 솔루션의 비교

생화학 반응을 설명하는 표준 시스템에는 Monod 용어의 긴 사슬이 포함되어 있습니다. FLOW-3D 의 ASM 모델은 WWTP에서 박테리아 종의 Monod 기반 성장 및 붕괴를 완벽하게 추적 할 수 있습니다. ASM 모델은 FLOW-3D 의 유체 역학 솔버 와 통합되어 속도 및 압력 장을 기반으로 한 박테리아의 움직임이 성장 및 붕괴 속도와 결합 될 수 있습니다.

FLOW-3D 의 ASM 솔버 결과가 벨기에 Zele의 폐수 처리장 (WWTP)에서 배출 될 때 다양한 유입수 종 농도의 붕괴 및 성장에 대해 보여줄 것 입니다. 종 및 유체 역학 계산을 정확하게 추적하면 폐수 처리 전문가가 정량적으로 뒷받침되는 설계 및 운영 결정을 내릴 수 있습니다.

Zele WWTP

Zele WWTP는 1983 년 50,000 명의 주민을 위해 벨기에에서 건설되었습니다. 일반적으로이 WWTP의 유입수는 가정용 폐수 40 %와 산업 폐수 60 %로 구성됩니다. 1 차 처리 공정 후 유입수는 생물학적 활성 슬러지 처리장으로 흘러 재활용 활성 슬러지와 혼합됩니다.

벨기에 Zele에서 나온 WWTP의 개략도 [2]. 녹색 상자는 2 차 처리 과정을 나타냅니다.
벨기에 Zele에서 나온 WWTP의 개략도 [2]. 녹색 상자는 2 차 처리 과정을 나타냅니다.

활성 오니 조 또는 폭기조는 약 400 m의 레인 6으로 분할되어 하나의 플러그 유동 폭기조 구성 3 각. 폭기조에서 나오는 유출 물은 각각 2050 m 3 용적의 2 개의 2 차 정화기 (SC1 및 SC2)로 이동합니다 . 최종 폐수는 인근 하천으로 배출됩니다. 2 차 정화기 아래에서 활성 슬러지의 일부는 폭기조로 다시 재활용되어 2 차 처리의 효율성을 높입니다.

우리는 2 차 처리 구성 요소의 기하학적 구조와 다양한 종의 유입 농도에 대한 자세한 정보를 이용할 수 있기 때문에 사례 연구를 위해이 WWTP를 선택했습니다. 정보는 상세하지만 완전하지는 않으며이 불완전한 정보는 폐수 농도에 중대한 영향을 미칠 것이며 나중에 논의 할 것입니다.

기하학, 메싱 및 물리학

지오메트리 생성 및 메싱은 간단했습니다. FLOW-3D 에는 완전한 WWTP를 완전히 정의하는 데 사용 된 기본 지오메트리 모양 모음이 있습니다. 이러한 모양은 생성하기 쉽고 외부 CAD 소프트웨어를 사용하여 생성 된 일부 지오메트리와 달리 오류가 없습니다. 마찬가지로, 구조화 된 그리드를 사용하면 구조화되지 않은 그리드 생성과 관련된 일반적인 오류를 처리하는 시간이 절약되었습니다.

폭기조 내부의 물리학은 복잡하며 질량 및 운동량 보존 방정식 (Navier-Stokes 방정식), 종 수송, 반응 역학, 산소 용해 및 연속 밀도 평가의 완전한 시스템을 해결해야합니다. FLOW-3D 는 가장 정확한 계산을 위해 완전히 결합 된 방식으로 이러한 모든 물리학을 설명합니다.

FLOW-3D의 Zele WWTP 설정. 화살표는 흐름 방향을 나타내며 유입수는 녹색 화살표의 시작 부분에서 도메인으로 들어갑니다.
FLOW-3D의 Zele WWTP 설정. 화살표는 흐름 방향을 나타내며 유입수는 녹색 화살표의 시작 부분에서 도메인으로 들어갑니다.

세 가지 표준 수학적 모델 인 ASM-1, ASM-2 및 ASM-3 중에서 연구자들은이 WWTP에서 ASM-1 수학적 모델을 사용합니다. 이는 간단하면서도 많은 중요한 생화학 과정을 다루기 때문입니다. ASM-1 모델은 일반적으로 폐수에서 발견되거나 처리 과정에서 생성되는 13 종의 진화를 고려합니다 [표 1].

종 IDZele의 초기 유입 농도 (mg / l)
가용성 불활성 유기물SI7.5
쉽게 생분해되는 기질SS400.0
미립자 불활성 유기물XI40.0
천천히 생분해되는 기질XS40.0
활성 종속 영양 바이오 매스XB, H120.0
활성 독립 영양 바이오 매스XB, A5.0
바이오 매스 붕괴로 인한 미립자 제품XP0.0
산소SO0.0
질산염 및 아질산염 질소SNO0.0
암모늄 질소SNH15.0
용해성 생분해 성 유기 질소SND8.2
미립자 생분해 성 유기 질소XND11.3
알칼리도SALKNot included

표 1. 표준 ASM-1 수학 시스템의 종 목록과 Zele WWTP에서 측정 된 초기 유입수 농도. 이러한 초기 농도 중 일부는 추론되며 큰 불확실성이 관련 될 수 있습니다. S와 X는 각각 용해성 물질과 미립자 물질을 나타냅니다.

이들 종 각각은 반응하지 않는 불활성 종 (SI 및 XI)을 제외하고 하나 이상의 생화학 적 과정에 의존합니다. 불활성 종의 유입 및 유출 농도는 XI의 경우와 같이 침전으로 인해 달라질 수 있습니다. SALK는 WWTP에서 측정되지 않았기 때문에이 사례 연구에서 무시되었습니다.

관심 유출량

폐수 엔지니어가 관심을 갖는 주요 유출량은 총 화학적 산소 요구량 (COD tot ), 암모늄 질소 (SNH) 농도, 아질산염 및 질산염 질소 (SNO) 및 총 킬달 질소 (TKN)입니다.

  • COD tot = SI + SS + XI + XS
  • TKN ~ XND + SND + SNH

이 양은 처리 된 물의 전반적인 품질을 나타냅니다.

유출량측정 된 유입 농도 (mg / l)FLOW-3D 유출 농도 (mg / l)
CODtot600264.04
SNH1530.34
SNO01.86
TKN3537.28

총 COD, SNH 및 TKN의 농도는 폐수가 폭기조를 통과하여 WWTP를 빠져 나 가면서 감소해야합니다. 이 동작은 총 COD [표 2]에 대해 올바르게 예측되지만 SNH 및 TKN에 대해서는 그렇지 않습니다. SNO의 농도는 증가 할 것으로 예상되며 이는 ASM 솔버에 의해 정확하게 예측됩니다. 모든 폐수 종의 농도는 아래 애니메이션에 표시됩니다.

Zele WWTP에 있는 모든 종의 진화에 대한 시뮬레이션 결과

애니메이션은 Zele WWTP에있는 모든 종의 진화에 대한 시뮬레이션 결과를 보여줍니다.

WWTP 데이터에 대한 결과의 민감도

나는 폐수에서 일부 종의 잘못된 진화를 모델링의 가정과 누락된 WWTP 데이터에 기인합니다. 유입수에서 측정 된 종 농도의 불확실성; 초기 농도에 대한 정보 누락; 그리고 입자상 물질의 침강 특성에 대한 누락 된 데이터는 폐수의 종 농도에 영향을 미쳤을 가능성이 있습니다.

마찬가지로 불완전한 지오메트리 사양은 WWTP 내부의 유체 역학 계산의 정확성에 부정적인 영향을 미칠 수 있습니다. 또한 폭기조에 산소를 살포하는 것에 대한 정보는 부분적으로 만있었습니다. 산소는 다른 종의 부패와 성장에 큰 영향을 미치는 중요한 구성 요소입니다.

WWTP의 모든 데이터를 항상 측정 할 수있는 것은 아닙니다. 이러한 경우 보정 된 수치 모델을 가상 실험실로 효과적으로 사용하여 다양한 WWTP 설계를 테스트 할 수 있습니다. 이 사례 연구는 특히 폭기조에서 WWTP의 2 차 처리 부분에서 종의 농도를 추적 할 수 있음을 보여줍니다. 그리고 이것은 유체 역학 효과를 고려하면서 할 수 있습니다. 완전한 WWTP 데이터와 문제 사양이 존재하는 경우 엔지니어와 설계자는 WWTP 플랜트 운영 및 설계 최적화에 대해 더 나은 정보를 바탕으로 결정을 내릴 수 있습니다.

우리는 활성 슬러지 모델을 추가로 개발하고 보정하기 위해 폐수 처리 산업의 연구원 및 전문가와 협력 할 수 있습니다. 귀하의 WWTP 프로젝트 및 연구에 대해 논의하려면 adwaith@flow3d.com 으로 이메일을 보내 주십시오 .

참고 문헌

[1] Henze M., Lossdrecht M.C.M., Ekama G.A., Brdjanovic D., Biological Wastewater Treatment, Principles, Modelling and Design, IWA publishing 2008.

[2] Peterson B., Vanrollenghem P.A., Gernaey K., Henze M. (2002) Evaluation of an ASM-1 model calibration procedure on a municipal–industrial wastewater treatment plant, Journal of Hydroinformatics, 4(1): 15-38.

[3] Henze, M., Grady, C. P. L. Jr., Gujer, W., Marais, G. v. R. & Matsuo, T. (1987) Activated Sludge Model No. 1. IAWPRC Scientific and Technical Reports No. 1. London, UK.

Simulation results from FLOW-3D highlighting the droplet formation and the input pressure pulse

연속 잉크젯 인쇄

Continuous Inkjets

연속 잉크젯 인쇄는 약 150 년 동안 축적 된 기술입니다. 간단히 말해, 프린트 헤드가 작동하면 연속적인 유체 흐름이있는 액적 생성 방법입니다. 이 개념은 1867 년 Lord Kelvin에 의해 처음 특허를 받았지만 80 년 이상 지난 1951 년 Siemens가 최초의 상용 장치를 선보였습니다. 처음에 이 기술은 만료일, 배치 코드, 이름 및 제품 로고와 같은 가변 정보의 비접촉식 고속 인쇄에 사용되었습니다.

물방울 생성

노즐 크기 선택

액적 생성을위한 시스템 매개 변수를 계산하기 위해 Rayleigh 제트 불안정성 이론을 사용할 수 있습니다. 이 이론에 따르면 물방울 형성으로 이어지는 제트 분리에 대한 자극의 최적 파장 (λ)은 대략 다음과 같습니다.

Nozzle size selection
Nozzle size selection

작동 주파수 선택

최적의 드롭 생성 주파수는 최적의 파장에서 직접 계산할 수 있습니다. 위의 이론과 알려진 산업 매개 변수를 사용하여 FLOW-3D 에서 계산 모델을 설정하는 동안 125μm의 노즐 반경과 10kHz의 주파수가 사용되었습니다

FLOW-3D 결과 검증

FLOW-3D 는 강력하고 정확한 표면 장력 모델로 인해 연속 잉크젯 인쇄와 같은 액적 기반 공정을 시뮬레이션하는 데 적합합니다.

아래 시뮬레이션 결과에서 10kHz의 주파수에서 진동하는 입력 압력 펄스를 볼 수 있습니다. 평균 액적 크기는 약 240 μm이며 이론적으로 추정 된 액적 크기 약 250 μm와 잘 일치합니다.

Simulation results from FLOW-3D highlighting the droplet formation and the input pressure pulse
Simulation results from FLOW-3D highlighting the droplet formation and the input pressure pulse

OLED Mura Problem

이론적으로는 정확히 동일한 진폭으로 압력 펄스를 생성 할 수 있습니다. 그러나 OLED의 잉크젯 인쇄와 같은 산업 응용 분야에서 모든 노즐은 본질적으로 불완전한 제조 또는 작동 매개 변수로 인해 약간 다릅니다. 이러한 모든 결함은 액적 부피의 변동을 일으켜 OLED 패널의 각 하위 픽셀에 증착 된 유기 화합물의 부피를 변화시켜 증착 된 필름 두께의 비례적인 변화를 초래합니다. 이러한 두께 변화는 잉크젯 인쇄 OLED 디스플레이에서 패널 휘도 불균일의 가장 중요한 원인 중 하나입니다 (Madigan et al. ). 이러한 패널 휘도의 불균일성을 “무라 효과”라고합니다.

무라 문제를 해결하는 한 가지 접근 방식은 평균 법칙을 사용하는 것입니다. 이것이 의미하는 바는 서로 다른 노즐 (픽셀 내 혼합)의 방울을 무작위로 결합하여 방울 부피의 양 및 음 오류를 평균화하여 방울 부피 오류를 거의 0에 가깝게 만드는 것입니다.

FLOW-3D 에서 픽셀 내 혼합 과정을 시뮬레이션하기 위해 입력 압력 펄스 진폭에 약간의 임의성이 추가되었습니다. 최대 변동의 크기는 1.7MPa의 원래 압력 진폭에 더하여 200kPa로 설정되었습니다. 아래 애니메이션은 무작위성이있는 케이스와 무작위성이없는 초기 케이스의 비교를 보여줍니다.

압력 펄스의 무작위성 대 일정한 진폭의 경우를 비교하는 애니메이션.

예상대로 액적 생성은 액적 모양, 액적 크기, 액적 간 간격 및 비행 속도 측면에서 균일하지 않습니다. 그러나 오른쪽의 일정한 진폭 케이스는 균일 한 모양과 크기의 균일 한 간격의 물방울을 생성합니다.

연속 잉크젯 인쇄는 저장소에서 마이크로 미터 크기의 노즐 뱅크로 액체를 보내는 고압 펌프로 시작하여 진동하는 압전 결정의 진동에 의해 결정되는 주파수에서 연속적인 물방울 흐름을 생성합니다. 특히 인쇄 응용 분야의 경우, 잉크 방울은 외부 전기장의 존재로 인해 연속 흐름에서 편향됩니다. 이것은 인쇄 매체의 표면에 패턴을 생성합니다. 이 기술의 장점 중 일부는 높은 처리량, 높은 액적 속도, 프린트 헤드에서 기판까지의 거리 증가, 연속 작동으로 인한 노즐 막힘 없음입니다. 이러한 긍정적 인 특성 덕분에이 기술은 오늘날 종이에 일반 인쇄 잉크에서 다양한 재료 (생존 세포 포함)를 증착하는 것으로 발전했습니다.

Continuous inkjet animation

결론

FLOW-3D 는 연속 잉크젯 인쇄 프로세스와 관련된 물리학에 대한 이해를 촉진하는 데 사용되었습니다. 강력한 표면 장력 모델 덕분에 FLOW-3D 는 다양한 고급 액적 생성 및 증착 응용 분야에서도 유용 할 수 있습니다. 예를 들어 OLED 프린팅의 경우 FLOW-3D 를 사용하여 픽셀 내 혼합 중에 발생하는 액 적의 변화를 효과적으로 이해하여 OLED 패널의 품질을 높일 수 있습니다.

References

Madigan C. F., Hauf C. R., Barkley L. D., Harjee N., Vronsky E., Slyke S. A. V., Advancements in Inkjet Printing for OLED Mass Production. Kateeva, Inc.

Gravitational and sedimentation microfluidic technique (Huh et al. Anal Chem 2007)의 중력 회로도를 사용한 입자 분류

중력을 사용한 미세 유체 입자 분류

Microfluidics Particle Sorting Using Gravity

미세 유체 입자 분류는 진단, 화학적 및 생물학적 분석, 식품 및 화학 처리, 환경 평가에 적용됩니다. 이전 블로그에서 유체 역학을 사용한 미세 유체 입자 분류에 대해 이야기했습니다 . 같은 주제를 바탕으로 중력을 사용하여 미세 입자를 분류하는 또 다른 방법에 대해 논의하겠습니다. 아래 애니메이션에서 볼 수 있습니다.

유비쿼터스 중력(Ubiquitous gravity)은 미세 유체 장치에서 미세 입자를 분류하는 데 사용할 수 있습니다. 중력이 입자의 움직임에 수직으로 작용할 때 입자는 반경에 따른 속도로 안정됩니다. 또한 입자의 운동은 입자의 밀도, 유체의 밀도 및 유체의 점도 사이의 차이에서 비롯된 유체 역학적 효과의 영향을받습니다. 아래 이미지는 중력 분류 기술 회로도를 보여줍니다.

Gravitational and sedimentation microfluidic technique (Huh et al. Anal Chem 2007)의 중력 회로도를 사용한 입자 분류
Gravitational and sedimentation microfluidic technique (Huh et al. Anal Chem 2007)의 중력 회로도를 사용한 입자 분류

부력 대 항력

앞서 언급했듯이 중력은 서로 다른 입자가 서로 다른 속도로 침전되도록합니다. 모든 입자의 밀도가 같고 입자 밀도가 주변 유체의 밀도보다 낮 으면 부력 우세와 항력 우세라는 두 가지 유형의 분류를 사용할 수 있습니다. 반경이 더 큰 입자는 더 많은 부력을 경험하고 작은 입자 위의 경로를 따르는 경향이 있습니다. 그러나 외장 액체 (입자를 운반하는 용액)의 유입 속도가 충분히 높으면 항력 효과가 우세하기 시작하고 더 큰 입자가 더 작은 입자의 경로 아래로 이동하는 경향이 있습니다.

FLOW-3D 시뮬레이션 결과

경쟁하는 부력과 항력은 아래 FLOW-3D 에서 얻은 시뮬레이션 결과에서 명확하게 볼 수 있습니다 . 그림 1은 부력 지배적 인 입자 분류의 경우를 보여줍니다. 더 큰 (빨간색) 입자는 수평 채널의 상단을 향해 정렬됩니다. Fig. 2에 나타난 결과는 부력이 우세한 경우의 유입 초 속도를 20 배로 설정 한 후 얻은 것이다. 더 높은 입구 속도에서 더 큰 입자는 더 많은 운동량을 전달하므로 그 위치는 수직 부력의 영향을받지 않습니다. 따라서 입자는 수평 채널의 상단으로 올라가지 않습니다. 대신 그들은 계속해서 바닥으로 이동합니다.

부력

Buoyancy dominant sorting
Buoyancy dominant sorting

Drag

Figure 2. Drag dominant sorting
Figure 2. Drag dominant sorting

LOW-3D 의 입자 모델은입자 분류 또는 기타 입자 역학과 관련된 미세 유체 시뮬레이션에 성공적이고 쉽게 사용할 수 있습니다. 지금까지 우리는 FLOW-3D 의 입자 모델을사용하여 두 가지 입자 분류 기술을 보았습니다. 하나는 유체 역학을 사용하고 다른 하나는 중력을 사용합니다.

중간 단계, 땅은 여전히 ​​대부분 물에 젖은 상태임

폭우에 따른 홍수 시뮬레이션

Flash Flood Simulation

최근에는 우리나라에서도 국지성 폭우가 빈발하고, 기상 이변에 따라 단시간의 폭우에 의해 돌발 홍수가 발생하고 있습니다.

FLOW-3D를 이용한 수치해석으로 홍수 발생시 주요 피해지역이 어떻게 분포될지, 상류피해 영역과 하류피해 영역을 사전에 검토하여 피해를 최소화할 수 있는 의사결정에 도움을 줄 수 있습니다.

FLOW-3D 사용자들이 이미 알고 있듯이, 우리는 보통 극단적인 이벤트 모델링, 복잡한 자유 표면 등을 매우 잘 예측합니다. 이상하게도 첫 번째 질문은 강수량을 모형화할 수 있는가 하는 것이었습니다. 역사적으로 우리의 소프트웨어는 홍수 평야/방류형 문제에 특별히 사용되지 않았습니다. 2D 솔루션이 대부분 잘 작동하기 때문에 이것을 해결하기 위해 일주일 동안 아침을 먹었는데, 돌발 홍수의 경우, 내 결론은 우리가 그것을 매우 잘 한다는 것이다.

돌발홍수 연구에 사용된 지형
돌발홍수 연구에 사용된 지형, 상부 층은 다공성이고 하부 층은 기반암임

여기에 사용된 두 가지 주요 모델은 매스 소스 와 다공성 매체 모델 입니다.

강우전 초기상태
강우전 초기상태

이 시뮬레이션을 설정하는 방법은 다공성 매체인 상부층이있는 지형이 있다는 것입니다. 이 층은 불투과성 기반암과 겹칩니다. 또한 상부 층은 상부 표면의 유체 공급원으로 정의됩니다.

이 시뮬레이션에서 180분 동안 지속되는 강수와 함께 돌발 홍수를 모델링하고 있습니다. 하층은 기반암이며 초기에 예상 한대로 흐름은 상층, 투과성 층으로의 침투와 일반적인 이류 / 네비어스톡스 현상에 의해 제어됩니다. 투과 층이 포화되면 지표수 현상이 더 많이 나타나기 시작하며, 이는 차례로 협곡의 더 낮은 고도에서 급류 흐름으로 이어집니다.

중간 단계, 땅은 여전히 ​​대부분 물에 젖은 상태임
중간 단계, 땅은 여전히 ​​대부분 물에 젖은 상태임

노력의 결과로 중요하고 실행 가능한 정보를 쉽게 시각화할 수 있습니다. 이 경우 우리는 두 가지 상황을 동시에 시뮬레이션했습니다: 협곡(도메인 오른쪽)에서 형성된 런오프와 후속 급류 그리고 미니 캐니언에서 인접한 이미 존재하는 수역(도메인 왼쪽)으로의 런오프.

강수량 맵은 전체 영역에 걸쳐 공간적으로 또는 일시적으로 모두 완전히 정의될 수 있으며, 다공성 행동은 포화 매체와 비포화 매체 모두에 대해 모델링할 수 있습니다.

폭우 후 급류 형성
폭우 후 급류 형성. 오른쪽의 협곡에서 물이 고이고 주요 급류가 형성되고 왼쪽의 기존 수역으로 무거운 방류 줄기가 있습니다.

마지막으로 큰 부피의 흐름과 물이 고이는 것을 볼 것으로 예상했던 하단 (고도 측면) 근처에 위치한 프로브에 위치한 수심의 시계열을 보는 것은 매우 흥미 롭습니다. 

처음에는 다공성 층이 물을 적극적으로 흡수하는 동안 지표수 풀링이 보이지 않습니다. 그런 다음 층이 포화되면 강수 유입과 급류를 통한 유출 사이의 정상 상태 균형을 관찰 할 때까지 수심의 급격한 증가를 관찰됩니다.

마지막으로 영역을 다양한 하위 구성 요소로 타일링하여 강수율 매핑을 훨씬 더 정교하게 만들 수 있습니다.

CFD가 레이저 용접을 만나면 : 불꽃이 어떻게 날아갑니까?

Pareekshith Allu Senior CFD Engineer | Additive Manufacturing | Laser Welding | Business Development

When CFD meets laser welding: How sparks fly!

CFD 또는 전산 유체 역학은 수치적 방법을 사용하여 유체 흐름을 연구하는 것입니다. 유체 흐름의 기본 방정식에는 솔루션 해가 없으므로 컴퓨터를 사용하여 방정식을 반복적으로 계산하는 수치해석 방법으로 해결합니다. 일반적으로 CFD 도구는 공기 역학, 엔진 연소, 물 및 환경 흐름, 미세 유체 및 제조 공정에서 광범위한 연구 및 엔지니어링 문제에 적용될 수 있습니다. CFD가 개발에 중요한 역할을 한 기술을 매일 접할 가능성이 있습니다. FLOW-3D 소프트웨어 제품 제조업체인 Flow Science Inc.에서는 자유 표면 흐름 문제 라고하는 특수한 문제 해결에 중점을 둡니다 . 

자유 표면 흐름이란 무엇입니까? 밀도 차이가 큰 두 유체간에 인터페이스가 공유되는 분야는 자유 표면 흐름입니다. 예를 들어, 기체-액체 경계면이 제한되지 않고 시간에 따라 자유롭게 움직이고 변경할 수 있다는 점에서 강의 물과 주변 공기 사이에 자유 표면이 존재합니다. FLOW-3D 솔버의 기본 DNA 인 Volume of Fluid 또는 VoF 방법 은 자유 표면의 진화를 추적하는 강력한 계산 기술입니다. 우리는 지난 40 년 동안 이 문제에 거의 전적으로 집중했습니다.

자유 표면 흐름은 제조산업 분야에서도 널리 사용됩니다. 금속 주조에서는 용융 금속과 용융 금속이 채우는 금형 또는 다이의 공기 사이에 자유 표면이 존재합니다. L-PBF ( Laser Powder Bed fusion) 라고하는 적층 제조 공정에서 레이저를 사용하여 분말 입자를 녹이고 융합하여 공정에서 자유 표면 용융 풀을 만듭니다. 그리고 레이저 용접에서는 레이저 빔에 의해 녹아서 두 개의 금속 부품 / 부품을 함께 융합 할 때 형성되는 자유 표면 용융 풀이 있습니다. 

이 게시물에서는 레이저 용접 공정에 대한 CFD 시뮬레이션이 유용한 이유를 설명합니다.

레이저 기술은 지난 몇 년 동안 상당히 발전했으며 이제 다른 레이저 제조업체는 다양한 파장에서 펄싱 기능이 있는 고출력 레이저를 제공 할 수 있습니다. 레이저와 로봇 자동화 시스템, 컨트롤러 및 프로세스 센서의 통합은 다양한 제조 산업에서 사용을 확대하여 열 입력이 적고 열 영향 영역이 더 작은 레이저 용접 조인트를 가능하게합니다. 

레이저-재료 상호 작용은 복잡하며이를 정확하게 모델링하려면 이러한 시간적 및 공간적 규모와 관련된 물리학을 구현해야합니다. 레이저 열원은 표면에 에너지를 축적하여 기판을 녹이고 용융 금속 풀을 만듭니다. 용융 풀은 전력, 속도 및 스캔 경로와 같은 레이저 가공 매개 변수와 용융 풀의 자유 표면에 동적 증기압을 적용하는 차폐 가스의 영향을 더 많이받습니다. 또한 용접되는 기판의 재료 특성이 중요한 역할을합니다. 용융된 풀의 상 변화와 증발은 용융 풀을 더욱 압박하는 반동 압력을 유발할 수있는 반면 표면 장력은 풀 내의 유체 대류에 영향을줍니다. 키홀 링이있는 경우 레이저 광선이 키홀 내에 갇혀 추가 반사 영향을 받을 수 있습니다. 기판에 더 많은 에너지를 전달합니다. 불안정한 키홀이 붕괴되면 갇힌 공극이 진행되는 응고 경계에 의해 포착되는 다공성 형성으로 이어질 수 있습니다. 

분명히 많은 일이 진행되고 있습니다. 이것이 CFD 시뮬레이션이 강력 할 수있는 곳이며 FLOW-3D WELD를 개발할 때 레이저-재료 상호 작용을 이해하는 데 많은 노력을 기울이는 이유입니다. 자유 표면 추적 및 레이저 에너지 증착, 차폐 가스 역학, 상 변화, 반동 압력, 표면 장력, 레이저 광선 추적 및 응고와 함께 유체 및 열 흐름 방정식을 통합하는 물리 기반 모델은 레이저의 복잡한 상호 작용을 캡처하는 데 매우 정확합니다. 용접과정을 해석하는 기능은 용융 풀의 안정성에 대한 다양한 공정 매개 변수의 영향을 분리하고 엔지니어와 연구원이 용접 일정을 최적화하는 데 도움이 될 수 있습니다.

CFD 시뮬레이션은 레이저 용접 프로세스를 분석하고 개선하는데 도움이되는 프레임 워크를 제공 할 수 있습니다. 불안정한 용융 풀은 키홀 유발 다공성, 파열 및 스패 터와 같은 결함을 초래할 수 있기 때문에 용융 풀의 작동 방식을 이해하는 것은 조인트의 품질에 매우 중요합니다. 그 후, FLOW-3D WELD 모델의 출력인 응고된 용융 풀 데이터 및 열 구배와 같은 결과를 미세 구조 또는 유한 요소 분석 모델에 입력하여 각각 결정 성장 및 열 응력 진화를위한 길을 닦을 수 있습니다.

이 게시물이 CFD를 사용하여 레이저 용접 프로세스를 시뮬레이션하는 이점을 이해하는데 도움이 되기를 바랍니다.

레이저 용접 공정을 더 잘 이해하기 위해 CFD 시뮬레이션 적용을 고려해 보셨습니까? 어떤 특징 / 물리 현상이 모델링되기를 원하십니까? 질문과 의견이 있으면 언제든지 flow3d@stikorea.co.kr 또는 미국 본사의 paree.allu@flow3d.com에게 연락하십시오.

바람이 개방형 철광석 골재 저장소에 미치는 영향 분석 (비산먼지 배출 방지 연구)

다양한 구성에 대한 비산 먼지 배출

이 기사는 Dhananjay Sharma, EI, CFM, 유압 모델링 엔지니어, AECOM 에 의해 기고되었습니다  .

바람이 개방형 골재 저장소에 미치는 영향은 전 세계적으로 환경 문제가 되고 있습니다. 2.7km2 철골 저장소 부지에서 이런 문제가 관찰되었습니다. 이 시설은 철도 운송차량를 통해 광석을 공급받는데, 이 운송차량은 자동 덤프에 의해 비워집니다. 그런 다음 이 광석은 일련의 컨베이어와 이송 지점을 통과하여 저장 장소중 하나로 운송됩니다. 비산먼지 배출은 풍력이 비축된량에 미치는 영향의 결과로 관찰된 결과입니다.

두 가지 다른 구성(옵션 A와 B)을 FLOW-3D로 모델링하여 비산먼지 배출의 영향을 연구했습니다. 옵션 A에는 4줄에 9개 더미가 있는 36개의 비축량이 있고 옵션 B에는 1줄에 총 16개의 비축량이 있습니다.

또한 장벽이 있는 공기와 장벽이 없는 공기의 속도를 비교하기 위해 비축물 주변을 따라 30미터 높이의 장벽을 모델링할 수도 있습니다. 10m 높이에서 기준 초속 7.5m(m/s)의 풍속이 두 구성을 모두 모델링하는데 사용되었습니다. 비축 옵션 A와 B에 대해 네 가지 풍향 방향이 분석되었습니다.

물리적 및 수치 적 모델링

초기 모델 설정

FLOW-3D 에서 비산 먼지 배출을 모델링하기  위해, 공기 온도는 15 ° C로 가정되었습니다. 단일의 균일한 비압축성 유체 옵션이 선택되었습니다. z 방향에서 -9.81 m / s의 중력이 사용되었습니다. 유체는 점성과 난류로 간주되었습니다. 2- 방정식 (ke) 모델은 옵션 A 및 B 구성 모두에 대해 표면 마찰없이 난류를 계산하는데 사용되었습니다.

초기 조건

1/7 power 법칙 (pproximately a logrithmic law-of-the-wall distribution)에 기반한 속도 프로파일이 각 시뮬레이션에 대한 초기 조건으로 지정되었습니다. 비축 분석에서 가장 관심있는 기준 속도는 12 및 7.5m / s입니다. 풍속을 증가시키고 파일에 인접한 속도에 미치는 영향을 측정하여 분석을 수행했으며, 레이놀즈 스케일링이 이러한 속도에 대해 유지된다는 것을 확인했습니다 (즉, 들어오는 풍속 스케일링과 파일에 인접한 속도 스케일링 간의 선형 관계).   그런 다음 7.5m / s의 속도 만 사용하여 FLOW-3D 시뮬레이션을 구성했습니다. 이러한 시뮬레이션의 결과는 12m / s 조건을 충족하도록 확장 할 수 있습니다.

풍력 프로파일 power 법칙을 사용하여 10m에서 7.5m / s 이상 및 이하의 다양한 높이에 대한 속도를 추정했습니다. 경계에서 속도를 적용하는 이 방법은 경계를 따라 지형 변화를 허용하지 않습니다. 기준 속도는 서쪽, 남서부 및 남풍 방향에 대해 해발 10 미터에서 할당되었습니다. 동풍의 경우, 속도는 뒤쪽 (Y- 최대) 경계에서 경사 10 미터 위의 기준 높이에 할당되었습니다.

풍력 프로파일 power 법칙은 z 방향으로 최대 360m까지 모든 미터에서 계산되었습니다. 속도는 메쉬 크기와 동일한 간격으로 평균화되었습니다. 속도가 할당된 높이 간격은 2, 4, 6, 8, 10, 20, 70, 181, 270 및 360 미터입니다. 속도 프로파일을 설정 한 후 각 높이 간격에 대한 값은 네 가지 풍향 (서쪽, 남서쪽, 남쪽 및 동쪽) 각각에 대해 X 및 Y 구성 요소로 세분화되었습니다. 초기 조건은 메쉬 블록의 외부면에 할당되어 비축에 도달하기 전에 속도 프로파일이 개발 될 수있는 충분한 수평 공간을 남겼습니다.

풍력 프로필 power 법칙은 다음과 같습니다.

\ displaystyle {{u} _ {x}} = {{u} _ {r}} {{\ left ({\ frac {{{{z} _ {x}}}} {{{{z} _ { r}}}}} \ right)} ^ {\ propto}}, 여기서

U x  = 높이에서의 풍속 x
U r  = 기준 높이에서의 풍속
Z x  = 높이 x
Z r  = 기준 높이
α = 1/7 ‐ 대기 안정성 계수

지형

3 개의 지형파일인 스테레오리소그래피 (STL) 파일이 생성되어 모델에 통합되었습니다. 개별 파일은 지형, 창고 및 기둥에 해당합니다. 옵션 A와 B에 대해 다른 STL 파일이 생성되었습니다.

메싱

모델 도메인은 각 풍향에 대해 조정되었습니다. 메쉬 크기는 옵션 A의 경우 240 만에서 330 만 셀, 옵션 B의 경우 130 만 셀입니다. 정확하게 기둥 근처에 높이 2m, 길이 4m, 너비 4m의 셀 크기를 사용했습니다. 해당 지역의 속도를 계산합니다.

경계 조건

비축 시뮬레이션에는 네 가지 경계 유형이 사용되었습니다. 모든 풍향에 대해 상단 경계 (Z-max)가 정체 압력으로 지정되었습니다. 바람의 방향에 따라 두 개의 측벽이 유출 경계 조건으로 지정되었습니다. 나머지 두 측벽에는 그리드 오버레이 경계가 지정되었습니다. 그리드 오버레이를 사용하면 초기 조건의 속도를 모델에 입력 할 수 있습니다. 중첩 된 블록을 사용하여 원하는 메시 해상도와 배율을 만들었습니다. 내포된 블록 사이의 경계면에서 대칭 경계 조건이 사용되었습니다. 대칭을 사용하면 블록간에 정보를 전송할 수 있습니다. 그림 1은 서쪽 풍향 (y 방향)에 대한 경계 조건 설정을 보여줍니다. 다른 풍향의 경우 경계 조건을 적용하는 데 유사한 방법이 사용되었습니다.

경계 조건 서쪽 풍향
그림 1. 서쪽 풍향의 경계 조건

장벽

FLOW-3D 의 배플 기능은 비축된 곳의 주변에 바람 장벽을 만드는데 사용되었습니다. 옵션 A와 B의 배플은 높이가 30 미터였으며 지형을 따라 여러 부분으로 구성되었습니다. 모델링된 장벽은 본질적으로 다공성입니다. 34 %의 다공성 값 (즉, 34 % 개방 면적) 및 해당 속도 대 압력 강하 값은 장벽 제조업체에서 얻었습니다. FLOW-3D의  모델과 연관된 흐름 다공성 손실이 지정될 수있는 배플 알고리즘을 사용합니다. 배플은 무한히 얇고 부피를 차지하지 않습니다.

시뮬레이션 결과

옵션 A

옵션 A의 경우 풍속 7.5m / s에 대한 장벽이 있거나없는 4 가지 풍향을 분석하고 시뮬레이션했습니다.

바람의 방향배리어없는 최대 속도 (m / s)배리어가있는 최대 속도 (m / s)최대 속도 감소
서부13.58611.27817 %
남서부13.04510.79617 %
남쪽12.35212.122 %
동쪽9.768.59712 %

각 시뮬레이션의 최대 속도와 장벽과 장벽이 없는 경우 사이의 최대 속도 감소는 위의 표 1에 나와 있습니다. 장벽은 남풍의 최대 속도에 가장 적은 영향을 미칩니다. 옵션 A에 대한 장벽 추가로 최대 속도가 2 % 감소했습니다. 장벽은 서풍 또는 남서풍이있는 전체 파일 케이스의 속도에 가장 큰 영향을 미쳤습니다. 최대 속도는 서풍과 남서풍 모두에서 17 % 감소했습니다.

옵션 B

옵션 B의 경우 풍속 7.5m / s에 대한 장벽이 있거나없는 네 가지 풍향을 분석하고 시뮬레이션했습니다.

그림 2. 옵션 A : 장벽이없는 서풍의 비축량에서 계산 된 속도 크기
장벽이있는 속도 크기 서풍
그림 3. 옵션 A : 장벽이있는 서풍 방향의 비축에서 계산 된 속도 크기
바람의 방향배리어없는 최대 속도 (m / s)배리어 포함 최대 속도 (m / s)최대 속도 감소
서부15.9711.3629 %
남서부15.149.2139 %
남쪽13.410.124 %
동쪽12.787.1544 %
그림 4. 옵션 B. 장벽이없는 동풍의 비축량에서 계산 된 속도 크기
그림 5. 옵션 B : 장벽이있는 동풍의 비축에서 계산 된 속도 크기

결론

모델 결과는 비축물 주변에 장벽을 추가하는 것이 속도를 줄이고 비산먼지 배출을 방지하는데 도움이 된다는 것을 분명히 보여주었습니다. 현장 주변의 장벽 추가와 관련된 비용이 있지만, 이 옵션은 먼지 배출량을 줄임으로써 환경 규범을 달성하는 데 도움이 될 것입니다. 모델 결과를 보면 FLOW-3D가 비산먼지 방출을 연구하기 위한 정확하고 신뢰할 수 있는 도구로 사용될 수 있다는 것이 분명합니다. 추가 설계 변경과 철골 배치의 새로운 옵션이 제안될 경우 FLOW-3D에서 쉽게 모델링하여 비용 및 환경적으로 효과적인 최적의 구성을 결정할 수 있습니다.

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

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

선택적 레이저 용융법(SLM: Selective Laser Melting)은 3D 프린팅 기술의 하나로 최근 주목 받고 있습니다. SLM에서는 레이저 조사 중 높은 온도 구배로 인해, 용융과 재응고 현상이 일어나므로 용융금속 유체의 거동이 중요한 역할을 담당하고 있어, 구성 부품의 최종 구조를 결정합니다.

FLOW-3D@ WELD를 이용하여 T-6Al-4V(64티타늄 합금)에 대한 선택적 레이저 용융법 (SLM) 시뮬레이션이 가능합니다.

SLM 개념도
SLM 개념도

금속 분말을 얇게 깔아 생긴 분말층에 레이저를 조사하면 조사된 부분만 용융, 응고 됩니다. 이 공정을 반복하면서 적층하여 3차원 형상을 만듭니다. 금속을 재료로 하여 고강도 제품을 만들수 있으므로, 기존의 시작 제품(Rapid Prototyping)뿐만 아니라, 짧은 납기일, 저비용, 고기능 등을 목적으로 한 Additive Manufactuing 기술로서 주목받고 있습니다.

FLOW-3D@ WELD를 이용한 해석을 통해서, 표면의 경사에 따라 용융지의 형상과 온도 분포가 결정된다는 것을 알 수 있습니다.

용융 풀의 최대 깊이는 SLM의 형태학적 변화에 따라 달라지며 평균 깊이는 42μm입니다.

선택적 레이저 용융법 (SLM) 해석 결과
선택적 레이저 용융법 (SLM) 해석 결과

 *Source: National Cheng Kung University, Department of Materials Science and Engineering, Taiwan YC Wu, WS Hwang

Figure 4 A view of the ogee spillway and Type 2 piers in the 3D CFD model

NUMERICAL ANALYSIS AND THE REAL WORLD : IT LOOKS PRETTY BUT IS IT RIGHT?

D. K. H. Ho, S. M. Donohoo, K. M. Boyes and C. C. Lock
Advanced Analysis, Worley Pty Limited
L7, 116 Miller Street, North Sydney, NSW 2060 Australia
Tel: +61 2 8923 6817 e-mail: david.ho@worley.com.au

Abstract

엔지니어링 설계에서 유한 요소, 유한 차분 및 전산 유체 역학 분석 소프트웨어와 같은 수치 도구의 일상적인 사용이 최근 몇 년 동안 증가했습니다. 소프트웨어 및 하드웨어 기술의 발전은보다 비선형적이고 복잡한 3 차원 분석이 수행되고 있음을 의미합니다.

그러나 본질적으로 “블랙 박스”인 이러한 강력한 소프트웨어는 “컴퓨팅”기술을 보유하고 있지만 광범위한 엔지니어링 경험이 필요하지 않은 분석가의 손에 “컴퓨터 보조 재해”로 이어질 수 있습니다. 품질 보증 절차의 엄격한 구현은 수치 모델이나 분석 기법이 정확한지 확인할 필요가 없을 수 있습니다.

이 백서에서는 복잡성이 증가하는 세 가지 실제 토목 공학 응용 프로그램에서 수치 분석 결과를 검증하는 방법을 설명합니다. 여기에는 유한 요소법을 이용한 수조 탱크의 구조 해석, 전산 유체 역학법을 이용한 수력 구조물 위의 홍수 조사, 유한 ​​차분법을 이용한 안벽 시공 시뮬레이션 등이 있습니다. 입력 데이터의 불확실성 수준과 각 사례에 대한 계산 결과의 신뢰성에 대해 논의합니다. 분석 과정에서 몇 가지 흥미로운 결과가 발견되었습니다.

첫 번째 사례 연구는 시공의 질이 구조물의 성능에 상당한 영향을 미친다는 것을 보여주었습니다. 그러나 설계자는 설계 단계에서 이러한 상황을 수량화하고 분석하지 못할 수도 있습니다. 필요할 경우 향후 역분석은 물론 설계 검증의 기준점이 될 수 있도록 공사 종료 시 모니터링의 중요성이 필수적입니다. 유한 요소 분석은 복잡한 문제를 분석할 수 있는 강력한 수치 도구이지만, 분석가들은 문제의 행동이 단순하고 잘 이해되는 것처럼 보일 수 있는 상황에서 예상치 못한 결과를 만날 수 있도록 준비해야 합니다.

두 번째 사례 연구에서는 중요한 배수로 구조에 전산 유체 역학 분석이 처음으로 적용 되었기 때문에 엄격한 검증 프로세스가 강조됩니다. 그것은 2D ogee 방수로 프로파일로 시작하여 문제의 방수로의 3D 모델을 분석하기 위해 진행되는 방식으로 수행되었습니다.
계산된 결과를 각 단계에서 이론 및 물리적 테스트 데이터와 비교했습니다. 유체 흐름 문제의 비선형적 특성에도 불구하고, 분석은 확신을 가지고 실제 설계 목적에 적합한 결과를 제공할 수 있었습니다.

최종 사례 연구에서는 안벽의 거동이 시공 이력과 매립 방식에 영향을 받은 것으로 나타났습니다. 벽의 움직임은 매우 가변적인 토양 속성에도 불구하고 질적으로도 단순한 비선형 토양 모델을 사용하여 정확하게 예측되었습니다. 지속적인 모니터링 기록이 없기 때문에 검증은 어려웠습니다. 계산된 결과를 검증하는 열쇠는 수치 소프트웨어 도구를 사용하지 않는 독립적인 계산을 찾는 것입니다. 대부분의 경우 이러한 솔루션을 사용할 수 있습니다. 그러나 다른 경우에는 실험실 또는 현장 관찰에만 의존할 수 있습니다.

Introduction

오늘날 수치 해석은 대부분의 엔지니어링 설계에서 필수적인 부분을 형성합니다. 따라서 결과 검증의 필요성은 분석 기술 / 방법론을 신뢰할 수 있고 설계자가 계산 된 결과에 대한 확신을 가질 수 있도록 설계 프로세스 전반에 걸쳐 매우 중요합니다.

일반적인 관행은 고전 이론, 실험 데이터, 게시 된 데이터, 유사한 구조의 성능 및 다른 사람이 수행 한 수치 계산에 대해 결과를 검증하는 것입니다. 때때로 소프트웨어 개발자가 제공 한 벤치 마크 또는 검증 예제가 이러한 목적으로 사용될 수 있지만 전체 범위의 문제를 포괄 할만큼 포괄적 인 경우는 거의 없습니다.

수치 해석을 시작하기 전에 분석가는 입력 데이터의 신뢰성, 소프트웨어 도구가 문제의 문제를 해결할 수 있는지 여부 및 결과를 검증하는 방법을 결정해야합니다. 검증 프로세스가 많은 실무자들에 의해 품질 보증 절차의 일부로 채택되었지만 비용이 많이 드는 실패가 여전히 발생했습니다 [1].

Validation

결과 검증의 필요성은 수치 분석의 사용 (남용)에서 일부 나쁜 업계 관행을 관찰함으로써 강화 될 수 있습니다. 수치 계산을 수행하기 위해 고용 된 일부 엔지니어 / 분석가는 계산 뒤에있는 기본 이론을 완전히 이해하지 못하거나 숨겨진 함정을 처리 할 수있는 실제 엔지니어링 경험이 충분하지 않을 수 있습니다.

일부 소프트웨어가 “CAD와 유사”해지고 많은 사람들이 작동하기 쉽다고 주장하기 때문에 엔지니어링 회사가 대학원 엔지니어 대신 초보를 고용하여 수치 모델링 및 분석을 수행하는 경향이 점차 증가하고 있습니다.

사용자는 복잡한 지오메트리 모델을 생성하고, 적절한 요소와 메시를 만들고, 각 하중 케이스에 대한 경계 조건 (접촉, 하중 및 고정)을 적용하고, 속성을 할당하고, 제출에 필요한 모든 플래그 / 스위치 / 버튼을 설정하는 데 상당한 노력을 기울일 것입니다.

분석이 실행됩니다. 자체 검사를위한 일부 품질 보증 절차는 전처리 단계에서 따를 수 있지만 계산이 완료되고 결과가 후 처리 될 때까지 많은 사용자는 출력이 어느 정도 정확하다고 쉽게 믿을 것입니다. 지오메트리 생성은 수치 모델링 프로세스의 일부일뿐입니다. 가장 어려운 문제 중 하나는 전체 설계 프로세스에서 불확실성을 다루는 것입니다. 재료 속성 및 로딩 순서와 같은 입력과 관련된 불확실성이 있습니다.

예를 들어 모델이 선형 또는 비선형 방식으로 동작하는지 여부와 같이 솔루션 유형의 적절성과 관련된 불확실성이 있습니다. 마지막으로 결과 해석과 관련된 불확실성이 있습니다. 수치 분석에서 결과를 검증하고 문제를 발견하는 데있어 분석가를위한 좋은 방법에 대한 간단한 지침은 없습니다. 그러나 다음 방법을 통해 점차적으로 달성 할 수 있습니다.

• 수치 적 방법 과정에 대한 좋은 이해 – 이것은 학부 및 / 또는 대학원 수준의 공식 교육을 통해 얻을 수 있으며 지속적인 전문성 개발의 일환으로 자습을 통해 더욱 향상 될 수 있습니다.
• 특정 유형의 문제에 대한 기본 이론과 해결책의 범위를 잘 이해합니다. 이 역시 위와 같은 교육을 통해 이루어질 수 있습니다.
• 실제 문제를 해결하는 데 공학적 판단을 사용하고 수치 분석을 수행 한 경험이 있습니다. 이는 숙련 된 엔지니어가 분석가를 적절하게 감독하는 환경에서 작업함으로써 얻을 수 있습니다.

품질 보증 시스템의 구현이 실행 가능한 솔루션으로 이어지는 엔지니어링 판단을 대체하는 것은 아니라는 점에 유의해야합니다. 복잡한 대규모 모델을 분석하기 전에 시뮬레이션 기술과 문제의 근본적인 동작을 완전히 이해하기 위해 간단한 테스트 모델을 사용하여 수치 “실험”을 수행해야하는 경우가 매우 많습니다.

경험에 따르면 때때로 테스트 모델 자체가 분석가가 최종 설계 솔루션에 도달 할 수있는 충분한 정보를 제공 할 수 있습니다. 해당 대형 복합 모델의 분석은 설계 기대치를 확인하는 것입니다. 다음 사례 연구는 결과 검증이 수행 된 방법과 신뢰 수준 및 불확실성이 해결된 방법을 보여줍니다.

Applications

일반적인 토목 공학 프로젝트에서 수치 분석은 구조 역학, 기하학 및 유체 역학의 세 가지 기본 분야 중 하나 또는 조합을 포함 할 수 있습니다. 문제의 성격은 토양-구조 상호 작용, 유체-구조 상호 작용 또는 토양-유체 상호 작용 중 하나로 분류 될 수 있습니다.

어떤 경우에는 세 가지 모두를 포함 할 수 있습니다. 잠재적 인 복잡성을 고려하여, 정확도를 잃지 않고 실제 목적을 위해 중요한 동작을 캡처하지 않고 문제를 단순화하기 위해 몇 가지 가정과 이상화가 이루어져야합니다. 이러한 문제를 해결할 수있는 범용 및 특수 수치 분석 소프트웨어가 있습니다. 두 가지 유형의 소프트웨어가 사례 연구에 사용되었습니다.

Case 1 – Deflection of a steel water tank

직경 약 90m의 대형 원형 강철 물 탱크는 처음 채울 때 큰 벽면이 휘어지면서 탱크의 장기적인 구조적 무결성에 대한 우려를 불러 일으켰습니다.

물의 높이는 전체 저장 용량에서 약 10m였습니다. 지붕 구조는 탱크 내부에있는 기둥으로 거의 전적으로지지되었습니다. 스트레이크(strakes)는 벽의 바닥 1/3이 더 두꺼운 고급 강판으로 구성되었습니다. 1 차 윈드 거더는 탱크 상단 주위에 용접되었고 2 차 윈드 거더는베이스 위 2/3에 위치했습니다. 하단 스트레이 크는 환형베이스 플레이트에 필렛 용접되었습니다. 내부 기둥의 기초를 제외한 전체 바닥은 용접 된 강판으로 덮여있었습니다.

이 탱크는 유능한 중간층 사암과 미사암 기반암 위에 압축된 채움물 위에 세워졌습니다. 일련의 축 대칭 유한 요소 분석 (FEA)을 수행하여 관찰된 처짐을 예측할 수 있는지 여부를 결정하고 매일 물을 채우고 비울 때 피로 파괴가 발생할 가능성으로 인해 벽 바닥의 응력 상태를 계산했습니다.

내부 기둥과 지붕 빔을 포함하는 탱크의 12 분의 1 섹터에 대한 3 차원 모델을 처음에 분석하여 벽이 얼마나 많은 지붕 자중을지지하고 축 대칭 가정의 타당성을 조사했는지 조사했습니다. 이 분석의 결과는 지붕 구조의 강성 기여도가 중요하지 않아 후속 축 대칭 모델에 포함되지 않았 음을 보여주었습니다.

그러나 지붕 자체 무게의 작은 부분이 벽에 적용됩니다. 축 대칭 모델은 모든 강철 섹션, 필렛 및 맞대기 용접 및 기초로 구성되었습니다 (그림 1). 그것들은 몇 개의 3 노드 삼각형 축 대칭 요소가있는 4 노드 비 호환 모드 사변형으로 이산화되었습니다.

용접 재료를 통해서만 하중 전달이 허용되도록 용접이 모델링되었습니다. 용접 연결부에 미세한 메시를 사용하여 응력 상태를 정확하게 포착했습니다. 롤러 지지대는 모델의 측면 및 하단 경계에 적용되었습니다. 다음과 같은 하중이 적용되었습니다 :

철골 구조물의 자중, 지붕 자중, 벽의 정수압, 수위에 따른 바닥의 균일 한 압력. 한 모델은 용접 또는베이스의 강판이 플라스틱 힌지를 형성하기 위해 항복되었다고 가정했습니다. 이 경우 벽 바닥에서 핀 연결이 모델링되었습니다.

Partial FE mesh of tank/foundation. Insert shows mesh and stress distribution at wall base
그림2 Partial FE mesh of tank/foundation. Insert shows mesh and stress distribution at wall base

벽 처짐은 그림 2에 나와 있습니다. 측정 범위와 계산 된 결과는 비교 목적으로 표시됩니다. 계산 된 벽 처짐을 검증하기 위해 두 벽 두께에 대한 Timoshenko 및 Woinowsky-Krieger [2]에 기반한 고전 이론도 그림에 표시되었습니다. 계산 된 편향은 이론적 계산에 의해 제한됨을 관찰 할 수 있습니다.

벽 두께의 변화로 인한 전이가 분석에서 포착되었습니다. 이것은 유한 요소 모델에 대한 확신을 제공했습니다. 윈드 거더와 구속 된베이스의 영향도 볼 수 있습니다. 윈드 거더 설치로 인해 초기 변형이 발생하여 공사가 끝날 때 벽 상단이 안쪽으로 당겨질 수 있습니다. 굽힘 동작이 발생한베이스 근처를 제외하고는 후프 동작이 벽 동작을 지배했습니다.

계산된 최대 처짐이 측정된 순서와 동일하더라도 최대 돌출이 발생한 높이는 예측되지 않았습니다. 실제로 조사 데이터는 몇 가지 가능한 시나리오를 제안했습니다.베이스에 플라스틱 힌지 형성 (그러나이 영역에서 계산 된 응력은 항복 강도를 초과하지 않았습니다). 지반 재료의 국부적 인 베어링 고장 (다시 현장에서 균열과 같은 명백한 지시 신호가 보이지 않음); 또는 탱크 건설이 끝날 때 내장 된 기하학적 결함이있었습니다. 사전 변형 된 탱크에서 역 분석을 수행하여 측정 된 처짐이 정수압 하에서 “회복”되었습니다. 그러나 계산된 응력은 수율을 훨씬 초과했습니다. 불행히도 탱크는 완성 후 첫 번째 충전 전에 즉시 조사되지 않았습니다.

Figure 2 Wall deflection of water tank
Figure 2 Wall deflection of water tank

탱크의 원래 디자인과 건설이 2000 년대 초에 수행되었다는 점은 흥미 롭습니다. 설계 계산에 관련 표준 [3]을 사용했습니다. 이 표준은 탱크 벽이 후프 동작만으로 작용한다고 가정하고이 구조의 경우가 아닌베이스의 제약 조건을 무시합니다. 벽 처짐의 크기는 기초 강성을 고려한 Rish [4]가 개발 한 고전 이론 [2] 또는 FEA와 같은 수치 분석에 의해 결정될 수 있습니다. 고급 강철을 사용하면 설계자는 강도에는 적합하지만 서비스 가능성에는 필요하지 않은 더 얇은 섹션을 선택해야합니다. 굽힘 강성은 큐브 두께에 의해 결정됩니다. 수중 부하에서 후속 벽 변형 프로파일은 제작 품질에 영향을받습니다. 이것은 설계 단계에서 추정하기 어려웠을 것입니다.

사례 2 – 배수로 배출

호주의 많은 댐 구조는 제한된 수 문학적 정보로 1950 년대와 60 년대에 설계 및 건설되었습니다. 이러한 기존 방수로 구조는 수정 된 가능한 최대 홍수 수준에 대처하기 위해 크기가 작습니다. 증가 된 홍수 조건 하에서 방수로 꼭대기에 대한 음압 생성과 같은 잠재적 인 문제가 발생할 수 있습니다. 이는 방수로 및 게이트 구조에 불안정성 또는 캐비테이션 손상을 유발할 수 있습니다. 역사적으로 스케일링 된 물리적 모델은 이러한 동작을 연구하기 위해 수력 학 실험실에서 구성되었지만 비용이 많이 들고 시간이 많이 걸리며 스케일링 효과와 관련된 많은 어려움이 있습니다. 오늘날 고성능 컴퓨터와보다 효율적인 전산 유체 역학 (CFD) 코드를 사용하여 수리적 구조의 동작을 합리적인 시간과 비용으로 수치 적으로 조사 할 수 있습니다. 이 분석 기법은 대도시 지역에 주요 상수원을 제공하는 가장 큰 콘크리트 중력 댐에 호주에서 처음으로 적용 되었기 때문에 검증을 수행 할 필요가있었습니다. 이것은 그림 3과 같이 조사 프로세스에 통합되었습니다. 순서도는 간단한 2D에서 상세한 3D 방수로 모델로 어떻게 발전했는지 보여줍니다.

Figure 3 Flowchart showing the validation process
Figure 3 Flowchart showing the validation process

미 육군 공병대 [5]에서 발표 한 광범위한 데이터가 있기 때문에 검증을 위해 ogee 방수로 프로필 (그림 4 참조)이 선택되었습니다. 계산 된 결과는 조사의 각 단계에서 검토되었습니다. 게시 된 데이터에서 크게 벗어나면 프로젝트가 중단됩니다. 이것은 프로젝트가 시작되기 전에 고객과 상호 합의되었습니다.

Figure 4 A view of the ogee spillway and Type 2 piers in the 3D CFD model
Figure 4 A view of the ogee spillway and Type 2 piers in the 3D CFD model

이러한 종류의 분석의 초기 어려움 중 하나는 개방 채널 중력 흐름 문제에서 자유 표면의 정확한 계산이었습니다. 자유 표면을 추적하는 데 적응 형 메싱 및 반복 방법을 사용하는 것은 일부 유한 체적 CFD 코드에서 사용되었지만 성공은 제한적이었습니다. 본 연구에 사용 된 코드는 SOLA-VOF 방법으로 Navier-Stokes 방정식을 해결합니다. 유체 운동의 과도 동작을 해결하기 위해 유한 차분 방법이 사용되었습니다. 유체의 부피 (VOF) 함수는 자유 표면 운동을 계산하는 데 사용됩니다 [6].

분석에 대한 자세한 내용은 [7]에 설명되어 있습니다. 계산 된 파고 압력 분포, 자유 표면 프로파일 및 정상 상태에서의 배출 속도는 검증 목적으로 사용되었습니다. 다른 상류 수두 (H) 아래의 배수로 꼭대기를 따라 압력 분포가 그림 5에 나와 있습니다. 일부 압력 진동은 코드가 일반 메시와 곡선 배수로 장애물 사이의 인터페이스에서 계산을 처리하는 방식에 기인 할 수 있습니다. 훨씬 더 미세한 메쉬는 이러한 불규칙성을 부드럽게 만들었습니다. 압력 분포에 대한 교각의 영향은 3D 모델에서 올바르게 예측되었습니다 (그림 6).

계산된 자유 표면 프로파일 (그림 7)도 게시 된 데이터와 잘 일치했습니다. Savage와 Johnson [8]은 분석 기법에 대한 신뢰도를 높이는 동일한 CFD 코드를 사용하여 유사한 유효성 검사를 수행했습니다. 문제의 배수로에 대한 후속 분석은 스케일링 된 물리적 모델 테스트에서 얻은 결과와 비교할 때 상당히 좋은 결과를 제공했습니다.

Figure 5 Comparison of crest pressure for various heads (2D model), Hd is the design head
Figure 5 Comparison of crest pressure for various heads (2D model), Hd is the design head
Figure 6 Comparison of crest pressure next to pier (3D model)
Figure 6 Comparison of crest pressure next to pier (3D model)
Figure 7 Upper nappe profile next to pier
Figure 7 Upper nappe profile next to pier

분석에서 배수로의 기하학적 구조와 물 속성이 잘 정의되었습니다. 물은 비압축성이며 고정 된 온도에서 일정한 특성을 가지고 있다고 가정했습니다. 실제로 좋은 품질의 콘크리트 표면 마감을 얻을 수 있기 때문에 배수로 경계는 매끄럽다 고 가정했습니다. 불확실성은 메쉬 밀도와 적절한 난류 모델의 선택이라는 두 가지 소스에서 비롯됩니다. 메쉬 크기는 메모리 양과 컴퓨터의 클럭 속도에 의해 제한됩니다.

높은 레이놀즈 수의 난류 흐름은 소용돌이와 소용돌이의 형성을 포착 할 수있는 매우 미세한 메시로 계산할 수 있지만 현재 메시 밀도는 검증 및 설계 목적에 필요한 변수를 예측하기에 충분히 미세했습니다. 조사 결과는 큰 와류, k-ε 및 RNG 모델과 같은 난류 모델의 선택에 의해 크게 영향을받지 않는 것으로 나타났습니다. 분명히 벽 거칠기와 난류 모델의 도입은 방전율을 감소시킬 것입니다. 그러나 다시 분석 결과는 사용 된 메시에 거의 영향을 미치지 않음을 보여줍니다. 향후 분석은 다른 메쉬 밀도로 인한 이산화 오류를 조사 할 것입니다.

사례 3 – 안벽 건설
주요 컨테이너 항구 시설은 설계 단계에서 최소한의 수치 분석을 수행하여 약 25 년 전에 건설되었습니다. 당시에는 이러한 분석 도구를 사용하는 것이 비용 효율적이지 않은 것으로 간주되었습니다. 다수의 컨테이너 크레인이 측면을 따라 이어지는 2km 길이의 안벽을 건설하기 위해 광범위한 준설 및 매립 작업이 수행되었습니다.

시설이 완공 된 이후 일련의 콘크리트 카운터 포트 유닛으로 구성된 안벽과 후방 크레인 빔은 크레인이 할 수 있도록 후방 빔에 대한 레벨 조정 작업이 수행 될 정도로 지속적으로 이동하고 있습니다. 정상적으로 작동합니다. 그러나 영향을받는 두 구조물의 움직임을 저지하기 위해보다 영구적 인 해결책을 모색했습니다. 토양-구조 상호 작용 및 시공 시뮬레이션을 처리 할 수있는 명시 적 유한 차이 분석을 사용하여 다양한 교정 옵션의 순위를 지정했습니다.

그라우트 기둥, 타이백 앵커 및 말뚝 지지대와 같은 다양한 제안 된 개선을 분석하기 전에, 토양 및 구조적 특성과 시공 과정의 선택이 적절하도록 계산 모델을 관찰에 대해 보정해야한다고 결정했습니다. 지질 및 지질 공학 정보는 현장 및 실험실 테스트 데이터를 포함하는 현장 조사 보고서에서 평가되었습니다. 시설의 범위를 고려할 때 현장에서 만나는 특정 토양 유형에 대해 상당한 분산 테스트 데이터가 예상됩니다. 수력 모래 충전재에 대한 표준 침투 테스트 (SPT) 블로우 횟수 (N) 및 콘 침투 테스트 (CPT) 저항 (qc)에 대한 몇 가지 일반적인 기록이 그림 8과 9에 나와 있습니다.

Figure 8 SPT ‘N’ profiles
Figure 8 SPT ‘N’ profiles
Figure 9 CPT profiles
Figure 9 CPT profiles

이 결과로부터 평균 해수면 위와 아래에있는 모래 채우기의 강도와 강성의 대비를 관찰 할 수 있습니다. 이 현상은 배치 방법에 기인한다고 제안되었다 [9]. 또한 기초 수준에서 진동 압축 된 모래의 특성에도 변동이있었습니다. 분석을 위해 선택된 토양 특성은 테스트 데이터, 인근 사이트의 경험 및 유사한 토양 조건에 대한 발표 된 데이터를 기반으로합니다. 그것들은 표 1에 요약되어 있습니다. 일반적으로 시설의 건설 순서는 다음과 같습니다.

  1. Removal of pockets of soft marine clay by dredging
  2. Dredging of sand to the required level
  3. Vibro-compaction of the sand on which the counterfort units were to be founded
  4. Placement of gravel for the quay wall foundation.
  5. Placement of concrete counterfort units weighing 360 tonne each
  6. Placement of hydraulic sand fill behind the units
  7. Surcharging the fill just behind the capping beam
  8. Construct capping beam and place more sand fill to the finished level
  9. Additional surcharge prior to the operation of container cranes.

Table 1 Soil properties used in the construction
simulation of the quay wall

Table 1 Soil properties used in the construction simulation of the quay wal
Table 1 Soil properties used in the construction simulation of the quay wal

2D 평면 변형 모델의 수치 시뮬레이션에서 구성 순서 (그림 10)와 하중은 다음 단계에 따라 단순화 / 이상적입니다.

  1. The starting condition of the seabed consisted of the vibrocompacted sand, gravel bed, native sand, clay and fissured clay at depth. The “in-situ” stresses were also switched on in this step.
  2. Placement of counterfort unit (using equivalent linear elastic beam elements) with a vertical force applied through the centre of gravity of the unit to represent the buoyant self-weight.
  3. Sequentially placing hydraulic sand fill behind the unit to the level prior to surcharging.
  4. Apply an equivalent trapezoidal pressure to represent the surcharge.
  5. Placement of capping beam and the sand fill to the required level.
  6. Apply additional surcharge.
  7. Application of repeated loads from the crane seaward and landward legs.
Figure 10 Construction sequence
Figure 10 Construction sequence

분석에서는 침수 된 물질과 평균 해수면 위에있는 물질을 나타 내기 위해 적절한 밀도를 사용했습니다. 안벽의 장기적인 움직임이 중요했기 때문에 배수 된 토양 매개 변수가 사용되었습니다. 토양은 분석에서 Mohr-Coulomb 실패 기준을 따르는 것으로 가정되었습니다. 단순한 탄성-완전 소성 응력-변형 거동이 가정되었습니다. 일련의 강체 다이어그램으로 표현 된 안벽 이동의 역사는 그림 11에 나와 있습니다. 벽의 상단과 바닥에서 계산 된 수직 및 수평 이동은 그림 12와 13에 표시됩니다. 수치는 모니터링 된 데이터와 해당 상한 및 하한 (해당 상자에 표시됨)입니다. 측정에서 산란의 양에도 불구하고 벽 건설에 대해 계산 된 움직임은 합리적으로 잘 비교되었습니다. 조사 데이터와 예측을 일치시키기 위해 분석에서 토양 속성을 변경하려는 시도가 없었습니다. 반복되는 크레인 하중의 래칫 효과를 관찰 할 수 있습니다. 불행히도 반복적 인 크레인 하중 하에서 벽 이동에 대한 기준이 없었기 때문에 이러한 예상 이동을 비교할 수 없었습니다. 문제의 복잡성과 고도로 가변적 인 토양 특성을 고려할 때 계산 된 결과는 매우 고무적입니다.

Figure 11 Wall deformations
Figure 11 Wall deformations

토양에서 플라스틱 구역의 발달도 분석에서 계산되었습니다. 벽의 발가락 아래의 토양이 여러 번 과도하게 압박을받는 것으로 밝혀졌습니다. 접촉 압력은 경사 하중으로 인한 베어링 고장에 대한 안전 지표 (FOS)를 결정하는 데 사용되었습니다. 지지력은 계산 방법에 의해 크게 영향을 받았다고보고되었습니다 [10]. 원래의 기초 디자인은 덴마크 코드 [11]를 기반으로했기 때문에이 경우 일관성을 위해 사용되었습니다. 편심의 함수로서 FOS의 발전과 수평 대 수직 추력 (H / V)의 비율이 각각 그림 14와 15에 나와 있습니다.

Figure 12 Wall top movements
Figure 12 Wall top movements
Figure 13 Wall base movements
Figure 13 Wall base movements
Figure 14 ‘FOS’ vs. eccentricity
Figure 14 ‘FOS’ vs. eccentricity
Figure 15 ‘FOS’ vs. H/V ratio
Figure 15 ‘FOS’ vs. H/V ratio

그림은 벽이 추가 요금과 반복적 인 적재 단계 동안 국부적 인 베어링 고장에 가까웠음을 보여줍니다. 크레인 하중 하에서 FOS의 명백한 증가는 벽에 대한 수직 하중이 증가하는 반면 유지된 토양의 수평 압력이 다소 일정하게 유지됨에 따라 편심이 감소했기 때문입니다.

끝 맺는 말
세 가지 매우 다른 실제 응용 프로그램의 유효성 검사 프로세스가 설명되었습니다. 각 사례의 주요 특징과 결과는 표 2에 요약되어 있습니다. 재료 및 하중 불확도 및 예상 결과가 강조 표시됩니다. 건설 품질은 구조의 성능에 상당한 영향을 미치는 것으로 나타났습니다.

이는 분석가가 프로젝트의 설계 단계에서 정량화하고 정확하게 분석하지 못할 수도 있습니다. 구조가 완료된 직후 모니터링의 중요성을 간과해서는 안됩니다. 이것은 미래의 역 분석을위한 유용한 자료가 될 것입니다. 수치 도구가 이러한 복잡한 문제를 분석 할 수 있다는 사실에도 불구하고 분석가는 어떤 매개 변수가 중요하거나 중요하지 않은지 식별 할 준비가되어 있어야합니다.

익숙하지 않은 문제를 분석 할 때 유효성 검사 프로세스를 점진적으로 수행해야합니다. 아마도 검증 방법을 찾는 핵심은 수치 분석 도구를 사용하지 않고 솔루션에 도달 할 수있는 다른 방법이 있는지 묻는 것입니다. 많은 경우 이러한 솔루션은 광범위한 문헌 검색 후에 존재합니다. 그러나 다른 경우에는 실험실 테스트와 현장 관찰이 유일한 대안이 될 것입니다.

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

References
[1] Puri, S.P.S. (1998) “Avoiding Engineering Failures Caused by Computer-Related Errors”, J. Comp. in Civil Engineering, ASCE, 12(4), 170-172.
[2] Timoshenko, S.P. and Woinowsky-Krieger, S. (1959) Theory of Plates and Shells, 2nd edition, McGraw-Hill Kogakusha. p.580.
[3] BS2654 (1989) Manufacturing of vertical steel welded non-refrigerated storage tanks with butt-welded shells for the petroleum industry.
[4] Rish, R.F. (1977) “Design of Cylindrical Tanks on Elastic Foundations”, Civil Engineering Transactions, The Institution of Engineers, Australia, 192-195.
[5] US Army Corps of Engineers (1990) Hydraulic Design of Spillways, Engineer Manual No. 1110-2-1603.
[6] Hirt, C.W. and Nichols, B.D. (1981) “Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries”, J. Comp. Phys. 39, 201- 225.
[7] Ho, D.K.H., Boyes, K.M and Donohoo, S.M. (2001) “Investigation of Spillway Behaviour under Increased Maximum Flood by Computational Fluid Dynamics Technique”, Proc. Conf. 14th Australasian Fluid Mechanics, Adelaide, December, 577-580.
[8] Savage, B.M. and Johnson, M.C. (2001) “Flow over Ogee Spillway: Physical and Numerical Model Case Study”, J. Hydraulic Engineering, ASCE, 127(8), 640-649.
[9] Lee, K.M., Shen, C.K., Leung, D.H.K. and Mitchell, J.K. (1999) “Effects of placement method on geotechnical behaviour of hydraulic fill sands” J. Geotech. and Geoenviron. Engineering, ASCE, 125(10), 832-846.
[10] Sieffert, J.G. and Bay-Gress, Ch. (2000) “Comparison of European bearing capacity calculation methods for shallow foundations”, Proceedings of the Institution of Civil Engineers, Geotechnical Engineering, 143, April, 65-74.
[11] DS 415 (1984) Code of Practice for Foundation Engineering. Table 2 Summary of findings for the three case studies

2 Fluid, 2 Temperature 모델

2 Fluid, 2 Temperature 모델

우주선 및 자동차 연료 탱크 및 특정 미세 유체 장치는 안전하고 효율적인 작동을 위해 정확한 액체 및 기체 상태 모델링이 필요합니다. 이러한 시스템에 유체 계면이 존재하는 것 외에도, 열 전달 및 상 변화의 물리학도 정확하게 포착해야합니다. 얼마나 복잡합니까!

이러한 복잡한 시나리오를 시뮬레이션하기 위해 FLOW-3D v12.0에는 2 Fluid, 2 Temperature 모델이 도입되었습니다.

 

단순화 된 모델 : 2 Fluid, 1 Temperature

FLOW-3D 의 인터페이스 추적 방법인 TruVOF는 열 전달 및 위상 변화를 포함하여 2 Fluid 모델과 함께 작동합니다. 그러나,이 모델의 단순화 중 하나는, 인터페이스를 갖는 메쉬 셀의 온도가 다음의 개략도에 도시 된 바와 같이 혼합물 온도 (따라서 단순화 된 모델) Tmix로 표현된다는 것입니다.

온도가 경계면을 가로 질러 연속적이고 매끄러 울 때 혼합물 근사치가 적절하지만, 열-물리적 특성의 큰 차이로 인해 액체 및 가스가 있는 경우에는 이를 추정 할 수 없습니다. 이러한 시스템에서 용액의 정확도는 액체-기체 혼합물을 함유하는 셀에서 유체 에너지 및 온도의 평균으로부터 발생하는 과도한 수치 확산에 의해 압도 될 수 있습니다. 단순화 된 온도 슬립 모델은 이러한 경우 부분적인 솔루션만 제공합니다.

단순화 된 모델-2 Fluid, 1 Temperature

종합 모델 : 2 Fluid, 2 Temperature

1 Temperature 접근 방식의 결함을 극복하기 위해 2 Fluid 솔루션에 대한 2 Temperature 모델이 버전 11.3에 도입되었습니다. 여기에는 아래 회로도에 표시된 것처럼 각 유체에 대한 에너지 전달 방정식을 해결하고 각 상의 온도를 저장하는 작업이 포함됩니다. 자유 표면이 있는 메쉬 셀은 이제 액체 (T1)와 가스 (T2) 온도를 모두 나타냅니다.

종합 모델 : 2 유체, 2 온도

탱크 슬로싱(Tank sloshing)

탱크 슬로싱에 대한 이 사례 연구에서, 액체는 초기 온도 300K이고 가스는 400K입니다. 단순화 된 모델과 포괄적인 모델 사이의 수치 확산 정도의 차이는 아래 애니메이션에 나와 있습니다. 온도 윤곽에서 시간이 지남에 따라 용액의 수치 확산은 1 Temperature 접근 방식으로 보여지고 계면 물리를 완전히 가리게 됩니다.

단순화 된 모델 : 2 Fluid, 1 Temperature

종합 모델 : 2 Fluid, 2 Temperature

공기중 드롭 용접(Drop welding in air)

이 낙하 용접 사례 연구에서 액체 금속은 중력 하에서 2300K에서 공기를 통해 고체화 된 금속 베드로 떨어집니다. 공기 및 베드 초기 온도는 293K입니다. simplified model에서는 수치 확산으로 인해 액체 금속 낙하 온도가 베드에 도달하기 전에도 급격히 감소하기 시작합니다. 반면에 comprehensive model에서는 방울이 초기 온도를 유지하여 훨씬 더 나은 솔루션을 제공합니다.

단순화 된 모델을 사용한 온도 필드 진화

종합 모델의 온도 필드

FLOW-3D의 2 Fluid, 2 Temperature 모델과 유체 인터페이스 추적을 결합하면 사용자는 특히 연료 슬로싱 시스템과 같이 복잡한 열전달 및 위상 변화 문제를 정확하게 모델링 할 수 있습니다.

이 새로운 모델에 대한 제안이나 의견은 adwaith@flow3d.com에 문의하십시오.

FLOW-3D 용접해석 개요

FLOW-3D 용접해석 개요

자료 제공: FLOW Science Japan

용접은 금속의 상변화, 용융시의 유동, 방열에 의한 응고 등을 포함한 복잡한 물리 현상입니다. FLOW-3D@에는 그 현상을 정밀하게 모델링하는 기능이 있고, 용접 현상을 충실하게 재현할 수 있습니다.

특히 용융 금속의 유동은 표면 장력의 영향이 강하고, 그 해석은 정확한 자유 계면의 추적이 필요합니다.

FLOW-3D@의 정확한 계면 추적 기능인 TruVOF®는 그 현상을 파악하기에 매우 적합합니다.

해석 조건으로 레이저의 power, spot size, distance, 움직임(분사방향) 등을 입력할 수 있으며, 보호가스는 밀도와 유속으로 설정 할 수 있습니다. 스패터의 분출 해석은 현재 지원되지 않습니다.

해석 결과로 용접 비드의 폭, 깊이, 기공의 유무 등 관찰할 수 있습니다.

열응력 해석은 FLOW-3D 유동 해석 결과를 Abaqus 등의 구조해석 프로그램에서 불러와서 별도로 열응력 해석을 수행해야 합니다.

해석 필요성

FLOW-3D 를 이용한 용접해석은

  • 높은 방사 강도와 고온으로 직접 관찰하기 어려운 내부 현상의 상세 내용을 가시화
  • 온도, 열, 용접 속도, 위치 관계, 재료 물성 등의 파라미터 연구 검토
  • 결함 예측 (공기연행, 응고 수축, 금속 산화) 등의 필요성

해석을 통해 얻는 이점

금속의 상변화, 용융시의 유체의 힘, 방열에 의한 응고 등의 물리모델 용접 현상을 분석할 수 있습니다.

또한 용융시에는 표면장력의 영향이 강하고, 자유계면을 추적하는 수치해석 방법에 대해 높은 정밀도가 요구됩니다. FLOW-3D@는 이러한 요구사항을 잘 처리할 수 있는 장점을 가지고 있어, 용접 용융-응고의 연속적인 현상을 정확하게 파악합니다.

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

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

레이저 Soldering의 용융지 형성에 미치는 유체 대류의 영향

레이저 Soldering의 용융지 형성에 미치는 유체 대류의 영향

자료 제공: 오하이오 주립대학교
자료 제공: FLOW Science Japan

Laser Soldering
Laser Soldering

레이저 Soldering에서는 금속 분말은 레이저 빔 주위에 고리 모양으로 배치된 여러 종류의 분말 공급노즐을 통해 불활성 캐리어 가스 중으로 분사됩니다. 용융지의 형상은 표면장력으로 유도되는 말랑고니 대류 패턴의 영향을 크게 받는 것으로 알려져 있습니다. 따라서 용융지 내의 유체 흐름을 정밀하게 예측하기 위해서 FOW-3D@ WELD를 이용한 레이저 Soldering 프로세스를 개략화하여 시뮬레이션 할 수 있습니다.

해석과 시뮬레이션 비교1
해석과 시뮬레이션 비교1
해석과 시뮬레이션 비교2

다른 레이저 출력의 예측된 용접폭, 높이 및 용융지의 깊이는 실험 결과 실험 측정값과 동등한 결과를 얻을 수 있는 것을 확인할 수 있습니다. 용융지의 열유동에 의해 가장 깊은 용융 영역은 마랑고니 대류에 의해 유도된 2개의 마주보는 표면 흐름의 충돌로 인해 형성되는 것을 확인할 수 있습니다.

열전달(Heat Transfer)

열전달(Heat Transfer)

열전달은 전도, 대류 및 복사를 통한 열 에너지의 전달입니다. 일반적이지만 매우 중요한 물리적 현상입니다. 재료 특성 및 기타 물리적 현상은 온도 (또는 열에너지)에 매우 민감합니다. FLOW-3DFLOW-3D  CAST의 열전달 모델은 전도, 대류 및 복사를 통해 유체 내, 고체 및 공극 내에서 열전달을 처리하는 완전 복합 열전달 방정식을 해석합니다.

또한, 이 모델은 사용자가 다양한 애플리케이션을 모델링 할 수 있도록 유연하고 편리한 옵션을 제공합니다.

  • 명시적 및 암시적인 열전달 옵션을 모두 사용할 수 있습니다. 암시적 방법을 사용하여 명시적 접근과 관련된 시간 단계별 크기의 안정성 제한을 제거 할 수 있습니다. 전도성 또는 열전달의 안정성 제한이 시뮬레이션에서 다른 안정성 제한보다 실질적으로 작을 때, 암시적 방법을 사용하면 계산 효율성이 크게 향상 될 수 있습니다.

  • 각기 다른 매체 사이의 열전달 계수는 흐름 유형에 따라 사용자 정의되거나 자동으로 계산 될 수 있습니다.


  • 1차 및 2차 열에너지의 이류 알고리즘을 모두 사용할 수 있습니다. 1차 옵션은 효율적이고 견고하며 대부분의 열전달 문제에 적합하지만 높은 열 구배가 예상되는 시뮬레이션의 경우 인공적인 열 확산으로 이어질 수 있습니다. 2차 옵션은 가령, 부력 중심의 흐름에서 온도 구배를 해결하는 것이 중요한 상황에 적합합니다.


  • 유체와 고체 사이의 열전달을 모델링하기 위해 여러 가지 옵션을 사용할 수 있습니다 (지정된 열유속에서 전원, 규정 온도까지). 이러한 옵션은 다양한 프로세스 및 응용 프로그램을 모델링 할 수 있는 유연성과 성능을 제공합니다.

다른 물리 모델과 함께 FLOW-3DFLOW-3D  CAST의 열전달 모델은 고급 모델링 기능을 위한 견고한 토대가 됩니다. 예를 들어, 액체 / 고체 및 액체 / 증기 상 변화 모델을 사용하여 금속 응고, 물의 건조 및 비등, 분무 냉각을 시뮬레이션 할 수 있습니다. 점성 가열은 고속 점성 흐름에도 포함될 수 있습니다.

Advances in Magnetohydrodynamic Liquid Metal Jet Printing

Advances in Magnetohydrodynamic Liquid Metal Jet Printing

Scott Vader1, Zachary Vader1, Ioannis H. Karampelas2 and Edward P. Furlani2, 3
1Vader Systems, Buffalo, NY
2Dept. of Chemical and Biological Engineering, 3 Dept. of Electrical Engineering,
University at Buffalo SUNY, NY 14260, Office: (716) 645-1194, Fax: (716) 645-3822, efurlani@buffalo.edu

ABSTRACT

자기유체역학적 액체 금속 제트 프린팅

우리는 용해된 금속 방울을 3D 물체로 만드는 새로운 주문형 DOD(Drop-on-Demand) 인쇄 방법을 제안합니다. 이 접근 방식에서는 단단한 금속 와이어가 인쇄 헤드 내에서 용해된 다음 펄스 자기장에 노출됩니다.

적용된 필드가 챔버에 침투하여 액상 금속 내에 자기 유압(MHD) 기반 압력 펄스를 유도하여 금속 일부가 노즐 챔버를 통해 이동된 후 배출됩니다. 표면 장력은 분출된 금속 위에 작용하여 가해진 압력에 따라 초 당 수 미터 범위의 속도로 구형 방울을 형성합니다.

잠시 비행한 후 방울이 기질에 충돌하여 냉각되어 고체 덩어리를 형성합니다. 따라서 패턴이 있는 증착 및 드롭 방식의 고형화를 통해 3D 솔리드 구조를 인쇄할 수 있습니다.

현재 연구에서는 샘플 프린팅 구조와 함께 시제품 MHD 프린팅 시스템 개발에 대한 발전된 점을 제시합니다. 또한 드롭 생성을 관리하는 기본 물리학에 대해 논의하고 장치 성능을 예측하기 위한 새로운 컴퓨팅 모델을 소개합니다.

Computational model of magnetohydrodynamic-based drop generation
Computational model of magnetohydrodynamic-based drop generation (printhead reservoir and ejection chamber
not shown): (a) the magnetic field generated by a pulsed coil is shown

INTRODUCTION

주문형 드롭온 잉크젯 프린팅은 상업 및 소비자 이미지 재현을 위한 잘 확립된 방법입니다. 이 기술을 추진하는 원리와 동일한 원리가 기능 인쇄 및 적층 제조 분야에도 적용될 수 있습니다.

Early stage prototype of a single nozzle printhead
Early stage prototype of a single nozzle printhead

기존의 잉크젯 기술은 폴리머에서 살아있는 세포에 이르는 다양한 재료를 증착하고 패터링하여 다양한 기능성 매체, 조직 및 장치를 프린팅하는 데 사용되어 왔습니다. 현재 진행 중인 작업을 통해 잉크젯 인쇄를 3D 금속 부품으로 확장하려고 시도하고 있습니다.

현재, 대부분의 3D 금속 인쇄 애플리케이션은 고체 물체를 형성하기 위해 레이저(예: 선택적 레이저 소거 [1] 및 직접 금속 소거[2]) 또는 전자 빔(예: 전자 빔 용해 [3])과 같은 외부 유도 에너지원에 의해 소거 또는 녹는 퇴적 금속 분말을 포함합니다.

그러나 이러한 방법은 비용과 복잡성, 즉 3D 프린팅 공정에 앞서 금속을 분쇄해야 한다는 점에서 일정한 단점이 있을 수 있습니다.

이 프레젠테이션에서는 자기 유압 역학 원리를 기반으로 하는 금속 적층 제조의 근본적으로 다른 접근 방식을 제안합니다. 이 방법은 스풀링된 고체 금속 와이어를 인쇄 헤드에 공급하고 노즐에서 업스트림을 예열하여 노즐 챔버에 공급되는 액체 금속 저장소를 형성하는 것입니다. 챔버가 채워지면 액체 금속 내에서 과도 전류를 유도하는 펄스 자기장이 인가됩니다. 유도 전류가 인가된 필드에 결합되어 로렌츠 힘 밀도를 생성하여, 인가된 압력에 따라 속도가 달라지는 용융 금속 방울을 배출하는 작용을 하는 챔버 내의 유사 압력을 제공합니다.

방울은 냉각된 기질에 투영되어 고체 덩어리를 형성합니다. 3D 솔리드 구조를 패터닝으로 인쇄할 수 있습니다. 방울의 침적과 방울의 현명한 응고입니다. 이 유망한 신기술은 낮은 재료 비용, 높은 제조율 및 매력적인 재료 특성 때문에 적층 제조 애플리케이션에 광범위한 영향을 미칠 수 있습니다.

현재 작업에서는 새로운 3D 인쇄 시스템을 도입하고 기기 개발의 진보를 설명하고 샘플 인쇄 구조를 시연합니다. 또한 드롭 생성-배출 메커니즘에 대해 설명하고 인쇄 성능을 예측하기 위한 일련의 새로운 컴퓨팅 모델을 제시합니다.

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

[FLOW-3D 이론] 1. 개요

  1. 개요

FLOW-3D는 범용 전산 유체 역학(CFD) 소프트웨어입니다. 유체의 운동 방정식을 계산하기 위해 특별히 개발된 수치 기법을 사용하여 다중 스케일, 다중 물리 흐름 문제에 대해 과도적 3차원 해결책을 얻습니다. 다양한 물리적 및 수치 옵션을 통해 사용자는 다양한 유체 흐름 및 열 전달 현상 분석을 위해 FLOW-3D를 적용할 수 있습니다.

유체 운동은 비선형, 과도, 2차 미분 방정식으로 설명됩니다. 이러한 방정식을 풀기 위해 유체 운동 방정식을 사용해야합니다. 이러한 방법을 개발하는 과학을 전산 유체 역학이라고 합니다. 이 방정식의 수치해는 대수적 표현으로 다양한 항을 근사화 합니다. 그런 다음 결과 방정식을 해결하여 원래 문제에 대한 대략적인 해결책을 제시합니다. 이 과정을 시뮬레이션이라고 합니다. FLOW-3D에서 사용할 수 있는 수치해석 알고리즘의 개요는 운동 방정식에 대한 섹션에 나옵니다.

일반적으로 수치 모델은 계산 Mesh 또는 그리드로 시작합니다. 이것은 여러 개의 서로 연결된 요소 또는 셀로 구성됩니다. 이러한 셀은 물리적 공간을 해당 볼륨과 관련된 여러 노드가 있는 작은 볼륨으로 세분화합니다. 노드는 압력, 온도 및 속도와 같은 미지수의 값을 저장하는데 사용됩니다. Mesh는 사실상 원래의 물리적 공간을 대체하는 숫자 공간입니다. 또한 별도의 위치에서 흐름 파라미터를 정의하고, 경계 조건을 설정하고, 유체 운동 방정식의 수치 근사치를 개발하는 방법을 제공합니다. FLOW-3D 접근 방식은 흐름 영역을 직사각형 셀의 격자로 세분하는 것입니다. 이 격자는 brick elements라고도 합니다.

계산 Mesh는 물리적 공간을 효과적으로 이산화 시킵니다. 각 유체 매개 변수는 불연속 지점에서 값 배열에 의해 Mesh로 표시됩니다. 실제 물리적 파라미터는 공간에서 연속적으로 변하기 때문에 노드 사이의 간격이 미세한 Mesh는 더 거친 Mesh보다 현실을 더욱 잘 표현해줍니다. 그런 다음 수치 근사치의 기본 속성에 도달합니다. 그리드 간격이 줄어들면 유효한 모든 유효한 수치 근사가 원래 방정식에 접근합니다. 근사치가 이 조건을 만족하지 않으면 올바르지 않은 것으로 간주해야 합니다.

동일한 물리적 공간에 대해 격자 간격을 줄이거나 Mesh를 조정하면 더 많은 요소와 노드가 생겨 수치 모델의 크기가 커집니다. 그러나 유체 흐름 및 열 전달의 실제 현실과는 별도로, 시뮬레이션 엔지니어들이 적절한 크기의 Mesh를 선택하도록 하는것과 밀접한 관계에 있는 설계 주기, 컴퓨터 하드웨어 및 마감일의 현실적인 문제도 있습니다. 이러한 제약 조건을 만족시키는 것과 사용자가 정확한 결과를 얻는 것 사이에서 타협점을 찾는 것은, CFD 모델 개발 못지않은 중요한 균형 잡힌 행위입니다.

직사각형 그리드는 규칙적이거나 구조적인 특성 때문에 생성 및 저장이 매우 쉽습니다. 균일하지 않은 그리드 간격은 복잡한 흐름 도메인을 매칭할 때 유연성을 더합니다. 연산 셀은 세 개의 지수를 사용하여 연속적으로 번호가 매겨집니다. 즉, x 방향은 i, y 방향은 j, z 방향은 k입니다. 이 방법으로 3차원 Mesh의 각 셀은 물리적 공간의 점의 좌표와 유사한 고유한 주소(i, j, k)로 식별할 수 있습니다.

구조화된 직사각형 그리드는 수치적 방법의 개발의 상대적 용이성, 원래의 물리적 문제와의 관계에 대한 후자의 투명성, 그리고 마지막으로 수치적 해결의 정확성과 안정성의 추가적인 이점을 가지고 있습니다. 유한 차분법과 유한 체적법에 기초한 가장 오래된 수치 알고리즘은 원래 이러한 Mesh에서 개발되었습니다. 이것은 FLOW-3D에서 수치적 접근방식의 핵심을 형성합니다. 유한차분법은 테일러 확장의 특성과 파생된 정의의 직접적인 적용에 기초합니다. 미분 방정식에 대한 수치적 해결책을 얻기 위해 적용된 방법 중 가장 오래된 방법이며, 첫 번째 적용은 1768년 오일러에 의해 개발된 것으로 간주됩니다. 유한체적법은 유체 운동을 위한 보존법의 일체형태에서 직접 파생되므로 자연적으로 보존 특성을 보유합니다.

FLOW-3D는 일반적인 유체 방정식의 다른 제한 사례에 해당하는 여러 모드에서 작동할 수 있습니다. 예를 들어, 하나의 모드는 압축 가능한 흐름을 위한 것이고 다른 하나는 압축할 수 없는 흐름 상황을 위한 것입니다. 후자의 경우 유체의 밀도와 에너지가 일정하다고 가정할 수 있으므로 계산할 필요가 없습니다. 또한 1유체 모드와 2유체 모드가 있습니다. 자유 표면은 단일 유체 비압축 모드에 포함될 수 있습니다. 이러한 작동 모드는 동작 방정식에 대한 다양한 선택에 해당합니다.

자유 표면은 FLOW-3D로 수행된 많은 시뮬레이션에서 존재합니다. 유량 매개변수와 재료 특성(밀도, 속도, 압력 등)이 불연속성을 경험하기 때문에 모든 계산 환경에서 자유 표면을 모델링하는 것은 어렵습니다. FLOW-3D에서는, 액체에 인접한 가스의 관성이 무시되고, 가스에 의해 점유되는 부피는 균일한 압력과 온도로만 표현되는 빈 공간, 질량의 공백으로 대체됩니다. 대부분의 경우 가스 모션의 세부 사항은 훨씬 무거운 액체의 움직임에 중요하지 않기 때문에 이 접근 방식은 계산 노력을 줄이는 이점이 있습니다. 자유 표면은 액체의 외부 경계 중 하나가 됩니다. 자유 표면의 경계 조건에 대한 적절한 정의는 자유 표면 역학을 정확하게 포착하기 위해 중요합니다.

VOF(Volume of Fluid) 방법은 이러한 목적으로 FLOW-3D에 사용됩니다. 유체 함수의 볼륨 정의, VOF 전송 방정식 해결 방법, 자유 표면의 경계 조건 설정 등 세 가지 주요 구성요소로 구성됩니다.

일부 물리 및 수치 모델은 Flow Science의 기술 노트: http://users.flow3d.com/technical-notes/ 에 자세히 설명되어 있으며, 여기에는 예제도 포함되어 있습니다.

Cavitation | 캐비테이션

캐비테이션이란 무엇입니까?

The spillways of the Glen Canyon dam in 1983 (Lee and Hoopes, 1996).

캐비테이션은 유체 흐름의 매우 낮은 압력 또는 포화 압력을 높이는 온도 상승으로 인해 유체 내에서 증기 또는 기포가 빠르게 발생하는 것입니다. 기포의 갑작스런 출현 (및 후속 붕괴)은 비압축성 유체 내에서 압력의 급격한 변화를 일으켜 심각한 기계적 손상을 일으킬 수 있습니다. 캐비테이션에 의해 유도 된 힘은 1983 년 Glen Canyon 댐의 배수로에서 경험 한 손상에서 볼 수 있듯이 며칠 내에 수 피트의 암석을 침식 할 가능성이 있습니다 (Lee and Hoopes, 1996).

또한 고압 다이 캐스팅에서 캐비테이션이 발생할 수 있습니다. 다이의 수축 및 곡선을 통한 용융 합금의 빠른 이동은 급속한 압력 강하를 초래하고 후속 캐비테이션으로 이어질 수 있습니다. 생성된 증기 기포는 최종 주조에서 다공성을 유발하거나 더 나쁜 경우 다이에 손상을 일으켜 주조품을 훼손시키고 다이 수명을 감소시킬 수 있습니다.

캐비테이션은 터빈과 파이프에 손상을 줄 수 있고, 댐의 배수로에서 콘크리트를 침식하는 등의 원인이 될 수 있습니다. 아래 이미지는 댐의 배수로 바닥 근처의 콘크리트 침식을 보여줍니다. 댐에 사용되는 콘크리트는 일반적으로 강도가 높지만 캐비테이션은 여전히 그것을 부식시킬 수 있습니다.

Eroded concrete due to cavitation on the spillway of a dam

캐비테이션은 때때로 오염 물질과 유기 분자를 분해하고, 소수성 화학 물질을 결합하고, 캐비테이션 기포의 파열로 인해 생성 된 충격파를 통해 신장 결석을 파괴하고, 혼합을위한 난류를 증가시켜 수질 정화와 같은 특정 산업 응용 분야에서 의도적으로 유도됩니다.

따라서 캐비테이션이 발생할 가능성이있는 위치와 그 강도를 이해하는 것이 중요합니다. 캐비테이션을 실험을 수행하거나 실험 결과의 현상을 시각화하는 것이 어렵고, 잠재적으로 손상 될 수 있으므로 수치해석 시뮬레이션으로 검토하는 것이 매우 필요하고, 유용합니다.

Real-World Applications | 실제 응용 분야

  • 물 및 환경 구조 내에서 손상을 주는 캐비테이션 시뮬레이션
  • 다이 손상 및 주조 다공성을 유발할 수 있는 고압 다이 캐스팅 중 캐비테이션 시뮬레이션
  • MEMS 장치 내의 열 거품 형성 시뮬레이션
  • 열 전달 표면의 비등 거동 예측
  • 캐비테이션 역학으로 인한 혼합 예측

Modeling Cavitation in FLOW-3D

FLOW-3D의 캐비테이션 모델은 thermal bubble jets 와 MEMS devices를 시뮬레이션하는데 성공적으로 사용되었습니다. FLOW-3D는 “active”또는 “passive” 모델 옵션을 제공합니다. Active 모델은 기포 영역을 열고 수동 모델은 흐름을 통해 캐비테이션 기포의 존재를 추적하고 전파하지만, 기포 영역의 형성을 시작하지는 않습니다.

Active모델은 더 큰 캐비테이션 영역이 예상되고 유동장에 영향을 미치는 경우에 가장 적합하며, Passive모델은 작은 기포의 간단한 모양이 예상되는 시뮬레이션에 가장 적합합니다. 활성 모델과 에너지 전송 계산을 통해 위상 변화도 옵션입니다. 기포는 계면에서의 증발 또는 응축으로 인해 추가로 팽창하거나 수축 할 수 있습니다.

Sample Results

아래 시뮬레이션은 수축 노즐을 보여줍니다. 애니메이션은 매우 일시적인 진동 동작을 보여주는 캐비테이션 버블의 진화를 보여줍니다. 캐비테이션 부피 분율은 초기 연속 액체에서 캐비테이션의 시작을 시각화하기 위해 플롯됩니다.

아래 애니메이션은 진입 속도가 8m/s이고 수렴 기울기가 18 °이고 발산 기울기가 8 ° 인 벤츄리 내의 캐비테이션을 보여줍니다. 다시 말하지만, 캐비테이션의 과도 동작은 잘 모델링되어 있으며, 모델은 22ms의 실험 결과와 비교하여 17.4ms의 캐비테이션주기 기간을 예측합니다 (Stutz and Reboud 1997).

Cavitation in a venturi

물 탱크를 통해 이동하는 고속 발사체를 시뮬레이션하여 발사체 후류에서 생성 된 저압 영역의 공동 기둥을 보여줍니다. 발사체의 초기 속도는 600m / s입니다. 아래는 탱크의 움직임과 후행하는 캐비테이션 유체의 애니메이션입니다. 발사체가 감속함에 따라 캐비테이션 기둥의 반경이 좁아집니다.

@

High-speed bullet

References

Lee, W., Hoopes, J.A., 1996, Prediction of Cavitation Damage for Spillways, Journal of Hydraulic Engineering, 122(9): 481-488.

Plesset, M.S., Prosperetti, A., 1977, Bubble Dynamics and Cavitation, Annual Revue of Fluid Mech, 9: 145-185.

Rouse, H., 1946. Elementary Mechanics of Fluids, New York: Dover Publications, Inc.

Stutz, B., Reboud, J.L., 1997, Experiments on unsteady cavitation, Experiments in Fluids, 22: 191-198.

The realm of operations of FLOW-3D

ADDITIVE MANUFACTURING SIMULATIONS

Capabilities of FLOW-3D

FLOW-3D는 자유 표면 유체 흐름 시뮬레이션을 전문으로하는 다중 물리 CFD 소프트웨어입니다. 자유 표면의 동적 진화를 추적하는 소프트웨어의 알고리즘인 VOF (Volume of Fluid) 방법은 Flow Science의 설립자인 Tony Hirt 박사가 개척했습니다.

또한 FLOW-3D에는 금속 주조, 잉크젯 인쇄, 레이저 용접 및 적층 제조 (AM)와 같은 광범위한 응용 분야를 시뮬레이션하기위한 물리 모델이 내장되어 있습니다.
적층 제조 시뮬레이션 소프트웨어, 특히 L-PBF (레이저 파우더 베드 융합 공정)의 현상 유지는 열 왜곡, 잔류 응력 및지지 구조 생성과 같은 부분 규모 모델링에 도움이되는 열 기계 시뮬레이션에 초점을 맞추고 있습니다.

유용하지만 용융 풀 역학 및 볼링 및 다공성과 같은 관련 결함에 대한 정보는 일반적으로 이러한 접근 방식의 영역 밖에 있습니다. 용융 풀 내의 유체 흐름, 열 전달 및 표면 장력이 열 구배 및 냉각 속도에 영향을 미치며 이는 다시 미세 구조 진화에 영향을 미친다는 점을 명심하는 것도 중요합니다.

FLOW-3D와 이산 요소법 (DEM) 및 WELD 모듈을 사용하여 분말 및 용융 풀 규모에서 시뮬레이션 할 수 있습니다.
구현되는 관련 물리학에는 점성 흐름, 열 전달, 응고, 상 변화, 반동 압력, 차폐 가스 압력, 표면 장력, 움직이는 물체 및 분말 / 입자 역학이 포함됩니다. 이러한 접근 방식은 합금에 대한 공정을 성공적으로 개발할 수 있게 하고, AM 기계 제조업체와 AM 기술의 최종 사용자 모두에게 관심있는 미세 구조 진화에 대한 통찰력을 제공하는데 도움이 됩니다.

The realm of operations of FLOW-3D
The realm of operations of FLOW-3D

FLOW-3D는 레이저 분말 베드 융합 (L-PBF), 직접 에너지 증착 (DED) 및 바인더 제트 공정으로 확장되는 기능을 가지고 있습니다.
FLOW-3D를 사용하면 분말 확산 및 패킹, 레이저 / 입자 상호 작용, 용융 풀 역학, 표면 형태 및 후속 미세 구조 진화를 정확하게 시뮬레이션 할 수 있습니다. 이러한 기능은 FLOW-3D에 고유하며 계산 효율성이 높은 방식으로 달성됩니다.

예를 들어 1.0mm x 0.4mm x 0.3mm 크기의 계산 영역에서 레이저 빔의 단일 트랙을 시뮬레이션하기 위해 레이저 용융 모델은 단 8 개의 물리적 코어에서 약 2 시간이 걸립니다.
FLOW-3D는 모든 관련 물리 구현 간의 격차를 해소하는 동시에 업계 및 연구 표준에서 허용하는 시간 프레임으로 결과를 생성합니다. 분말 패킹, 롤러를 통한 파워 확산, 분말의 레이저 용융, 용융 풀 형성 및 응고를 고려하고 다층 분말 베드 융합 공정을 위해 이러한 단계를 순차적으로 반복하여 FLOW-3D에서 전체 AM 공정을 시뮬레이션 할 수 있습니다.

FLOW-3D의 다층 시뮬레이션은 이전에 응고된 층의 열 이력을 저장한다는 점에서 독특하며, 열 전달을 고려하여 이전에 응고된 층에 확산된 새로운 분말 입자 세트에 대해 시뮬레이션이 수행됩니다.
또한, 응고 된 베드의 열 왜곡 및 잔류 응력은 FLOW-3D를 사용하여 평가할 수 있으며, 보다 복잡한 분석을 수행하기 위해 FLOW-3D의 압력 및 온도 데이터를 Abaqus 및 MSC Nastran과 같은 FEA 소프트웨어로 내보낼 수 있습니다.

Sequence of a multi-layer L-PBF simulation setup in FLOW-3D

Ease of Use

FLOW-3D는 다양한 응용 분야에서 거의 40 년 동안 사용되어 왔습니다. 사용자 피드백을 기반으로 UI 개발자는 소프트웨어를 사용하기 매우 직관적으로 만들었으며 새로운 사용자는 시뮬레이션 설정의 순서를 거의 또는 전혀 어려움없이 이해합니다.
사용자는 FLOW3D에서 구현 된 다양한 모델의 이론에 정통하며 새로운 실험을 설계 할 수 있습니다. 실습 튜토리얼, 비디오 강의, 예제 시뮬레이션 및 기술 노트의 저장소도 사용할 수 있습니다.
사용자가 특정 수준의 경험에 도달하면 고급 수치 교육 및 소프트웨어 사용자 지정 교육을 사용할 수 있습니다.

Available Literature

실험 데이터에 대해 FLOW-3D 모델을 검증하는 몇 가지 독립적으로 발표된 연구가 있습니다. 여기에서 수록된 저널 논문은 레이저 용접 및 적층 제조 공정으로 제한됩니다. 더 많은 참조는 당사 웹 사이트에서 확인할 수 있습니다.

Laser Welding

  1. L.J.Zhang, J.X.Zhang, A.Gumenyuk, M.Rethmeier, S.J.Na, Numerical simulation of full penetration laser welding of thick steel plate with high power high brightness laser, Journal of Materials Processing Technology, Volume 214, Issue 8, 2014.
    A study by researchers from BAM in Germany, KAIST in Korea, and State Key Laboratory of Mechanical Behavior of Materials in China that focuses on keyhole dynamics and full penetration laser welding of steel plates.
  2. Runqi Lin, Hui-ping Wang, Fenggui Lu, Joshua Solomon, Blair E.
    Carlson, Numerical study of keyhole dynamics and keyhole-induced porosity formation in remote laser welding of Al alloys, International Journal of Heat and Mass Transfer, Volume 108, Part A, 2017.
    General Motors (GM) and Shangai University collaborated on a study on the influence of welding speed and weld angle of inclination on porosity occurrence in laser keyhole welding.
  3. Koji Tsukimoto, Masashi Kitamura, Shuji Tanigawa, Sachio Shimohata, and Masahiko Mega, Laser Welding Repair for Single Crystal Blades, International Gas Turbine Congress, Tokyo, 2015.
    Mitsubishi Heavy Industry’s study on laser welding repair using laser cladding for single Ni crystal alloys used in gas turbine blades.

Additive Manufacturing

  1. Yu-Che Wu, Cheng-Hung San, Chih-Hsiang Chang, Huey-Jiuan Lin, Raed Marwan, Shuhei Baba, Weng-Sing Hwang, Numerical modeling of melt-pool behavior in selective laser melting with random powder distribution and experimental validation, Journal of Materials Processing Technology, Volume 254, 2018
    This paper discusses powder bed compaction with random packing for different powder-size distributions, and the importance of considering evaporation effects in the melting process to validate the melt pool dimensions.
  2. Lee, Y.S., and W.Zhang, Mesoscopic simulation of heat transfer and fluid flow in laser powder bed additive manufacturing, Proceedings of the Annual International Solid Freeform Fabrication Symposium, Austin, TX, USA. 2015
    A study conducted by Ohio State University researchers to understand the influence of process parameters in formation of balling defects.
  3. Y.S. Lee, W.Zhang, Modeling of heat transfer, fluid flow and solidification microstructure of nickel-base superalloy fabricated by laser powder bed fusion, Additive Manufacturing, Volume 12, Part B, 2016
    A study conducted by Ohio State University researchers to understand the influence of solidification parameters, calculated from the temperature fields, on solidification morphology and grain size using existing theoretical models in laser powder bed fusion processes.

 

 

자유 표면 모델링 방법

본 자료는 국내 사용자들의 편의를 위해 원문 번역을 해서 제공하기 때문에 일부 오역이 있을 수 있어서 원문과 함께 수록합니다. 자료를 이용하실 때 참고하시기 바랍니다.

Free Surface Modeling Methods

An interface between a gas and liquid is often referred to as a free surface. The reason for the “free” designation arises from the large difference in the densities of the gas and liquid (e.g., the ratio of density for water to air is 1000). A low gas density means that its inertia can generally be ignored compared to that of the liquid. In this sense the liquid moves independently, or freely, with respect to the gas. The only influence of the gas is the pressure it exerts on the liquid surface. In other words, the gas-liquid surface is not constrained, but free.

자유 표면 모델링 방법

기체와 액체 사이의 계면은 종종 자유 표면이라고합니다.  ‘자유’라는 호칭이 된 것은 기체와 액체의 밀도가 크게 다르기 때문입니다 (예를 들어, 물 공기에 대한 밀도 비는 1000입니다).  기체의 밀도가 낮다는 것은 액체의 관성에 비해 기체의 관성은 일반적으로 무시할 수 있다는 것을 의미합니다.  이러한 의미에서, 액체는 기체에 대해 독립적으로, 즉 자유롭게 움직입니다.  기체의 유일한 효과는 액체의 표면에 대한 압력입니다.  즉, 기체와 액체의 표면은 제약되어있는 것이 아니라 자유롭다는 것입니다.

In heat-transfer texts the term ‘Stephen Problem’ is often used to describe free boundary problems. In this case, however, the boundaries are phase boundaries, e.g., the boundary between ice and water that changes in response to the heat supplied from convective fluid currents.

열전달에 관한 문서는 자유 경계 문제를 묘사할 때 “Stephen Problem’”라는 용어가 자주 사용됩니다.  그러나 여기에서 경계는 상(phase) 경계, 즉 대류적인 유체의 흐름에 의해 공급된 열에 반응하여 변화하는 얼음과 물 사이의 경계 등을 말합니다.

Whatever the name, it should be obvious that the presence of a free or moving boundary introduces serious complications for any type of analysis. For all but the simplest of problems, it is necessary to resort to numerical solutions. Even then, free surfaces require the introduction of special methods to define their location, their movement, and their influence on a flow.

이름이 무엇이든, 자유 또는 이동 경계가 존재한다는 것은 어떤 유형의 분석에도 복잡한 문제를 야기한다는 것은 분명합니다. 가장 간단한 문제를 제외한 모든 문제에 대해서는 수치 해석에 의존할 필요가 있습니다. 그 경우에도 자유 표면은 위치, 이동 및 흐름에 미치는 영향을 정의하기 위한 특별한 방법이 필요합니다.

In the following discussion we will briefly review the types of numerical approaches that have been used to model free surfaces, indicating the advantages and disadvantages of each method. Regardless of the method employed, there are three essential features needed to properly model free surfaces:

  1. A scheme is needed to describe the shape and location of a surface,
  2. An algorithm is required to evolve the shape and location with time, and
  3. Free-surface boundary conditions must be applied at the surface.

다음 설명에서는 자유 표면 모델링에 사용되어 온 다양한 유형의 수치적 접근에 대해 간략하게 검토하고 각 방법의 장단점을 설명합니다. 어떤 방법을 사용하는지에 관계없이 자유롭게 표면을 적절히 모델화하는 다음의 3 가지 기능이 필요합니다.

  1. 표면의 형상과 위치를 설명하는 방식
  2. 시간에 따라 모양과 위치를 업데이트 하는 알고리즘
  3. 표면에 적용할 자유 표면 경계 조건

Lagrangian Grid Methods

Conceptually, the simplest means of defining and tracking a free surface is to construct a Lagrangian grid that is imbedded in and moves with the fluid. Many finite-element methods use this approach. Because the grid and fluid move together, the grid automatically tracks free surfaces.

라그랑주 격자 법

개념적으로 자유 표면을 정의하고 추적하는 가장 간단한 방법은 유체와 함께 이동하는 라그랑주 격자를 구성하는 것입니다. 많은 유한 요소 방법이 이 접근 방식을 사용합니다. 격자와 유체가 함께 움직이기 때문에 격자는 자동으로 자유 표면을 추적합니다.

At a surface it is necessary to modify the approximating equations to include the proper boundary conditions and to account for the fact that fluid exists only on one side of the boundary. If this is not done, asymmetries develop that eventually destroy the accuracy of a simulation.

표면에서 적절한 경계 조건을 포함하고 유체가 경계의 한면에만 존재한다는 사실을 설명하기 위해 근사 방정식을 수정해야합니다. 이것이 수행되지 않으면 결국 시뮬레이션의 정확도를 훼손하는 비대칭이 발생합니다.

The principal limitation of Lagrangian methods is that they cannot track surfaces that break apart or intersect. Even large amplitude surface motions can be difficult to track without introducing regridding techniques such as the Arbitrary-Lagrangian-Eulerian (ALE) method. References 1970 and 1974 may be consulted for early examples of these approaches.

라그랑지안 방법의 주요 제한은 분리되거나 교차하는 표면을 추적 할 수 없다는 것입니다. ALE (Arbitrary-Lagrangian-Eulerian) 방법과 같은 격자 재생성 기법을 도입하지 않으면 진폭이 큰 표면 움직임도 추적하기 어려울 수 있습니다. 이러한 접근법의 초기 예를 보려면 참고 문헌 1970 및 1974를 참조하십시오.

The remaining free-surface methods discussed here use a fixed, Eulerian grid as the basis for computations so that more complicated surface motions may be treated.

여기에서 논의된 나머지 자유 표면 방법은 보다 복잡한 표면 움직임을 처리할 수 있도록 고정된 오일러 그리드를 계산의 기준으로 사용합니다.

Surface Height Method

Low amplitude sloshing, shallow water waves, and other free-surface motions in which the surface does not deviate too far from horizontal, can be described by the height, H, of the surface relative to some reference elevation. Time evolution of the height is governed by the kinematic equation, where (u,v,w) are fluid velocities in the (x,y,z) directions. This equation is a mathematical expression of the fact that the surface must move with the fluid:

표면 높이 법

낮은 진폭의 슬로 싱, 얕은 물결 및 표면이 수평에서 너무 멀리 벗어나지 않는 기타 자유 표면 운동은 일부 기준 고도에 대한 표면의 높이 H로 설명 할 수 있습니다. 높이의 시간 진화는 운동학 방정식에 의해 제어되며, 여기서 (u, v, w)는 (x, y, z) 방향의 유체 속도입니다. 이 방정식은 표면이 유체와 함께 움직여야한다는 사실을 수학적으로 표현한 것입니다.

Finite-difference approximations to this equation are easy to implement. Further, only the height values at a set of horizontal locations must be recorded so the memory requirements for a three-dimensional numerical solution are extremely small. Finally, the application of free-surface boundary conditions is also simplified by the condition on the surface that it remains nearly horizontal. Examples of this technique can be found in References 1971 and 1975.

이 방정식의 유한 차분 근사를 쉽게 실행할 수 있습니다.  또한 3 차원 수치 해법의 메모리 요구 사항이 극도로 작아지도록 같은 높이의 위치 값만을 기록해야합니다.  마지막으로 자유 표면 경계 조건의 적용도 거의 수평을 유지하는 표면의 조건에 의해 간소화됩니다.  이 방법의 예는 참고 문헌의 1971 및 1975을 참조하십시오.

Marker-and-Cell (MAC) Method

The earliest numerical method devised for time-dependent, free-surface, flow problems was the Marker-and-Cell (MAC) method (see Ref. 1965). This scheme is based on a fixed, Eulerian grid of control volumes. The location of fluid within the grid is determined by a set of marker particles that move with the fluid, but otherwise have no volume, mass or other properties.

MAC 방법

시간 의존성을 가지는 자유 표면 흐름의 문제에 대해 처음 고안된 수치 법이 MAC (Marker-and-Cell) 법입니다 (참고 문헌 1965 참조).  이 구조는 컨트롤 볼륨 고정 오일러 격자를 기반으로합니다.  격자 내의 유체의 위치는 유체와 함께 움직이고, 그 이외는 부피, 질량, 기타 특성을 갖지 않는 일련의 마커 입자에 의해 결정됩니다.

Grid cells containing markers are considered occupied by fluid, while those without markers are empty (or void). A free surface is defined to exist in any grid cell that contains particles and that also has at least one neighboring grid cell that is void. The location and orientation of the surface within the cell was not part of the original MAC method.

마커를 포함한 격자 셀은 유체로 채워져있는 것으로 간주되며 마커가 없는 격자 셀은 빈(무효)것입니다.  입자를 포함하고, 적어도 하나의 인접 격자 셀이 무효인 격자의 자유 표면은 존재하는 것으로 정의됩니다.  셀 표면의 위치와 방향은 원래의 MAC 법에 포함되지 않았습니다.

Evolution of surfaces was computed by moving the markers with locally interpolated fluid velocities. Some special treatments were required to define the fluid properties in newly filled grid cells and to cancel values in cells that are emptied.

표면의 발전(개선)은 국소적으로 보간된 유체 속도로 마커를 이동하여 계산되었습니다.  새롭게 충전된 격자 셀의 유체 특성을 정의하거나 비어있는 셀의 값을 취소하거나 하려면 특별한 처리가 필요했습니다.

The application of free-surface boundary conditions consisted of assigning the gas pressure to all surface cells. Also, velocity components were assigned to all locations on or immediately outside the surface in such a way as to approximate conditions of incompressibility and zero-surface shear stress.

자유 표면 경계 조건의 적용은 모든 표면 셀에 가스 압력을 할당하는 것으로 구성되었습니다. 또한 속도 성분은 비압축성 및 제로 표면 전단 응력의 조건을 근사화하는 방식으로 표면 위 또는 외부의 모든 위치에 할당되었습니다.

The extraordinary success of the MAC method in solving a wide range of complicated free-surface flow problems is well documented in numerous publications. One reason for this success is that the markers do not track surfaces directly, but instead track fluid volumes. Surfaces are simply the boundaries of the volumes, and in this sense surfaces may appear, merge or disappear as volumes break apart or coalesce.

폭넓게 복잡한 자유 표면 흐름 문제 해결에 MAC 법이 놀라운 성공을 거두고 있는 것은 수많은 문헌에서 충분히 입증되고 있습니다.  이 성공 이유 중 하나는 마커가 표면을 직접 추적하는 것이 아니라 유체의 체적을 추적하는 것입니다.  표면은 체적의 경계에 불과하며, 그러한 의미에서 표면은 분할 또는 합체된 부피로 출현(appear), 병합, 소멸 할 가능성이 있습니다.

A variety of improvements have contributed to an increase in the accuracy and applicability of the original MAC method. For example, applying gas pressures at interpolated surface locations within cells improves the accuracy in problems driven by hydrostatic forces, while the inclusion of surface tension forces extends the method to a wider class of problems (see Refs. 1969, 1975).

다양한 개선으로 인해 원래 MAC 방법의 정확성과 적용 가능성이 증가했습니다. 예를 들어, 셀 내 보간 된 표면 위치에 가스 압력을 적용하면 정 수력으로 인한 문제의 정확도가 향상되는 반면 표면 장력의 포함은 방법을 더 광범위한 문제로 확장합니다 (참조 문헌. 1969, 1975).

In spite of its successes, the MAC method has been used primarily for two-dimensional simulations because it requires considerable memory and CPU time to accommodate the necessary number of marker particles. Typically, an average of about 16 markers in each grid cell is needed to ensure an accurate tracking of surfaces undergoing large deformations.

수많은 성공에도 불구하고 MAC 방법은 필요한 수의 마커 입자를 수용하기 위해 상당한 메모리와 CPU 시간이 필요하기 때문에 주로 2 차원 시뮬레이션에 사용되었습니다. 일반적으로 큰 변형을 겪는 표면의 정확한 추적을 보장하려면 각 그리드 셀에 평균 약 16 개의 마커가 필요합니다.

Another limitation of marker particles is that they don’t do a very good job of following flow processes in regions involving converging/diverging flows. Markers are usually interpreted as tracking the centroids of small fluid elements. However, when those fluid elements get pulled into long convoluted strands, the markers may no longer be good indicators of the fluid configuration. This can be seen, for example, at flow stagnation points where markers pile up in one direction, but are drawn apart in a perpendicular direction. If they are pulled apart enough (i.e., further than one grid cell width) unphysical voids may develop in the flow.

마커 입자의 또 다른 한계는 수렴 / 발산 흐름이 포함된 영역에서 흐름 프로세스를 따라가는 작업을 잘 수행하지 못한다는 것입니다. 마커는 일반적으로 작은 유체 요소의 중심을 추적하는 것으로 해석됩니다. 그러나 이러한 유체 요소가 길고 복잡한 가닥으로 당겨지면 마커가 더 이상 유체 구성의 좋은 지표가 될 수 없습니다. 예를 들어 마커가 한 방향으로 쌓여 있지만 수직 방향으로 떨어져 있는 흐름 정체 지점에서 볼 수 있습니다. 충분히 분리되면 (즉, 하나의 그리드 셀 너비 이상) 비 물리적 공극이 흐름에서 발생할 수 있습니다.

Surface Marker Method

One way to limit the memory and CPU time consumption of markers is to keep marker particles only on surfaces and not in the interior of fluid regions. Of course, this removes the volume tracking property of the MAC method and requires additional logic to determine when and how surfaces break apart or coalesce.

표면 마커 법

마커의 메모리 및 CPU 시간의 소비를 제한하는 방법 중 하나는 마커 입자를 유체 영역의 내부가 아니라 표면에만 보존하는 것입니다.  물론 이는 MAC 법의 체적 추적 특성이 배제되기 때문에 표면이 분할 또는 합체하는 방식과 시기를 특정하기위한 논리를 추가해야합니다.

In two dimensions the marker particles on a surface can be arranged in a linear order along the surface. This arrangement introduces several advantages, such as being able to maintain a uniform particle spacing and simplifying the computation of intersections between different surfaces. Surface markers also provide a convenient way to locate the surface within a grid cell for the application of boundary conditions.

2 차원의 경우 표면 마커 입자는 표면을 따라 선형으로 배치 할 수 있습니다.  이 배열은 입자의 간격을 균일하게 유지할 수있는 별도의 표면이 교차하는 부분의 계산이 쉽다는 등 몇 가지 장점이 있습니다.  또한 표면 마커를 사용하여 경계 조건을 적용하면 격자 셀의 표면을 간단한 방법으로 찾을 수 있습니다.

Unfortunately, in three-dimensions there is no simple way to order particles on surfaces, and this leads to a major failing of the surface marker technique. Regions may exist where surfaces are expanding and no markers fill the space. Without markers the configuration of the surface is unknown, consequently there is no way to add markers. Reference 1975 contains examples that show the advantages and limitations of this method.

불행히도 3 차원에서는 표면에 입자를 정렬하는 간단한 방법이 없으며 이로 인해 표면 마커 기술이 크게 실패합니다. 표면이 확장되고 마커가 공간을 채우지 않는 영역이 존재할 수 있습니다. 마커가 없으면 표면의 구성을 알 수 없으므로 마커를 추가 할 방법이 없습니다.
참고 문헌 1975이 방법의 장점과 한계를 보여주는 예제가 포함되어 있습니다.

Volume-of-Fluid (VOF) Method

The last method to be discussed is based on the concept of a fluid volume fraction. The idea for this approach originated as a way to have the powerful volume-tracking feature of the MAC method without its large memory and CPU costs.

VOF (Volume-of-Fluid) 법

마지막으로 설명하는 방법은 유체 부피 분율의 개념을 기반으로합니다. 이 접근 방식에 대한 아이디어는 대용량 메모리 및 CPU 비용없이 MAC 방식의 강력한 볼륨 추적 기능을 갖는 방법에서 시작되었습니다.

Within each grid cell (control volume) it is customary to retain only one value for each flow quantity (e.g., pressure, velocity, temperature, etc.) For this reason it makes little sense to retain more information for locating a free surface. Following this reasoning, the use of a single quantity, the fluid volume fraction in each grid cell, is consistent with the resolution of the other flow quantities.

각 격자 셀 (제어 체적) 내에서 각 유량 (예 : 압력, 속도, 온도 등)에 대해 하나의 값만 유지하는 것이 일반적입니다. 이러한 이유로 자유 표면을 찾기 위해 더 많은 정보를 유지하는 것은 거의 의미가 없습니다. 이러한 추론에 따라 각 격자 셀의 유체 부피 분율인 단일 수량의 사용은 다른 유량의 해상도와 일치합니다.

If we know the amount of fluid in each cell it is possible to locate surfaces, as well as determine surface slopes and surface curvatures. Surfaces are easy to locate because they lie in cells partially filled with fluid or between cells full of fluid and cells that have no fluid.

각 셀 내의 유체의 양을 알고 있는 경우, 표면의 위치 뿐만 아니라  표면 경사와 표면 곡률을 결정하는 것이 가능합니다.  표면은 유체 가 부분 충전 된 셀 또는 유체가 전체에 충전 된 셀과 유체가 전혀없는 셀 사이에 존재하기 때문에 쉽게 찾을 수 있습니다.

Slopes and curvatures are computed by using the fluid volume fractions in neighboring cells. It is essential to remember that the volume fraction should be a step function, i.e., having a value of either one or zero. Knowing this, the volume fractions in neighboring cells can then be used to locate the position of fluid (and its slope and curvature) within a particular cell.

경사와 곡률은 인접 셀의 유체 체적 점유율을 사용하여 계산됩니다.  체적 점유율은 계단 함수(step function)이어야 합니다, 즉, 값이 1 또는 0 인 것을 기억하는 것이 중요합니다.  이 것을 안다면, 인접 셀의 부피 점유율을 사용하여 특정 셀 내의 유체의 위치 (및 그 경사와 곡률)을 찾을 수 있습니다.

Free-surface boundary conditions must be applied as in the MAC method, i.e., assigning the proper gas pressure (plus equivalent surface tension pressure) as well as determining what velocity components outside the surface should be used to satisfy a zero shear-stress condition at the surface. In practice, it is sometimes simpler to assign velocity gradients instead of velocity components at surfaces.

자유 표면 경계 조건을 MAC 법과 동일하게 적용해야 합니다.  즉, 적절한 기체 압력 (및 대응하는 표면 장력)을 할당하고, 또한 표면에서 제로 전단 응력을 충족 시키려면 표면 외부의 어떤 속도 성분을 사용할 필요가 있는지를 확인합니다.  사실, 표면에서의 속도 성분 대신 속도 구배를 지정하는 것이보다 쉬울 수 있습니다.

Finally, to compute the time evolution of surfaces, a technique is needed to move volume fractions through a grid in such a way that the step-function nature of the distribution is retained. The basic kinematic equation for fluid fractions is similar to that for the height-function method, where F is the fraction of fluid function:

마지막으로, 표면의 시간 변화를 계산하려면 분포의 계단 함수의 성질이 유지되는 방법으로 격자를 통과하고 부피 점유율을 이동하는 방법이 필요합니다.  유체 점유율의 기본적인 운동학방정식은 높이 함수(height-function) 법과 유사합니다.  F는 유체 점유율 함수입니다.

A straightforward numerical approximation cannot be used to model this equation because numerical diffusion and dispersion errors destroy the sharp, step-function nature of the F distribution.

이 방정식을 모델링 할 때 간단한 수치 근사는 사용할 수 없습니다.  수치의 확산과 분산 오류는 F 분포의 명확한 계단 함수(step-function)의 성질이 손상되기 때문입니다.

It is easy to accurately model the solution to this equation in one dimension such that the F distribution retains its zero or one values. Imagine fluid is filling a column of cells from bottom to top. At some instant the fluid interface is in the middle region of a cell whose neighbor below is filled and whose neighbor above is empty. The fluid orientation in the neighboring cells means the interface must be located above the bottom of the cell by an amount equal to the fluid fraction in the cell. Then the computation of how much fluid to move into the empty cell above can be modified to first allow the empty region of the surface-containing cell to fill before transmitting fluid on to the next cell.

F 분포가 0 또는 1의 값을 유지하는 같은 1 차원에서이 방정식의 해를 정확하게 모델링하는 것은 간단합니다.  1 열의 셀에 위에서 아래까지 유체가 충전되는 경우를 상상해보십시오.  어느 순간에 액체 계면은 셀의 중간 영역에 있고, 그 아래쪽의 인접 셀은 충전되어 있고, 상단 인접 셀은 비어 있습니다.  인접 셀 내의 유체의 방향은 계면과 셀의 하단과의 거리가 셀 내의 유체 점유율과 같아야 한다는 것을 의미합니다.  그 다음 먼저 표면을 포함하는 셀의 빈 공간을 충전 한 후 다음 셀로 유체를 보내도록 위쪽의 빈 셀에 이동하는 유체의 양의 계산을 변경할 수 있습니다.

In two or three dimensions a similar procedure of using information from neighboring cells can be used, but it is not possible to be as accurate as in the one-dimensional case. The problem with more than one dimension is that an exact determination of the shape and location of the surface cannot be made. Nevertheless, this technique can be made to work well as evidenced by the large number of successful applications that have been completed using the VOF method. References 1975, 1980, and 1981 should be consulted for the original work on this technique.

2 차원과 3 차원에서 인접 셀의 정보를 사용하는 유사한 절차를 사용할 수 있지만, 1 차원의 경우만큼 정확하게 하는 것은 불가능합니다.  2 차원 이상의 경우의 문제는 표면의 모양과 위치를 정확히 알 수없는 것입니다.  그래도 VOF 법을 사용하여 달성 된 다수의 성공 사례에서 알 수 있듯이 이 방법을 잘 작동시킬 수 있습니다.  이 기법에 관한 초기의 연구 내용은 참고 문헌 1975,1980,1981를 참조하십시오.

The VOF method has lived up to its goal of providing a method that is as powerful as the MAC method without the overhead of that method. Its use of volume tracking as opposed to surface-tracking function means that it is robust enough to handle the breakup and coalescence of fluid masses. Further, because it uses a continuous function it does not suffer from the lack of divisibility that discrete particles exhibit.

VOF 법은 MAC 법만큼 강력한 기술을 오버 헤드없이 제공한다는 목표를 달성 해 왔습니다.  표면 추적이 아닌 부피 추적 기능을 사용하는 것은 유체 질량의 분할과 합체를 처리하는 데 충분한 내구성을 가지고 있다는 것을 의미합니다.  또한 연속 함수를 사용하기 때문에 이산된 입자에서 발생하는 숫자를 나눌 수 없는 문제를 겪지 않게 됩니다.

Variable-Density Approximation to the VOF Method

One feature of the VOF method that requires special treatment is the application of boundary conditions. As a surface moves through a grid, the cells containing fluid continually change, which means that the solution region is also changing. At the free boundaries of this changing region the proper free surface stress conditions must also be applied.

VOF 법의 가변 밀도 근사

VOF 법의 특수 처리가 필요한 기능 중 하나는 경계 조건의 적용입니다.  표면이 격자를 통과하여 이동할 때 유체를 포함하는 셀은 끊임없이 변화합니다.  즉, 계산 영역도 변화하고 있다는 것입니다.  이 변화하고있는 영역의 자유 경계에는 적절한 자유 표면 응력 조건도 적용해야합니다.

Updating the flow region and applying boundary conditions is not a trivial task. For this reason some approximations to the VOF method have been used in which flow is computed in both liquid and gas regions. Typically, this is done by treating the flow as a single fluid having a variable density. The F function is used to define the density. An argument is then made that because the flow equations are solved in both liquid and gas regions there is no need to set interfacial boundary conditions.

유체 영역의 업데이트 및 경계 조건의 적용은 중요한 작업입니다.  따라서 액체와 기체의 두 영역에서 흐름이 계산되는 VOF 법에 약간의 근사가 사용되어 왔습니다.  일반적으로 가변 밀도를 가진 단일 유체로 흐름을 처리함으로써 이루어집니다.  밀도를 정의하려면 F 함수를 사용합니다.  그리고, 흐름 방정식은 액체와 기체의 두 영역에서 계산되기 때문에 계면의 경계 조건을 설정할 필요가 없다는 논증이 이루어집니다.

Unfortunately, this approach does not work very well in practice for two reasons. First, the sensitivity of a gas region to pressure changes is generally much greater than that in liquid regions. This makes it difficult to achieve convergence in the coupled pressure-velocity solution. Sometimes very large CPU times are required with this technique.

공교롭게도 이 방법은 두 가지 이유로 인해 실제로는 그다지 잘 작동하지 않습니다.  하나는 압력의 변화에 대한 기체 영역의 감도가 일반적으로 액체 영역보다 훨씬 큰 것입니다.  따라서 압력 – 속도 결합 해법 수렴을 달성하는 것은 어렵습니다.  이 기술은 필요한 CPU 시간이 매우 커질 수 있습니다.

The second, and more significant, reason is associated with the possibility of a tangential velocity discontinuity at interfaces. Because of their different responses to pressure, gas and liquid velocities at an interface are usually quite different. In the Variable-Density model interfaces are moved with an average velocity, but this often leads to unrealistic movement of the interfaces.

두 번째 더 중요한 이유는 계면에서 접선 속도가 불연속이되는 가능성에 관련이 있습니다.  압력에 대한 반응이 다르기 때문에 계면에서 기체와 액체의 속도는 일반적으로 크게 다릅니다.  가변 밀도 모델은 계면은 평균 속도로 동작하지만, 이는 계면의 움직임이 비현실적으로 되는 경우가 많습니다.

Even though the Variable-Density method is sometimes referred to as a VOF method, because is uses a fraction-of-fluid function, this designation is incorrect. For accurately tracking sharp liquid-gas interfaces it is necessary to actually treat the interface as a discontinuity. This means it is necessary to have a technique to define an interface discontinuity, as well as a way to impose the proper boundary conditions at that interface. It is also necessary to use a special numerical method to track interface motions though a grid without destroying its character as a discontinuity.

가변 밀도 방법은 유체 분율 함수를 사용하기 때문에 VOF 방법이라고도하지만 이것은 올바르지 않습니다. 날카로운 액체-가스 인터페이스를 정확하게 추적하려면 인터페이스를 실제로 불연속으로 처리해야합니다. 즉, 인터페이스 불연속성을 정의하는 기술과 해당 인터페이스에서 적절한 경계 조건을 적용하는 방법이 필요합니다. 또한 불연속성으로 특성을 훼손하지 않고 격자를 통해 인터페이스 동작을 추적하기 위해 특수한 수치 방법을 사용해야합니다.

Summary

A brief discussion of the various techniques used to numerically model free surfaces has been given here with some comments about their relative advantages and disadvantages. Readers should not be surprised to learn that there have been numerous variations of these basic techniques proposed over the years. Probably the most successful of the methods is the VOF technique because of its simplicity and robustness. It is this method, with some refinement, that is used in the FLOW-3D program.

여기에서는 자유 표면을 수치적으로 모델링 할 때 사용하는 다양한 방법에 대해 상대적인 장점과 단점에 대한 설명을 포함하여 쉽게 설명하였습니다.  오랜 세월에 걸쳐 이러한 기본적인 방법이 많이 제안되어 온 것을 알고도 독자 여러분은 놀라지 않을 것입니다.  아마도 가장 성과를 거둔 방법은 간결하고 강력한 VOF 법 입니다.  이 방법에 일부 개량을 더한 것이 현재 FLOW-3D 프로그램에서 사용되고 있습니다.

Attempts to improve the VOF method have centered on better, more accurate, ways to move fluid fractions through a grid. Other developments have attempted to apply the method in connection with body-fitted grids and to employ more than one fluid fraction function in order to model more than one fluid component. A discussion of these developments is beyond the scope of this introduction.

VOF 법의 개선은 더 나은, 더 정확한 방법으로 유체 점유율을 격자를 통과하여 이동하는 것에 중점을 두어 왔습니다.  기타 개발은 물체 적합 격자(body-fitted grids) 관련 기법을 적용하거나 여러 유체 성분을 모델링하기 위해 여러 유체 점유율 함수를 채용하기도 했습니다.  이러한 개발에 대한 논의는 여기에서의 설명 범위를 벗어납니다.

References

1965 Harlow, F.H. and Welch, J.E., Numerical Calculation of Time-Dependent Viscous Incompressible Flow, Phys. Fluids 8, 2182.

1969 Daly, B.J., Numerical Study of the Effect of Surface Tension on Interface Instability, Phys. Fluids 12, 1340.

1970 Hirt, C.W., Cook, J.L. and Butler, T.D., A Lagrangian Method for Calculating the Dynamics of an Incompressible Fluid with Free Surface, J. Comp. Phys. 5, 103.

1971 Nichols, B.D. and Hirt, C.W.,Calculating Three-Dimensional Free Surface Flows in the Vicinity of Submerged and Exposed Structures, J. Comp. Phys. 12, 234.

1974 Hirt, C.W., Amsden, A.A., and Cook, J.L.,An Arbitrary Lagrangian-Eulerian Computing Method for all Flow Speeds, J. Comp. Phys., 14, 227.

1975 Nichols, B.D. and Hirt, C.W., Methods for Calculating Multidimensional, Transient Free Surface Flows Past Bodies, Proc. of the First International Conf. On Num. Ship Hydrodynamics, Gaithersburg, ML, Oct. 20-23.

1980 Nichols, B.D. and Hirt, C.W., Numerical Simulation of BWR Vent-Clearing Hydrodynamics, Nucl. Sci. Eng. 73, 196.

1981 Hirt, C.W. and Nichols, B.D., Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries, J. Comp. Phys. 39, 201.

업무에 적합한 올바른 CFD 소프트웨어 선택 방법

업무에 적합한 올바른 CFD 소프트웨어 선택 방법

많은 제품들이 모두 자신의 소프트웨어가 가장 적합하다고 말하기 떄문에, 사람들은 자신의 업무에 적합한 CFD 소프트웨어 선택에 어려움을 겪습니다. 그 이유는 유체 흐름 및 열 전달 분석을 위한 소프트웨어 패키지는 다양한 형태로 제공됩니다. 이러한 패키지는 물리적 근사치와 수치적 솔루션 기법이 크게 다르기 때문에 적합한 패키지를 선택하는 것이 어렵습니다.

아래 내용에서 올바른 CFD 소프트웨어를 선택할 때 고려해야 할 중요한 항목을 설명합니다.

Spillway’s tailrace over natural rock

1. 메싱 및 지오메트리

유한 요소 또는 “바디 맞춤 좌표”를 사용하는 솔루션 방법은 유동 영역의 기하학적 구조를 준수하는 해석용 그리드를 생성해야합니다. 정확한 수치 근사를 위해 허용 가능한 요소 크기와 모양으로 이러한 그리드를 생성하는 것은 쉽지 않은 작업입니다. 복잡한 경우 이러한 유형의 그리드 생성에는 며칠 또는 몇주의 노력이 소요될 수 있습니다. 일부 프로그램은 직사각형 그리드 요소만 사용하여 이러한 생성 문제를 제거하려고 시도하지만 흐름 및 열 전달 특성을 변경하는 “계단현상” 경계 문제를 해결해야 합니다. FLOW-3D는 FAVOR ™ (분수 면적 / 체적) 방법을 사용하여 기하학적 특성이 매끄럽게 포함된 생성하기 쉬운 직사각형 그리드를 사용하여 두 문제를 모두 해결합니다. 간단하고 강력한 솔리드 모델러가 FLOW-3D와 함께 패키지로 제공되거나 사용자가 CAD 프로그램에서 기하학적 데이터를 가져올 수 있습니다.

2. 운동량 방정식과 대략적인 흐름 모델

유체 운동량의 정확한 처리는 여러 가지 이유로 중요합니다. 첫째, 복잡한 지오메트리를 통해 유체가 어떻게 흐를지 예측할 수 있는 유일한 방법입니다. 둘째, 유체에 의해 가해지는 동적 힘 (즉, 압력)은 모멘텀을고려하여야만 계산할 수 있습니다. 마지막으로, 열 에너지의 대류 이동을 계산하려면 개별 유체 입자가 다른 유체 입자 및 제한 경계와 관련하여 어떻게 움직이는지를 정확하게 파악할 수 있어야 합니다.

이것은 운동량의 정확한 처리를 의미합니다. 모멘텀의 보존을 대략적으로만 하는 단순화된 흐름 모델은 실제적인 유체 구성과 온도 분포를 예측하는데 사용할 수 없기 때문에 FLOW-3D에서는 사용되지 않습니다.

3. 액체-고체 열 전달 영역

액체와 고체 (예 : 금속-금형) 사이의 열 전달에는 계면 영역의 정확한 추정이 필요합니다. 계단 경계는 이 영역을 과대 평가합니다. 예를 들어, 실린더의 표면적은 27 %의 비율로 과대 평가됩니다. FLOW-3D 전 처리기의 각 제어 볼륨에 대해 FAVOR ™ 방법에 의해 정확한 계면 영역이 자동으로 계산됩니다.

4. 액체-고체 열 전달에 대한 볼륨 효과 제어

제어 볼륨의 크기는 액체 / 고체 인터페이스를 포함하는 제어 볼륨에서도 열이 흐르기 때문에 액체와 고체 사이에서 교환되는 열의 속도와 양에 영향을 미칠 수 있습니다. FLOW-3D에서는 액체-고체 인터페이스에서 열 전달 속도를 계산할 때 체적 크기와 전도도가 고려됩니다.

5. 암시성(Implicitness)과 정확성

비선형 및 결합 방정식에 대한 암시적 방법에는 각 반복에서 under-relaxation 특성이 있는 반복 솔루션 방법이 필요합니다. 이 동작은 일부 상황에서 심각한 오류 (또는 매우 느린 수렴)를 일으킬 수 있습니다 (예 : 큰 종횡비로 제어 볼륨을 사용하거나 실제로 중요하지 않은 효과를 예상하여 암시성이 사용되는 경우).

FLOW-3D에서는 계산 노력FLOW-3D에서는 계산 작업이 덜 필요하기 때문에 가능한 경우 언제나 명시적 수치 방법을 사용하며, 수치 안정성 요구 사항은 정확도 요구 사항과 동일합니다. Implicit vs. Explicit Numerical Methods 문서에서 자세히 알아보세요.

6. 대류 전송을 위한 암시적 수치 방법 (Implicit Numerical Methods)

임의적으로 큰 시간 단계 크기를 계산에 사용할 수 있는 암시적 수치 기법은 CPU 시간을 줄이는데 널리 사용되는 방법입니다. 불행히도 이러한 방법은 대류 해석에 정확하지 않습니다. 암시적 방법은 근사 방정식에 확산 효과를 도입하여 시간 단계 독립성을 얻습니다. 물리적 확산(예 : 열전도)에 수치적 확산을 추가하는 것은 확산 속도만 수정하기 때문에 심각한 문제를 일으키지 않을 수 있습니다. 그러나 대류 과정에 수치 확산을 추가하면 모델링되는 물리적 현상의 특성이 완전히 바뀝니다. FLOW-3D에서 시간 단계는 프로그램에 의해 자동으로 제어되어 정확한 시간 근사치를 보장합니다.

7. 이완 및 수렴 매개 변수 (Relaxation and Convergence Parameters)

암시적 근사를 사용하는 수치 방법은 하나 이상의 수렴 및 이완 매개 변수를 선택해야합니다. 이러한 매개 변수를 잘못 선택하면 발산 또는 수렴 속도가 느려질 수 있습니다. FLOW-3D에서는 하나의 수렴 및 하나의 이완 매개 변수만 사용되며, 두 매개 변수는 프로그램에 의해 동적으로 선택됩니다. 사용자는 수치해석 솔버를 제어하는 ​​매개 변수를 설정할 필요가 없습니다.

8. 자유 표면 추적

액체-가스 인터페이스 (즉, 자유 표면)를 모델링하는 데 사용되는 두 가지 방법이 있습니다. 그 중 하나는 액체 및 가스 영역의 흐름을 계산하고 계면을 유체 밀도의 급격한 변화로 처리하는 것입니다. 일반적으로 밀도 불연속성은 고차 수치 근사를 사용하여 모델링됩니다.

불행히도, 이 치료는 몇몇 그리드 셀에 걸쳐 인터페이스가 매끄럽게 진행되도록 해주며, 그러한 인터페이스에 일반적으로 존재하는 접선 유속의 급격한 변화는 설명하지 않습니다. 또한 이 기법은 가스가 계산 영역으로 유입되는 액체로 대체될 경우 탈출 포트 또는 가스의 싱크로도 보완해야 합니다. 또한 이러한 방법은 일반적으로 유체의 비압축성을 만족시키기 위해 더 많은 노력을 기울여야 합니다.  가스 영역은 거의 균일한 압력 조정을 통해 솔루션 수렴 속도를 늦추는 경향이 있기 때문에 이러한 현상이 발생합니다.

FLOW-3D에서는 다른 기술인 VOF (Volume-of-Fluid) 방법이 사용됩니다. 이것은 인터페이스가 단계 불연속으로 긴밀하게 유지되는 진정한 3 차원 인터페이스 추적 체계입니다. 또한 선택적 표면 장력을 포함하여 수직 및 접선 응력 경계 조건이 인터페이스에 적용됩니다. 가스 영역은 사용자가 모델에 포함되도록 요청하지 않는 한 계산되지 않습니다.

FLOW-3D – CFD Software Simulation Gallery


FLOW-3D – CFD Software Simulation Gallery

FLOW-3D는 광범위한 산업 응용 분야 및 물리적 공정에서 액체 및 가스의 동적 거동을 연구하는 엔지니어를 위한 완전하고 다양한 CFD 시뮬레이션 플랫폼을 제공합니다. FLOW-3D는 자유 표면 및 다상 응용 분야에 매우 큰 강점을 가지고 있으며, 미세 유체 공학, 생물 의학 장치, 수자원 사회 기반 시설, 항공 우주, 소비재, 적층 제조, 잉크젯 인쇄, 레이저 용접, 자동차, 해양, 에너지 등 광범위한 산업에 사용됩니다.
https://www.flow3d.co.kr에서 FLOW-3D를 살펴봐 주시기 바랍니다.

 






Moving Boundaries: An Eulerian Approach

Moving Boundaries: An Eulerian Approach

많은 문제에서, 유체 및 고체 영역의 내부 경계가 그 안에서 이동할 수 있도록하면서 공간에 고정 된 그리드를 유지하는 것이 유리합니다. 이는 리 메싱의 필요성을 피할 수 있으므로 이러한 경계의 형태에 급격한 변화가 발생할 때마다 적절합니다. 메시 생성도 크게 단순화되었습니다.

고정 그리드 내에서 유체 인터페이스, 침전물, 응고 된 유체 및 탄성 재료의 경계 이동을 모델링하기위한 다양한 접근 방식이 표시됩니다. 유체 경계의 이동은 VOF (Volume-of-Fluid) 방법의 변형으로 수행되며, 각 계산 셀에서 유체의 양을 나타내는 양이 고정 메시를 통해 조정됩니다.

퇴적물의 침식 및 퇴적은 퇴적물 수색 모델을 사용하여 계산됩니다. 국부적 인 침식 속도는 패킹 된 퇴적물 / 유체 경계면에 존재하는 국부적 인 전단 응력을 기반으로하며, 증착은 Stokes 유동 근사치로 예측됩니다.

Emptying of gravure cell (same cell dimensions as filling case); a
three-dimensional perspective is shown. The transfer roll surface
(block at top) is moving away from the gravure roll at 0.5m/s. The
static contact of the fluid with all surfaces is 30°. The elapsed time
is 150

충진 층 경계면은 퇴적물 농도와 퇴적물의 포장 분율에 따라 달라집니다. 용융 금속은 온도가 빙점 아래로 떨어지면 굳을 수 있습니다. 응고 된 “유체”는 동결 및 용융을 유발하는 열유속의 양으로부터 결정된대로 표면이 증가하거나 수축하는 고체처럼 처리됩니다.

탄성 응력은 응고 된 재료 / 공기 인터페이스를 예측하는 VOF 방법을 사용하여 동일한 고정 그리드 내의 운동량 균형에 탄성 응력 계산을 추가하여 응고 된 영역에서 계산됩니다.

매우 일시적인 흐름 문제의 경우 유체와 공극 공간 사이 또는 두 개의 혼합 불가능한 유체 사이에있는 유체 인터페이스는 문제의 역학에 따라 자유롭게 움직여야합니다.

한 가지 해결책은 인터페이스와 함께 변형되는 메시를 만드는 것입니다. 이것은 시뮬레이션 중에 인터페이스의 형태가 거의 변경되지 않는 상황에서 잘 작동합니다. 그러나보다 일반적인 경우에는 시뮬레이션 중에 새 메시를 반복적으로 생성해야하거나 변경되지 않은 메시 내에서 자유 표면 경계를 생성하는 방법이 필요합니다. 이 작업은 후자를 제시합니다. VOF (Vol-of-fluid) 함수는 자유 표면의 위치를 추적하는 데 사용됩니다. 또한이 함수는 곡률을 계산하여 표면 장력의 영향을 예측하는 데 사용됩니다.

<원문보기> Moving-Boundaries-an-Eularian-Approach.pdf

접촉선의 고정(Contact Line Pinning)

접촉선의 고정(Contact Line Pinning)

증발하는 빗방울에서 남은 잔류의 물은 새로 씻은 자동차에서 좋지 못할 수 있습니다. 그러나, 동일한 증발 공정은, 예를 들어, 드롭 잔류 물이 인쇄 된 이미지 또는 텍스트의 일부가되는 잉크젯 인쇄에서 유리할 수있다. 그러나 동일한 증발 과정이 어떤 경우엔 도움이 될 수 있습니다 예를 들면, 잉크 찌꺼기가 인쇄 된 이미지나 텍스트의 일부가 되는 잉크젯 인쇄가 그렇습니다.

액체 방울의 증발로 인한 잔류의 물이 예상치 못한 방식으로 나타날 수 있습니다. 커피 링 얼룩이 잘 알려진 예이며, 커피의 잔류의 물이 물방울의 바깥 쪽 가장자리에 모여 얇은 원형 링 얼룩이 남습니다. 이 현상은 흥미로운 유체역학적인 과정의 결과입니다. 커피 링 얼룩이 형성 되려면 액체가 증착 된 고체 표면에 고정 된 접촉선이 있어야합니다. 고정 된 접촉선은 액체 방울이 고체 기판과 교차하는 액체 방울의 외부의 가장자리가 방울이 증발함에 따라 정지 상태를 유지함을 의미합니다. 증발은 기판의 열에 의해 발생하며 방울의 얇은 외부의 가장자리에서 가장 크게 생깁니다. 표면 장력은 액체가 증발하면서 손실 된 액체를 대체하기 위해 가장자리를 향해 발생하게 됩니다. 이는 결국 더 많은 용질을 가장자리로 운반하며 모든 액체가 증발 한 후, 결과적으로 커피 링 얼룩을 형성하게하는 더 높은 농도의 용질 잔류 물을 생성합니다.

모델링 접근법

FLOW-3D v12.0의 최신 업데이트로 인해 ‘접촉선의 고정’ 모델이 개발되었으며, 소프트웨어의 기능이 표면 장력 중심의 애플리케이션으로도 광범위하게 확장되었습니다. 표면 접촉의 고정 및 비고정 특성은 잉크젯 인쇄, 코팅 및 스프레이 냉각에서 중요한 역할을 합니다. 습윤 특성에 대한 표면 공법은 미세 유체 장치에서 액체 샘플의 이동을 제어하는 ​​데 사용될 수 있습니다. 모델의 주요 특징은 방울의 가장자리를 고정 위치에 고정하는 수단을 제공하는 것입니다. 형상 구성 요소 및 하위 구성 요소중에 표면에 ‘고정’ 속성을 지정할 수 있습니다. 유체의 접촉선은 처음 표면과 접촉하는 곳에 고정됩니다. 전방 속도를 0으로 유지하면 고정이 적용됩니다. 유체는 접촉선과 표면을 따라 이동하는 것이 아니라 롤오버하여 접촉점을 지나야만 이동할 수 있습니다.

커피 링 얼룩 검증

그림 1은 평평한 수평 표면에 놓인 원형 물방울의 결과를 보여줍니다. 표면은 30 ℃의 일정한 온도로 유지됩니다. 초기 유체 온도는 20 ℃이고 주변 공극의 온도는 일정한 20 ℃입니다. 유체는 밀도 0.967 g/cm3, 점도 0.02022 poise, 비열 1.645e+07 cm2/s/K, 열전도도 1.2964e+4 g*cm/s3/K, 표면 장력 계수 33.15 g/cm2의 일반적인 잉크를 나타냅니다.

그림 1. 고정 된 접촉선을 사용하여 건조 공정 중의 물방울 모양의 변화.

액적 표면의 초기 곡률 반경은 7.5e-03 cm이고, 차지하는 공간은 반경 4.5e-03 cm의 원이며, 겉보기의 초기 접촉각은 37.87 도입니다. 그림 1-a를 참조하시기 바랍니다. 지정된 정적 접촉각은 0 도입니다.

정압에 의한 상변화 모델이 활성화됩니다. 공극 내의 증기 분압은 0이고 상변화 수용 계수는 Rsize = 0.01 입니다.

잉크가 건조될 때 기판 상에 고체가 잔류하는 물이 형성되는 것을 포착하기 위해 잔류 물 모델도 켜집니다. 유체에 용해 된 안료의 농도는 초기 농도 0.01 g/cm3 이고 최대 농도 rmax = 1.1625 g/cm3 에서 운반이 가능한 스칼라로 표시됩니다. 용해 된 안료는 질량 평균을 기준으로 안료의 단위질량당 0.05 poise의 속도로 유체의 순 점도를 향상시킵니다.

이 공정은 3.0 도의 방위 방향으로 하나의 셀에 걸쳐있는 축 대칭 원통형 메쉬로 모델링됩니다. (x 간격 = 6e-05 cm, z 간격 = 4e-05 cm.)

그림 1은 유체가 증발함에 따라 접촉선이 고정 된 상태를 유지하고 있음을 보여줍니다. 0 도의 정적 접촉각 조건은 액적의 중심을 향한 압력 구배를 가져오고, 이는 접촉선 방향으로의 유동을 생성합니다. 용해 된 안료의 농도는 증발로 인해 자유 표면 근처에서 증가하며, 흐름을 따라 농도는 접촉선을 향해 더욱 재분배합니다. (그림 2). 액체가 계속 증발함에 따라, 남아있는 액체의 안료 농도는 증가합니다. 농도가 최대 rmax에 도달하면, 과잉된 안료는 고체가 잔류하는 물로 전환됩니다.

그림2. g / cm3 단위의 안료 농도 및 t = 2.0ms에서의 흐름 패턴. 흐름은 고정 된 접촉선을 향하여 안료 농도가 증가합니다.

접촉선 근처의 유체가 먼저 건조되어 고체가 잔류하는 물이 남습니다. 해당 영역의 유체에 안료 농도가 높기 때문에 고체가 잔류하는 물의 특징인 ‘커피 링’ 패턴이 기판 표면에 생성됩니다. (그림 3 및 4). 안료의 총 질량(용해 + 건조 잔류 물)은 초기 질량의 0.025 % 이내로 보존됩니다.

그림 3. 모든 유체가 증발 된 후 기판 표면에 건조된 잔류 물의 분포 (단위 : g / cm3) .
가장 높은 농도는 고정 된 접촉선의 위치에 있으며, 이는 ‘커피 링’ 효과를 만들어냅니다.
그림 4. 유체가 완전히 증발 한 후 초기 액적의 반경을 따라 건조된 잔류 물의 예상 분포.

물방울 벽의 검증

그림5. 수직 벽에 고정 된 물방울의 변형 : t = 0 ms (파란색), t = 4e-02 ms (연한 파랑) t = 0.2 ms (빨간색).
해당 이미지는 “Effects of microscale topography”, Y.V.Kalinin, V.Berejnov and R. E. Thorne, Langmuir 25, 5391-5397. (2009). 에서의 이미지입니다.

접촉선 고정 응용의 두 번째 예는 수직의 벽에 고정 된 한 방울의 액체 알루미늄의 거동입니다. 유체 밀도는 2.7 g / cm3, 표면 장력 계수 200 g / cm2 및 점도 0.27 poise입니다. 정적 접촉각은 0 도입니다.

초기의 겉보기의 접촉각이 90도가 되도록 반경 0.5cm의 물방울을 수직 벽에 놓습니다 (그림 5). 7e+06 cm/s2의 중력 크기는 표면 장력의 복원 작용을 없애고 액적이 눈에 띄도록 변형시키기 위하여 인위적으로 향상되었습니다. 결과들은 비슷한 크기의 물방울에 대한 실험 결과와의 질적 비교를 포함하여 그림 5에서 보여줍니다.

요약

FLOW-3D의 접촉선 고정 모델은 표면 장력 및 벽의 접착 기능을 확장하여 표면 공법에서 복잡한 상호 작용을 모델링합니다. 접촉선 고정이 실제로 응용되는 분야에 관하여 더 많은 예시와 추가적인 참조를 찾으신다면 여기에서 찾을 수 있습니다.

Spraying

Spraying

Dispensing Liquids with Swirl Spray Nozzles

소용돌이-스프레이(Swirl-spray) 노즐은 화학 청소기, 의약품 및 연료에 액체를 분사하는 일반적인 방법입니다. 액체의 성공적인 분무화를 위해서는 노즐에 침투하는 공기 코어(air core)를 형성해야 합니다. CFD는 최적의 스프레이 콘을 위한 기하학, 스월 속도(swirl velocity) 및 유체 특성의 영향을 탐구하는 효과적인 방법입니다.

FLOW-3D simulation of a swirl spray nozzle

이 예에서는 2차원 축 대칭 스월 흐름(axisymmetric swirl flow)이 시뮬레이션되었습니다. 대칭 축을 따른 공기 코어가 노즐 전체 길이까지 거의 침투했습니다. 왼쪽 그림은 벡터가 평면에서 속도 분포를 나타내는 압력 분포입니다. 오른쪽 그림은 속도의 스월 성분에 의해 색상이 지정되고 빨간색은 더 높은 값을 나타냅니다.

스프레이 콘과 방울 크기의 눈금이 너무 넓기 때문에 스프레이의 전체 분자를 직접 계산할 수 없습니다. 또한, 원자화는 외부 동요, 노즐의 미세한 결함 및 기타 영향과 밀접하게 연관된 혼란스러운 과정입니다. 그러나 스프레이 콘이 노즐을 빠져나갈 때(즉, 벽 두께, 원뿔 각도, 축 속도 및 방위각 속도) 스프레이 콘의 특성을 예측할 수 있다는 것은 이러한 유형의 흐름 장치를 최적화하는 데 큰 도움이 됩니다.

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)과 포화도 사이의 비선형 관계를 모형화할 수 있습니다.

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

Non-Newtonian Fluids

Non-Newtonian Fluids

혈액, 케첩, 치약, 샴푸, 페인트 및 로션과 같은 비 뉴턴 유체는 점도가 다양한 복잡한 유변학을 가지고 있습니다. FLOW-3D는 변형 및 온도에 따라 달라지는 비 뉴턴 점도를 가진 유체를 모델링합니다. 전단 및 온도 의존 점도는 Carreau, 거듭 제곱 법칙 함수 또는 단순히 표 형식 입력을 통해 설명됩니다. 일부 폴리머, 세라믹 및 반고체 금속의 특성인 시간 의존적 또는 요 변성 거동(thixotropic behavior)도 시뮬레이션 할 수 있습니다.

Hand Lotion Pump

핸드 로션 펌프는 종종 몇 가지 설계 문제와 관련이 있습니다. 펌프가 공극을 막지 않고 효과적으로 작동하고 로션을 연속적으로 생성하는 것이 중요합니다. 좋은 디자인은 노력을 덜 필요로하며 이상적으로는 로션을 원하는 위치로 향하게합니다. FLOW-3D의 움직이는 물체 모델은 노즐이 아래로 밀리는 것을 시뮬레이션하여 저장소의 로션을 가압하는 데 사용됩니다. 로션의 압력과 로션을 추출하는 데 필요한 힘을 연구 할 수 있습니다. 동일한 고정 구조화 된 메시 내에서 여러 설계 변수를 쉽게 분석 할 수 있습니다.

FLOW-3D’s TruVOF method accurately captures the pulsating lotion as the ball regulates the frequency of dispensing lotion.

Food Processing

Food Processing

식품 가공 산업은 분사 장비를 최적으로 설계 및 제조하기 위해 일반적으로 뉴턴이 아닌 유체, 슬러리, 고형물과 유체의 혼합물 등 복잡한 유체를 관리하기 위한 다양한 요구사항을 가지고 있습니다. 이는 상업용 등급 장비의 일관성 및 내구성과 품질에 필수적입니다. 또한, 포장 디자인의 혁신은 한 제품을 다른 제품과 명확하게 구별할 수 있습니다. 예를 들어, 꿀, 케첩 또는 크리머(creamer)를 깨끗하고 정확하게 분사하는 것은 소비자가 상점에서 선택하는 것일 수 있습니다. 운송 및 보관 요구 사항은 모양을 보다 효율적으로 엔지니어링하고 컨테이너 자재를 더 많이 선택해야 합니다. 1.5리터의 물 또는 세탁 세제를 이동 또는 떨어뜨리는 동안 유체 하중은 업스트림 설계의 중요한 부분이 될 수 있습니다.

Velocity contours of a non-Newtonian fluid such as ketchup out of a slit nozzle.

꿀, 옥수수 시럽, 치약과 같은 점성 액체는 보통 단단한 표면에 접촉할 때 코일을 형성하는 경향을 보입니다. 이러한 효과는 관찰하기에 흥미롭고 재미있지만, 제품 안에 공기가 들어가 포장을 어렵게 만들 수 있는 포장 공정에서는 달갑지 않을 수 있습니다.

코일링이 발생하는 조건은 유체의 점성, 유체가 떨어지는 거리 및 유체의 속도에 따라 달라집니다. FLOW-3D는 다양한 물리적 프로세스 매개변수를 연구하기 위한 정확한 도구를 제공하여 효율적인 프로세스를 설계하는 데 도움이 됩니다.

Flushing

Flushing

화장실이 어떻게 작동하는지 궁금한 적이 있습니까? 사실 꽤 복잡합니다. 손잡이를 밀면 물이 용기를 채우기 시작합니다. 용기의 유체 레벨이 트랩 상단 (보울 뒤) 위로 올라가면 위어 유형의 흐름이 시작됩니다. 흐름이 충분히 빠르면 트랩 상단에 거품이 형성되어 사이펀이 생성됩니다. 그 시점에서 사이펀은 용기에서 물을 빼내고 변기가 내립니다.

많은 지역에서 물 절약은 중요한 문제이며 가정과 상업용 모두에 저 유량 화장실이 필요합니다. 그러나 화장실이 첫 번째 시도에서 작업을 완료하지 못하면 물 절약 목표가 실패합니다. FLOW-3D는 최적의 결과를 얻기 위해 다양한 설계를 모델링하는 데 사용할 수 있습니다.

Toilet Flushing Examples

아래 3D 애니메이션에서 FLOW-3D는 물 동작의 세척 순서를 보여줍니다. 물의 두 영역이 공과 함께 초기화됩니다. 공은 6 자유도의 완전 결합 유체-고체 모션을 시뮬레이션하기 위해 움직이는 물체 모델(GMO)을 사용하여 모델링됩니다. 중력은 수세식 탱크에서 물을 용기로 밀어 넣습니다. 분석은 정체 영역과 공이 영역을 벗어나는 기간을 나타내는 흐름 프로파일과 압력 윤곽을 보여줍니다. 공 대신 다른 질량과 모양을 사용할 수 있습니다. 플러싱 과정에서 잔여 물도 분석 할 수 있습니다.

아래의 횡단면 플롯은 수조의 흐름 재순환과 상세한 흐름 프로필을 보여줍니다. Collision 모델은 규정된 반발 및 마찰 계수를 기반으로 바운싱을 예측하는 공을 시뮬레이션하는 데 사용되었습니다. 물과 공기 사이의 일시적인 예리한 경계면은 FLOW-3D의 TruVOF 방법을 사용하여 잘 유지됩니다.

Filling / 충진

Filling / 충진

재료 비용을 줄이고 사이클 시간을 개선하기 위해 소비재 회사는 슬로 싱, 튀기 및 공기 혼입을 포함한 많은 자유 표면 유체 문제를 처리해야합니다.

Predicting Entrained Air in a Bottle Filling Example

혼입된 공기는 생산 라인에서 컨테이너가 채워질 때 액체의 부피를 증가시킬 수 있습니다. 아래 왼쪽 이미지는 높이가 약 20cm인 병에 1.2 초 동안 채우는 것을 보여줍니다. 색상 음영은 액체에서 공기의 부피 비율을 나타냅니다. 병에서 짧은 시간과 높은 수준의 혼합으로 인해 공기가 표면으로 올라와 빠져 나갈 시간이 없었습니다. 그러나 오른쪽 이미지에서 볼 수 있듯이 약 1.7 초의 추가 시간이 지나면 표면으로 상승하는 공기로 인한 액체 부피 감소가 명확하게 보입니다. FLOW-3D의 드리프트 플럭스 모델을 사용하면 액체의 기포와 같은 성분을 분리하여 분리 할 수 있습니다.

Air entrainment (left) and separation of air and liquid (right)

In by 9, out by 5 – Rapid evaluation of Tide® bottle filling

FLOW-3D를 사용하여 새로운 Tide 병 디자인의 채우기를 모델링하는 방법을 설명하는이 기사는 The Procter and Gamble Company의 기술 부문 책임자 인 John McKibben이 기고했습니다.

오전 9시에 긴급한 이메일을 받았다고 상상해보십시오.

새로운 Tide® 병 디자인 중 하나가 핸들을 채우고 충전 장비에 문제가있을 수 있음을 방금 깨달았습니다. 프로토 타입 병도없고 몇 주 동안도 없을 것입니다. 디자이너와 소비자는 디자인의 모양을 좋아하지만 그것이 채우는 방식은 우리 생산 시설의 쇼 스토퍼가 될 수 있습니다.
이 상황을 접했을 때 저는 3D 지오메트리 (그림 1)의 스테레오 리소그래피 (.stl) 파일을 요청하여 응답을 시작했고 제가 할 수있는 일을 확인했습니다. FLOW-3D는 .stl 파일을 사용하여 지오메트리를 입력 할 수 있으며 채우기에 대한 자유 표면 문제를 해결할 수 있어야한다는 것을 알고있었습니다. 나는 이것이 잠재적 인 문제에 대한 좋은 질적 이해를 제공 할 것으로 기대했지만,이 응용 프로그램에 대해 얼마나 정확한지에 대해서는 약간 불확실했습니다.

Setting up and Running the Simulation

오후 1 시경에 지오메트리 파일, 유량 및 유체 속성을 받았습니다. 몇 시간 내에 시뮬레이션이 실행되어 예비 결과를 제공했습니다. 저는 제 고객을 초대하여 결과를 간단히 살펴 보았고 그는 “보스의 상사”도 함께 살펴 보았습니다. 그래서 저녁 5 시까 지 예비 결과를보고 원래 우려 사항이 문제가 아니라고 판단했습니다.

그러나 결과는 몇 가지 다른 질문을 제기했습니다. 핸들을 채우면 유입되는 유체 분사가 많이 분리되었습니다. 나는 이것이 동반 된 공기와 거품의 양을 증가시킬 것이라는 것을 알고 있었다 (우리는 결국 세탁 세제를 채우고있다). FLOW-3D 공기 혼입 모델을 테스트하기로 결정했습니다. 이 모델은 원래 난류 제트 용으로 개발되었으며,이 층류 문제를 볼 때 얼마나 잘 수행 될지 확신 할 수 없었습니다.

Figure 2: Filled results
Experimental comparison of bottle filling model with and without the air entrainment model, courtesy of The Procter & Gamble Company.

그림 2는 공기 유입 모델이 있거나없는 병 충전 모델의 결과를 보여줍니다. 혼입 된 공기가 포함되면 충전 레벨이 크게 증가합니다. 혼입 된 공기가 병 상단에서 액체를 밀어 내지는 않지만 공기 혼입 정확도를 확인해야 할만큼 충분히 가깝습니다.

그림 3은 몇 주 후에 실행 된 실험의 이미지와 공기 혼입 수준을 비교합니다 (시제품 병이 제공되었을 때). 제트 분리 및 충진 수준의 정 성적 일치는 우수하며 시뮬레이션이 병 설계를 선별하기에 충분히 정확하다는 것을 확인했습니다.

Energy

Energy

전 세계 에너지 부문의 엔지니어는 전산 유체 역학(CFD)을 통해 해결책을 찾기 위해 광범위한 프로세스에서 매일 복잡한 설계 문제에 직면합니다. 특히 자유 표면 흐름과 관련이 높은 이러한 문제의 대부분은 FLOW-3D가 매우 정확한 분석을 제공하여 문제 해결에 적합합니다.

  • Fuel or cargo sloshing inside containers on the high seas
  • Wave effects on offshore platforms
  • Performance optimization for separation devices undergoing 6 DOF motion
  • Design of devices to capture energy from waves

Energy Case Studies

천연자원이 계속 감소함에 따라, 대체 자원과 방법을 탐구하고 가능한 한 효과적으로 현재 공급량을 사용하고 있습니다. 엔지니어는 사고를 예방하고 채굴 및 기타 에너지 수확 기법으로 인한 환경적 영향을 평가하기 위해 FLOW-3D를 사용합니다.

잔해물 수송 테스트를 CFD 모델링과 결합하면 ECCS 스트레이너가 견딜 수 있어야 하는 잔해물 부하를 다른 방법으로는 가정해야하는 지나치게 보수적 인 값에서 크게 줄일 수 있습니다. CFD는 또한 수두 손실 테스트를 지원하기 위해 ECCS 스트레이너 주변의 흐름 패턴, 수두 손실 테스트 및 플랜트 설계 수정을 식별하는데 있어 containment pool 수위 변화를 식별하는데 유용함이 입증되었습니다.
( By combining debris transport testing with CFD modeling, the debris loads that the ECCS strainers must be able to withstand can be significantly reduced from the overly conservative values that must otherwise be assumed. CFD has also proved to be valuable in identifying containment pool water level changes, flow patterns in the vicinity of the ECCS strainer to support head loss testing, and plant design modifications.
-Tim Sande & Joe Tezak, Alion Science and Technology)


Architects Achieve LEED Certification in Sustainable Buildings

Comparing HVAC System Designs

Debris Transport in a Nuclear Reactor Containment Building

Evaluation of the Wind Effects on the Iron-Ore Stock Pile

Mining

Mining

점성이 높은 유체, 비 뉴턴 흐름, 슬러리 또는 심지어 입상 흐름의 형태를 취할 수있는 많은 채광 응용 프로그램에서, 남은 잔여물인 테일링은 까다로운 시뮬레이션 문제를 제공합니다. FLOW-3D는 비 뉴턴 유체, 슬러리 및 입상 흐름에 대한 특수 모델을 포함하여 이러한 분석을 수행하는 데 필요한 모든 도구를 제공합니다. FLOW-3D의 자유 표면 유동 모델링 기능과 결합되어, 이러한 어렵고 환경 적으로 민감한 문제에 대한 탁월한 모델링 솔루션을 제공합니다.

풍력 강제 분석 하에서 광석 비축더미 먼지 표류 현상은 밀접하게 관련되어 있습니다. 여기서 FLOW-3D의 드리프트 플럭스 모델을 사용하면 엔지니어가 광석 침착 및 혼입 패턴과 개선 솔루션의 효과를 연구 할 수 있습니다.

Mining Simulations

Example of a tailing breach simulation using FLOW-3D‘s granular flow model.

액화 작용과 기계적 방해 모두 균열의 동적 특징의 일부이며, 이는 결국 물과 같은 뉴턴의 흐름과는 대조적으로 미세한 흐름의 매우 독특한 특성으로 자신을 멈추게 됩니다.

Oil & Gas Separators

Oil & Gas Separators

FLOW-3D는 오일 및 물과 같은 혼합 불가능한 유체를 모델링할 수 있으며, 개방된 환경(주변 공기)에 관련된 구성 요소 간의 고유한 인터페이스를 정확하게 추적할 수 있습니다.

전체 영역에 영향을 미치는 역학 관계로 유체가 자유롭게 혼합될 수 있습니다. 시간이 지남에 따라 유체는 연속 위상과 분산 위상 사이의 표류 관계에 따라 다시 분리됩니다.

중력 분리기의 성능은 CFD 모델링을 통해 향상 될 수 있습니다 :

  • Develop vessel inlet configurations to improve uniformity of gas and liquid flows and avoid mixing of oil & water due to wave-induced sloshing
  • Determine the influence of internal equipment on hydraulic efficiency and separation performance
  • Gauge the effects of changes in operating conditions
  • Accurately model small-scale phenomena (multi-phase flows, drops, particles, bubbles)

Sloshing & Slamming

Sloshing & Slamming

LNG 선과 같은 전 세계 항해 선박의 내부 슬로싱 운동에 의해 유발되는 가속 및 하중은 이러한 선박의 안전 설계에 중요한 요소입니다. 선박은 이러한 파도 슬로싱 역학으로 인해 상당한 내부 강제력을 경험할 수 있으며,화물 탱크는 격리 시스템 및 지지 구조물에 추가 하중을 경험할 수 있습니다. FLOW-3D의 비관성 기준 프레임 모델을 사용하면 컨테이너에서 유체 모션을 정확하게 추적하기 위해 복잡한 모션 매개 변수를 규정 할 수 있습니다. FLOW-3D는 연료 탱크의 액체화물과 추진체 움직임을 정확하게 예측합니다.

슬래밍 분석과 “그린 워터”(파동 스프레이 또는 선체 오버 토핑의 결과로 데크를 덮는 물)는 FLOW-3D의 일반적인 응용 분야이며, 이러한 문제에 대한 FLOW-3D의 TruVOF 접근 방식은 엔지니어에게 고유 한 기능을 제공합니다.

이 시뮬레이션은 FLOW-3D로 캡처한 것처럼 파도로 움직이는 LNG 유조선의 복잡하고 완전히 결합된 동작을 보여줍니다. 탱크에서 유체의 슬로싱은 선박의 역학에 영향을 미치며, 여기에는 bilge keels 과 bulbous bow가 포함되어 입사파에서 안정성을 유지하는데 도움이됩니다. 또한 탱크의 유체 슬로싱은 예상대로 선박의 움직임에서 위상 이동이 있습니다. FLOW-3D의 움직이는 물체 모델, 밀도 평가 모델 및 파동 생성 경계 조건은 해양 엔지니어가 유사한 문제를 분석하는 데 도움이 될 수 있습니다.

Optofluidics

Optofluidics

광유체학(Optofluidics)은 광학 분야와 미세유체학 분야를 합친 것입니다. 대부분의 광유체 응용은 렌즈를 만들기 위해 다양한 굴절률의 유체를 사용하는 것을 포함합니다. 이 접근 방식의 주요 장점은 렌즈의 동적 재구성 가능성입니다. 또한, 미세 유체 흐름은 이전에는 달성할 수 없었던 해상도를 달성하는 현미경에 이 기술을 쉽게 통합할 수 있도록 합니다. 마이크로 유체 소자에 빛을 집중시키는 마이크로 유체 소자 분야에서는 광학적 특성화가 필요한 랩온 어 칩 애플리케이션에 중요합니다. 마이크로 채널에서 나오는 빛의 효율적인 조명과 반사는 장치의 성능에 매우 중요합니다.

FLOW-3D는 L2 (Liquid core-liquid cladding) 렌즈와 같은 미세 유체 렌즈의 형성과 관련된 유체 역학을 시뮬레이션하는 데 사용됩니다.

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.

Thermocapillary Actuation

Thermocapillary Actuation

표면 장력의 온도 의존성은 유체 방울을 패턴 있는 표면 위로 전달하는 데 사용될 수 있습니다. 표면은 유체 방울이 친수성-수소성 인터페이스에 의해 형성된 채널을 따르도록 제한되도록 친수성 또는 친수성 접촉 각도로 패턴화됩니다. 또한 프로그램 가능한 방식으로 가열된 마이크로 히터의 배열은 열전압 작동을 유발하여 유체 방울을 뜨거운 영역에서 차가운 지역으로 유도합니다. 아래 이미지는 문제 설정의 상단 및 단면 뷰(Anton A)를 보여줍니다. Darhuber 외.) 다음에 Flow-3D를 설정합니다.

Liquid droplet moving along hydrophilic microstripe
Top-view of a liquid droplet moving along a hydrophilic microstripe. The array of Ti-resistors (shown in light gray) beneath the hydrophilic stripes locally heat the droplet thereby modifying the surface tension and propelling the liquid toward the colder regions of the device surface. The dark gray stripes represent the leads and contacts (Au) for the heating resistors.
Cross sectional view of device
Cross-sectional view of a portion of the device containing two micro-heaters and an overlying droplet.

더 차가운 표면 온도 영역은 인접한 따뜻한 지점보다 더 높은 표면 장력을 유지하여 액체를 당기는 접선 표면 힘을가합니다. 부분적 습윤 (접촉각> 0) 표면은 전체 습윤 표면 (접촉각 = 0)에 비해 부피 손실이 적은 유체 수송을 허용하기 때문에 바람직한 옵션입니다.

FLOW-3D setup of three microheaters

Top view of the setup in FLOW-3D showing three microheaters in pink, yellow and blue respectively. The central hydrophilic strip is shown in black with a fluid (water) droplet in sky blue.

아래 애니메이션은 완전히 젖은 표면과 부분적으로 젖은 표면의 비교를 보여줍니다. 예상대로 완전히 젖은 표면은 부분적으로 젖은 표면보다 액적을 더 평평하게 (그리고 더 많이 퍼지게) 만듭니다. 히터가 한 번에 하나씩 활성화되면 물방울이 더 차가운 영역으로 이동됩니다. 더 많은 유체가 남겨질수록 시뮬레이션이 끝날 때까지 완전히 젖은 표면은 더 많은 유체 볼륨을 잃는 것을 볼 수 있습니다. 따라서 부분적으로 젖은 표면은 유체 손실을 줄이기위한 더 바람직한 옵션입니다. 두 경우 모두 소수성 표면으로 둘러싸인 중앙 친수성 스트립으로 인해 물방울이 중앙에 머물러야합니다.

Animation of the results post-processed in FlowSight.

References

Anton A. Darhuber, Joseph P. Valentino, Sandra M. Trian and Sigurd Wagner, Thermocapillary Actuation of Droplets on Chemically Patterned Surfaces by Programmable Microheater Arrays, Journal of Microelectrochemical Systems, Vol. 12, No. 6, December 2003

Electrokinetics

Dielectrophoresis

유전 영동은 분극성 입자에 힘을 생성하여 균일하지 않은 전기장 (일반적으로 AC 전기장)에서 움직임을 유도합니다. 유전 영동력은 마이크로스케일 및 나노스케일 바이오 입자를 특성화, 처리 또는 조작하는 데 사용할 수 있습니다. 여기에는 세포, 바이러스, 박테리아, DNA 등의 분류, 포획 및 분리가 포함될 수 있습니다. 유전 영동은 FLOW-3D에서 완전히 설명 할 수 있으며 날카로운 인터페이스가 있거나 없는 단일 유체 또는 2 유체 흐름과 같이 코드에서 사용할 수있는 다른 모든 유체 흐름 옵션과 함께 활성화 될 수 있습니다.

Electro-wetting

전도성 액적에서 액체와 전극 사이에 인가되는 얇은 유전체 코팅 전위를 갖는 전극 상에 배치되면, 드롭 평면화와 전극 표면 확산이 일어납니다. 이 현상은 종종electro-wetting라 부릅니다. 현상은 전하 층의 발달과 관련되어 있으므로, 외부 전기장을 그들을 이동, 합체, 깨지거나 하는 원인을 조작하기 위해 사용될 수 있습니다.

 

Lab-On-Chip Electro-wetting Applications

Lab-on-chip 기반electro-wetting 은 분리된 물방울을 조절할 수 있어 설계자들이 복잡한 절차를 전통적인 실험실 장치를 달지만 훨씬 작은 volumes 으로 비슷한 실험을 수행할 수 있습니다. 이러한 기기는 효율적으로 운송, 병합되어 있으며 분리된 물방울들이 요구합니다. FLOW-3D는 사용자가이 장치를 조작하는 데 사용되는 기하학적 파라미터들 및 전압의 영향을 시뮬레이션 할 수 있도록 하여 설계 프로세스에 유용한 도구가 될 수 있습니다.

아래의 애니메이션은 수송 시뮬레이션 병합 및 분할 방울에 FLOW-3D의 기능을 보여줍니다. Lab-on-chip은 약 300 ㎛로 분리 된 두 개의 평행 한 플레이트로 구성됩니다. 바닥 판은 방울을 조작하기 위해 사용되는 그 안에 삽입 된 전극을 보유하고 있습니다. 액 적은 물 (약간 도전성) 실리콘 오일에 의해 둘러싸여 있습니다. 액체 방울의 부피가 800nl 관한 것입니다.

This lab-on-a-chip electrowetting simulation demonstrates an electric field being applied in order to split a small droplet.

Here an electric field is being applied in order to merge two small droplets.

This simulation shows an electric field being applied to a small droplet to control its motion.

Droplet Based Microfluidics

Droplet Based Microfluidics

연속 미세 유체와 달리 액적 기반 미세 유체는 개별 볼륨의 유체만 조작합니다. 마이크로 스케일 액적 시스템은 액적 역학에 대한 깊은 이해가 있는 경우 높은 처리량을 허용 할 수도 있습니다. 전산 유체 역학은 이러한 시스템의 동작을 이해하고 예측하는데 매우 유용한 도구입니다. 수치 시뮬레이션을 통한 액적 역학 연구는 잉크젯 기술의 확산에 중요한 역할을했습니다. 지난 10 년 동안 FLOW-3D는 상업 및 학술 응용 분야 모두에서 이러한 연구에 선호되는 분석 도구로 확립되었습니다. FLOW-3D는 음향 유도 잉크젯, 피에조 잉크젯, 열 거품 잉크젯 및 기타 여러 유형을 연구하는데 사용되었습니다.

강력한 표면 장력 모델과 성공적인 잉크젯 모델링 이력을 갖춘 FLOW-3D는 자연스럽게 액적 기반 미세 유체 공정 모델링으로 확장됩니다. FLOW-3D는 오늘날 다른 응용 분야 중에서도 정밀 액적 생성, 정밀 액적 증착, 액적 병합, 액적 분리, 액적 움직임, 흐름 집중, 잉크젯 인쇄 및 액적 T 접합을 시뮬레이션하는 데 사용됩니다.

Co flow fluid dynamics
Co-Flow
Droplet based microfluidics
Flow Focusing
Piston driven inkjets simulation
Inkjets
T-junction device multi-phase flow simulation
Multi-phase Flows
Computational analysis drop formation low viscosity
Precision Droplet Creation

Digital Microfluidics

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

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

접촉선의 이해(Contact Line Insights)

접촉선의 이해(Contact Line Insights)

FLOW-3D는 코팅 성능 향상에 관심이있는 엔지니어에게 이상적인 수치 모델링 기능을 많이 갖추고 있습니다. 전산 시뮬레이션은 코팅 흐름에 영향을 미치는 여러 물리적 과정의 상대적 중요성과 효과를 연구 할 수있는 훌륭한 방법입니다. 물리적인 테스트에서 항상 프로세스를 분리하거나 해당 프로세스의 크기를 임의로 조정할 수있는 것은 아닙니다. 여기에서는 리 볼렛 형성(rivulet formation), 핑거링(fingering), 증발, 거친 표면에서의 접촉선 이동 및 유체 흡수와  관련하여 정적 및 동적 접촉각에 대하여 FLOW-3D의 처리에 대해 설명합니다.

 

정적 및 동적 접촉각(Static and Dynamic Contact Angles)

FLOW-3D는 정적 접촉각의 함수로 동적 접촉각을 정확하게 계산하고 입력으로 설정하며 자유 표면 인터페이스에서 작용하는 관련된 힘을 정확하게 계산하여 유체의 소수성을 캡처 할 수 있습니다. 아래 시뮬레이션은 물방울이 경사를 따라 내려갈 때 정적 접촉각이 동적 접촉각에 미치는 영향을 보여줍니다.

 

흡수(Absorption)

종이 기판에 액 적의 충격 및 흡수는 전산 유체 역학 소프트웨어를 사용하여 연구 할 수 있습니다. 여기서 FLOW-3D는 섬유층에서 물방울 충돌을 시뮬레이션하는데 사용되며 표면 장력, 접촉각 및 점도와 관련된 유체 전면의 전파를 살펴 봅니다.

 

 

아래의 FLOW-3D 시뮬레이션에서, 낙하는 직경이 40 미크론이며 초기 하향 속도는 300 cm / s입니다. 기재는 종이이고, 기공률이 30 % 인 20 미크론 두께입니다.

 

 

액체 필름의 핑거링(Fingering in Liquid Films)

FLOW-3D에서 동적 접촉선은 동적 접촉각이나 접촉선의 위치를 ​​지정할 필요없이 직접 모델링됩니다. 이는 소량의 유체에서 유체에 영향을 미치는 모든 동적 힘을 포함하는 수치 모델을 사용하여 수행됩니다. 정적 접촉각은 액체-고체 접착력을 특성화 하는데 사용됩니다.

액체 시트의 핑거링. 왼쪽은 0 °, 오른쪽은 70 °

여기서, 이러한 접근법의 힘의 적용은 경사 표면 아래로 흐르는 액체 필름에서 관찰 된 핑거링에 의해 제공됩니다. 실험적 관찰에 따르면 두 가지 뚜렷한 핑거링 패턴이 발생합니다. 첫 번째 패턴은 작은 정적 접촉각(즉, 습윤 조건)이며 상하한이 모두 하향으로 움직이는 쐐기형 핑거를 나타냅니다. 두 번째 패턴은 큰 정적 접촉각(즉, 습윤 조건이 열악함)이며 가장 균일한 폭을 가진 긴 핑거이고 가장 큰 한계점은 하향으로 움직이지 않는 것이 특징입니다.

 

 

증발 효과(Evaporative Effects)

퇴적(Deposit)

분산 된 고체 물질을 함유하는 액 적은 고체 표면에서 건조 될 때, 함유하고 있는 고체 물질을 침전물로서 남깁니다. 이 침전물의 형상이 많은 인쇄 공정, 청소 및 코팅 공정에 중요한 영향을 미칩니다. 한 종류의 퇴적물의 전형적인 예는 위의 이미지와 같이 엎질러 진 커피 패치의 둘레를 따라 링 얼룩이 형성되는 “커피 링” 문제입니다. 이 유형의 링 침전물은 액체의 증발로 인한 표면 장력 구동 흐름의 결과로, 특히 낙하 둘레에서 발생합니다.

 

건조(Drying)

FLOW-3D의 증발 잔류 액체 모델은 건조 후 톨루엔으로 형성된 잔류된 물의 3D형상을 시뮬레이션합니다. (30 배 확대)

건조는 코팅 공정의 중요한 부분입니다. 하지만 건조의 결함으로 잘 도포 된 코팅을 완전히 취소 할 수도 있습니다. 건조 중에 온도 및 용질 구배는 밀도 및 표면 장력 구배로 인해 코팅 내 유동을 유도 할 수 있으며, 이는 코팅 품질을 잠재적으로 파괴 할 수 있습니다. FLOW-3D의 증발 잔류 물 모델을 사용하면 건조로 인한 흐름을 시뮬레이션하고 값 비싼 물리적 실험에 소요되는 시간을 줄일 수 있습니다.

 

모델링 링 형성(Modeling Ring Formation)

증발에 의해 접촉 라인에서 생성 된 흐름 시뮬레이션

윗쪽 그림에서 FLOW-3D는 증발이 가장 큰 접촉선에서의 증착으로 인해 에지 피닝(edge pinning)이 발생함을 보여줍니다. 증발은 증발로 인한 열 손실로 인해 액체를 냉각시킵니다 (색상은 온도를 나타냄). 동시에 고체 표면은 전도에 의해 액체를 가열합니다. 접촉선 주변에서 증발이 가장 커서, 액체가 접촉선을 향해 흘러 정적 조건을 재설정합니다. 최종 결과는 액체가 완전히 증발하는 액체 가장자리에 현탁 된 고체의 증착입니다.

 

 

참고
[1] Deegan, R., Bakajin, O., Dupont, T. et al. Capillary flow as the cause of ring stains from dried liquid drops, Nature 389, 827–829 (1997).

 

Microfluidic Circuits

Microfluidic Circuits

생물학에서 물질을 한 장소에서 다른 장소로 운반하거나 수백 개의 검사를 병렬로 수행하기 위해 사용하는 미세 유체 회로 장치 분야에서 최근 발전하고 있습니다. 일반적으로 이러한 회로는 특정 논리(AND, OR, XOR 등) 또는 여러 로직의 조합을 기반으로 합니다. 따라서 이러한 회로를 마이크로유체 논리 회로라고도 합니다. 전자 회로와 유사하게 오일은 채널과 공압 밸브를 통과하며 압력 디퍼렌셜에 의해 구동됩니다(전자 회로의 기존 전위/전압 디퍼렌셜과는 대조적으로). FLOW-3D의 움직이는 물체 모델은 유체 흐름과 결합되어 공압 밸브의 움직임을 시뮬레이션할 수 있습니다.

Simulation of a pneumatic latching valve used in microfluidic demultiplexer. The animation starts at stage 3 – the open stage, and finally evolves to stage 7 – the closed stage.

Read the Microfluidic Circuit – Pneumatic Latching Valve blog.

Magnetic Fields

Magnetic Fields

균일 한 자기장에서 자성 입자는 자화되어 쌍극자-쌍극자 상호 작용으로 인해 사슬 모양의 미세 구조로 조립됩니다. 조립된 체인은 외부 필드의 방향에 맞춰 정렬되는 경향이 있습니다. FLOW-3D를 사용한 이 분석에서는 superparamagnetic beads의 초기 무작위 분포와 베이스에 내장 된 spherical (gold colored) magnetic dipole elements 배열을 포함하는 마이크로 채널을 통해 균일한 필드가 z 방향으로 위쪽으로 적용됩니다. 적용된 필드가있는 경우 비드는 자화되어 개별 체인과 같은 구조로 조립됩니다. 이러한 구조는 고정된 쌍극자 요소에 끌립니다. 분석은 입자 체인의 자체 조립과 내장 된 쌍극자 요소에 체인의 후속 부착을 보여줍니다. 계산 모델은 유체가 입자 운동에 점성 항력을 제공하고 움직이는 입자가 차례로 유체 흐름을 변경하는 완전 결합 입자 유체 상호 작용을 고려합니다. 모델링 결과는 University of Buffalo에서 제공합니다. University at Buffalo의 FLOW-3D 작업에 대한 자세한 내용을 보려면 여기로 이동하십시오.
( In a uniform magnetic field, magnetic particles become magnetized and assemble into chain-like microstructures due to dipole-dipole interactions. The assembled chains tend to align with the direction of the external field. In this analysis using FLOW-3D, a uniform field is applied upward in the z-direction through a micro-channel that contains an initial random distribution of superparamagnetic beads and an array of spherical (gold colored) magnetic dipole elements embedded in its base. In the presence of an applied field, the beads become magnetized and assemble into discrete chain-like structures.  These structures in turn, are attracted to the anchored dipole elements. The analysis shows the self-assembly of particle chains and the subsequent attachment of the chains onto the embedded dipole elements. The computational model takes into account fully-coupled particle fluid interaction where the fluid provides a viscous drag on particle motion and the moving particles, in turn, alter the fluid flow. Modeling results courtesy of the University of Buffalo. Go here for more information about the University at Buffalo’s work with FLOW-3D. )

Lab-on-a-chip

다양한 표면 장력을 사용하는 패턴화된 표면

마이크로 채널의 패턴화된 표면은 액체 사이의 실제 물리적 벽 없이도 여러 액체가 나란히 흐르는 특정 경로를 따라 한 저장소에서 다른 저장소로 액체를 운반하는 데 사용할 수 있습니다. 패턴화된 표면은 랩 온어 칩 (lab-on-a-chip), 바이오어세이, 마이크로 리액터 및 화학적 및 생물학적 감지를 통해 유체를 운반하는 데 사용됩니다. 이 경우 표면 장력은 패턴화된 흐름을 생성하기 위해 마이크로 채널의 유체 흐름을 조작하는데 사용됩니다. 고체 표면에서 유체의 친수성 또는 소수성 거동을 이용하여 마이크로 채널을 통한 여러 유체의 움직임을 제어합니다. 마이크로 채널 내부의 유체 흐름은 층상이므로 여러 유체 흐름 (이 경우 2 개)이 난류 혼합없이 나란히 흐를 수 있습니다. 유체 흐름의 측면에는 물리적 벽이 없기 때문에 흐름은 소위 가상 벽에 의해 제한됩니다. 이 벽은 기본적으로 두 유체 사이의 친수성 경계입니다.

Patterned surfaces in micro channels
Experimental results showing the three phases – A, B and C (left to right), Bin Zhao et al.

위 그림은 마이크로 채널의 실험을 보여줍니다. 중앙 수평 채널의 중간 스트립은 친수성이지만 상부 및 하부 수직 채널과 함께 나머지 채널은 소수성의 정도가 다릅니다. 소수성은 접촉각의 몇도 정도만 다릅니다. 상부 채널의 접촉각은 118o이고 하부 채널의 접촉각은 112o입니다. 그러나 접촉각의 작은 차이는 유체가 이러한 영역으로 흐르기 위해 상당히 다른 압력을 필요로합니다.

Numerical Simulation

처음에는 모든 채널이 다른 유체(투명)로 채워집니다. 분홍색 액체가 수평 채널로 밀리면 중앙 영역(단계 A)의 친수성 경로를 사용합니다. 압력이 증가하면 유체는 하부 친수성-수성 장벽을 깨고 하부 친수성 영역(단계 B)으로 흐르기 시작합니다. 압력을 더 높이면 마침내 유체가 상부 친수성-수소성 장벽을 부수고 상부 영역에서도 흐르기 시작합니다(Phase C).

Numerical results - patterned surfaces using varied surface tension
Numerical results showing the three phases – A, B and C.

위의 수치 결과는 둘 사이에 중요한 차이가 있다는 점을 고려할 때 실험에서 패턴화된 표면 연구의 전반적인 아이디어와 합리적인 비교 가능성을 보여줍니다. 위에 표시된 수치 결과는 과도 상태 (압력이 지속적으로 증가)이므로 유체 경계가 실험 결과와 정확히 유사하지 않습니다. 마찬가지로 유체 특성은 실험에 사용 된 특성과 정확히 유사하지 않습니다. 그럼에도 불구하고 유체 1은 실험에서와 같이 압력이 증가함에 따라 단계 A, B 및 C를 통과합니다. 단계 B에서 투명한 유체는 계속해서 위쪽 채널을 통해 흐르지 만 분홍색 유체만 아래쪽 영역으로 흐릅니다. 이것은 실험과 일치합니다. 흥미로운 것은 C 단계에서 나타난 기포 형성입니다. C 단계에서 기포 형성과 같은 흥미로운 물리학에 대한 계시와 연구는 미세 유체 장치의 설계 및 제작 과정에 중요 할 수 있습니다.

FLOW-3D Results

아래 애니메이션은 위의 실험에 대한 FLOW-3D의 시뮬레이션 결과를 보여줍니다. 유체 1 (하늘색)은 실험의 분홍색 유체와 동일합니다. 처음에는 전체 도메인이 Fluid 2 (투명 유체)로 채워집니다. 압력은 단계적으로 증가하고 시뮬레이션이 진행됨에 따라 세 단계를 모두 볼 수 있습니다.

Evolution of fluid flow with increasing pressure in patterned micro channels created by varying contact angles.

Ref: Bin Zhao, Jeffrey S. Moore, David J. Beebe, Surface-Directed Liquid Flow Inside Microchannels, Science 291, 1023 (2001)

Learn more about the power and versatility of modeling microfluidic applications with 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

Cell Behavior

Cell Behavior

정밀하고 신중하게 제어되는 화학 반응성 구배를 생성 할 수있는 능력은 미세 유체학을 운동성, 화학성 및 소수의 미생물 집단에서 항생제에 대한 내성을 단기간에 진화시키고 개발하는 능력을 연구하는 이상적인 도구가 됩니다. FLOW-3D는 연구자들이 아래 예제에 표시된 것처럼 새롭고 더 나은 gradient generators를 고안하는 데 도움이 될 수 있습니다.

1-D Gradient generator with de-coupled convection and diffusion

FLOW-3D를 사용한 이 1-D 미세유체 팔레트 시뮬레이션에서는 표시된 흐름선을 통해 주 중앙 마이크로 채널에서 대류 셀의 깨끗한 디커플링을 확인할 수 있습니다. 이 흐름은 모두 대류 단위로만 제한되며 마이크로 채널로 유출되는 단 한 개의 흐름도 없어 대류 및 확산의 디커플링이 우수합니다. 소스 농도의 진화는 그림에서 볼 수 있으며, 애니메이션이 끝날 때쯤이면 눈에 띄게 일정해집니다.

This FLOW-3D simulation of a 2-D microfluidic palette demonstrates a spatio-temporal control on the generated gradients. The source and sink are rotated at an angular velocity. Also, after every t seconds, the active access port is deactivated and the next port is turned on. To see the live status of the diffusion inside the chamber, three line probes are placed in the simulation (marked in red, blue and black, respectively, in the bottom right window of the simulation).2-D 마이크로 유체 팔레트의 이  FLOW-3D 시뮬레이션은 생성된 그라데이션에 대한 spatio-temporal 제어를 보여줍니다. 소스 및 sink는 각 속도로 회전합니다. 또한 t초마다 활성 액세스 포트가 비활성화되고 다음 포트가 켜집니다. 챔버 내부의 확산 상태를 확인하기 위해 시뮬레이션에 세 개의 라인 프로브가 배치됩니다(시뮬레이션의 오른쪽 하단 창에 각각 빨간색, 파란색 및 검은색 표시).

Read the Microfluidic Palette – A Gradient Generator blog.