CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

FEM heat transfer problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 14, 2012, 16:29
Default FEM heat transfer problem
  #1
New Member
 
Rafal
Join Date: Apr 2012
Posts: 1
Rep Power: 0
Reload is on a distinguished road
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.
Reload is offline   Reply With Quote

Old   April 16, 2012, 06:45
Default
  #2
Senior Member
 
Rami Ben-Zvi
Join Date: Mar 2009
Posts: 155
Rep Power: 17
Rami is on a distinguished road
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:
\rho C_{p}\partial T/\partial t=-q_{i,i}+G
with
q_{i}=-kT_{,i}

Assume T distribution within an element is
T=\sum_{I}N{_I}T{_I}
with N_I and T_I 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 N_J over the element volume.

The LHS will be
\partial T_{J}/\partial t\int_{V}\rho C_{p}N_{I}N_{J}dV

The RHS, after using intergation by parts and the Gauss theorem, yields
-\int_{\partial V}N_{J}q_{n}dA-T_{J}\int_{V}N_{I,i}kN_{J,i}dV+\int_{V}N_{J}GdV
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
M_{IJ}\dot{T}_{J}=BC \: term -K_{IJ}T_{J}+F_I

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 05:35.
Rami is offline   Reply With Quote

Reply

Tags
fem, heat flow

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Heat transfer problem at nucleate boiling MrStuebb FLUENT 12 March 14, 2017 14:31
Natural Convection heat transfer problem srinivasa FLUENT 21 November 11, 2016 07:08
simple heat transfer problem openfoam1 OpenFOAM 2 January 27, 2010 15:17
Wall heat transfer coefficient (HTC) problem Mohamed khamis CFX 1 January 16, 2010 00:12
How can I increase Heat Transfer at Domain Interf? B.Simon CFX 3 October 28, 2008 19:53


All times are GMT -4. The time now is 06:09.