rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Merge branch 'devel_vignesh' of https://fsunuc.physics.fsu.edu/git/rtang/ANASEN_analysis into devel_vignesh

Browse Source
This commit is contained in:
Vignesh Sitaraman 2026-04-27 18:24:06 -04:00
parent 4a7dc22acd e41f2e14ec
commit 15e60a0cc1
5 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
.gitignore vendored
View File

@ -34,3 +34,4 @@ anasen_fem/scalars.dat.names
myenv/ myenv/
eloss_calculations/ eloss_calculations/
anasen_fem/He96_CO2_4_260Torr.gas anasen_fem/He96_CO2_4_260Torr.gas
anasen_fem/heavy_ion_track.csv

5
anasen_fem/garfield_sim.py
View File

@ -42,7 +42,7 @@ else:
gas.LoadGasFile(gas_file) gas.LoadGasFile(gas_file)
# --- 3. FIELD MAP SETUP --- # --- 3. FIELD MAP SETUP ---
fm = ROOT.Garfield.ComponentElmer() fm = ROOT.Garfield.ComponentElmer2d()
fm.Initialise("wires2d/mesh.header", fm.Initialise("wires2d/mesh.header",
@ -58,6 +58,7 @@ fm.SetMedium(0, gas)
# --- 4. SENSOR AND DRIFT SETUP --- # --- 4. SENSOR AND DRIFT SETUP ---
sensor = ROOT.Garfield.Sensor() sensor = ROOT.Garfield.Sensor()
sensor.AddComponent(fm) sensor.AddComponent(fm)
sensor.SetArea(-50.0, -50.0, -5.0, 50.0, 50.0, 5.0) #hardcoding the sesnsor area to define a psuedo 3d geometry
# Heavy Ion Drift (RKF) - Best for the general track # Heavy Ion Drift (RKF) - Best for the general track
drift = ROOT.Garfield.DriftLineRKF() drift = ROOT.Garfield.DriftLineRKF()
@ -69,7 +70,7 @@ aval.SetSensor(sensor)
# --- 5. EXECUTION --- # --- 5. EXECUTION ---
# Starting position (e.g., near the IC wires at r=23mm or closer to Anodes) # Starting position (e.g., near the IC wires at r=23mm or closer to Anodes)
x0, y0, z0, t0 = 35.0, 0.0, 0.0, 0.0 x0, y0, z0, t0 = 3.50, 0.0, 0.0, 0.0
print(f"Simulating heavy ion drift from r={x0}...") print(f"Simulating heavy ion drift from r={x0}...")
drift.DriftIon(x0, y0, z0, t0) drift.DriftIon(x0, y0, z0, t0)

2
anasen_fem/run.py
View File

@ -10,7 +10,7 @@ while val<178.3+0.1:
os.system("ElmerGrid 14 2 wires2d.msh -2d") os.system("ElmerGrid 14 2 wires2d.msh -2d")
os.system("ElmerSolver wires2d.sif") os.system("ElmerSolver wires2d.sif")
os.system("./paraview_plotter.py") os.system("./paraview_plotter.py")
# os.system("python3 garfield_sim.py") os.system("python3 garfield_sim.py")
os.system("cp wires2d.msh wires2d/mesh_files/wires2d%02d_%1.4f.msh"%(count,val)) os.system("cp wires2d.msh wires2d/mesh_files/wires2d%02d_%1.4f.msh"%(count,val))
os.system("cp wires2d.sif wires2d/sif_files/wires2d_%02d_%1.4f.sif"%(count,val)) os.system("cp wires2d.sif wires2d/sif_files/wires2d_%02d_%1.4f.sif"%(count,val))
os.system("cp wires2d/elfield_anasen_t0001.vtu wires2d/vtu_files/elfield_anasen_%02d_%1.4f.vtu"%(count,val)) os.system("cp wires2d/elfield_anasen_t0001.vtu wires2d/vtu_files/elfield_anasen_%02d_%1.4f.vtu"%(count,val))

2
anasen_fem/scalars.dat.names
View File

@ -2,7 +2,7 @@ Metadata for SaveScalars file: ./scalars.dat
Elmer version: 26.1 Elmer version: 26.1
Elmer compilation date: 2026-03-15 Elmer compilation date: 2026-03-15
Solver input file: wires2d.sif Solver input file: wires2d.sif
File started at: 2026/04/02 15:10:11 File started at: 2026/04/27 17:44:16
Variables in columns of matrix: Variables in columns of matrix:
1: res: potential difference 1: res: potential difference

2
anasen_fem/wires_gmsh2d_bc.py
View File

@ -224,7 +224,7 @@ gmsh.option.setNumber("Mesh.Algorithm", 6)
gmsh.model.mesh.generate(dim=2) gmsh.model.mesh.generate(dim=2)
# gmsh.model.mesh.refine() # gmsh.model.mesh.refine()
# gmsh.model.mesh.refine() # gmsh.model.mesh.refine()
gmsh.model.mesh.setOrder(3)
gmsh.write("wires2d.msh") gmsh.write("wires2d.msh")
gmsh.model.mesh.setOrder(1)
#gmsh.fltk.run() #gmsh.fltk.run()
gmsh.finalize() gmsh.finalize()