stefan-problem.nb

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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     29152,        734]
NotebookOptionsPosition[     27006,        659]
NotebookOutlinePosition[     27473,        676]
CellTagsIndexPosition[     27430,        673]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 "Two Phase Stefan Problem ",
 StyleBox["Mathematica",
  FontSlant->"Italic"],
 " Code"
}], "Title",
 CellChangeTimes->{{3.658907084181588*^9, 3.6589071026296434`*^9}}],

Cell[CellGroupData[{

Cell["Muhammad Zeeshan Khalid", "Author",
 CellChangeTimes->{{3.6589071064798636`*^9, 3.658907111647159*^9}, {
  3.797172466993026*^9, 3.797172468657581*^9}}],

Cell["\<\
University of Engineering and Technology Taxila\
\>", "Institution",
 CellChangeTimes->{{3.6589071142613087`*^9, 3.658907124844914*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["", "Author",
 CellChangeTimes->{3.6589071297351933`*^9}],

Cell["", "Institution",
 CellChangeTimes->{3.658907135839543*^9}],

Cell[CellGroupData[{

Cell["Abstract", "AbstractSection",
 CellChangeTimes->{{3.658907140062784*^9, 3.65890714156487*^9}}],

Cell["\<\
Two Phase Stefan Problem was solved using eigen function expansion method. \
Transcendental equation was solved and eigen values obtained were graphically \
validated. In order to find the coefficients Numerical integration scheme was \
applied and integrated from 0 to 1. Numerical differentiation scheme was \
applied to obtain the interface position of the melted material. It was found \
numerical scheme reduced the computational time significantly and made the \
problem less time consuming and tracable.\
\>", "Abstract",
 CellChangeTimes->{{3.6589072237655716`*^9, 3.6589074397699265`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["Mathematica",
  FontSlant->"Italic"],
 " code"
}], "Section",
 CellChangeTimes->{{3.6589074443301873`*^9, 3.6589074821593513`*^9}}],

Cell[CellGroupData[{

Cell["\<\
Transcendental equation solution for Phase two\
\>", "Subsection",
 CellChangeTimes->{{3.6589074853535337`*^9, 3.658907495698126*^9}, {
  3.658907749358634*^9, 3.6589077520467877`*^9}}],

Cell["\<\
Transcendental equation was solved by assuming normalized initial interface \
location and and outer cylindrical boundary at 2.\
\>", "Text",
 CellChangeTimes->{{3.6589075225546618`*^9, 3.658907585014234*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{"eqn1", ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"BesselJ", "[", 
      RowBox[{"0", ",", "\[Beta]"}], "]"}], 
     RowBox[{"BesselY", "[", 
      RowBox[{"1", ",", 
       RowBox[{"\[Beta]", "*", "2"}]}], "]"}]}], "-", 
    RowBox[{
     RowBox[{"BesselJ", "[", 
      RowBox[{"1", ",", 
       RowBox[{"\[Beta]", "*", "2"}]}], "]"}], 
     RowBox[{"BesselY", "[", 
      RowBox[{"0", ",", "\[Beta]"}], "]"}]}]}]}], 
  TraditionalForm]], "Equation",
 CellChangeTimes->{
  3.6589075900325212`*^9, {3.658907640665417*^9, 3.6589076440316095`*^9}}],

Cell["\<\
Plot[Evaluate[eqn1],{\[Beta],0,30},Frame\[Rule]False, \
FrameLabel\[Rule]{Style[\[OpenCurlyDoubleQuote]\[Beta]\[CloseCurlyDoubleQuote]\
,16],Style[ \[OpenCurlyDoubleQuote]Transdental Equation\
\[CloseCurlyDoubleQuote], 16]},PlotStyle\[Rule]{Thick,Red} ]\
\>", "Text",
 CellChangeTimes->{{3.6589076514410334`*^9, 3.6589076820117817`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Roots of the equation", "Subsection",
 CellChangeTimes->{{3.658907685524983*^9, 3.6589077047680836`*^9}}],

Cell["\<\
Clear[root, guess, rootlist1];
guess = 1; rootlist1 = {}; n = 20;
Do[{root = FindRoot[Evaluate[eqn1 == 0], {\[Beta], guess}];
   AppendTo[rootlist1, root[[1, 2]]];
   guess = 3 + root[[1, 2]]}, {i, 1, n}];
rootlist1\
\>", "Text",
 CellChangeTimes->{
  3.658907728954467*^9, {3.6589078888556128`*^9, 3.658907888860613*^9}}],

Cell[TextData[Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "2.404825557695773`", ",", "5.520078110286308`", ",", "8.653727912911013`", 
   ",", "11.791534439014281`", ",", "14.930917708487785`", ",", 
   "18.071063967910924`", ",", "21.21163662987926`", ",", 
   "24.352471530749302`", ",", "27.493479132040253`", ",", 
   "30.634606468431976`", ",", "33.77582021357357`", ",", 
   "36.917098353664045`", ",", "40.05842576462824`", ",", 
   "43.19979171317673`", ",", "46.341188371661815`", ",", 
   "49.482609897397815`", ",", "52.624051841115`", ",", "55.76551075501998`", 
   ",", "58.90698392608094`", ",", "62.048469190227166`"}], "}"}]], "Output",
 CellChangeTimes->{3.654595579034833*^9, 3.654598357149732*^9, 
  3.654599912530695*^9, 3.6546000003007154`*^9, 3.6546023089597626`*^9, 
  3.6546027838029222`*^9, 3.6546028235291944`*^9, 3.6546029242859573`*^9, 
  3.6546033291271133`*^9, 3.654607334177189*^9, 3.654607379766796*^9, 
  3.65556312581509*^9, 3.6556317998372645`*^9, 3.655632709566298*^9, 
  3.655800397559953*^9, 3.6559923000432205`*^9}]], "Text",
 CellChangeTimes->{
  3.658907728954467*^9, {3.6589078888556128`*^9, 3.658907900599284*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Transcendental equation for Phase One\
\>", "Subsection",
 CellChangeTimes->{{3.658907737884978*^9, 3.6589077456324205`*^9}}],

Cell[TextData[{
 "eqn2:=BesselJ[0,g]\nPlot[Evaluate[eqn2],{g,0,30},Frame\[Rule]False, \
FrameLabel\[Rule]{Style[\"g\",16],Style[ \"Transdental Equation\", \
16]},PlotStyle\[Rule]{Thick,Green} ]\n\n",
 StyleBox["Roots of the equation\n",
  FontSize->16,
  FontWeight->"Bold"],
 "Clear[root,guess,rootlist2];\nguess=1;rootlist2={};n=20;\n\
Do[{root=FindRoot[Evaluate[eqn2==0],{g,guess}];\n\
AppendTo[rootlist2,root[[1,2]]];\nguess=3+root[[1,2]]},{i,1,n}];\nrootlist2\n",
 Cell[BoxData[
  RowBox[{"{", 
   RowBox[{
   "2.404825557695773`", ",", "5.520078110286308`", ",", "8.653727912911013`",
     ",", "11.791534439014281`", ",", "14.930917708487785`", ",", 
    "18.071063967910924`", ",", "21.21163662987926`", ",", 
    "24.352471530749302`", ",", "27.493479132040253`", ",", 
    "30.634606468431976`", ",", "33.77582021357357`", ",", 
    "36.917098353664045`", ",", "40.05842576462824`", ",", 
    "43.19979171317673`", ",", "46.341188371661815`", ",", 
    "49.482609897397815`", ",", "52.624051841115`", ",", "55.76551075501998`",
     ",", "58.90698392608094`", ",", "62.048469190227166`"}], "}"}]], "Output",
  CellChangeTimes->{3.654595579034833*^9, 3.654598357149732*^9, 
   3.654599912530695*^9, 3.6546000003007154`*^9, 3.6546023089597626`*^9, 
   3.6546027838029222`*^9, 3.6546028235291944`*^9, 3.6546029242859573`*^9, 
   3.6546033291271133`*^9, 3.654607334177189*^9, 3.654607379766796*^9, 
   3.65556312581509*^9, 3.6556317998372645`*^9, 3.655632709566298*^9, 
   3.655800397559953*^9, 3.6559923000432205`*^9}]
}], "Program",
 CellChangeTimes->{{3.658907762025358*^9, 3.658907853067566*^9}, {
  3.658907923141574*^9, 3.658907924058626*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Coefficients for Phase One and Two\
\>", "Section",
 CellChangeTimes->{{3.6589079806688643`*^9, 3.6589079842550697`*^9}, {
  3.658908018057003*^9, 3.65890802535142*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Beta]", "[", "m_", "]"}], ":=", 
  RowBox[{"rootlist1", "[", 
   RowBox[{"[", "m", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"g", "[", "n_", "]"}], ":=", 
  RowBox[{"rootlist2", "[", 
   RowBox[{"[", "n", "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.6545982786752434`*^9, 3.6545983083589416`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"b1", "[", "n_", "]"}], ":=", 
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    RowBox[{"R", "*", 
     RowBox[{"BesselJ", "[", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"g", "[", "n", "]"}], "R"}]}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b2", "[", "n_", "]"}], ":=", 
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"BesselJ", "[", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"g", "[", "n", "]"}], "R"}]}], "]"}], "2"], ",", 
    RowBox[{"{", 
     RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"B", "[", "n_", "]"}], ":=", 
  RowBox[{
   RowBox[{"B", "[", "n", "]"}], "=", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{"b1", "[", "n", "]"}], 
     RowBox[{"b2", "[", "n", "]"}]], ")"}]}]}]}], "Input",
 CellChangeTimes->{{3.6502935714909334`*^9, 3.650293582703575*^9}, {
   3.650293646493223*^9, 3.650293739497543*^9}, {3.6502937907564745`*^9, 
   3.650293837068124*^9}, {3.6502939736499357`*^9, 3.650293974194967*^9}, {
   3.6502943947910233`*^9, 3.650294396619128*^9}, {3.65029444910213*^9, 
   3.650294454294427*^9}, {3.650294515928952*^9, 3.650294761487997*^9}, {
   3.650294809484743*^9, 3.6502949107605352`*^9}, {3.650294984164734*^9, 
   3.6502949908921185`*^9}, {3.650295023127962*^9, 3.6502950670514746`*^9}, {
   3.6502951048266354`*^9, 3.650295160247805*^9}, {3.6502956559971604`*^9, 
   3.65029570319086*^9}, {3.650297242844923*^9, 3.650297249749318*^9}, {
   3.6502974427213554`*^9, 3.6502974871858983`*^9}, {3.6502982090111847`*^9, 
   3.650298238338862*^9}, {3.650298292761975*^9, 3.650298442316529*^9}, {
   3.6502985751641273`*^9, 3.650298604495805*^9}, {3.65030046972749*^9, 
   3.6503004862164335`*^9}, {3.650300526654746*^9, 3.650300597758813*^9}, {
   3.650300735741705*^9, 3.6503008119050617`*^9}, {3.650302470906951*^9, 
   3.6503025188916955`*^9}, {3.6503027200372005`*^9, 
   3.6503027379092226`*^9}, {3.650302789772189*^9, 3.6503029224067755`*^9}, {
   3.650303047448927*^9, 3.6503030617667465`*^9}, {3.6503032353046722`*^9, 
   3.650303236232725*^9}, {3.6503033058227053`*^9, 3.6503033495912085`*^9}, {
   3.6503039254591465`*^9, 3.6503040234327507`*^9}, {3.650304137819293*^9, 
   3.6503042016669445`*^9}, {3.650423057336937*^9, 3.650423180268968*^9}, {
   3.6504232628046885`*^9, 3.6504232851589675`*^9}, {3.6504233313206077`*^9, 
   3.650423337075937*^9}, {3.6504233699028144`*^9, 3.650423373857041*^9}, {
   3.6504234454721365`*^9, 3.6504234493323574`*^9}, {3.650423538990486*^9, 
   3.6504235486490383`*^9}, {3.650423581250903*^9, 3.650423653845055*^9}, {
   3.6504236941393595`*^9, 3.650423764874406*^9}, {3.650424187296567*^9, 
   3.650424198685218*^9}, {3.65046427645275*^9, 3.650464277234795*^9}, {
   3.650464349799946*^9, 3.6504643504729843`*^9}, {3.6504644931101427`*^9, 
   3.650464522683834*^9}, {3.6504660820510244`*^9, 3.6504662406400957`*^9}, {
   3.6504667498572206`*^9, 3.6504667789278836`*^9}, {3.6506364724534473`*^9, 
   3.650636519084114*^9}, {3.650636648296505*^9, 3.6506366689226847`*^9}, {
   3.65063868067675*^9, 3.6506386818998203`*^9}, {3.650639128642372*^9, 
   3.650639144081256*^9}, {3.6506405787873163`*^9, 3.6506406444320707`*^9}, {
   3.650800921147909*^9, 3.650800933436612*^9}, {3.6509942772842827`*^9, 
   3.6509942863508015`*^9}, {3.6509943648392906`*^9, 
   3.6509943977421727`*^9}, {3.6509945054093313`*^9, 3.650994510147602*^9}, {
   3.6509945720891447`*^9, 3.6509945758603606`*^9}, {3.650994832421035*^9, 
   3.6509948403284874`*^9}, {3.65099525823139*^9, 3.6509952584704037`*^9}, {
   3.650996531783633*^9, 3.6509965373939543`*^9}, 3.6509967346342354`*^9, {
   3.650996791829507*^9, 3.650996846539636*^9}, {3.6509970596138234`*^9, 
   3.650997163808783*^9}, {3.6509974245076942`*^9, 3.6509974326971626`*^9}, 
   3.651571837420647*^9, {3.651571873045685*^9, 3.6515718907016945`*^9}, {
   3.651571932715098*^9, 3.651571933251128*^9}, {3.6515719789187403`*^9, 
   3.651572015361825*^9}, {3.6515721094632072`*^9, 3.6515721100462403`*^9}, {
   3.652179406904944*^9, 3.6521794085190363`*^9}, {3.653653149593973*^9, 
   3.6536531576294327`*^9}, {3.653739740647565*^9, 3.6537398415953393`*^9}, {
   3.654063638066285*^9, 3.654063655371275*^9}, {3.6540758647816133`*^9, 
   3.654075884733754*^9}, {3.6541611435993176`*^9, 3.654161158747184*^9}, {
   3.654163185067083*^9, 3.654163193748579*^9}, {3.654163236399019*^9, 
   3.6541632469666233`*^9}, {3.654163304660923*^9, 3.6541633322975035`*^9}, {
   3.654170005428185*^9, 3.654170008808378*^9}, {3.6542549791457205`*^9, 
   3.654254998275814*^9}, {3.6542551789881506`*^9, 3.6542552015224395`*^9}, {
   3.6542573170614414`*^9, 3.6542573246148734`*^9}, {3.654259782545459*^9, 
   3.654259790045888*^9}, {3.654598771346423*^9, 3.6545989309685526`*^9}, {
   3.6546073551533885`*^9, 3.654607362381802*^9}, {3.655800403590298*^9, 
   3.6558004039833207`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"B", "[", "n", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "1", ",", "10"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.654598948543558*^9, 3.6545989598662057`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.4559469741516875`", ",", 
   RowBox[{"-", "0.24435287574709988`"}], ",", "0.1769751628633535`", ",", 
   RowBox[{"-", "0.1424615373290745`"}], ",", "0.12102716239544418`", ",", 
   RowBox[{"-", "0.10622764923071092`"}], ",", "0.09529771187453798`", ",", 
   RowBox[{"-", "0.08684041206885848`"}], ",", "0.08006865070402151`", ",", 
   RowBox[{"-", "0.07450269020010696`"}]}], "}"}]], "Output",
 CellChangeTimes->{3.654598963330404*^9, 3.654599920889173*^9, 
  3.654600009112219*^9, 3.6546023206554317`*^9, 3.6546027900472794`*^9, 
  3.6546028307076054`*^9, 3.654602931431366*^9, 3.654603336006507*^9, 
  3.654607341330598*^9, 3.6546073950696716`*^9, 3.6555631308043756`*^9, 
  3.6556318036134806`*^9, 3.6556327161116724`*^9, 3.6558004076315293`*^9, 
  3.655992305265519*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"c1", "[", "m_", "]"}], ":=", 
  RowBox[{
   RowBox[{"BesselY", "[", 
    RowBox[{"0", ",", 
     RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"R", "*", 
      RowBox[{"BesselJ", "[", 
       RowBox[{"0", ",", 
        RowBox[{
         RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"c2", "[", "m_", "]"}], ":=", 
  RowBox[{
   RowBox[{"BesselJ", "[", 
    RowBox[{"0", ",", 
     RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"R", "*", 
      RowBox[{"BesselY", "[", 
       RowBox[{"0", ",", 
        RowBox[{
         RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b3", "[", "m_", "]"}], ":=", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"BesselY", "[", 
     RowBox[{"0", ",", 
      RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], "2"], 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"R", "*", 
      SuperscriptBox[
       RowBox[{"BesselJ", "[", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}], "2"]}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b4", "[", "m_", "]"}], ":=", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"BesselJ", "[", 
     RowBox[{"0", ",", 
      RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], "2"], 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"R", "*", 
      SuperscriptBox[
       RowBox[{"BesselY", "[", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}], "2"]}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b5", "[", "m_", "]"}], ":=", 
  RowBox[{"2", 
   RowBox[{"BesselJ", "[", 
    RowBox[{"0", ",", 
     RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], 
   RowBox[{"BesselY", "[", 
    RowBox[{"0", ",", 
     RowBox[{"\[Beta]", "[", "m", "]"}]}], "]"}], 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"R", "*", 
      RowBox[{"BesselJ", "[", 
       RowBox[{"0", ",", 
        RowBox[{
         RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}], 
      RowBox[{"BesselY", "[", 
       RowBox[{"0", ",", 
        RowBox[{
         RowBox[{"\[Beta]", "[", "m", "]"}], "R"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0", ",", "1"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"A", "[", "m_", "]"}], ":=", 
  RowBox[{
   RowBox[{"A", "[", "m", "]"}], "=", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"c1", "[", "m", "]"}], "-", 
       RowBox[{"c2", "[", "m", "]"}]}], ")"}], 
     RowBox[{
      RowBox[{"b3", "[", "m", "]"}], "+", 
      RowBox[{"b4", "[", "m", "]"}], "-", 
      RowBox[{"b5", "[", "m", "]"}]}]], ")"}]}]}]}], "Input",
 CellChangeTimes->{{3.6502935714909334`*^9, 3.650293582703575*^9}, {
   3.650293646493223*^9, 3.650293739497543*^9}, {3.6502937907564745`*^9, 
   3.650293837068124*^9}, {3.6502939736499357`*^9, 3.650293974194967*^9}, {
   3.6502943947910233`*^9, 3.650294396619128*^9}, {3.65029444910213*^9, 
   3.650294454294427*^9}, {3.650294515928952*^9, 3.650294761487997*^9}, {
   3.650294809484743*^9, 3.6502949107605352`*^9}, {3.650294984164734*^9, 
   3.6502949908921185`*^9}, {3.650295023127962*^9, 3.6502950670514746`*^9}, {
   3.6502951048266354`*^9, 3.650295160247805*^9}, {3.6502956559971604`*^9, 
   3.65029570319086*^9}, {3.650297242844923*^9, 3.650297249749318*^9}, {
   3.6502974427213554`*^9, 3.6502974871858983`*^9}, {3.6502982090111847`*^9, 
   3.650298238338862*^9}, {3.650298292761975*^9, 3.650298442316529*^9}, {
   3.6502985751641273`*^9, 3.650298604495805*^9}, {3.65030046972749*^9, 
   3.6503004862164335`*^9}, {3.650300526654746*^9, 3.650300597758813*^9}, {
   3.650300735741705*^9, 3.6503008119050617`*^9}, {3.650302470906951*^9, 
   3.6503025188916955`*^9}, {3.6503027200372005`*^9, 
   3.6503027379092226`*^9}, {3.650302789772189*^9, 3.6503029224067755`*^9}, {
   3.650303047448927*^9, 3.6503030617667465`*^9}, {3.6503032353046722`*^9, 
   3.650303236232725*^9}, {3.6503033058227053`*^9, 3.6503033495912085`*^9}, {
   3.6503039254591465`*^9, 3.6503040234327507`*^9}, {3.650304137819293*^9, 
   3.6503042016669445`*^9}, {3.650423057336937*^9, 3.650423180268968*^9}, {
   3.6504232628046885`*^9, 3.6504232851589675`*^9}, {3.6504233313206077`*^9, 
   3.650423337075937*^9}, {3.6504233699028144`*^9, 3.650423373857041*^9}, {
   3.6504234454721365`*^9, 3.6504234493323574`*^9}, {3.650423538990486*^9, 
   3.6504235486490383`*^9}, {3.650423581250903*^9, 3.650423653845055*^9}, {
   3.6504236941393595`*^9, 3.650423764874406*^9}, {3.650424187296567*^9, 
   3.650424198685218*^9}, {3.65046427645275*^9, 3.650464277234795*^9}, {
   3.650464349799946*^9, 3.6504643504729843`*^9}, {3.6504644931101427`*^9, 
   3.650464522683834*^9}, {3.6504660820510244`*^9, 3.6504662406400957`*^9}, {
   3.6504667498572206`*^9, 3.6504667789278836`*^9}, {3.6506364724534473`*^9, 
   3.650636519084114*^9}, {3.650636648296505*^9, 3.6506366689226847`*^9}, {
   3.65063868067675*^9, 3.6506386818998203`*^9}, {3.650639128642372*^9, 
   3.650639144081256*^9}, {3.6506405787873163`*^9, 3.6506406444320707`*^9}, {
   3.650800921147909*^9, 3.650800933436612*^9}, {3.6509942772842827`*^9, 
   3.6509942863508015`*^9}, {3.6509943648392906`*^9, 
   3.6509943977421727`*^9}, {3.6509945054093313`*^9, 3.650994510147602*^9}, {
   3.6509945720891447`*^9, 3.6509945758603606`*^9}, {3.650994832421035*^9, 
   3.6509948403284874`*^9}, {3.65099525823139*^9, 3.6509952584704037`*^9}, {
   3.650996531783633*^9, 3.6509965373939543`*^9}, 3.6509967346342354`*^9, {
   3.650996791829507*^9, 3.650996846539636*^9}, {3.6509970596138234`*^9, 
   3.650997163808783*^9}, {3.6509974245076942`*^9, 3.6509974326971626`*^9}, 
   3.651571837420647*^9, {3.651571873045685*^9, 3.6515718907016945`*^9}, {
   3.651571932715098*^9, 3.651571933251128*^9}, {3.6515719789187403`*^9, 
   3.651572015361825*^9}, {3.6515721094632072`*^9, 3.6515721100462403`*^9}, {
   3.652179406904944*^9, 3.6521794085190363`*^9}, {3.653653149593973*^9, 
   3.6536531576294327`*^9}, {3.653739740647565*^9, 3.6537398415953393`*^9}, {
   3.654063638066285*^9, 3.654063655371275*^9}, {3.6540758647816133`*^9, 
   3.654075884733754*^9}, {3.6541611435993176`*^9, 3.654161158747184*^9}, {
   3.654163185067083*^9, 3.654163193748579*^9}, {3.654163236399019*^9, 
   3.6541632469666233`*^9}, {3.654163304660923*^9, 3.6541633322975035`*^9}, {
   3.654170005428185*^9, 3.654170008808378*^9}, {3.6542549791457205`*^9, 
   3.654254998275814*^9}, {3.6542551789881506`*^9, 3.6542552015224395`*^9}, {
   3.6542573170614414`*^9, 3.6542573246148734`*^9}, {3.654259782545459*^9, 
   3.654259790045888*^9}, {3.6545989821604805`*^9, 3.654598998998444*^9}, {
   3.6546027011961975`*^9, 3.6546027415275044`*^9}, {3.6546028814985104`*^9, 
   3.6546029013656464`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"A", "[", "m", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"m", ",", "1", ",", "20"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.65407588798094*^9, 3.6540759055459447`*^9}, {
  3.6556327279283485`*^9, 3.6556327281923637`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1.9783708873997186`", ",", "4.388395971065906`", ",", "2.590924905872302`",
    ",", "3.816846060833698`", ",", "2.720709629865935`", ",", 
   "3.649467369906522`", ",", "2.7862706924374003`", ",", 
   "3.563465082302065`", ",", "2.8277977594771744`", ",", 
   "3.509324552993323`", ",", "2.857200406365075`", ",", "3.471377004638419`",
    ",", "2.879462064083754`", ",", "3.4429415244331154`", ",", 
   "2.8970924191940233`", ",", "3.4206385804049284`", ",", 
   "2.911513743238947`", ",", "3.4025553226900063`", ",", 
   "2.923601357913211`", ",", "3.3875188949159516`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.6545990232648315`*^9, 3.654599924211363*^9, 3.654600017478698*^9, 
   3.6546023274148183`*^9, 3.654602797487705*^9, 3.6546028415152235`*^9, 
   3.654602904561829*^9, 3.6546029388707914`*^9, 3.6546033437209477`*^9, 
   3.654607351485179*^9, 3.6546074039011765`*^9, 3.6555631375517616`*^9, 
   3.655631810648883*^9, {3.6556327260382404`*^9, 3.655632736453836*^9}, 
   3.655992321466446*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Solution to the Stefan Problem", "Subsection",
 CellChangeTimes->{{3.658908083377739*^9, 3.6589080967165017`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"sol1", "=", 
   RowBox[{"S", "/.", 
    RowBox[{
     RowBox[{"NDSolve", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"S", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            UnderoverscriptBox["\[Sum]", 
             RowBox[{"n", "=", "1"}], "10"], 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"B", "[", "n", "]"}], 
               RowBox[{"Exp", "[", 
                RowBox[{
                 RowBox[{"-", "t"}], "*", 
                 RowBox[{
                  RowBox[{"g", "[", "n", "]"}], "^", "2"}]}], "]"}], 
               RowBox[{"(", 
                RowBox[{"-", 
                 RowBox[{"BesselJ", "[", 
                  RowBox[{"1", ",", 
                   RowBox[{
                    RowBox[{"g", "[", "n", "]"}], 
                    RowBox[{"S", "[", "t", "]"}]}]}], "]"}]}], ")"}]}], ")"}],
              "0.5"}]}], "-", 
           RowBox[{"Sum", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"A", "[", "m", "]"}], 
              RowBox[{"E", "^", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "t"}], "*", 
                 RowBox[{
                  RowBox[{"\[Beta]", "[", "m", "]"}], "^", "2"}]}], " ", 
                ")"}]}], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", 
                  RowBox[{"BesselJ", "[", 
                   RowBox[{"1", ",", 
                    RowBox[{
                    RowBox[{"\[Beta]", "[", "m", "]"}], 
                    RowBox[{"S", "[", "t", "]"}]}]}], "]"}]}], " ", 
                 RowBox[{"BesselY", "[", 
                  RowBox[{"0", ",", 
                   RowBox[{
                    RowBox[{"\[Beta]", "[", "m", "]"}], 
                    RowBox[{"S", "[", "t", "]"}]}]}], "]"}]}], "+", 
                RowBox[{
                 RowBox[{"BesselJ", "[", 
                  RowBox[{"0", ",", 
                   RowBox[{
                    RowBox[{"\[Beta]", "[", "m", "]"}], 
                    RowBox[{"S", "[", "t", "]"}]}]}], "]"}], " ", 
                 RowBox[{"BesselY", "[", 
                  RowBox[{"1", ",", 
                   RowBox[{
                    RowBox[{"\[Beta]", "[", "m", "]"}], " ", 
                    RowBox[{"S", "[", "t", "]"}]}]}], "]"}]}]}], ")"}], 
              "0.5"}], ",", 
             RowBox[{"{", 
              RowBox[{"m", ",", "1", ",", "10"}], "}"}]}], "]"}]}]}], ",", 
         RowBox[{
          RowBox[{"S", "[", "0.01", "]"}], "\[Equal]", "0.01"}]}], "}"}], ",",
        "S", ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0.01", ",", "3"}], "}"}], ",", 
       RowBox[{"MaxSteps", "\[Rule]", 
        RowBox[{"10", "^", "8"}]}]}], "]"}], "[", 
     RowBox[{"[", 
      RowBox[{"1", ",", "1"}], "]"}], "]"}]}]}], ";"}]], "Text",
 CellChangeTimes->{3.6589081115533504`*^9}],

Cell["\<\
p1 = Plot[Evaluate[Table[sol1[t], {S, 0, 2}]], {t, 0.01, 1}, FrameLabel -> \
{Style[\"Dimensionless time\", 16], Style[ \"Interface Position\", 16]}, \
Frame -> True, PlotStyle -> {Black, Thick}]\
\>", "Text",
 CellChangeTimes->{3.658908124105068*^9}],

Cell["\<\
Using this general code we can get graphs by changing dimensionless \
parameters and values.\
\>", "Text",
 CellChangeTimes->{{3.6589081467163615`*^9, 3.6589082022055354`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[" acknowledgments ", "AcknowledgmentsSection",
 CellChangeTimes->{{3.658908451841814*^9, 3.6589084593902454`*^9}}],

Cell["\<\
I would like to acknowledge B.G Higgens for his useful suggestions.\
\>", "Acknowledgments",
 CellChangeTimes->{{3.6589082161333323`*^9, 3.6589082239547796`*^9}, {
  3.6589083010631895`*^9, 3.6589083172621164`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["REFERENCES", "ReferenceSection",
 CellChangeTimes->{{3.6589084638364997`*^9, 3.658908466723665*^9}}],

Cell["\<\
[1] Wolfram Research, Inc., Mathematica, Version 5.2, Champaign, IL, 2005.
[2] [1]\tB.G Higgens, \[OpenCurlyDoubleQuote]Transient Heat Conduction in a \
Finite Length\[CloseCurlyDoubleQuote], \
sites.google.com/site/chemengwithmathematica/home/ech51-material-balances/\
class-notes.
\
\>", "Text",
 CellChangeTimes->{{3.658908403893071*^9, 3.6589084261443443`*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1904, 943},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->FrontEnd`FileName[{"Article"}, "Preprint.nb", 
  CharacterEncoding -> "WindowsANSI"],
PrivateNotebookOptions -> {"NotebookAuthor" -> ""}
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 184, 6, 83, "Title"],
Cell[CellGroupData[{
Cell[788, 32, 158, 2, 39, "Author"],
Cell[949, 36, 146, 3, 15, "Institution"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1132, 44, 62, 1, 39, "Author"],
Cell[1197, 47, 65, 1, 15, "Institution"],
Cell[CellGroupData[{
Cell[1287, 52, 100, 1, 37, "AbstractSection"],
Cell[1390, 55, 608, 9, 40, "Abstract"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[2047, 70, 158, 5, 70, "Section"],
Cell[CellGroupData[{
Cell[2230, 79, 195, 4, 36, "Subsection"],
Cell[2428, 85, 219, 4, 29, "Text"],
Cell[2650, 91, 583, 18, 23, "Equation"],
Cell[3236, 111, 347, 6, 29, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3620, 122, 111, 1, 36, "Subsection"],
Cell[3734, 125, 332, 9, 119, "Text"],
Cell[4069, 136, 1155, 19, 29, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5261, 160, 135, 3, 36, "Subsection"],
Cell[5399, 165, 1655, 30, 202, "Program"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[7103, 201, 179, 4, 70, "Section"],
Cell[7285, 207, 357, 9, 52, "Input"],
Cell[7645, 218, 5015, 87, 91, "Input"],
Cell[CellGroupData[{
Cell[12685, 309, 238, 6, 31, "Input"],
Cell[12926, 317, 816, 13, 30, "Output"]
}, Open  ]],
Cell[13757, 333, 7287, 160, 154, "Input"],
Cell[CellGroupData[{
Cell[21069, 497, 290, 7, 31, "Input"],
Cell[21362, 506, 1049, 18, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22448, 529, 120, 1, 36, "Subsection"],
Cell[22571, 532, 3037, 81, 84, "Text"],
Cell[25611, 615, 261, 5, 29, "Text"],
Cell[25875, 622, 186, 4, 29, "Text"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[26110, 632, 119, 1, 69, "AcknowledgmentsSection"],
Cell[26232, 635, 225, 4, 29, "Acknowledgments"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26494, 644, 106, 1, 69, "ReferenceSection"],
Cell[26603, 647, 375, 8, 65, "Text"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)