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

Question related to fill in boundary conditions in comsol

  1. Oct 1, 2008 #1
    Firstly I really feel so lucky to find this forum. Since I don't have a strong physics background but now dealing with many problems directly related to physics.

    I'm now doing some simulation in comsol and need to solve some PDEs. I'm using this PDE coefficient form in comsol. The equations need to be solved are:

    -[tex]\rho[/tex][tex]\omega[/tex]2[tex]\vec{u}[/tex]- [tex]\nabla[/tex]T=0
    [tex]\nabla \vec{D}[/tex] = 0
    [tex]\vec{E}[/tex] = -[tex]\nabla V[/tex]

    The equations specify the behavior of a piezoelectric substrate when subjected to an electric field. In these equations the variable need to be solved is [tex]\vec{u}[/tex]=[tex]\vec{u}[/tex](u,v,w,V).

    Then PDE coefficient form in comsol is: -[tex]\nabla[/tex][tex]\cdot[/tex](C[tex]\nabla[/tex][tex]\vec{u}[/tex]+[tex]\alpha[/tex][tex]\vec{u}[/tex]-[tex]\gamma[/tex])+a[tex]\vec{u}[/tex]+[tex]\beta[/tex][tex]\cdot[/tex][tex]\nabla[/tex][tex]\vec{u}[/tex]=[tex]\vec{f}[/tex]

    in this step I need to transfer my equations to this PDE coefficient form.

    Then I need to specify the boundary condition:

    The neumann boundary condition in comsol coefficient form specifies:

    The mixed boundary condition in comsol coefficient form specifies:

    q, g,h,r are coefficients I need to fill in according to my specicial case.

    one of the boundary conditions writes V=Vp, n[tex]\cdot[/tex]=0. And I really not sure how to define this in the mixed boundary condition and decide the value of q,g,h,r. Could anyone with comsol experience give me some hint?
  2. jcsd
  3. Oct 3, 2008 #2


    User Avatar
    Science Advisor
    Gold Member

    Hi overgift,

    I was first wondering is there some specific reason you're using the PDE form rather than the MEMS/piezoelectric application modes? Those would build you the coupling "automatically" and specifying for example the electric potential boundary condition would be a fairly straightforward task? Or perhaps you're doing something which is beyond the capabilities of those implementations.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook