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

zero flux bc's on FDM 9-noded patch

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 14, 2001, 12:35
Default zero flux bc's on FDM 9-noded patch
  #1
E Siebrits
Guest
 
Posts: n/a
Given equation:

d/dx(X1*du/dx + X2*du/dy) + d/dy(Y2*du/dx + Y1*du/dy) = ...

where X1, Y1, etc are constants.

2nd-order FDM discretization of this equation system on a regular mesh (with rectangular elements) yields a 9-noded molecule in variable 'u' with the usual W, E, S, N, and C (center) coeffs, but also NW, NE, SW, and SE.

My question is how to correctly apply a zero flux bc, especially when acute or obtuse corners are present on the mesh boundaries of an irregular-shaped 2D domain?

Any pointers will be welcome!

Thanks
  Reply With Quote

Old   August 14, 2001, 14:37
Default Re: zero flux bc's on FDM 9-noded patch
  #2
John C. Chien
Guest
 
Posts: n/a
(1). Through the coordinate transformation of the original equations and the boundary conditions.
  Reply With Quote

Old   August 14, 2001, 21:09
Default Re: zero flux bc's on FDM 9-noded patch
  #3
joe
Guest
 
Posts: n/a
your answer is like this ...

A : How can I cook chinese food ?

B : Through your hand.

  Reply With Quote

Old   August 15, 2001, 00:07
Default Re: zero flux bc's on FDM 9-noded patch
  #4
John C. Chien
Guest
 
Posts: n/a
(1). I guess there is the missing part. That is the mesh. (2). By using the coordinate transformation, you will have the governing equation written in the general coordinates, which is yet unknown. (3). But if you then create a mesh which is normal to the boundary, then the zero flux boundary condition is just: t(Xi_at wall, Eta, Zeta)= t(Xi_at wall - 1, Eta, Zeta), if Xi is the coordinate normal to the local wall. (4). If you are serious about the zero flux condition, then the numerical grid generation is a good one to try. (5). By the way, I normally by-pass the cooking step, and buy the Chinese food (cooked) at the local mall's PANDA EXPRESS fast food place. (6). So, zero flux means the normal gradient is zero, that is you need to define the local wall point and a point on the line normal to the wall point. That means a local mesh line must be normal to the wall. When you have this local behavior of the mesh, the zero flux condition is very simple.
  Reply With Quote

Old   August 15, 2001, 00:19
Default Re: zero flux bc's on FDM 9-noded patch
  #5
John C. Chien
Guest
 
Posts: n/a
(1). If on the other hand, you insist on using the regular coordinates (Cartesian or the like), then you can find a point of the mesh near the wall first, then draw a normal to the wall, find the wall point, then assign the value on the wall equal to the near the wall point. (2). After you have these wall point values updated, then use interpolation to compute the values at other wall points if needed.
  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
Using starToFoam clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04
StarToFoam error Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38
CheckMeshbs errors ivanyao OpenFOAM Running, Solving & CFD 2 March 11, 2009 03:34
Import gmsh msh to Foam adorean Open Source Meshers: Gmsh, Netgen, CGNS, ... 24 April 27, 2005 08:19
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


All times are GMT -4. The time now is 17:53.