Mercurial > fem-fenics-eugenio

annotate devel/example/Ficticious_Domain/Unsteady/L2_penalization/VelocityUpdate.ufl @ 192:945c19831a16

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Unsteady NS with L2 penalization.
author gedeone-octave <marcovass89@hotmail.it>
date Sun, 10 Nov 2013 19:33:44 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
192
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
1 ## Copyright (C) 2013 Marco Vassallo <gedeone-octave@users.sourceforge.net>
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
2 ##
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
3 ## This program is free software; you can redistribute it and/or modify it under
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
4 ## the terms of the GNU General Public License as published by the Free Software
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
5 ## Foundation; either version 3 of the License, or (at your option) any later
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
6 ## version.
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
7 ##
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
8 ## This program is distributed in the hope that it will be useful, but WITHOUT
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
9 ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
10 ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
11 ## details.
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
12 ##
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
13 ## You should have received a copy of the GNU General Public License along with
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
14 ## this program; if not, see <http://www.gnu.org/licenses/>.
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
15
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
16 # Define function spaces (P2-P1)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
17 V = VectorElement("CG", triangle, 2)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
18 Q = FiniteElement("CG", triangle, 1)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
19
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
20 # Define trial and test functions
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
21 u = TrialFunction(V)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
22 v = TestFunction(V)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
23
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
24 # Define coefficients
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
25 k = Constant(triangle)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
26 u1 = Coefficient(V)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
27 p1 = Coefficient(Q)
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
28
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
29 # Define bilinear and linear forms
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
30 a = inner(u, v)*dx
945c19831a16 Unsteady NS with L2 penalization.
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
31 L = inner(u1, v)*dx - k*inner(grad(p1), v)*dx