Zunächst weiß ich leider nicht, wie man im Forum mathematische Formeln erzeugt,
deswegen leider hier als Text.
Integral _ (Neumann-Rand) u^2 * ds
Also das Integral von u^2 auf dem Neumannrand ausgewertet.
Dies kann man als L2 Norm auffassen und somit
|| u^2 ||^2 _ Neumannrand
Ist es möglich dies in Python umzusetzen?
Wenn ja, wie ?
Ich habe bisher folgendes gefunden.
dolfin.fem.norms.norm(v, norm_type='L2', mesh=None)
Es ist also möglich beispielsweise norm(u, 'L2', mesh) zu schreiben.
Allerdings möchte ich die Norm ja nur auf dem Rand und nicht auf meinem gesamten mesh.
Code: Alles auswählen
n = 50
mesh = UnitSquareMesh(n,n)
V = FunctionSpace(mesh, "Lagrange, 1")
p = Expression('5.0')
u = Function(V)
u.interpolate(p)
Vielen Dank
MatheMagie