def Update_Moose_after_Init(Sim_name, Mesh, MooseFile, Porosity_File): SatInitName = 'dataFile = PrimaryFiles/' + Sim_name + '/MOOSEValues_sat_updated.txt' PressInitName = 'dataFile = PrimaryFiles/' + Sim_name + '/MOOSEValues_press_updated.txt' PorosityInitName = 'dataFile = PrimaryFiles/' + Sim_name + '/' + Porosity_File + '.txt' lines_to_change = [168, 173, 178] new_lines = [SatInitName, PressInitName, PorosityInitName] change_input(MooseFile, lines_to_change, new_lines)
def Update_Porosity(Sim_name, LammpsFileDir, Porosity_file, Mesh, nparticles, radius, nnodes, nelements, ann_r, dom_r, NN): input_path = "string input_path = \"/home/crhea/Dropbox/Thesis/Mesh/\";" input_path_lammps_Data = "string input_path_lammps_Data = \"/home/crhea/Dropbox/Thesis/PrimaryFiles/" + Sim_name + "/" + LammpsFileDir + "/\";" output_path = "string output_path = \"/home/crhea/Dropbox/Thesis/PrimaryFiles/" + Sim_name + "/\";" mesh_to_read = "string mesh_to_read = \"" + Mesh + "_nodes\";" mesh_connectivity = "string mesh_connectivity = \"" + Mesh + "_connectivity\";" file_name = "string file_name = \"lammps_pos_out_simple\";" output = "string outputname = \"Porosity\";" domain_radius = 'double domain_radius = ' + str(dom_r) + ';' annulus_radius = 'double annulus_radius = ' + str(ann_r) + ';' nnodes = 'int nnodes =' + str(nnodes) + ';' number_cell_elements = 'int number_cell_elements = ' + str(nelements) + ';' num_of_NN = 'int num_of_NN = ' + str(NN) + ';' num_partic_lines = 'int nparticles = ' + str(nparticles) + ';' rad = 'double particle_rad = ' + str(radius) + ';' lines_to_change = [ 638, 639, 640, 641, 642, 643, 644, 645, 646, 647, 649, 651, 653, 654 ] new_lines = [ input_path, input_path_lammps_Data, output_path, mesh_to_read, mesh_connectivity, file_name, output, domain_radius, annulus_radius, num_partic_lines, nnodes, number_cell_elements, num_of_NN, rad ] change_input( '/home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' + Porosity_file + '.cpp', lines_to_change, new_lines)
def Update_Porosity(Sim_name, Porosity_file, Mesh, nparticles, radius, nnodes, nelements, ann_r, dom_r, NN): input_path = "string input_path = \"/home/crhea/Dropbox/Thesis/Mesh/\";" input_path_lammps_Data = "string input_path_lammps_Data = \"/home/crhea/Dropbox/Thesis/PrimaryFiles/" + Sim_name + "/LAMMPSFILES/\";" output_path = "string output_path = \"/home/crhea/Dropbox/Thesis/PrimaryFiles/" + Sim_name + "/\";" mesh_to_read = "string mesh_to_read = \"" + Mesh + "_nodes\";" mesh_connectivity = "string mesh_connectivity = \"" + Mesh + "_connectivity\";" file_name = "string file_name = \"lammps_pos_out_simple\";" output = "string outputname = \"Porosity\";" domain_radius = 'double domain_radius = ' + str(dom_r) + ';' annulus_radius = 'double annulus_radius = ' + str(ann_r) + ';' nnodes = 'int nnodes =' + str(nnodes) + ';' number_cell_elements = 'int number_cell_elements = ' + str(nelements) + ';' num_of_NN = 'int num_of_NN = ' + str(NN) + ';' num_partic_lines = 'int nparticles = ' + str(nparticles) + ';' rad = 'double particle_rad = ' + str(radius) + ';' lines_to_change = [ 686, 687, 688, 689, 690, 691, 692, 693, 694, 695, 697, 699, 701, 702 ] new_lines = [ input_path, input_path_lammps_Data, output_path, mesh_to_read, mesh_connectivity, file_name, output, domain_radius, annulus_radius, num_partic_lines, nnodes, number_cell_elements, num_of_NN, rad ] change_input(Porosity_file + '.cpp', lines_to_change, new_lines)
def Update_MOOSE(Sim_name, MooseFile, Porosity_File, time_Step): WorkingDirectory = '/home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' PorosityInitName = 'dataFile = ' + WorkingDirectory + Porosity_File + '.txt' PorosityInitNameOld = 'dataFile = ' + WorkingDirectory + Porosity_File + '_old.txt' #OutputName = 'file_base = PrimaryFiles/'+Sim_name+'/MOOSEFILES/MOOSEOutput'+str(time_Step) lines_to_change = [120, 125] #,237] new_lines = [PorosityInitName, PorosityInitNameOld] #,OutputName] change_input(MooseFile, lines_to_change, new_lines)
def Update_Lammps(Sim_name, LammpsFile): RestartRead = 'read_restart /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/LAMMPSFILES/restart1' Dump = 'dump pos all custom 100 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/LAMMPSFILES/pos_lammps_out.txt id type x y z vx vy vz' Restart = 'restart 100 PrimaryFiles/' + Sim_name + '/LAMMPSFILES/restart1 PrimaryFiles/' + Sim_name + '/LAMMPSFILES/restart2' Fix1 = 'fix 1 all viscous/field 0.5 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/VelForLammpsX.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/VelForLammpsY.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/Viscosity.txt' lines_to_change = [1, 31, 32, 35] new_line_lammps2 = [RestartRead, Dump, Restart, Fix1] change_input(LammpsFile, lines_to_change, new_line_lammps2)
def Update_Lammps(Sim_name, LammpsFileDir, LammpsFile): RestartRead = 'read_restart /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' + LammpsFileDir + '/restart1' Dump = 'dump pos all custom 1 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' + LammpsFileDir + '/pos_lammps_out.txt id type x y z vx vy vz' Restart = 'restart 1 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' + LammpsFileDir + '/restart1 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' + LammpsFileDir + '/restart2' Fix1 = 'fix 1 all viscous/field 50.0 /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/VelForLammpsX.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/VelForLammpsY.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/Viscosity.txt' Hookean = 'pair_style hooke/cap ${kn} ${kt} ${gamma_n} ${gamma_t} ${coeffFric} 0 10.0 10.0 0.78 0.99 50.0 ' + '/home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/SaturationInterpolated.txt' lines_to_change = [1, 34, 35, 40, 44] new_line_lammps2 = [RestartRead, Dump, Restart, Hookean, Fix1] change_input(LammpsFile, lines_to_change, new_line_lammps2)
def Update_Moose_after_Init(Sim_name, Mesh, MooseFile, Porosity_File): WorkingDirectory = '/home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' SatInitName = 'dataFile = ' + WorkingDirectory + 'MOOSEValues_sat_updated.txt' PressInitName = 'dataFile = ' + WorkingDirectory + 'MOOSEValues_press_updated.txt' PorosityInitName = 'dataFile = ' + WorkingDirectory + Porosity_File + '.txt' PorosityInitNameOld = 'dataFile = ' + WorkingDirectory + Porosity_File + '.txt' lines_to_change = [110, 115, 120, 125] new_lines = [ SatInitName, PressInitName, PorosityInitName, PorosityInitNameOld ] change_input(MooseFile, lines_to_change, new_lines)
def Update_Lammps_Init(Sim_name,LammpsFile,ParticlesInput,Diameter,DomainVals): DiameterData = 'variable Diam equal '+str(Diameter) InitData = 'read_data '+ParticlesInput Dump = 'dump pos all custom 100 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/LAMMPSFILES/pos_lammps_out.txt id type x y z vx vy vz' Restart = 'restart 100 PrimaryFiles/'+Sim_name+'/LAMMPSFILES/restart1 PrimaryFiles/'+Sim_name+'/LAMMPSFILES/restart2' Fix1 = 'fix 1 all viscous/field 0.5 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/VelForLammpsX.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/VelForLammpsY.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/Viscosity.txt' if len(Domain) == 4: #Rectangle DomainSet = 'region box block '+str(DomainVals[0])+' '+str(DomainVals[1]) + ' '+str(DomainVals[2])+' '+str(DomainVals[3])+' -3 3 units box' lines_to_change = [2,29,30,43,49,53] else: #cirlce pass new_line_lammps = [DiameterData,InitData,DomainSet,Dump,Restart,Fix1] change_input(LammpsFile+"_init",lines_to_change,new_line_lammps)
def Update_Lammps_Init(Sim_name,LammpsFileDir,LammpsFile,ParticlesInput,Diameter,DomainVals): DiameterData = 'variable Diam equal '+str(Diameter) InitData = 'read_data '+ParticlesInput Dump = 'dump pos all custom 100 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/'+LammpsFileDir+'/pos_lammps_out.txt id type x y z vx vy vz' Restart = 'restart 100 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/'+LammpsFileDir+'/restart1 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/'+LammpsFileDir+'/restart2' Fix1 = 'fix 1 all viscous/field 0.5 /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/VelForLammpsX.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/VelForLammpsY.txt /home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/Viscosity.txt' Hookean = 'pair_style hooke/cap ${kn} ${kt} ${gamma_n} ${gamma_t} ${coeffFric} 0 1.01 1.1 1.0 0.8 50.0 '+'/home/crhea/Dropbox/Thesis/PrimaryFiles/'+Sim_name+'/SaturationInterpolated.txt' if len(Domain) == 4: #Rectangle DomainSet = 'region box block '+str(DomainVals[0])+' '+str(DomainVals[1]) + ' '+str(DomainVals[2])+' '+str(DomainVals[3])+' -3 3 units box' lines_to_change = [2,29,30,41,42,43,46] else: #cirlce pass new_line_lammps = [DiameterData,InitData,DomainSet,Dump,Restart,Hookean,Fix1] change_input(LammpsFile+"_init",lines_to_change,new_line_lammps)
def Init_FF(Sim_name,MooseFile,Mesh,Porosity_File,Init_Press): Meshchange = 'file = Mesh/'+Mesh+'.e' outputinitSatPress = "PrimaryFiles/"+Sim_name+"/MOOSEValues" InitConst("Mesh/"+Mesh+".e",outputinitSatPress,Init_Press) SatInitName = 'dataFile = PrimaryFiles/'+Sim_name+'/MOOSEValues_sat_init.txt' PressInitName = 'dataFile = PrimaryFiles/'+Sim_name+'/MOOSEValues_press_init.txt' PorosityInitName = 'dataFile = PrimaryFiles/'+Sim_name+'/'+Porosity_File+'.txt' XvelName = ' output = PrimaryFiles/'+Sim_name+'/velocitiesX' YvelName = ' output = PrimaryFiles/'+Sim_name+'/velocitiesY' Sat_Out = ' output = PrimaryFiles/'+Sim_name+'/MOOSEValues_sat_updated' Press_Out = ' output = PrimaryFiles/'+Sim_name+'/MOOSEValues_press_updated' OutputName = 'file_base = PrimaryFiles/'+Sim_name+'/MOOSEFILES/MOOSEOutput' lines_to_change = [3,110,115,120,141,146,151,156,179] new_lines = [Meshchange,SatInitName,PressInitName,PorosityInitName,XvelName,YvelName,Sat_Out,Press_Out,OutputName] change_input(MooseFile,lines_to_change,new_lines)
def Init_FF(Sim_name, MooseFileDir, MooseFile, Mesh, Porosity_File, Init_Press): WorkingDirectory = '/home/crhea/Dropbox/Thesis/PrimaryFiles/' + Sim_name + '/' Meshchange = 'file = ' + Mesh + '.e' outputinitSatPress = "MOOSEValues" InitConst(Mesh + ".e", outputinitSatPress, Init_Press) SatInitName = 'dataFile = ' + WorkingDirectory + 'MOOSEValues_sat_init.txt' PressInitName = 'dataFile = ' + WorkingDirectory + 'MOOSEValues_press_init.txt' PorosityInitName = 'dataFile = ' + WorkingDirectory + Porosity_File + '.txt' PorosityInitNameOld = 'dataFile = ' + WorkingDirectory + Porosity_File + '.txt' XvelName = ' output =' + WorkingDirectory + 'velocitiesX' YvelName = ' output = ' + WorkingDirectory + 'velocitiesY' Sat_Out = ' output = ' + WorkingDirectory + 'MOOSEValues_sat_updated' Press_Out = ' output = ' + WorkingDirectory + 'MOOSEValues_press_updated' OutputName = 'file_base = ' + WorkingDirectory + MooseFileDir + '/MOOSEOutput' lines_to_change = [3, 110, 115, 120, 125, 146, 151, 156, 161, 184] new_lines = [ Meshchange, SatInitName, PressInitName, PorosityInitName, PorosityInitNameOld, XvelName, YvelName, Sat_Out, Press_Out, OutputName ] change_input(MooseFile, lines_to_change, new_lines)
def Update_MOOSE(Sim_name, MooseFile, Porosity_File, time_Step): PorosityInitName = 'dataFile = PrimaryFiles/' + Sim_name + '/' + Porosity_File + '.txt' #OutputName = 'file_base = PrimaryFiles/'+Sim_name+'/MOOSEFILES/MOOSEOutput'+str(time_Step) lines_to_change = [178] #,237] new_lines = [PorosityInitName] #,OutputName] change_input(MooseFile, lines_to_change, new_lines)