PDE TOOLBOX: PLOT HEAT FLUX AT 3D NODAL POSITIONS

Dear Matlab community,

I am trying to plot the heat flux data on a 3d geometry. I am using the pde matlab toolbox, specifically the thermalmodel suite.

Here is my code in a nutshell:

thermalmodel = createpde('thermal','transient');

importGeometry(thermalmodel,'sink.stl');

THEN I SET MESH, BCS, ICS AND THE SOLVE. SOLUTION IS RIGHT

thermalsink = solve(thermalmodel,tlist);

PLOT OF TEMPERATURES IS RIGHT BUT WHEN I TRY TO PLOT THE HEAT FLUX:

figure

[qx,qy,qz] = evaluateHeatFlux(thermalsink);

pdeplot3D(thermalmodel,'FlowData',[qx qy qz]);

I GET THIS ERROR:

Error using pdeplot3D (line 243)

Length of flowdata must be 3*number of points

I followed what is written in many worked examples, but it seems not to work. I guess there is some reshaping needed, although this wasn't apparently needed in the worked examples that can be found on the Internet. Can you help me?

Many thanks

Answer 1

[Giulia Ulpiani]

Got the answer by myself.

I didn't consider it was a transient solution and I had to isolate the time step I wished to display. Thanks anyway