FLOW-3D CAST 사양

FLOW-3D CAST Feature

CAST virtual foundry conference banner

Active Simulation Control

실행중인 해석의 제어 파라미터는 History probes에서 사용자가 정의한 조건에 따라, 런타임 동안에 자동으로 변경 될 수 있습니다. History probes에 의해 기록된 시뮬레이션 변수는 경계 조건, mass source 및 General Moving Object 기능을 이용하여, 시간에 따른 개체의 동작을 제어하기 위해 사용될 수있습니다. 예를 들어, 고압다이캐스팅 해석에서 게이트에 설정한 History probes에 유체가 도달하면, 그 정보를 캡처하는 데이터 출력 주파수를 증가시켜 플런저의 속도를 고속으로 자동 전환 될 수있습니다. 고압다이캐스팅 해석은 유체가 게이트에 도달 할 때 자동으로 고속 전환됩니다. 이 프로세스는 새로운 실행 시뮬레이션 제어 기능을 통해 자동으로 진행됩니다. 저속 구간에서 플런저의 움직임은 trigger 슬리브의 용융물에 혼입되는 공기의 양을 최소화하기 위해 Barkhudarov 방법 1을 사용하여 계산됩니다. 이 결과는 훨씬 더 높은 품질의 주조품이 나올수 있도록 설계하는데 도움이 될 수 있습니다. Read the development note > Read the blog post >

Batch Postprocessing & Report Generation

Batch 후처리 및 보고서 생성은 해석 결과 분석시 사용자의 해석 처리 시간을 절약하기 위해 개발되었습니다. Batch 후처리는, 해석이 완료된 후, 사용자가 애니메이션, 시나리오, 그래프, 텍스트 데이터 시리즈를 정의하여 자동으로 생성되도록 할 수 있습니다. 그래픽 요청은 백그라운드에서 FlowSight를 실행하여 처리되도록 FLOW-3D Cast에 정의되어 있습니다. 원하는 해석 결과를 생성할 수 있는 컨텍스트 파일을 사용하면 Batch 후처리 기능을 사용할 수 있습니다. Batch 후처리가 완료되면, 사용자는 쉽게 자신의 관리자, 동료, 또는 클라이언트에 보낼 수있는 HTML5 형식의 완벽한 기능을 갖춘 보고서를 만들 수 있습니다. 이미지 및 동영상도 보고서에 포함 할 수 있고, 사용자는 텍스트, 캡션, 참고 문헌의 형식을 완벽하게 제어 하고 유지할 수 있습니다. Read the blog post >

Metal Casting Models

Squeeze Pin Model

스퀴즈 핀은 주조시 주입 공급이 어려운 영역에서, 응고하는 동안 금속 수축을 보상하기 위해 사용되는 실제의 다이 캐스팅 머신의 동작을 모델링하는 해석을 할 수 있습니다. 스퀴즈 핀은 선택된 표면에 cylinderical squeeze pin을 추가하여, STL 파일 또는 대화식으로 생성 될 수 있습니다. Read the development note >

Intensification Pressure Model

새로운 플런저 타입 형상이 추가 되었습니다. 강화된 압력 조건으로 macro-shrinkage 와 micro-porosity 제거를 지정할 수 있습니다.

Thermal Die Cycling model

FLOW-3D Cast v4.1's full process thermal die cycling model

다이싸이클링 (Thermal die cycling, TDC) 모델에 새로운 두 가지의 단계가 추가되었습니다. 금형이 열린 상태에서 제품이 여전히 금형 내부에 있는 ejection 단계와, 금형이 닫혔지만 사출 바로전의 preparation 단계가 추가되었습니다. 또한, 마지막 싸이클만이 아닌 모든 금형 싸이클 모두 수렴된 결과를 전달하기 위해 TDC 솔버가 성능 손실 없이 최적화 되었습니다. Read the blog post >

Valves and Vents

Modeling valves and vents in FLOW-3D Cast v4.1

밸브와 밴트의 외부 압력과 온도는 이제 사용자가 다이 캐스팅 공정에서 충진중에 보다 실제적인 동작을 정의 할 수 있도록, 시간의 표 함수로서 정의 할 수있습니다. 밸브 및 벤트의 압력 및 온도는 프로세스 설계 단계에서 유용한 제품 내부에 설정된 프로브에 의해 제어 될 수 있습니다.

PQ2 Diagram

PQ2다이어그램의 사용은 사용자가 더 나은 슬리브의 플런저 실제 움직임과 유사하게 적용 할 수 있습니다. 새로운 기능은 실제 공정 변수가 아직 알려져 있지 않았을 때 다이캐스팅 설계 단계 중에 특히 유용합니다. Read the blog post >

Cooling Channels

냉각 채널은 금형 각각의 냉각 유로에 의해 제거되거나 추가된 열의 총량에 의해 제어 될 수 있습니다. Read the development note >

Air Entrainment Model

Air entrainment 모델에 compressibility를 입력하는 새로운 옵션이 추가되었습니다. 고압 다이캐스팅의 충진 공정과 같은 경우, 공기 압축성은 유체 압력의 변화로 인한 유체의 흐름에 중요한 인자가 됩니다.
 

Cavitation Model

캐비테이션 모델은 유동 조건의 더 넓은 범위에 걸쳐 유체의 캐비테이션 거동을 나타내도록 개선되었습니다. 캐비테이션 생성에 대한 새로운 옵션은 경험적 관계를 기반으로, 기존의 일정한 속도로 생성되는 방식에서 보완되었습니다. 새로운 passive gas model 옵션은 open bubbles이 아닌 유체내에 cavitationg gas를 추적하여, 계산에 필요한 격자와 계산시간을 줄일 수 있습니다. Read the development note >

Two-fluid Phase Change Model

Two-fluid phase change model 은 과냉각을 포함하도록 확장되었습니다. 일정한 과냉각 온도를 정의하고 가스 온도가 응축이 일어나기 전에 포화점 이하로 내려갈 수 있게 함으로써 구현됩니다.

Simulation Results and Analysis

Simulation Results File Editor

사용자가 FLOW-3D Cast v4.1 결과 파일들을 병합 및 제거 할 수 있는 편집 유틸리티

Linking flsgrf.* files

Restart 해석 결과 파일들(flsgrf.*)은 FlowSight 에서 하나의 연속적인 애니메이션 결과를 표시하기 위해 restart source 결과로 링크될 수 있습니다.

Fluid/wall Contact Time

A new spatial quantity has been added to the solution output that stores the time that metal spent in contact with each geometric component, as well as the time spent by each component with metal.

용탕이 각 geometry 컴포넌트를 접촉한 시간과 각 컴포넌트가 용탕과의 접촉 시간을 나타내는 새로운 공간적 양이 해석 아웃풋에 추가 되었습니다.

Performance and Usability

Calculators

열전달 계수, 열 침투 깊이, 밸브 손실 계수, 슬리브에 용탕량(깊이), 플런저의 속도를 계산할 수 있는 Calculators 기능이 Model Setup 창에서 바로 가능해졌습니다. 또한 유틸리티 메뉴에서도 가능합니다.

Thermal Die Cycling

Heat transfer database in FLOW-3D Cast v4.1

열전달 계수 데이터베이스와 각 싸이클 단계들이 입력되어있어 간편하게 다이싸이클링 해석을 하실 수 있습니다.

GMRES Pressure Solver

GMRES pressure solver의 속도가 솔버 데이터 구조의 최적화로 인해 2배까지 향상되었습니다. 이로 인해 메모리 사용량이 20% 미만으로 증가할 수 있습니다. Read the blog post >

Sampling Volumes

Sampling volume 기능은 STL로 정의할 수 있습니다. 각 sampling volume에 의해 계산된 양들의 목록은 유체의 부피, 최대/최소 온도, 파티클의 갯수와 같은 전체 해석 영역에 대해 모두 같은 양이 되도록 확장되었습니다.

 

FSI/TSE Model

구조분석 모델의 성능이 부분적인 coupling으로 해석 솔버의 병렬화와 최적화를 통해 향상되었습니다.

Workspaces

Workspaces 를 이전에 설치된 FLOW-3D에서 가져올 수 있습니다. Workspaces 와 사용자가 선택한 시뮬레이션들을 복사할 수 있습니다.

Expanded Simulation Pre-check

Simulation pre-check 기능은 preprocessor checks를 포함하고, 문제가 발생하는 경우 링크됩니다.

Improved Transparency

Depth-peeling 옵션은 transparent geometries 를 좀 더 잘 표현하고, v4.0보다 10배 빨라졌습니다.

Interactive Tools

Baffles, history probes, void/fluid pointers, valves, mass-momentum sources, squeeze pins에 대한 새로운 대화형 생성 기능이 추가되었습니다. 또한 probing과 clipping 도구들이 대화형으로 개선되었습니다.

General Enable/Disable

모든 objects (e.g., mesh blocks)은 활성화/비활성화 할 수 있습니다.

Estimated Remaining Simulation Time

솔버 메세지 파일에 short-print로 추정된 잔여 해석 시간이 추가 되었습니다.

Tabular Data

테이블 형식의 데이터에서 선택된 데이터를 마우스 오른쪽 버튼을 클릭하여 csv파일 또는 외부 파일에 복사, 저장할 수 있습니다.

1 23-10 Michael R. Barkhudarov, Minimizing Air Entrainment, The Canadian Die Caster, June 2010

[FLOW-3D 이론] Numerical Approximations 수치근사 – 압력 솔루션 알고리즘

Numerical Approximations 수치근사

Pressure Solution Algorithm 압력 해 알고리즘

질량보존의 수치 처리는 압축성과 비압축성에서 상당히 다르다. 그러나 어느 경우든지 적합한 질량 방정식은 셀내의 압력을 결정하고 속도를 갱신하는 알고리즘에 이르게한다. 압축 유동이나 제한적 압축 유동(내재적 압력 해법을 사용해야 하는)에서 연속방정식(10.6) 또는 (10.8)은 셀에서의 압력과 속도의 타원조건으로 직접 해석될 수있다.  압축유동에서는 연속방정식(10.1)이 포물선 방정식으로 즉, 시간에 대해 전진하는 알고리즘에 의해 해석된다. 이 때 압력은 상태방정식 밀도가 갱신되는 셀 밀도와 같게함으로써 결정된다.  이 경우 시간단계 크기가 음파 전달에관한 안정성을 확실하게할 만큼 충분히 작다면 속도는 갱신될 필요가 없다. 압축유동에 대해 내재적 선택을 사용하면 더 큰 시간단계를 허용하나 충격파나 저밀도 파형에는 덜 정확할 수 있다.

때때로 인위적 제한 압축성을 유체에 추가하는 것이 해에 상당한 에러를 일으키지 않고 수렴을 증진시킬 수 있다. FLOW-3D에서 이는IMP = 2(디폴트는 1)로 지정함으로써 자동적으로 이루어진다. 상세한 모델 내용을 위해 https://users.flow3d.com/tech-notes/default.asp의 사용자 주소에서Flow Science Technical Note #55를 참조하라.

Incompressible SOR Method  비압축 SOR방식

식(10.315)으로부터 계산된 속도는 제한적 압축 연속 방정식(10.8)에 대한 다음의 이산화 근사식을 만족시켜야하며,

(10.352)

여기서XCi는 셀i의 중심의x위치이다. 원통좌표계에서는CYL = 1.0 및 Ri = XIM1/XCi 이며XIM1는 망내 마지막 실재 셀의 바깥 가장자리의 반경(x위치)이다. 데카르트 형상에서는 모든i 에 대해Ri = 1.0 이고CYL = 0.0이다. 항RSOR 은 셀내의 유체 체적 소스를 뜻한다. 압축성계수1/(ρc)2는 2유체 문제의 공식에 의해 계산된다.

 (10.353)

1유체 문제에서는 단지 두번째 항만 존재한다. 첨자l v는 각기 유체 1과 2를 뜻한다. 제한적 압축성은 두가지 목적, 물리적및 수치적, 으로 사용될 수 있다. 이 두 입력변수RCSQL = 1/(︀ρc2)︀ 와 RCSQV = 1/(︀ρc2)︀를 지정함으로써 유체의 제한적인 물리적 압축성을 모델할 수 있다. 또한 RCSQL 와 RCSQV를 적절히 지정함으로써 물체 경계상에서의 자유표면의 붕괴에 의해 종종 발생하는 수치적 압력 파동의 효과를 완화시킬 수 있다.

속도가 식 (10.352)을 만족시키기 위해 압력 그러므로 유체가 차지하고 있는 셀 내의 속도를 조절하는 것이 필요하다. 이는 둘 중 하나의 방식으로 행해진다. 가장 간단한 방법은 successive over-relaxation(연속가속완화) (SOR)반복 과정이다. 계산망을 망내의 첫번쩨 비 경계셀에서 시작하여 하나씩 쓸어나간다. 쓸림은 먼저i에 대해 시행되고 다음에j그리고 마지막으로k값에 대해서 되어진다. 계산은 단지 유체를 포함는 유체가 빈 이웃이 없는 셀들에 대해 실행된다.  셀(i, j, k)에서의 속도가 식(10.352)을 만족시키기 위해 필요한 압력변화는 아래와 같으며, 여기서 S는 식(10.352)의 좌측이다.

(10.354)

식(10.354)은 단순히S = 0를 이루기 위해 필요한p값을 생성하는 완화과정의Newton형태이다. 각 셀에서S를 평가하기 위해 사용된 속도 값은 반복과정에서 사용 가능한 가장 최신의 값이다.  식(10.354)으로부터의 결과를 이용하여 셀압력의 세 추정치는

(10.355)

이며, 셀의 면들에 위치한 속도들의 새 추정치는 다음 식과 같다.

      (10.355)

여기서 여기에 나타나는 속도들 또한 반복중에 가장 최신의 값들이다. 반복과정을 시작하기 위해 식(10.315)으로부터의 새 추정 속도들은 전 시간단계로부터 남아있는 압력값과 함께 이용된다. 물론 면적이0인 곳에서의 속도는 이 단계에 수정되지 않는다.

자유표면을 가지는 셀들에서, 즉 유체가있지만 하나 또는 더 많은 빈 이웃 셀들을 가지는 셀에서는 다른 과정이 이용된다. 이러한 셀들에서 요구되는 경계조건은 압력이 표면에서 지정된 값, 즉 표면에서의 ps 이다. 표면압력은 이웃한 void영역 압력, PR, 과 표면장력 압력, PS,의 합과 같도록 지정되며

(10.357)

여기서n은 인접한 빈 공간의 색인이다. PS의 평가는Surface Tension with Wall Adhesion에서 기술된다. 표면압력, ps 은 셀내의 정압분포를 가정하여 표면셀의 중심에서의 압력pi,j,k 으로 외삽하여 해석에 적용된다.정압변화는NF에의해 정의된 바와같은 표면에 수직한 방향에서의 순수가속도에 의존한다. 이 표면 셀 압력은 압력 반복동안에 변하지않으며 고정 경계값으로 취급된다. 이런 방법으로 셀내의 자유 표면의 실제 위치가 확실하게 고려된다.

고압, 단열 기포가 존재할 때[Hir92], [BC94] 수치 불안정성이 공간지역 압력의 외재적 근사로 인해 발생할 수 있다고 알려져 있다(단열 기포모델은Variable Pressure (Adiabatic) Void Region에서 기술되어 있다). 이러한 어려움을 없애기 위해FLOW-3D에 내재적 기포모델이 추가되어 있다. 이의 목적은 기포압력 변화가 한 사이클의 마지막에서 계산되도록 기대하고 이를 통상적인 압력-속도 반복과정에 포함하는 것이다. 내재적 기포 모델은 기포의 “강성도”가 너무 크지 않다면 잘 작동한다. 다른 말로, 강성 기포는 기대되고 그리고 실제압력 변화가 한시간 단계내에서 너무 다른 기포를 뜻한다. 이런 강성 기포들이 발생하면 해석은 매 사이클 마다 큰 압력변화와속도와 다른 양 들에서 이에 상응하는 커다란 변동을 가질 수 있다(상세 내용을 위해 Ref. [Hir92] 참조하라).

2 유체문제에서 모든 셀들은 유체로 가득 차 있다고 간주된다; 즉, 압력과 속도는 반복하는 동안에 모든 유체 셀 내에서 조절된다. 표면장력이 작용되면 압력PSi,j, 은 두 유체중에 하나에만 작용해서 표면장력으로 의한 경계면을 통과하는 압력에서의 불연속성이 유지된다.

완전한 반복은 식(10.354), (10.355) 및 (10.356)에따라 모든 유체가 가득찬 셀내의 압력및 속도를 조절하는 것으로 이루어져 있다. 반복시의 수렴은 모든 셀들이 어떤 일정 작은 수인EPSI ·VFi,j,k 보다 작은S값을 가질 때 이루어진다.

EPSI의 값은 자동적으로FLOW-3D에 의해 각 시간 사이클에서 시간 단계 크기의 함수로 계산된다. 이 알고리즘은 입력변수EPSADJ의 값이 양수이면 원용된다. 선택적으로 한 EPSI의 상수가 한 계산 과정에서 사용 가능하다.

어떤 경우에는 반복의 수렴이 식(10.355)에서의δp를 완화인자OMEGA로 곱함으로써 가속화될 수 있다.

OMEGA 는 1.7 또는 1.8이 최적값이다. 어떤경우에도 는 2.0이 넘어서는 안되는데 이는 이럴 경우 불안정한 반복이 발생하기 때문이다. 압축성유동에서OMEGA는 1.0으로 지정된다. 또한 시간 단계 크기가 상당히 대류 안정성 한계보다 작을 때 비압축성 유동에 대해OMEGA는 1.0으로 사용하는것이 권장된다. 이는 해에서의 잠재적 압력 잡음을 감소시킬 수 있다.

Incompressible Line Implicit SADI Method 비압축성 선 내재적 SADI 방식

앞에 언급된 압력을 계산하기 위한SOR반복법은 간단하고 많은 문제에서 잘 작동한다. 그러나SOR방식의 수렴이 상당히 느려지는 경우들이 있다. 예를들면, 한방향으로의 셀 크기가 다른 방향으로 보다 훨씬 큰 망은SOR압력 완화가 횡방향의 작은 셀크기에 의해 제한되므로 큰 셀방향으로 느린 수렴성을 보여줄 것이다.

이런 형태의 더딘 수렴에 대한 보완은 더 작은 셀크기의 방향에서 더 내재적인 해석 방식을 사용하는 것이다. 이런 목적으로 수정된Alternating-Direction-Implicit (SADI)방식이 개발되었다. SADI는 망 셀의 한 i, j, 또는 k열을 따라 압력에 대한 표준 3중 대각해법에 근거한다.   이 해법은 주기적 경계를 포함하는 모든 경계조건에 적용 가능하다.

주기적 경계가 원통좌표계에서 방위각의 방향에 사용될 때IADIY=1로 지정함으로써-단지 이 방향으로만- ADI압력 해법을 사용하는 것이 권장된다. 이는 가끔 발생할 수 있는 압력과 속도해에서의 수치 잡음을 제거하는데 도움이 될 것이다. 다른 방법으로는OMEGA= 1.0을 지정하여 상향 완화를 잠금으로써 잡음을 감소시키는데 도움이될 수 있다.

SADI방식이 z방향으로사용되면 반복은 모든 i j색인을 거치고 각(i, j)쌍에 대해k-색인 방향에서 압력에 대해 내재적으로 해석하는 것으로 이루어진다. 이웃 열들에서 필요한 압력 값들은 표준ADI에서는 반드시 항상 그렇게 실행되지는 않지만 최신의 반복값을 항상 취하는데 이는 수렴을 향상시킨다.

SADI방식은 방향의 어떤 조합으로도 사용될 수있다: 어느 하나 또는 둘 또는 셋 모두. 이는 더 비용이 드는 내재적 소해가 단지 전체 반복의 수렴을 향상시키는데 필요한 방향에서만으로 제한될 수 있다는 것을 뜻한다.

단지 한 또는 두 방향으로만 내재적으로 처리될 때 열간의 상향 완화는SOR상향 완화에서 사용되는 변수OMEGA에 의해 조절된다. 일반적으로SADI는 이 변수에 그렇게 민감하지 않으며 디폴트 값OMEGA=1.7은 보통 만족스럽다.  SADI 가 세 방향 모두에 사용될 때 이 경우 상향 완화가 최대값에서 고정되므로 변수OMEGA 는 영향을 미치지 않는다.

Compressible Solution Method 압축성 해 방식

압축성 유동에서 셀 압력은 연속방정식 밀도를 상태방정식으로부터 결정되는 밀도와 동일시함으로써 결정된다. 이 방정식에는 SOR 와 SADI 방식 둘 다 이용 가능하다. 두 알고리즘에서 식(10.354)에서 사용되는 반복 함수는 아래와 같이 정의된다.

  (10.358)

인자∇ · (uAΦ)는 시간n+1에서의 속도에 의거하여 셀의 압축및 팽창을 조절하며 여기서

 (10.359)

C1은 유체 1에서의 음속이다. 외재적 해석 방식에서 외재적 모멘텀 방정식으로부터의 속도는 이며, 이는 가 셀압력에 대한 반복에 독립적이라는 것에 주목한다. 상태 방정식밀도는

 (10.360)

로 정의된다.

상태방정식 밀도는 압축이나 팽창에 의한 시간 정도 n에서 n+1사이에 발생하는 에너지변화에 대해 조절되지 않는 것에 주목한다. 반복이 매시간 단계에서 되풀이 되므로 에너지를 갱신하지 않는데서 비롯되는 오차가 기껏해야 한시간 단계 쳐지고 시간단계 안정성을 이루는 목적을 위해서도 의미가 없다. 식(10.354)에서 사용된 양δS/δp은 각 사이클에서 한 번씩 평가되고 저장되는 항DSDPU = δDTi,j,k/δpi,j,k을 필요로 한다.

SADI 해법은 비압축성에서와 같이 진행한다. 반복함수는 아직 3중 대각 시스템(횡방향에서의p대한 최신의 반복값을 유지하며)으로 처리될 수 있다. 모든 셀에 대해 어느 경우에도 수렵에 도달한다.

  (10.361)

SOR해석 알고리즘에서 상태방정식에서의 밀도에 대한 압력의 의존도가 비선형(코드의 사용자 수정에 의해 가능)일 경우에 유용한 선택이 주어진다. 다음셀로 전진하기 전에 식(10.361)의 조건을 만족하기 위해 셀내의 압력을 변화시키도록내부 반복이 실행된다. 사용자는 내부 반복의 최대수(IITMX)와  완화인자 OMEGA를 지정해도 된다. IITMX > 1이면OMEGA = 1.0이 권장된다.

지정된 속도 및 지정된 압력경계 조건은 계산영역 내에 이용 가능한 반복법중에 어느 방법으로도 계산하기 어려운 균일한 압력변화를 형성할 수가 있다.  이런 경우에 전반적으로 균일한 압력조절을 주기 위해 추가 알고리즘이 압축 해 과정에 주어진다. 이 선택은 변수IPUN를 1로 지정함으로써 활성화된다. If IPUN = 0이면 균일 압력변화는 평가되지 않는다. 이 균일 압력조절은 셀들에 대해 Si,j,k와 셀 체적의 곱을 합하고 이 결과를δS/δp와 셀 체적의 곱의 셀들에 대한 합으로 나눔으로써 계산된다.

GMRES Pressure-Velocity Solvers GMRES 압력-속도 해법

새로운 압력-속도해법이FLOW-3D [AMS90], [BBC+94] [Saa96] 에서 실행되고 있다. GMRES는  일반화된 최소 잔류 방법을 뜻한다. GMRES솔버에 추가하여 또 새로운 선택적 알고리즘- 일반화된 짝 구배(GCG)알고리즘-이 새 GMRES솔버에서 점성항을 위해 실행되고 있다. 이 새 솔버는 많은 범주의 문제들에 대해 아주 정확하고 효과적이다. 좋은 수렴성, 대칭 및 속도성을 갖는다; 그러나SOR 이나 SADI방법보다 더 많은 메모리를 사용한다. GMRES 솔버는  과소 또는 상향 완화를 사용하지 않는다.

사용된 방법들에 대한 상세한 내용은http://users.flow3d.com/tech-notes/default.asp에 있는 사용자 주소상의 Flow  Science TN68에서 찾을 수 있다.

Improved Generalized Minimal Residual (GMRES) Solver in FLOW-3D—How it works and when to use it

Overview
In FLOW-3D1 there are three linear solvers (as of Version 9.2) with which to solve for pressure from the linear system of continuity equations throughout the domain: the successive over-relaxation (SOR) algorithm, the alternating direction implicit (ADI) algorithm and the generalized minimal residual (GMRES) solver. With GMRES, the system of equations is solved simultaneously throughout the domain by an iterative technique. This is quite unlike the approach used in the SOR algorithm, which adjusts the pressure on a cell-by-cell basis to enforce the continuity equation, or the ADI algorithm, which adjusts the pressures along each mesh column, whose direction alternates through the chosen directions. Although the SOR algorithm is very simple and memory efficient—no extra information from each cell needs to be stored—it can require a large number of iterations to converge, especially for problems where nearly uniform pressure adjustments must be made over a large region of the domain. However, it will always eventually converge so long as the relaxation parameter Ω is less than or equal to 1.
With GMRES, the number of iterations required for convergence is typically much smaller (less than 10) than for SOR; however, it does not always converge. Therefore, for many problems the GMRES algorithm is computationally much more efficient because it is able to converge with far fewer iterations than other solver schemes. This is especially true for problems where the pressures over a large region of the domain are intimately coupled; examples are incompressible flow through a piping network, compressing a gas within a storage tank and most confined flow problems. SOR still may be more computationally efficient for fairly shallow free-surface problems where the pressure is more or less controlled by the location of the free surface—the pressure between two points in the liquid are not strongly coupled.

Development of New Pressure-Velocity Solvers in FLOW-3D [FLOW-3D의 새로운 압력-속도 해법의 개발]

1 Introduction
The purpose of this note is to document the development of new pressure-velocity solvers in FLOW-3Dr. In the following section, new solvers are described ¯rst followed by a section of sample problems where typical simulations were performed to illustrate the application of these new solvers. Two appendices are added to describe in detail the general minimum residual (GMRES) and generalized conjugate gradient (GCG) algorithms used in the new solvers.

FLOW-3D/MP Features List

FLOW-3D/MP Features

FLOW-3D/MP v6.1 은 FLOW-3D v11.1 솔버에 기초하여 물리 모델, 특징 및 그래픽 사용자 인터페이스가 동일합니다. FLOW-3D v11.1의 새로운 기능은 아래 파란색으로 표시되어 있으며 FLOW-3D/MP v6.1 에서 사용할 수 있습니다. 새로운 개발 기능에 대한 자세한 설명은 FLOW-3D v11.1에서 새로운 기능을 참조하십시오.

Meshing & Geometry

  • Structured finite difference/control volume meshes for fluid and thermal solutions
  • Finite element meshes in Cartesian and cylindrical coordinates for structural analysis
  • Multi-Block gridding with nested, linked, partially overlapping and conforming mesh blocks
  • Fractional areas/volumes (FAVOR™) for efficient & accurate geometry definition
  • Mesh quality checking
  • Basic Solids Modeler
  • Import CAD data
  • Import/export finite element meshes via Exodus-II file format
  • Grid & geometry independence
  • Cartesian or cylindrical coordinates
Flow Type Options
  • Internal, external & free-surface flows
  • 3D, 2D & 1D problems
  • Transient flows
  • Inviscid, viscous laminar & turbulent flows
  • Hybrid shallow water/3D flows
  • Non-inertial reference frame motion
  • Multiple scalar species
  • Two-phase flows
  • Heat transfer with phase change
  • Saturated & unsaturated porous media
Physical Modeling Options
  • Fluid structure interaction
  • Thermally-induced stresses
  • Plastic deformation of solids
  • Granular flow
  • Moisture drying
  • Solid solute dissolution
  • Sediment transport and scour
  • Cavitation (potential, passive tracking, active tracking)
  • Phase change (liquid-vapor, liquid-solid)
  • Surface tension
  • Thermocapillary effects
  • Wall adhesion
  • Wall roughness
  • Vapor & gas bubbles
  • Solidification & melting
  • Mass/momentum/energy sources
  • Shear, density & temperature-dependent viscosity
  • Thixotropic viscosity
  • Visco-elastic-plastic fluids
  • Elastic membranes & walls
  • Evaporation residue
  • Electro-mechanical effects
  • Dielectric phenomena
  • Electro-osmosis
  • Electrostatic particles
  • Joule heating
  • Air entrainment
  • Molecular & turbulent diffusion
  • Temperature-dependent material properties
  • Spray cooling
Flow Definition Options
  • General boundary conditions
    • Symmetry
    • Rigid and flexible walls
    • Continuative
    • Periodic
    • Specified pressure
    • Specified velocity
    • Outflow
    • Grid overlay
    • Hydrostatic pressure
    • Volume flow rate
    • Non-linear periodic and solitary surface waves
    • Rating curve and natural hydraulics
    • Wave absorbing layer
  • Restart from previous simulation
  • Continuation of a simulation
  • Overlay boundary conditions
  • Change mesh and modeling options
  • Change model parameters
Thermal Modeling Options
  • Natural convection
  • Forced convection
  • Conduction in fluid & solid
  • Fluid-solid heat transfer
  • Distributed energy sources/sinks in fluids and solids
  • Radiation
  • Viscous heating
  • Orthotropic thermal conductivity
  • Thermally-induced stresses
Turbulence Models
  • RNG model
  • Two-equation k-epsilon model
  • Two-equation k-omega model
  • Large eddy simulation
Metal Casting Models
  • Thermal stress & deformations
  • Iron solidification
  • Sand core blowing
  • Sand core drying
  • Permeable molds
  • Solidification & melting
  • Solidification shrinkage with interdendritic feeding
  • Micro & macro porosity
  • Binary alloy segregation
  • Thermal die cycling
  • Surface oxide defects
  • Cavitation potential
  • Lost-foam casting
  • Semi-solid material
  • Core gas generation
  • Back pressure & vents
  • Shot sleeves
  • PQ2 diagram
  • Squeeze pins
  • Filters
  • Air entrainment
  • Temperature-dependent material properties
  • Cooling channels
  • Fluid/wall contact time
Numerical Modeling Options
  • TruVOF Volume-of-Fluid (VOF) method for fluid interfaces
  • First and second order advection
  • Sharp and diffuse interface tracking
  • Implicit & explicit numerical methods
  • GMRES, point and line relaxation pressure solvers
  • User-defined variables, subroutines & output
  • Utilities for runtime interaction during execution
Fluid Modeling Options
  • One incompressible fluid – confined or with free surfaces
  • Two incompressible fluids – miscible or with sharp interfaces
  • Compressible fluid – subsonic, transonic, supersonic
  • Stratified fluid
  • Acoustic phenomena
  • Mass particles with variable density or diameter
Shallow Flow Models
  • General topography
  • Raster data interface
  • Subcomponent-specific surface roughness
  • Wind shear
  • Ground roughness effects
  • Laminar & turbulent flow
  • Sediment transport and scour
  • Surface tension
  • Heat transfer
  • Wetting & drying
Advanced Physical Models
  • General Moving Object model with 6 DOF–prescribed and fully-coupled motion
  • Rotating/spinning objects
  • Collision model
  • Tethered moving objects (springs, ropes, mooring lines)
  • Flexing membranes and walls
  • Porosity
  • Finite element based elastic-plastic deformation
  • Finite element based thermal stress evolution due to thermal changes in a solidifying fluid
  • Combusting solid components
Chemistry Models
  • Stiff equation solver for chemical rate equations
  • Stationary or advected species
Porous Media Models
  • Saturated and unsaturated flow
  • Variable porosity
  • Directional porosity
  • General flow losses (linear & quadratic)
  • Capillary pressure
  • Heat transfer in porous media
  • Van Genunchten model for unsaturated flow
Discrete Particle Models
  • Massless marker particles
  • Mass particles of variable size/mass
  • Linear & quadratic fluid-dynamic drag
  • Monte-Carlo diffusion
  • Particle-Fluid momentum coupling
  • Coefficient of restitution or sticky particles
  • Point or volumetric particle sources
  • Charged particles
  • Probe particles
Two-Phase & Two-Component Models
  • Liquid/liquid & gas/liquid interfaces
  • Variable density mixtures
  • Compressible fluid with a dispersed incompressible component
  • Drift flux
  • Two-component, vapor/non-condensable gases
  • Phase transformations for gas-liquid & liquid-solid
  • Adiabatic bubbles
  • Bubbles with phase change
  • Continuum fluid with discrete particles
  • Scalar transport
  • Homogeneous bubbles
  • Super-cooling
Coupling with Other Programs
  • Geometry input from Stereolithography (STL) files – binary or ASCII
  • Direct interfaces with EnSight®, FieldView® & Tecplot® visualization software
  • Finite element solution import/export via Exodus-II file format
  • PLOT3D output
  • Neutral file output
  • Extensive customization possibilities
  • Solid Properties Materials Database
Data Processing Options
  • State-of-the-art post-processing tool, FlowSight™
  • Batch post-processing
  • Report generation
  • Automatic or custom results analysis
  • High-quality OpenGL-based graphics
  • Color or B/W vector, contour, 3D surface & particle plots
  • Moving and stationary probes
  • Measurement baffles
  • Arbitrary sampling volumes
  • Force & moment output
  • Animation output
  • PostScript, JPEG & Bitmap output
  • Streamlines
  • Flow tracers
User Conveniences
  • Active simulation control (based on measurement of probes)
  • Mesh generators
  • Mesh quality checking
  • Tabular time-dependent input using external files
  • Automatic time-step control for accuracy & stability
  • Automatic convergence control
  • Mentor help to optimize efficiency
  • Change simulation parameters while solver runs
  • Launch and manage multiple simulations
  • Automatic simulation termination based on user-defined criteria
  • Run simulation on remote servers using remote solving
Multi-Processor Computing

FLOW-3D/MP Software Overview

FLOW-3D/MP Overview

FLOW-3D/MP 는 엔지니어가 계산할 도메인이 매우 크거나 시뮬레이션 실행 시간이 너무 많이 소요될 것 같은 문제를 해결 하기 위해 고성능 컴퓨팅 클러스터에서 실행되도록 설계한 FLOW-3D의 분산 메모리 버전입니다. FLOW-3D/MP는 클러스터의 컴퓨팅 노드에서 여러 CPU 코어에 계산 속도를 높이기 위해 병렬화하는 하이브리드 MPI-OpenMP의 방법을 사용합니다. 시뮬레이션 도메인에 따라서 그들 사이의 연산 작업을 분할 클러스터의 연산 노드에 분산된 여러 서브 도메인으로 분할됩니다. 다른 서브 도메인의 결과의 동기화는 메시지 교환 인터페이스 (MPI) 라이브러리를 통 이용하여 노드 사이에서 데이터를 교환함으로써 수행됩니다. 각각의 하위 도메인 내에서의 OpenMP 스레드는 계산을 더욱 병렬화하게 됩니다. 솔버의 성능을 강화하는 MPI와 OpenMP 병렬 처리 결과의 조합은 매우 오래 걸리는 시뮬레이션의 실행 시간을 줄이는 효과가 큽니다.

Why use FLOW-3D/MP?

현재 하드웨어가 멀티 코어, 멀티 CPU 노드 (즉 ccNUMA 공유 메모리)인 고성능 컴퓨팅 (HPC)인 경우 구성은 Infiniband와 같은 고속 네트워크 인프라 스트럭처를 통해 연결됩니다.
더 좋은 연산 성능과 효율의 장점, 전력소비 절감과 비용 감소 및 우수한 유연성을 위해 멀티코어 클러스터 시스템은 과학분야와 같은 고성능 컴퓨팅이 필요한 분야에서 널리사용되고 있습니다.

사용자 지원을 강화하고 멀티 코어 클러스터 솔루션의 정확성을 향상시키기 위해 그리드 해상도를 높이는 등 더 많은 기능을 강화시키고 있습니다.
FLOW-3D/MP는 설계 및 솔루션 정확도를 유지하고 실행시간을 크게 감소시키는 등 클러스터 시스템에서 최고의 기능을 발휘할 수 있도록 최적화되었습니다.
마지막으로, 독립형 스테이션의 메모리 제한은 FLOW-3D / MP의 분산 메모리 접근방식으로 해결 될 수 있습니다.

What kind of performance can I expect?

물론, FLOW-3D/MP의 실제 성능은 시뮬레이션에 따라 다르지만, 솔버는 금속 주조, 물, 환경, 미세 유체 및 항공 우주 등 다양한 애플리케이션을 위해 512 코어까지 확장하여 보여 주었습니다. 여러가지 경우에 대한 성능 그래프와 함께 세부 사항은 벤치 마크의 페이지에 제시했습니다.

How to use FLOW-3D/MP?

FLOW-3D/MP는 일반적으로 클러스터 컴퓨터에 설치되고 실행됩니다. 클러스터 계산은 슈퍼 컴퓨팅 시설의 독립 실행 형 클러스터 또는 일부가 될 수 있습니다. FLOW-3D/ MP와 함께 제공되는 그래픽 사용자 인터페이스는 사용자가 쉽게 설정하고 시뮬레이션을 실행할 수 있습니다. PBS, Torque, SGE와 같은 작업 스케줄러를 사용하여 실행되는 대규모 클러스터 시뮬레이션의 경우, 사용자는 highly configurable 및 독립적인 작업 스케줄러 제출 유틸리티에 액세스 할 수 있습니다.

What’s in FLOW-3D/MP v6.1?

FLOW-3D/MP V6.1은 FLOW-3D v11.1을 기반으로합니다. 일부 주요 기능은 새로운 입자 모델(particle model), 스퀴즈 핀 모델( squeeze pins model), 계류 라인(mooring lines)과 활성화된 시뮬레이션 제어를 포함합니다. 모든 모델은 FLOW-3D/MP에 대해 하이브리드 MPI-OpenMP의 방법론과 호환됩니다. FLOW-3D / MP의 계산 부하 균형은 매우 중요하고 솔버의 성능에 크게 영향을 미칩니다. 로드 밸런싱은 정적(시뮬레이션이 시작되기 전) 및 동적(시뮬레이션 진행중)으로 분류 될 수있습니다.
정적로드 밸런싱을 달성하기 위해, FLOW-3D/MP는 여러 하위 도메인 (MPI 도메인)에 있는 하나의 계산 도메인을 분할하는 자동 분해 도구를 제공하여, 그들 사이를 균등하게 활성화된 cells을 배분합니다. 서브 도메인 사이의 동기화 시간을 최소화하는 것은 성능을 향상시킵니다.
V6.1에서 상기 분해 단계는 사용자의 경험을 반영하여 셋업에서 중단을 피하기 위해 해석 단계와 결합되었습니다. 동적로드 균형을 달성하기 위해, 동적 스레드 밸런싱 기능은 시뮬레이션 과정 동안의 OpenMP 스레드를 조정하는데 사용될 수있습니다. one-fluid에서, 자유 표면 시뮬레이션은 최대 20 %의 성능 향상이 이 기능을 사용하여 달성되었습니다.
V6.1의 다른 중요한 개선은 복잡한 지형 모델링 홍수 이벤트에 사용하여 GMRES 압력 솔버, 일괄 처리 및 보고서 생성의 최적화, 래스터 데이터 인터페이스를 포함합니다. 새로운 모델과 기능에 대한 자세한 내용은 FLOW-3D의 v11.1 페이지를 참조하십시오. * 성능 메트릭은 시뮬레이션은 24 시간에서 실행될 수있는 횟수로 정의된다. 높은 막대가 더 나은 성능을 나타냅니다.

FLOW-3D Features

The features in blue are newly-released in FLOW-3D v12.0.

Meshing & Geometry

  • Structured finite difference/control volume meshes for fluid and thermal solutions
  • Finite element meshes in Cartesian and cylindrical coordinates for structural analysis
  • Multi-Block gridding with nested, linked, partially overlapping and conforming mesh blocks
  • Conforming meshes extended to arbitrary shapes
  • Fractional areas/volumes (FAVOR™) for efficient & accurate geometry definition
  • Closing gaps in geometry
  • Mesh quality checking
  • Basic Solids Modeler
  • Import CAD data
  • Import/export finite element meshes via Exodus-II file format
  • Grid & geometry independence
  • Cartesian or cylindrical coordinates

Flow Type Options

  • Internal, external & free-surface flows
  • 3D, 2D & 1D problems
  • Transient flows
  • Inviscid, viscous laminar & turbulent flows
  • Hybrid shallow water/3D flows
  • Non-inertial reference frame motion
  • Multiple scalar species
  • Two-phase flows
  • Heat transfer with phase change
  • Saturated & unsaturated porous media

Physical Modeling Options

  • Fluid structure interaction
  • Thermally-induced stresses
  • Plastic deformation of solids
  • Granular flow
  • Moisture drying
  • Solid solute dissolution
  • Sediment transport and scour
  • Sludge settling
  • Cavitation (potential, passive tracking, active tracking)
  • Phase change (liquid-vapor, liquid-solid)
  • Surface tension
  • Thermocapillary effects
  • Wall adhesion
  • Wall roughness
  • Vapor & gas bubbles
  • Solidification & melting
  • Mass/momentum/energy sources
  • Shear, density & temperature-dependent viscosity
  • Thixotropic viscosity
  • Visco-elastic-plastic fluids
  • Elastic membranes & walls
  • Evaporation residue
  • Electro-mechanical effects
  • Dielectric phenomena
  • Electro-osmosis
  • Electrostatic particles
  • Joule heating
  • Air entrainment
  • Molecular & turbulent diffusion
  • Temperature-dependent material properties
  • Spray cooling

Flow Definition Options

  • General boundary conditions
    • Symmetry
    • Rigid and flexible walls
    • Continuative
    • Periodic
    • Specified pressure
    • Specified velocity
    • Outflow
    • Outflow pressure
    • Outflow boundaries with wave absorbing layers
    • Grid overlay
    • Hydrostatic pressure
    • Volume flow rate
    • Non-linear periodic and solitary surface waves
    • Rating curve and natural hydraulics
    • Wave absorbing layer
  • Restart from previous simulation
  • Continuation of a simulation
  • Overlay boundary conditions
  • Change mesh and modeling options
  • Change model parameters

Thermal Modeling Options

  • Natural convection
  • Forced convection
  • Conduction in fluid & solid
  • Fluid-solid heat transfer
  • Distributed energy sources/sinks in fluids and solids
  • Radiation
  • Viscous heating
  • Orthotropic thermal conductivity
  • Thermally-induced stresses

Numerical Modeling Options

  • TruVOF Volume-of-Fluid (VOF) method for fluid interfaces
  • Steady state accelerator for free-surface flows
  • First and second order advection
  • Sharp and diffuse interface tracking
  • Implicit & explicit numerical methods
  • Immersed boundary method
  • GMRES, point and line relaxation pressure solvers
  • User-defined variables, subroutines & output
  • Utilities for runtime interaction during execution

Fluid Modeling Options

  • One incompressible fluid – confined or with free surfaces
  • Two incompressible fluids – miscible or with sharp interfaces
  • Compressible fluid – subsonic, transonic, supersonic
  • Stratified fluid
  • Acoustic phenomena
  • Mass particles with variable density or diameter

Shallow Flow Models

  • General topography
  • Raster data interface
  • Subcomponent-specific surface roughness
  • Wind shear
  • Ground roughness effects
  • Manning’s roughness
  • Laminar & turbulent flow
  • Sediment transport and scour
  • Surface tension
  • Heat transfer
  • Wetting & drying

Turbulence Models

  • RNG model
  • Two-equation k-epsilon model
  • Two-equation k-omega model
  • Large eddy simulation

Advanced Physical Models

  • General Moving Object model with 6 DOF–prescribed and fully-coupled motion
  • Rotating/spinning objects
  • Collision model
  • Tethered moving objects (springs, ropes, breaking mooring lines)
  • Flexing membranes and walls
  • Porosity
  • Finite element based elastic-plastic deformation
  • Finite element based thermal stress evolution due to thermal changes in a solidifying fluid
  • Combusting solid components

Chemistry Models

  • Stiff equation solver for chemical rate equations
  • Stationary or advected species

Porous Media Models

  • Saturated and unsaturated flow
  • Variable porosity
  • Directional porosity
  • General flow losses (linear & quadratic)
  • Capillary pressure
  • Heat transfer in porous media
  • Van Genunchten model for unsaturated flow

Discrete Particle Models

  • Massless marker particles
  • Multi-species material particles of variable size and mass
  • Solid, fluid, gas particles
  • Void particles tracking collapsed void regions
  • Non-linear fluid-dynamic drag
  • Added mass effects
  • Monte-Carlo diffusion
  • Particle-fluid momentum coupling
  • Coefficient of restitution or sticky particles
  • Point or volumetric particle sources
  • Initial particle blocks
  • Heat transfer with fluid
  • Evaporation and condensation
  • Solidification and melting
  • Coulomb and dielectric forces
  • Probe particles

Two-Phase & Two-Component Models

  • Liquid/liquid & gas/liquid interfaces
  • Variable density mixtures
  • Compressible fluid with a dispersed incompressible component
  • Drift flux with dynamic droplet size
  • Two-component, vapor/non-condensable gases
  • Phase transformations for gas-liquid & liquid-solid
  • Adiabatic bubbles
  • Bubbles with phase change
  • Continuum fluid with discrete particles
  • Scalar transport
  • Homogeneous bubbles
  • Super-cooling
  • Two-field temperature

Coupling with Other Programs

  • Geometry input from Stereolithography (STL) files – binary or ASCII
  • Direct interfaces with EnSight®, FieldView® & Tecplot® visualization software
  • Finite element solution import/export via Exodus-II file format
  • PLOT3D output
  • Neutral file output
  • Extensive customization possibilities
  • Solid Properties Materials Database

Data Processing Options

  • State-of-the-art post-processing tool, FlowSight™
  • Batch post-processing
  • Report generation
  • Automatic or custom results analysis
  • High-quality OpenGL-based graphics
  • Color or B/W vector, contour, 3D surface & particle plots
  • Moving and stationary probes
  • Visualization of non-inertial reference frame motion
  • Measurement baffles
  • Arbitrary sampling volumes
  • Force & moment output
  • Animation output
  • PostScript, JPEG & Bitmap output
  • Streamlines
  • Flow tracers

User Conveniences

  • Active simulation control (based on measurement of probes)
  • Mesh generators
  • Mesh quality checking
  • Tabular time-dependent input using external files
  • Automatic time-step control for accuracy & stability
  • Automatic convergence control
  • Mentor help to optimize efficiency
  • Units on all variables
  • Custom units
  • Component transformations
  • Moving particle sources
  • Change simulation parameters while solver runs
  • Launch and manage multiple simulations
  • Automatic simulation termination based on user-defined criteria
  • Run simulation on remote servers using remote solving
  • Copy boundary conditions to other mesh blocks

Multi-Processor Computing

  • Shared memory computers
  • Distributed memory clusters

FlowSight

  • Particle visualization
  • Velocity vector fields
  • Streamlines & pathlines
  • Iso-surfaces
  • 2D, 3D and arbitrary clips
  • Volume render
  • Probe data
  • History data
  • Vortex cores
  • Link multiple results
  • Multiple data views
  • Non-inertial reference frame
  • Spline clip