accretionpy 0.0.1__py3-none-any.whl

accretionpy/__init__.py

@@ -0,0 +1,2 @@
+
+

accretionpy/accretion.py

@@ -0,0 +1,65 @@
+"""this code uses user input and Eq 1 from Kummer (2025) to determine what
+ kind of BBH accretion will occur in a given system."""
+
+import numpy as np
+
+with open('initial_conditions.txt', 'r') as file:
+    variables = [int(line.strip()) for line in file if line.strip()]
+    #mstar, mbh1, mbh2, aout, eout, ain, ein = [int(i) for i in line]
+    #print(mstar, mbh1, mbh2, aout, eout, ain, ein)
+
+mstar = variables[0]
+mbh1 = variables[1]
+mbh2 = variables[2]
+aout = variables[3]
+eout = variables[4]
+ain = variables[5]
+ein = variables[6]
+
+#unit conversions
+#mstar = mstari * 1.989*10**30
+#mbh1 = mbh1i * 1.989*10**30
+#mbh2 = mbh2i * 1.989*10**30
+#aout = aouti * 6.957*10**8
+#ain = aini * 6.957*10**8
+
+qout = mstar / (mbh1 + mbh2)
+
+print(qout)
+if qout < 1:
+    rmin = 0.425 * aout * (1 - eout)*((qout)*(1+qout))
+else:
+    rmin = 0.425 * aout * (1 - eout)*((1/qout)*(1+1/qout))
+
+print(rmin)
+
+aapo = ain * (1 + ein)
+
+if aapo > rmin:
+    print(f"{aapo:.2f} is greater than {rmin:.2f}")
+    print("Ballistic accretion will occur in this system.")
+
+if rmin > aapo:
+    print(f"{aapo:.2f} is less than {rmin:.2f}")
+    print("A circumbinary disk will form in this system.")
+
+#acirc = (aout*(0.5 - 0.0986*np.log(qout))**4)*(1+qout)
+#print(acirc)
+#acav = 1.6*acirc
+#print(acav)
+##since acav / abin = 1.68, solve for abin
+#
+#abin = acav / 1.68
+#print(abin)
+#
+#if abin != rmin:
+#    print(f"binary sparations not equal: {abin:.2f} Rsun and {rmin} Rsun")
+#
+
+
+
+# eventually add plot with projected accretion path
+# use more robust calculations with modules and functions
+# calculate the critical value for BA accretion
+
+

accretionpy-0.0.1.dist-info/METADATA

@@ -0,0 +1,6 @@
+Metadata-Version: 2.4
+Name: accretionpy
+Version: 0.0.1
+License-File: LICENSE
+Requires-Dist: numpy
+Dynamic: license-file

accretionpy-0.0.1.dist-info/RECORD

@@ -0,0 +1,7 @@
+accretionpy/__init__.py,sha256=daEdpEyAJIa8b2VkCqSKcw8PaExcB6Qro80XNes_sHA,2
+accretionpy/accretion.py,sha256=qbNTGc7oa_QvhOjdITtIJ_Qc9_7EGsuNBUKoluzQDfI,1603
+accretionpy-0.0.1.dist-info/licenses/LICENSE,sha256=Qh1pDusI_asSyP76QAiMNptyRkw5dPmyN-h2egumZQ4,1074
+accretionpy-0.0.1.dist-info/METADATA,sha256=dK-DQWxl7ICRAQgJUQ56jVFhWhaq94q7k_x6KLnRLos,120
+accretionpy-0.0.1.dist-info/WHEEL,sha256=aeYiig01lYGDzBgS8HxWXOg3uV61G9ijOsup-k9o1sk,91
+accretionpy-0.0.1.dist-info/top_level.txt,sha256=dbwmkdPAwKqzhlNf_NcXD_IOV9poQmYg18njUXXv_1w,12
+accretionpy-0.0.1.dist-info/RECORD,,

accretionpy-0.0.1.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (82.0.1)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

accretionpy-0.0.1.dist-info/licenses/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) [2026] [Emma Chambers]
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

accretionpy-0.0.1.dist-info/top_level.txt

@@ -0,0 +1 @@
+accretionpy