I want to simulate the propagation of a sinusoidal plane wave in a rectangular domain using Finite Elements Method. First, the wave should propagate through a fluid medium, then it will encounter a wall of solid material with some thickness, and then continue its propagation in the fluid. I am aware I have to consider some layer of absorbing material to avoid the spurious effects of wave reflection. What I found in the literature to be most suitable for this type of problem is the PML.

I can illustrate my objective with the following image

enter image description here

Ideally, I'd use an open-source Finite Element software. I have been trying to make this work with FEniCS, and I managed to get some good results in a time-domain simulation when there's no PML. However, I couldn't find a way to implement the PML. Many examples use a frequency-domain approach, which I'm not sure if can be used in this type of problem.

After some research I found a Thesis in which the author developed a PML mixed-formulation of time-domain equations suitable for nonlinear problems. The author also gave the weak form of these equations. But as far as I know, FEniCS demands the equations to be in the bilinear-linear form $a(u,v) = L(v)$, or $F(u,v)=0$ - no mixed-formulation. Also, the coding seems to be a bit complicated for a non-expert.

I've been trying to rewrite his equations into this different Weak Form ever since, but I haven't a formal background in FE analysis, so up to now I've been unsuccessful. I can share my progress if you want me to.

The point is, I'm persuaded that there must be a simpler way to get this done. So I'd like to get from you some advice of methodology I should be following, or maybe point me to some meaningful reference that will help me achieve my goal.

Thanks in advance!


The general wave equation is $\partial ^2 u / \partial t ^2 = c^2 \nabla ^2 u + f$.

The weak form is

$a(u,v) = \int_{\Omega} (vu^{n+1} + c^2 \Delta t^2 \nabla v \nabla u^{n+1})d \Omega$

$L_{n+1}(v) = \int_{\Omega}(2u^{n} - u^{n-1} + \Delta t^2 f^{n+1})vd \Omega$

where the exponent $n$ denotes the time step.

I have successfully implemented this into FEniCS by assembling a linear system with $a$, $L$ and the boundary conditions.

The next step is to implement the PML. I've seen lots of work using frequency-domain formulations, but as I mentioned, I'm not sure how to work with them. So I found this Thesis, where the author developed a set of mixed-formulation equations that account for the PML and are in time-domain. Consider the following for the fluid domain:

enter image description here (eq. 5.9)

and its weak form

enter image description here (eq. 5.10a)

enter image description here (eq. 5.10c)

The problem is I was unable to translate those into FEniCS syntax. There are two equations and FEniCS can only solve one of the following:

  • Linear systems $Ax=b$
  • Linear variational problems $a=L$
  • Nonlinear variational problems $F=0$

That means I must, somehow, rewrite the full set of governing equations into FEniCS syntax. I couldn't find a way to this. I tried adding eqs. 5.10a and c together and making them equal 0, in order to use the nonlinear variational solver.

It didn't work though. Not only the math became over complicated, but also the coding got messy. I felt I had to do a PhD in FEA just to know how code that (I'm a grad student). Instead I'm looking for something easier to understand and to code.

  • $\begingroup$ I think you should at least write down your PDE. Maybe it's obvious to those with CFD backgrounds, but it took me a while before i even realized it was a CFD problem :) Also, a link to the thesis you reference would be helpful, if possible. $\endgroup$ – LedHead May 31 '20 at 17:46
  • $\begingroup$ Thanks for your time @LedHead. I added some information and I hope this will clarify my question. $\endgroup$ – Lucas Vieira May 31 '20 at 21:38
  • $\begingroup$ Can you type your equations using MathJax please? $\endgroup$ – nicoguaro Jun 1 '20 at 16:13
  • $\begingroup$ @LedHead, it seems that it is not a CFD problem. It appears to me that the fluid is an "acoustic fluid", that is, a fluid that satisfies the wave equation. $\endgroup$ – nicoguaro Jun 1 '20 at 16:15
  • $\begingroup$ I think that the OP presents two things: (1) Related to modeling/methods for the problem; and (2) related to how to do this in FEniCS. For the second, I suggest that you ask in FEniCS' Q&A. There is a question related to coupled problems: fenicsproject.org/qa/12931/solving-a-coupled-pdes-using-fenics $\endgroup$ – nicoguaro Jun 1 '20 at 16:20

It looks to me like you want to look at the step-62 tutorial program of deal.II, written by Daniel Garcia. It covers pretty much exactly the application you have in mind.

  • $\begingroup$ That looks promising (I'd upvote your answer, but I don't have enough reputation in this forum). I'll take a closer look as soon as I can. Thanks! $\endgroup$ – Lucas Vieira Jun 2 '20 at 0:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.