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 (11 years, 1 month ago) by kstocke1
File size:	14685 byte(s)
Log Message:
#	User	Rev	Content
1	kstocke1	3933	(* 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 *)