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

flux/value for UDS

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 30, 2003, 10:12
Default flux/value for UDS
  #1
Andrew Garrard
Guest
 
Posts: n/a
Hello CFD People, I have a question: When creating a UDS one may specify the flux at one boundary and the value at the other and the solver will produce the distribution in between. One can then extract the value of the UDS at the boundary where the flux is specified using a UDF or using the GUI. Is it possiable to do this the other way around? Whereby one sets the value of the UDS on two walls, the distribution is calculated and a UDF extract the value of the flux that would be needed to produce the same value on either of the walls. Any thoughts on this would be apprecited
  Reply With Quote

Old   September 1, 2003, 03:12
Default Re: flux/value for UDS
  #2
Greg Perkins
Guest
 
Posts: n/a
I can't see why you cn't do this so long as your problem is well-posed.

Most times I do use a combination of flux and specified bcs. But there's no reason you can't use just fixed if that makes sense.

Have you tried?? What problem did you have?

Greg
  Reply With Quote

Old   September 1, 2003, 08:21
Default Re: flux/value for UDS
  #3
Andrew Garrard
Guest
 
Posts: n/a
Sorry, I don't really understand what you mean. I have tried using the F_UDSI_G funciton, but then you have to specify the direction vectors and unit vectors and it all gets rather complicated. If there is a UDF command to do this I would be very happy: say, for example, you have a square box with two oppostie edges, say face 1 and face 2. You specity the boundary conditons as UDS value for both. Is there a UDF command that you can use to extract the FLUX at the walls where the value boundary condition has been speicfied.

If you specify the flux at a wall and solve, you can extract the UDS value with F_UDSI. I want to know if you can do it the other way around?
  Reply With Quote

Old   September 1, 2003, 09:31
Default Re: flux/value for UDS
  #4
Greg Perkins
Guest
 
Posts: n/a
I see what you mean - by flux you really mean the gradient f the scalar at the boundary (since the flux may be this multiplied by some other scalar, ie density).

OK, off the top of my head I'd do somethign like (if its a bit worng it cos I just typed in here and not my C compiler...)

give face,f_thread for face and thread for face at boundary...for the flux perpendicular to the boundary: uds is uds you are using...

#define Distance_h(a,b,n) (pow(a[n]-b[n],2.0)) #if ND_ND==2 #define Distance(a,b) (pow(Distance_h(a,b,0)+Distance_h(a,b,1),0.5)) #else #define Distance(a,b) (pow(Distance_h(a,b,0)+Distance_h(a,b,1)+Distance_ h(a,b,2),0.5)) #endif

real A[ND_ND] cell_t cf; Thread *tf; real grad,dn; real xc[ND_ND],xf[ND_ND];

cf = F_C0(face,f_thread); tf = F_THREAD_C0(face,f_thread);

C_CENTROID(xc,cf,tf) F_CENTROID(xf,face,f_thread);

dn = Distance(xc,xf);

grad = (F_UDSI(face,f_thread,uds)-C_UDSI(cf,tf))/dn;

You can work out right the sign.

Greg

  Reply With Quote

Old   September 1, 2003, 09:33
Default Re: flux/value for UDS
  #5
Greg Perkins
Guest
 
Posts: n/a
There's a few typos above but I'll assume you can work those out without any probs.

Greg
  Reply With Quote

Old   September 4, 2003, 10:42
Default Re: flux/value for UDS
  #6
Andrew Garrard
Guest
 
Posts: n/a
Hi, thanks for your help with this, I have been using a C_UDSI fuction so far which has not been producing very helpful resutls, so I am very keen to give your suggestions a try. I am a little busy at the moment but I will get onto it in a few weeks. I have never used a #if and #else function before. Did the forum reproduce the UDF accuratly and the # lines should all be on one line or should these be seperated? I think that this forum has a few bugs when it comes to displaying line breaks. If it would be easier, I would appreciate a .txt file containg the UDF. my actual e-mail address is a.garrard@sheffield.ac.uk Thanks again.
  Reply With Quote

Old   September 4, 2003, 20:32
Default Re: flux/value for UDS
  #7
Greg Perkins
Guest
 
Posts: n/a
Hi

yeah the formatting gets scrambled in the forum, but you can just paste it inot your compiler and reformat.

The #if defines are useful for doing different things and specifying them on the command line at compile time.

In the code I check for either a 2D or 3D version of FLUENT and calculate the distance using either 2 or 3 components of the vector.

I wrote that code in the forum so it might have some typos or small bugs. See how it goes...

Greg
  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



All times are GMT -4. The time now is 13:22.