[daenerys]

Hi all,

I'm trying to visualize plots and contours in matlab, same type of plots Ansys Fluent can do automatically.

Lets say I want to plot the static pressure contour. I exported a data file from Fluent after solving the field.
I got an ascii file with 4 columns in it:
column 1 - node numbers
column 2 - x coordinate
column 3 - y coordinate
column 4 - static pressure value

I wrote the following matlab code, but it's not working.

clc
close all
clear all

data=importdata('static_pressure');
index=data.data(:,1);
x=data.data(:,2);
y=data.data(:,3);
p=data.data(:,4);
interval=20;
xplot=x(1:interval:end);
yplot=y(1:interval:end);
pplot=p(1:interval:end);
[xi,yi]=meshgrid(xplot,yplot);
zi=griddata(x,y,p,xi,yi);

figure(1)
hcb=colorbar;
colorTitleHandle = get(hcb,'Title');
titleString = '[Pa]';
set(colorTitleHandle ,'String',titleString);
caxis ([0.5,1.8]);
title('Static Pressure Distribution');
xlabel('X coordinate');
ylabel('Y coordinate');
contourf(x,y,p,10000, 'edgecolor', 'none');

I get the following error:
Error using ==> contourf at 65
Z must be size 2x2 or greater.

Error in ==> load at 20
contourf(x,y,p,10000, 'edgecolor', 'none');


Any help would be appreciated. Thx.

[daenerys]

Hi everybody,

I managed to get it to work, so for future reference here is my code:

clc
close alll
clear all

%read data
data=importdata('static_pressure');
x=data.data(:,2);
y=data.data(:,3);
p=data.data(:,4);
%enforce y==0 where y<1e-9 (to avoid NaN at centerline)
y(y<1e-9)=0;
%make rectangular grid
Nx = 1000;
Ny = Nx;
xx = linspace(min(x),max(x),Nx);
yy = linspace(min(y),max(y),Ny);
[xi,yi]=meshgrid(xx,yy);
%interpolate results on rectangular grid
zi=griddata(x,y,p,xi,yi);
%mirror
xi = [ xi(end:-1:2,: ); xi];
yi = [-yi(end:-1:2,: ); yi];
zi = [ zi(end:-1:2,: ); zi];

%plot
figure()
hcb=colorbar;
colorTitleHandle = get(hcb,'Title');
titleString = '[Pa]';
set(colorTitleHandle ,'String',titleString);
title('Pressure Distribution');
xlabel('X coordinate');
ylabel('Y coordinate');
contourf(xi,yi,zi,100, 'edgecolor', 'none');
axis([0 0.0508 -0.015 0.015]);

Enjoy!

[Bodo1993]

Hi, I use the file 'surfaces' to extract 2D surfaces of the phase fraction from OpenFOAM transient results. I postprocess the data using MATLAB to plot contours of the phase fraction and make videos. I am wondering if I have to use the interpolation function 'griddata' to do so. The problem with 'griddata' is that it takes time to do the interpolation, e.g. a simulation with 100 time steps may take more than an hour to generate the video. Any suggestions?

[aadicfd]

Originally Posted by daenerys

Hi everybody,

I managed to get it to work, so for future reference here is my code:

clc
close alll
clear all

%read data
data=importdata('static_pressure');
x=data.data(:,2);
y=data.data(:,3);
p=data.data(:,4);
%enforce y==0 where y<1e-9 (to avoid NaN at centerline)
y(y<1e-9)=0;
%make rectangular grid
Nx = 1000;
Ny = Nx;
xx = linspace(min(x),max(x),Nx);
yy = linspace(min(y),max(y),Ny);
[xi,yi]=meshgrid(xx,yy);
%interpolate results on rectangular grid
zi=griddata(x,y,p,xi,yi);
%mirror
xi = [ xi(end:-1:2,: ); xi];
yi = [-yi(end:-1:2,: ); yi];
zi = [ zi(end:-1:2,: ); zi];

%plot
figure()
hcb=colorbar;
colorTitleHandle = get(hcb,'Title');
titleString = '[Pa]';
set(colorTitleHandle ,'String',titleString);
title('Pressure Distribution');
xlabel('X coordinate');
ylabel('Y coordinate');
contourf(xi,yi,zi,100, 'edgecolor', 'none');
axis([0 0.0508 -0.015 0.015]);

Enjoy!

Hi Daenerys,
I have been following the same method, however my matlab plots come out to be different from my Ansys plot. Did you also face such problem? I cannot figure out the reason for this problem.

Thanks!