trunk/nonperiodicVSS/friction.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     14487,        426]
NotebookOptionsPosition[     12885,        369]
NotebookOutlinePosition[     13239,        385]
CellTagsIndexPosition[     13196,        382]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"a", " ", "=", " ", "33."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b", " ", "=", " ", "13."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"r", " ", "=", " ", "40."}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"eta", " ", "=", " ", "1."}], ";"}], "\n", 
 RowBox[{
  RowBox[{"ab", " ", "=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"a", "^", "2"}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{"b", "^", "2"}], ")"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"v", " ", "=", " ", 
   RowBox[{"32", "*", "Pi", "*", 
    RowBox[{"eta", "/", "3"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.584807426161751*^9, 3.584807453166168*^9}, {
   3.5848098994271517`*^9, 3.58480990098132*^9}, {3.584890669262165*^9, 
   3.5848906732045937`*^9}, 3.5848942046023273`*^9, 3.5848942364960403`*^9, {
   3.584894394914665*^9, 3.584894407017727*^9}, {3.584894867879013*^9, 
   3.5848948703643208`*^9}, {3.584895906363385*^9, 3.584895906913582*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"RESISTANCE", " ", "TENSORS"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.5848081048411083`*^9, 3.584808114547517*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"sphere", " ", "rotation"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5848955452222424`*^9, 3.584895552580024*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"XiS", " ", "=", " ", 
  RowBox[{"8.", "*", "Pi", "*", "eta", "*", 
   RowBox[{"(", 
    RowBox[{"r", "^", "3"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.584895562690569*^9, 3.584895581538089*^9}}],

Cell[BoxData["1.608495438637974`*^6"], "Output",
 CellChangeTimes->{3.5848955833332977`*^9, 3.584895910460045*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"S", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"2", "/", "ab"}], ")"}], "*", 
    RowBox[{"Log", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"a", "+", "ab"}], ")"}], "/", "b"}], "]"}]}]}], ";"}]], "Input",\

 CellChangeTimes->{{3.584895591890572*^9, 3.584895628889739*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"ellipsoid", " ", "axial", " ", 
   RowBox[{"rotation", " ", "--"}], " ", "about", " ", "long", " ", "axis"}], 
  " ", "*)"}]], "Input",
 CellChangeTimes->{{3.584807845501628*^9, 3.584807881963066*^9}, 
   3.58480810065349*^9, {3.584810005255805*^9, 3.584810006597301*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"nXiA", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"ab", "^", "2"}], ")"}], "*", 
    RowBox[{"(", 
     RowBox[{"b", "^", "2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dXiA", " ", "=", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", "*", "a"}], " ", "-", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"b", "^", "2"}], ")"}], "*", "S"}]}], ")"}]}], ";"}]}], "Input",\

 CellChangeTimes->{{3.584894479745336*^9, 3.584894547193737*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"XiA", " ", "=", " ", 
  RowBox[{"v", " ", "*", 
   RowBox[{"(", " ", 
    RowBox[{"nXiA", "/", "dXiA"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.584807580926756*^9, 3.5848076774271517`*^9}, {
  3.584893446441492*^9, 3.584893456926993*^9}, {3.584893513019639*^9, 
  3.5848935374433117`*^9}, {3.5848937295629377`*^9, 3.5848937320361967`*^9}, {
  3.5848940598863373`*^9, 3.5848940827349033`*^9}, {3.584894124697816*^9, 
  3.584894131830501*^9}, {3.5848944118465767`*^9, 3.584894413689575*^9}, {
  3.584894553484161*^9, 3.5848945620889597`*^9}}],

Cell[BoxData["107748.7123755539`"], "Output",
 CellChangeTimes->{
  3.58480767838589*^9, 3.584812369332225*^9, 3.58489067890679*^9, 
   3.584893466947144*^9, {3.584893515809536*^9, 3.5848935380765*^9}, 
   3.584893667613534*^9, 3.584893745846806*^9, 3.584893991496749*^9, {
   3.5848940686270847`*^9, 3.584894083689426*^9}, 3.584894132885084*^9, 
   3.584894209263301*^9, 3.5848944220659027`*^9, 3.584894562584774*^9, 
   3.5848956426941223`*^9, 3.584895912816411*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"ellipsoid", " ", "equatorial", " ", 
    RowBox[{"rotation", " ", "--"}], " ", "about", " ", "short", " ", 
    "axes"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.5848078641082973`*^9, 3.584807887971993*^9}, {
  3.584810008477813*^9, 3.58481001069223*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"nXiB", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"a", "^", "4"}], ")"}], "-", 
    RowBox[{"(", 
     RowBox[{"b", "^", "4"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dXiB", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "*", 
        RowBox[{"(", 
         RowBox[{"a", "^", "2"}], ")"}]}], "-", 
       RowBox[{"(", 
        RowBox[{"b", "^", "2"}], ")"}]}], ")"}], "*", "S"}], " ", "-", " ", 
    RowBox[{"(", 
     RowBox[{"2", "*", "a"}], ")"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.584894566445219*^9, 3.5848945876718807`*^9}, {
  3.5848946888115788`*^9, 3.5848947150455017`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"XiB", " ", "=", " ", 
  RowBox[{"v", " ", "*", " ", 
   RowBox[{"(", 
    RowBox[{"nXiB", "/", "dXiB"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.584807667200397*^9, 3.5848076699441557`*^9}, {
   3.5848077499771643`*^9, 3.584807786900511*^9}, {3.584893547107358*^9, 
   3.584893579226317*^9}, {3.58489389632187*^9, 3.584893901176014*^9}, {
   3.584893943459343*^9, 3.584893946518909*^9}, {3.584894088000602*^9, 
   3.5848941203441153`*^9}, 3.5848944172115498`*^9, {3.584894722257359*^9, 
   3.584894733911747*^9}}],

Cell[BoxData["269785.72042901954`"], "Output",
 CellChangeTimes->{
  3.584807788832862*^9, 3.584812370076736*^9, 3.5848906801945047`*^9, 
   3.5848935801190357`*^9, 3.5848936688170233`*^9, 3.5848939033674097`*^9, {
   3.58489396555275*^9, 3.584893992584134*^9}, 3.5848940700482283`*^9, {
   3.584894103054511*^9, 3.584894120965992*^9}, 3.584894210383068*^9, 
   3.584894423252715*^9, 3.584894734499061*^9, 3.58489564491859*^9, 
   3.584895914221643*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.584887428382176*^9, 3.584891660811152*^9}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"FRICTION", " ", "FACTORS", " ", "via", " ", "Wikipedia"}], " ", 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.584887434985119*^9, 3.58488743942144*^9}, {
  3.5848907429833384`*^9, 3.584890746521543*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"xi", " ", "=", 
   RowBox[{
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"p", "^", "2"}], ")"}], "-", "1"}], "]"}], "/", "p"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"s", " ", "=", " ", 
   RowBox[{"2", "*", 
    RowBox[{
     RowBox[{"ArcTanh", "[", "xi", "]"}], "/", "xi"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p", " ", "=", " ", 
   RowBox[{"a", "/", "b"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"re", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"a", "*", 
      RowBox[{"(", 
       RowBox[{"b", "^", "2"}], ")"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"1", "/", "3"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"eqSphere", " ", "=", " ", 
   RowBox[{"8", "*", "Pi", "*", "eta", "*", 
    RowBox[{"(", 
     RowBox[{"re", "^", "3"}], ")"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.584887453255073*^9, 3.584887474611525*^9}, {
  3.584887518352754*^9, 3.5848875201183167`*^9}, {3.5848875705654163`*^9, 
  3.584887612684559*^9}, {3.584891207938551*^9, 3.5848912105877657`*^9}, {
  3.58489477575947*^9, 3.584894789206271*^9}, {3.5848948909071417`*^9, 
  3.5848949045983353`*^9}, {3.584895894748207*^9, 3.58489594161933*^9}, {
  3.584895989916945*^9, 3.584896035299241*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.5848876347988167`*^9, 3.584887638277546*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"ellipsoid", " ", "axial", " ", "friction", " ", 
   RowBox[{"factor", " ", "--"}], " ", "about", " ", "long", " ", "axis"}], 
  " ", "*)"}]], "Input",
 CellChangeTimes->{{3.58488764265452*^9, 3.584887652272901*^9}, {
  3.58489128903328*^9, 3.5848912914194727`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Fax", " ", "=", " ", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"4", "/", "3"}], ")"}], "*", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"xi", "^", "2"}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"2", "-", 
      RowBox[{"(", 
       RowBox[{"s", "/", 
        RowBox[{"(", 
         RowBox[{"p", "^", "2"}], ")"}]}], ")"}]}], ")"}]}]}]}]], "Input",
 CellChangeTimes->{{3.58488765759788*^9, 3.584887700018955*^9}, {
   3.5848879174727783`*^9, 3.584887938981941*^9}, {3.584890603482909*^9, 
   3.5848906063233852`*^9}, 3.584890946638309*^9, {3.584894926277274*^9, 
   3.584894983989938*^9}}],

Cell[BoxData["0.7687260259259584`"], "Output",
 CellChangeTimes->{
  3.58488770586084*^9, {3.58488792206218*^9, 3.584887939492504*^9}, 
   3.584890608261299*^9, 3.5848906924156*^9, 3.584890948261269*^9, 
   3.584891214003551*^9, 3.584893997414321*^9, 3.584894794095375*^9, 
   3.5848948750854387`*^9, 3.584894906946406*^9, {3.584894984628289*^9, 
   3.5848950094024963`*^9}, 3.584895649579417*^9, 3.5848960393110647`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{"(*", " ", 
   RowBox[{"ellipsoid", " ", "equatorial", " ", "friction", " ", 
    RowBox[{"factor", " ", "--"}], " ", "about", " ", "short", " ", "axes"}], 
   " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.584887710741127*^9, 3.584887713271037*^9}, 
   3.584887813383315*^9, {3.584891294243541*^9, 3.5848912960435057`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Feq", " ", "=", " ", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"4", "/", "3"}], ")"}], " ", "*", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "/", "p"}], ")"}], "^", "2"}], "  ", "-", " ", 
      RowBox[{"p", "^", "2"}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"2", " ", "-", " ", 
      RowBox[{"s", "*", 
       RowBox[{"(", 
        RowBox[{"2", " ", "-", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "/", "p"}], ")"}], "^", "2"}]}], ")"}]}]}], 
     ")"}]}]}]}]], "Input",
 CellChangeTimes->{{3.584887727157157*^9, 3.584887769116691*^9}, {
   3.584887926989279*^9, 3.584887942885796*^9}, {3.584890609957733*^9, 
   3.584890611532776*^9}, 3.584890951791182*^9}],

Cell[BoxData["1.9247682885909356`"], "Output",
 CellChangeTimes->{
  3.5848877700852623`*^9, {3.584887929629698*^9, 3.5848879434378977`*^9}, 
   3.584890612056086*^9, 3.5848906935323477`*^9, 3.584890952137286*^9, 
   3.584891216311064*^9, 3.584893998534431*^9, 3.584894876538658*^9, 
   3.5848950109231567`*^9, 3.584895650617403*^9, 3.584896040714982*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"ellipsoid", " ", "axial", " ", "rotational", " ", "friction"}], 
   " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.58489104521352*^9, 3.584891053491618*^9}, {
  3.584891299813236*^9, 3.5848913004309483`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fAx", " ", "=", " ", 
  RowBox[{"Fax", " ", "*", " ", "eqSphere"}]}]], "Input",
 CellChangeTimes->{{3.584891056631095*^9, 3.584891064865149*^9}}],

Cell[BoxData["107748.71237555385`"], "Output",
 CellChangeTimes->{3.5848910656263733`*^9, 3.584891217398337*^9, 
  3.58489400012298*^9, 3.584895019483296*^9, 3.584895652455336*^9, 
  3.584896042503642*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "ellipsoid", " ", "equatorial", " ", "rotational", " ", "friction"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.5848910141246853`*^9, 3.584891022193733*^9}, {
  3.5848913017358427`*^9, 3.5848913022309103`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fEq", " ", "=", " ", 
  RowBox[{"Feq", " ", "*", " ", "eqSphere"}]}]], "Input",
 CellChangeTimes->{{3.584891025068652*^9, 3.584891029821259*^9}}],

Cell[BoxData["269785.7204290194`"], "Output",
 CellChangeTimes->{3.584891030572219*^9, 3.584891218233328*^9, 
  3.584894001208784*^9, 3.584895020435354*^9, 3.584895653375164*^9, 
  3.58489604404213*^9}]
}, Open  ]]
},
WindowSize->{740, 876},
WindowMargins->{{Automatic, 190}, {Automatic, 15}},
FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 22, \
2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1053, 25, 103, "Input"],
Cell[1613, 47, 204, 4, 43, "Input"],
Cell[1820, 53, 164, 3, 27, "Input"],
Cell[CellGroupData[{
Cell[2009, 60, 223, 5, 27, "Input"],
Cell[2235, 67, 114, 1, 30, "Output"]
}, Open  ]],
Cell[2364, 71, 335, 11, 27, "Input"],
Cell[2702, 84, 318, 6, 27, "Input"],
Cell[3023, 92, 524, 17, 43, "Input"],
Cell[CellGroupData[{
Cell[3572, 113, 568, 10, 27, "Input"],
Cell[4143, 125, 468, 7, 27, "Output"]
}, Open  ]],
Cell[4626, 135, 343, 7, 43, "Input"],
Cell[4972, 144, 717, 22, 43, "Input"],
Cell[CellGroupData[{
Cell[5714, 170, 540, 10, 27, "Input"],
Cell[6257, 182, 453, 7, 27, "Output"]
}, Open  ]],
Cell[6725, 192, 109, 1, 43, "Input"],
Cell[6837, 195, 284, 6, 43, "Input"],
Cell[7124, 203, 1342, 37, 88, "Input"],
Cell[8469, 242, 94, 1, 27, "Input"],
Cell[8566, 245, 311, 6, 27, "Input"],
Cell[CellGroupData[{
Cell[8902, 255, 615, 17, 27, "Input"],
Cell[9520, 274, 421, 6, 27, "Output"]
}, Open  ]],
Cell[9956, 283, 364, 7, 43, "Input"],
Cell[CellGroupData[{
Cell[10345, 294, 782, 23, 27, "Input"],
Cell[11130, 319, 355, 5, 27, "Output"]
}, Open  ]],
Cell[11500, 327, 289, 6, 43, "Input"],
Cell[CellGroupData[{
Cell[11814, 337, 169, 3, 27, "Input"],
Cell[11986, 342, 205, 3, 27, "Output"]
}, Open  ]],
Cell[12206, 348, 264, 6, 27, "Input"],
Cell[CellGroupData[{
Cell[12495, 358, 169, 3, 27, "Input"],
Cell[12667, 363, 202, 3, 27, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
Revision:	3933
Committed:	Thu Aug 8 16:15:12 2013 UTC (12 years, 1 month ago) by kstocke1
File size:	14685 byte(s)
Log Message:
#	Content
1	(* Content-type: application/vnd.wolfram.mathematica *)
2
3	(* Wolfram Notebook File *)
4	(* http://www.wolfram.com/nb *)
5
6	(* CreatedBy='Mathematica 8.0' *)
7
8	(CacheID: 234)
9	(* Internal cache information:
10	NotebookFileLineBreakTest
11	NotebookFileLineBreakTest
12	NotebookDataPosition[ 157, 7]
13	NotebookDataLength[ 14487, 426]
14	NotebookOptionsPosition[ 12885, 369]
15	NotebookOutlinePosition[ 13239, 385]
16	CellTagsIndexPosition[ 13196, 382]
17	WindowFrame->Normal*)
18
19	(* Beginning of Notebook Content *)
20	Notebook[{
21	Cell[BoxData[{
22	RowBox[{
23	RowBox[{"a", " ", "=", " ", "33."}], ";"}], "\[IndentingNewLine]",
24	RowBox[{
25	RowBox[{"b", " ", "=", " ", "13."}], ";"}], "\[IndentingNewLine]",
26	RowBox[{
27	RowBox[{"r", " ", "=", " ", "40."}], ";"}], "\[IndentingNewLine]",
28	RowBox[{
29	RowBox[{"eta", " ", "=", " ", "1."}], ";"}], "\n",
30	RowBox[{
31	RowBox[{"ab", " ", "=", " ",
32	RowBox[{"Sqrt", "[",
33	RowBox[{
34	RowBox[{"(",
35	RowBox[{"a", "^", "2"}], ")"}], "-",
36	RowBox[{"(",
37	RowBox[{"b", "^", "2"}], ")"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
38	RowBox[{
39	RowBox[{"v", " ", "=", " ",
40	RowBox[{"32", "", "Pi", "",
41	RowBox[{"eta", "/", "3"}]}]}], ";"}]}], "Input",
42	CellChangeTimes->{{3.584807426161751^9, 3.584807453166168^9}, {
43	3.5848098994271517`^9, 3.58480990098132^9}, {3.584890669262165*^9,
44	3.5848906732045937`^9}, 3.5848942046023273`^9, 3.5848942364960403`*^9, {
45	3.584894394914665^9, 3.584894407017727^9}, {3.584894867879013*^9,
46	3.5848948703643208`^9}, {3.584895906363385^9, 3.584895906913582*^9}}],
47
48	Cell[BoxData[
49	RowBox[{"\[IndentingNewLine]",
50	RowBox[{"(*", " ",
51	RowBox[{"RESISTANCE", " ", "TENSORS"}], " ", "*)"}]}]], "Input",
52	CellChangeTimes->{{3.5848081048411083`^9, 3.584808114547517^9}}],
53
54	Cell[BoxData[
55	RowBox[{"(*", " ",
56	RowBox[{"sphere", " ", "rotation"}], " ", "*)"}]], "Input",
57	CellChangeTimes->{{3.5848955452222424`^9, 3.584895552580024^9}}],
58
59	Cell[CellGroupData[{
60
61	Cell[BoxData[
62	RowBox[{"XiS", " ", "=", " ",
63	RowBox[{"8.", "", "Pi", "", "eta", "*",
64	RowBox[{"(",
65	RowBox[{"r", "^", "3"}], ")"}]}]}]], "Input",
66	CellChangeTimes->{{3.584895562690569^9, 3.584895581538089^9}}],
67
68	Cell[BoxData["1.608495438637974`*^6"], "Output",
69	CellChangeTimes->{3.5848955833332977`^9, 3.584895910460045^9}]
70	}, Open ]],
71
72	Cell[BoxData[
73	RowBox[{
74	RowBox[{"S", " ", "=", " ",
75	RowBox[{
76	RowBox[{"(",
77	RowBox[{"2", "/", "ab"}], ")"}], "*",
78	RowBox[{"Log", "[",
79	RowBox[{
80	RowBox[{"(",
81	RowBox[{"a", "+", "ab"}], ")"}], "/", "b"}], "]"}]}]}], ";"}]], "Input",\
82
83	CellChangeTimes->{{3.584895591890572^9, 3.584895628889739^9}}],
84
85	Cell[BoxData[
86	RowBox[{"(*", " ",
87	RowBox[{"ellipsoid", " ", "axial", " ",
88	RowBox[{"rotation", " ", "--"}], " ", "about", " ", "long", " ", "axis"}],
89	" ", "*)"}]], "Input",
90	CellChangeTimes->{{3.584807845501628^9, 3.584807881963066^9},
91	3.58480810065349^9, {3.584810005255805^9, 3.584810006597301*^9}}],
92
93	Cell[BoxData[{
94	RowBox[{
95	RowBox[{"nXiA", " ", "=", " ",
96	RowBox[{
97	RowBox[{"(",
98	RowBox[{"ab", "^", "2"}], ")"}], "*",
99	RowBox[{"(",
100	RowBox[{"b", "^", "2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
101	RowBox[{
102	RowBox[{"dXiA", " ", "=", " ",
103	RowBox[{"(",
104	RowBox[{
105	RowBox[{"2", "*", "a"}], " ", "-", " ",
106	RowBox[{
107	RowBox[{"(",
108	RowBox[{"b", "^", "2"}], ")"}], "*", "S"}]}], ")"}]}], ";"}]}], "Input",\
109
110	CellChangeTimes->{{3.584894479745336^9, 3.584894547193737^9}}],
111
112	Cell[CellGroupData[{
113
114	Cell[BoxData[
115	RowBox[{"XiA", " ", "=", " ",
116	RowBox[{"v", " ", "*",
117	RowBox[{"(", " ",
118	RowBox[{"nXiA", "/", "dXiA"}], ")"}]}]}]], "Input",
119	CellChangeTimes->{{3.584807580926756^9, 3.5848076774271517`^9}, {
120	3.584893446441492^9, 3.584893456926993^9}, {3.584893513019639*^9,
121	3.5848935374433117`^9}, {3.5848937295629377`^9, 3.5848937320361967`*^9}, {
122	3.5848940598863373`^9, 3.5848940827349033`^9}, {3.584894124697816*^9,
123	3.584894131830501^9}, {3.5848944118465767`^9, 3.584894413689575*^9}, {
124	3.584894553484161^9, 3.5848945620889597`^9}}],
125
126	Cell[BoxData["107748.7123755539`"], "Output",
127	CellChangeTimes->{
128	3.58480767838589^9, 3.584812369332225^9, 3.58489067890679*^9,
129	3.584893466947144^9, {3.584893515809536^9, 3.5848935380765*^9},
130	3.584893667613534^9, 3.584893745846806^9, 3.584893991496749*^9, {
131	3.5848940686270847`^9, 3.584894083689426^9}, 3.584894132885084*^9,
132	3.584894209263301^9, 3.5848944220659027`^9, 3.584894562584774*^9,
133	3.5848956426941223`^9, 3.584895912816411^9}]
134	}, Open ]],
135
136	Cell[BoxData[
137	RowBox[{"\[IndentingNewLine]",
138	RowBox[{"(*", " ",
139	RowBox[{"ellipsoid", " ", "equatorial", " ",
140	RowBox[{"rotation", " ", "--"}], " ", "about", " ", "short", " ",
141	"axes"}], " ", "*)"}]}]], "Input",
142	CellChangeTimes->{{3.5848078641082973`^9, 3.584807887971993^9}, {
143	3.584810008477813^9, 3.58481001069223^9}}],
144
145	Cell[BoxData[{
146	RowBox[{
147	RowBox[{"nXiB", " ", "=", " ",
148	RowBox[{
149	RowBox[{"(",
150	RowBox[{"a", "^", "4"}], ")"}], "-",
151	RowBox[{"(",
152	RowBox[{"b", "^", "4"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
153	RowBox[{
154	RowBox[{"dXiB", " ", "=", " ",
155	RowBox[{
156	RowBox[{
157	RowBox[{"(",
158	RowBox[{
159	RowBox[{"2", "*",
160	RowBox[{"(",
161	RowBox[{"a", "^", "2"}], ")"}]}], "-",
162	RowBox[{"(",
163	RowBox[{"b", "^", "2"}], ")"}]}], ")"}], "*", "S"}], " ", "-", " ",
164	RowBox[{"(",
165	RowBox[{"2", "*", "a"}], ")"}]}]}], ";"}]}], "Input",
166	CellChangeTimes->{{3.584894566445219^9, 3.5848945876718807`^9}, {
167	3.5848946888115788`^9, 3.5848947150455017`^9}}],
168
169	Cell[CellGroupData[{
170
171	Cell[BoxData[
172	RowBox[{"XiB", " ", "=", " ",
173	RowBox[{"v", " ", "*", " ",
174	RowBox[{"(",
175	RowBox[{"nXiB", "/", "dXiB"}], ")"}]}]}]], "Input",
176	CellChangeTimes->{{3.584807667200397^9, 3.5848076699441557`^9}, {
177	3.5848077499771643`^9, 3.584807786900511^9}, {3.584893547107358*^9,
178	3.584893579226317^9}, {3.58489389632187^9, 3.584893901176014*^9}, {
179	3.584893943459343^9, 3.584893946518909^9}, {3.584894088000602*^9,
180	3.5848941203441153`^9}, 3.5848944172115498`^9, {3.584894722257359*^9,
181	3.584894733911747*^9}}],
182
183	Cell[BoxData["269785.72042901954`"], "Output",
184	CellChangeTimes->{
185	3.584807788832862^9, 3.584812370076736^9, 3.5848906801945047`*^9,
186	3.5848935801190357`^9, 3.5848936688170233`^9, 3.5848939033674097`*^9, {
187	3.58489396555275^9, 3.584893992584134^9}, 3.5848940700482283`*^9, {
188	3.584894103054511^9, 3.584894120965992^9}, 3.584894210383068*^9,
189	3.584894423252715^9, 3.584894734499061^9, 3.58489564491859*^9,
190	3.584895914221643*^9}]
191	}, Open ]],
192
193	Cell[BoxData["\[IndentingNewLine]"], "Input",
194	CellChangeTimes->{3.584887428382176^9, 3.584891660811152^9}],
195
196	Cell[BoxData[
197	RowBox[{"\[IndentingNewLine]",
198	RowBox[{"(*", " ",
199	RowBox[{"FRICTION", " ", "FACTORS", " ", "via", " ", "Wikipedia"}], " ",
200	"*)"}]}]], "Input",
201	CellChangeTimes->{{3.584887434985119^9, 3.58488743942144^9}, {
202	3.5848907429833384`^9, 3.584890746521543^9}}],
203
204	Cell[BoxData[{
205	RowBox[{
206	RowBox[{"xi", " ", "=",
207	RowBox[{
208	RowBox[{"Sqrt", "[",
209	RowBox[{
210	RowBox[{"(",
211	RowBox[{"p", "^", "2"}], ")"}], "-", "1"}], "]"}], "/", "p"}]}],
212	";"}], "\[IndentingNewLine]",
213	RowBox[{
214	RowBox[{"s", " ", "=", " ",
215	RowBox[{"2", "*",
216	RowBox[{
217	RowBox[{"ArcTanh", "[", "xi", "]"}], "/", "xi"}]}]}],
218	";"}], "\[IndentingNewLine]",
219	RowBox[{
220	RowBox[{"p", " ", "=", " ",
221	RowBox[{"a", "/", "b"}]}], ";"}], "\[IndentingNewLine]",
222	RowBox[{
223	RowBox[{"re", " ", "=", " ",
224	RowBox[{
225	RowBox[{"(",
226	RowBox[{"a", "*",
227	RowBox[{"(",
228	RowBox[{"b", "^", "2"}], ")"}]}], ")"}], "^",
229	RowBox[{"(",
230	RowBox[{"1", "/", "3"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
231	RowBox[{
232	RowBox[{"eqSphere", " ", "=", " ",
233	RowBox[{"8", "", "Pi", "", "eta", "*",
234	RowBox[{"(",
235	RowBox[{"re", "^", "3"}], ")"}]}]}], ";"}]}], "Input",
236	CellChangeTimes->{{3.584887453255073^9, 3.584887474611525^9}, {
237	3.584887518352754^9, 3.5848875201183167`^9}, {3.5848875705654163`*^9,
238	3.584887612684559^9}, {3.584891207938551^9, 3.5848912105877657`*^9}, {
239	3.58489477575947^9, 3.584894789206271^9}, {3.5848948909071417`*^9,
240	3.5848949045983353`^9}, {3.584895894748207^9, 3.58489594161933*^9}, {
241	3.584895989916945^9, 3.584896035299241^9}}],
242
243	Cell[BoxData[""], "Input",
244	CellChangeTimes->{{3.5848876347988167`^9, 3.584887638277546^9}}],
245
246	Cell[BoxData[
247	RowBox[{"(*", " ",
248	RowBox[{"ellipsoid", " ", "axial", " ", "friction", " ",
249	RowBox[{"factor", " ", "--"}], " ", "about", " ", "long", " ", "axis"}],
250	" ", "*)"}]], "Input",
251	CellChangeTimes->{{3.58488764265452^9, 3.584887652272901^9}, {
252	3.58489128903328^9, 3.5848912914194727`^9}}],
253
254	Cell[CellGroupData[{
255
256	Cell[BoxData[
257	RowBox[{"Fax", " ", "=", " ",
258	RowBox[{
259	RowBox[{"(",
260	RowBox[{"4", "/", "3"}], ")"}], "*",
261	RowBox[{
262	RowBox[{"(",
263	RowBox[{"xi", "^", "2"}], ")"}], "/",
264	RowBox[{"(",
265	RowBox[{"2", "-",
266	RowBox[{"(",
267	RowBox[{"s", "/",
268	RowBox[{"(",
269	RowBox[{"p", "^", "2"}], ")"}]}], ")"}]}], ")"}]}]}]}]], "Input",
270	CellChangeTimes->{{3.58488765759788^9, 3.584887700018955^9}, {
271	3.5848879174727783`^9, 3.584887938981941^9}, {3.584890603482909*^9,
272	3.5848906063233852`^9}, 3.584890946638309^9, {3.584894926277274*^9,
273	3.584894983989938*^9}}],
274
275	Cell[BoxData["0.7687260259259584`"], "Output",
276	CellChangeTimes->{
277	3.58488770586084^9, {3.58488792206218^9, 3.584887939492504*^9},
278	3.584890608261299^9, 3.5848906924156^9, 3.584890948261269*^9,
279	3.584891214003551^9, 3.584893997414321^9, 3.584894794095375*^9,
280	3.5848948750854387`^9, 3.584894906946406^9, {3.584894984628289*^9,
281	3.5848950094024963`^9}, 3.584895649579417^9, 3.5848960393110647`*^9}]
282	}, Open ]],
283
284	Cell[BoxData[
285	RowBox[{"\n",
286	RowBox[{"(*", " ",
287	RowBox[{"ellipsoid", " ", "equatorial", " ", "friction", " ",
288	RowBox[{"factor", " ", "--"}], " ", "about", " ", "short", " ", "axes"}],
289	" ", "*)"}]}]], "Input",
290	CellChangeTimes->{{3.584887710741127^9, 3.584887713271037^9},
291	3.584887813383315^9, {3.584891294243541^9, 3.5848912960435057`*^9}}],
292
293	Cell[CellGroupData[{
294
295	Cell[BoxData[
296	RowBox[{"Feq", " ", "=", " ",
297	RowBox[{
298	RowBox[{"(",
299	RowBox[{"4", "/", "3"}], ")"}], " ", "*", " ",
300	RowBox[{
301	RowBox[{"(",
302	RowBox[{
303	RowBox[{
304	RowBox[{"(",
305	RowBox[{"1", "/", "p"}], ")"}], "^", "2"}], " ", "-", " ",
306	RowBox[{"p", "^", "2"}]}], ")"}], "/",
307	RowBox[{"(",
308	RowBox[{"2", " ", "-", " ",
309	RowBox[{"s", "*",
310	RowBox[{"(",
311	RowBox[{"2", " ", "-", " ",
312	RowBox[{
313	RowBox[{"(",
314	RowBox[{"1", "/", "p"}], ")"}], "^", "2"}]}], ")"}]}]}],
315	")"}]}]}]}]], "Input",
316	CellChangeTimes->{{3.584887727157157^9, 3.584887769116691^9}, {
317	3.584887926989279^9, 3.584887942885796^9}, {3.584890609957733*^9,
318	3.584890611532776^9}, 3.584890951791182^9}],
319
320	Cell[BoxData["1.9247682885909356`"], "Output",
321	CellChangeTimes->{
322	3.5848877700852623`^9, {3.584887929629698^9, 3.5848879434378977`*^9},
323	3.584890612056086^9, 3.5848906935323477`^9, 3.584890952137286*^9,
324	3.584891216311064^9, 3.584893998534431^9, 3.584894876538658*^9,
325	3.5848950109231567`^9, 3.584895650617403^9, 3.584896040714982*^9}]
326	}, Open ]],
327
328	Cell[BoxData[
329	RowBox[{"\[IndentingNewLine]",
330	RowBox[{"(*", " ",
331	RowBox[{"ellipsoid", " ", "axial", " ", "rotational", " ", "friction"}],
332	" ", "*)"}]}]], "Input",
333	CellChangeTimes->{{3.58489104521352^9, 3.584891053491618^9}, {
334	3.584891299813236^9, 3.5848913004309483`^9}}],
335
336	Cell[CellGroupData[{
337
338	Cell[BoxData[
339	RowBox[{"fAx", " ", "=", " ",
340	RowBox[{"Fax", " ", "*", " ", "eqSphere"}]}]], "Input",
341	CellChangeTimes->{{3.584891056631095^9, 3.584891064865149^9}}],
342
343	Cell[BoxData["107748.71237555385`"], "Output",
344	CellChangeTimes->{3.5848910656263733`^9, 3.584891217398337^9,
345	3.58489400012298^9, 3.584895019483296^9, 3.584895652455336*^9,
346	3.584896042503642*^9}]
347	}, Open ]],
348
349	Cell[BoxData[
350	RowBox[{"(*", " ",
351	RowBox[{
352	"ellipsoid", " ", "equatorial", " ", "rotational", " ", "friction"}], " ",
353	"*)"}]], "Input",
354	CellChangeTimes->{{3.5848910141246853`^9, 3.584891022193733^9}, {
355	3.5848913017358427`^9, 3.5848913022309103`^9}}],
356
357	Cell[CellGroupData[{
358
359	Cell[BoxData[
360	RowBox[{"fEq", " ", "=", " ",
361	RowBox[{"Feq", " ", "*", " ", "eqSphere"}]}]], "Input",
362	CellChangeTimes->{{3.584891025068652^9, 3.584891029821259^9}}],
363
364	Cell[BoxData["269785.7204290194`"], "Output",
365	CellChangeTimes->{3.584891030572219^9, 3.584891218233328^9,
366	3.584894001208784^9, 3.584895020435354^9, 3.584895653375164*^9,
367	3.58489604404213*^9}]
368	}, Open ]]
369	},
370	WindowSize->{740, 876},
371	WindowMargins->{{Automatic, 190}, {Automatic, 15}},
372	FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 22, \
373	2012)",
374	StyleDefinitions->"Default.nb"
375	]
376	(* End of Notebook Content *)
377
378	(* Internal cache information *)
379	(*CellTagsOutline
380	CellTagsIndex->{}
381	*)
382	(*CellTagsIndex
383	CellTagsIndex->{}
384	*)
385	(*NotebookFileOutline
386	Notebook[{
387	Cell[557, 20, 1053, 25, 103, "Input"],
388	Cell[1613, 47, 204, 4, 43, "Input"],
389	Cell[1820, 53, 164, 3, 27, "Input"],
390	Cell[CellGroupData[{
391	Cell[2009, 60, 223, 5, 27, "Input"],
392	Cell[2235, 67, 114, 1, 30, "Output"]
393	}, Open ]],
394	Cell[2364, 71, 335, 11, 27, "Input"],
395	Cell[2702, 84, 318, 6, 27, "Input"],
396	Cell[3023, 92, 524, 17, 43, "Input"],
397	Cell[CellGroupData[{
398	Cell[3572, 113, 568, 10, 27, "Input"],
399	Cell[4143, 125, 468, 7, 27, "Output"]
400	}, Open ]],
401	Cell[4626, 135, 343, 7, 43, "Input"],
402	Cell[4972, 144, 717, 22, 43, "Input"],
403	Cell[CellGroupData[{
404	Cell[5714, 170, 540, 10, 27, "Input"],
405	Cell[6257, 182, 453, 7, 27, "Output"]
406	}, Open ]],
407	Cell[6725, 192, 109, 1, 43, "Input"],
408	Cell[6837, 195, 284, 6, 43, "Input"],
409	Cell[7124, 203, 1342, 37, 88, "Input"],
410	Cell[8469, 242, 94, 1, 27, "Input"],
411	Cell[8566, 245, 311, 6, 27, "Input"],
412	Cell[CellGroupData[{
413	Cell[8902, 255, 615, 17, 27, "Input"],
414	Cell[9520, 274, 421, 6, 27, "Output"]
415	}, Open ]],
416	Cell[9956, 283, 364, 7, 43, "Input"],
417	Cell[CellGroupData[{
418	Cell[10345, 294, 782, 23, 27, "Input"],
419	Cell[11130, 319, 355, 5, 27, "Output"]
420	}, Open ]],
421	Cell[11500, 327, 289, 6, 43, "Input"],
422	Cell[CellGroupData[{
423	Cell[11814, 337, 169, 3, 27, "Input"],
424	Cell[11986, 342, 205, 3, 27, "Output"]
425	}, Open ]],
426	Cell[12206, 348, 264, 6, 27, "Input"],
427	Cell[CellGroupData[{
428	Cell[12495, 358, 169, 3, 27, "Input"],
429	Cell[12667, 363, 202, 3, 27, "Output"]
430	}, Open ]]
431	}
432	]
433	*)
434
435	(* End of internal cache information *)