 |
KratosMultiphysics
KRATOS Multiphysics (Kratos) is a framework for building parallel, multi-disciplinary simulation software, aiming at modularity, extensibility, and high performance. Kratos is written in C++, and counts with an extensive Python interface.
|
Variables | |
| string | file_name = 'sp_data.hdf5' |
| f = h5py.File(file_name, 'r') | |
| test_id = f.attrs['test_id'] | |
| internal_radius = f.attrs['internal_radius'] | |
| external_radius = f.attrs['external_radius'] | |
| interface_radius = f.attrs['interface_radius'] | |
| thickness = f.attrs['thickness'] | |
| volume = f.attrs['volume'] | |
| real_probe_height = f.attrs['real_probe_height'] | |
| target_porosity = f.attrs['target_porosity'] | |
| porosity = f.attrs['porosity'] | |
| density = f.attrs['density'] | |
| height_factor = real_probe_height / thickness | |
| tuple | porosity_factor = (1.0-target_porosity)/(1.0-porosity) |
| float | gram_factor = 1000.0 |
| initial_radii = np.array(f['0'].get('radius')) | |
| initial_radii_3 = np.power(initial_radii,3) | |
| float | initial_masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * initial_radii_3 |
| initial_mass = np.sum(initial_masses) | |
| list | all_times = [] |
| list | all_pressures = [] |
| list | all_sps = [] |
| int | failure_step = 175 |
| int | max_step = 2*failure_step |
| times = np.zeros(max_step+1) | |
| int | p_rate = 5e10 |
| float | psi_factor = 0.000145038 |
| list | pressures = [p_rate*t*psi_factor for t in times] |
| t_f = times[failure_step] | |
| float | max_num_bonds = 6.0 |
| sps = np.zeros(max_step+1) | |
| all_radii = np.array(f[str(i)].get('radius')) | |
| continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds')) | |
| xs = np.array(f[str(i)].get('x')) | |
| ys = np.array(f[str(i)].get('y')) | |
| xs_2 = np.power(xs,2) | |
| ys_2 = np.power(ys,2) | |
| distance_2 = xs_2 + ys_2 | |
| weak_radius = interface_radius | |
| internal_radii = np.where(distance_2<weak_radius**2,all_radii,0.0) | |
| external_radii = np.where(distance_2>=weak_radius**2,all_radii,0.0) | |
| cont_internal_radii = np.where(continuum_bonds>numbonds,internal_radii,0.0) | |
| cont_external_radii = np.where(continuum_bonds>0.0,external_radii,0.0) | |
| total_radii_3 = np.power(cont_internal_radii+cont_external_radii,3) | |
| float | masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * total_radii_3 |
| current_mass = np.sum(masses) | |
| current_sp = initial_mass - current_mass | |
| list | exp_times = [] |
| list | exp_pressures = [] |
| list | exp_sps = [] |
| fields = line.strip().split() | |
| string | graph_name = 'sp_bonds_t.pdf' |
| list | graph_labels |
| productions | |
| label | |
| loc | |
| prop | |
| bbox_inches | |
| list sp_statistics.all_pressures = [] |
| sp_statistics.all_radii = np.array(f[str(i)].get('radius')) |
| list sp_statistics.all_sps = [] |
| list sp_statistics.all_times = [] |
| sp_statistics.bbox_inches |
| sp_statistics.cont_external_radii = np.where(continuum_bonds>0.0,external_radii,0.0) |
| sp_statistics.cont_internal_radii = np.where(continuum_bonds>numbonds,internal_radii,0.0) |
| sp_statistics.continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds')) |
| sp_statistics.current_mass = np.sum(masses) |
| sp_statistics.current_sp = initial_mass - current_mass |
| sp_statistics.density = f.attrs['density'] |
| list sp_statistics.exp_pressures = [] |
| list sp_statistics.exp_sps = [] |
| list sp_statistics.exp_times = [] |
| sp_statistics.external_radii = np.where(distance_2>=weak_radius**2,all_radii,0.0) |
| sp_statistics.external_radius = f.attrs['external_radius'] |
| sp_statistics.f = h5py.File(file_name, 'r') |
| int sp_statistics.failure_step = 175 |
| sp_statistics.fields = line.strip().split() |
| string sp_statistics.file_name = 'sp_data.hdf5' |
| float sp_statistics.gram_factor = 1000.0 |
| list sp_statistics.graph_labels |
| string sp_statistics.graph_name = 'sp_bonds_t.pdf' |
| sp_statistics.height_factor = real_probe_height / thickness |
| sp_statistics.initial_mass = np.sum(initial_masses) |
| float sp_statistics.initial_masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * initial_radii_3 |
| sp_statistics.initial_radii = np.array(f['0'].get('radius')) |
| sp_statistics.initial_radii_3 = np.power(initial_radii,3) |
| sp_statistics.interface_radius = f.attrs['interface_radius'] |
| sp_statistics.internal_radii = np.where(distance_2<weak_radius**2,all_radii,0.0) |
| sp_statistics.internal_radius = f.attrs['internal_radius'] |
| sp_statistics.label |
| sp_statistics.loc |
| float sp_statistics.masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * total_radii_3 |
| float sp_statistics.max_num_bonds = 6.0 |
| int sp_statistics.max_step = 2*failure_step |
| int sp_statistics.p_rate = 5e10 |
| sp_statistics.porosity = f.attrs['porosity'] |
| tuple sp_statistics.porosity_factor = (1.0-target_porosity)/(1.0-porosity) |
| sp_statistics.pressures = [p_rate*t*psi_factor for t in times] |
| sp_statistics.productions |
| sp_statistics.prop |
| float sp_statistics.psi_factor = 0.000145038 |
| sp_statistics.real_probe_height = f.attrs['real_probe_height'] |
| sp_statistics.sps = np.zeros(max_step+1) |
| sp_statistics.t_f = times[failure_step] |
| sp_statistics.target_porosity = f.attrs['target_porosity'] |
| sp_statistics.test_id = f.attrs['test_id'] |
| sp_statistics.thickness = f.attrs['thickness'] |
| sp_statistics.times = np.zeros(max_step+1) |
| sp_statistics.total_radii_3 = np.power(cont_internal_radii+cont_external_radii,3) |
| sp_statistics.volume = f.attrs['volume'] |
| sp_statistics.weak_radius = interface_radius |
| sp_statistics.xs = np.array(f[str(i)].get('x')) |
| sp_statistics.xs_2 = np.power(xs,2) |
| sp_statistics.ys = np.array(f[str(i)].get('y')) |
| sp_statistics.ys_2 = np.power(ys,2) |
1.9.1