zero flux bc's on FDM 9-noded patch

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 14, 2001, 12:35 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

 August 14, 2001, 14:37 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.

 August 14, 2001, 21:09 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.

 August 15, 2001, 00:07 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.

 August 15, 2001, 00:19 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04 Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38 ivanyao OpenFOAM Running, Solving & CFD 2 March 11, 2009 03:34 adorean Open Source Meshers: Gmsh, Netgen, CGNS, ... 24 April 27, 2005 08:19 Andrea CFX 2 October 11, 2004 05:12

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

 Contact Us - CFD Online - Privacy Statement - Top