PDE modeler plotting method

I made the 2D simulation model with PED tool box.

Here I attached my PDE model.

It simulate the glass vial and liquid product in the glass vial freezing in the lyophilizer chamber.

Lyophilizer consists of total 3 shelfs and there are total 4 vials.

It has the difference in the shape of the vial at the center shelf and the bottom shelf.

This is for verifying how the vial geometry impact the thermal history of the liquid in the vial.

Although I succeed in making 2D modeler, I do not know how to plot the thermal history of the liquid.

For example, I want to plotting the temperature for top layer of the liquid in accordance with the time.

I want to make the plot for all 4 vials and want to compare how the temperature differ from each other.

But I do not know how to plotting it.

Please help

Answer 1

[Ravi Kumar]

1. You can export solution from Solve dropdown menu. You will have the solution for all time-steps as variable 'u' in your workspace. Each column in 'u' is solution corresponding to each time step.
2. Export the mesh data from Mesh dropdown menu. You will have 'p', 'e', and 't' in your workspace.
3. You can now use the pdeplot function and variable in your workspace to plot the desired temperature For example, plot solution at 5000th time-step as:

   pdeplot(p,e,t,'XYData',u(:,5000))

[Pyo Kwang Jeon]

This is very helpful information! Thanks a lot.

But what i wanted is not pdeplotting.

I want to plot the temperature of the point according to the time.

And i don't know what mesh mean the point of the top of the liquid.

Can i find which row of data 'u' is the very point that i want to plotting?

[Pyo Kwang Jeon]

This is exactly what i want to plot

[Ravi Kumar]

You can use the pdeInterpolant function to interpolate solution to arbitrary point within domain as shown below.

Export the solution and mesh data as mentioned in my previous post, hence, having ''p', 't', and 'u' in the workspace, construct the interpolant:

V = pdeInterpolant(p,t,u);

Then evaluate at a desired location, say x = 460 and y = 350, as:

Tinterp  = evaluate(V,460,350)

Then you can plot Tinterp vs time to get the plot you want.

If you are solving a heat transfer problem, consider using the heat transfer vertical workflow. This is a programmatic workflow and you should be able to simplify your setup and postprocessing using the function of ThermalModel. Refer to the example.

Regards,

Ravi