Paradigm Change in CFD Led by Open-Source Codes 이신형교수서울대학교공과대학조선해양공학과 CAE Conference 년 11 월 15 일 ( 금 ) COEX

Similar documents
82 제 1 발표장 (2 일금 ) CFD 응용 [V] 이남훈 1*, 류태광 2 NUMERICAL VERIFICATION OF SHAKE TABLE TEST FOR THE LIQUID STORAGE TANK N. Lee and T. Yoo 1.,.,,,.,. Baek e

00-CD..

03_전체논문취합_전체(추계)_수정본.hwp

PowerPoint 프레젠테이션


歯J PDF

Getting Started

, Yard Bottom Slamming, Slamming,, 10-8 Probability Level Bottom Slamming., Bottom Slamming,, Evaluation, Allowable Criteria, Ballast Reduction, Botto

PowerPoint 프레젠테이션


소성해석

NX Nastran brochure (Korean)

<31325FB1E8B0E6BCBA2E687770>

Precipitation prediction of numerical analysis for Mg-Al alloys

Microsoft Word - KSR2012A103.doc

2. 수치시뮤레이션 2.1 기본방정식과수치조건 기본방정식은 Navier-Stokes 방정식이며 FEM 수치기법으로이산화하여구조격자를만들어계산을수행하였다. k- 을사용한수송방정식은 t (ρε)+ (ρεu x i )= i x j [( μ+ μ t σ ε ) ε + C 1ε

7 LAMPS For use on a flat surface of a type 1 enclosure File No. E Pilot Lamp File No. E Type Classification Diagram - BULB Type Part Mate

¿ÀǼҽº°¡À̵å1 -new

김우전5.hwp

femap brochure (Korean)

Automation for STAR-CCM+ Kim Hak-Soon

Your title goes here

슬라이드 1

< DC1A4C3A5B5BFC7E22E666D>

Turbine Digital Flowmeter SEMI U+ 특징 PVC, PTFE, P.P, PVDF 등 다양한 재질 Size, 유량, Connection별 주문제작 정밀성, 내화학성이 우수 4~20mA, Alarm, 통신(RS485) 등 출력 제품과 Controll

ETL_project_best_practice1.ppt

untitled


Microsoft Word - KSR2013A320

Microsoft Word - KSR2012A021.doc

(2 일금 ) 제 1 발표장 Opensource CFD toolkit OpenFOAM, SNUFOAM. (1) (2).,,,,. (3) (1), (2) RANS (Reynolds-averaged Navier-Stokes) SST. - PIMPLE (merg

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 10, Oct ,,. 0.5 %.., cm mm FR4 (ε r =4.4)

(2) : :, α. α (3)., (3). α α (4) (4). (3). (1) (2) Antoine. (5) (6) 80, α =181.08kPa, =47.38kPa.. Figure 1.

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Feb.; 28(2),

2 Journal of Disaster Prevention


.4 편파 편파 전파방향에수직인평면의주어진점에서시간의함수로 벡터의모양과궤적을나타냄. 편파상태 polriion s 타원편파 llipill polrid: 가장일반적인경우 의궤적은타원 원형편파 irulr polrid 선형편파 linr polrid k k 복소량 편파는 와 의

Microsoft Word Á¤Áظð.doc

Vertical Probe Card Technology Pin Technology 1) Probe Pin Testable Pitch:03 (Matrix) Minimum Pin Length:2.67 High Speed Test Application:Test Socket

03.hwp

<C1A4C3A5BAB8B0EDBCAD2D D30355F33B1B32E687770>

< FBEC8B3BBB9AE2E6169>

GNU/Linux 1, GNU/Linux MS-DOS LOADLIN DOS-MBR LILO DOS-MBR LILO... 6

Microsoft PowerPoint - chap01-C언어개요.pptx

00829A_SHR-6164-KOR.indb

À±È°Ãß°è½ÉÆ÷ÆÊÇ÷¿ 59ȸ-ÃÖÁ¾

Introduction Capillarity( ) (flow ceased) Capillary effect ( ) surface and colloid science, coalescence process,


88 KOREA INSTITUTE OF LOCAL FINANCE

108 KOREA INSTITUTE OF LOCAL FINANCE

82-01.fm

Microsoft Word - KSR2012A172.doc


untitled

[ 물리 ] 과학고 R&E 결과보고서 유체내에서물체의마찰력에미치는 표면무늬에대한연구 연구기간 : ~ 연구책임자 : 홍순철 ( 울산대학교 ) 지도교사 : 김영미 ( 울산과학고 ) 참여학생 : 김형규 ( 울산과학고 ) 노준영 (

한약재품질표준화연구사업단 단삼 ( 丹參 ) Salviae Miltiorrhizae Radix 생약연구과

2017 년 6 월한국소프트웨어감정평가학회논문지제 13 권제 1 호 Abstract

I I-1 I-2 I-3 I-4 I-5 I-6 GIS II II-1 II-2 II-3 III III-1 III-2 III-3 III-4 III-5 III-6 IV GIS IV-1 IV-2 (Complement) IV-3 IV-4 V References * 2012.

삼성955_965_09

03_ok.hwp

M&S. EDISON. e-science EDISON_CFD[2] EDISON. (geometry kernel software),,, EDISON_CFD.,.. 4. e-science EDISON M&S.., EDISON , b f. Find to min

14(4) 09.fm

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 28, no. 3, Mar guidance system: MGS), MGS. MGS, (pulse repetit

¾Ö´º¾ó¸®Æ÷Æ®(2010)1.ps

Microsoft PowerPoint - 3.공영DBM_최동욱_본부장-중소기업의_실용주의_CRM


논단 : 제조업 고부가가치화를 통한 산업 경쟁력 강화방안 입지동향 정책동향 <그림 1> ICT융합 시장 전망 , 년 2015년 2020년 <세계 ICT융합 시장(조 달러)> 2010년 2015년 2020년 <국내 ICT


백서2011표지

VOL /2 Technical SmartPlant Materials - Document Management SmartPlant Materials에서 기본적인 Document를 관리하고자 할 때 필요한 세팅, 파일 업로드 방법 그리고 Path Type인 Ph

(1 일목 ) 제 3 발표장 47 수치기법 [I] 이은택 1, 안형택 2* SIMULATION ON FLOW PAST A CIRCULAR CYLINDER USING UNSTRUCTURED MESH BASED INCOMPRESSIBLE FLUID SOLVER(ULSAN3

Microsoft Word - 1-차우창.doc

ePapyrus PDF Document

Sun ONE Portal Server, Mobile Access, h 6.2 ƒe 1 Ï Û ( Ñ ) l d e ' f el d ' f f i. 2 f CPU d Šf th. l hh Š Š h h Š. l hh Š f f hšš. l hh j j l

슬라이드 1

Microsoft PowerPoint - 7-Work and Energy.ppt

hapter_ i i 8 // // 8 8 J i 9K i? 9 i > A i A i 8 8 KW i i i W hapter_ a x y x y x y a /()/()=[W] b a b // // // x x L A r L A A L L A G // // // // /

주요국 에너지 Profile 분석_아랍에미리트

What's new in NX 8 (Korean)

67~81.HWP

<303520B1E2BABBBBE7BEF7BAB8B0EDBCAD2E687770>

Microsoft Word - KSR2012A038.doc

PowerPoint 프레젠테이션

<C3CA3520B0FAC7D0B1B3BBE7BFEB202E687770>

SEMI U+

LEICA C-LUX

Microsoft Word - KSR2013A299

슬라이드 1

example code are examined in this stage The low pressure pressurizer reactor trip module of the Plant Protection System was programmed as subject for

H122_E.book

슬라이드 1

ESP1ºÎ-04

13 Who am I? R&D, Product Development Manager / Smart Worker Visualization SW SW KAIST Software Engineering Computer Engineering 3

1. Features IR-Compact non-contact infrared thermometer measures the infrared wavelength emitted from the target spot and converts it to standard curr

Oracle Apps Day_SEM

<4D F736F F F696E74202D204D5249B1E2BCFAB5BFC7E2B9D7BDC3C0E5C7F6C8B25F F325B315D205BC8A3C8AF20B8F0B5E55D>

TTA Verified : HomeGateway :, : (NEtwork Testing Team)

12(4) 10.fm

< C6AFC1FD28C3E0B1B8292E687770>

<353020B9DAC3E1BDC42DC5ACB6F3BFECB5E520C4C4C7BBC6C3BFA1BCADC0C720BAB8BEC820B0EDB7C1BBE7C7D7BFA120B0FCC7D120BFACB1B82E687770>

Transcription:

Paradigm Change in CFD Led by Open-Source Codes 이신형교수서울대학교공과대학조선해양공학과 CAE Conference 2013 2013 년 11 월 15 일 ( 금 ) COEX

CFD 패러다임의변화 In-house codes - Unstructured mesh - Pre/post integration - CAD import - GUI - General purpose soler Commercial codes - Community deelopment - OOP - Customization capability CFD for open innoation - Huge computing resources requirements - License fees - Monopoly - Difficulty in customization/integration Open-source codes 2

OpenFOAM 이란? OpenFOAM? Open Field Operation And Manipulation Open Source CFD toolbox 1990 년대초영국의 Imperial College 에서개발시작 주저자 Dr. Hroje Jasak : Wikki Ltd. Dr. Henry Weller : OpenCFD Ltd. 2004 년말 OpenCFD Ltd. 에서 1.0 이공개 S/W 로발표 3

OpenFOAM 특징 GNU GPL 을적용한 open source software 어플리케이션 (Application) 이라고칭하는실행파일을만들기위한 C++ 라이브러리들의모음 라이브러리들은해석하고자하는편미분방정식과최대한유사하게만들어져있음 소스와함께미리컴파일된다양한어플리케이션이함께제공됨 4

OpenFOAM 구성 OpenFOAM 은 Application 이라는실행파일을만들기위한 C++ 라이브러리모음 Open Source CFD Toolbox 각종 Standard Application 이포함 CFD, Molecular dynamics, Electromagnetics, Solid dynamics, Finance Application 은 Soler 와 Utility 로구분 Soler 는연속체역학계산을위한프로그램 Utility 는데이터처리등의작업에필요한프로그램

OpenFOAM Library FiniteVolume Mesh GenericPatchFields DynamicMesh FMotionSoler Parallel postprocessing OpenFOAM TurbulenceModels TransportModels ThermoPhysicalModels lagrangian CombustionModels Engine ODE regionmodels

OpenFOAM Applications Solers Basic Combustion Compressible Incompressible HeatTransfer Multiphase DNS Utilities Mesh ParallelProcessing PostProcessing PreProcessing Surface Thermophysical Miscellaneous

OpenFOAM 호환성 전처리프로그램 OpenFOAM 포맷지원 ü Harpoon, Gridgen, Pointwise, Engrid 격자파일호환 StarCCM+, Fluent, Plot3D 후처리프로그램 FieldView, Ensight, Tecplot MpCCI 8

OpenFOAM 라이선스 OpenFOAM 은등록된 trademark OpenFOAM 의라이선스는 GNU GPL GNU : GNU is Not Unix GPL : General Public License Http://www.gnu.org/licenses/gpl.html 9

OpenFOAM 의장점 사용자및 CPU 수에따른라이선스비용부담감소 하드웨어인프라대규모확장가능 계산시간감소, CFD 적용분야확대 전용프로그램제작 계산과정의단순화 CFD에대한접근성확장, 해의정확성보장 기능추가용이 CFD personnel 의효율극대화 10 Experts lead/guide deelopment Field engineers use codes/packages

OpenFOAM 의단점 OS 의제한 ( 리눅스 ) 계산용클러스터장비로문제해결 윈도우용제작가능 TUI 작업환경 GUI 제작으로극복가능 어려운사용방법 매뉴얼의한계, 기술지원조직의부재 정보및기술공유의어려움 저변확대, 커뮤니티형성, 국내행사마련 11

국내현황 교육 넥스트폼 : 기본교육및고급교육 저자초청세미나 2009 년 Dr. Hroje Jasak 초청세미나 ( 서울대조선해양공학과 ) 2011 년 Dr. Hroje Jasak 초청세미나, 150 명참석 ( 넥스트폼, 서울대조선해양공학과 ) 2012 전산유체공학회춘계학술대회 오픈폼세션 중소기업청산학연공동기술개발사업 12 선박해석용다상유동시뮬레이션프로그램개발 넥스트폼-서울대조선해양공학과

Audi, Volkswagen, Volo 해외적용사례 External aerodynamics of car NIST (National Institute of Standard and Technology) / DOE firefoam / FDS NREL (National Renewable Energy Laboratory) / IEA (International Energy Agency) LES wind turbine / wind farm analysis code DLR, NASA, US Nay Bombardier, ABB, Basf, Hitachi, Man, PBMR, Shell, Siemens, DNV, GL 13

주요행사 국외 International OpenFOAM Workshop ü 2006년부터매년개최 ü 크로아티아, 이태리, 카나다, 스웨덴, 미국, 독일 ü 8th Workshop : 2013년 6월 11 ~ 14, 제주, 한국 Open Source CFD Conference 매년개최 2012년 10월 29 ~ 30, 런던, 영국 국내 14 OKUCC (OpenFOAM Korea Users Community Conference) ü 2012년대전 KISTI : 약 70명참가 ü 2013년서울대학교호암교수회관 : 약 70명참가 ü 2014년부산개최예정

SNUFOAM 소개 SNUFOAM OpenFOAM 을기반으로함 조선해양유체분야에특화되어있는맞춤형코드 ü 난류유동 ü 자유수면유동 ü 캐비테이션유동 ü 수치수조 ü 부유체 6자유도운동 여러모듈들의지속적인개발및업그레이드 15

SNUFOAM-벽함수 벽함수에따른난류유동장의차이점분석 WF-1: 난류의생성과소멸이같다는가정으로부터유도 ü OpenFOAM 에서제공하는라이브러리 G k = t w C 1/ 4 m k ky P 1/ 2 WF-2: 난류의생성과소멸이같다는가정을포함하지않음 ü OpenFOAM 에서제공하지않는라이브러리로새로개발하여코드에삽입 G k =t w t / r ky C k w 1/ 4 P m 1/ 2 16

SNUFOAM- 벽함수 KCS 선수부주위의난류유동장비교 p k epsilon WF-1 WF-2 17

SNUFOAM- 벽함수 상용코드인 Fluent 와비교 상용코드와유사한결과를보임 p k epsilon SNUFOAM with WF-2 Fluent 18

자유수면기법의추가 SNUFOAM- 자유수면 보통자유수면을처리하기위하여 VOF (olume of fluid) 기법을사용 t t ò W ò W ò r udw + n ruuds = n ( m Ñu) ds - ò S ò OpenFOAM 에서는 compressie 기법의일종인 interface compression 만을제공 S ü 스킴의특성상비물리적인자유수면의형상이나타날수있음 eff ò W ÑpdW + adw + a n uds = 0 0 a 1 S ò skñadw W 19 C=0 C=0.5 C=1.0 CICSAM 이나 HRIC 등의자유수면을처리할수있는기법이필요

SNUFOAM- 자유수면 SNUFOAM 의고해상도자유수면포착기법 Compressie Interface Capturing Scheme for Arbitrary Meshes (CICSAM) - Ubbink (1997) ~ α ~ ~ f = γf αfhyperc + (1- γf ) αfuq, 0 γf 1 ì 1+ cos 2q f ü = miní, 1ý î 2 þ High-Resolution Interface Capturing (HRIC) - Muzaferija and Peric (1999) ~ a ~ α f f = g ~ * (1 ) ~ f a + - g f a D, 0 g f ì ~ αd ï * = ~ í2α ï î1 D : 0 > ~ α ~ D, αd > 1 : 0 ~ αd 0.5 : 0.5 ~ α 1 D 1 γ f γ f = cosθ f 1/2 20 Modified CICSAM (M-CICSAM) - Waclawczyk and Koronowicz (2008) ~ α = γ ~ ~ f αfcbc + (1- γf ) αffromm, 0 γf 1 1/4 γ f = cosθ f f Modified HRIC - Rhee (2009) ~ a f = g ~ * (1 ) ~ f a f + - g f a fuq, 0 g f 1 γ f = cosθ f 1/2

SNUFOAM- 자유수면 댐붕괴문제 시간에따라물기둥이붕괴하면서자유수면의형상이변함 실험과잘일치하고있음 21 Shape of free-surface interface Vertical displacement of free-surface interface (Top) Horizontal displacement of free-surface interface (Bottom)

SNUFOAM- 자유수면 슬로싱문제 슬로싱으로인한바닥과벽면에서의시간에따른압력변화를실험과비교 계산결과가실제물리적현상을잘표현하고있음 Pressure at P1 22 Pressure at P2 Pressure at P3

SNUFOAM- 캐비테이션 함정이나선박의타혹은추진기에서발생하는캐비테이션유동 일반적으로물은섭씨 100 도이상에서끓기시작하지만압력이증기압보다낮아질때에도발생 함정이나선박의타혹은추진기에서가장많이발생 캐비테이션에의한방사소음을줄이기위해서는캐비테이션해석이필수적 23

24 SNUFOAM- 캐비테이션 캐비테이션모델 Two phase mixture model s eap cond R R R E t Q + + = + Ñ ) ( Q E R cond R eap R s Merkle et al. 0 Kunz et al. 0 Schnerr and Sauer 0 Singhal et al. f m f r ø ö ç ç è æ ø ö ç è æ - - U L U f P P C ch l cond 2 2 1 r ( ) ø ö ç ç è æ ø ö ç è æ - - U L U f P P C ch l eap 2 2 1 1 r a m a r ( ) ø ö ç ç è æ - - U L C ch cond 2 1 a a ( ) ø ö ç ç è æ ø ö ç è æ - - - U L U P P MIN C ch l eap 2 2 1 1,0) ( r a r a m r r a ( ) l m l cond P P R C r a a r r r - - - 3 2 3 1 ( ) l m l eap P P R C r a a r r r - - 3 2 3 1 m f r m m f r r l l l ch cond f P P C r r r g - - 3 2 ( ) l l ch eap f P P C - - 1 3 2 r r r g ø ö ç ç è æ Ñ Ñ t f s m

SNUFOAM- 캐비테이션 Sheet caitation Hemispherical headform body

SNUFOAM- 캐비테이션 Super caitation Wedge

SNUFOAM- 캐비테이션 Cloud caitation Delft twisted foil time

캐비테이션의압축성효과 SNUFOAM- 캐비테이션

SNUFOAM- 캐비테이션 캐비테이션의압축성효과 : sheet caitation

SNUFOAM- 캐비테이션 Propeller caitation and erosion KP505 J = 0.3

SNUFOAM- 수치수조 파를생성, 소멸시키기위한기법 파랑하중에의한선박이나해양구조물의영향을알아보기위하여필요 조파와소파를위하여운동량생성항기법을사용 t ò W rudw + ò S n ruuds = ò S 생성된파가해석해와잘일치하고있음을확인 ò n ( m Ñu) ds - ÑpdW + skñadw + b ( u - u eff W ò W ò W I ) dw Wae profile (T: 15, 10, 5(s), H:2m) Wae profile (T: 15, 10(s)), H:4m)

SNUFOAM-6DOF 선박이나해양구조물의 6 자유도움직임을위한라이브러리 격자변형기법등을사용하여부유체의운동을해석할수있음 Pitch motion only Surge, heae and pitch motion

적용사례 - KCS KCS 의기본제원 3,600 TEU container carrier KCS specification Designation Prototype Model Scale ratio 1 1/31.6 Speed (m/s) 12.3467 2.1964 Froude number (Fr) 0.26 0.26 Reynolds number (Re) 2.4 10 9 1.4 10 7 Length (m) 230.0 7.2786 Breadth (m) 32.2 1.0190 Depth (m) 19.0 0.6013 Draft (m) 10.8 0.3418 Wetted surface area (m) 9,498.0 9.5121 Displacement (m) 52,030.0 1.6490 Block coefficient (C B ) 0.6505 0.6505

적용사례 - KCS 저항및자항계산 Soler : SNUFOAM KCS resistance test KCS self-propulsion test with KP505 propeller

적용사례 - KCS 자유수면분포해석 경계면포착기법에따른선체표면에서자유수면의분포를비교 ü 선체표면에서자유수면의높이가실험과비교적잘일치하고있음 Wae pattern 비교 ü 상용코드와비교하여유사한결과를보임 0.02 0.01 Data Interface Compression SuperBee anleer modified HRIC modified CICSAM Fluent (RKE) z/lpp 0-0.01-0.02-1 -0.8-0.6-0.4-0.2 0 x/lpp SNUFOAM (RKE) Hull wae profile Wae pattern

적용사례 DTMB 5415 DTMB 5415 (DDG-51) 의기본제원 US Nay Combatant, DTMB 5415 DTMB 5415 specification Scale Ratio 24.832 Length (L) Draft (T) 5.72 m 0.248 m Wet Surface Area (S) 4.861 m 2 Adance Velocity (U 0 ) 2.09 m/s Froude Number (Fn) 0.28 Reynolds Number (Re) 1.26 10 7 Sonar dome과 transom stern을가짐 Twin open-water propeller를가짐 Shaft는 struts에의해지지됨 Full scale 함정은존재하지않음 (DDG-51) Sinkage at FP Sinkage at AP -0.0028L -0.0009L

적용사례 DTMB 5415 선체표면에서의압력분포 Sonar dome 의전면에서정체압이걸리는것을확인 Sonar dome 의측면에서는 sonar dome 을타고흐르는빠른유동으로인하여압력이급격히낮아짐 Pressure coefficient

적용사례 DTMB 5415 자유수면형상의비교 상용코드인 FLUENT* 와비교 함정주위의파형을잘나타내고있음 * Zdraistch, F., Hydrodynamics Analysis of Ships Side by Side in Waes using AQWA and Resistance and Diffraction Simulation oer a Ship Hull using ANSYS-CFD, The Society of Naal Architects & Marine Engineers, Report of Texas Section Technical Meeting, Houston, 2008. Wae pattern

적용사례 활주선 활주면주위의유동해석과양항력비교 활주선을단순화시킨형상 양력, 항력, 모멘트를실험및 Saitsky s formular 와비교 Lift Drag Moment

적용사례 활주선 고속활주선에적용하여계산 활주선의속도에따라자세가달라짐 15knot 30knot

감사합니다