Question related to fill in boundary conditions in comsol

Click For Summary
SUMMARY

This discussion focuses on solving partial differential equations (PDEs) related to piezoelectric substrates using COMSOL Multiphysics. The user is specifically working with the PDE coefficient form to model the behavior of the substrate under an electric field, defined by equations such as -\rho\omega2\vec{u}- \nablaT=0 and \nabla \vec{D} = 0. The user seeks guidance on specifying boundary conditions, particularly Neumann and mixed boundary conditions, and how to determine coefficients q, g, h, and r for their specific case. The conversation also suggests considering the MEMS/piezoelectric application modes in COMSOL for a more streamlined approach.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with COMSOL Multiphysics, specifically the PDE coefficient form
  • Knowledge of piezoelectric materials and their behavior under electric fields
  • Experience with boundary condition definitions in numerical simulations
NEXT STEPS
  • Explore the MEMS/piezoelectric application modes in COMSOL for simplified modeling
  • Research how to define Neumann and mixed boundary conditions in COMSOL
  • Learn about the implications of different boundary condition coefficients (q, g, h, r) in simulations
  • Study the mathematical formulation of piezoelectric equations and their implementation in COMSOL
USEFUL FOR

This discussion is beneficial for engineers, researchers, and students working with COMSOL Multiphysics, particularly those involved in simulating piezoelectric materials and solving PDEs in physics-related applications.

overgift
Messages
7
Reaction score
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?
 
Physics news on Phys.org
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.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 0 ·
Replies
0
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
7
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 0 ·
Replies
0
Views
1K