{"id":284150,"date":"2023-05-16T07:19:53","date_gmt":"2023-05-16T07:19:53","guid":{"rendered":"\/forum\/forums\/topic\/wind-profiles\/"},"modified":"2023-05-16T07:19:53","modified_gmt":"2023-05-16T07:19:53","slug":"wind-profiles","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/wind-profiles\/","title":{"rendered":"Wind profiles"},"content":{"rendered":"

Hi everyone,<\/p>\n

I used the UDF file for the Wind direction and speed in transient mode. But the result looks not true while the upper wind speed is much slower than the building surrounding.  Below is my script for the UDF.<\/p>\n

Looking forward for your suggestions.<\/p>\n

Thanks a lot<\/p>\n

\"\"<\/p>\n

\n
#include “udf.h”<\/div>\n

<\/p>\n

DEFINE_PROFILE(velocity,t,i)<\/div>\n
{<\/div>\n

<\/p>\n

  face_t f;<\/div>\n
  real flow_time;<\/div>\n
  flow_time=RP_Get_Real(“flow-time”);<\/div>\n
  real x[ND_ND];<\/div>\n
  real z;<\/div>\n

<\/p>\n

  if(flow_time >=0 && flow_time < 3600){<\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
    F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1-1)*flow_time\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >=3600 && flow_time < 7200){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1-1)*(flow_time-3600)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 7200 && flow_time < 10800){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1-1)*(flow_time-7200)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 10800 && flow_time < 14400){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1.6-1)*(flow_time-10800)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 14400 && flow_time < 18000){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1.6+(1.9-1)*(flow_time-14400)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 18000 && flow_time < 21600){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1.9+(1-1.9)*(flow_time-18000)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 21600 && flow_time < 25200){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(2.5-1)*(flow_time-21600)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 25200 && flow_time < 28800){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (2.5+(4.1-2.5)*(flow_time-25200)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.14);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 28800 && flow_time < 32400){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (4.1+(3.1-4.1)*(flow_time-28800)\/3600)*pow(270\/10,0.22)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 32400 && flow_time < 36000){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (3.1+(1.8-3.1)*(flow_time-32400)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 36000 && flow_time < 39600){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1.8+(1-1.8)*(flow_time-36000)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 39600 && flow_time < 43200){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1.6-1)*(flow_time-39600)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.14);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 43200 && flow_time < 46800){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1.6+(1-1)*(flow_time-43200)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 46800 && flow_time < 50400){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1-1)*(flow_time-46800)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
  if(flow_time >= 50400 && flow_time < 54000){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = (1+(1-1)*(flow_time-50400)\/3600)*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
   if(flow_time == 54000){    <\/div>\n
    begin_f_loop(f,t)<\/div>\n
    {<\/div>\n
  F_CENTROID(x,f,t);<\/div>\n
    z=x[2];<\/div>\n
    F_PROFILE(f,t,i) = 1*pow(270\/10,0.14)*pow(z\/370,0.22);<\/div>\n
    }<\/div>\n
    end_f_loop(f,t)<\/div>\n
  }<\/div>\n
 <\/div>\n

<\/p>\n

}<\/div>\n<\/div>\n","protected":false},"template":"","class_list":["post-284150","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t