# 3D color map of F=F(x,y,z) where data is given as a 2D array of the form [x y z F]

9 views (last 30 days)
Saeid on 3 Jan 2023
Commented: Saeid on 5 Jan 2023
The output of a CFD calculation is usually given in the form of a 2D array [x y z F] where F is a function such as pressure or velocity that is calculated for the given points xyz in the 3D space. The result is then given as a color map, as shown below as an example for a relatively simple channel geomtry. How can this be done in MATLAB giventhe [x y z F] output format? ##### 2 CommentsShow 1 older commentHide 1 older comment
Torsten on 3 Jan 2023
Use "slice" for the three visible planes in one plot.

Karim on 3 Jan 2023
Edited: Karim on 3 Jan 2023
Hello , you can create those plots both for both 2D and 3D elements using the patch command. Essentially you need a connecitivty matrix (i.e. the faces of the patch) which you have since these are the faces of the elements used in the mesh. Hense next to the coordindates of the nodes, you also need the connecivity. If you can share a set of sample data (i.e. the mesh, the nodes and the value for plotting) i can adjust the example to meet you data. You can use the paperclip symbol to attach data (of a small example) to you question.
See below for a demonstration of the concept:
% the data here will differ depending on the used FE package, these are
% outputs from ansys i used as a demo in another matlab answer
% this file contains the xyz coordinates of the grid points
% this file contains the deformation values, which we will use to create the 'F' variable
% let's use the magnitude of the deformation as color
F = sqrt( sum(displacements(:,2:end).^2,2));
% this file containts the node conectivity (i.e. the nodes that make up the elements)
% now create the plot
figure
patch('Faces',nodes(:,2:end),'Vertices',coordinates(:,2:end),'FaceVertexCData',F,'FaceColor','interp',...
'EdgeAlpha',0.15) % put the edge alpha to 0 if you went to hide the element edges
view([-15 45])
axis equal
colormap jet
colorbar
title("Deformed") Saeid on 5 Jan 2023
Hi Karim,
this still looks very promising, thank you for the help!
I may still have sent you two files for mesh and results that correspond to slightly different geometries. Let me check this in my files and, before bothering you again, I will also try to run your routine on a simpler geometry (channel with just a round or slit profile and NOT both of them). Also keep in mind that this is a CFD problem where a viscoelastic (LPTT) fluid enters form the round end of this channel and exits from the narrow end. The BCs are no fluid movement at the outside borders + a certain flow rate (or fluid velocity) at the entrance and the solution is given in the form of fluid velocity at every node, which will then be converted to pressure values at each node using the said viscoelastic constitutive equaion. I don't know whether this then leads to the differences you have observed, but thought it would be helpful to mention.
As to the mesh, I have just used a coarse mesh for speed and to reduce the size of the files I am sending here, and that is probably the reason why the number of nodes is so low, but I can assure you that COMSOL can go up with node numbers.
One final questions about the colormap of the solution: is it not possible because of the discrepancy you mentioned above, or is it generally not posible?

### Categories

Find more on Computational Fluid Dynamics (CFD) in Help Center and File Exchange

R2022b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!