|
[Sponsors] |
December 21, 2005, 09:10 |
convective stiffness matrix for FEM
|
#1 |
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, 16:54 |
Re: convective stiffness matrix for FEM
|
#2 |
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, 03:03 |
Re: convective stiffness matrix for FEM
|
#3 |
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, 03:07 |
Re: convective stiffness matrix for FEM
|
#4 |
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, 04:22 |
Re: convective stiffness matrix for FEM
|
#5 |
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, 06:55 |
Re: convective stiffness matrix for FEM
|
#6 |
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. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Force can not converge | colopolo | CFX | 13 | October 4, 2011 23:03 |
Efficient Conversion of Local Stiffness Matrix to Global Stiffness | nlane | Main CFD Forum | 2 | September 14, 2010 19:26 |
Extracted stiffness matrix (K) !? | mechanicaldesign | ANSYS | 0 | January 22, 2010 05:04 |
OpenFOAM version 1.6 details | lakeat | OpenFOAM Running, Solving & CFD | 42 | August 26, 2009 22:47 |
stiffness matrix | happy_clp | Main CFD Forum | 0 | March 25, 2009 04:31 |