pumapy.physics_models.finite_volume

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.assign_prescribed_bc_cy(not_dir_x, not_dir_y, not_dir_z, Td, dirichlet_bc_xfaces, dirichlet_bc_yfaces, dirichlet_bc_zfaces, len_x, len_y, len_z, i_iv)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.create_Ab_indices_cy(I_A, J_A, I_b, counter_A, counter_b, i_cv, len_x, len_y, len_z, len_xyz, side_bc)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.create_T_ivs_cy(T, tf, i_cv, len_x, len_y, len_z, len_xyz, side_bc, t_sw, t_se, t_nw, t_ne, t_tsw, t_tse, t_tnw, t_tne)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.index_at_p(index, size)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.index_at_s(index, size)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.assign_prescribed_bc_cy(not_dir_x, not_dir_y, not_dir_z, Td, dirichlet_bc_xfaces, dirichlet_bc_yfaces, dirichlet_bc_zfaces, len_x, len_y, len_z, i_iv)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.create_Ab_indices_cy(I_A, J_A, I_b, counter_A, counter_b, i_cv, len_x, len_y, len_z, len_xyz, side_bc)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.create_T_ivs_cy(T, tf, i_cv, len_x, len_y, len_z, len_xyz, side_bc, t_sw, t_se, t_nw, t_ne, t_tsw, t_tse, t_tnw, t_tne)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.index_at_p(index, size)

pumapy.physics_models.finite_volume.anisotropic_conductivity_utils.index_at_s(index, size)

pumapy.physics_models.finite_volume.elasticity_utils

pumapy.physics_models.finite_volume.elasticity_utils.assign_prescribed_bc_cy(not_dir_x, not_dir_y, not_dir_z, Dd, dirichlet_bc_xfaces, dirichlet_bc_yfaces, dirichlet_bc_zfaces, len_x, len_y, len_z, i_iv)

pumapy.physics_models.finite_volume.elasticity_utils.create_Ab_indices_cy(I_A, J_A, I_b, counter_A, counter_b, i_cv, len_x, len_y, len_z, len_xyz, side_bc)

pumapy.physics_models.finite_volume.elasticity_utils.create_u_ivs_cy(u, uf, i_cv, len_x, len_y, len_z, len_xyz, side_bc, u_sw, u_se, u_nw, u_ne, u_tsw, u_tse, u_tnw, u_tne)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_p(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_p_pad(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_s(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_s_pad(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.pad_domain_cy(ws_pad, orient_pad, need_to_orient, len_x, len_y, len_z, side_bc)

pumapy.physics_models.finite_volume.elasticity_utils

pumapy.physics_models.finite_volume.elasticity_utils.assign_prescribed_bc_cy(not_dir_x, not_dir_y, not_dir_z, Dd, dirichlet_bc_xfaces, dirichlet_bc_yfaces, dirichlet_bc_zfaces, len_x, len_y, len_z, i_iv)

pumapy.physics_models.finite_volume.elasticity_utils.create_Ab_indices_cy(I_A, J_A, I_b, counter_A, counter_b, i_cv, len_x, len_y, len_z, len_xyz, side_bc)

pumapy.physics_models.finite_volume.elasticity_utils.create_u_ivs_cy(u, uf, i_cv, len_x, len_y, len_z, len_xyz, side_bc, u_sw, u_se, u_nw, u_ne, u_tsw, u_tse, u_tnw, u_tne)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_p(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_p_pad(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_s(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.index_at_s_pad(index, size)

pumapy.physics_models.finite_volume.elasticity_utils.pad_domain_cy(ws_pad, orient_pad, need_to_orient, len_x, len_y, len_z, side_bc)

pumapy.physics_models.finite_volume.isotropic_conductivity

class pumapy.physics_models.finite_volume.isotropic_conductivity.IsotropicConductivity(workspace, cond_map, direction, side_bc, prescribed_bc, tolerance, maxiter, solver_type, display_iter, matrix_free)[source]

Bases: pumapy.physics_models.utils.linear_solvers.PropertySolver

assemble_Amatrix()[source]

assemble_bvector()[source]

compute()[source]

compute_effective_coefficient()[source]

error_check()[source]

initialize()[source]

log_input()[source]

log_output()[source]

setup_cond_matrices()[source]

pumapy.physics_models.finite_volume.isotropic_conductivity_utils

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.compute_flux(T, cond, l_x, l_y, l_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.index_at_p(i, j, k, len_x, len_y, len_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.index_at_s(i, j, k, len_x, len_y, len_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.setup_matrices_cy(_kf, l_x, l_y, l_z, domain_bc_check, bc_check, prescribed_bc, side_bc)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.compute_flux(T, cond, l_x, l_y, l_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.index_at_p(i, j, k, len_x, len_y, len_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.index_at_s(i, j, k, len_x, len_y, len_z)

pumapy.physics_models.finite_volume.isotropic_conductivity_utils.setup_matrices_cy(_kf, l_x, l_y, l_z, domain_bc_check, bc_check, prescribed_bc, side_bc)

pumapy.physics_models.finite_volume.mpfa_conductivity

Further explained in publication: Semeraro, F., Ferguson, J.C., Acin, M., Panerai, F. and Mansour, N.N., 2021. Anisotropic analysis of fibrous and woven materials part 2: Computation of effective conductivity. Computational Materials Science, 186, p.109956. https://www.sciencedirect.com/science/article/abs/pii/S092702562030447X

class pumapy.physics_models.finite_volume.mpfa_conductivity.AnisotropicConductivity(workspace, cond_map, direction, side_bc, dirichlet_bc, tolerance, maxiter, solver_type, display_iter)[source]

Bases: pumapy.physics_models.utils.linear_solvers.PropertySolver

assemble_matrices()[source]

compute()[source]

compute_Kmat(i, i_cv)[source]

compute_effective_coefficient()[source]

compute_fluxes()[source]

compute_transmissibility(i, i_cv_in)[source]

error_check()[source]

initialize()[source]

initialize_mpfa()[source]

log_input()[source]

log_output()[source]

pumapy.physics_models.finite_volume.mpsa_elasticity

class pumapy.physics_models.finite_volume.mpsa_elasticity.Elasticity(workspace, elast_map, direction, side_bc, tolerance, maxiter, solver_type, display_iter, dirichlet_bc)[source]

Bases: pumapy.physics_models.utils.linear_solvers.PropertySolver

assemble_matrices()[source]

compute()[source]

compute_Cmat(i, i_cv)[source]

compute_effective_coefficient()[source]

compute_stresses()[source]

compute_transmissibility(i, i_cv_in)[source]

error_check()[source]

initialize()[source]

initialize_mpsa()[source]

log_input()[source]

log_output()[source]

pumapy.physics_models.finite_volume.mpxa_matrices

matrices used in MPFA and MPSA functions

pumapy.physics_models.finite_volume.mpxa_matrices.create_d1_mpfa(k, d)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.create_d1_mpsa(c, d)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.create_d2_mpfa(k, d)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.create_d2_mpsa(c, d)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.create_mpfa_indices()[source]

pumapy.physics_models.finite_volume.mpxa_matrices.create_mpsa_indices()[source]

pumapy.physics_models.finite_volume.mpxa_matrices.div_Ed_mpfa(q_sw, q_se, q_nw, q_ne, q_tsw, q_tse, q_tnw, q_tne)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.div_Ed_mpsa(s_sw, s_se, s_nw, s_ne, s_tsw, s_tse, s_tnw, s_tne)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.div_Eu_mpfa(E_sw, E_se, E_nw, E_ne, E_tsw, E_tse, E_tnw, E_tne)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.div_Eu_mpsa(E_sw, E_se, E_nw, E_ne, E_tsw, E_tse, E_tnw, E_tne)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Ampfa(k, ndx, ndy, ndz)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Ampsa(c, ndx, ndy, ndz)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Bmpfa(k)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Bmpsa(c)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Cmpfa(k, ndx, ndy, ndz)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Cmpsa(c, ndx, ndy, ndz)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Dmpfa(k)[source]

pumapy.physics_models.finite_volume.mpxa_matrices.fill_Dmpsa(c)[source]