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

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

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2010, 16:46
Default 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: 6
HectorRedal is on a distinguished road
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.
HectorRedal is offline   Reply With Quote

Old   September 13, 2010, 05:25
Default
  #2
Senior Member
 
Rami Ben-Zvi
Join Date: Mar 2009
Posts: 148
Rep Power: 8
Rami is on a distinguished road
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<s<1, -1<t<1).

Good luck,
Rami
Rami is offline   Reply With Quote

Old   September 13, 2010, 15:58
Default
  #3
Member
 
Hector Redal
Join Date: Aug 2010
Location: Madrid, Spain
Posts: 67
Rep Power: 6
HectorRedal is on a distinguished road
Hello Rami,

Thank you for your answer.
I will try to follow your instructions.

Thanks,

Hector.
HectorRedal is offline   Reply With Quote

Old   March 13, 2015, 09:25
Default
  #4
New Member
 
Samir
Join Date: Jun 2014
Posts: 4
Rep Power: 3
samkam786 is on a distinguished road
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 ?
samkam786 is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Comparison: Finite Volume Method vs. Analytic Method m-fry Main CFD Forum 1 April 20, 2010 14:40
Derivation of Momentum Equation in Integral Form Demonwolf Main CFD Forum 2 October 29, 2009 20:53
Calculate velocity inflow on a 2D surface in 3D quarkz Main CFD Forum 4 May 10, 2009 05:54
Solid-Liquid Two Phase Flow Numerical Simulation with the Finite Element Method Qing Hao Main CFD Forum 2 January 9, 1999 00:31
What is C.V. based finite element method C-H Kuo Main CFD Forum 3 November 5, 1998 10:07


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