Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

FEM: periodic boundary conditions (1D)

  1. Jun 12, 2015 #1
    I am trying to set up the mass matrix for a 1D system which I want to solve using finite elements. So the mass matrix is defined as

    [tex]
    M = \int{NN^T}dL,
    [/tex]
    where N is the finite element linear basis functions. I use hat functions.

    Say I have 10 elements, corresponding to 11 nodes running from -5 to 5 so the spacing is 1. Node 1 is equal to node 11 since I want to employ periodic boundary conditions.

    My issue is that I am not sure how to construct the mass matrix for the 10th node. As shown here, the elements for the 10th node will be (I use periodic boundary conditions, so [itex]x_{N+1}=x_1[/itex])

    [tex]
    M_{10,10} = \frac{x_{1}-x_{10}}{3} = -10/3\\
    M_{10,1} = \frac{x_{1}-x_{10}}{6} = -10/6
    [/tex]
    All other elements have positive values given by 1/3 and 1/6, respectively.

    Are my values for [itex]M_{10,10}[/itex] and [itex]M_{10,1}[/itex] correct? I find it odd that their values are so much different than the values in the "bulk".
     
  2. jcsd
  3. Jun 12, 2015 #2

    RUber

    User Avatar
    Homework Helper

    This is 1D, so is it a line?
    Is x_11 actually x_1 or is it x_1 + P where P is your period in X?
    Thus f(x_11) = f(x_1) but x_11 is not actually x_1.
    If this is the case, use x_11 for x_11 and in the matrix assign it the position of x_1 since it will be multiplied by that node.
     
  4. Jun 15, 2015 #3
    Thanks, that is also what I thought. But then my matrix has dimensions 10x10, but my field will have 11 values since x_11 and x_1 are not the same. But that wont work when I multiply them together (?)

    Am I missing something?
     
  5. Jun 15, 2015 #4

    RUber

    User Avatar
    Homework Helper

    x_11 and x_1 are not the same in reality, but in the matrix they have the same role.
    The main thing is that you are keeping with the periodic nature in your x, if all other entries are 1/3 and 1/6, then so should the ones at the edge.
    You aren't actually putting the 11th value in, but when you plot it, you can manually put the 11th point in.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook