# convective stiffness matrix for FEM

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Search this Thread Display Modes
 December 21, 2005, 08:10 convective stiffness matrix for FEM #1 Juergen Kertz Guest   Posts: n/a Hello! After some more or less successful experiments I came to the "segregated" methods, which would finally meet my needs considering memory consumption. For previous tries used only diffusive terms (e.g. penalty method), and the stiffness matrix is well done with: Aij=mue*INTEGRAL{dNi/dx*dNj/dx + dNi/dy*dNj/dy}dOmega where i and j are matrix indices. However, as I tried to add the convective term, I found an additional k-index, and I do not know what to do with this. The meaning is not described in any paper I found. Cij=INTEGRAL{wi*(Nk*uk*dNj/dx + Nk*vk*dNj/dy)}dOmega I know that the uk and vk are taken from the previous iteration to avoid the nonlinearity,but I do not know where to put the k-index (if it even means an index) into the algorithm and my code. For further explanation or my set-up: I want to calculate steady viscous flow (plastics extrusion die) in 2D (for the first try) with linear 3-node-triangle-elements. Can somebody please tell me what the k's mean, or point me to a document where the matrix assembly algorithm is described? Thank you! Juergen

 December 21, 2005, 15:54 Re: convective stiffness matrix for FEM #2 Renato. Guest   Posts: n/a Hi Juergen, take a look in page 21 of http://www.nacad.ufrj.br/~rnelias/El...._Tese_Msc.pdf Although written in portuguese, I think it can help you... Cheers Renato.

 December 22, 2005, 02:03 Re: convective stiffness matrix for FEM #3 Rami Guest   Posts: n/a Juergen, According to the Einstein convention a repeated index means summation. Therefore k is a dummy index to be summed on the relevant range (probably k=1,2,3 for 3D). I hope this helps, Rami

 December 22, 2005, 02:07 Re: convective stiffness matrix for FEM #4 Juergen Kertz Guest   Posts: n/a Hi Renato, Thank you very much for that document! This is a very good work! The things are described "for the programmer", and this is what I needed. I think I can handle my problem with this now, I will try today. Did you ever think of translating your work (it is your work, isn't it?) to english (or even better german)? I think A LOT of people could use it very well!! As far as I could read the titles with my little italian knowledge, you describe highly sophiticated concepts in an easy-to-read way. Thanks again! Juergen

 December 22, 2005, 03:22 Re: convective stiffness matrix for FEM #5 Juergen Kertz Guest   Posts: n/a Thank you Rami for you answer! This is exactly what Renato's document says (which I'm trying to read just in this moment) Looking at my problem with this knowledge makes my further work much easier! In fact, I feel the answer is so obvious, I must be some kind of stupid that I didn't get it on my own. Thanks, Juergen

 December 23, 2005, 05:55 Re: convective stiffness matrix for FEM #6 Renato. Guest   Posts: n/a Hi Juergen, When I was stating my CFD studies I had the same kind of doubt about Einstein's notation and after that, I decided to write all FEM matrices to incompressible flow in my MSc's dissertation. Regarding the translation: unfortunately, here in Brazil, we're not encouraged to write our dissertations and thesis in any other language rather than Portuguese. We only publish articles in English. In fact, most of the Brazilians don't need to know a second language. We are a huge country where everybody speak Portuguese, besides being surrounded by Latin countries. It brings me to a curiosity: "How many north Americans learn a second language? Maybe Spanish... They are also a huge country speaking an international language" Lucky guys... Regards and Merry Christmas Renato.

 Thread Tools Search this Thread Search this Thread: Advanced Search 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 Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post colopolo CFX 13 October 4, 2011 22:03 nlane Main CFD Forum 2 September 14, 2010 18:26 mechanicaldesign ANSYS 0 January 22, 2010 04:04 lakeat OpenFOAM Running, Solving & CFD 42 August 26, 2009 21:47 happy_clp Main CFD Forum 0 March 25, 2009 03:31

All times are GMT -4. The time now is 17:24.

 Contact Us - CFD Online - Privacy Statement - Top