lanl / dfnWorks

dfnWorks is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. If you download the software please fill out our interest form to stay up to date on releases https://goo.gl/forms/VE39oKsyp4LVC6Gj2 and join our google group https://groups.google.com/d/forum/dfnworks-users . Precompiled Docker Container https://hub.docker.com/r/ees16/dfnworks
https://dfnworks.lanl.gov/
Other
73 stars 42 forks source link

How to define triple intersections? #49

Closed Beniuy closed 1 year ago

Beniuy commented 1 year ago

Hi, I want to define a triple intersection, but it shows that

Rejected User Defined Polygon Fracture (Defined By Coordinates) 2 rejectCode = -6: Fracture too close to another fracture's edge. Or Rejected User Defined Polygon Fracture (Defined By Coordinates) 3 rejectCode = -12: Fracture created a triple intersection with an angle too small.

How to avoid the two aforementioned error ? My polygons.dat is

nPolygons: 3
4 {-4,0,-1} {-4,0,1}  {4,0,1}  {4,0,-1} 
4 {-4,-1,-1} {-4,-1,1} {-2.9,0.1,1} {-2.9,0.1,-1}
4 {-4,1,-1} {-4,1,1} {-2.9,-0.1,1} {-2.9,-0.1,-1}

with domain size {8,10,4}. [This is modified from the example user_polygons]

Best regards

hymanjd commented 1 year ago

It appears that your fractures are not acceptable with your fram parameter h, what is that value?

Beniuy commented 1 year ago

Thank you ! I set the 'h' small enough compared with the fracture size, but it seems that the two edges can't be collinear. Here is a picture for reference. Can I make the three edges collinear? image

hymanjd commented 1 year ago

I"m guessing that your top edge is the issue. Try running with disableFram : True and visualmode to look at it closer in paraview

Beniuy commented 1 year ago

Thanks for your advice! I set the disableFram true, and the multiple fractures are not rejected now. But it says cannot find uge file ? image Tips : My purpose is construct a model like this, it has triple intersections. image

My polygon file is as follow

nPolygons: 5 4 {-1.0,0,-1} {-1.0,0,1} {+1.0,0,1} {+1.0,0,-1} 4 {-2,-1,-1} {-2,-1,1} {-1,0.0,1} {-1,0.0,-1} 4 {-2,1,-1} {-2,1,1} {-1.0,-0.0,1} {-1.0,-0.0,-1} 4 {+2,-1,-1} {+2,-1,1} {+1, 0,1} {+1,0,-1} 4 {+2,1,-1} {+2,1,1} {+1.,-0.0,1} {+1.,-0.0,-1}

hymanjd commented 1 year ago

If you run with disable fram, it will not create a uge file, and you can’t run flow. This is just for visualization.

I’d suggest extending your polygon edges a little past the intersections. This might fix your problem.

[signature_61897647]

Jeffrey D. Hyman Staff Scientist Energy and Natural Resources Security (EES-16) Los Alamos National Laboratory https://dfnworks.lanl.gov/ Pronouns: he, him, his Office: +1 (505) 665-8737

From: Xu Yue @.> Reply-To: lanl/dfnWorks @.> Date: Wednesday, February 8, 2023 at 8:54 PM To: lanl/dfnWorks @.> Cc: Jeffrey Hyman @.>, Comment @.***> Subject: [EXTERNAL] Re: [lanl/dfnWorks] How to define triple intersections? (Issue #49)

Thanks for your advice! I set the disableFram true, and the multiple fractures are not rejected now. But it says cannot find uge file ? [Image removed by sender. image]https://urldefense.com/v3/__https:/user-images.githubusercontent.com/50405929/217713291-17187393-5687-4be7-84f9-3aad7ce42012.png__;!!Bt8fGhp8LhKGRg!D6_qcnLnPcuq5I-kSv1NOuv-RhX2VFivv2CEq_ZbulUHTvy0STr-AhAyO_uWmhizpmCfukNb4t-_oispq5-vCzpw$ Tips : My purpose is construct a model like this, it has triple intersections. [Image removed by sender. image]https://urldefense.com/v3/__https:/user-images.githubusercontent.com/50405929/217713854-f38590aa-269b-40a0-8123-282fa4ca87fb.png__;!!Bt8fGhp8LhKGRg!D6_qcnLnPcuq5I-kSv1NOuv-RhX2VFivv2CEq_ZbulUHTvy0STr-AhAyO_uWmhizpmCfukNb4t-_oispq0g3NOnt$

My polygon file is as follow

nPolygons: 5 4 {-1.0,0,-1} {-1.0,0,1} {+1.0,0,1} {+1.0,0,-1} 4 {-2,-1,-1} {-2,-1,1} {-1,0.0,1} {-1,0.0,-1} 4 {-2,1,-1} {-2,1,1} {-1.0,-0.0,1} {-1.0,-0.0,-1} 4 {+2,-1,-1} {+2,-1,1} {+1, 0,1} {+1,0,-1} 4 {+2,1,-1} {+2,1,1} {+1.,-0.0,1} {+1.,-0.0,-1}

— Reply to this email directly, view it on GitHubhttps://urldefense.com/v3/__https:/github.com/lanl/dfnWorks/issues/49*issuecomment-1423588565__;Iw!!Bt8fGhp8LhKGRg!D6_qcnLnPcuq5I-kSv1NOuv-RhX2VFivv2CEq_ZbulUHTvy0STr-AhAyO_uWmhizpmCfukNb4t-_oispqwrWVk3W$, or unsubscribehttps://urldefense.com/v3/__https:/github.com/notifications/unsubscribe-auth/AEKOGDEPTR4WKK47PEP5KLTWWRS7JANCNFSM6AAAAAAUS3MP3M__;!!Bt8fGhp8LhKGRg!D6_qcnLnPcuq5I-kSv1NOuv-RhX2VFivv2CEq_ZbulUHTvy0STr-AhAyO_uWmhizpmCfukNb4t-_oispq0jbc9Bi$. You are receiving this because you commented.Message ID: @.***>

Beniuy commented 1 year ago

Thank you very much ! I am really grateful for your prompt advice !