Help with finite element program I developed for Diff Eqns

Click For Summary
SUMMARY

The discussion centers on a finite element program (MFEM) developed in Java for solving boundary value problems (BVP) and initial value problems (IVP) related to partial differential equations, specifically the scalar Helmholtz equation in wave propagation through a rectangular waveguide in TE mode. The user encounters issues with applying Dirichlet boundary conditions, resulting in discrepancies between expected and actual field plots. The program utilizes nodal-based rectangular elements and a custom mesh generator called Mmesh. The user seeks assistance in resolving these boundary condition application issues to achieve accurate field representations.

PREREQUISITES
  • Understanding of finite element methods (FEM)
  • Familiarity with boundary value problems (BVP) and initial value problems (IVP)
  • Knowledge of the scalar Helmholtz equation
  • Experience with Java programming and custom mesh generation tools
NEXT STEPS
  • Investigate the implementation of Dirichlet boundary conditions in finite element analysis
  • Explore the effects of boundary conditions on eigenvalue problems in wave propagation
  • Learn about the characteristics of nodal-based rectangular elements in FEM
  • Review techniques for visualizing field plots in electromagnetic simulations
USEFUL FOR

This discussion is beneficial for developers and researchers working on finite element analysis, particularly those focused on electromagnetic wave propagation, boundary value problems, and Java-based simulation tools.

mdn
Messages
49
Reaction score
0
I developed finite element program (MFEM) in java for BVP &IVP to compute partial differential equation.
I am facing one problem and description is as follows
my problem is on generalized eigenvalue problem generated in wave propagation through rectangular wave guide in TE mode. (Differential equation is scalar Helmholtz equation).

I am getting correct answer but without specifying Dirichlet boundary condition! (Tangential component of electric field should be zero on boundary element).
and if i applied boundary condition, getting same answer but field plot is different .
one thing i used nodal based rectangular element and mesh generated in my own flash program Mmesh.

please help me to solve above problem.
 
Engineering news on Phys.org
mdn said:
I developed finite element program (MFEM) in java for BVP &IVP to compute partial differential equation.
I am facing one problem and description is as follows
my problem is on generalized eigenvalue problem generated in wave propagation through rectangular wave guide in TE mode. (Differential equation is scalar Helmholtz equation).

I am getting correct answer but without specifying Dirichlet boundary condition! (Tangential component of electric field should be zero on boundary element).
and if i applied boundary condition, getting same answer but field plot is different .
one thing i used nodal based rectangular element and mesh generated in my own flash program Mmesh.

please help me to solve above problem.

Welcome to the PF.

We probably will need a lot more information to be of help on this question. Can you post some plots that show the output from your program?
 
Dear Berkeman,
I have attached field plots from my program for TE10(without bc),TE20 (without bc),TE11(without bc) and TM11 (with bc).
Thanks in advance.
 

Attachments

  • te10wbc.jpg
    te10wbc.jpg
    25.5 KB · Views: 480
  • te20wbc.jpg
    te20wbc.jpg
    25.9 KB · Views: 454
  • te11wbc.jpg
    te11wbc.jpg
    26.9 KB · Views: 476
  • tm11bc.jpg
    tm11bc.jpg
    27.5 KB · Views: 461

Similar threads

Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
3K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K