PDE solver on 2D sphere

7 次查看(过去 30 天)
Mark vanRossum
Mark vanRossum 2018-8-29
Hi,
I would like to solve a spatio-temporal non-linear PDE on a 2D sphere (not a 'shell' with a finite thickness).
I see there is a number of mesh generators (cubed sphere, icosahedral grid) for Matlab, but how to integrate them with a PDE solver?
Is it possible with PDEtoolbox or other Matlab tools? Any examples?
Climate modelers probably do this all the time...
  5 个评论
Mark vanRossum
Mark vanRossum 2018-9-3
That's a good idea Torsten. Will the poles give trouble, and do I need special boundary conditions there?
Torsten
Torsten 2018-9-3
I don't know. My guess is that the periodic boundary condition will be the most difficult part to deal with when using the PDE toolbox.
Best wishes
Torsten.

请先登录,再进行评论。

回答(1 个)

Ravi Kumar
Ravi Kumar 2018-9-4
If you are referring to thin walled hollow sphere, a 3-D volume as your problem domain, then you can create such geometry in PDE Toolbox as:
radius = 1;
thickness = radius/10;
gm = multisphere([radius-thickness, radius], 'Void',[1,0]);
model = createpde;
model.Geometry = gm;
pdegplot(model,'FaceLabels','on','FaceAlpha',0.3,'CellLabels','on')
print('Hollow_Sphere','-dpng','-r200')

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by