- #1
overgift
- 8
- 0
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]
T=cES-eEi
Di=[tex]\epsilon[/tex]S
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:
n[tex]\cdot[/tex](C[tex]\nabla[/tex][tex]\vec{u}[/tex]+[tex]\alpha[/tex][tex]\vec{u}[/tex]-[tex]\gamma[/tex])+q[tex]\vec{u}[/tex]=[tex]\vec{g}[/tex]
The mixed boundary condition in comsol coefficient form specifies:
n[tex]\cdot[/tex](C[tex]\nabla[/tex][tex]\vec{u}[/tex]+[tex]\alpha[/tex][tex]\vec{u}[/tex]-[tex]\gamma[/tex])+q[tex]\vec{u}[/tex]=[tex]\vec{g}[/tex]-hT[tex]\mu[/tex]
h[tex]\vec{u}[/tex]=[tex]\vec{r}[/tex]
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?
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]
T=cES-eEi
Di=[tex]\epsilon[/tex]S
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:
n[tex]\cdot[/tex](C[tex]\nabla[/tex][tex]\vec{u}[/tex]+[tex]\alpha[/tex][tex]\vec{u}[/tex]-[tex]\gamma[/tex])+q[tex]\vec{u}[/tex]=[tex]\vec{g}[/tex]
The mixed boundary condition in comsol coefficient form specifies:
n[tex]\cdot[/tex](C[tex]\nabla[/tex][tex]\vec{u}[/tex]+[tex]\alpha[/tex][tex]\vec{u}[/tex]-[tex]\gamma[/tex])+q[tex]\vec{u}[/tex]=[tex]\vec{g}[/tex]-hT[tex]\mu[/tex]
h[tex]\vec{u}[/tex]=[tex]\vec{r}[/tex]
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?