CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Fluent UDF and Scheme Programming

Discretization of gradient vector - unstructured mesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 14, 2012, 04:28
Default Discretization of gradient vector - unstructured mesh
  #1
New Member
 
Join Date: Jun 2012
Posts: 5
Rep Power: 6
Dapb is on a distinguished road
Hello,

I have an unstructured triangular mesh that has two different zones (on bottom: zone1, on top: zone2). To impose the boundary condition at the interface, I have these equations:

a [grad(T)*n]1 + b = c [grad(T)*n]2 + d
Tf1 = Tf2 = Tf (interface values)

where the indices 1 and 2 designate the zone and a, b, c and d are constant.

I would like to know how can I descretize the scalar product grad(T)*n at the interface (if it was a structured mesh, I could use (Tf-T1)/ds and (T2-Tf)/ds, but this is not the case) to can impose the value of Tf at the interface.

Thanks in advance

Last edited by Dapb; June 15, 2012 at 03:55.
Dapb is offline   Reply With Quote

Old   June 15, 2012, 03:54
Default
  #2
New Member
 
Join Date: Jun 2012
Posts: 5
Rep Power: 6
Dapb is on a distinguished road
Another question:

Instead of imposing the value of Tf, is there a way to impose directly the flux at both interfaces (interface and interface_SHADOW) even in this situation:

flux(interface) = a [grad(T)*n]1 = c [grad(T)*n]2 + d - b

flux(interface_SHADOW) = c [grad(T)*n]2 = a [grad(T)*n]1 + b - d

I tried to do it using C_T_G(c,t) in the function DEFINE_PROFILE(name,t,i) to calculate grad(T) in the zones 1 and 2, but Fluent pointed out an error.
Dapb 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
[ICEM] Negative volume error in hybrid mesh siw ANSYS Meshing & Geometry 4 September 3, 2014 05:25
[ICEM] Problem making structured mesh on a surface froztbear ANSYS Meshing & Geometry 4 November 10, 2011 09:52
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10


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