Check Keywords Warn ! Set drift distance, E fields, and voltage between GEMs ! dist: distance in cm ! Edrift: drift field in V/cm ! Etrans: transfer field in V/cm ! Eind: induction field in V/cm ! deltaV1: voltage (in V) between 1st GEM ! deltaV2: voltage (in V) between 2nd GEM $ldist = 0.2 $udist = 0.5 $mdist = 0.1 $Edrift = 50 $Etrans = 500 $Eind = 1000 $deltaV1= 850 $deltaV2= 950 $WTuel = 0 $WTucp = 0 $WTlcp = 0 $WTucp2 = 0 $WTlcp2 = 0 $WTlel = 0 Header Mesh DB "." "Gemcell" End Simulation Coordinate System = Cartesian 3D Simulation Type = Steady State Steady State Max Iterations = 1 Output File = "Gemcell.result" Post File = "Gemcell.ep" End Constants Permittivity Of Vacuum = 8.8542e-12 End Body 1 Equation = 1 Material = 1 End Body 2 Equation = 1 Material = 2 End Body 3 Equation = 1 Material = 3 End Body 4 Equation = 1 Material = 3 End Body 5 Equation = 1 Material = 2 End Body 6 Equation = 1 Material = 3 End Body 7 Equation = 1 Material = 3 End Equation 1 Active Solvers(1) = 1 Calculate Electric Energy = True End Solver 1 Equation = Stat Elec Solver Variable = Potential Variable DOFs = 1 Procedure = "StatElecSolve" "StatElecSolver" Calculate Electric Field = True Calculate Electric Flux = False Linear System Solver = Iterative Linear System Iterative Method = BiCGStab Linear System Max Iterations = 1000 Linear System Abort Not Converged = True Linear System Convergence Tolerance = 1.0e-10 Linear System Preconditioning = ILU1 Steady State Convergence Tolerance = 5.0e-7 ! Adaptive Mesh Refinement = True ! Adaptive Remesh = True ! Adaptive Save Mesh = True ! Adaptive Error Limit = 1.0e-12 End ! Gas for both 1st and 2nd GEM Material 1 Relative Permittivity = 1 Density = 1 End ! Dielectric for both 1st and 2nd GEM Material 2 Relative Permittivity = 3.23 Density = 2 End ! Copper Material 3 for both 1st and 2nd GEM Relative Permittivity = 1.0e10 Density = 3 End ! Upper copper plate of 1st GEM Boundary Condition 1 Target Boundaries = 1 Potential = $-1*Etrans*mdist*Eind*ldist - deltaV1 + WTucp End ! Lower copper plate Boundary Condition 2 Target Boundaries = 2 Potential = $-1*Etrans*mdist*Eind*ldist + WTlcp End ! Upper copper plate of 2nd GEM Boundary Condition 3 Target Boundaries = 3 Potential = $-1*Eind*ldist - deltaV2 + WTucp2 End ! Lower copper plate of 2nd GEM Boundary Condition 4 Target Boundaries = 4 Potential = $-1*Eind*ldist + WTlcp2 End ! Upper electrode Boundary Condition 5 Target Boundaries = 5 Potential = $-1*Eind*ldist*Etrans*mdist - deltaV1 - deltaV2 - Edrift*udist + WTuel End ! Lower electrode Boundary Condition 6 Target Boundaries = 6 Potential = $WTlel End ! Set up boundary A for hole 1 for both 1st and 2nd GEM Boundary Condition 7 Target Boundaries = 7 End ! Link to half A of hole 2 for both 1st and 2nd GEM Boundary Condition 8 Target Boundaries = 8 Periodic BC = 7 Periodic BC Potential = Logical True End ! Set up boundary B for hole 3 for both 1st and 2nd GEM Boundary Condition 9 Target Boundaries = 9 End ! Link to half B of hole 2 for both 1st and 2nd GEM Boundary Condition 10 Target Boundaries = 10 Periodic BC = 9 Periodic BC Potential = Logical True End ! Set up boundary C for hole 1 side for both 1st and 2nd GEM Boundary Condition 11 Target Boundaries = 11 End ! Link to the side containing hole for 3 both 1st and 2nd GEM Boundary Condition 12 Target Boundaries = 12 Periodic BC = 11 Periodic BC Potential = Logical True End