PDE modeler plotting method

5 次查看(过去 30 天)
Pyo Kwang Jeon
Pyo Kwang Jeon 2018-11-27
评论: Ravi Kumar 2018-11-28
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

回答(1 个)

Ravi Kumar
Ravi Kumar 2018-11-27
  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))
junk2.png
  3 个评论
Pyo Kwang Jeon
Pyo Kwang Jeon 2018-11-28
wanted.PNGThis is exactly what i want to plot
Ravi Kumar
Ravi Kumar 2018-11-28
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

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Geometry and Mesh 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by