mikaem / Oasis

55 stars 48 forks source link

Issue with pre-boundary_conditions with Oasis+VaMPy #55

Open perezrmaria opened 9 months ago

perezrmaria commented 9 months ago

Hi! We're a team trying to run a simulation on the left atrium and have encountered an error that we're not able to solve.

We've used VaMPy to preprocess our mesh, and previously tried with a mesh provided by them for testing (as we assume that mesh is okay) and after running oasis NSfracStep problem=Artery mesh_path=models/artery/artery.xml.gz save_solution_after_cycle=0 we've run into this error:

=== Mesh information === X range: 25.1176 to 42.5953 (delta: 17.4777) Y range: 24.3107 to 37.2822 (delta: 12.9715) Z range: 22.3903 to 44.4398 (delta: 22.0495) Number of cells: 158380 Number of cells per processor: 158380 Number of edges: 0 Number of faces: 320834 Number of facets: 320834 Number of vertices: 27982 Volume: 249.2360 Number of cells per volume: 635.4619 Creating initial folders === Initial pressure and area fraction === Boundary ID=2, pressure: 0.32641, area fraction: 0.33467 Boundary ID=3, pressure: 0.67359, area fraction: 0.66533 Traceback (most recent call last): File "/Users/maria/miniconda3_86/envs/vampy_env/bin/oasis", line 8, in <module> sys.exit(main()) ^^^^^^ File "/Users/maria/miniconda3_86/envs/vampy_env/lib/python3.11/site-packages/oasis/run_oasis.py", line 11, in main from oasis import NSfracStep File "/Users/maria/miniconda3_86/envs/vampy_env/lib/python3.11/site-packages/oasis/NSfracStep.py", line 172, in <module> vars().update(pre_solve_hook(**vars())) ^^^^^^^^^^^^^^^^^^^^^^^^ File "/Users/maria/miniconda3_86/envs/vampy_env/lib/python3.11/site-packages/oasis/problems/NSfracStep/Artery.py", line 177, in pre_solve_hook eval_dict["centerline_u_x_probes"] = Probes(probe_points.flatten(), V) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ TypeError: Probes() takes no arguments

We've tried with our own mesh and got the exact same error, for which we haven't found anything in the documentation. If someone has experienced this error before or has any knowledge on how to fix it we would be really happy to use some advice!

Thank you so much!