
[Sponsors] 
September 21, 2011, 03:04 
Gauss Quadrature points

#1 
Member
supercommandodhruv
Join Date: Sep 2011
Posts: 57
Rep Power: 6 
Hi All,
I am calculating volume flow in a channel, and want to use the Gauss quadrature to calculate the volume flow rate. I searched on the net, and found some tables that have the points ( abscissa) and the corresponding weighting functions. As far as i understood, I can select any of the points in my domain as a point corresponding to one of the Gauss point (say of 5X5 matrix points), and then can use coordinate transformation to select the other points. What happens however, if the other points are outside the domain? How can I select the points then? Is there some rule to fix the points according to my domain? Thanks, Dhruv 

September 21, 2011, 16:47 

#2 
Senior Member
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 8 
You can set those values to zero. But that might compromise the accuracy of the Gaussian quadrature: it is exact for a polynomial of degree 2n and below. How well is your data represented by a polynomial if some of it is identically zero? Probably not so good.


September 22, 2011, 02:42 
Gaussian Quadrature

#3  
Member
supercommandodhruv
Join Date: Sep 2011
Posts: 57
Rep Power: 6 
Quote:
Thanks for the reply. However, I did not fully understand your explanation. What do you mean when you say some of it is identically zero? I am attaching here the geometry, the cut and the region where I have the computation results. I want to set some points (say a 4 X 4 array) and find out the volume flow. How can I find a polynomial which can represent my data points? 

September 22, 2011, 08:56 

#4  
Senior Member
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 8 
Quote:
As for your problem, I have two suggestions:  use an OF postprocessing utility. I have never done this one, but search the forum for integrating fluxes over patches.  Extract the points, and use a trapezoidal or a Romberg integration. I don't think Gaussian quadratures are of much help here, as they require nodes at specific locations. That will force you to interpolate the flux to those locations, and then you face the question of how do you interpolate. Finally, when I think about it, since this is a finite volume method, the flux is constant across a cell (at least in orthogonal meshes). Thus the selfconsistent way of getting the total flux is to add the fluxes of all the cells on a patch. Mirko 

September 26, 2011, 04:59 
Gau

#5  
Member
supercommandodhruv
Join Date: Sep 2011
Posts: 57
Rep Power: 6 
Quote:


Tags 
gauss quadrature, weighting functions 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Problem in fvschemes divSchemes cannot use Gauss linearUpwind  nico765  OpenFOAM Bugs  27  June 3, 2014 05:29 
please help me which solver is better for my application  Ger_US  OpenFOAM  8  January 8, 2013 13:29 
Submarine with SimpleFoam  alex_rubel  OpenFOAM  20  February 25, 2011 07:34 
solution diverges when linear upwind interpolation scheme is used  subash  OpenFOAM  0  May 29, 2010 01:23 
BlockMeshmergePatchPairs  hjasak  OpenFOAM Native Meshers: blockMesh  11  August 15, 2008 07:36 