Page History

{list-pages:direction=siblings}
h3. A complete example for a lossless cavity
{code}
  ModelInfo: {
    File: dds3.ncdf                             //mesh file. It is the file converted using acdtool
    BoundaryCondition: {                   //specify boundary conditions. The numbers here are sideset in cubit
      Magnetic: 1, 2                           //reference surfaces 1 and 2 are symmetric planes
      Electric: 3 4                               //set reference surfaces 3 and 4 to be electric boundary condition
      Exterior: 6                                 //surface group 6 (maybe many surfaces) is metal
   }
   SurfaceMaterial: {                          //for each metal (exterior) surface group, list the sigma values
      ReferenceNumber: 6
      Sigma: 5.8e7
    }
  }

  FiniteElement: {
     Order: 2                                     //set the finite element basis function order to be used.
     CurvedSurfaces: on
  }

  EigenSolver: {
      NumEigenvalues:     1                //want to compute 1 mode
      FrequencyShift: 		10.e9     //the eigenfrequency of the mode should be above 10GHz
  }

  CheckPoint: {
    Action: save
    Directory: eigens                         //eigenvectors are saved out into this directory
  }

  PostProcess: {
    Toggle: off                                   //postprocess switch
    ModeFile: dds                              //The prefix of the mode filename.
  }

  Log: thisrun.log                              //If you want more printout logged into the file
{code}

with lossy materials
{code}
 ModelInfo: {
  File: ./pillbox.ncdf
  BoundaryCondition: {
    Electric: 1,2,3,4
    Exterior: 6
  }
   Material : {
     Attribute: 1
     Epsilon:   1.0
     Mu:        1.0
   }
   Material : {
      Attribute: 2
      Epsilon:   1.0
      Mu:        1.0
      EpsilonImag: -0.2        //lossy material
   }
 }

  FiniteElement: {
   Order: 1
   Curved Surfaces: off
  }

  PostProcess: {
    Toggle: off
    ModeFile: mode
    SymmetryFactor: 2
  }

  EigenSolver: {
    NumEigenvalues:     2
    FrequencyShift: 		5e9
  }
{code}
h3. A complete example with waveguide loaded cavity

{code}
  

ModelInfo: {
    File: cell1fourthc026ds-pbc.ncdf
    BoundaryCondition: {
      Magnetic: 1, 2,3,4
    Periodic_M: 3 Exterior: 6
 //master surface
    WaveguidePeriodic_S: 7 4   //forslave eachsurface, numberthe appearedmesh here,should itbe shouldexactly havesame atas leastthose oneon Portthe containermaster later.surface
    }
  }

  FiniteElementExterior: {6
    OrderTheta: -150          1//phase
    Curved Surfaces: on
  }
}

  PostProcessFiniteElement: {
  Order: 2
  ToggleCurvedSurfaces: on
    ModeFileScalarPotential: test
1 //use A-V formulation
}

PostProcess:
{
  Toggle: on
  ModeFile: mode
  SymmetryFactor:  8.
}
EigenSolver: {
    NumEigenvalues:     1
    FrequencyShift: 	10e9
}

	9.e9
  }

  CheckPointModelInfo: {
    ActionFile: savecell1fourth.ncdf
    DirectoryBoundaryCondition: eigens{
   }


   PortMagnetic: {1,2,3,4
      Exterior:  ReferenceNumber6
      Waveguide: 7     //thisfor each number appeared here, it should matchhave surfaceat groupsleast inone waveguidePort boundarycontainer conditionlater.
 Absorbing and Waveguide have the same effects. Origin:Omega3P internally will figure  0.0, 0.0415, 0.0out which BC to use.
    }
 //the origin of the 2D port in the 3D coordinate system
 }

  FiniteElement: {
    Order:          XDirection: 1.0, 0.0,
    Curved 0.0 Surfaces: on
  }

 //the xPostProcess: axis{
 of  the 2DToggle: porton
 in the 3D coordinateModeFile: systemtest
  }

  EigenSolver: {
    YDirectionNumEigenvalues: 0.0, 0.0,   -1.0
   //the y axis of the 2D port in the 3D coordinate systemFrequencyShift: 		9.e9
  }

  CheckPoint: {
    Action: save
    Directory: eigens
  }


  ESolverPort: {
        ReferenceNumber: 7     //this number Type:should Analyticmatch surface groups in waveguide boundary condition.
        //analyticOrigin: expression is used
  0.0, 0.0415, 0.0    //the origin of the 2D port in the Mode: {3D coordinate system
        XDirection: 1.0, 0.0,    0.0    //the x axis of the 2D WaveguideType:port Rectangularin the 3D coordinate system
 //it is a rectangular waveguide
   YDirection: 0.0, 0.0,   -1.0   //the y axis of the 2D port in the 3D coordinate system
  ModeType      ESolver: TE{
 1 0              Type: Analytic     //load the TE10 mode
      //analytic expression is used
                AMode: 0.028499{
                        WaveguideType: Rectangular     //dimensionit is ofa therectangular waveguide in x
                        BModeType: 0.0134053TE 1 0                    //load the TE10 mode
  //dimension  of the waveguide in y
                 }
A: 0.028499          }
   }
{code}
h3. Load TEM mode in a coax waveguide
{code}
   Port: {
    //dimension of the waveguide ReferenceNumber:in 2x
         Origin:     0.0, 0.0, 0.011
        ESolverB: {
0.0134053                Type: Analytic
          //dimension of the waveguide in  Mode: {y
                        WaveguideType: Coax}
         }
   }

   Port: {
        ModeTypeReferenceNumber: TEM2
        Origin:        0.0, 0.0, 0.011
        AESolver: 0.0011{
   //smaller radius
            Type: Analytic
           B: 0.0033  //larger radius
 Mode: {
              }
        }
   }

{code}

h3. Load TE11 mode in a circular waveguide
{code}
   Port: {
WaveguideType: Coax
                   ReferenceNumber: 2
    ModeType: TEM
   Origin:     0.0, 0.0, 0.1
        XDirection: 1.0,  0.0,   A: 0.0
0011  //smaller radius
     YDirection: 0.0,  1.0,    0.0
        ESolver: {
   B: 0.0033  //larger radius
         Type: Analytic
      }
        }
   }

Mode:  {
   Port: {
        ReferenceNumber: 2
        Origin:   Waveguide type: Circular
   0.0, 0.0, 0.1
        XDirection: 1.0,  0.0,    0.0
       Mode typeYDirection: 0.0, TE 1 1.0,    0.0
        ESolver: {
                AType: 0.03Analytic
                }
Mode:  {
      }
   }
{code}

h3. Load two TE modes in the same rectangular waveguide
{code}
  Port: {
   Waveguide   Reference numbertype: 9  // FPCCircular
      Origin: 0.0, 0.198907, -0.4479152585
      XDirection:  -1.0, 0.0, 0.0
     Mode YDirectiontype: 0.0, 0.0,TE 1 1.0
           
      ESolver:  {
      A: 0.03
      Type: Analytic
             Mode: {}
        }
   }

  Port: {
      Reference WaveguideTypenumber: Rectangular
9  // FPC
      Origin: 0.0, 0.198907, -0.4479152585
         ModeTypeXDirection: TE -1 1.0, 0.0, 0.0
      YDirection: 0.0, 0.0, 1.0
      ESolver:  {
    A: 0.1348935946  
      Type: Analytic
              BMode: 0.024973714999999970 {
            }
     }
  }

  PortWaveguideType: { Rectangular
      Reference number: 9  // FPC
      Origin: 0.0, 0.198907, -0.4479152585
 ModeType: TE 1   XDirection:  -1.0, 0.0, 0.0
1
        YDirection: 0.0, 0.0, 1.0          A: 0.1348935946
      ESolver:  {
             TypeB: Analytic0.024973714999999970
            }
 Mode:  {
  }
  }

  Port: {
      Reference number: 9      WaveguideType: Rectangular// FPC
      Origin: 0.0, 0.198907, -0.4479152585
            ModeTypeXDirection: TE 2 -1.0, 0.0, 0.0
      YDirection: 0.0, 0.0, 1.0
      ESolver: {
     A: 0.1348935946  
     Type: Analytic
               BMode: 0.024973714999999970  {
            }
       }
  WaveguideType: }
{code}

h3. Make a non-planar surface absorbing boundary
{code}
  Port: {
Rectangular
             ReferenceNumber: 5       ModeType: TE 2  0
   //reference surface ID
    Origin:     0.0, 0.0, 0.0     A:  //not used
0.1348935946
      XDirection: 1.0, 0.0,    0.0    //not used
    YDirectionB: 0.0, 1.0,024973714999999970
     0.0    //not used
    ESolver: {}
      Type: Analytic}
   }

  Port: Mode:{
    ReferenceNumber: 5              //reference    Mode number: 1surface ID
    Origin:     0.0, 0.0, 0.0       //not used
    Waveguide typeXDirection: ABC
 1.0, 0.0,    0.0    //not used
    YDirection: 0.0, 1.0,    0.0     Mode type//not used
    ESolver: ABC{
      Type: }Analytic
    }
  }
{code}

h3. LinearSolver options in EigenSolver container
* The first option is that user does not provide anything. The EigenSolver container in the input file looks like:
{code}
  EigenSolver: {
    NumEigenvalues:Mode:{
                       Mode number: 1
            1
    FrequencyShift:       Waveguide type: ABC
    10.e9
  }
{code} 
 In this case, Omega3P will use the default option for linear solver for solving shifted linear systems

* The second option is to use Krylov subspace method with different preconditioner.
{code}Mode type: ABC
      }
    }
  }