Skip to content
Snippets Groups Projects
Commit 37095c13 authored by Kenan Cengiz's avatar Kenan Cengiz
Browse files

Update Liutex.py

parent 5826df57
No related tags found
1 merge request!1output vtu instead
Pipeline #18147 failed
......@@ -200,8 +200,8 @@ def main():
mesh = update_mesh_with_omega_liutex(mesh, a, b)
# Save the mesh with vorticity and omega data to a VTK file
mesh.save("Liu_Rortex_T106C_Re80K_19.09.vtk")
print("Dataset saved as Liu_Rortex_T106C_Re80K_19.09.vtk")
mesh.save("Liu_Rortex_T106C_Re80K_19.09.vtu")
print("Dataset saved as Liu_Rortex_T106C_Re80K_19.09.vtu")
# Selecting seed points in high vorticity regions
vorticity_threshold = np.percentile(vorticity_magnitude, 90) # Top 10% vorticity magnitude
......@@ -224,8 +224,8 @@ def main():
# Save vorticity lines
vorticity_mesh = mesh.streamlines_from_source(seed_polydata, vectors='Velocity')
vorticity_mesh.save("Vorticity_Lines_Blade_19.09.vtk")
print("Vorticity lines saved as Vorticity_Lines_Blade_19.09.vtk")
vorticity_mesh.save("Vorticity_Lines_Blade_19.09.vtu")
print("Vorticity lines saved as Vorticity_Lines_Blade_19.09.vtu")
plotter.add_axes()
plotter.view_isometric()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment