A Structure-Preserving Partitioned Finite Element Method For The 2D Wave Equation
A Structure-Preserving Partitioned Finite Element Method For The 2D Wave Equation
A Structure-Preserving Partitioned Finite Element Method For The 2D Wave Equation
4
Hamiltonian Methods for Nonlinear Control, Universidad Técnica
Federico Santa María, Valparaíso, Chile, May 1-4, 2018
Abstract: Discretizing open systems of conservation laws while preserving the power-balance
at the discrete level can be achieved using a new Partitioned Finite Element Method (PFEM),
where an integration by parts is performed only on a subset of the variables in the weak
formulation. Moreover, since boundary control and observation appear naturally in this
formulation, the method is suitable both for simulation and control of infinite-dimensional port-
Hamiltonian systems. The method can be applied using FEM software, and comes along with
worked-out test cases on the 2D wave equation in different geometries and coordinate systems.
1 1
0.8 0.8
0 0
0 0
1.19 s 1.59 s
1 1
0.8 0.8
0
0
0
0
As we did with our own codes in § 4.1, a convergence
analysis of the first mode was also performed for the square
domain using FreeFem++. The results for the error of the
Fig. 2. Snapshots of simulation for a harmonic excitation first natural frequency are presented in Fig. 5 and follow
at two of the boundaries of the domain. The variable the same trend that we obtained with our code.
α1 is shown.
2
10
P1P0
Finally, the convergence of the time domain simulation P1P0P0
P1P0P0
P1P1
P1P1P1
P1P2P2
P1P2
with boundary conditions given by (32) was verified. Fig. 1
P1P1P1
P1P2P2
10
3 shows the L2 norm of the error, taking into account the
state of the system after one second of simulation. The
numerical result obtained using 1600 elements was taken 0
10
error (%)
10 -2 −2
10
−3
10
1
0 1 2 3 4
10 10 10 10 10
Number of states
L 2 norm error of variable
10 -3
Fig. 5. Convergence of the first natural frequency of
the 2D linear wave equation, discretization using
FreeFem++.
10 -4
5. PERSPECTIVES AND CONCLUSIONS
10 0 10 1 10 2
Two extensions of the PFEM are proposed before conclud-
Number of DOF ing.
Fig. 3. Convergence of the L2 error as a function of the 5.1 Extension to other coordinate systems
number of degrees of freedom of the discretization
scheme. The error is computed after 1 second of
The method presented above was written in cartesian
simulation. Polynomials of order 1, 0 and 0 were
coordinates. For other geometries, it can be easier to deal
used as approximation functions for the variables
with other curvilinear systems, like the polar coordinates.
(P1 P0 P0 ). A first-order convergence is observed.
In this latter case, the methodology to obtain the dis-
cretized equations remain the same: for example, in (26),
4.3 Using FEM software for more complex geometries the same kind of algebraic structure is obtain for J, with
matrices Dρ and Dθ and their transpose. For the choice
One of the main advantages of the method proposed in of the basis functions, polynomials are suitable for the ρ
this paper compared to previous work on power-preserving variable, but care must be taken to use periodic functions
discretization methods is that it proves compatible with for the θ variable: trigonometric polynomials will solve this
classic FEM software. problem, as detailed in (Boyd, 2001, chap. 18).