Questions on well model

[Hammond, Glenn E]

Michael/Jenn,

Is there a reason why circle_center(2) is not referenced instead of cur_location(2) in 

temp_trajectory(i,2) = cur_location(2)....

at lines:

https://bitbucket.org/pflotran/pflotran/src/2c2909ed88b2262b360e5872dbc685cd22d4a2d7/src/pflotran/pm_well.F90#lines-1224

https://bitbucket.org/pflotran/pflotran/src/2c2909ed88b2262b360e5872dbc685cd22d4a2d7/src/pflotran/pm_well.F90#lines-1262

https://bitbucket.org/pflotran/pflotran/src/2c2909ed88b2262b360e5872dbc685cd22d4a2d7/src/pflotran/pm_well.F90#lines-1309

 I understand that circle_center(2) and cur_location(2) are identical at line 1224, but they may not be for 1262 and 1309.

Thanks,

Glenn

[Frederick, Jennifer Mary]

Glenn, 

This was a portion of the code that Michael wrote. I defer to him to answer this one.

Jenn

 

Jennifer M. Frederick, Ph.D. Computational Geoscientist & Engineer Principal Member of the Technical Staff 08931  |  Earth Systems Analysis Center for Energy and Earth Systems Sandia National Laboratories   Mobile 505.313.3233 Email jennifer....@sandia.gov Web https://www.sandia.gov/-jmfrede Exceptional service in the national interest.

---

From: 'Hammond, Glenn E' via pflotran-dev <pflotr...@googlegroups.com>
Sent: Thursday, July 10, 2025 11:25:34 AM
To: pflotran-dev <pflotr...@googlegroups.com>
Cc: Bayat, Elyce Jane <eba...@sandia.gov>; de Toledo Camargo, Julia <julia....@pnnl.gov>
Subject: [EXTERNAL] [pflotran-dev: 6341] Questions on well model

 

--
You received this message because you are subscribed to the Google Groups "pflotran-dev" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-dev...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/pflotran-dev/PH0PR09MB7836A86334C697EBB1EB04319A48A%40PH0PR09MB7836.namprd09.prod.outlook.com.

[Michael Nole]

This seems to me like it might be a bug.

At line 1224, your well is kicking off to horizontal in the x-direction. When you specify a kickoff radius, the code draws a well along the edge of a circle with that radius until the well points in the x-direction. circle_center is the center of that circle. The reason temp_trajectory(i,2) uses cur_location at line 1224 is because the well trajectory is not changing in the y-direction.

At line 1262, you’re kicking off to horizontal in the y-direction and the trajectory should not be changing in the x-direction. I would think that in this case, temp_trajectory(i,1) should use cur_location(1).

And since at line 1309 you’re moving in along the circle in x,y, and z temp_trajectory(i,2) should use circle_center(2).

Michael

To view this discussion visit https://groups.google.com/d/msgid/pflotran-dev/SJ0PR09MB61124EB113E5283A21563172CD55A%40SJ0PR09MB6112.namprd09.prod.outlook.com.