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

K-E model

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 4, 2009, 00:08
Default K-E model
  #1
New Member
 
MARK
Join Date: Apr 2009
Posts: 9
Rep Power: 16
shineone is on a distinguished road
hello everyone

i am using K-E model for simulation in downtown area.
could you please tell me how to set wind profile in this software.

Thank you very much for your help
shineone is offline   Reply With Quote

Old   August 11, 2009, 08:20
Default Wind profile
  #2
Member
 
michael malin
Join Date: Mar 2009
Location: london
Posts: 33
Rep Power: 18
phoenics_cfd is on a distinguished road
The wind object is described here:

http://www.cham.co.uk/phoenics/d_pol...m#Wind_Profile

and the wind_profile object here:

http://www.cham.co.uk/phoenics/d_pol...-type.htm#Wind

The wind object is very closely related to the wind_profile object. The latter specifies only the wind inflow boundaries, but later versions of PHOENICS are eqipped with the wind object, which is a single object that always fills the entire domain. Based on the wind direction it creates inflow boundaries at the domain edges using logarithmic or power law profiles on the upwind faces, and fixed pressure boundaries on the downwind faces. In addition, the upper, or sky boundary can also be made into a fixed pressure boundary.
phoenics_cfd is offline   Reply With Quote

Old   August 12, 2009, 05:01
Default
  #3
New Member
 
Song Zhu
Join Date: Mar 2009
Posts: 1
Rep Power: 0
migao is on a distinguished road
nice!PHOENICS is becoming easiler to use!
migao is offline   Reply With Quote

Old   August 12, 2009, 10:02
Default thanks
  #4
New Member
 
MARK
Join Date: Apr 2009
Posts: 9
Rep Power: 16
shineone is on a distinguished road
thank you very much
it is really helpful
i will try follow your instructions
shineone is offline   Reply With Quote

Old   August 14, 2009, 11:00
Default
  #5
New Member
 
MARK
Join Date: Apr 2009
Posts: 9
Rep Power: 16
shineone is on a distinguished road
Thank you, Michael
the dialog you showed to me is exactly what i want to set for wind profile
but I am using phoenics 3.5, it doesn't show the settings as you showed. could you please tell me How can i set these parameter such as roughness height in Phoenics 3.5.

and also i have some problem with relaxation control.
could you please tell me how can i get convergence.
thank you very much.
shineone is offline   Reply With Quote

Old   August 14, 2009, 12:40
Default Wind profile
  #6
Member
 
michael malin
Join Date: Mar 2009
Location: london
Posts: 33
Rep Power: 18
phoenics_cfd is on a distinguished road
There are three alternative methods of introducing a wind inlet boundary condition, as follows: (a) a PIL Do loop over the inlet height in the Q1 file with appropriate PATCH and COVAL statements (see library case T306 as an example ); (b) a single PATCH per inlet plus COVAL statements, with the inlet values supplied via FORTRAN coding in Group 13 Section 12 of GROUND; or (c) a single PATCH per inlet plus COVAL statements, but with INFORM commands to define the wind inlet values, as in the example below, which will need two INLET PATCHes (named BLINL and BLINLW) in the Q1 file plus COVAL statements:

************************************************** **********
Echo InForm settings for Group 13
inform13begin
REAL(VREF,REFHEIGHT,RHIN,ROUGHH,AKAP,ZRFDRH)
REAL(VTAU,VTAUX,VTAUZ)
REAL(VELX,VELZ,IN_ANG,PI)
** constants
AKAP=0.41;PI=3.1415926
** wind speed, reference height, roughness height
VREF=1.5;REFHEIGHT=10.0;ROUGHH=0.022
** wind incidence angle
IN_ANG=30.0*PI/180.

VELZ=VREF*COS(IN_ANG)
VELZ

VELX=VREF*SIN(IN_ANG)
VELX

ZRFDRH=REFHEIGHT/ROUGHH

VTAU=VREF*AKAP/LOG(ZRFDRH)
VTAU

VTAUX=VELX*AKAP/LOG(ZRFDRH)
VTAUX

VTAUZ=VELZ*AKAP/LOG(ZRFDRH)
VTAUZ


RHIN=1.189
CHAR(FORM1,FORMX,FORMZ,FORMK,FORME)

FORM1=VTAU/AKAP
FORMX=VTAUX/AKAP
FORMZ=VTAUZ/AKAP

FORMK=VTAU^2/0.3
FORME=VTAU^3/AKAP

(SOURCE OF P1 AT BLINL IS RHIN*(:FORMZ:*LOGE(YG/:ROUGHH WITH IMA$
T<100)
(SOURCE OF U1 AT BLINL IS RHIN*(:FORMZ:*LOGE(YG/:ROUGHH*:FORMX:*$
LOGE(YG/:ROUGHH: WITH IMAT<100)
(SOURCE OF W1 AT BLINL IS RHIN*(:FORMZ:*LOGE(YG/:ROUGHH*:FORMZ:*$
LOGE(YG/:ROUGHH: WITH IMAT<100)
(SOURCE OF KE AT BLINL IS RHIN*(:FORMZ:*LOGE(YG/:ROUGHH*:FORMK: $
WITH IMAT<100)
(SOURCE OF EP AT BLINL IS RHIN*(:FORMZ:*LOGE(YG/:ROUGHH*:FORME:/$
YG WITH IMAT<100)

(SOURCE OF P1 AT BLINW IS RHIN*(:FORMX:*LOGE(YG/:ROUGHH WITH IMA$
T<100)
(SOURCE OF U1 AT BLINW IS RHIN*(:FORMX:*LOGE(YG/:ROUGHH*:FORMX:*$
LOGE(YG/:ROUGHH WITH IMAT<100)
(SOURCE OF W1 AT BLINW IS RHIN*(:FORMX:*LOGE(YG/:ROUGHH*:FORMZ:*$
LOGE(YG/:ROUGHH WITH IMAT<100)
(SOURCE OF KE AT BLINW IS RHIN*(:FORMX:*LOGE(YG/:ROUGHH*:FORMK: $
WITH IMAT<100)
(SOURCE OF EP AT BLINW IS RHIN*(:FORMX:*LOGE(YG/:ROUGHH*:FORME:/$
YG WITH IMAT<100)

inform13end
phoenics_cfd is offline   Reply With Quote

Old   August 15, 2009, 23:23
Default Thank you
  #7
New Member
 
MARK
Join Date: Apr 2009
Posts: 9
Rep Power: 16
shineone is on a distinguished road
Thanks a lot Michael

I am trying this programming.
and seems it can get better result.
thanks again.
shineone is offline   Reply With Quote

Reply

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
Reynolds Stress model in CFX vs Fluent Tim CFX 1 October 7, 2009 07:19
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 1 January 8, 2001 16:49
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 0 December 21, 2000 04:06
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 1 December 17, 2000 08:01
references about the fan/radiator model Mihai ARGHIR FLUENT 0 December 17, 2000 07:40


All times are GMT -4. The time now is 10:00.