def test_mesh_biquadratic():
from pygeoiga.nurb.cad import make_surface_biquadratic
k, B = make_surface_biquadratic()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=data_path + "biquadratic.geo",
save_msh=data_path + "biquadratic.msh")
print(script)
plot_mesh(mesh)
def test_mesh_square():
from pygeoiga.nurb.cad import make_surface_square
U, V, B = make_surface_square()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=10,
save_geo=data_path + "square_surf.geo",
save_msh=data_path + "square_surf.msh")
print(script)
plot_mesh(mesh)
def test_mesh_quarter_disk():
from pygeoiga.nurb.cad import quarter_disk
k, B = quarter_disk()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=data_path + "quarter_disk.geo",
save_msh=data_path + "quarter_disk.msh")
print(script)
plot_mesh(mesh)
def test_create_script_biquadratic():
from pygeoiga.FE_solvers.run_moose import create_script
from pygeoiga.nurb.nrb_to_gmsh import convert_single_NURB_to_gmsh
from pygeoiga.nurb.cad import make_surface_biquadratic
k, B = make_surface_biquadratic()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=datapath + "biquadratic.geo",
save_msh=datapath + "biquadratic.msh")
input = datapath + "biquadratic.msh"
create_script(input,
topology_info=physical_tag_id,
bot_bc=40,
top_bc=10,
geometry=None,
kappa=4)
def test_create_script_quarter_disk():
from pygeoiga.FE_solvers.run_moose import create_script
from pygeoiga.nurb.nrb_to_gmsh import convert_single_NURB_to_gmsh
from pygeoiga.nurb.cad import quarter_disk
k, B = quarter_disk()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=datapath + "quarter_disk.geo",
save_msh=datapath + "quarter_disk.msh")
input = datapath + "quarter_disk.msh"
create_script(input,
topology_info=physical_tag_id,
right_bc=10,
left_bc=40,
geometry=None,
kappa=4)
def test_run_workflow_biquadratic():
from pygeoiga.FE_solvers.run_fenics import convert_msh_to_xdmf, run_simulation
from pygeoiga.nurb.nrb_to_gmsh import convert_single_NURB_to_gmsh
from pygeoiga.nurb.cad import make_surface_biquadratic
k, B = make_surface_biquadratic()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=datapath + "biquadratic.geo",
save_msh=datapath + "biquadratic.msh")
input = datapath + "biquadratic.msh"
plt.pause(3) # needs a bit of time to create the files
convert_msh_to_xdmf(input)
run_simulation(input,
topology_info=physical_tag_id,
bot_bc=40,
top_bc=10,
geometry=None,
kappa=4)
def test_run_workflow_square():
from pygeoiga.FE_solvers.run_fenics import convert_msh_to_xdmf, run_simulation
from pygeoiga.nurb.nrb_to_gmsh import convert_single_NURB_to_gmsh
from pygeoiga.nurb.cad import make_surface_square
U, V, B = make_surface_square()
mesh, script, physical_tag_id = convert_single_NURB_to_gmsh(
B,
size=0.5,
save_geo=datapath + "square_surf.geo",
save_msh=datapath + "square_surf.msh")
input = datapath + "square_surf.msh"
plt.pause(3)
convert_msh_to_xdmf(input)
run_simulation(input,
topology_info=physical_tag_id,
right_bc=10,
left_bc=20,
geometry=None,
kappa=4)