Input:

results from solver

output:

circumferential averge physical quantities( such as streamwise velocity)

How to?

1. draw an structured 2D axial meshin Pointwise

a. similar nodes distribution as the mesh in the solver ( normally unstructured)

b. export as plot3D format

2. interpolate solution results on the 2D plane in Tecplot

a. import the CFD results in Tecplot

b. import the 2D mesh

file /load data/

c. interpolate data on the 2D mesh

d. rotate the 2D mesh and then interpolate again ( using for loop)

tecplot macro

#!MC 1410

$!Varset |NumLoop|=360

$!Loop |NumLoop|

$!Varset |num|=(|Loop|*1+0)

$!RotateData

  ZoneList =  [4]

  Angle = |num|

  XVar = 1

  YVar = 2

  ZVar = 3

  NormalX = 1

  NormalY = 0

  NormalZ = 0

$!LinearInterpolate

  SourceZones =  [1]

  DestinationZone = 4

  VarList =  [4-9]

  LinearInterPConst = 0

  LinearInterpMode = DontChange

$!WriteDataSet  "C:\Users\kaiming\Documents\ZJU\output_|num|.dat"

  IncludeText = No

  IncludeGeom = No

  IncludeDataShareLinkage = Yes

  ZoneList =  [4]

  Binary = No

  UsePointFormat = Yes

  Precision = 9

  TecplotVersionToWrite = TecplotCurrent

$!EndLoop

output like this

average the result , reshape the data into an array, then get the circumferential average results with 2D matrix distribution -- matlab

matlab code is as follows:

n=359;

a=[];

b=[];

c=[];

for i=1:n

    filename=sprintf('output_%d.dat',i);

    fileinfo = importdata(filename);

    ux=fileinfo.data(:,7);

    uy=fileinfo.data(:,8);

    uz=fileinfo.data(:,9);

    a=[a,ux];

    b=[b,uy];

    c=[c,uz];

end

aver_Ux=sum(a,2)/n;

aver_Uy=sum(b,2)/n;

aver_Uz=sum(c,2)/n;
aver_Ux=aver_Ux.'  % transpose

aver_Uy=aver_Uy.'  % transpose

aver_Uz=aver_Uz.'  % transpose
aver_Ux_array=reshape(aver_Ux, 499,[]);

aver_Uy_array=reshape(aver_Uy, 499,[]);

aver_Uz_array=reshape(aver_Uz, 499, []);

save -ascii aver_ux_array.dat aver_Ux_array;

save -ascii aver_uy_array.dat aver_Uy_array;

save -ascii aver_uz_array.dat aver_Uz_array;

巴特西

circumferential averge streamwise velocity using Tecplot and Matlab

Input:

output:

How to?

最新文章

热门文章