Universit¨at T ¨ubingen T ¨ubingen, den 12.06.2018 Mathematisches Institut
Dr. Bal´azs Kov´acs
7. Exercise sheet for Numerik f ¨ur Differentialgleichungen auf Oberfl¨achen
Exercise 16. An evolving surfaceΓ(t)is given directly by the mappingX(·,t):Γ0 →R3: Γ(t) =
x ∈R3 with x1= p1+max{0,p1},x2 = g(p,t)p2 q
p22+p23
,x3 = g(p,t)p3 q
p22+p23
p∈Γ0
,
where the initial surfaceΓ0is the unit sphere, and the function g(p,t) =e−2t
q
p22+p23+ (1−e−2t)
(1−p21)(p21+0.05) +p21 q
1−p21
.
Write a short code which visualises the surface evolution using the direct mapping. As an initial triangulation use a mesh from a previous programming exercise.
Hint.The code is indeed short. Apart from the information given above, nothing else is needed.
Exercise 17. Visualise the evolving surface Γ(t)from Exercise 12, using the computed velocity, and the pseudo-code from part (c).
Please, download the codehttps://na.uni-tuebingen.de/ex/surfPDE_ss18/func_v.m, and bring your laptop along!
Exercise 18. (a) Derive the fully discrete scheme ak-step BDF method for the numerical solution of the heat equation on an evolving surface. Recall that the corresponding matrix vector formulation is
d dt
M(t)u(t)+A(t)u(t) =b(t), with given initial data.
(b) Write a pseudo-code.
Discussed on the tutorials on 26.06.2018.