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

FEM num integration: quadrature rules for a cube

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 11, 2008, 13:53
Default FEM num integration: quadrature rules for a cube
  #1
Anthony
Guest
 
Posts: n/a
Hi all, I am looking for a web site where I can find some high order quadrature rules (i.e. integration formula) for a cube, which are exact up to order 7. I am using Q2 cubic elements in my FEM code.

I don't want to use any 3D generalization of a 1D n-integration points (IP) Gauss formula because the number of IP varies like n.n.n, which is huge. I know some cheaper (less IP) integration formula exist but did not manage to find them on the web.

Seems to be quite tough to find a web site that gather these formula.

Thanks for your help, Anthony
  Reply With Quote

Old   December 11, 2008, 14:51
Default Re: FEM num integration: quadrature rules for a cu
  #2
Jed
Guest
 
Posts: n/a
Encyclopaedia of Cubature Formulas:

http://www.cs.kuleuven.ac.be/~nines/...h/ecf/ecf.html

I don't know what you mean by Q2 cubic elements. Q2 elements usually mean tensor product quadratic which means the maximum order is 6 and you end up not beating the tensor product rule. If you are using a truncated space then the cubature rules above will help.

Note that if the basis and quadrature rule are tensor products, you can exploit the tensor products for evaluation, never forming the <code>Q\times P</code> basis and derivative matrices. This operation is sufficiently fast that I don't even form the Jacobian for Q2 and higher order, I just form a Q1 preconditioning matrix. This uses half the memory and the linear solver (preconditioned with ML's parallel algebraic multigrid) converges in half the time of an assembled Q2 matrix, but I get arbitrary spectral order (Q2-Q8, then the tensor product operation breaks out of L1 cache so further p-refinement is less efficient). If you're going to bother assembling the matrix, by all means use truncated spaces and the cubature formulas, but tensor products are nice too.
  Reply With Quote

Old   December 12, 2008, 06:07
Default Re: FEM num integration: quadrature rules for a cu
  #3
Anthony
Guest
 
Posts: n/a
A Q2 cubic element is a finite element the shape of which is a cube with 27 nodes (8 corners + 12 middles of edges + 6 centers of face + 1 center of the cube = 27). The basis functions of each node are quadratic (2d order polynomial) i.e. Q2.

What I am looking for is a Gauss or Newton-Cotes like integration formula to compute the sum of a polynomial function of order 6 or 7 on the basic cube [-1:1] x [-1:1] x [-1:1] or alternatively [0:1] x [0:1] x [0:1].

With this type of element, the mass matrix terms are polynomes of order 6 (w_i*w_j*det J = 2+2+2 = 6 ), that's why to integrate these terms exactly, I need a quadrature rule or order 6 or 7.

Hope it is clearer right now, sorry to have been cryptic. Anthony

  Reply With Quote

Old   December 13, 2008, 15:48
Default Re: FEM num integration: quadrature rules for a cu
  #4
Jed
Guest
 
Posts: n/a
I thought cubic possibly referred to order 3 polynomials which didn't make sense. I don't think you can beat the tensor product quadrature rule here, just use the 3*3*3 Gauss rule. If you use a truncated space (like corners + edges) then there are better rules (see the link in my last post).

Note that there <code>f = (1-x^2) (1-y^2) (1-z^2)</code> is the center basis function so your mass matrix needs to integrate <code>f*f*det(J)</code> which has order at least 12 (depending on the element mapping type). There is a big difference between getting all polynomials of order 12 (which requires many points) and just tensor products of quartics (which the 27-point rule will do).

As my last post indicated, assembling a matrix based on Q2 or higher elements is usually counter-productive.
  Reply With Quote

Old   December 16, 2008, 08:36
Default Re: FEM num integration: quadrature rules for a cu
  #5
Anthony
Guest
 
Posts: n/a
Thx Jed, Anthony
  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
Quadrature rules on tetrahedron Jonas Holdeman Main CFD Forum 1 July 9, 2011 02:41
Gambit help: Cube inside cube Jack Martinez FLUENT 13 August 11, 2010 06:29
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM kenn Main CFD Forum 2 July 18, 2004 18:28
Integration Rules Khurram Main CFD Forum 1 January 21, 2003 21:38
How good is CFD? kai Main CFD Forum 39 April 7, 2000 12:48


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