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

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

(* CreatedBy='Wolfram 14.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       154,          7]
NotebookDataLength[     44176,       1204]
NotebookOptionsPosition[     39374,       1130]
NotebookOutlinePosition[     39767,       1146]
CellTagsIndexPosition[     39724,       1143]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
3D flow; solution from Melcher and Schwarz, Physics of Fluids 11,2604 \
(1968);doi:10.1063/1.1691866 \
\>", "Text",
 CellChangeTimes->{{3.9305509570634203`*^9, 3.9305509613024607`*^9}, {
  3.993276519648038*^9, 3.993276520406001*^9}, {3.993276681025387*^9, 
  3.9932766949125137`*^9}, {3.993276744348081*^9, 3.9932767802970133`*^9}},
 FormatType->
  "TextForm",ExpressionUUID->"f4ad87b8-89db-4a0b-95e1-cf736cdfddfa"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"vz1", "[", 
    RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"B1", "+", " ", 
      RowBox[{"C1", " ", "z"}]}], ")"}], " ", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "n"}], " ", "z"}], "]"}], " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vx1", "[", 
    RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "n"}], " ", "z"}]], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"F1", " ", 
         SuperscriptBox["n", "2"]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "q", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"B1", " ", "n"}], "+", 
           RowBox[{"C1", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", 
              RowBox[{"n", " ", "z"}]}], ")"}]}]}], ")"}]}]}], ")"}]}], 
     SuperscriptBox["n", "2"]], 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vy1", "[", 
    RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"\[ImaginaryI]", " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "n"}], " ", "z"}]], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"B1", " ", 
            SuperscriptBox["k", "2"]}], "+", 
           RowBox[{"\[ImaginaryI]", " ", "F1", " ", "n", " ", "q"}]}], 
          ")"}]}], "+", 
        RowBox[{"C1", " ", 
         SuperscriptBox["k", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"n", " ", "z"}]}], ")"}]}]}], ")"}]}], 
     RowBox[{"k", " ", 
      SuperscriptBox["n", "2"]}]], 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"p1", "[", 
     RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], "=", 
    RowBox[{"2", " ", "C1", " ", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"-", "n"}], " ", "z"}]], 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"s", " ", "t"}], "-", 
       RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
       RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vz2", "[", 
   RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"B2", "+", " ", 
     RowBox[{"C2", " ", "z"}]}], ")"}], " ", 
   RowBox[{"Exp", "[", 
    RowBox[{"n", " ", "z"}], "]"}], 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vx2", "[", 
    RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"-", 
     FractionBox[
      RowBox[{"\[ImaginaryI]", " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"n", " ", "z"}]], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"\[ImaginaryI]", " ", "F2", " ", 
          SuperscriptBox["n", "2"]}], "+", 
         RowBox[{"q", " ", 
          RowBox[{"(", 
           RowBox[{"C2", "+", 
            RowBox[{"B2", " ", "n"}], "+", 
            RowBox[{"C2", " ", "n", " ", "z"}]}], ")"}]}]}], ")"}]}], 
      SuperscriptBox["n", "2"]]}], 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vy2", "[", 
    RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"n", " ", "z"}]], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "F2"}], " ", "q"}], "-", 
        FractionBox[
         RowBox[{"\[ImaginaryI]", " ", 
          SuperscriptBox["k", "2"], " ", 
          RowBox[{"(", 
           RowBox[{"C2", "+", 
            RowBox[{"B2", " ", "n"}], "+", 
            RowBox[{"C2", " ", "n", " ", "z"}]}], ")"}]}], 
         SuperscriptBox["n", "2"]]}], ")"}]}], "k"], 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"p2", "[", 
     RowBox[{"z_", ",", "x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
    RowBox[{"2", " ", "C2", " ", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"n", " ", "z"}]], 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"s", " ", "t"}], "-", 
       RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
       RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]]}]}], ";"}], 
  "\[IndentingNewLine]", "\n"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.734944381831294*^9, 3.734944457893009*^9}, {
   3.734945406611662*^9, 3.734945456679431*^9}, {3.7349455195552607`*^9, 
   3.734945580158185*^9}, {3.734945619955332*^9, 3.7349457502224827`*^9}, 
   3.7349459060317383`*^9, {3.7372248711391478`*^9, 3.737224878620694*^9}, {
   3.737225004791161*^9, 3.737225026739876*^9}, {3.737375888579014*^9, 
   3.7373759092180967`*^9}, {3.737377385505518*^9, 3.737377407896426*^9}, 
   3.753370641494532*^9, {3.753541797606564*^9, 3.753541819271274*^9}, {
   3.993153818781959*^9, 3.993153829419084*^9}, {3.993153904242517*^9, 
   3.993153906377001*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"d014ecdd-0fe9-4c45-81b3-9b083ac8104c"],

Cell["HD tractions at z = 0 ", "Text",
 CellChangeTimes->{{3.754161132770631*^9, 3.754161137721384*^9}, {
  3.993276985532673*^9, 
  3.993276997168865*^9}},ExpressionUUID->"77d8c928-0909-4865-90b9-\
d40e39777755"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Tzz1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], "-", 
       RowBox[{"p1", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Txx1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}]}], "-", 
       RowBox[{"p1", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Tyy1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}]}], "-", 
       RowBox[{"p1", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Tzx1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Tzy1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Txy1hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy1", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Tzz2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], "-", 
       RowBox[{"p2", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Txx2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}]}], "-", 
       RowBox[{"p2", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Tyy2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}]}], "-", 
       RowBox[{"p2", "[", 
        RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Tzx2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Tzy2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vz2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "z"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Txy2hd", "=", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vx2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
         "]"}], "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"vy2", "[", 
           RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
         "]"}]}], ")"}], "/.", 
      RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.734946148700568*^9, 3.734946161318375*^9}, {
  3.734954733836834*^9, 3.73495473491609*^9}, {3.7350421781463423`*^9, 
  3.73504226338925*^9}, {3.737126783083868*^9, 3.737126790125691*^9}, {
  3.737385122106275*^9, 3.73738516958351*^9}, {3.738078553874476*^9, 
  3.7380785605668297`*^9}, {3.738078602747737*^9, 3.7380786031981487`*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"144d5637-17b8-463e-89ef-78ba20700fb6"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "2"}], " ", "B1", " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    RowBox[{"s", " ", "t"}], "-", 
    RowBox[{"\[ImaginaryI]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "x"}], "+", 
       RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", "n"}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658828782*^9},
 CellLabel->"Out[9]=",ExpressionUUID->"2655ee4a-5981-484d-9116-de266ba8b719"],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "x"}], "+", 
         RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "q", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"\[ImaginaryI]", " ", "F1", " ", "n"}], "-", 
         RowBox[{"B1", " ", "q"}]}], ")"}]}], "+", 
      RowBox[{"C1", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["n", "2"], "+", 
         SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], 
   SuperscriptBox["n", "2"]]}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658863372*^9},
 CellLabel->"Out[10]=",ExpressionUUID->"2155c204-3a39-408b-b45c-ac3a99aa0c58"],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "x"}], "+", 
         RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"C1", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["k", "2"], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "+", 
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "B1"}], " ", 
          SuperscriptBox["k", "2"]}], "-", 
         RowBox[{"\[ImaginaryI]", " ", "F1", " ", "n", " ", "q"}]}], 
        ")"}]}]}], ")"}]}], 
   SuperscriptBox["n", "2"]]}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658894821*^9},
 CellLabel->"Out[11]=",ExpressionUUID->"3284adeb-b5f8-4ac5-9d0a-8728b4780eb0"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "F1"}], " ", 
      SuperscriptBox["n", "2"]}], "+", 
     RowBox[{"2", " ", "\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{"C1", "-", 
        RowBox[{"B1", " ", "n"}]}], ")"}], " ", "q"}]}], ")"}]}], 
  "n"]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658925027*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"81c43482-c369-4ea3-aed3-d7097986663a"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", "\[ImaginaryI]", " ", 
      SuperscriptBox["k", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"C1", "-", 
        RowBox[{"B1", " ", "n"}]}], ")"}]}], "+", 
     RowBox[{"F1", " ", 
      SuperscriptBox["n", "2"], " ", "q"}]}], ")"}]}], 
  RowBox[{"k", " ", "n"}]]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.99327865895411*^9},
 CellLabel->"Out[13]=",ExpressionUUID->"f93427ab-17bc-4ed0-9140-9879e80c1d43"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "2"}], " ", 
      SuperscriptBox["k", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"C1", "-", 
        RowBox[{"B1", " ", "n"}]}], ")"}], " ", "q"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "F1", " ", 
      SuperscriptBox["n", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"k", "-", "q"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{"k", "+", "q"}], ")"}]}]}], ")"}]}], 
  RowBox[{"k", " ", 
   SuperscriptBox["n", "2"]}]]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658990033*^9},
 CellLabel->"Out[14]=",ExpressionUUID->"2e104a09-6122-4088-9a56-fbc1ac906972"],

Cell[BoxData[
 RowBox[{"2", " ", "B2", " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    RowBox[{"s", " ", "t"}], "-", 
    RowBox[{"\[ImaginaryI]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "x"}], "+", 
       RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", "n"}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278658996139*^9},
 CellLabel->"Out[15]=",ExpressionUUID->"af823b5f-5a24-4df9-8d19-16cfd3705c69"],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "x"}], "+", 
         RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "q", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"\[ImaginaryI]", " ", "F2", " ", "n"}], "+", 
         RowBox[{"B2", " ", "q"}]}], ")"}]}], "+", 
      RowBox[{"C2", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["n", "2"], "+", 
         SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], 
   SuperscriptBox["n", "2"]]}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278659021353*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"eab5b21c-c65a-462e-9f6d-426aeb3716c7"],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"2", " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "x"}], "+", 
         RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"C2", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["k", "2"], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "+", 
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"B2", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         RowBox[{"\[ImaginaryI]", " ", "F2", " ", "n", " ", "q"}]}], 
        ")"}]}]}], ")"}]}], 
   SuperscriptBox["n", "2"]]}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278659049836*^9},
 CellLabel->"Out[17]=",ExpressionUUID->"e5060396-3dc6-44fc-b840-eebca675a00a"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"F2", " ", 
      SuperscriptBox["n", "2"]}], "-", 
     RowBox[{"2", " ", "\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{"C2", "+", 
        RowBox[{"B2", " ", "n"}]}], ")"}], " ", "q"}]}], ")"}]}], 
  "n"]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278659072143*^9},
 CellLabel->"Out[18]=",ExpressionUUID->"9531ebdb-1f93-4525-bd39-a53dbc2b99b5"],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"\[ImaginaryI]", " ", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{
      RowBox[{"s", " ", "t"}], "-", 
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "x"}], "+", 
         RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", " ", 
       SuperscriptBox["k", "2"], " ", 
       RowBox[{"(", 
        RowBox[{"C2", "+", 
         RowBox[{"B2", " ", "n"}]}], ")"}]}], "-", 
      RowBox[{"\[ImaginaryI]", " ", "F2", " ", 
       SuperscriptBox["n", "2"], " ", "q"}]}], ")"}]}], 
   RowBox[{"k", " ", "n"}]]}]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278659100403*^9},
 CellLabel->"Out[19]=",ExpressionUUID->"0ca7ea9e-da3d-484f-babc-2443f1d4b07e"],

Cell[BoxData[
 FractionBox[
  RowBox[{"\[ImaginaryI]", " ", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", "\[ImaginaryI]", " ", 
      SuperscriptBox["k", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"C2", "+", 
        RowBox[{"B2", " ", "n"}]}], ")"}], " ", "q"}], "+", 
     RowBox[{"F2", " ", 
      SuperscriptBox["n", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], 
  RowBox[{"k", " ", 
   SuperscriptBox["n", "2"]}]]], "Output",
 CellChangeTimes->{3.9932770085417833`*^9, 3.993278659139118*^9},
 CellLabel->"Out[20]=",ExpressionUUID->"fb8e4089-334e-4436-9a24-b13b7feb8dcf"]
}, Open  ]],

Cell["\<\
for simplicity consider the 2D case only - zx plane , membrane z=h[x] \
fluctuates about the  planar state z=0.  F1 and F2 are rotational solutions \
with no interface deformation.
\
\>", "Text",
 CellChangeTimes->{{3.9931526275914507`*^9, 3.99315263481493*^9}, {
  3.9931533799706287`*^9, 3.993153388985091*^9}, {3.993276814605913*^9, 
  3.993276857970573*^9}, {3.993276897178973*^9, 3.993276926825905*^9}, {
  3.993277467679253*^9, 
  3.993277485133049*^9}},ExpressionUUID->"794c5fae-86af-4abd-95b1-\
0f458ac0b812"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{
  3.993278556535132*^9},ExpressionUUID->"15cf0a66-61ca-42f9-a8b8-\
621fca3d947b"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", " ", "=", " ", "q"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"F1", " ", "=", " ", "0"}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{"F2", " ", "=", " ", "0"}], ";"}]}], "Input",
 CellChangeTimes->{{3.993276929252819*^9, 3.9932769292560167`*^9}},
 CellLabel->"In[21]:=",ExpressionUUID->"30c0eb2e-8f8e-4cee-83f2-da662636d79f"],

Cell["\<\
 BC 1 and BC 2: velocity continuity z- normal velocity, x -tangential velocity\
\>", "Text",
 CellChangeTimes->{{3.7541617037858887`*^9, 3.75416171724622*^9}, 
   3.754161768502489*^9, {3.99327692043736*^9, 3.9932769240130577`*^9}, {
   3.993277054725178*^9, 
   3.9932770795715303`*^9}},ExpressionUUID->"82039b8e-582a-4976-95d1-\
92fbaaa27d9d"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"eqVz", "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"vz1", "[", 
       RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], "-", 
      RowBox[{"vz2", "[", 
       RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], ")"}], "/.", 
    RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"eqVx", "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"vx1", "[", 
       RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], "-", 
      RowBox[{"vx2", "[", 
       RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}]}], ")"}], "/.", 
    RowBox[{"z", "\[Rule]", "0"}]}], "]"}]}], "\[IndentingNewLine]"}], "Input",\

 CellChangeTimes->{
  3.480516376018009*^9, {3.480516512332141*^9, 3.480516522202848*^9}, {
   3.481396795653274*^9, 3.4813968073576813`*^9}, {3.683295478729539*^9, 
   3.6832954829356537`*^9}, {3.734664476838166*^9, 3.734664485965081*^9}, {
   3.734664546823557*^9, 3.734664554460435*^9}, {3.736085369330634*^9, 
   3.736085382154909*^9}, 3.7360855877639837`*^9, {3.754160174664154*^9, 
   3.7541601754493227`*^9}, {3.7541604656355467`*^9, 
   3.7541604660451307`*^9}, {3.754160869258175*^9, 3.754160869457418*^9}, {
   3.754161723706852*^9, 3.754161822720517*^9}, 3.993276944364698*^9},
 CellLabel->"In[24]:=",ExpressionUUID->"667d2e43-0210-43ac-95fc-5673eb7fca00"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{"B1", "-", "B2"}], ")"}], " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    RowBox[{"s", " ", "t"}], "-", 
    RowBox[{"\[ImaginaryI]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "x"}], "+", 
       RowBox[{"k", " ", "y"}]}], ")"}]}]}]]}]], "Output",
 CellChangeTimes->{3.930550250806239*^9, 3.993152637448516*^9, 
  3.993153460437995*^9, 3.993276968065548*^9, 3.993277013518012*^9, 
  3.993278663235083*^9},
 CellLabel->"Out[24]=",ExpressionUUID->"821f285a-fe28-4a05-8524-6985721907e5"],

Cell[BoxData[
 FractionBox[
  RowBox[{"\[ImaginaryI]", " ", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "x"}], "+", 
        RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "C1"}], "+", "C2", "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"B1", "+", "B2"}], ")"}], " ", "q"}]}], ")"}]}], 
  "q"]], "Output",
 CellChangeTimes->{3.930550250806239*^9, 3.993152637448516*^9, 
  3.993153460437995*^9, 3.993276968065548*^9, 3.993277013518012*^9, 
  3.9932786632655983`*^9},
 CellLabel->"Out[25]=",ExpressionUUID->"bc141bed-d128-4179-8adf-726f56cf2d6f"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"eqVx", "\[Equal]", "0"}], ",", " ", 
      RowBox[{"eqVz", "\[Equal]", "0"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"B2", ",", "C2"}], "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.993152664565734*^9, 3.9931526678649273`*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"e06cbfd1-897f-450f-8ed9-95bb78db5e73"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"B2", "\[Rule]", "B1"}], ",", 
    RowBox[{"C2", "\[Rule]", 
     RowBox[{"C1", "-", 
      RowBox[{"2", " ", "B1", " ", "q"}]}]}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.993152668904715*^9, 3.993153479272196*^9, 
  3.9932770175156593`*^9, 3.9932786650627193`*^9},
 CellLabel->"Out[26]=",ExpressionUUID->"930d14a1-7539-471a-8129-db05d9ee068c"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"B2", "=", "B1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"C2", "=", 
   RowBox[{"C1", "-", 
    RowBox[{"2", " ", "B1", " ", "n"}]}]}], 
  ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.754161862172842*^9, 3.754161877542512*^9}, {
  3.930550278760919*^9, 3.9305502805962257`*^9}, {3.993277034931377*^9, 
  3.993277038347972*^9}},
 CellLabel->"In[27]:=",ExpressionUUID->"dd740573-23fc-4ab6-978b-d334eda628a2"],

Cell["BC 3 area incompressibility ", "Text",
 CellChangeTimes->{{3.993152922262754*^9, 3.993152929735931*^9}, {
  3.993277047021332*^9, 3.993277048618205*^9}, {3.993277083792521*^9, 
  3.993277095101892*^9}},ExpressionUUID->"3f177db1-0807-4205-851a-\
2518cd8e60e6"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"vx1", "[", 
          RowBox[{"z", ",", "x", ",", "y", ",", "t"}], "]"}], ",", "x"}], 
        "]"}], "/.", 
       RowBox[{"z", "->", "0"}]}], ")"}], "==", "0"}], ",", " ", "C1"}], 
   "]"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.993152931405468*^9, 3.993152969908237*^9}, {
  3.9931530206087093`*^9, 3.993153047284206*^9}, {3.993277100532433*^9, 
  3.993277106809614*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"a2b85397-2795-4e72-a949-8ab8e13fbb8a"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"C1", "\[Rule]", 
    RowBox[{"B1", " ", "q"}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.99315294897042*^9, 3.9931529704246407`*^9}, {
   3.993153038462901*^9, 3.9931530477654753`*^9}, 3.993153499148468*^9, 
   3.9932771086924667`*^9},
 CellLabel->"Out[29]=",ExpressionUUID->"00115711-ff0f-4fcc-bb0b-82524ea00793"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"C1", "=", 
   RowBox[{"B1", " ", "q"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.993277115221258*^9, 3.993277117348399*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"324cbd60-b070-4dd6-8756-131625050405"],

Cell["\<\
BC 4 : tractions balance,  planar interface with normal in z direction 
normal traction\
\>", "Text",
 CellChangeTimes->{{3.735041924605235*^9, 3.7350419382189283`*^9}, {
   3.7350419717958927`*^9, 3.735042039720652*^9}, {3.735042135440206*^9, 
   3.735042167894506*^9}, {3.735042235019622*^9, 3.7350422362037354`*^9}, {
   3.735042269765376*^9, 3.735042437133864*^9}, {3.73504246741686*^9, 
   3.7350425383847437`*^9}, {3.735042571821682*^9, 3.7350425805731564`*^9}, {
   3.7350426684153023`*^9, 3.735042673363592*^9}, {3.73504272040926*^9, 
   3.735042729160207*^9}, {3.7350433117733727`*^9, 3.7350433207358503`*^9}, {
   3.735043455266927*^9, 3.735043520164506*^9}, {3.7360866869822807`*^9, 
   3.736086714545104*^9}, 3.7360867851951733`*^9, {3.737385247087268*^9, 
   3.7373852677256927`*^9}, {3.7373872156876793`*^9, 3.737387230710207*^9}, {
   3.7380788028203278`*^9, 3.738078825590032*^9}, {3.7533725107690563`*^9, 
   3.7533725639615917`*^9}, {3.754161367375412*^9, 3.7541614441169777`*^9}, {
   3.755006568873129*^9, 3.7550065731095457`*^9}, {3.7550066638971033`*^9, 
   3.755006664225905*^9}, {3.930550384233891*^9, 3.930550395701104*^9}, {
   3.993277119841569*^9, 3.993277143948267*^9}, {3.993277276869624*^9, 
   3.993277291326063*^9}, {3.993278616577636*^9, 
   3.993278617636589*^9}},ExpressionUUID->"1b02b2cd-5e5e-411e-aa70-\
0c3665029878"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{"tHDnZ", " ", "=", " ", 
    RowBox[{"Simplify", "[", 
     RowBox[{
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{"\[Eta]1", " ", "Tzz1hd"}], " ", "-", " ", 
        RowBox[{"\[Eta]2", " ", "Tzz2hd"}]}], ")"}], " ", "/.", " ", 
      RowBox[{"z", " ", "->", " ", "0"}]}], "]"}]}], "\n"}]}]], "Input",
 CellChangeTimes->{{3.993277146434667*^9, 3.99327717322436*^9}, {
  3.9932772820680447`*^9, 3.993277296839823*^9}, {3.993278629917941*^9, 
  3.993278637466114*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"2e175bf6-df95-4eb4-811c-1e4328e3dc27"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "2"}], " ", "B1", " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    RowBox[{"s", " ", "t"}], "-", 
    RowBox[{"\[ImaginaryI]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "x"}], "+", 
       RowBox[{"k", " ", "y"}]}], ")"}]}]}]], " ", "q", " ", 
  RowBox[{"(", 
   RowBox[{"\[Eta]1", "+", "\[Eta]2"}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.993153533469091*^9, {3.99327715352577*^9, 3.993277175746769*^9}, 
   3.9932772947293386`*^9, 3.993278671087823*^9},
 CellLabel->"Out[30]=",ExpressionUUID->"44f59191-5a58-4b38-88bd-8653792d23ba"]
}, Open  ]],

Cell["membrane normal stress : (* H=-d^2h/dx^2 *)", "Text",
 CellChangeTimes->{{3.993277259683009*^9, 3.9932772686693773`*^9}, 
   3.993277304510478*^9},ExpressionUUID->"d2cd57de-a3b9-4033-a8b9-\
1c622c23ba59"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"tmem", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "\[Kappa]"}], " ", 
    RowBox[{"(", 
     RowBox[{"-", 
      RowBox[{"D", "[", 
       RowBox[{
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{
          RowBox[{"s", " ", "t"}], "-", 
          RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
          RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]], ",", " ", 
        RowBox[{"{", 
         RowBox[{"x", ",", "4"}], "}"}]}], "]"}]}], ")"}]}], " ", "-", 
   RowBox[{"\[Sigma]", " ", 
    RowBox[{"D", "[", 
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"s", " ", "t"}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]], ",", " ", 
      RowBox[{"{", 
       RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.993153561329705*^9, 3.993153571810046*^9}, {
  3.993153707076737*^9, 3.9931537168775253`*^9}, {3.993155041442789*^9, 
  3.993155062226645*^9}, {3.993277243376012*^9, 3.993277270436764*^9}, {
  3.993278592368479*^9, 3.993278610031033*^9}},
 CellLabel->"In[31]:=",ExpressionUUID->"09b83750-1f45-48c6-9159-7664cb9aa847"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]], " ", 
   SuperscriptBox["q", "4"], " ", "\[Kappa]"}], "+", 
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"s", " ", "t"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "q", " ", "x"}], "-", 
     RowBox[{"\[ImaginaryI]", " ", "k", " ", "y"}]}]], " ", 
   SuperscriptBox["q", "2"], " ", "\[Sigma]"}]}]], "Output",
 CellChangeTimes->{3.9931535722643003`*^9, 3.9931537173450747`*^9, 
  3.993155046604031*^9, 3.9932773090720654`*^9, 3.9932786738526573`*^9},
 CellLabel->"Out[31]=",ExpressionUUID->"4c71351d-83be-49b9-ab66-3ebc7aede4e8"]
}, Open  ]],

Cell["\<\
find the normal velocity amplitude, this is the relaxtaion rate s\
\>", "Text",
 CellChangeTimes->{{3.993277346881728*^9, 3.993277383547002*^9}, {
  3.993277494541038*^9, 
  3.99327749488457*^9}},ExpressionUUID->"532f3af4-3a6a-48ac-8530-\
6e589ebc2a4b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sol", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"tHDnZ", "==", "tmem"}], ",", " ", "B1"}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.993277311691657*^9, 3.99327733694508*^9}, {
  3.993277508548937*^9, 3.9932775093166103`*^9}, {3.993278688294367*^9, 
  3.993278692699733*^9}},
 CellLabel->"In[34]:=",ExpressionUUID->"8ea5f504-f47a-4318-b12f-a13ed046c534"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"B1", "\[Rule]", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{"q", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SuperscriptBox["q", "2"], " ", "\[Kappa]"}], "+", "\[Sigma]"}], 
        ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"\[Eta]1", "+", "\[Eta]2"}], ")"}]}]]}]}], "}"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.9932773272084417`*^9, 3.993277337390935*^9}, 
   3.993277509874714*^9, {3.993278676964429*^9, 3.993278693112803*^9}},
 CellLabel->"Out[34]=",ExpressionUUID->"ad3fb18c-6126-414f-ab4f-db1d7a9e46f6"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.9932775039871283`*^9, 
  3.993277506508168*^9}},ExpressionUUID->"15ced8d9-f358-4780-b9b9-\
4f8252cc583a"],

Cell["kinematic condition", "Text",
 CellChangeTimes->{{3.993277387898533*^9, 
  3.99327741054736*^9}},ExpressionUUID->"3dd3793d-74c8-4ba8-be1f-\
326674267310"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"s", "=", 
  RowBox[{"B1", "/.", "sol"}]}]], "Input",
 CellChangeTimes->{{3.993277398081143*^9, 3.993277404040779*^9}, {
  3.9932775117796392`*^9, 3.9932775147705307`*^9}},
 CellLabel->"In[33]:=",ExpressionUUID->"f17be489-f526-4b92-b618-f48af60f970a"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"-", 
   FractionBox[
    RowBox[{"q", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        SuperscriptBox["q", "2"], " ", "\[Kappa]"}], "+", "\[Sigma]"}], 
      ")"}]}], 
    RowBox[{"2", " ", 
     RowBox[{"(", 
      RowBox[{"\[Eta]1", "+", "\[Eta]2"}], ")"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{{3.993277497772623*^9, 3.993277515174111*^9}, 
   3.993278678888824*^9},
 CellLabel->"Out[33]=",ExpressionUUID->"e6be9f0e-57ef-463b-b51b-c77439fbec67"]
}, Open  ]],

Cell["\<\
 Note : the tangential traction is not used in this case, for a membrane it \
serves to determine the tension gradients\
\>", "Text",
 CellChangeTimes->{{3.9932771870293503`*^9, 3.993277239908988*^9}, {
  3.993277415721508*^9, 
  3.993277418872403*^9}},ExpressionUUID->"645fbd0f-4ff9-4c38-b819-\
ced17fce600d"],

Cell[BoxData[
 RowBox[{"tHDtX", " ", "=", " ", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"(", " ", 
     RowBox[{"Tzx1hd", " ", "-", " ", 
      RowBox[{"lam", " ", "Tzx2hd"}]}], ")"}], " ", "/.", " ", 
    RowBox[{"z", " ", "->", " ", "0"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.993277236089229*^9, 
  3.993277236091898*^9}},ExpressionUUID->"f37d088a-e47e-4d8f-86a9-\
470e2e198581"]
},
WindowSize->{808, 732},
WindowMargins->{{155, Automatic}, {25, Automatic}},
FrontEndVersion->"14.3 for Mac OS X ARM (64-bit) (July 8, 2025)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"5c6208e1-19d6-48c0-a912-23b1732fb4ee"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[554, 20, 426, 8, 33, "Text",ExpressionUUID->"f4ad87b8-89db-4a0b-95e1-cf736cdfddfa"],
Cell[983, 30, 6625, 181, 378, "Input",ExpressionUUID->"d014ecdd-0fe9-4c45-81b3-9b083ac8104c"],
Cell[7611, 213, 213, 4, 33, "Text",ExpressionUUID->"77d8c928-0909-4865-90b9-d40e39777755"],
Cell[CellGroupData[{
Cell[7849, 221, 6684, 185, 437, "Input",ExpressionUUID->"144d5637-17b8-463e-89ef-78ba20700fb6"],
Cell[14536, 408, 466, 12, 33, "Output",ExpressionUUID->"2655ee4a-5981-484d-9116-de266ba8b719"],
Cell[15005, 422, 895, 26, 52, "Output",ExpressionUUID->"2155c204-3a39-408b-b45c-ac3a99aa0c58"],
Cell[15903, 450, 962, 29, 52, "Output",ExpressionUUID->"3284adeb-b5f8-4ac5-9d0a-8728b4780eb0"],
Cell[16868, 481, 730, 22, 51, "Output",ExpressionUUID->"81c43482-c369-4ea3-aed3-d7097986663a"],
Cell[17601, 505, 765, 22, 51, "Output",ExpressionUUID->"f93427ab-17bc-4ed0-9140-9879e80c1d43"],
Cell[18369, 529, 943, 28, 52, "Output",ExpressionUUID->"2e104a09-6122-4088-9a56-fbc1ac906972"],
Cell[19315, 559, 449, 11, 33, "Output",ExpressionUUID->"af823b5f-5a24-4df9-8d19-16cfd3705c69"],
Cell[19767, 572, 895, 26, 52, "Output",ExpressionUUID->"eab5b21c-c65a-462e-9f6d-426aeb3716c7"],
Cell[20665, 600, 936, 28, 52, "Output",ExpressionUUID->"e5060396-3dc6-44fc-b840-eebca675a00a"],
Cell[21604, 630, 708, 21, 51, "Output",ExpressionUUID->"9531ebdb-1f93-4525-bd39-a53dbc2b99b5"],
Cell[22315, 653, 825, 23, 51, "Output",ExpressionUUID->"0ca7ea9e-da3d-484f-babc-2443f1d4b07e"],
Cell[23143, 678, 963, 28, 52, "Output",ExpressionUUID->"fb8e4089-334e-4436-9a24-b13b7feb8dcf"]
}, Open  ]],
Cell[24121, 709, 527, 11, 77, "Text",ExpressionUUID->"794c5fae-86af-4abd-95b1-0f458ac0b812"],
Cell[24651, 722, 128, 3, 29, "Input",ExpressionUUID->"15cf0a66-61ca-42f9-a8b8-621fca3d947b"],
Cell[24782, 727, 369, 9, 70, "Input",ExpressionUUID->"30c0eb2e-8f8e-4cee-83f2-da662636d79f"],
Cell[25154, 738, 355, 7, 33, "Text",ExpressionUUID->"82039b8e-582a-4976-95d1-92fbaaa27d9d"],
Cell[CellGroupData[{
Cell[25534, 749, 1420, 31, 70, "Input",ExpressionUUID->"667d2e43-0210-43ac-95fc-5673eb7fca00"],
Cell[26957, 782, 569, 15, 33, "Output",ExpressionUUID->"821f285a-fe28-4a05-8524-6985721907e5"],
Cell[27529, 799, 741, 21, 51, "Output",ExpressionUUID->"bc141bed-d128-4179-8adf-726f56cf2d6f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28307, 825, 447, 11, 29, "Input",ExpressionUUID->"e06cbfd1-897f-450f-8ed9-95bb78db5e73"],
Cell[28757, 838, 421, 10, 33, "Output",ExpressionUUID->"930d14a1-7539-471a-8129-db05d9ee068c"]
}, Open  ]],
Cell[29193, 851, 472, 11, 70, "Input",ExpressionUUID->"dd740573-23fc-4ab6-978b-d334eda628a2"],
Cell[29668, 864, 265, 4, 33, "Text",ExpressionUUID->"3f177db1-0807-4205-851a-2518cd8e60e6"],
Cell[CellGroupData[{
Cell[29958, 872, 641, 17, 49, "Input",ExpressionUUID->"a2b85397-2795-4e72-a949-8ab8e13fbb8a"],
Cell[30602, 891, 378, 8, 33, "Output",ExpressionUUID->"00115711-ff0f-4fcc-bb0b-82524ea00793"]
}, Open  ]],
Cell[30995, 902, 238, 5, 29, "Input",ExpressionUUID->"324cbd60-b070-4dd6-8756-131625050405"],
Cell[31236, 909, 1366, 21, 55, "Text",ExpressionUUID->"1b02b2cd-5e5e-411e-aa70-0c3665029878"],
Cell[CellGroupData[{
Cell[32627, 934, 608, 14, 70, "Input",ExpressionUUID->"2e175bf6-df95-4eb4-811c-1e4328e3dc27"],
Cell[33238, 950, 610, 16, 33, "Output",ExpressionUUID->"44f59191-5a58-4b38-88bd-8653792d23ba"]
}, Open  ]],
Cell[33863, 969, 210, 3, 33, "Text",ExpressionUUID->"d2cd57de-a3b9-4033-a8b9-1c622c23ba59"],
Cell[CellGroupData[{
Cell[34098, 976, 1220, 30, 31, "Input",ExpressionUUID->"09b83750-1f45-48c6-9159-7664cb9aa847"],
Cell[35321, 1008, 792, 18, 33, "Output",ExpressionUUID->"4c71351d-83be-49b9-ab66-3ebc7aede4e8"]
}, Open  ]],
Cell[36128, 1029, 263, 6, 33, "Text",ExpressionUUID->"532f3af4-3a6a-48ac-8530-6e589ebc2a4b"],
Cell[CellGroupData[{
Cell[36416, 1039, 432, 10, 29, "Input",ExpressionUUID->"8ea5f504-f47a-4318-b12f-a13ed046c534"],
Cell[36851, 1051, 635, 18, 53, "Output",ExpressionUUID->"ad3fb18c-6126-414f-ab4f-db1d7a9e46f6"]
}, Open  ]],
Cell[37501, 1072, 154, 3, 29, "Input",ExpressionUUID->"15ced8d9-f358-4780-b9b9-4f8252cc583a"],
Cell[37658, 1077, 160, 3, 33, "Text",ExpressionUUID->"3dd3793d-74c8-4ba8-be1f-326674267310"],
Cell[CellGroupData[{
Cell[37843, 1084, 274, 5, 29, "Input",ExpressionUUID->"f17be489-f526-4b92-b618-f48af60f970a"],
Cell[38120, 1091, 508, 15, 53, "Output",ExpressionUUID->"e6be9f0e-57ef-463b-b51b-c77439fbec67"]
}, Open  ]],
Cell[38643, 1109, 320, 7, 55, "Text",ExpressionUUID->"645fbd0f-4ff9-4c38-b819-ced17fce600d"],
Cell[38966, 1118, 404, 10, 29, "Input",ExpressionUUID->"f37d088a-e47e-4d8f-86a9-470e2e198581"]
}
]
*)