FLOW-3D 는 선박설계, 슬로싱 동역학, 파도에 미치는 영향 및 환기를 포함하여 해안 및 해양 관련 분야에 사용할 수 있는 이상적인 소프트웨어입니다.
자유 표면 유체 역학, 파동 생성, 움직이는 물체, 계선 및 용접 공정과 관련한 FLOW-3D 의 기능은 해양 및 해양 산업에서 CFD 공정을 모델링하는 데 매우 적합한 도구입니다. 해안 응용 분야의 경우 FLOW-3D 해안 응용 분야의 경우 FLOW-3D 는 해안 구조물에 대한 심한 폭풍 및 쓰나미 파동의 세부 사항을 정확하게 예측하고 돌발 홍수 및 중요 구조물 홍수 및 피해 분석에 사용됩니다. 기능은 다음과 같습니다.
자유 표면 – 파동 유체 역학 및 오버 토핑 : 규칙 및 불규칙파 및 파동 스펙트럼 (Pierson Moskowitz, JONSWAP)
Seakeeping – slamming, planing, porpoising 및 선체 선체 변위 : 완전히 결합된 선박 및 수중 차량 유체 역학
선체 – Resistance, stability and dynamics: surging, heaving, pitching and rolling motion (response amplitude operators or RAOs)
슬로싱 – LNG / 밸러스트 탱크
해양 공학 – 파동 에너지 변환기
해안 응용 분야의 경우, FLOW-3D 는 강력한 폭풍과 쓰나미 현상에 의한 해안 구조물이 받는 영향에 대한 세부 사항 예측, 돌발 홍수에 의한 중요한 시설물에 대한 정확한 피해 분석 등을 위해 사용됩니다.
Mooring Lines, Springs and Ropes
FLOW-3D (계류선 및 스프링 등)의 특수 물체를 다른 움직이는 물체에 부착하면 엔지니어가 선박 런칭, 부유 장애물 역학, 부표, 파도에너지 변환기 등을 정확하게 포착할 수 있습니다.
Welding
FLOW-3D 용접 모듈이 추가되면서 조선업계의 용접분야에서는 다공성 등 용접 결함을 최소화할 수 있어 선체의 품질을 크게 높이는 동시에 생산 시간을 최적화할 수 있습니다.
Coastal & Maritime Case Studies
FLOW-3D 사용자들은 연약한 해안선 보호, 구조물에 대한 파장 시뮬레이션, 선체 설계 최적화, 선박 내 환기 연구 등 해안 및 해양 애플리케이션에 FLOW-3D를 사용합니다.
우리는 보트가 세계 항해를 하면서 마주칠 것 같은 다양한 조건에서 항해를 할 수 있는지를 볼 수 있었습니다. 그리고 속도뿐만 아니라 연료 효율과 안전도 고려하도록 설계를 수정할 수 있었습니다. – Pete Bethune, skipper of Earthrace
Validation of Sloshing Simulations in Narrow Tanks / Aerial Landslide Generated Wave Simulations / Earthrace: Speed, Fuel Efficiency and Safety
Mary Kathryn WalkerFlorida Institute of Technology, mwalker2022@my.fit.edu Robert J. Weaver, Ph.D.Associate ProfessorOcean Engineering and Marine SciencesMajor Advisor Chungkuk Jin, Ph.D.Assistant ...
Waqed H. Hassan | Zahraa Mohammad Fadhe* | Rifqa F. Thiab | Karrar MahdiCivil Engineering Department, Faculty of Engineering, University of Warith Al-Anbiyaa, Kerbala 56001, IraqCivil ...
응고 모델은 열전달이 활성화되고(Physics → Heat Transfer → Fluid internal energy advection) 유체비열(Fluids → Fluid 1 → Thermal Properties → Specific heat)과 전도도(Fluids → Fluid 1 → Thermal Properties → Thermal Conductivity) 이 지정될 때 사용될 수 있다. 단지 유체 1만 상 변화를 겪을 수 있다.
응고모델을 활성화하기 위해 Fluids → Fluid 1 → Solidification Model 을 체크하고 물성 Fluids → Fluid 1 → Solidification Model 가지에서 Liquidus temperature, Solidus temperature, 그리고 Latent heat of fusion 를 지정한다. 가장 간단한 모델(Latent Heat Release Definition 에 펼쳐지는 메뉴에서 Linearly with constant 를 선택)에서, 잠열은 물체가 Liquidus 에서 Solidus 온도로 냉각될 때 선형적으로 방출된다. 고상에서의 상변화열을 포함하는, 잠열 방출의 더 자세한 모델을 위해 온도의 함수로 잠열방출을 정의하기 위해 Specific energy vs. temperature 또는 Solid fraction vs. temperature 선택을 사용한다. 이 지정에 대한 더 자세한 내용은 이론 매뉴얼의 Heat of Transformation 를 참조한다.
응고는 유체의 강직성 및 유동저항을 뜻한다. 이 강직성은 두 가지로 모델링 된다. 낮은 고상율에 대해 즉 Fluids → Fluid 1 → Solidification Model → Solidified Fluid 1 Properties → Coherent Solid Fraction 의 coherency 점 밑에서는 점도는 고상율의 함수이다. 간섭 고상율보다 큰 고상율에 대해서는 고상율의 함수에 비례하는 항력계수를 갖는 Darcy 형태의 항력이 이용된다. 이 항력은 모멘텀 방정식에 (bx,by,bz) 로써 추가된다- Momentum Equations 를 보라. 이 항력의 계산은 Solidification Drag Model 에서 기술된다. 항력계수는 사용자가 유동저항에 양을 조절할 수 있는 Coefficient of Solidification Drag 인자를 포함한다. 항력계수는 FLOW-3D 출력에서 기록된 속도에 상응하는 지역 상 평균 속도에 의해 곱해진다.
Fluid 1 Properties)을 지나면 항력은 무한대가 되고 계산격자 관련하여 유동이 있을 수 없다(단 예외로 Moving Solid Phase를 참조).
Note
모든 유체가 완전히 응고하면 모사를 정지시키기 위해 General → Finish condition → Solidified fluid fraction 를 이용한다. General → Finish condition → Finish fraction 은 모사를 중지하기 위한 고상율 값을 정한다.
주조 시 mushy zone 은 액상과 고상이 혼합물로 존재하는 지역이다. 이 지역 혼합점도는 동축의 수지상 조직(과냉각된 액체 안에서 방사상으로 자라는 결정으로 된 구조) 이 액체 안에서 자유롭게 부유할 때 영향을 미친다.
일단 수지상 조직의 간섭성이 발생하여 고정된 고상 망이 형성되면 액상이 고정된 다공 수지상 구조를 통과해야 하므로 추가의 유동손실이 발생한다. 다른 방법으로는 간섭점을 지난 액/고상 혼합물은 다공물질을 통한 유동 대신에 고점도의 유체로 간주될 수 있다. 점성유체로 간주하는 접근은 예를 들면 연속 이중 롤 주조 과정같이 고상이 계속 이동 및 변형할 때 유용하다.
Solidification Drag Models in FLOW-3D, FLOW-3D 내 응고 항력모델
응고에 의한 항력계수를 정의하기 위해 사용자는 우선 열전달 및 응고모델을 활성화 해야 한다. 이들은 Model Setup → Physics 탭 에서 활성화될 수 있다. 수축모델 또한 응고모델 창에서 활성화될 수 있다.
일단 Solidification 모델이 활성화되면 항력의 공식이 지정될 필요가 있다. Solidification대화의 밑 좌측 모퉁이에서 Porous media drag-based 와 Viscosity-based 의 항력공식 중의 선택을 한다.
Viscosity-based 공식은 점성 유체로 취급하며 Viscosity 영역 내Flow model for solidified metal 입력 밑에서 지정되는 순수 고상 점성을 갖는 고상화된 유체로 간주된다. 이 접근법은 경직성의 항력모델(즉, 응고 금속이 롤러 사이로 압착될 때)을 사용할 수 없는 경우의 모사에 이용된다. 이 점성은 고상율에 따라 선형으로 변한다.고상율이0일 때 점도는 유체1의 점도이다.고상율이1이면 점도는 Solidification 패널에서 지정된 값과 같다.
Porous media drag-based 공식은 응고상태를 결정하기 위해 고상율을 사용한다. 고상율이 Critical Solid Fraction 이거나 초과하면 이때 항력은 무한대가 된다-즉, 액상/고상 혼합물은 고체같이 거동한다. 고상율이 Coherent Solid Fraction 보다 작으면 항력은 0이다. 이 두 값 사이에서 유동은 mushy 지역에 있고 이를 통한 유동은 마치 다공질 내에서의 유동같이 처리된다. 또한 모델은 고상율이 Coherent Solid Fraction 보다 작을 때 자동적으로 용융 금속의 점도를 조절한다. 이 상태에서 고상결정은 점도를 올리지만 결합하지는 않는다(즉, 간섭 없음). 일단 유체가 Coherent Solid Fraction 에 도달하면 항력방정식이 고려되고 점도는 간섭성에 도달하기 전의 값으로 일정하게 된다. 임계 및 간섭 고상율은 사용자가 정의하며 논문이나 책 등에서 찾을 수 있다. 이 식에서는 Coefficient of Solidification Drag 가 정의되어야 한다. 이는 Solidification 창 또는 Fluid 1 → Solidification Model→Solidified Fluid 1 Properties tree → Other 트리를열어 Model Setup →Fluids 탭에서 될 수 있다.
How to Calculate Permeability 투과성 계산법
밑에 주어진 Darcy법칙은 수지상 구조를 위한 다공매질내의 수학적 유동기술이다.[Poi87].
(19)
여기서 u 는 수지상 구조 내 유동의 속도이고 ∇P 는 지역 압력구배, 그리고 K 는 mushy 구역의 특정 투수성이다. 이 방정식은 단지 유동이 거의 정상 상태이고, 관성효과가 없으며 유체의 체적율이 일정하고 균일하며 액체-액체의 상호작용 힘이 없을 때 유효하다. 투수성을 정의하는데 이용될 수 있는 대 여섯 개의 모델이 있으나 FLOW-3D 는 밑에 보여주는 Blake-Kozeny 을 이용한다. 다른 모델들은 코드와 함께 제공되는 소스코드를 사용자 사양에 맞게 수정하여 추가할 수 있다.
(20)
여기서
C2 는 전형적으로 와 같은 비틀림
fs 는 고상율이고
λ1는 유동을 위한 특정 치수
이 응용에서 수지상 가지 간격(DAS)이 이용된다.
식 (11.19) 을 식(11.20) 에 적용하면 투수성을 위한 다음 식을 얻는다.
(21)
수지상 가지 간격(DAS)에 대한 일반적인 값들은 밑에 주어져 있다.
Range of Cooling Rates in Solidification Processes¶
COOLING RATE, K/s
PRODUCTION PROCESSES
DENDRITE ARM SPACING,
to
large castings
5000 to 200
to
small castings, continuous castings, die castings, strip castings, coarse powder atomization
200 to 5
to
fine powder atomization, melt spinning, spray deposition, electron beam or laser surface melting
5 to 0.05
Range of cooling rates in solidification processes [CF85]
How FLOW-3D Defines the Coefficient of Solidification Drag FLOW-3D 가 응고 항력계수를 결정하는법
FLOW-3D 는 액고상 변화를 모델링하기 위해 다공매질항력을 이용한다. 항력은 고상율의 함수이다. 사용자에게 두 수축모델이 이용 가능하다; 급속 수축 모델 과 완전 유동모델. 급속 수축 모델은 상변화와 연관된 체적변화를 고려하지 않으며 유체는 정지해 있다고 가정한다. 완전 유동모델은 상변화가 관련된 체적변화를 고려한다. 항력은 투수성에 역으로 비례하므로 다음과 같이 표현될 수 있다.
(22)
여기서, Fd 는 FLOW-3D 에서 사용된 항력계수이다. 이 항력계수는 지역 속도에 의해 곱해지고 모멘텀 방정식의 오른쪽에서 차감된다 (Momentum Equations 참조). 식 (11.22) 를 재정리하고 식 (11.21) 로부터의 투수성에 치환하면 다음을 얻는다.
Macro-Segregation during Alloy Solidification 합금응고시 거시적 편절
편절 모델은 대류와 확산에 의한 용질 이동에 따른 이원합금 요소에서의 변화를 모델링 하도록 되어 있다. 이 모델링은 Physics → Solidification 로 부터 될 수 있다.
Activate binary alloy segregation model 을 체크하고 편절 모델을 활성화한다.
여러 온도에서 평형에 있는2원합금 요소농도를 정의하는 상태도는 직선의 고상선 및 액상선을 가진다고 가정된다. 상태도는 입력데이터에 의해 구성되고 전처리 그림파일 prpplt 에 포함된다. Analyze → Existing 에서 이용 가능하다
Macro-Segregation Model (under Fluids → Fluid 1 → Solidification Model)에 관련된 일부 유체물성 트리가 밑에 보여진다. 상태도는 Reference Solute Concentration 에서의 the Solidus 와 Liquidus Temperatures 값들에 의해 정의된다. 추가로 Concentration Variables 밑의 Partition coefficient 도 정의되어야 한다. 그렇지 않으면 Pure Solvent Melting Temperature 가 정의될 수 있다. Partition coefficient 와 Pure Solvent Melting Temperature 둘 다가 지정되면 용매 용융 온도는 상태도로부터 재 정의된다.
Eutectic Temperature 또는 Eutectic Concentration 는 융해작용을 정의하기 위해 지정될 수 있다. 또 이 두 변수가 다 지정되면 Eutectic Concentration 은 상태도에서 재 정의된다.
Diffusion Coefficients 는 고상과 액상 사이의 용질의 확산계수 비율을 정의한다. 액체 내의 용질의 분자 확산계수는 Physics → Solidification 에서 specifying Solute diffusion coefficient 를 지정함으로써 정해진다. RMSEG 는 용질의 난류 확산계수 승수를 정의한다; 이는 입력파일에서 직접 지정된다.
용질 재 분배에 의한 농도변화가 중요하면 Physics → Density evaluation → Density evaluated as a function of other quantities를 정하고 용질농도의 선형함수로써 금속농도를 정의하기 위해 Fluids → Segregation model 밑의 Solutal Expansion Coefficient 를 용질 확장계수로 지정한다. 이 경우 Reference Solute Concentration 이 기준농도로 사용될 것이다. 추가로 Fluids → Fluid 1 → Density Properties → Volumetric Thermal Expansion 은 액체 내 열부력 효과를 참작하기 위해 지정될 수 있다(또한 Buoyant Flow참조).
초기 용질농도는 Meshing & Geometry → Initial → Global → Uniform alloy solute concentration 에서 지정될 수 있다. 불 균일한 초기 분포는 Alloy solute concentration 밑의 초기유체 구역 안에서 정의될 수 있다. 추가로 농도는 Initial Conditions: Region Values 에서 기술된 바와 같이 2차함수를 사용하는 부분을 편집하여 공간상의2차함수로 변화할 수 있다. 압력과속도 경계에서 용질 경계조건을 정하기 위해 Boundaries → Boundary face → Solute concentration 를 이용한다.
액상 및 고상 구성은 후처리에서 데이터 변환을 이용하여 그려질 수 있다. 용융 응고금속은 금속 내 용융의 질량 분율을 저장하는 SLDEUT 를 그림으로써 가시화될 수 있다.
액상 내 열구배가 크면 Physics → Heat Transfer → Second order monotonicity preserving 를 지정함으로써 더 나은 정확성을 위해 고차원 이류법을 사용한다.
mushy 지역에서의 유동손실은 수지상 가지 간격(DAS)의 함수인 Fluids → Fluid 1 → Solidification Model → Solidified Fluid 1 Properties → Coefficient of Solidification Drag 에 의해 조절된다. 후자는 이 모델에 의해 계산되지 않으므로 사용자는 Coefficient of Solidification Drag 를 지정해야 한다
Note
표준 응고모델 과는 달리 상태도상의 용융점을 지나 고상선을 외삽하여 정의되므로 여기서 응고선의 값은 음수일 수 있다.
Microporosity Formation 미세다공형성
미세다공모델은 단지 응고(Solidification참조)를 모델링할 때 사용될 수 있고 Physics → Solidification → Activate micro-porosity model 에서 활성화된다. 필요한 입력은 Fluids → Densities → Fluid 1 and Fluids → Solidification Properties → Solidified Fluid 1 Properties → Density 에서 정의되는 액체와 고상 유체밀도이며 고상유체밀도는 액체밀도보다 크다. 또한 Fluids → Fluid 1 → Solidification Model → Solidified Fluid 1 Properties 안에 있는 Critical Solid Fraction 은 1.0보다작게 설정되어야 한다.
Square of the speed of sound at critical solid fraction 값이 정의될 수 있다. 이는 수축에 의해 mushy 지역에서 전개되는 커다란 음압에서의 응고유체의 압축성을 기술한다. Critical pressure at which gas pores can form 값은 모델이 Initial tab 탭에서 또는 재 시작 데이터에서 정의되는 유체내의 초기 압력과 결합되도록 한다.
Intensification pressure 또한 다공 생성을 지연시키기 위해 응고 시 shot sleeve plunger 에 의해 형성되는 추가압력을 고려하기 위한 고압 주조모사를 위해 정의될 수 있다. Intensification pressure 가 클수록 더 적은 양의 다공이 주조 시 응고 과정에서 발생할 것이다.
미세 다공 모델은 응고 모델의 활성화 이외의 어떤 다른 설정을 필요로 하지 않는다. 이는 완전 유동방정식이나 속도장이 0인 경우, 즉 순수한 열 문제에서도 함께 사용될 수 있다.
이 모델은 후처리 과정의 공간 및 이력에서 사용 가능한 Percent micro-porosity 라고 불리는 추가 출력 양을 생성한다.
이동고상 선택은 연속주조 모델링을 가능하게 한다. Continuous Casting Phantom 요소는 응고된 이동 유체가 있는 지역에서 정의된다. 이는 지정된 영역을 차지하지만 정의에만 존재하므로 환영요소라고 한다. 이는 실제로 면적이나 체적을 차지하지 않으므로 체적이 없고 결과에서도 고체요소로 보이지 않는다. 이는 Meshing & Geometry → Geometry → Component → Component Type 옆 펼쳐지는 메뉴에서 정의된다.
다른 방법으로는 입력파일(prepin.*)에서 IFOB(N) 변수가 4로 지정되고 N 은 요소 번호이다. 이 파일은 File → Edit Simulation…. 을 선택하여 이용될 수 있다. 또한 입력파일에서 시간의 함수(TOBS(t) 에 의해 지정되는)일 수 있는 가상 요소의 속도성분 UTOBS(t,N), VTOBS(t,N) 그리고 WTOBS(t,N) 이 지정된다.
Fluids → Fluid 1 → Solidification Properties → Solidified Fluid 1 Properties → Coherent Solid Fraction 에 의해 정의된 간섭 고상율 보다 큰 고상율에 대해서는 Darcy 형태의 항력 이 유체를 가상 요소의 속도로 움직이게 하는데 사용된다. 고상율이 Fluids → Fluid 1 → Solidification Properties → Solidified Fluid 1 Properties → Critical Solid Fraction 에서 지정된 경직점을 능가하게 되면 가상 요소의 속도를 따라 움직일 것이다.
Note
가상 요소는 요소 그림에 안 나타나나 Component number 를 그릴 때는 보여진다.가상 요소는 균일속도가 요소의 전체에 적용되므로 평평해야 한다.
Solidification Shrinkage 응고수축
체적 수축은 소재가 응고하고 응고소재의 밀도가 액체소재의 밀도보다 클 때 나타난다(즉, Fluids → Fluid 1 → Solidification Model → Solidified Fluid 1 Properties → Density > Fluids → Fluid 1 → Density Properties → Density). 수축모델은 그러므로 Solidification 모델이 활성화되어야 하고 고상/액상의 두 밀도가 정의되어야 한다. 수축은 단지 1유체의 뚜렷한 경계면 문제에서만 모델링 될 수 있다.
두 가지 수축모델이 있다. Shrinkage model with flow effects 를 선택하면 완전 열 유체방정식을 해석한다(이론 매뉴얼의Solidification Shrinkage and Porosity Models 참조). 그러나 이 모델은 특히 장시간의 응고가 고려되면 컴퓨터 계산시간이 많이 소요된다. 다른 방법으로 사용자 Interface 에 Shrinkage model 이라고 불리는 단순모델이 있다.
이 모델은 단지 열전달 방정식의 해석에 의존하며 특히 내재적 열전달 모델 (Numerics → Explicit/implicit options → Heat transfer → Implicit –Thermal solution 참조)과 사용시에 빨리 해석할 수 있다. 액체 체적 내로의 유동 통로가 없을 때 내부공동이 발생한다.
이 두 모델에서 유입은 mushy 지역 유동에 대한 항력계수를 계산함으로써 정의된다. 격자 내 모든 점에서의 항력함수는 상수승수 Fluids → Solidification properties → Other → Coefficient of Solidification Drag (Solidification Drag Model 참조)를 가지는 지역 고상율의 함수로 계산된다. 항력함수의 역의 값은 공간 그림에서 가시화 될 수 있다: 이 그림을 위한 변수이름은 ‘drag coefficient’ 이다.
Mushy 지역에서의 커다란 유동 손실에 따른 부분적 유입이 Shrinkage model with flow effects 에서 발생할 수 있지만 단순화된 Shrinkage model 은 완전 유입이 아니면 유입이 없게 된다. 후자는 유입 통로를 따라 지역 고상율이 Fluids → Fluid 1 → Solidification Model → Solidified Fluid 1 Properties → Critical Solid Fraction (디폴트는1.0)에서 정의된 임계값보다 커질 때 발생한다. 추가로 고립된 액체 내의 금속의 고상율이 Coherent Solid Fraction 에 도달할 때까지 단순모델에서의 유입은 고립부 상부로부터 발생한다. 그 후로는 유입이 고립부의 가장 뜨거운 부분에서부터 발생한다.
모든 유체가 완전히 응고되면 모사가 정지하도록 General → Additional finish condition → Solidified fluid fraction 를 사용한다. 변수 Finish fraction 는 유체가 지정된 고상율에 도달할 때 모사가 정지하도록 하는데 사용될 수 있다.
Note
이송 방향을 결정하기 위해 단순 수축 모델에서 중력이 필요하며 좌표축 중 하나를 따라야합니다. 둘 이상의 중력 구성 요소가 0이 아닌 경우, 가장 큰 중력 구성 요소가 공급 방향을 결정하는 데 사용됩니다.
FLOW-3D 시뮬레이션에서 일반적인 움직이는 물체 (GMO)는 사용자가 규정하거나 유체 흐름과 동적으로 결합되는 모든 종류의 모션을 가진 강체입니다. 고정된 축 / 포인트와 같은 6 자유도 또는 모션 구속 조건을 가질 수 있습니다. 규정된 힘과 토크는 결합된 동작 하에서 GMO에 적용될 수 있습니다. GMO 모델은 충돌 및 연속 접촉을 포함하여 강체 상호 작용뿐만 아니라 독립적인 동작 유형에서 여러 개의 강체를 허용합니다. 이 모델은 견고하고 효율적이며 강력하고 상업용 전산 유체 역학 소프트웨어에서 FLOW-3D가 유일합니다.
FLOW-3D 자동차 차동 부분의 3D 시뮬레이션
모델링 기능
6 개의 DOF를 소유 할 수있는 최대 500 개의 움직이는 물체를 허용하거나 고정된 축 또는 고정된 점을 중심으로 회전 할 수 있습니다. 다른 모션 제약 패턴도 가능합니다.
물체는 유체 흐름과 완벽하게 결합되거나 사용자가 모션을 처리할 수 있습니다.
물체는 밀도로 특징 지어지는 여러 가지 재료로 만들 수 있습니다.
객체의 기하학 및 동작의 복잡성에 대한 제한이 없습니다.
지정된 시간에 따른 힘과 토크를 대상에 적용 할 수 있습니다.
모델 충돌 및 움직이는 물체와 움직이지 않는 물체 사이의 지속적인 접촉
스프링과 로프는 물체에 닿을 수 있습니다.
개체에서 다공성 허용
열전도 및 대류가 허용됩니다.
유압식, 중력 식, 비 관성식, 스프링 식, 사용자 정의 제어력 및 토크는 결합 된 모션이있는 물체에 대해 고려됩니다.
시각 및 수치 출력을 포함한 완벽한 후 처리 기능
움직이는 물체 시뮬레이션
FLOW-3D고객은 처음 사용하는 방법보다 더욱 효율적으로 움직이는 물체 모델의 적용을 사용했습니다. waterwheels에서 shot sleeve, 에너지 디바이스 파동에 이르기까지 우리는 복잡한 메쉬 및 집중적인 컴퓨팅 리소스에 의존하지 않고도 엔지니어링 문제를 해결하는데 모델을 사용하는 방법에 깊은 인상을 받았습니다. 이 모델을 사용하는 진지하고 상상력이 좋은 예를 모두 검토하려면 YouTube 재생 목록 을 방문하십시오.
본 자료는 국내 사용자들의 편의를 위해 원문 번역을 해서 제공하기 때문에 일부 오역이 있을 수 있어서 원문과 함께 수록합니다. 자료를 이용하실 때 참고하시기 바랍니다.
FAVOR™ vs. Body-Fitted Coordinates
The simplicity of the fractional area/volume method FAVOR™ for modeling complex geometric regions is very attractive. But, can it compete in terms of accuracy with deformed grids such as those employed by finite-element or body-fitted coordinate methods? A comparison between these methods shows that there are only small differences between the capabilities of the two approaches.
복잡한 기하 영역을 모델링하기 위해 면적/체적 점유율을 측정하는 FAVOR TM 법의 간결함은 매우 매력적입니다. 그러나 유한 요소법과 BFC (body-fitted coordinate : 물체 적합 좌표) 법에서 사용되는 변형 격자와 비교할 때 정확도면에서 경쟁 할 수 있을까요? 이러한 기술을 비교함으로써 이 두 가지 방법의 기능의 차이는 매우 적다는 것을 알 수 있습니다.
In the absence of solid boundaries there is little fundamental difference between grids of differently shaped control volumes. Some methods require users to store more information (e.g., node locations and various geometric factors) and some exhibit differing levels of accuracy depending on the amount of element distortion. In all cases, however, the underlying idea is a discrete approximation in which fluid forces and fluxes are computed for each element in the grid.
고체 경계가 존재하지 않는 경우, 형상이 다른 컨트롤 볼륨 격자 사이에는 근본적인 차이는 거의 없습니다. 어떤 방법이 더 많은 정보를 저장하도록 요구할 경우 (노드 위치, 다양한 기하학적 형상 요소 등)나 요소의 왜곡의 양에 따라 정확도 수준은 다를 수 있습니다. 그러나 어떤 경우에도 기초가 되는 개념은 격자의 각 요소마다 유체의 힘과 플럭스가 계산된 근사치입니다.
It is the issue of obstacle boundaries that is most often raised as an advantage of deformable grids because they can be constructed to fit the geometry. Two consequences come with this flexibility. One consequence is that these grids must be unstructured for general use. This is because structured grids can only undergo limited distortion before elements are turned inside out. The other consequence is that the distortion of an element makes it more difficult to construct accurate numerical approximations.
이것은 기하학적 형상에 적합하도록 구축 될 수있는 변형 가능한 격자의 장점으로 거론되는 가장 흔한 장애물 경계의 문제입니다. 이러한 유연성은 2 개의 결과가 포함됩니다. 하나는 이러한 격자는 범용으로 비구조이어야한다는 것입니다. 구조 격자의 경우 조금 왜곡을 더한 것만으로 요소가 뒤집어지기 때문입니다. 또 하나는 요소의 왜곡에 의해 정확한 수치 근사를 구축하는 것이 더욱 어려워 질 것입니다.
Structured FAVOR™ Grids vs. Unstructured Grids
The FAVOR™ concept can be used in connection with any type of grid including grids consisting of rectangular or distorted elements and whether the grid is structured or unstructured. Structured grids are best because they are easy to generate and the indices for neighboring elements are known. Rectangular grid elements make it easy to compute the fractional areas and volumes of elements used by the FAVOR™ method.
FAVOR TM의 개념은 사각형 요소와 왜곡 요소로 이루어진 격자를 포함하여, 또한 구조 격자 또는 비구조 격자 여부에 관계없이 모든 종류의 격자와 관련시켜 사용할 수 있습니다. 그 중에서도 구조 격자가 좋습니다. 구조 격자는 쉽게 생성 할 수 인접한 요소의 인덱스를 알고 있기 때문입니다. 직사각형의 격자 요소를 사용하면 FAVOR TM 법에서 사용되는 요소의 면적 점유율과 체적 점유율을 쉽게 계산할 수 있습니다.
The ease with which structured, rectangular grids can be generated makes this an obvious choice for the FAVOR™ method. However, the specification of open and closed grid regions implicit in the FAVOR™ method introduces a type of unstructured computing environment, because only those elements having a finite amount of open volume are actually computed. This results in a computation that is analogous to an unstructured grid computation.
직사각형 구조 격자를 쉽게 생성 할 수 있기 때문에 FAVOR TM 법에서는이를 선택하는 것이 당연하게되어 있습니다. 그러나 FAVOR TM 법에서는 열린 격자 공간과 닫힌 격자 영역의 지정을 잠재적하는 일종의 비구조 계산 환경입니다. 이것은 열린 체적의 크기가 유한한 요소만 실제로 계산되어 그 결과, 비 구조 격자의 계산과 비슷한 계산이 되기 때문입니다.
The analogy with a true unstructured grid is not perfect. For example, the FAVOR™ method requires storage for all elements whether they are blocked or not, while unstructured grids require storage of neighbor lists. Storing all elements is the price paid for automatically knowing which elements are neighbors. Of course, if heat conduction is to be computed in the solid regions surrounding a fluid, then blocked elements must be kept anyway.
진짜 비구조 격자와 유사하지만 완전하지 않습니다. 예를 들어, FAVOR TM 법은 차단되어 있는지 여부에 관계없이 모든 요소의 스토리지가 필요하지만, 비 구조 격자는 인접리스트의 스토리지가 필요합니다. 어떤 요소가 인접해 있는지를 자동으로 알 수 있는 대신 모든 요소를 저장해야 합니다. 물론, 유체를 둘러싼 고체 영역에서 열전도가 계산되는 경우 차단되는 요소는 어느쪽으로도 유지해야합니다.
One disadvantage associated with rectangular, structured grids is that they cannot be distorted to increase resolution in a localized region. Grid lines can be constructed closer together in a particular region for increased resolution, but these grid lines then extend across the entire grid. As a result, the number of elements in a grid may become large. Offsetting this disadvantage is the simplicity of grid generation. There is also the possibility of introducing multiple grid blocks joined at their boundaries to provide increased local resolution without a large increase in the number of grid cells.
직사각형 구조 격자와 관련된 단점 중 하나는 국소화된 영역에서 변형시켜 해상도를 올릴 수 없는 것입니다. 격자선은 특정 영역에서 해상도를 높이기 위해 세밀하게 만들 수 있지만, 이러한 격자는 그 격자 전체에 퍼집니다. 그 결과 하나의 격자 내의 요소의 수가 증가 할 수 있습니다. 이 단점을 상쇄하는 것이 격자 생성의 간단함입니다. 또한 경계에 결합된 복수의 격자 블록을 도입하여 격자 셀의 수를 크게 증가시키지 않고 국소적인 해상도를 올릴 가능성도 있습니다.
The FAVOR™ method is seen to have properties analogous to unstructured grids but without the overhead associated with the construction of unstructured body-fitted or finite-element grids. Since grid generation often requires a major investment of time and effort, the elimination of this task makes FAVOR™ a very desirable alternative.
FAVOR TM 법은 비 구조 격자와 유사한 특성을 가진 것으로 볼 수 있지만, 비 구조의 물체 적합 격자와 유한 요소 격자의 구축과 관련된 오버 헤드는 없습니다. 격자 생성에 시간과 노력을 많이 투자해야하는 경우가 많기 때문에이 작업이 불필요하게 됨으로써 FAVOR TM은 매우 바람직한 선택이 되고 있습니다.
FAVOR™ Elements vs. Body-Fitted Coordinate Elements
A second issue with deformed grids is that their deformation can be used to fit a bounding solid surface. This is accomplished by moving the nodes of elements closest to the surface onto the surface. In contrast, in the FAVOR™ method a surface is allowed to cut through an element and its location is recorded not by moving the edges of the element but in terms of the fractional face areas and fractional volume of the element that are not covered by the solid.
변형 격자의 두 번째 문제는 경계의 고체 표면에 대한 적합성에 변형을 사용할 수 있는 것입니다. 이것은 표면에 가장 가까운 요소 노드를 표면으로 이동함으로써 달성됩니다. 반대로 FAVOR TM 법에서는 표면에서 요소를 자를 수 있도록 허용되고 그 위치는 요소의 끝을 이동함으로써가 아니라 고체로 덮여 있지 않은 요소의 표면적 점유율 및 볼륨 점유율로 기록됩니다.
What we wish to show is that this fractional area/volume technique for defining solid boundaries has the same consequences for numerical approximation as does a deformed (i.e., body-fitted coordinates) grid technique.
우리가 보여주고 싶은 것은 경계를 정의하는 면적/체적 점유율 기법의 수치 근사의 결과가 변형 (BFC) 격자 기법과 동일한 결과를 가지고 있다는 것 입니다.
The most important point to recognize about the FAVOR™ method is that approximations of fluid-dynamic quantities are restricted to the open regions of elements. This restriction introduces fractional areas and volumes of elements as factors directly into the discrete approximations. For example, the flux of a quantity from one element to another has the fractional area of the fluxing boundary that is open to flow as a multiplier.
FAVOR ™ 법에 대해 이해하는 가장 중요한 점은 유체 역학의 유량의 근사치는 요소가 열려있는 영역에 한정되는 것입니다. 이 제한은 요소의 면적/체적 점유율은 인자로 분리된 사치에 직접 도입됩니다. 예를 들어, 요소 요소에 유량의 플럭스는 승수(multiplier)로 흐름에 대해 열려있는 플럭스 경계의 면적 점유율을 가지고 있습니다.
In this way FAVOR™ and body-fitted coordinates both compute fluxes across the faces of elements that employ the same areas. In FAVOR™ the areas are stored as fractions of the original element face areas. In a BFC method the areas are computed from the coordinates of the nodes defining the faces, and often times they are also stored so they don’t have to be recomputed.
이처럼 FAVOR TM과 BFC는 모두 동일한 면적을 사용하는 요소의 면을 통과하는 유량을 계산합니다. FAVOR TM는 면적은 원래 요소 표면적의 비율로 저장됩니다. BFC 법에서는 면적은 표면을 정의하는 노드의 좌표에서 계산됩니다. 많은 경우 저장도 되기 때문에 재 계산의 필요가 없습니다.
Figure 1: FAVOR™ blocked cell left (a) and BFC cell right (b). Solid region is shaded.
When constructing difference approximations in a grid of non-uniform elements it is necessary to know the effective element widths in different directions. A simple example is given in Fig. 1, which shows an element with a solid boundary. In this case the solid boundary is parallel with the vertical sides of the element. For a deformed element, Fig. 1b, the width of the cell would be h and this width would be used for computing differences in the horizontal direction.
불균일 요소 격자에서 차분 근사를 구축 할 때 다양한 방향의 유효 성분 폭을 알고 있어야 합니다. 간단한 예를 그림 1에 나타냅니다. 이것은 고체 경계를 가지는 요소를 보여줍니다. 이 경우 고체 경계 요소의 수직 측면과 평행하고 있습니다. 변형 된 요소 (그림 1b)에 대해 셀의 폭은 h되어,이 폭을 사용하여 가로 방향의 차이가 계산됩니다.
In FAVOR™ the width of the open portion of the cell is equal to the product of the open volume fraction and the original cell width. It is this product that is used in FAVOR™ for difference approximations in the horizontal direction and it is the same as the width of the deformed cell, h.
FAVOR TM에서 셀의 열 부분의 폭은 열린 체적 점유율과 원래의 셀 폭의 곱 같습니다. 이것은 수평방향에서의 근사 차이에 대한 FAVOR TM이 사용된 제품이며 이는 변형 된 셀의 폭 h와 같습니다.
When differences are computed in a vertical direction, for example, across the top of the element (i.e., parallel to the solid boundary in Fig. 1) the FAVOR™ method approximation involves a ratio of the fractional area at the top divided by the fractional volume. This ratio has a value of unity because the solid boundary blocks the same fraction of area horizontally as it does volume, which again makes the approximation similar to that of a BFC cell. Thus, there is no difference between the effective cell widths used in making difference approximations in FAVOR™ and body-fitted coordinate cells.
세로 방향, 예를 들어 요소의 표면과 교차 (그림 1의 고체 경계와 평행) 방향으로 차이를 계산하는 경우 FAVOR TM 법의 근사치는 윗면의 면적 점유율을 체적 점유율로 나눈 비율이 관계 합니다. 이 비율의 값은 1입니다. 이것은 고체 경계에 따라 부피와 마찬가지로 동일한 면적 점유율이 옆으로 차단되기 때문입니다. 따라서 여기서도 근사 BFC 셀과 비슷합니다. 이처럼 FAVOR TM 셀과 BFC 셀에서 차등 근사 할 때 사용하는 유효한 셀 폭에 차이는 없습니다.
Figure 2: FAVOR™ cell left (a) and body-fitted coordinate cell right (b). Solid region is shaded.
If the solid surface is slanted as shown in Fig. 2, then the top and bottom face areas are different but are still known quantities in either case. Further, the average horizontal width of the open portion of the cell (h) is still the same because the volume fraction does not change when the solid surface is rotated about a midpoint within the cell (see dashed line). Here again there is little to distinguish between a FAVOR™ and a body-fitted coordinate cell.
그림 2와 같이 고체 표면이 경사져있는 경우에는 위 아래의 표면적은 다르지만 두 경우 모두 알려진 양입니다. 또한 셀의 열 부분의 가로 방향의 평균 폭 (h)도 마찬가지입니다. 이것은 셀의 중간 지점 (점선 참조)을 중심으로 단단한 표면을 회전해도, 체적 점유율은 변하지 않기 때문입니다. 여기서도 FAVOR TM 셀과 BFC 셀을 구별 할 수있는 것은 거의 없습니다.
Summary
It is well known that body-fitted coordinate grids can be difficult to construct, and there has been a huge effort devoted to the development of “automatic” grid generators. It is also well known that even with the best of grid generators it still takes a significant amount of time to establish a workable and well-behaved grid.
The simple rectangular construction of FAVOR™ grids makes them extremely easy to generate. Fractional areas and fractional volumes must be computed to define obstacles placed within a grid, but theses computations are well defined and easy to automate using simple algorithms.
Numerical accuracy is not sacrificed when selecting FAVOR™ over a body-fitted coordinate gridding method. The two approaches simply represent different ways to approximate bounding surfaces.
Finally, the numerical advantages inherent in the structured, smoothly varying, strictly orthogonal grids used by the FAVOR™ method should not be overlooked. These advantages, as well as the ability to automatically represent porous media (i.e., another example of a fractional area/volume region) are additional reasons why the FAVOR™ method was selected as the basis of FLOW-3D.
BFC 격자의 구축이 어려운 것은 잘 알려져 있으며 “자동”격자 생성 기능의 개발에 많은 노력을 쏟고 있습니다. 또한 최고 수준의 격자 생성 기능을 사용하더라도 유효하고 적절한 기능을 수행 격자를 조합하려면 상당한 시간이 걸릴 수 있다는 것이 잘 알려져 있습니다.
FAVOR TM 격자는 간단한 직사각형 구조이기 때문에 매우 쉽게 생성 할 수 있습니다. 격자에 배치 장애물을 정의하는 면적 점유율과 체적 점유율을 계산해야 하지만, 이러한 계산은 잘 정의되어 있으며, 간단한 알고리즘을 사용하여 쉽게 자동화 할 수 있습니다.
격자 생성 기법으로 BFC 대신 FAVOR TM을 선택하여 수치적 정확성이 희생되는 것은 아닙니다. 이 두 가지 방법은 경계 표면의 근사 방법이 다를 뿐입니다.
마지막으로, FAVOR TM 법에서 사용되는 구조화 된 변동이 부드럽고 엄격한 직교 격자 고유의 수많은 장점은 간과해서는 안됩니다. 이러한 장점과 다공성 매체 (면적 / 체적 점유율 영역의 또 다른 예)를 자동으로 표현하는 기능은 FAVOR TM 법이 FLOW-3D의 기초로 뽑힌 또 하나의 이유입니다.
FLOW-3D의 수치 모델링 기능은 코팅 성능 향상에 관심이 있는 엔지니어에게 이상적입니다. 계산 시뮬레이션은 코팅 흐름에 영향을 미치는 다양한 물리적 공정의 상대적 중요성과 효과를 연구하는 훌륭한 방법입니다. 물리적 테스트에서 프로세스를 분리하거나 해당 프로세스의 규모를 임의로 조정하는 것이 항상 가능한 것은 아닙니다. 이 섹션에서는 리 블릿 형성(rivulet formation), 핑거링(fingering), 증발, 거친 표면 위의 접촉선 이동 및 유체 흡수와 관련하여 FLOW-3D의 정적 및 동적 접촉각 처리에 대해 설명합니다.
Static and Dynamic Contact Angles
FLOW-3D는 입력으로 설정된 정적 접촉각의 함수로 동적 접촉각과 자유 표면 인터페이스에서 작용하는 관련 힘을 정확하게 계산하여 유체의 소수성을 캡처 할 수 있습니다. 아래 시뮬레이션은 물방울이 경사면 아래로 이동함에 따라 정적 접촉각이 동적 접촉각에 미치는 영향을 보여줍니다.
L.M. Hocking 박사는 그의 저서 [“A moving fluid interface on a rough surface,” J. Fluid Mech., 76, 801, (1976)]에서 표면에 미세한 요철이 흐름 구조를 유도하기 때문에 Contact line이 고체 표면을 통해 이동할 수 있으며 이는 거시적 관점에서 “velocity slip”로 해석 될 수 있다고 했습니다.
이 가설에 대한 전산 해석은 FLOW-3D를 이용하여 쉽게 수행됩니다. 선택된 테스트는 가로, 규칙적으로 이격 된 직사각형 슬롯 패턴 이차원 고체 표면 구성됩니다. 슬롯은 2mm 깊이 10mm 폭, 그리고 그들 사이 폭 10mm 고체 조각을 갖고 이격 됩니다. 이 크기는 전형적으로 상대적으로 부드러운 표면에 긁힌 모양입니다. 액체와 고체 사이의 정적인 접촉각이 60 °가 되도록 선택 하였습니다. 작동 유체는 물로 선정되었고 시험은 채널을 통해 속도30cm / s의 평균 물높이 15mm의 채널의 바닥에 있는 거친 표면을 두고 구동 이루어져 있습니다. 채널의 상단은 free-slip boundary로 정해집니다.
Hocking의 주장대로 micro-scale 교란이 Large scale 관점에서 보았을 때 계산된 속도장으로 보면 velocity slip의 한 종류로서 해석 될 수 있습니다. 아래는 계산된 수평 속도 분포를 나타내고 있습니다. 이것은 표면 바로 위에 제어 볼륨 층의 계산 된 수평 속도 분포를 제공하는 X-Y 플롯에 그래픽으로 보여 주고 있습니다. 격자 미세화에 의해 표면의 고체 부분의 윗쪽 속도가 영이 되는 경향이 있지만, 슬롯들 위에 있는 속도는 영이 안되게 유지됩니다. 많은 요철 위의 이러한 속도의 평균은 효과적인 슬립으로 해석 될 수 있는 non-zero 수평 이송 속도를 일으킵니다.
Evaporative Effects
분산된 고체 물질을 포함하는 액체 방울이 고체 표면에서 건조되면 고체 물질이 침전물로 남습니다. 이 퇴적물의 패턴은 많은 인쇄, 청소 및 코팅 공정에 중요한 의미를 갖습니다. 한 가지 유형의 침전물의 전형적인 예는 왼쪽 이미지와 같이 유출 된 커피 조각의 둘레를 따라 링 얼룩이 형성되는 “커피 링”문제입니다. 이러한 유형의 링 침전물은 액체의 증발로 인한 표면 장력 구동 흐름의 결과로 발생하며, 특히 방울 주변에서 발생합니다 [1].
Drying
건조는 코팅 공정의 중요한 부분입니다. 잘 도포된 코팅은 건조 결함으로 인해 완전히 손상될 수 있습니다. 건조 중에 온도 및 용질 구배는 밀도 및 표면 장력 구배로 인해 코팅 내 흐름을 유도 할 수 있으며, 이로 인해 잠재적으로 코팅 품질이 손상 될 수 있습니다. FLOW-3D의 증발 잔류물 모델을 통해 사용자는 건조로 인한 흐름을 시뮬레이션하고 값 비싼 물리적 실험에 소요되는 시간을 줄일 수 있습니다.
FLOW-3D’s evaporation residue model simulates a 3D view of residue formed from toluene after drying (magnified 30x)
Modeling Ring Formation
FLOW-3D는 증발이 가장 큰 접촉 라인에서의 증착으로 인해 에지 고정이 발생 함을 보여줍니다.
링 형성 모델링 증발에 의해 접촉 라인에서 생성 된 흐름 시뮬레이션 증발은 증발로 인한 열 손실로 인해 액체를 냉각시킵니다 (색상은 온도를 나타냄). 동시에 고체 표면은 전도에 의해 액체를 가열합니다. 증발은 접촉 라인 근처에서 가장 크므로 액체가 접촉 라인을 향해 흐르게하여 정적 상태를 다시 설정합니다. 최종 결과는 액체가 완전히 증발하는 액체 가장자리에 부유 고체가 증착됩니다.
FLOW-3D의 접촉 선 고정 모델에 대해 자세히 알아보십시오.
Simulation of flow generated at a contact line by evaporation