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

Estimation of a diffusion coefficient at a face

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 23, 2009, 16:17
Default Estimation of a diffusion coefficient at a face
  #1
New Member
 
Peter
Join Date: Mar 2009
Posts: 15
Rep Power: 8
pifou74 is on a distinguished road
Hello,

I would like to calculate precisely a term such as div(Gf.grad(phi)) in order to create a source term for an equation.
To achieve this, I need to know how fluent evaluate the value of the coefficient Gf at faces, knowing Gc at cell centers around the face.
I have tried to estimate Gf from surrounding cells c0 and c1 like:
Gf=(2Gc0*Gc1/(Gc0+Gc1)) but it doesn't work very well. I have also tried simple linear interpolation with bad results.

If someone has an idea or an information it would help me a lot.

Thanks
pifou74 is offline   Reply With Quote

Old   September 26, 2009, 06:51
Default
  #2
New Member
 
Reza
Join Date: Sep 2009
Location: Iran
Posts: 7
Rep Power: 7
arhn is on a distinguished road
Send a message via Yahoo to arhn
hi
according to Patankar book (p. 45) you should use harmonic mean with respect to distance between center points and interface.
for interface that palace in meddle your formula is correct but in general:
Gf = G0*G1/( G1(1-fe)+G0(fe) )
fe = (distance between G0 and interface)/( distance between G0 and G1 )
arhn is offline   Reply With Quote

Old   September 26, 2009, 07:47
Default
  #3
New Member
 
Peter
Join Date: Mar 2009
Posts: 15
Rep Power: 8
pifou74 is on a distinguished road
Thank you Arhn for your answer. Do you think that, in fluent, it is made like you explain?
Of course I know the very good Patankar's book. You're suggestion is good, nevertheless hard to apply in unstructured mesh, due to the fact that it is quite hard to estimate precisely fe in this case especially when grid is not orthogonal .
Have you any experience about this?
pifou74 is offline   Reply With Quote

Old   September 26, 2009, 08:19
Default
  #4
New Member
 
Reza
Join Date: Sep 2009
Location: Iran
Posts: 7
Rep Power: 7
arhn is on a distinguished road
Send a message via Yahoo to arhn
I write some MATLAB codes that use unstructured meshes like Delaunay or Voronoi and use that formulation. I since don't use FLUENT, sorry.
arhn is offline   Reply With Quote

Old   September 26, 2009, 08:39
Default
  #5
New Member
 
Peter
Join Date: Mar 2009
Posts: 15
Rep Power: 8
pifou74 is on a distinguished road
Ok thank you arhn.

Have a nice day
pifou74 is offline   Reply With Quote

Old   September 26, 2009, 16:56
Default
  #6
New Member
 
Peter
Join Date: Mar 2009
Posts: 15
Rep Power: 8
pifou74 is on a distinguished road
Re Hello Arhn,

I have think about your message. In a Delaunay grid, I think that you can have some triangles which are not orthogonal?
It would help my understanding if you can tell me, in this case how you estimate fe?
Indeed, in this case (non orthogonal grid) the face center, and the cells centers are not on a same line. Moreover I don't know for you but in my case I can only access to the coordinates of the centers of the face and cells. Do you make a kind of approximation?

Thank you for your help

Last edited by pifou74; September 26, 2009 at 17:26.
pifou74 is offline   Reply With Quote

Old   September 26, 2009, 17:40
Default
  #7
New Member
 
Reza
Join Date: Sep 2009
Location: Iran
Posts: 7
Rep Power: 7
arhn is on a distinguished road
Send a message via Yahoo to arhn
Hi
In my code I had coordinate of cell centroids and vertexes of each triangle.
With command "polyxpoly" by using of coordinate of 4 point (2 cell centroids and 2 vertexes) I find the cross point of line between cell centroids and interface line then with respect to the point's coordinates I calculate Fe.
arhn is offline   Reply With Quote

Old   September 27, 2009, 04:08
Default
  #8
New Member
 
Peter
Join Date: Mar 2009
Posts: 15
Rep Power: 8
pifou74 is on a distinguished road
Thank you for all, I will try
pifou74 is offline   Reply With Quote

Old   September 27, 2009, 04:51
Default
  #9
New Member
 
Reza
Join Date: Sep 2009
Location: Iran
Posts: 7
Rep Power: 7
arhn is on a distinguished road
Send a message via Yahoo to arhn
You are welcome.
arhn 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Import netgen mesh to OpenFOAM hsieh Open Source Meshers: Gmsh, Netgen, CGNS, ... 32 September 13, 2011 05:50
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00


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