FEM heat transfer problem

 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 14, 2012, 15:29 FEM heat transfer problem #1 New Member   Rafal Join Date: Apr 2012 Posts: 1 Rep Power: 0 Hello everyone! I'm trying to make a program that uses finite element method (FEM) to transient heat transfer. I'm using Matlab and already have a meshing algorithm. I'm trying to understand how to go further. The best for me would be to understand a simple example. So I've made one – a simple cube. I've put everything in a *.pdf file to make everything clear by proper text format. I marked all my (5) questions in red. They are also numbered so you can write answers to "question no #". Everything else is written in pdf. Link to download: http://wyslijto.pl/plik/g66y5hcqtt - click in FEM [v0.1].pdf In shortcut: I don't understand solving a FEM heat transient problem. I try to make it in Matlab. I made a simple cube example, it's made from 6 tetrahedra. I'm looking for someone to help me with my problems by solving this example or showing me a text algorithm how to do this - step by step. I will be really grateful for any help. Best regards, Reload.

 April 16, 2012, 05:45 #2 Senior Member   Rami Ben-Zvi Join Date: Mar 2009 Posts: 148 Rep Power: 9 Reload, It's simple, and I'll just sketch the FEM procedure for the case in hand (conduction, transient solution). Let us write the equation (using tensor notation) as: with Assume T distribution within an element is with and being the nodal shape functions and temperatures, and the sum is on the element nodes. Now, apply Galerkin's method, i.e., integrate both hands multiplied by over the element volume. The LHS will be The RHS, after using intergation by parts and the Gauss theorem, yields The first term on the RHS is the boundary normal flux, integrated over the element surface. Now sum over all the elements and assemble to get the global matrix equation You have to specify appropriate initial conditions, the shape functions (apparently you used linear functions) and a time scheme (simple differencing will do). You also have to specify the boundary conditions. You may wish to choose G distribution (e.g., use the same shape function as for T), and - voila! That's the FEM in a nutshell. I hope I have clarified your doubts, Rami Last edited by Rami; April 17, 2012 at 04:35.

 Tags fem, heat flow

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post MrStuebb FLUENT 11 March 27, 2016 11:41 srinivasa FLUENT 20 February 20, 2012 13:49 openfoam1 OpenFOAM 2 January 27, 2010 15:17 Mohamed khamis CFX 1 January 16, 2010 00:12 B.Simon CFX 3 October 28, 2008 19:53

All times are GMT -4. The time now is 01:33.