Elmer solver input file problem plz someone help me! i'm beginner in Garfield++

[Dk_ntm]

i am encountering with this type of error when i creating & running a .sif file for double layer GEM detector.

ELMER SOLVER (v 9.0) STARTED AT: 2024/10/27 15:26:38
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
…

MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
ERROR:: LoadInputFile:
ERROR:: LoadInputFile: Unknown specifier: [target boundaries]
ERROR:: LoadInputFile: In section: [boundary condition 5]
ERROR:: LoadInputFile: For property name:[boundary condition 6]
STOP 1
here is the sif file. plz suggest me!
siffile.ipynb (12.4 KB)

[hschindl]

This is probably more a question for the Elmer forum, but perhaps @jrenner can help.
Could you also upload your .sif file?

[Dk_ntm]

Here is the .sif file
sif.ipynb (6.2 KB)

[hschindl]

Isn’t that a python notebook that you use to create the .sif file, not the .sif itself?

[Dk_ntm]

It is a .sif file

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 = {dist_lower_electrode}
$udist = {dist_upper_electrode}
$mdist = {dist_mid_electrode}
$Edrift = {E_drift}
$Etrans = {E_trans}
$Eind = {E_ind}
$deltaV1= {V_GEM1}
$deltaV2= {V_GEM2}
$WTuel = {wt_uel}
$WTucp = {wt_ucp}
$WTlcp = {wt_lcp}
$WTucp2 = {wt_ucp2}
$WTlcp2 = {wt_lcp2}
$WTlel = {wt_lel}
“”"

gemcell_sif_text = “”"
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
Material 1
Relative Permittivity = 1
Density = 1
End

! Dielectric
Material 2
Relative Permittivity = 3.23
Density = 2
End

! Copper
Material 3
Relative Permittivity = 1.0e10
Density = 3
End

! Upper copper plate of 1st GEM
Boundary Condition 1
Target Boundaries = 1
Potential = $-1Etransmdist - deltaV1 + WTucp
End

! Lower copper plate
Boundary Condition 2
Target Boundaries = 2
Potential = $-1Etransmdist + WTlcp
End

! Upper copper plate of 2nd GEM
Boundary Condition 3
Target Boundaries = 3
Potential = $-1Eindldist - deltaV2 + WTucp2
End

! Lower copper plate of 2nd GEM
Boundary Condition 4
Target Boundaries = 4
Potential = $-1Eindldist + WTlcp2
End

! Upper electrode
Boundary Condition 5
Target Boundaries = 5
Potential = $-1EindldistEtransmdist - deltaV1 + deltaV2 - Edrift*udist + WTuel
End
! Lower electrode
Boundary Condition 6
Target Boundaries = 6
Potential = $WTlel
End

! Set up boundary A for hole 1
Boundary Condition 7
Target Boundaries = 7
End
! Link to half A of hole 2
Boundary Condition 8
Target Boundaries = 8
Periodic BC = 7
Periodic BC Potential = Logical True
End

! Set up boundary B for hole 3
Boundary Condition 9
Target Boundaries = 9
End
! Link to half B of hole 2
Boundary Condition 10
Target Boundaries = 10
Periodic BC = 9
Periodic BC Potential = Logical True
End

! Set up boundary C for hole 1 side
Boundary Condition 11
Target Boundaries = 11
End
! Link to the side containing hole 3
Boundary Condition 12
Target Boundaries = 12
Periodic BC = 11
Periodic BC Potential = Logical True
End

[jrenner]

In the Jupyter notebook that was sent, the keyword End appears to be missing after the declaration of Boundary Condition 5:

! Upper electrode
Boundary Condition 5
  Target Boundaries = 5
  Potential = $-1*Eind*ldist*Etrans*mdist  - deltaV1  + deltaV2 - Edrift*udist + WTuel
  
! Lower electrode
Boundary Condition 6
  Target Boundaries = 6
  Potential = $WTlel
End

This seems to have been corrected in the version you sent in the previous message. Did that solve the problem or change the error output?

[Dk_ntm]

yeah thank you! i corrected it but when i running this .sif file, i’m encountering with this problem. And also tell me is my .sif file is correct for double Layer GEM Detector?
ELMER SOLVER (v 9.0) STARTED AT: 2024/10/29 16:28:56

ParCommInit: Initialize #PEs: 1

MAIN:

MAIN: =============================================================

MAIN: ElmerSolver finite element software, Welcome!

MAIN: This program is free software licensed under (L)GPL

MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.

MAIN: Webpage Elmer - Elmer - CSC Company Site, Email elmeradm@csc.fi

MAIN: Version: 9.0 (Rev: unknown, Compiled: 2024-08-07)

MAIN: Running one task without MPI parallelization.

MAIN: Running with just one thread per task.

MAIN: HYPRE library linked in.

MAIN: MUMPS library linked in.

MAIN: Lua interpreter linked in.

MAIN: =============================================================

MAIN:

MAIN:

MAIN: -------------------------------------

MAIN: Reading Model: case.sif

LoadInputFile: Scanning input file: case.sif

LoadInputFile: Scanning only size info

LoadInputFile: First time visiting

LoadInputFile: Reading base load of sif file

LoadInputFile: Loading input file: case.sif

LoadInputFile: Reading base load of sif file

CheckKeyword: Unlisted keyword: [periodic bc potential] in section: [boundary condition 8]

CheckKeyword: Unlisted keyword: [periodic bc potential] in section: [boundary condition 10]

CheckKeyword: Unlisted keyword: [periodic bc potential] in section: [boundary condition 12]

LoadInputFile: Number of BCs: 12

LoadInputFile: Number of Body Forces: 0

LoadInputFile: Number of Initial Conditions: 0

LoadInputFile: Number of Materials: 3

LoadInputFile: Number of Equations: 1

LoadInputFile: Number of Solvers: 1

LoadInputFile: Number of Bodies: 7

ListTagKeywords: Setting weight for keywords!

ListTagKeywords: No parameters width suffix: normalize by area

ListTagKeywords: Setting weight for keywords!

ListTagKeywords: No parameters width suffix: normalize by volume

Loading user function library: [StatElecSolve]…[StatElecSolver_Init0]

ERROR:: LoadMesh: Requested mesh > ./Gemcell < does not exist!

STOP 1

i have attached the .sif file plz someone help me
sif.txt (3.5 KB)

[jrenner]

It looks like it can’t find the mesh Gemcell. Have you created the mesh in Gmsh and converted it to Elmer format?

The .sif file alone will not simulate a double layer GEM if the corresponding geometry has not been created and meshed correctly.