#OpenFOAM #우분투 #Linux #window 코드 작업 환경 설정 자 일단 Setup하기 위해서 몇가지가 필요한데요. 저 같은 경우는 Window 환경에서 작업을 진행하고 있습니다. OpenFoam은 Linux 환경에서 쓰는게 편하기 때문에 Ubuntu or PowerShell을 설치해주는게 필요합니다. 위 두 프로그램은 우리가 Window 환경에서 Linux 프로그램을 쓸 수 있도록 하게 됩니다. 그 다음으로는 어떠한 편집기를 쓸건지를 결정해야합니다. --> 코드 작업용 저 같은 경우는 그냥 Linux 환경이 익숙하지 않아서 Window 환경에서 사용하는 편집기를 그대로 이용합니다. 그러기 위해서는 그냥 Linux 환경이 설치 된 파일 경로만 알고 있으면 됩니다. 다음과 같이 파일 경로창에 입력을 해주세요 \\wsl.localhost\ 그러면 저는 Ubuntu를 설치 했기 때문에 Ubuntu 폴더가 보입니다. 들어가서 home/user 파일 경로에 들어가시면 제가 설치 해

#OpenFoam #문법 #wmake C++을 이용하게 된다면 Build라는 작업은 필수적으로 발생합니다. (파이썬, 매트랩과 다른부분이죠) 이후에 우리가 작성한 코드를 실행할 수 있게됩니다. Build 먼저 빌드라는 것은? 소스 코드로 부터 실행 파일을 만드는 작업입니다. 이때 필요한 것이 컴파일러라는 도구입니다! 자 코드를 짜다 보면 우리는 이전에 작업 했던 것을 가지고 와서 내 코드에 적용하고 싶죠.. 이때 코드에는 # include header를 하게 됩니다. 그렇게 되다보면 Header 위치에 대한 정보를 줘야합니다. 다만, 보통의 header에는 내용을 포함하고 있지는 않습니다. --> 내부 내용을 비공개하기 위해서 일부는 Header만 제공하고, 내부 내용은 Library라는 것을 통해 제공합니다. 그러면 이제 Library에 대한 위치를 알려줘야 하죠!! 즉, Header와 Libary 위치에 대하여 정보를 컴파일러에 줘야합니다!! 그렇게되면 실행코드를 만들 수 있는

Open Foam은 C++ 기반으로 작성해서 객체를 생성하고, 그 객체에 어떠한 Data를 넣어서 Solver를 세팅하는 구조 입니다. 먼저 C++ 에서 객체 생성은 다음의 법칙을 따릅니다. BaseClass* instance = new DerivedClass(); 여기서 *은 Pointer로 C++에서의 Source Code의 용량및 메모리 관리를 하기 위해 사용하는 방법입니다. 다음으로는 C++에서 autoPtr 또는 smatPtr을 사용하는데요, 이는 기존 Pointer가 가지고 있는 단점을 개선하기 위해서 C++에서 내놓은 좋은 방법입니다. --> 코딩을 하지 않는 다면 굳이 알필요는 없습니다. 뭐 다 몰라도 되지만, 쉽게 말씀 드리면 그냥 어떠한 Solver의 객체를 쉽게 생성할 수 있는 방법이 존재한다. 그럼 Turbulence의 객체를 만들어 볼까요? autoPtr tbrptr ( incompressible::t

#대류 #Convection #Prandtl #Nusselt #경계층 열의 대류(Convection)와 경계층 이론 이번 포스팅에서는 열전달의 두 번째 메커니즘 — 대류(Convection) 에 대해 다룹니다. 전도(Conduction)가 물질 내부의 분자 운동에 의한 열전달이라면, 대류는 유체의 거동(속도)에 의해 에너지가 이동하는 현상이에요. 1. 대류(Convection)란? 즉, 유체의 운동과 함께 열전달이 일어나는 것이 이 대류에 해당합니다. 이 이동현상에는 2가지가 존재하는데요. 구분 설명 예시 자연 대류 (Natural Convection) 밀도 차이에 의한 부력(Buoyancy)이 원인 공기 중의 히터 주변 상승 기류 강제 대류 (Forced Convection) 외부 힘(팬, 펌프 등)에 의해 유체가 이동 열교환기, 팬 냉각 장치 등 바로 자연대류와 강제 대류입니다. 자연 대류로는 밀도차 --> 압력차라 보시면 되겠죠. 차가운 공기는(밀도가 큰) 아래로 가라 앉고

#전도 #CFD #StarCCM+ #열해석 #Conduction #푸리에수 #Biot 수 #FO 전도(Conduction) 열전달 이론 Simcenter STAR-CCM+ Heat Transfer 중 전도(Conduction) 파트를 정리해봅니다. 전도는 열이 온도 구배(Temperature Gradient)에 의해 이동하는 대표적인 열전달 메커니즘입니다. 1. 전도(Conduction)의 기본 이론 전도는 고체 내부나 정지된 유체에서 분자 진동이나 자유전자 이동으로 열이 전달되는 현상입니다. 전도의 지배방정식은 Fourier의 열전도 법칙으로 표현됩니다. q = -k ∇T q : 열유속 벡터 (W/m²) k : 열전도율 (W/m·K) ∇T : 온도 구배 이 경우는 등방성에 해당 되는데요, Star CCM+에서는 k를 Tensor로 제공하여 이방성 까지 해석이 가능하도록 제공하고 있습니다. 2. 무차원 수 (Fourier, Biot)와 Time Step 결정 비정상(Transien

#파이썬 #가상환경 #설정 #Python으로만 해당 글은 window 환경에서 Python만 가지고 가상환경을 만들기, 왜 만들어야 하는지 다루는 글입니다. 먼저 가상환경을 만드는 방법은 간단합니다. Python을 이용하여 가상환경 구축하기 일단 venv 라는 모듈이 깔려 있어야 합니다. 없다면 pip install venv를 통해 간단히 설치해주면 됩니다. python -m은 module을 실행하기 위한 명령어입니다. 그렇다면 test라는 가상환경을 만들기 위해서는 간단히 다음과 같이 명령어를 입력해주시면 됩니다. python -m venv test 그러면, 우리는 test라는 파일이 생긴것을 볼 수 있습니다. 그러면 그안에 activate라는 batch 파일이 생긴것을 볼 수 있죠? Window에서는 이 부분만 실행해주면 됩니다. cd test/Scripts activate 이렇게 해주게 되면 스크립트 폴더 내에 있는 activate를 실행시킬 수 있게됩니다. 그러면, 가상환경이

#유동층 #반응기 #Geldart C #거동 #Channeling #응집체 #프랙탈 #구조 Geldart Class 분류 Geldart Class는 밀도, 직경에 따라서 분리를 합니다. A(Aeratable)는 공기처럼 유동이 잘되는 Class, B(Bubbling)는 유동화 할때이 같이 일어나는 Class, D는 Spout (터지는) 현상이 발견되는 Class로 구분이 되죠. C(Cohesvie)는 응집력인 강한 Class에 해당해서 유동이 잘 안되는 특징을 가지고 있습니다. Geldart Class 분 이 중 C 클래스는 Cohesive Group으로 분류되며 이는 유동화 될때, 내부 정전기적 인력이 강하여 유동화가 잘 안일어 나죠. 대표 현상으로는 다음과 같습니다. 현상명 설명 Channelling 유체가 저항이 적은 통로로만 흐름, 다수의 입자 비활성화 Slugging 응집 덩어리가 유동화되지 않은 채 덩어리째 상하로 움직임 Defluidization 유동화가 유지되지 못하고

#Data #Driven #Fluid #POD #Proper #Orthogonal #Decomposition 안녕하세요. 화공 공대생입니다. Data Driven Fluid Mechanis (Miguel A. Mendenz) 책을 리뷰를 시작하고자 합니다. Data-Driven Fluid Mechanics - 예스24 Big data and machine learning are driving profound technological progress across nearly every industry, and are rapidly shaping fluid mechanics r... www.yes24.com 이 책은 유체역학에서의 데이터 기반 접근법과 머신러닝 기법을 결합하여, 모델 차수 감소, 시스템 식별, 유동 제어, 난류 모델링 등 다양한 주제를 다루고 있습니다. 특히, 전통적인 물리 기반 모델링과 최신 데이터 기반 기법의 통합을 강조하며, 이론과 실습을 균형 있게 제공하고 있어서

#화공공대생 #OpenFoaM #BuoyantSimpleFoam 오늘은 온도 변화에 따른 유동에 대해서 알아볼 예정입니다. 이전 시간에는 온도 변화가 기체의 밀도에 영향을 준다고 했죠? PV=nRT로 부터 쉽게 알 수 있을 것입니다. 이로 인해서 뜨거운 공기는 위로 차가운 공기는 아래로 이동하게 됩니다. 이를 자연대류라고 하죠. 자 OpenFoam에서는 Simple Foam을 제공해서 이를 해석이 가능합니다. <- 이게 뭐가 문제죠? 라고 생각 하실 수 있는데, Simple Foam은 비압축성 해석에 사용되는 거죠! 아 그러면 밀도 변화가 없는데...? 어떻게 구현한다고? >> 사실 우리는 이러한 해답을 알고 있습니다. 대충 온도 변화에 따라서 유동을 주면 되겠구나! 라고요!! 자 먼저 Simpe Foam에 대한 수식을 볼까요? ∇∙u=0 ∇∙(u⊗u)−∇∙R=−∇p+Su 전형적인 압력 값을 통한 유속 변화를 예측하는 수식이 됩니다. 여기서 u는 유속, R은 Stress Tensor

#유체역학 #부력 #대류 #Boussinesq #부시네스크 대류현상 공기의 대류는 차가운 공기는 아래로 뜨거운 공기는 위로 올라 갑니다. 이를 대류 현상이라고 합니다. 이 원리는 간단합니다. 바로 밀도 차에 의해서 이러한 혼합이 발생하는 것이죠. 하지만 우리는 온도차에 의해서 뜨거운 공기는 위로 찬공기는 아래로 가라앉는다고 이해하고 있죠. 이 부분에 대해서 오늘 좀 알아보려고 합니다. 기체가 이상기체 방정식을 따른다고 보죠. 그러면 우리는 밀도는 온도에 반비례한다고 알 수 있습니다. 아 Ok! 뜨거운거는 가벼운 공기고, 차가운 공기는 무거운 공기야!! 이 부분이 온도차에 의해서 밀도차가 발생하게 되고 이에 의해 혼합이 발생하는 것입니다. 자 그러면? 밀도차에 의해서 왜 혼합이 발생하는가를 보죠. 이를 이해하려면 부력을 이해해주어야 합니다. 부력 이란 우리 근처에서 쉽게 발견할 수 있습니다. 소금물에 몸이 뜨거나 아니면 배가 바다에 뜨는 현상이 바로 부력입니다. 이에 대한 정의를 해보

#난류 #모델링 #RANS #LES #DNS #openFoam #Cylinder #화공공대생 #k-epsilon #k-omega #SST 이전 시간에 이어서 오늘은 Cylinder 주변에서 발생하는 난류와 박리에 대한 시뮬레이션을 해보겠습니다. 이에 대한 이론적인 내용은 중요하니 따로 다루겠고, 이번에는 어떻게 모델링하고 시뮬레이션 하는 방법에 대하여 알아보겠습니다. 난류 모델링 개요 먼저 난류모델링은 3가지고 구분해서 볼 수 있습니다. 1. 직접푸는 DNS( Direct Numerical Simulation) 2. 공간 평균 방법 LES(Large Eddy Simulation) 3. 시간 평균 방법 RANS(Reynolds Averaged Numerical Simulation) 1 -> 3으로 가는 순서로 컴퓨팅 시간이 감소하며, 디테일 적인 요소는 줄어듭니다. DNS는 직접 푸는 방법에 해당한다 정도만 알고 넘어가면됩니다. LES는 시간에 따른 축적이 필요 없으며, 동일 시간대

#Solid #Stress #Shear Stress 해당 포스팅은 Stress, Shear stress, Solid Stress in mulitphase에 대하여 다루고 있습니다. Stress란? deformation (변형)이 발생하는 동안의 힘을 의미합니다. 힘은 단위 면적당 힘을 의미하며 압력과 동일한 단위 차원을 쓰죠! (stress=F/A) 변형을 보게되면 아래와 같게 되는데 tension Compression의 경우 면에 수직에 작용하는 것이 되고 해당 shear의 경우 면에 tangential 하게 작용하는 힘이 됩니다. shear와 관련 shear stress란? shear라는게 비틀림을 의미합니다. 지구과학 시간에 보면 막 지진처럼 끊어졌던 것들 있죠? 그런 것이 바로 shear stress에 의한 현상입니다. 이제 연속성 관점에서 보게 된다면 유체에서 정의하는 전단 응력이 됩니다. 이때 유체는 연속적으로 변형이 이뤄지는 것으로 보는것이고, 이 변형에 대한 기준은 축

#난류 #CFD #모델링 #전달현상 해당 내용은 아래 강의 내용을 정리한 내용입니다. 자세한 내용은 아래 내용을 참고 부탁드립니다. https://elearning.edison.re.kr/cfd/courseware/c4/index.html?id=drg_1_1 난류이론 층류유동 및 난류유동 강의 내용 소개 층류-천이-난류:시각적 이해 경계층의 정의 난류 경계층의 발달 난류 경계층의 속도분포 난류유동 특성 유동의 전달현상 골프공 dimple 익형 항력 경계층 이론 적분해석 압력분포의 영향 경계층 이론의 결과 난류 경계층에 대한 몇 가지 사실 난류의 생성 과정 난류 운동에너지와 소산율 난류경계층 기초 이론-벽법칙 난류경계층 기초 이론-몇 가지 유용한 결과 난류경계층 해석 Tip: 격자 생성 난류 이론 이해를 위한 핵심 기초난류이론 마무리 퀴즈 퀴즈 1 퀴즈 2 퀴즈 3 난류모델링 강의 내용 ... elearning.edison.re.kr 전달현상이란? 유동이 확산되어 혼합하는 과정이라고

#OpenFoam #Mixer #MRF #Sliding Method #Coriolis force. OpenFoam에서 교반기를 모사하기 위해서 MRF 방법을 이용한접근을 해보고자 한다. 교반기란 먼저 교반기에 대해서 알아보자. 교반기 액체를 휘저어 물질을 섞는 장치이다. 교반기 구성 Motor, Gear Box, Mount Plate, Shaft, Impeller, Brushing 여기서 Impeller가 우리가 중점적으로 봐야할 것에 해당한다. Impleller Type 임펠러는 Axial 방향(수직전단)에 대한 것과 Radial 방향(분사형)에 대한 것으로 나뉠 수가 있다. 수직 전단의 경우 유체를 아래로 밀어내주어 발생하는 유체의 흐름으로 Down Up Flow가 발생하게 된다. 수평 전단의 경우는 유체를 옆으로 밀어내어 방사형태의 유체 흐름을 만들어 내는 타입을 말한다. 즉, 유체를 어떻게 밀어내서 유동을 유도하는 것이냐에 따라서 Impeller Type을 결정할 수 있다.

#OpenFoam #튜토리얼 #코드 #수정하기 #방정식 #추가 #식 #화공공대생 안녕하세요. 화공공대생입니다. 이번시간은 식을 추가하여서 원하는 방정식을 만족하도록 수정하는 방법을 알아보겠습니다. 목표 : 비압축성 유동에 에너지 전달 방정식을 추가 지배방정식 코딩의 시작은 지배방정식 부터 시작하죠? 비압축성 + 에너지 방정식 먼저 상위 식 2개는 icoFoam에 구현되어 있습니다. 먼저 우분투를 실행을 해보도록하죠~ 저는 우분투 대신 cmd에서 bash를 통해 진입했어요! 폴더복사 복사해야할 icoFoam 폴더 위치는 아래에 있습니다. $FOAM_SOLVERS/incompressible/icoFoam 해당 파일들을 모두 복사해서 원하는 폴더 안에 복사해주죠. 저같은 경우는 AddTeq 폴더를 만들어서 해당 파일들을 전부다 복사해왔네요. mkdir AddTeq cd AddTeq cp -r $FOAM_SOLVERS/incompressible/icoFoam . .은 현재 디렉토리를 말하는

#오픈폼 #OpenFoam #튜토리얼 #화공공대생 #Tutorial 안녕하세요. 화공공대생입니다. 만약 파일구조가 익숙하지 않는다면 여기 링크를 확인해주세요. [OpenFoam] Window Openfoam 튜토리얼 3 파일구조 익히기 #오픈폼 #OpenFoam #CFD #Window #윈도우 #튜토리얼 #Tutorial 안녕하세요. 화공공대생입니다. 오늘... blog.naver.com 오늘은 파일 중 ControlDict에 대해서 알아보겠습니다. 해당 파일은 지배방정식을 결정하는 주요한 파일이며, 해석 시간, 해석 간격과 같은 부분을 설정하는 파일입니다. 여기에서 주요 옵션은 applicaiton이라고 생각하는데, 이 부분에 지배방정식을 포함하고 있습니다. 이에 따라서 어떤 파일들을 읽을지도 결정이되죠. 위 그림을 보자면 경우 inCompressible 솔버 중 icoFoam을 선택하였네요. icoFoam.C 파일을 열어보니 fvVectorMatrix를 Ueqn에 대한 Matri

#CFD #incompressible #Pressure #Solver #PISO #Algorithm 이전 시간에는 간단하게 Projection Method와 SIMPLE Algorithm에 대해 알아보았습니다. [CFD] Projection Method (투영법), Decoupling! 이론편 #CFD #Projection #Method #Chorin #영사 #유동해석 안녕하세요. 화공 공대생입니다. Projection M... blog.naver.com PISO algorithm (Pressure - Implicit with Splitting of Opertators) 은 SIMPLE algorithm의 확장된 버전에 해당합니다. 당연히 동일하게 비압축성 유체Solver로 사용하고, 비정상상태 솔버에 사용합니다. SIMPLE과 똑같이 Projection Method를 이용합니다. SIMPLE하고 거의 똑같지만 압력 보정을 여러번 거치는 것이 PISO 알고리즘입니다. SIMPLE보다

#오픈폼 #OpenFoam #CFD #Window #윈도우 #튜토리얼 #Tutorial 안녕하세요. 화공공대생입니다. 오늘은 OpenFoam을 이용하여 CFD Simulation을 하기 위한 Input File 구조들에 대해서 알아보겠습니다. 먼저 파일 파일의 일반적인 구조를 익히기 위해서 OpenFoam 프로그래밍에 대해서 간단하게 이해를 하고 있어야 합니다. OpenFoam은 Runtime에 결정될 수 있도록 프로그래밍이 짜여져있습니다. 이를 위해서 C++을 기반으로 객체지향(OOP)* 형태의 프로그래밍으로 구현이 되어있습니다. 여기서 객체 생성을 위해서 다양한 Input File들을 넣어주게 되는거죠.. >> "간단히 생각하면, Input File로 (객체) 프로그래밍한다. " 라고 생각하면 됩니다. 1. OpenFoam Case 파일의 일반적구조 뭔가 복잡해보이죠?? 하나씩 만들다가는 골머리 아플거 같습니다... 하지만 OpenFoam는 다양한 예제들을 제공해주고 있습니다.

#그리핏 #테니스엘보 #보호대 #엘보보호대 #테린이 #화공공대생 #외측상과염 안녕하세요 화공공대생입니다. 오랜만 포스팅인데 그리핏 엘보보호대 상품 후기로 왔습니다. ️️️ 한줄평을 먼저 해보자면 편하고 흘러내림 없고 제기능을 하는 제품입니다. ️️️ 내돈주고 사라고한다면? 2만원돈 OK! 요새 테니스에 미친듯이 빠져서 포스팅을 미루고 올리지 못하고 있네요.. 그만큼 저는 다양한 부상들과 함께 하고 있는 테린이 입니다. 오늘은 운좋게도 그리핏에서 어찌 알았는지 엘보 보호대를 협찬을 해주셔서 사용 후기를 알려드리겠습니다. 테니스엘보우는 외측상과염이라고 팔꿈치 바깥쪽에 생기는 염증입니다. 테니스 엘보 증상 [1] 발생원인은 힘줄이 다치면서 손상부위에 염증이 생기는데 테린이들의 이 부상을 가지는 주 원인은 볼을 테두리에 맞추면서 라켓에서 전달되는 진동이라고 생각하는데요.. 테니스 엘보 방지[1] 이 부상을 방지하기 위해서 많은 테니스인들이 보호대 착용을 하고있습니다. 그리핏 엘보 보호대는

5월의 다낭은 건기로 아름다운 해변과 멋진 산을 구경하기 좋다. 다만, 식당과 여행사들의 장사치들이 눈탱이 치는 도시이다. 다낭을 준비하면서 다들 다낭 XX카페들이 많다. 자유여행을 생각한다면 여기서는 구글 지도에서 필요한 정보만 긁어오자. 맛집, 렌트, 투어는 정보는 믿을만하지 못하다. (눈탱이가 절반이다.) 웬만하면 구글 Map 리뷰를 생각하자. 환전, 일정, 호텔 순으로 정리하고자 한다. 현지 돈 준비하기 (환전) 카드와 현금이 모두 필요한 곳이다. 다낭의 대부분은 현금으로 받기 때문에 현금을 준비하는 것을 추천한다. 카드의 경우 롯데마트, grab에서 유용하게 사용할 수 있으니 미리 신청하자. (신청일 기준 5일 전) 환전 가성비 순으로 보면 다음과 같다. (2024.05.11 기준). 달러 -> 동 (금은방 기준) 100 달러 기준 (256만 동) 카드 (트레블 솔, 로그, 월렛). (솔>로그>월렛) 순이다. 나는 시간이 없어서 월렛을 가져감. 5만 원 기준 (92만 ~천동