# How can be calculate a surface integral in the finite element method?

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

 September 9, 2010, 16:46 How can be calculate a surface integral in the finite element method? #1 Member   Hector Redal Join Date: Aug 2010 Location: Madrid, Spain Posts: 67 Rep Power: 7 Hi, I would like to post the following question about how to calculate the integral of a surface in the finite element method. As far as I know, when calculating, for example the stiff matrix, in the finite element formulation, the integral for such calculation is extended over the whole domain, for example a cube. I am developing a software for cfd that uses the finite element formulation, and there are some integrals that should be performed over the surface of the domain (These integrals are surface integrals). For integrals that are applied over the domain, the jacobi is used for differential transformation. My doubt arises when trying to calculate the surface integrals. As far as I know, the jacobi could not be used. Which kind of differential transformation should be used? I would appreciate any kind of help / suggestion. Thanks in advance, Hector.

 September 13, 2010, 05:25 #2 Senior Member   Rami Ben-Zvi Join Date: Mar 2009 Posts: 147 Rep Power: 8 Hello Hector, I assume (for simplicity) you use brick elements of any order. You integrate as usual (see below for details), except the following: Suppose you use the isoparametric coordinates (r,s,t) with the appropriate transformation to the (x,y,z) coordinates, and suppose your surface is defined by, say, r=1. The you simply use the shape functions (and/or their derivatives) as appearing in the integrand with the substitution of r=1, and instead of the Jacobian relating dV (the volume differential) to dr*ds*dt you will had a relation between dA (the area differential) and ds*dt. That's it. Now integrate either exactly (if possible) or numerically (e.g., Gauss quadrature on the square -1

 September 13, 2010, 15:58 #3 Member   Hector Redal Join Date: Aug 2010 Location: Madrid, Spain Posts: 67 Rep Power: 7 Hello Rami, Thank you for your answer. I will try to follow your instructions. Thanks, Hector.

 March 13, 2015, 09:25 #4 New Member   Samir Join Date: Jun 2014 Posts: 4 Rep Power: 3 hey ! I have been trying to formulate the FEM equation to calculate the temperature across the 2D cross section.. but the extended volume is not in a form of cube, but in the form of a quadrant. I can't use axisymmetric elements... Any hints ?

 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 m-fry Main CFD Forum 1 April 20, 2010 14:40 Demonwolf Main CFD Forum 2 October 29, 2009 20:53 quarkz Main CFD Forum 4 May 10, 2009 05:54 Qing Hao Main CFD Forum 2 January 9, 1999 00:31 C-H Kuo Main CFD Forum 3 November 5, 1998 10:07

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