Dear Cubit team,
I am using cubit 2025.8 in windows to set up a (preferably Hex) mesh for a complex geometry. I have hundreds of small surfaces/fractures in the doman which also cross-cut each other. For the mesh my first wish is to embed the fractures inside the volume. By embed I mean that the mesh in the fracture surfaces should share its nodes, edges, and surfaces with the volume. Secondly, I would like to have a bit of refinement on the surfaces. I very much appreciate any help. This is how I make the geometry in the cubit:
#!cubit
reset
#!python
import numpy as np
# Read Data from file
nested_fractures = np.load(r"...\DFN_Mesh\data\DFN_Frac_corners.npy")
for frac in range (len(nested_fractures)):
for point in nested_fractures[frac]:
cubit.cmd(f"create vertex location {point[0]} {point[1]} {point[2]}")
n_groups = len (nested_fractures) * 4
frac_ids = np.arange (1, n_groups+1).reshape(-1, 4)
for frac in frac_ids:
cubit.cmd(f"create surface vertex {frac[0]} {frac[1]} {frac[2]} {frac[3]}")
n_fracs = len (nested_fractures)
cubit.cmd ("brick x 100 y 50 z 520")
cubit.cmd (f"move Volume {n_fracs+1} z -260 include_merged")
#!cubit
imprint surface all
merge surface all
Here is a view into the geometry with 571 surfaces:
I have attahed my journal file and also data (a compressed .npy file).
DFN_mesh.jou (733 Bytes)
DFN_Frac_corners.zip (40.9 KB)
I very much appreciate any help.
Cheers