(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 19181, 645]*) (*NotebookOutlinePosition[ 20052, 673]*) (* CellTagsIndexPosition[ 20008, 669]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Problem #1: Kelvin-Helmholtz Instability\ \>", "Title", TextAlignment->Center], Cell["\<\ Two layer unbounded system with density contrast and shear. See \ Kundu, page 375. \ \>", "Text"], Cell[TextData[ "For z>0, the basic state is u=U1 and \[Rho]=\[Rho]1. For z<0, the basic \ state is u=U2 and \[Rho]=\[Rho]2. "], "Text"], Cell[TextData[{ "Your tasks:\n", StyleBox[" \[HappySmiley] ", FontSize->24], "For disturbances of form ", Cell[BoxData[ \(TraditionalForm\`e\^\(i(kx - \[Omega]t)\)\)]], " with k>0, find \[Omega](k).\n ", StyleBox["\[HappySmiley] ", FontSize->24], "For the special case of ", StyleBox["case1", FontFamily->"Courier", FontWeight->"Bold"], " , plot Im[\[Omega](k)]. Your\n graph might look like: \n ", Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.238095 0.0147151 0.327516 [ [.2619 .00222 -3 -9 ] [.2619 .00222 3 0 ] [.5 .00222 -3 -9 ] [.5 .00222 3 0 ] [.7381 .00222 -3 -9 ] [.7381 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 10 6 ] [.01131 .09659 -24 -4.5 ] [.01131 .09659 0 4.5 ] [.01131 .17847 -18 -4.5 ] [.01131 .17847 0 4.5 ] [.01131 .26035 -24 -4.5 ] [.01131 .26035 0 4.5 ] [.01131 .34223 -6 -4.5 ] [.01131 .34223 0 4.5 ] [.01131 .42411 -24 -4.5 ] [.01131 .42411 0 4.5 ] [.01131 .50599 -18 -4.5 ] [.01131 .50599 0 4.5 ] [.01131 .58787 -24 -4.5 ] [.01131 .58787 0 4.5 ] [.02381 .64303 -17.0625 0 ] [.02381 .64303 17.0625 12 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .01472 m .2619 .02097 L s [(1)] .2619 .00222 0 1 Mshowa .5 .01472 m .5 .02097 L s [(2)] .5 .00222 0 1 Mshowa .7381 .01472 m .7381 .02097 L s [(3)] .7381 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(4)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .21429 .01472 m .21429 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .40476 .01472 m .40476 .01847 L s .45238 .01472 m .45238 .01847 L s .54762 .01472 m .54762 .01847 L s .59524 .01472 m .59524 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .78571 .01472 m .78571 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (k) show 69.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .02381 .09659 m .03006 .09659 L s [(0.25)] .01131 .09659 1 0 Mshowa .02381 .17847 m .03006 .17847 L s [(0.5)] .01131 .17847 1 0 Mshowa .02381 .26035 m .03006 .26035 L s [(0.75)] .01131 .26035 1 0 Mshowa .02381 .34223 m .03006 .34223 L s [(1)] .01131 .34223 1 0 Mshowa .02381 .42411 m .03006 .42411 L s [(1.25)] .01131 .42411 1 0 Mshowa .02381 .50599 m .03006 .50599 L s [(1.5)] .01131 .50599 1 0 Mshowa .02381 .58787 m .03006 .58787 L s [(1.75)] .01131 .58787 1 0 Mshowa .125 Mabswid .02381 .03109 m .02756 .03109 L s .02381 .04747 m .02756 .04747 L s .02381 .06384 m .02756 .06384 L s .02381 .08022 m .02756 .08022 L s .02381 .11297 m .02756 .11297 L s .02381 .12935 m .02756 .12935 L s .02381 .14572 m .02756 .14572 L s .02381 .1621 m .02756 .1621 L s .02381 .19485 m .02756 .19485 L s .02381 .21122 m .02756 .21122 L s .02381 .2276 m .02756 .2276 L s .02381 .24398 m .02756 .24398 L s .02381 .27673 m .02756 .27673 L s .02381 .2931 m .02756 .2931 L s .02381 .30948 m .02756 .30948 L s .02381 .32586 m .02756 .32586 L s .02381 .35861 m .02756 .35861 L s .02381 .37498 m .02756 .37498 L s .02381 .39136 m .02756 .39136 L s .02381 .40773 m .02756 .40773 L s .02381 .44049 m .02756 .44049 L s .02381 .45686 m .02756 .45686 L s .02381 .47324 m .02756 .47324 L s .02381 .48961 m .02756 .48961 L s .02381 .52237 m .02756 .52237 L s .02381 .53874 m .02756 .53874 L s .02381 .55512 m .02756 .55512 L s .02381 .57149 m .02756 .57149 L s .02381 .60424 m .02756 .60424 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .02381 .64303 -78.0625 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (Im) show %%IncludeResource: font Math2Mono %%IncludeFont: Math2Mono /Math2Mono findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 75.125000 13.000000 moveto (@) show 81.125000 13.000000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (w) show %%IncludeResource: font Math2Mono %%IncludeFont: Math2Mono /Math2Mono findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 87.125000 13.000000 moveto (D) show 93.125000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Math1Mono %%+ font Math2Mono %%+ font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Math1Mono %%+ Math2Mono %%+ Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Math1Mono %%+ font Math2Mono %%+ font Courier grestore 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .01472 m .06244 .01472 L .10458 .01472 L .14415 .01472 L .18221 .01472 L .19227 .01472 L .19721 .01472 L .20178 .01472 L .20433 .01472 L .20562 .01472 L .20635 .01472 L .20704 .01832 L .20771 .02318 L .20845 .02638 L .20979 .03068 L .21133 .03456 L .21274 .03759 L .22272 .0533 L .24239 .07529 L .26336 .0947 L .30494 .12889 L .34501 .15956 L .38356 .18808 L .42456 .21782 L .46405 .24609 L .50202 .27306 L .54245 .30159 L .58136 .32892 L .62272 .35788 L .66256 .38571 L .70089 .41242 L .74167 .44079 L .78094 .46806 L .81869 .49426 L .8589 .52213 L .89758 .54892 L .93872 .5774 L .97619 .60332 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgGoo 00<007ooOol0=Goo00<007ooOol067oo001GOol00`00Oomoo`0hOol00`00Oomoo`0hOol00`00Oomo o`0cOol4000IOol005Moo`03001oogoo03Uoo`03001oogoo03Eoo`8003Ioo`03001oo`0001Yoo`00 Egoo00<007ooOol0=Woo00@007ooOol003Uoo`03001oogoo03=oo`03001oo`0001Yoo`00EWoo0P00 >Goo0P00=goo0`00=goo0P006Woo000NOol00`00Oomoo`3oOol001ioo`03001oogoo0?moo`007Woo 00<007ooOol0lgoo00@007oo000000Qoo`007Woo00<007ooOol0lgoo00<007oo00002Goo000HOood 0008Ool30009Ool001ioo`03001oogoo00Qoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo 00Moo`03001oo`0000]oo`03001oogoo00Qoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo 00Uoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo 00Uoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo00Uoo`03001oogoo00Qoo`03001oogoo 00Uoo`03001oogoo00Qoo`03001oogoo00aoo`03001oo`0000Uoo`007Woo00<007ooOol0:Goo0P00 b7oo00@007oo000000Qoo`007Woo00<007ooOol0:Woo00<007ooOol0aWoo00<007ooOol02Goo000N Ool2000[Ool00`00Oomoo`35Ool2000;Ool001ioo`03001oogoo02]oo`03001oogoo0=5oo`007Woo 00<007ooOol0:goo00<007ooOol0dGoo000NOol00`00Oomoo`0/Ool00`00Oomoo`3@Ool001ioo`80 02eoo`03001oogoo0=1oo`007Woo00<007ooOol0;Goo00<007ooOol0cgoo000NOol00`00Oomoo`0^ Ool00`00Oomoo`3>Ool001ioo`03001oogoo02moo`03001oogoo0Woo00<007ooOol0`goo0004Ool01000Oomoo`002goo00<007ooOol00`001Goo00<0 07ooOol0>Woo00<007ooOol0`Woo0004Ool01000Oomoo`0027oo00@007ooOol0009oo`03001oogoo 00Eoo`03001oogoo03]oo`03001oogoo0<5oo`001Goo0P002Woo0P000goo100017oo00<007ooOol0 ?7oo0P00`Goo000NOol2000oOol00`00Oomoo`2nOol001ioo`03001oogoo03moo`03001oogoo0;eo o`007Woo00<007ooOol0@7oo00<007ooOol0_7oo000NOol20012Ool2002lOol001ioo`03001oogoo 04=oo`03001oogoo0;Uoo`007Woo00<007ooOol0A7oo00<007ooOol0^7oo000NOol00`00Oomoo`15 Ool00`00Oomoo`2gOol001ioo`8004Moo`800;Moo`007Woo00<007ooOol0B7oo00<007ooOol0]7oo 000NOol00`00Oomoo`19Ool00`00Oomoo`2cOol001ioo`03001oogoo04Yoo`03001oogoo0;9oo`00 7Woo0P00C7oo0P00/Woo000NOol00`00Oomoo`1=Ool00`00Oomoo`2_Ool000]oo`8000Aoo`8000=o o`<000Eoo`03001oogoo04ioo`03001oogoo0:ioo`002Woo00@007ooOol000=oo`8000Ioo`03001o ogoo009oo`03001oogoo04moo`800:ioo`002Woo00@007ooOol000]oo`03001oogoo009oo`80059o o`03001oogoo0:]oo`002Woo00@007ooOol000Qoo`<000Eoo`03001oogoo059oo`03001oogoo0:Yo o`002Woo00@007ooOol000Qoo`03001oogoo00Eoo`03001oogoo05=oo`03001oogoo0:Uoo`002goo 0P002Goo100017oo00<007ooOol0E7oo0P00ZGoo000NOol2001GOol00`00Oomoo`2VOol001ioo`03 001oogoo05Moo`800:Ioo`007Woo00<007ooOol0FGoo00<007ooOol0Xgoo000NOol00`00Oomoo`1J Ool00`00Oomoo`2ROol001ioo`8005aoo`800:9oo`007Woo00<007ooOol0GGoo00<007ooOol0Wgoo 000NOol00`00Oomoo`1NOol2002OOol001ioo`03001oogoo061oo`03001oogoo09aoo`007Woo0P00 HWoo0P00W7oo000NOol00`00Oomoo`1SOol00`00Oomoo`2IOol001ioo`03001oogoo06Aoo`03001o ogoo09Qoo`007Woo00<007ooOol0IGoo0P00V7oo000NOol2001XOol00`00Oomoo`2EOol001ioo`03 001oogoo06Qoo`8009Eoo`001Goo0P0017oo0P0017oo00<007ooOol00Woo0`001Goo00<007ooOol0 JWoo00<007ooOol0TWoo0004Ool01000Oomoo`000goo0P0017oo00<007ooOol01Goo00<007ooOol0 0Woo00<007ooOol0Jgoo0P00TWoo0004Ool01000Oomoo`002Woo00<007ooOol017oo00<007ooOol0 0Woo0P00KWoo00<007ooOol0Sgoo0004Ool01000Oomoo`002Woo00@007ooOomoo`<000Eoo`03001o ogoo06ioo`8008moo`0017oo00@007ooOol000]oo`04001oogoo0007Ool00`00Oomoo`1`Ool00`00 Oomoo`27oo000:Ool2000:Ool40004Ool00`00Oomoo`37Ool00`00Oomoo`0eOol0 01ioo`800Eoo`03001oogoo01Qoo`007Woo00<007ooOol0ogoo 000NOol00`00Oomoo`3oOol001ioo`03001oogoo0?moo`00ogoo8Goo003oOolQOol00?moob5oo`00 ogoo8Goo003oOolQOol00?moob5oo`00ogoo8Goo000NOol20008Ool2003fOol001ioo`03001oogoo 00Qoo`03001oogoo0?Aoo`0047oo1@0000Ioo`00Ool007oo0003Ool01000Oomoogoo10000goo00<0 07ooOol0m7oo000BOol02@00Oomoogoo001oo`00Ool00003Ool01000Oomoo`000Woo00<007oo0000 0Woo00<007ooOol0m7oo000BOol02@00Oomoogoo001oo`00Ool00003Ool01000Oomoo`000Woo00<0 07oo00000Woo00<007ooOol0m7oo000BOol02@00Oomoogoo001oo`00Ool00003Ool01000Oomoo`00 17oo00@007ooOol00?Ioo`004Woo00<007ooOol01P000goo00D007ooOomoo`0000=oo`04001oogoo 003fOol0011oo`D000Uoo`8000Qoo`800?Ioo`00ogoo8Goo003oOolQOol00?moob5oo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0, 0, 0, 0}}], "\n ", StyleBox["case1", FontFamily->"Courier", FontWeight->"Bold"], " could correspond to a layer of fresh water over salt water with 1 m/s \ shear across the interface.\n ", StyleBox["\[HappySmiley] ", FontSize->24], "Do something else or make some intelligent comments about the problem\n \ you have solved. Make a nice notebook. Put it on your web page." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(case1 = {U1 -> 0, U2 -> 1, \[Rho]1 -> 1, \[Rho]2 -> 1.04, g -> 9.8}\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ \(U1 \[Rule] 0\), ",", \(U2 \[Rule] 1\), ",", \(\[Rho]1 \[Rule] 1\), ",", RowBox[{"\[Rho]2", "\[Rule]", StyleBox["1.04000000000000003`", StyleBoxAutoDelete->True, PrintPrecision->3]}], ",", RowBox[{"g", "\[Rule]", StyleBox["9.80000000000000071`", StyleBoxAutoDelete->True, PrintPrecision->2]}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[TextData[ "exampleofuse = (g*(-\[Rho]1^2 + \[Rho]2^2))/((U1 - \ U2)^2*\[Rho]1*\[Rho]2)/.case1"], "Input"], Cell[BoxData[ \(0.768923076923078063`\)], "Output"] }, Open ]] }, Open ]] }, FrontEndVersion->"Macintosh 3.0", ScreenRectangle->{{0, 1152}, {0, 850}}, WindowSize->{589, 757}, WindowMargins->{{150, Automatic}, {Automatic, 33}}, StyleDefinitions -> "ok.nb", MacintoshSystemPageSetup->"\<\ 00<0001804P000000]P2:?oQon82n@960dL5:0?l0080001804P000000]P2:001 0000I00000400`<300000BL?00400@00000000000000060801T1T000000@0000 00000000004000000000000000000000\>" ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1731, 51, 90, 4, 78, "Title"], Cell[1824, 57, 109, 3, 70, "Text"], Cell[1936, 62, 136, 2, 70, "Text"], Cell[2075, 66, 16278, 545, 381, "Text"], Cell[CellGroupData[{ Cell[18378, 615, 103, 2, 47, "Input"], Cell[18484, 619, 464, 12, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18985, 636, 110, 2, 63, "Input"], Cell[19098, 640, 55, 1, 52, "Output"] }, Open ]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)