Lumped mass on triangular unstructured mesh

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

 August 4, 2011, 23:08 Lumped mass on triangular unstructured mesh #1 New Member   Felipe Hernandez Join Date: Apr 2011 Posts: 13 Rep Power: 6 I was wondering what method is used for lumped mass matrices on a triangular mesh. Guass-Lobatto integration rules fail (in my experience) because the jacobian of the transformation from the bi-unit square to the triangle goes to zero on one corner. Therefore, zeros sometimes appear on the diagonals. (I'm not sure if I was clear in that explanation!) I've been using a row-sum technique because I haven't had much luck with the "special lumping technique". Is that what is usually done? Thanks!

 August 5, 2011, 03:09 #2 Senior Member   cfdnewbie Join Date: Mar 2010 Posts: 551 Rep Power: 11 Petrov, the transformation as you would describe it would indeed give a singularity when two nodes are projected onto one. Usually, what's done to prevent this is to define the transformation in another way: A |\ | \D |__ \ B C So D is chosen as the midpoint of the side. The Jacobi in that case behaves well, no singular point. You might read the book by Karniadakis about spectral methods for more on this! Hope this helps!

 August 5, 2011, 12:49 #3 New Member   Felipe Hernandez Join Date: Apr 2011 Posts: 13 Rep Power: 6 Great! One more question: Is this the correct generalization to a quadratic element? Code: ```D--G--C A | | | H H J F ==> E D | | | G A--E--B B-F--C``` (With J in the middle of the triangle)

 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 aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52 sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11 Remy Main CFD Forum 1 December 22, 2008 05:49 meaton OpenFOAM Running, Solving & CFD 6 April 27, 2008 07:56 Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09

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