MATLAB: It’s possible to update nodes value in a mesh that was created with pde toolbox

pdepdetoolbox

Hello,

i'm trying to update my mesh nodes values with the result of a simulation. What i like to do is something like this:

results = solvepde(model);
u = results.NodalSolutions;
mesh.Nodes = mesh.Nodes + u;

i know that this last line is wrong, but exist a function that do this?

Best Answer

Hello Gabriel,

You can achieve this using geometryFromMesh function. Following your example:

results = solvepde(model);
u = results.NodalSolutions;
oldNodes = model.Mesh.Nodes;
newNodes = oldNodes + u';
elements = model.Mesh.Elements;
% Delete the geometry. As you will create new geometry with deformed mesh.
model.Geometry = []
newMesh = geometryFromMesh(model,newNodes,elements);

Regards,

Ravi

Related Solutions

MATLAB: How to set space-varying boundary conditions

You can use applyBoundaryCondition. You just have to create a mesh for the geometry.

% Create a 3-dimensional PDE Model
model = createpde(3);
% Construct the Geometry
[x,y,z] = meshgrid([-3:0.25:3],[-2:0.25:2],[-1:0.25:1]);
% Create the convex hull.
x = x(:);
y = y(:);
z = z(:);
K = convhull(x,y,z);
% Put the data in the correct shape for |geometryFromMesh|.
nodes = [x';y';z'];
elements = K';
% Create a PDE model and import the mesh.
geometryFromMesh(model,nodes,elements);
% Plot geometry
figure
pdegplot(model,'FaceLabels','on')
title('Bracket with Face Labels')
% Now the new stuff
generateMesh(model)
ufun = @(region,state)[region.x.^2;region.x.^2 + region.y.^2;region.x.^2 + region.y.^2 - region.z]
applyBoundaryCondition(model,'dirichlet','Face',6,'u',ufun,'Vectorized','on')
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',[1;2;3]);
results = solvepde(model);
pdeplot3D(model,'ColorMapData',results.NodalSolution(:,1))

Alan Weiss

MATLAB mathematical toolbox documentation

MATLAB: Plotting the Electric Field of a 3D Parallel Plate

To get the gradient over multiple dimensions, get multiple outputs from it. See specifically Contour Plot of Vector Field for a 2D example. A 3D version would require quiver3. (I don’t have the geometry, so I can’t run your code.)

Best Answer

Related Solutions

MATLAB: How to set space-varying boundary conditions

MATLAB: Plotting the Electric Field of a 3D Parallel Plate

Related Question