Middle East Technical University
Mechanical Engineering Deparment
ME 582 Finite Element Analysis in Thermofluids
|Announcements | Files | Links|
13-05-2018Term Porject is available at the Files page.
12-05-2018Handouts 15, 16 and 17, steady Navier-Stokes solver code with sample input files and input file generators are available at the Files page.
02-05-2018Study Set 7 (including Hw 7) and Handout 14 are available at the Files page.
25-04-2018Handout 13 is updated with more examples.
18-04-2018Handout 13 is available at the Files page.
17-04-2018unsteady1D.m MATLAB code, Handouts 11 and 12 and the animation generated in COMSOL Tutorial 3 are available at the Files page.
11-04-2018Handouts 9 and 10 are available at the Files page.
04-04-2018generate2DinputFile.m MATLAB code, Handout 8 and Study Set 5 (including Hw 5) are available at the Files page.
28-03-2018steady2D.m MATLAB code and sample input files, Handouts 5, 6 and 7, and Study Set 4 (including Hw 4) are available at the Files page.
14-03-2018steady1D.m MATLAB code and sample input files, and Handout 4 that explains them are available at the Files page.
Study Set 3 (including Homework 3) is available at the Files page.
13-03-2018In Q4 of Hw 2 there should be "k" in the differential equation as shown below. You need it because "k" is different in each material. Sorry for noticing this mistake too late.
-d/dr(k r dT/dr)= 0
And a small warning for Q1. In the Adv.-Diff. equation we used in Chapter 2, there was a minus sign in front of the diffusion term. With that minus sign the boundary terms of the weak form were like + (w SV)_a + (w SV)_b . But in this question there is no minus sign in front of the second order term. So be careful not to make sign mistakes for NBC and MBC.
11-03-2018In part (c) of Q3 of Hw2, length of each element is given as 0.01 m. It should be 0.015 m.
08-03-2018I got some questions about Q1 of Hw2. Right hand side function f is not constant, but a function of x. When you change your integration variable from x to ksi, you also need to express all the x's in f(x) in terms of ksi. Just replacing x with ksi will not work. You need to use the x-ksi relation, which is different for each element. Therefore you need to calculate separate Fe's for each element. See Example 2.7 of Handout 3, where a DE with a non-constant coefficient is solved. As mentioned in other questions of the homework, you can use computer to calculate the integrals and automate or semi-automate the solution.
Last lecture, I might have made a mistake in writing the quadratic shape functions on the board. They should be as follows. Please check and correct them in case there is a mistake.
S1 = 0.5*ksi*(ksi-1)
S2 = (1-ksi*ksi)
S3 = 0.5*ksi*(ksi+1)
07-03-2018"Handout 3 - Examples of Chapter 2" and the second study set (including the second homework) are available at the Files page.
15-02-2018Spring 2018 Syllabus is available at the Files page.
Here is the list of students who can take the course. Any student whose name is not on the list can attend the lectures as an audience.
05-02-2018ME 582 web site is getting ready for the new semester. Hope to see you in the first lecture.