From 5b34251638f8cac5f798d80e96ad3be8cac20ba4 Mon Sep 17 00:00:00 2001 From: vsitaraman Date: Tue, 31 Mar 2026 10:16:28 -0400 Subject: [PATCH] modified: anasen_fem/run.py modified: anasen_fem/scalars.dat.names modified: anasen_fem/wires_gmsh2d_bc.py --- anasen_fem/run.py | 20 +++++++++++++------- anasen_fem/scalars.dat.names | 4 ++-- anasen_fem/wires_gmsh2d_bc.py | 4 ++-- 3 files changed, 17 insertions(+), 11 deletions(-) diff --git a/anasen_fem/run.py b/anasen_fem/run.py index 8034e1c..24c9ff2 100755 --- a/anasen_fem/run.py +++ b/anasen_fem/run.py @@ -1,15 +1,21 @@ +import code import os -#val=-178.3 -val=17.83 -count=11 -while val<178.3+0.1: +val=-178.3 +# val=17.83 +count=0 +while val<17.83+0.1: print(val) os.system("python3 wires_gmsh2d_bc.py "+str(val)) - os.system("ElmerGrid 14 2 wires2d.msh") - os.system("ElmerSolver wires2d.sif") + os.system("cp wires2d/mesh_files/wires2d.msh wires2d/mesh_files/%02d_%1.4f.msh"%(count,val)) + os.system("ElmerGrid 14 2 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("ElmerSolver wires2d/sif_files/wires2d_%02d_%1.4f.sif"%(count,val)) os.system("./paraview_plotter.py") - os.system("cp contour_output.png contour_output_z_%02d_%1.4f.png"%(count,val)) + os.system("cp contour_output.png png/contour_output_z_%02d_%1.4f.png"%(count,val)) + # os.system("cp field_output.png field_ouput_z_%02d_%1.4f.png"%(count,val)) + val=val+17.83 count = count + 1 + break diff --git a/anasen_fem/scalars.dat.names b/anasen_fem/scalars.dat.names index d1fd517..f2aeb7a 100755 --- a/anasen_fem/scalars.dat.names +++ b/anasen_fem/scalars.dat.names @@ -1,8 +1,8 @@ Metadata for SaveScalars file: ./scalars.dat Elmer version: 26.1 Elmer compilation date: 2026-03-15 -Solver input file: wires2d.sif -File started at: 2026/03/23 20:55:05 +Solver input file: wires2d/sif_files/wires2d_00_-178.3000.sif +File started at: 2026/03/31 10:12:11 Variables in columns of matrix: 1: res: potential difference diff --git a/anasen_fem/wires_gmsh2d_bc.py b/anasen_fem/wires_gmsh2d_bc.py index 7ec4144..8753a1f 100755 --- a/anasen_fem/wires_gmsh2d_bc.py +++ b/anasen_fem/wires_gmsh2d_bc.py @@ -147,10 +147,10 @@ gmsh.model.setPhysicalName(1,30,"cathode_wires") gmsh.model.addPhysicalGroup(2,[anasen_barrel],tag=13) gmsh.model.setPhysicalName(1,13,"gas") -gmsh.option.setNumber("Mesh.Algorithm", 6) +gmsh.option.setNumber("Mesh.Algorithm", 5) gmsh.model.mesh.generate(dim=2) gmsh.model.mesh.refine() -gmsh.write("wires2d.msh") +gmsh.write("wires2d/mesh_files/wires2d.msh") #gmsh.fltk.run() gmsh.finalize()