test(well): check well intersection feature for InternalWell andVtkmesh

Author: lucas-givord

geos-trame/tests/data/geosDeck/P4_launcher.sbatch

@@ -0,0 +1,31 @@
+#!/bin/sh
+ 
+#SBATCH --job-name="geos"
+#SBATCH --ntasks=8
+#SBATCH --nodes=1
+#SBATCH --time=3:00:00
+#SBATCH --partition=p4_general
+##SBATCH --exclusive
+#SBATCH --output=job_GEOS_P4_%j.out
+#SBATCH --error=job_GEOS_P4_%j.err
+##SBATCH --mem=734GB
+ 
+# do not change
+ulimit -s unlimited
+ulimit -c unlimited
+ 
+# loading of the module 
+module purge
+module use /workrd/SCR/NUM/GEOS_environment/p4/modulefiles/app
+module load geos/develop_3bf12d2/pangea4-gcc12.1-hpcxompi2.17.1-onemkl2023.2.0-Release-2024-06-20-10-05
+ 
+export HDF5_USE_FILE_LOCKING=FALSE
+export OMP_NUM_THREADS=1
+ 
+#----- Set standard Output and standrard error base file name. -----
+OUT_NAME=job_GEOS_P4_${SLURM_JOBID}
+ 
+srun --mpi=pmix_v3 --hint=nomultithread \
+    -n ${SLURM_NTASKS} geos \
+    -o modelA \
+    -i simulationA.xml

geos-trame/tests/data/geosDeck/co2flash.txt

@@ -0,0 +1 @@
+FlashModel CO2Solubility  1e5 1.5e8 5e4 283.15 383.15 1 0

geos-trame/tests/data/geosDeck/geosDeck.xml

@@ -98,7 +98,7 @@
                 polylineNodeCoords="{ { 800, 1450, 395.646 }, 
                                       { 800, 1450, -200.0 },
                                       { 800, 1450, 0.0 },
-                                      { 800, 1450, -554.354 } }" 
+                                      { 800, 1450, -554.354 } }"
                 polylineSegmentConn="{ { 0, 2 }, { 2, 1 }, { 1, 3 } }" 
                 radius="0.050"
                 numElementsPerSegment="5">

geos-trame/tests/data/geosDeck/pvtgas.txt

@@ -0,0 +1,2 @@
+DensityFun SpanWagnerCO2Density 1e5 1.5e8 5e4 283.15 383.15 1
+ViscosityFun FenghourCO2Viscosity 1e5 1.5e8 5e4 283.15 383.15 1

geos-trame/tests/data/geosDeck/pvtliquid.txt

@@ -0,0 +1,2 @@
+DensityFun PhillipsBrineDensity 1e5 1.5e8 5e4 283.15 383.15 1 0
+ViscosityFun PhillipsBrineViscosity 0

geos-trame/tests/data/geosDeck/well_1.vtk

@@ -0,0 +1,74 @@
+# SPDX-License-Identifier: Apache-2.0
+# SPDX-FileCopyrightText: Copyright 2023-2024 TotalEnergies.
+# SPDX-FileContributor: Lucas Givord - Kitware
+from trame.app import get_server
+from trame_client.utils.testing import enable_testing
+from geos_trame.app.core import GeosTrame
+
+
+def test_internal_well_intersection():
+
+    server = enable_testing( get_server( client_type="vue3" ), "message" )
+    file_name = "geos-trame/tests/data/geosDeck/geosDeck.xml"
+
+    app = GeosTrame( server, file_name )
+    app.state.ready()
+
+    app.deckInspector.state.object_state = [ "Problem/Mesh/0/VTKMesh/0", True ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [
+        "Problem/Mesh/0/VTKMesh/0/InternalWell/0",
+        True,
+    ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [
+        "Problem/Mesh/0/VTKMesh/0/InternalWell/0/Perforation/0",
+        True,
+    ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [
+        "Problem/Mesh/0/VTKMesh/0/InternalWell/0/Perforation/1",
+        True,
+    ]
+    app.deckInspector.state.flush()
+
+    assert app.deckViewer.well_engine.get_number_of_wells() == 1
+    assert len( app.deckViewer._perforations ) == 2
+
+
+def test_vtk_well_intersection():
+
+    server = enable_testing( get_server( client_type="vue3" ), "message" )
+    file_name = "geos-trame/tests/data/geosDeck/geosDeck.xml"
+
+    app = GeosTrame( server, file_name )
+    app.state.ready()
+
+    app.deckInspector.state.object_state = [ "Problem/Mesh/0/VTKMesh/0", True ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [ "Problem/Mesh/0/VTKMesh/0/VTKWell/0", True ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [
+        "Problem/Mesh/0/VTKMesh/0/VTKWell/0/Perforation/0",
+        True,
+    ]
+    app.deckInspector.state.flush()
+
+    app.deckInspector.state.object_state = [
+        "Problem/Mesh/0/VTKMesh/0/VTKWell/0/Perforation/1",
+        True,
+    ]
+    app.deckInspector.state.flush()
+
+    assert app.deckViewer.well_engine.get_number_of_wells() == 1
+    assert len( app.deckViewer._perforations ) == 2
+
+    app.deckInspector.state.object_state = [ "Problem/Mesh/0/VTKMesh/0/VTKWell/0", False ]
+    app.deckInspector.state.flush()
+
+    assert app.deckViewer.well_engine.get_number_of_wells() == 0