Blog: Update 6/4: rad_press_blowoff.nb

File rad_press_blowoff.nb, 33.0 KB (added by adebrech, 7 years ago)
Line 
1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 11.3' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 158, 7]
13NotebookDataLength[ 33626, 640]
14NotebookOptionsPosition[ 31721, 607]
15NotebookOutlinePosition[ 32056, 622]
16CellTagsIndexPosition[ 32013, 619]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell[BoxData[{
24 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]",
25 RowBox[{
26 RowBox[{"MS", "=",
27 RowBox[{"2", "*",
28 SuperscriptBox["10", "33"]}]}], ";"}], "\[IndentingNewLine]",
29 RowBox[{
30 RowBox[{"AU", "=",
31 RowBox[{"1.5", "*",
32 SuperscriptBox["10", "13"]}]}], ";"}], "\[IndentingNewLine]",
33 RowBox[{
34 RowBox[{"Rp", "=",
35 RowBox[{"1.5", "*",
36 SuperscriptBox["10", "10"]}]}], ";"}], "\[IndentingNewLine]",
37 RowBox[{
38 RowBox[{"M", "=",
39 RowBox[{"\[Pi]", " ",
40 SuperscriptBox["R", "2"],
41 RowBox[{"(",
42 RowBox[{"\[Theta]", " ", "a0"}], ")"}], "\[Rho]"}]}],
43 ";"}], "\[IndentingNewLine]",
44 RowBox[{
45 RowBox[{"EG0", "=",
46 FractionBox[
47 RowBox[{"G", " ", "M", " ", "Mstar"}], "a0"]}],
48 ";"}], "\[IndentingNewLine]",
49 RowBox[{
50 RowBox[{"A0", "=",
51 RowBox[{"\[Theta]", " ", "a0", " ",
52 RowBox[{"(",
53 RowBox[{"2", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
54 RowBox[{
55 RowBox[{"a", "=",
56 RowBox[{"a0", "+", "\[CapitalDelta]a"}]}], ";"}], "\[IndentingNewLine]",
57 RowBox[{
58 RowBox[{"\[CapitalDelta]F", "=",
59 FractionBox[
60 SuperscriptBox["a0", "3"],
61 SuperscriptBox["a", "3"]]}], ";"}], "\[IndentingNewLine]",
62 RowBox[{
63 RowBox[{"EG", "=",
64 FractionBox[
65 RowBox[{"G", " ", "M", " ", "Mstar"}], "a"]}],
66 ";"}], "\[IndentingNewLine]",
67 RowBox[{
68 RowBox[{"A", "=",
69 RowBox[{"\[Theta]", " ", "a", " ",
70 RowBox[{"(",
71 RowBox[{"2", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
72 RowBox[{
73 RowBox[{"F", "=",
74 RowBox[{"F0", " ", "\[CapitalDelta]F"}]}], ";"}], "\[IndentingNewLine]",
75 RowBox[{
76 RowBox[{"Pphot", "=",
77 RowBox[{"F", " ", "A", " ", "e\[Gamma]"}]}],
78 ";"}], "\[IndentingNewLine]",
79 RowBox[{
80 RowBox[{"\[Tau]blowoff", "=",
81 FractionBox["EG0", "Pphot"]}], ";"}], "\[IndentingNewLine]",
82 RowBox[{
83 RowBox[{"\[Tau]move", "=",
84 FractionBox[
85 RowBox[{"EG0", "-", "EG"}], "Pphot"]}], ";"}], "\[IndentingNewLine]",
86 RowBox[{
87 RowBox[{"\[Tau]rep", "=",
88 FractionBox["M", "Mdot"]}], ";"}], "\[IndentingNewLine]",
89 RowBox[{"escaperatio", "=",
90 FractionBox["\[Tau]blowoff", "\[Tau]rep"]}], "\[IndentingNewLine]",
91 RowBox[{
92 RowBox[{"moveratio", "=",
93 RowBox[{"FullSimplify", "[",
94 FractionBox["\[Tau]move", "\[Tau]rep"], "]"}]}],
95 RowBox[{"(*",
96 RowBox[{
97 RowBox[{"Valid", " ", "if", " ", "\[CapitalDelta]a"}], " ", ">", " ",
98 RowBox[{"2", "R"}]}], "*)"}]}], "\[IndentingNewLine]",
99 RowBox[{"Mdotbal", "=",
100 RowBox[{"Extract", "[",
101 RowBox[{
102 RowBox[{"Mdot", "/.",
103 RowBox[{"Solve", "[",
104 RowBox[{
105 RowBox[{"\[Tau]move", "\[Equal]", "\[Tau]rep"}], ",", "Mdot"}],
106 "]"}]}], ",", "1"}], "]"}]}], "\[IndentingNewLine]",
107 RowBox[{"abal", "=",
108 RowBox[{"FullSimplify", "[",
109 RowBox[{
110 RowBox[{"Extract", "[",
111 RowBox[{
112 RowBox[{"\[CapitalDelta]a", "/.",
113 RowBox[{"Solve", "[",
114 RowBox[{
115 RowBox[{"\[Tau]move", "\[Equal]", "\[Tau]rep"}], ",",
116 "\[CapitalDelta]a"}], "]"}]}], ",", "2"}], "]"}], "+", "a0"}],
117 "]"}]}], "\[IndentingNewLine]",
118 RowBox[{"Fbal", "=",
119 RowBox[{"FullSimplify", "[",
120 RowBox[{"Extract", "[",
121 RowBox[{
122 RowBox[{"F0", "/.",
123 RowBox[{"Solve", "[",
124 RowBox[{
125 RowBox[{"abal", "\[Equal]", "a"}], ",", "F0"}], "]"}]}], ",", "1"}],
126 "]"}], "]"}]}]}], "Input",
127 CellChangeTimes->{{3.734084179998208*^9, 3.7340842051947517`*^9}, {
128 3.734084439082707*^9, 3.7340844524289417`*^9}, {3.73408453758842*^9,
129 3.734084540807675*^9}, {3.734084627721187*^9, 3.734084627909923*^9},
130 3.735387911518156*^9, {3.735388085892416*^9, 3.735388088724641*^9}, {
131 3.735388156723371*^9, 3.73538817212212*^9}, {3.735388281843964*^9,
132 3.735388306033917*^9}, {3.735388537396063*^9, 3.735388539412677*^9}, {
133 3.73538886634972*^9, 3.7353888697316647`*^9}, {3.7353906109295263`*^9,
134 3.735390611577392*^9}, {3.735390650991634*^9, 3.735390660828896*^9}, {
135 3.735390712085197*^9, 3.7353907529753733`*^9}, {3.735390859849416*^9,
136 3.735390870290489*^9}, {3.735390913717967*^9, 3.735390949844298*^9}, {
137 3.735390987486464*^9, 3.73539101336615*^9}, {3.735391093071527*^9,
138 3.735391170710884*^9}, {3.7353912359378033`*^9, 3.735391252122828*^9}, {
139 3.735391299029881*^9, 3.7353913341225967`*^9}, {3.735391473422173*^9,
140 3.73539147540281*^9}, {3.735393388139484*^9, 3.7353933947435904`*^9}, {
141 3.735393486762179*^9, 3.73539349641508*^9}, {3.735393589367202*^9,
142 3.7353936384912443`*^9}, {3.735393699252473*^9, 3.735393799925232*^9}, {
143 3.735394040752419*^9, 3.7353940824095993`*^9}, {3.7353942093461437`*^9,
144 3.735394228573021*^9}, {3.7353942957297487`*^9, 3.735394297950942*^9}, {
145 3.735394350227005*^9, 3.735394396028558*^9}, {3.7353956637942047`*^9,
146 3.735395690745967*^9}, 3.7353957264577627`*^9, {3.7353957708569183`*^9,
147 3.7353958289921083`*^9}, {3.7353959562795362`*^9,
148 3.7353959747314997`*^9}, {3.735396349851812*^9, 3.735396352256922*^9}, {
149 3.735396526527966*^9, 3.735396654827194*^9}, 3.735397116014626*^9,
150 3.735397156720448*^9, {3.735397273422062*^9, 3.735397293726501*^9},
151 3.7353973647239943`*^9, {3.735397465053506*^9, 3.735397533087738*^9}, {
152 3.735397564290709*^9, 3.735397616517623*^9}, {3.735397708070107*^9,
153 3.735397752160228*^9}, {3.735398235383628*^9, 3.735398238360449*^9}, {
154 3.735398407447701*^9, 3.73539840850095*^9}, {3.735398612643667*^9,
155 3.735398631710144*^9}, {3.73539867454572*^9, 3.735398710963562*^9}, {
156 3.735398760774446*^9, 3.735398771832532*^9}, {3.735398865532209*^9,
157 3.735398906018365*^9}, {3.735398958236966*^9, 3.735398995105311*^9}, {
158 3.735399281009182*^9, 3.735399284842266*^9}, {3.735399448228963*^9,
159 3.735399457224431*^9}, {3.735399491309353*^9, 3.7353994973359737`*^9}, {
160 3.735399763431449*^9, 3.7353997656870832`*^9}},
161 CellLabel->"In[1]:=",ExpressionUUID->"3bd75454-65ba-4ab8-aa6c-93d1faf35385"],
162
163Cell[BoxData[
164 FractionBox[
165 RowBox[{"G", " ", "Mdot", " ", "Mstar", " ",
166 SuperscriptBox[
167 RowBox[{"(",
168 RowBox[{"a0", "+", "\[CapitalDelta]a"}], ")"}], "2"]}],
169 RowBox[{"2", " ",
170 SuperscriptBox["a0", "4"], " ", "e\[Gamma]", " ", "F0", " ", "R", " ",
171 "\[Theta]"}]]], "Output",
172 CellChangeTimes->{
173 3.7340842058729258`*^9, {3.7340844428580017`*^9, 3.7340844527430696`*^9},
174 3.7340845411975927`*^9, 3.734084628718831*^9, 3.735387740619869*^9,
175 3.7353878048017797`*^9, 3.7353880966661587`*^9, {3.735388299710059*^9,
176 3.735388306605991*^9}, 3.735388540806287*^9, 3.735388693488846*^9,
177 3.735388907535069*^9, 3.735389371780724*^9, 3.735390616796953*^9, {
178 3.735390651631131*^9, 3.735390661354269*^9}, {3.735390727428893*^9,
179 3.735390770261911*^9}, 3.735390879011875*^9, {3.735390933344758*^9,
180 3.735390950161211*^9}, {3.735390990017276*^9, 3.735391013772992*^9}, {
181 3.735391096727951*^9, 3.735391172700852*^9}, 3.735391252680607*^9, {
182 3.735391299824855*^9, 3.735391334453298*^9}, 3.735393766415402*^9, {
183 3.735394026363492*^9, 3.735394106846942*^9}, 3.735394239893298*^9,
184 3.735394298313966*^9, {3.735394355243819*^9, 3.7353943968347483`*^9},
185 3.7353956931655684`*^9, {3.7353957856861143`*^9, 3.735395829541792*^9}, {
186 3.735395956617731*^9, 3.735395975066366*^9}, 3.7353962955414677`*^9, {
187 3.735396353110695*^9, 3.735396363629683*^9}, {3.735396550084609*^9,
188 3.735396625537466*^9}, 3.735396655791259*^9, 3.735397116378934*^9,
189 3.735397162262988*^9, 3.735397544310498*^9, {3.735397602739293*^9,
190 3.735397618129449*^9}, 3.735397728309565*^9, 3.735398240006583*^9,
191 3.7353984096910477`*^9, 3.735398541405702*^9, 3.735398711660425*^9,
192 3.7353987733349953`*^9, 3.735398870155141*^9, 3.7353989074662037`*^9,
193 3.73539895875029*^9, 3.735399003880084*^9, 3.73539928520887*^9,
194 3.735399458310413*^9, 3.735399497638406*^9, 3.735399765957768*^9,
195 3.737099300632057*^9},
196 CellLabel->"Out[17]=",ExpressionUUID->"38231f98-dde7-44db-8546-fa146c06346d"],
197
198Cell[BoxData[
199 FractionBox[
200 RowBox[{"G", " ", "Mdot", " ", "Mstar", " ", "\[CapitalDelta]a", " ",
201 RowBox[{"(",
202 RowBox[{"a0", "+", "\[CapitalDelta]a"}], ")"}]}],
203 RowBox[{"2", " ",
204 SuperscriptBox["a0", "4"], " ", "e\[Gamma]", " ", "F0", " ", "R", " ",
205 "\[Theta]"}]]], "Output",
206 CellChangeTimes->{
207 3.7340842058729258`*^9, {3.7340844428580017`*^9, 3.7340844527430696`*^9},
208 3.7340845411975927`*^9, 3.734084628718831*^9, 3.735387740619869*^9,
209 3.7353878048017797`*^9, 3.7353880966661587`*^9, {3.735388299710059*^9,
210 3.735388306605991*^9}, 3.735388540806287*^9, 3.735388693488846*^9,
211 3.735388907535069*^9, 3.735389371780724*^9, 3.735390616796953*^9, {
212 3.735390651631131*^9, 3.735390661354269*^9}, {3.735390727428893*^9,
213 3.735390770261911*^9}, 3.735390879011875*^9, {3.735390933344758*^9,
214 3.735390950161211*^9}, {3.735390990017276*^9, 3.735391013772992*^9}, {
215 3.735391096727951*^9, 3.735391172700852*^9}, 3.735391252680607*^9, {
216 3.735391299824855*^9, 3.735391334453298*^9}, 3.735393766415402*^9, {
217 3.735394026363492*^9, 3.735394106846942*^9}, 3.735394239893298*^9,
218 3.735394298313966*^9, {3.735394355243819*^9, 3.7353943968347483`*^9},
219 3.7353956931655684`*^9, {3.7353957856861143`*^9, 3.735395829541792*^9}, {
220 3.735395956617731*^9, 3.735395975066366*^9}, 3.7353962955414677`*^9, {
221 3.735396353110695*^9, 3.735396363629683*^9}, {3.735396550084609*^9,
222 3.735396625537466*^9}, 3.735396655791259*^9, 3.735397116378934*^9,
223 3.735397162262988*^9, 3.735397544310498*^9, {3.735397602739293*^9,
224 3.735397618129449*^9}, 3.735397728309565*^9, 3.735398240006583*^9,
225 3.7353984096910477`*^9, 3.735398541405702*^9, 3.735398711660425*^9,
226 3.7353987733349953`*^9, 3.735398870155141*^9, 3.7353989074662037`*^9,
227 3.73539895875029*^9, 3.735399003880084*^9, 3.73539928520887*^9,
228 3.735399458310413*^9, 3.735399497638406*^9, 3.735399765957768*^9,
229 3.737099300886467*^9},
230 CellLabel->"Out[18]=",ExpressionUUID->"eb303ce1-c470-4faf-844b-62f72bab0060"],
231
232Cell[BoxData[
233 FractionBox[
234 RowBox[{"2", " ",
235 SuperscriptBox["a0", "4"], " ", "e\[Gamma]", " ", "F0", " ", "R", " ",
236 "\[Theta]"}],
237 RowBox[{"G", " ", "Mstar", " ", "\[CapitalDelta]a", " ",
238 RowBox[{"(",
239 RowBox[{"a0", "+", "\[CapitalDelta]a"}], ")"}]}]]], "Output",
240 CellChangeTimes->{
241 3.7340842058729258`*^9, {3.7340844428580017`*^9, 3.7340844527430696`*^9},
242 3.7340845411975927`*^9, 3.734084628718831*^9, 3.735387740619869*^9,
243 3.7353878048017797`*^9, 3.7353880966661587`*^9, {3.735388299710059*^9,
244 3.735388306605991*^9}, 3.735388540806287*^9, 3.735388693488846*^9,
245 3.735388907535069*^9, 3.735389371780724*^9, 3.735390616796953*^9, {
246 3.735390651631131*^9, 3.735390661354269*^9}, {3.735390727428893*^9,
247 3.735390770261911*^9}, 3.735390879011875*^9, {3.735390933344758*^9,
248 3.735390950161211*^9}, {3.735390990017276*^9, 3.735391013772992*^9}, {
249 3.735391096727951*^9, 3.735391172700852*^9}, 3.735391252680607*^9, {
250 3.735391299824855*^9, 3.735391334453298*^9}, 3.735393766415402*^9, {
251 3.735394026363492*^9, 3.735394106846942*^9}, 3.735394239893298*^9,
252 3.735394298313966*^9, {3.735394355243819*^9, 3.7353943968347483`*^9},
253 3.7353956931655684`*^9, {3.7353957856861143`*^9, 3.735395829541792*^9}, {
254 3.735395956617731*^9, 3.735395975066366*^9}, 3.7353962955414677`*^9, {
255 3.735396353110695*^9, 3.735396363629683*^9}, {3.735396550084609*^9,
256 3.735396625537466*^9}, 3.735396655791259*^9, 3.735397116378934*^9,
257 3.735397162262988*^9, 3.735397544310498*^9, {3.735397602739293*^9,
258 3.735397618129449*^9}, 3.735397728309565*^9, 3.735398240006583*^9,
259 3.7353984096910477`*^9, 3.735398541405702*^9, 3.735398711660425*^9,
260 3.7353987733349953`*^9, 3.735398870155141*^9, 3.7353989074662037`*^9,
261 3.73539895875029*^9, 3.735399003880084*^9, 3.73539928520887*^9,
262 3.735399458310413*^9, 3.735399497638406*^9, 3.735399765957768*^9,
263 3.7370993010815268`*^9},
264 CellLabel->"Out[19]=",ExpressionUUID->"c6abdde9-868b-4004-b8cb-6f4341ccf4dd"],
265
266Cell[BoxData[
267 RowBox[{
268 FractionBox["1", "2"], " ", "a0", " ",
269 RowBox[{"(",
270 RowBox[{"1", "+",
271 FractionBox[
272 SqrtBox[
273 RowBox[{
274 RowBox[{"G", " ", "Mdot", " ", "Mstar"}], "+",
275 RowBox[{"8", " ",
276 SuperscriptBox["a0", "2"], " ", "e\[Gamma]", " ", "F0", " ", "R", " ",
277 "\[Theta]"}]}]],
278 RowBox[{
279 SqrtBox["G"], " ",
280 SqrtBox["Mdot"], " ",
281 SqrtBox["Mstar"]}]]}], ")"}]}]], "Output",
282 CellChangeTimes->{
283 3.7340842058729258`*^9, {3.7340844428580017`*^9, 3.7340844527430696`*^9},
284 3.7340845411975927`*^9, 3.734084628718831*^9, 3.735387740619869*^9,
285 3.7353878048017797`*^9, 3.7353880966661587`*^9, {3.735388299710059*^9,
286 3.735388306605991*^9}, 3.735388540806287*^9, 3.735388693488846*^9,
287 3.735388907535069*^9, 3.735389371780724*^9, 3.735390616796953*^9, {
288 3.735390651631131*^9, 3.735390661354269*^9}, {3.735390727428893*^9,
289 3.735390770261911*^9}, 3.735390879011875*^9, {3.735390933344758*^9,
290 3.735390950161211*^9}, {3.735390990017276*^9, 3.735391013772992*^9}, {
291 3.735391096727951*^9, 3.735391172700852*^9}, 3.735391252680607*^9, {
292 3.735391299824855*^9, 3.735391334453298*^9}, 3.735393766415402*^9, {
293 3.735394026363492*^9, 3.735394106846942*^9}, 3.735394239893298*^9,
294 3.735394298313966*^9, {3.735394355243819*^9, 3.7353943968347483`*^9},
295 3.7353956931655684`*^9, {3.7353957856861143`*^9, 3.735395829541792*^9}, {
296 3.735395956617731*^9, 3.735395975066366*^9}, 3.7353962955414677`*^9, {
297 3.735396353110695*^9, 3.735396363629683*^9}, {3.735396550084609*^9,
298 3.735396625537466*^9}, 3.735396655791259*^9, 3.735397116378934*^9,
299 3.735397162262988*^9, 3.735397544310498*^9, {3.735397602739293*^9,
300 3.735397618129449*^9}, 3.735397728309565*^9, 3.735398240006583*^9,
301 3.7353984096910477`*^9, 3.735398541405702*^9, 3.735398711660425*^9,
302 3.7353987733349953`*^9, 3.735398870155141*^9, 3.7353989074662037`*^9,
303 3.73539895875029*^9, 3.735399003880084*^9, 3.73539928520887*^9,
304 3.735399458310413*^9, 3.735399497638406*^9, 3.735399765957768*^9,
305 3.737099301202676*^9},
306 CellLabel->"Out[20]=",ExpressionUUID->"c9ad2ac9-c0f6-42d3-a1ea-cc2f701ef7bf"],
307
308Cell[BoxData[
309 FractionBox[
310 RowBox[{"G", " ", "Mdot", " ", "Mstar", " ", "\[CapitalDelta]a", " ",
311 RowBox[{"(",
312 RowBox[{"a0", "+", "\[CapitalDelta]a"}], ")"}]}],
313 RowBox[{"2", " ",
314 SuperscriptBox["a0", "4"], " ", "e\[Gamma]", " ", "R", " ",
315 "\[Theta]"}]]], "Output",
316 CellChangeTimes->{
317 3.7340842058729258`*^9, {3.7340844428580017`*^9, 3.7340844527430696`*^9},
318 3.7340845411975927`*^9, 3.734084628718831*^9, 3.735387740619869*^9,
319 3.7353878048017797`*^9, 3.7353880966661587`*^9, {3.735388299710059*^9,
320 3.735388306605991*^9}, 3.735388540806287*^9, 3.735388693488846*^9,
321 3.735388907535069*^9, 3.735389371780724*^9, 3.735390616796953*^9, {
322 3.735390651631131*^9, 3.735390661354269*^9}, {3.735390727428893*^9,
323 3.735390770261911*^9}, 3.735390879011875*^9, {3.735390933344758*^9,
324 3.735390950161211*^9}, {3.735390990017276*^9, 3.735391013772992*^9}, {
325 3.735391096727951*^9, 3.735391172700852*^9}, 3.735391252680607*^9, {
326 3.735391299824855*^9, 3.735391334453298*^9}, 3.735393766415402*^9, {
327 3.735394026363492*^9, 3.735394106846942*^9}, 3.735394239893298*^9,
328 3.735394298313966*^9, {3.735394355243819*^9, 3.7353943968347483`*^9},
329 3.7353956931655684`*^9, {3.7353957856861143`*^9, 3.735395829541792*^9}, {
330 3.735395956617731*^9, 3.735395975066366*^9}, 3.7353962955414677`*^9, {
331 3.735396353110695*^9, 3.735396363629683*^9}, {3.735396550084609*^9,
332 3.735396625537466*^9}, 3.735396655791259*^9, 3.735397116378934*^9,
333 3.735397162262988*^9, 3.735397544310498*^9, {3.735397602739293*^9,
334 3.735397618129449*^9}, 3.735397728309565*^9, 3.735398240006583*^9,
335 3.7353984096910477`*^9, 3.735398541405702*^9, 3.735398711660425*^9,
336 3.7353987733349953`*^9, 3.735398870155141*^9, 3.7353989074662037`*^9,
337 3.73539895875029*^9, 3.735399003880084*^9, 3.73539928520887*^9,
338 3.735399458310413*^9, 3.735399497638406*^9, 3.735399765957768*^9,
339 3.7370993013183317`*^9},
340 CellLabel->"Out[21]=",ExpressionUUID->"1e0ab6d6-35e0-437f-b6a5-f20ee492f739"]
341}, Open ]],
342
343Cell[CellGroupData[{
344
345Cell[BoxData[{
346 RowBox[{
347 RowBox[{"G", "=",
348 RowBox[{"6.67", "*",
349 SuperscriptBox["10",
350 RowBox[{"-", "8"}]]}]}], ";"}], "\[IndentingNewLine]",
351 RowBox[{
352 RowBox[{"a0", "=",
353 RowBox[{"0.04747", " ", "AU"}]}], ";"}], "\[IndentingNewLine]",
354 RowBox[{
355 RowBox[{"Mstar", "=",
356 RowBox[{"1.35", "MS"}]}], ";"}], "\[IndentingNewLine]",
357 RowBox[{
358 RowBox[{"\[Rho]", "=",
359 RowBox[{"4.11", "*",
360 SuperscriptBox["10",
361 RowBox[{"-", "19"}]]}]}], ";"}], "\[IndentingNewLine]",
362 RowBox[{
363 RowBox[{"R", "=",
364 RowBox[{"6", " ", "Rp"}]}], ";"}], "\[IndentingNewLine]",
365 RowBox[{
366 RowBox[{"e\[Gamma]", "=",
367 RowBox[{"1.63", "*",
368 SuperscriptBox["10",
369 RowBox[{"-", "11"}]]}]}], ";"}], "\[IndentingNewLine]",
370 RowBox[{
371 RowBox[{"F0", "=",
372 RowBox[{"2", "*",
373 SuperscriptBox["10", "14"]}]}], ";"}], "\[IndentingNewLine]",
374 RowBox[{
375 RowBox[{"Mdot", "=",
376 RowBox[{"3", "*",
377 SuperscriptBox["10", "10"]}]}], ";"}], "\[IndentingNewLine]",
378 RowBox[{
379 RowBox[{"\[Theta]", "=",
380 FractionBox["\[Pi]", "6"]}], ";"}], "\[IndentingNewLine]",
381 RowBox[{
382 RowBox[{"\[CapitalDelta]a", "=",
383 RowBox[{"2", "R"}]}],
384 ";"}], "\[IndentingNewLine]", "escaperatio", "\[IndentingNewLine]", \
385"moveratio", "\[IndentingNewLine]",
386 FractionBox["Mdotbal", "Mdot"], "\[IndentingNewLine]",
387 FractionBox["abal", "a0"], "\[IndentingNewLine]",
388 FractionBox["Fbal",
389 "F0"], "\[IndentingNewLine]", "\[CapitalDelta]F", "\[IndentingNewLine]",
390 RowBox[{
391 RowBox[{"\[CapitalDelta]a", "=",
392 RowBox[{
393 RowBox[{"10", "Rp"}], "+", "R"}]}],
394 ";"}], "\[IndentingNewLine]", "\[Tau]move", "\[IndentingNewLine]", \
395"\[Tau]rep", "\[IndentingNewLine]", "Fbal", "\[IndentingNewLine]",
396 RowBox[{
397 RowBox[{"\[CapitalDelta]a", "=", "Rp"}],
398 ";"}], "\[IndentingNewLine]", "Fbal"}], "Input",
399 CellChangeTimes->{{3.735388346048518*^9, 3.7353883549933167`*^9}, {
400 3.735390570710659*^9, 3.735390594989675*^9}, {3.735390626684834*^9,
401 3.735390627245739*^9}, {3.735390708546039*^9, 3.7353907210949497`*^9}, {
402 3.735390758222827*^9, 3.735390767131804*^9}, {3.735390835339801*^9,
403 3.7353908382151527`*^9}, {3.735390872692635*^9, 3.735390875114319*^9}, {
404 3.7353912113522043`*^9, 3.735391212611041*^9}, {3.735391255756678*^9,
405 3.735391256396935*^9}, 3.735391289910912*^9, {3.735391342321745*^9,
406 3.7353913668086567`*^9}, 3.735394097903352*^9, 3.735394284227829*^9, {
407 3.735394424515667*^9, 3.735394440262532*^9}, {3.735395883629739*^9,
408 3.735395920316472*^9}, {3.735396191320887*^9, 3.7353962535027237`*^9}, {
409 3.7353966617450247`*^9, 3.735396666298485*^9}, {3.735397389390411*^9,
410 3.735397391477284*^9}, 3.7353975578896646`*^9, {3.7353976252083063`*^9,
411 3.735397643727426*^9}, {3.735398307504058*^9, 3.735398318456324*^9}, {
412 3.735398360541746*^9, 3.7353983631879053`*^9}, {3.735398725542017*^9,
413 3.735398725900507*^9}, {3.735398922183124*^9, 3.7353989254815207`*^9}, {
414 3.7353990424627657`*^9, 3.73539910581122*^9}},
415 CellLabel->"In[22]:=",ExpressionUUID->"3bc784c0-8c63-4e43-a23f-eba9a96acd0d"],
416
417Cell[BoxData["0.05443253460293959`"], "Output",
418 CellChangeTimes->{
419 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
420 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
421 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
422 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
423 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
424 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
425 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
426 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
427 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
428 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
429 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
430 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
431 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
432 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
433 3.737099305595427*^9},
434 CellLabel->"Out[32]=",ExpressionUUID->"43e3ce62-c147-4ca3-b83d-64afd276589e"],
435
436Cell[BoxData["0.010983528085341771`"], "Output",
437 CellChangeTimes->{
438 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
439 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
440 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
441 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
442 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
443 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
444 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
445 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
446 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
447 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
448 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
449 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
450 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
451 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
452 3.737099305596635*^9},
453 CellLabel->"Out[33]=",ExpressionUUID->"96631df3-19c7-4e60-bc3c-cb85935f2fff"],
454
455Cell[BoxData["91.04542659062027`"], "Output",
456 CellChangeTimes->{
457 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
458 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
459 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
460 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
461 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
462 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
463 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
464 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
465 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
466 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
467 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
468 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
469 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
470 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
471 3.737099305597931*^9},
472 CellLabel->"Out[34]=",ExpressionUUID->"81f7a2b9-83bf-4bc8-ae30-de04083557e6"],
473
474Cell[BoxData["5.892921205714251`"], "Output",
475 CellChangeTimes->{
476 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
477 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
478 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
479 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
480 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
481 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
482 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
483 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
484 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
485 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
486 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
487 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
488 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
489 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
490 3.73709930559914*^9},
491 CellLabel->"Out[35]=",ExpressionUUID->"9826d640-5953-418b-9ce5-218d6a88834c"],
492
493Cell[BoxData["0.010983528085341773`"], "Output",
494 CellChangeTimes->{
495 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
496 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
497 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
498 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
499 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
500 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
501 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
502 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
503 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
504 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
505 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
506 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
507 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
508 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
509 3.7370993056006927`*^9},
510 CellLabel->"Out[36]=",ExpressionUUID->"eca4b3bd-ee9d-4ad8-a995-cb395fa7a4b8"],
511
512Cell[BoxData["0.5085853894177522`"], "Output",
513 CellChangeTimes->{
514 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
515 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
516 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
517 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
518 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
519 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
520 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
521 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
522 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
523 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
524 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
525 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
526 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
527 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
528 3.737099305602434*^9},
529 CellLabel->"Out[37]=",ExpressionUUID->"625a7bcc-3b17-4254-911b-bf304d132f7e"],
530
531Cell[BoxData["2031.4948847050584`"], "Output",
532 CellChangeTimes->{
533 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
534 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
535 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
536 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
537 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
538 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
539 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
540 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
541 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
542 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
543 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
544 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
545 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
546 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
547 3.737099305604148*^9},
548 CellLabel->"Out[39]=",ExpressionUUID->"3662aa2e-3c3a-408d-bcca-7ddeedcdf24c"],
549
550Cell[BoxData["129976.42029615105`"], "Output",
551 CellChangeTimes->{
552 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
553 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
554 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
555 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
556 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
557 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
558 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
559 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
560 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
561 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
562 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
563 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
564 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
565 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
566 3.737099305605386*^9},
567 CellLabel->"Out[40]=",ExpressionUUID->"58638298-1159-4e4d-9dae-89c71a1faace"],
568
569Cell[BoxData["3.1259437366813174`*^12"], "Output",
570 CellChangeTimes->{
571 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
572 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
573 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
574 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
575 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
576 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
577 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
578 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
579 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
580 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
581 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
582 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
583 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
584 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
585 3.737099305606627*^9},
586 CellLabel->"Out[41]=",ExpressionUUID->"7490f399-ec44-4365-b703-0d1c81883e01"],
587
588Cell[BoxData["1.4919892559175934`*^11"], "Output",
589 CellChangeTimes->{
590 3.735388355501658*^9, 3.7353885517043962`*^9, {3.735388721591456*^9,
591 3.735388741704383*^9}, 3.73538894869025*^9, 3.735389376644249*^9,
592 3.7353906798142147`*^9, {3.735390760907627*^9, 3.735390767672577*^9},
593 3.7353908396225643`*^9, 3.735390876194805*^9, 3.7353912138121853`*^9,
594 3.735391257379287*^9, {3.7353913369783697`*^9, 3.735391367430542*^9}, {
595 3.735394063352459*^9, 3.735394098436739*^9}, 3.735394284726489*^9, {
596 3.735394414584363*^9, 3.735394440778715*^9}, {3.735395881202676*^9,
597 3.735395920899523*^9}, 3.735396175667062*^9, {3.735396221761755*^9,
598 3.735396254851386*^9}, 3.735396632536746*^9, 3.7353966669772787`*^9,
599 3.735397395123358*^9, {3.73539762820903*^9, 3.735397644155486*^9},
600 3.735398319761314*^9, 3.7353983638660994`*^9, 3.735398424862804*^9,
601 3.73539854368324*^9, {3.735398720293001*^9, 3.735398727362359*^9},
602 3.7353989280160646`*^9, 3.7353989603881903`*^9, {3.735399027455608*^9,
603 3.7353990793215*^9}, 3.7353991099512177`*^9, 3.735399913497563*^9,
604 3.737099305607904*^9},
605 CellLabel->"Out[43]=",ExpressionUUID->"2c1dd49e-40f7-42a1-8159-f7414912c2a0"]
606}, Open ]]
607},
608WindowSize->{808, 911},
609WindowMargins->{{542, Automatic}, {47, Automatic}},
610FrontEndVersion->"11.3 for Linux x86 (64-bit) (March 6, 2018)",
611StyleDefinitions->"Default.nb"
612]
613(* End of Notebook Content *)
614
615(* Internal cache information *)
616(*CellTagsOutline
617CellTagsIndex->{}
618*)
619(*CellTagsIndex
620CellTagsIndex->{}
621*)
622(*NotebookFileOutline
623Notebook[{
624Cell[CellGroupData[{
625Cell[580, 22, 5909, 138, 729, "Input",ExpressionUUID->"3bd75454-65ba-4ab8-aa6c-93d1faf35385"],
626Cell[6492, 162, 2017, 33, 61, "Output",ExpressionUUID->"38231f98-dde7-44db-8546-fa146c06346d"],
627Cell[8512, 197, 2015, 32, 57, "Output",ExpressionUUID->"eb303ce1-c470-4faf-844b-62f72bab0060"],
628Cell[10530, 231, 2004, 32, 61, "Output",ExpressionUUID->"c6abdde9-868b-4004-b8cb-6f4341ccf4dd"],
629Cell[12537, 265, 2172, 40, 73, "Output",ExpressionUUID->"c9ad2ac9-c0f6-42d3-a1ea-cc2f701ef7bf"],
630Cell[14712, 307, 2006, 32, 57, "Output",ExpressionUUID->"1e0ab6d6-35e0-437f-b6a5-f20ee492f739"]
631}, Open ]],
632Cell[CellGroupData[{
633Cell[16755, 344, 3058, 70, 627, "Input",ExpressionUUID->"3bc784c0-8c63-4e43-a23f-eba9a96acd0d"],
634Cell[19816, 416, 1186, 17, 35, "Output",ExpressionUUID->"43e3ce62-c147-4ca3-b83d-64afd276589e"],
635Cell[21005, 435, 1187, 17, 35, "Output",ExpressionUUID->"96631df3-19c7-4e60-bc3c-cb85935f2fff"],
636Cell[22195, 454, 1184, 17, 35, "Output",ExpressionUUID->"81f7a2b9-83bf-4bc8-ae30-de04083557e6"],
637Cell[23382, 473, 1183, 17, 35, "Output",ExpressionUUID->"9826d640-5953-418b-9ce5-218d6a88834c"],
638Cell[24568, 492, 1189, 17, 35, "Output",ExpressionUUID->"eca4b3bd-ee9d-4ad8-a995-cb395fa7a4b8"],
639Cell[25760, 511, 1185, 17, 35, "Output",ExpressionUUID->"625a7bcc-3b17-4254-911b-bf304d132f7e"],
640Cell[26948, 530, 1185, 17, 35, "Output",ExpressionUUID->"3662aa2e-3c3a-408d-bcca-7ddeedcdf24c"],
641Cell[28136, 549, 1185, 17, 35, "Output",ExpressionUUID->"58638298-1159-4e4d-9dae-89c71a1faace"],
642Cell[29324, 568, 1189, 17, 37, "Output",ExpressionUUID->"7490f399-ec44-4365-b703-0d1c81883e01"],
643Cell[30516, 587, 1189, 17, 94, "Output",ExpressionUUID->"2c1dd49e-40f7-42a1-8159-f7414912c2a0"]
644}, Open ]]
645}
646]
647*)
648