본문 바로가기

프로그램

ANSYS Fluent 연습 7번째 ANSYS Fluent 프로그램을 독학하고 있습니다... 오늘 연습해본 flunet 해석 진행캡처한 내용입니다. (1) Meshing 해석 진행한 모델은 제일 위에 그림과 같은 형태입니다. 자세한 형태는 언급하지 않겠지만, 좌측에서 우측으로 유체가 지나가는 모델입니다. 이때, mesh를 0.2mm로 hard하게 했더니 mesh짜는데만해도 1시간이 넘게 걸렸습니다. 그래도 혹시나 컴퓨터가 좋지 않아서 '응답없음'이 뜰까봐 기다렸다가 찍은 사진입니다. 그리고 해석 시간을 줄이기 위해서 아래 사진 처럼 모델을 1/4로 잘라서 해석을 진행했습니다. 이때, 1/4로 자른 이유는 해석 시간을 단축시키기 위한 것입니다. 총 nodes수는 위 사진에서 57만개, 아래 사진에서 21만개였습니다. skeness는 둘다 평..
【C언어】 평균과 표준편차를 구하는 코드 예시 평균과 표준편차를 구하는 C언어 코드 예시입니다. 평균과 표준편차를 main() 밖의 함수 하나로 계산하는 것입니다. 제가 직접 코딩한 것이고 visual studio 2010에서는 정상적으로 작동했습니다. 미처 캡처본을 찍지 못했는데, 바로 복붙하셔서 ctrl + F5로 확인해보시면 될 겁니다. C언어 공부하고있는 초반이시면 보시기 좋은 예시가 될 것 같습니다. 필요한 부분에 조금이나마 //설명문으로 적어놓았습니다. #define _CRT_SECURE_NO_WARNINGS #include #include //수학관련 함수를 사용합니다. #define SIZE 5 //사용자가 5개의 숫자를 입력하게 합니다. float p_calc_avr_stdev(float num[], int size, float *m..
【MATLAB】 분수 스칼라 연산 문제풀이 스칼라 연산에서 분수 계산문제 풀이입니다. [문제] 아래의 두 식 1)과 2)를 MATLAB으로 계산하시오. 1) 2) [문제풀이] 1) 수학적인 개념을 잘 이해하고 matlab에서 작성해서 계산해야함을 유의합니다. 우선 1)번식은 아래와 같이 작성하고 계산할 수 있습니다. 2) 위 문제보다 조금더 복잡하지만 아래와 같이 작성해서 계산하면 됩니다. 특별한 설명보다 자세히 보실 수 있도록 띄어쓰기로 구분해뒀으니 참고하시길 바랍니다. 이상!
【MATLAB】 화학반응 Arrhenius 방정식 문제풀이 - Chapra의 응용수치해석 3판 문제 2.20 [문제] 화학반응의 온도에 대한 의존성은 다음과 같은 Arrhenius 방정식으로 계산할 수 있다. 여기서 k는 반응속도(s^-1), A는 진동(또는 주파수) 인자, E는 활성에너지(J/mol), R은 기체상수[8.314 J/(mole*K)], 그리고 Ta는 절대온도(K)이다. 혼합물은 E = 1 x 10^5 J/mol 이며 A = 7 x 10^16 이다. MATLAB을 사용하여 온도 253에서 325K까지 반응속도의 값을 구하라. subplot을 사용하여 (a) Ta에 대한 k(녹색 선)와, (b) 1/Ta에 대한 log10k(붉은색 선)의 그래프를 나란히 그려라. semilogy 함수를 적용하여 (b)를 구하라. 두 subplot의 축에 라벨과 제목을..
【MATLAB】 문제 2.18 풀이 - Chapra의 응용수치해석 3판 문제 2.18 [문제] 다음의 MATLAB 명령어를 입력하면 무엇이 화면에 출력되는가? (a) (b) [문제풀이] (a) 순서대로 풀이하면 이렇습니다. 먼저, x에 5를 대입합니다. 그리고 세제곱을 해주어서 결과를 냅니다. 그리고 y = 8-x에 x를 넣어주어 결과를 8 - 5 이므로 y는 3이 됩니다. 여기서 중간에 쓰여진 x^3은 y=8-x 식에서 사용되지 않습니다. 그 이유는 변수가 이용되지 않았다는 것과 두번째줄의 결과는 세미콜론을 붙이지 않으면 알 수 있듯이 ans = 125로 나올 뿐입니다. 위 상태가 세미콜론(;)을 붙이지 않고 계산해본 결과입니다. 위에서 알 수 있듯이 x^3을 하면 ans변수에 125라는 값이 저장되는 것입니다. 만약에 y=8-ans라고..
【MATLAB】 매개변수 방정식, 3차원 원추형 나선 그래프 문제풀이 - Chapra의 응용수치해석 3판 문제 2.17 [문제] 다음의 매개변수 방정식은 원추형 나선(conical helix)을 나타낸다. 범위 t = 0에서 6π까지 증분 Δt = π/64로 x, y 및 z의 값을 계산하라. subplot을 사용하여 위쪽 구획에 (x,y)의 2차원 그래프(붉은색 실선)를, 아래쪽 구획에 (x,y,z)의 3차원 그래프(청록색 실선)를 그려라. 두 그래프 모두 축에 라벨을 표시하라. [문제풀이] 우선, t가 0부터 6π까지 증분을 Δt = π/64으로 아래와 같이 2가지 표현중 하나를 사용하면 됩니다. 1) 콜론 표기법 2) linspace() 저는 이번문제에서 좀더 표시하기 쉬운 1번 콜론 표기법으로 계속 풀어보겠습니다. 다만 2번 방법으로 하려면 마지막에 사이 갯수를 입력..
【MATLAB】 Maclaurim 급수전개로 cosine 문제풀이 - Chapra의 응용수치해석 3판 문제 2.15 [문제] Maclaurim 급수전개로 cosine을 나타내면 다음과 같다. MATLAB을 이용하여 의 항이 포함될 때까지의 급수전개에 대한 그래프(검정색 점선)와 함께 cosine의 그래프(실선)를 그려라. 급수전개를 계산하는데 내장함수 factorial을 사용하라. x축의 범위를 0에서 3π/2까지로 놓아라. [문제풀이] 우선, 쉬운 cosine 함수부터 그려보겠습니다. 먼저, 가장 위에 clear 명령으로 혹시모를 변수중첩을 방지했습니다. 그리고 x값은 linspace(x1,x2)를 이용해서 0부터 3π/2까지 범위로 100개의 인자를 사용했습니다. linspace의 인자 갯수의 defualt는 100개입니다. 그리고 5번째줄에서 y1 = cos(x)..
【MATLAB】 Manning 방정식 개수로 행렬문제 풀이 - chapra의 응용수치해석 3판 문제 2.10 [문제] Manning 방정식은 직사각형 단면의 개수로 내에서 물의 속도를 계산하는데 사용된다. 여기서 U는 속도(m/s), S는 수로의 기울기, n은 거칠기 계수, B는 폭(m), 그리고 H는 깊이(m)이다. 다음의 데이터는 5개의 수로에서 얻은 것이다. 거칠기 계수 n 수로의 기울기 S 수로의 폭 B 수로의 깊이 H 0.035 0.0001 10 2 0.020 0.0002 8 1 0.015 0.0010 20 1.5 0.030 0.0007 24 3 0.022 0.0003 15 2.5 이 값을 행렬에 저장하는데, 각 행은 각 수로를 표현하고 각 열은 각 매개변수를 표현하도록 하라. 매개변수 행렬 내의 값을 기초로 속도의 열벡터를 계산하는 MATLAB 문장을..
【MATLAB】 간단한 수학계산 문제풀이 [ MATLAB 간단한 수학계산 문제풀이 ] 1번. 2번. matlab에서 곱하기 기호는 * 입니다. 3번. 나누기 기호는 슬래쉬(/) 입니다. 역슬래쉬는 반대로 나눠줍니다. 4번. 괄호안에 4+3을 먼저 수행하여 7이되고 더하기보다 곱하기가 먼저이므로 2*7을 해서 14가 되고 마지막으로 3+14를 통해 답이 17이 나옵니다. 5번. 곱하기와 나누기는 우선순위가 같기때문에 좌측부터 수행이 됩니다ㅏ. 6번. 소수점계산도 가능합니다. 7번. 8번. 더하기보다 지수승이 먼저 계산됩니다. 9번. sqrt()는 루트(root)를 붙여주는 기능입니다. 괄호안에 숫자가 루트 내에 들어가게 됩니다. help sqrt하게 되면 위와 같이 'square root'라는 설명을 보실 수 있습니다. 10번. matlab에서는..
【MATLAB】 Hooke의 법칙, 포텐셜에너지 문제풀이 (for, max) [문제] 힘 F(N)이 스프링을 압축하기 위해 작용할 때, 변위 x(m)는 Hooke의 법칙에 의해 다음과 같이 표현할 수 있다. 여기서 k는 스프링 상수(N/m)이다. 스프링에 저장된 포텐셜에너지 U(J)는 다음의 식으로 계산할 수 있다. 5개의 스프링으로 실험을 하여 다음과 같은 데이터를 수집하였다. 힘 F(N) 14 18 8 9 13 변위 x(m) 0.013 0.020 0.009 0.010 0.012 MATLAB을 이용하여 F와 x를 벡터로 저장하라. 그리고 스프링 상수와 포텐셜에너지의 벡터를 계산하라. max 함수를 이용하여 최대 포텐셜에너지를 구하라. [문제풀이] 우선, 위 표를 벡터화해보면 다음과 같습니다. 다음으로 스프링 상수를 구해보겠습니다. 5개의 스프링을 사용했다고 했으니 위 F와 x의..