CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > ANSYS > CFX

implementing boundary conditions

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 5, 2000, 08:52
Default implementing boundary conditions
Andreas Abdon
Posts: n/a
Hello everybody!

I'm implementing boundary conditions on a user scalar equation in user fortran subroutine USRSRC. I want to assign the scalar a value in the near wall cells. The manual recommends the following procedure:

SP=-c SU=c*value AM=0 for all faces

where c is "set equal to approximately the maximum flux in the CV.. .. to stop potential errors due to the Rhie-Chow terms." (quote). I have tried this approach and also different constant values of c and found that it may distort the solution seriously, especially when the flux (velocity) is close to zero. Does anybody have a recommendation how to proceed? Is a large (1.0E10) value of c "safe" to use?


Andreas Abdon

Div. Heat Transfer Lund Institute of Technology Box 118 22100 Lund Sweden
  Reply With Quote

Old   March 14, 2000, 13:00
Default Re: implementing boundary conditions
Andreas Abdon
Posts: n/a
I'm stuck!

** Does ANYBODY have any info on this matter, please help! **

Even if you haven't come across these problems when implementing boundary conditions I would be grateful if you let me know how you did it. I have traced the problem to one of my equations, it's a diffusion eq. with NO CONVECTION. I have tried different c values and the high ones (like 1e10) gives a reasonable solution but slows down convergence seriously. A low value gives a very good convergence rate (final residual levels as for c=1e10) but the solution is not very good.

Andreas (...nervous breakdown coming soon)
  Reply With Quote

Old   March 15, 2000, 13:20
Default Re: implementing boundary conditions
Andreas Abdon
Posts: n/a
----> I FOUND IT! <------

Don't use AM(NCELL,0,IPHASE) when putting AM=0 for all six facial components! You have to loop over each individual face: AM(NCELL,1,IPHASE)=0, AM(NCELL,2,IPHASE)=0,AM(NCELL,3,IPHASE)=0, and so on....

Doing this, the solution is the same for c=1.0 nd c=1.0e10, i.e. dependence of c value is eliminated.


  Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
mesh file for flow over a circular cylinder Ardalan Main CFD Forum 7 December 15, 2020 13:06
Domain Imbalance HMR CFX 5 October 10, 2016 05:57
Impinging Jet Boundary Conditions Anindya Main CFD Forum 25 February 27, 2016 12:58
Fluent accuracy and boundary conditions Paolo Lampitella FLUENT 0 June 12, 2008 06:25
A problem about setting boundary conditions lyang Main CFD Forum 0 September 19, 1999 18:29

All times are GMT -4. The time now is 01:07.