Skip to content
Snippets Groups Projects
Commit 4d93c2bd authored by Nitschke, Ingo's avatar Nitschke, Ingo
Browse files

2024-01-24-16-56-14

parent 0cd58201
No related branches found
No related tags found
No related merge requests found
No preview for this file type
......@@ -169,7 +169,8 @@ For this purpose we define the flux potential
&:= \frac{\cAC}{2}\int_{\surf} \Tr(\Dlow\IbxiQ)\dS \label{eq:fluxAC_origin}
\end{align}
where $ \cAC\in\R $.
Since the Jaumann derivative is $ \R^3 $-metric compatible, \ie\ $ \Djau\Id=\nullb $, it holds $ \Dlow\IbxiQ = -\xi\Djau\Qb + \Sb[\Vb]\IbxiQ + \IbxiQ\Sb[\Vb] $ by \eqref{eq:Dlow_rel_to_Djau}.
Since the Jaumann derivative is $ \R^3 $-metric compatible, \ie\ $ \Djau\Id=\nullb $, it holds
$ \Dlow\IbxiQ = -\xi\Djau\Qb + \Sb[\Vb]\IbxiQ + \IbxiQ\Sb[\Vb] $ by \eqref{eq:Dlow_rel_to_Djau}. %MProved
The space of Q-tensor fields $ \tangentQR $ is closed by the Jaumann derivative \cite{NitschkeVoigt_2023}, \ie\ $ \Tr(\Djau\Qb) = 0 $ particularly.
Therefore the flux potential \eqref{eq:fluxAC_origin} results in
\begin{align}
......@@ -177,7 +178,7 @@ Therefore the flux potential \eqref{eq:fluxAC_origin} results in
&= \cAC \innerH{\tangentSymR}{\IbxiQ,\Sb[\Vb]}
= \cAC \innerH{\tangentR[^2]}{\IdS - \xi \IdS\Qb,\nablaC\Vb} \formPeriod
\label{eq:fluxAC}
\end{align}
\end{align}%MProved
Contrarily to the genuine flux potentials in \cite{Nitschke2023}, this potential is neither quadratic, nor generally positive at all.
Moreover it transfers one-to-one into energy exchange with an sink/source, \eg\ the surrounded environment, as we see in Section \ref{sec:energy_rate}.
Therefore, \eqref{eq:fluxAC} does not stipulate a passive mechanism, hence we call it the active flux potential.
......@@ -205,7 +206,7 @@ For a better classification we decompose \eqref{eq:fluxAC} into an isotropic and
\end{align}
with $ \cIA=\cAC $, $ \cNA = -\xi\cAC $, and $\energyAC = \energyIA + \energyNA$.
The active isotropic flux potential $ \energyIA $ solely depends on the surface extensibility $ \DivC\Vb $,
and the active nematic flux potential $ \energyNA $ on the deviatoric stress $ \projQR\Sb[\Vb] = \frac{1}{2}(\nablaC\Vb + \nablaC^T\Vb) - \frac{\DivC\Vb}{3}\Id $ in direction of $ \Qb $.
and the active nematic flux potential $ \energyNA $ on the deviatoric stress $ \projQR\Sb[\Vb] = \Sb[\Vb] - \frac{\DivC\Vb}{3}\Id $ in direction of $ \Qb $.
Hence the active isotropic flux potential vanishes for inextensible fluids.
However, it cannot be neglected if one ensure $ \DivC\Vb=0 $ by the Lagrange-multiplier method, as we did in \cite{Nitschke2023}.
Eventually, Section \ref{sec:active_forces} leads to an additional pressure-like force, which has not any impact on divergence-free solutions $ \Vb $.
......@@ -220,14 +221,41 @@ The variations $ \innerH{\tangentR}{\deltafrac{\energyIA}{\Vb},\Wb} =: -\innerH{
&= \cIA\DivC\IdS
= \cIA\GradC 1
= \cIA\meanc\normal \formComma \\
\FbNV
\FbNV \label{eq:active_nematic_force}
&= \cNA\DivC\left( \IdS\Qb\IdS \right) \formPeriod
\end{align}
\rednote{
\begin{itemize}
\item The isotropic force field, \resp\ its stress tensor field $ \cIA\IdS $, corresponds to the isotropic part of the ansatz in \cite{Al-IzziAlexander_PRR_2023}
\end{itemize}
}
The active flux potential does not depend on any Q-tensor rate, which could serves as a process variable associated to $ \Qb $.
As a consequence no active molecular forces $ \HbAC\in\tangentQR $ occurs in the Q-tensor equation, \resp\ we stipulate $ \Hb_{\IA} = \Hb_{\NV} = 0 $ in compliance with the setup in \cite{Nitschke2023}.
\subsubsection{Surface Conforming Q-Tensor}
Any Q-tensor $ \Qb\in\tangentQR $ field is uniquely orthogonal decomposable into a tangential (flat degenerated) Q-tensor field $ \qb\in\tangentQS $,
a scalar-valued normal eigenvector field $ \beta\in\tangentS[^0] $,
and a tangential vector-valued surface non-conforming field $ \etab\in\tangentS $, by
\begin{align}
\Qb &= \qb + \etab\otimes\normal + \normal\otimes\etab + \beta\left( \normal\otimes\normal -\frac{1}{2}\IdS \right) \formComma
\label{eq:qtensor_decomposition}
\end{align}
see \cite{NitschkeVoigt_2023}.
Substituting this decomposition into the active nematic force field \eqref{eq:active_nematic_force} yields
\begin{align*}
\FbNV
&= \cNA\left( \DivC\qb - \frac{1}{2}\GradC\beta \right)
= \cNA\left(\div\qb - \frac{1}{2}\nabla\beta + \left( \qb\dbdot\shop - \frac{\meanc}{2}\beta \right)\normal\right) \formComma
\end{align*}
which does not depend on any surface non-conforming parts of the Q-tensor field.
Therefore it holds $ \FbNV|_{\Qb\in\tangentQR} = \FbNV|_{\Qb\in\tangentCQR} $ \wrt\ the space of surface conforming Q-tensor fields
\begin{align*}
\tangentCQR
&:= \left\{ \Qb \in \tangentQR \mid \exists\lambda\in\tangentS[^0] : \Qb\normal = \lambda\normal \right\} \formPeriod
\end{align*}
\subsection{Active Beris-Edwards Models}
......
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 13.2' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 16875, 444]
NotebookOptionsPosition[ 14255, 392]
NotebookOutlinePosition[ 14677, 409]
CellTagsIndexPosition[ 14634, 406]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Once", "@",
RowBox[{"PrependTo", "[",
RowBox[{"$Path", ",", "\"\</home/ingo/git/surfsup\>\""}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"<<", "\"\<SurfsUp`\>\""}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"?", "\"\<SurfsUp`*\>\""}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Once", "@",
RowBox[{"LaunchKernels", "[", "8", "]"}]}], ";"}],
RowBox[{"(*",
RowBox[{"optional", ";", " ",
RowBox[{
"on", " ", "my", " ", "16", " ", "kernel", " ", "machine", " ", "it", " ",
"take", " ", "8", " ", "kernel"}]}], "*)"}]}], "\[IndentingNewLine]",
RowBox[{"$KernelCount", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Lagrange", " ", "Perspective"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{"InitMovingParameterization", "[",
RowBox[{
RowBox[{"{",
RowBox[{"y1", ",", "y2", ",",
RowBox[{"hf", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"y1", ",", "y2"}], "}"}]}], "]"}], ";"}],
"*)"}]}], "\[IndentingNewLine]",
RowBox[{"InitMovingParameterization", "[",
RowBox[{
RowBox[{"{",
RowBox[{"y1", ",", "y2", ",",
RowBox[{"hf", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"y1", ",", "y2"}], "}"}], ",", "t", ",",
RowBox[{"{",
RowBox[{
RowBox[{"vm1", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"vm2", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "}"}], ",", "\"\<S\>\""}],
"]"}]}], "Input",
CellChangeTimes->{{3.9120878485289516`*^9, 3.912087873391719*^9}, {
3.9120885662154903`*^9, 3.912088671128661*^9}, {3.912088823712516*^9,
3.912088843751397*^9}, {3.91208889556355*^9, 3.912088895775659*^9}, {
3.9120889641333847`*^9, 3.912088968687875*^9}, {3.912146071803026*^9,
3.912146085055224*^9}, {3.9121462265943623`*^9, 3.91214626151243*^9}, {
3.912146444785282*^9, 3.912146451409123*^9}, 3.9121465023360577`*^9, {
3.912166952414016*^9, 3.912166965550562*^9}, {3.912771282968028*^9,
3.912771312441022*^9}, {3.912771447127698*^9, 3.912771455780593*^9},
3.912771903197764*^9, 3.9129368097083387`*^9, {3.913282961111938*^9,
3.913282964139681*^9}, {3.9136979029913483`*^9, 3.913697924122053*^9}, {
3.913967512821052*^9, 3.9139675139075937`*^9}, {3.9139692792809887`*^9,
3.9139692838751802`*^9}, {3.913969355649913*^9, 3.913969361712606*^9}, {
3.9144949578987226`*^9, 3.914495165499811*^9}, {3.9144959814900007`*^9,
3.914496003417218*^9}, {3.914496102391836*^9, 3.914496129028253*^9}, {
3.914496385924666*^9, 3.914496389370826*^9}, 3.914496816014628*^9,
3.914497137215846*^9, {3.914497195796404*^9, 3.914497205443459*^9}, {
3.9145002758142567`*^9, 3.914500278613076*^9}, {3.9145003217368507`*^9,
3.914500333937392*^9}, {3.9145003800922203`*^9, 3.9145004723683977`*^9}, {
3.914568352151114*^9, 3.91456837502531*^9}, {3.9145686816083937`*^9,
3.914568695760436*^9}, {3.914568755046906*^9, 3.91456877979802*^9}, {
3.914568842681856*^9, 3.914568856162018*^9}, {3.914570792465082*^9,
3.914570934687232*^9}, 3.914572872952714*^9, {3.914575174014627*^9,
3.914575204325823*^9}, {3.9146576852779713`*^9, 3.914657694073289*^9}, {
3.914664584430881*^9, 3.914664585007571*^9}, {3.914664615095174*^9,
3.914664616087166*^9}, {3.914669504535619*^9, 3.914669504946134*^9}, {
3.914677936628172*^9, 3.914677937450904*^9}, {3.9149227495921307`*^9,
3.914922788249671*^9}, {3.9149251411560698`*^9, 3.91492520194254*^9}, {
3.9149252435433693`*^9, 3.914925271566181*^9}, {3.91492532661316*^9,
3.914925344732032*^9}, 3.914925472345007*^9, {3.914927041152981*^9,
3.914927086751133*^9}, {3.914927703431086*^9, 3.9149277122625437`*^9}, {
3.915001177706603*^9, 3.915001195497541*^9}, {3.915076816729044*^9,
3.915076863122754*^9}, {3.915079097344879*^9, 3.915079097839446*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"7e42b61f-2302-4df4-bfab-65ccf677b957"],
Cell[BoxData["8"], "Output",
CellChangeTimes->{3.915000815668765*^9, 3.915001208341139*^9,
3.9150065026266193`*^9, 3.915076866743194*^9, 3.915078934525136*^9,
3.915079101029883*^9, 3.915088554662814*^9},
CellLabel->"Out[5]=",ExpressionUUID->"3a4b31e0-b243-4ba1-a6d8-35ec64ebb492"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
TagBox["\<\"calculate frame once:\"\>",
"EchoLabel"], " ", "\<\"with Simplify\"\>"}]], "Echo",
CellChangeTimes->{
3.9150885546684637`*^9},ExpressionUUID->"a82544bc-c8dd-49c0-8cf3-\
b499a141301d"],
Cell[BoxData[
RowBox[{
TagBox["\<\"calculate normals once:\"\>",
"EchoLabel"], " ", "\<\"with Simplify\"\>"}]], "Echo",
CellChangeTimes->{
3.915088554709632*^9},ExpressionUUID->"4e5ef113-da38-40b6-bcaa-\
fb70a555151e"],
Cell[BoxData[
RowBox[{
TagBox["\<\"calculate surface identity once:\"\>",
"EchoLabel"], " ", "\<\"with Simplify\"\>"}]], "Echo",
CellChangeTimes->{
3.915088554762478*^9},ExpressionUUID->"837e0956-6de1-4fa0-a2d3-\
f59da11107f5"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Q", "=",
RowBox[{"DeployQ2FromComponentsE", "[",
RowBox[{
RowBox[{"Q1", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q2", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q3", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q4", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q5", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"IxiQ", "=",
RowBox[{
RowBox[{"IdentityMatrix", "[", "3", "]"}], "-",
RowBox[{"\[Xi]", "*", "Q"}]}]}], ";"}]}], "Input",
CellChangeTimes->{{3.9150779625357018`*^9, 3.915078086586111*^9}},
CellLabel->"In[7]:=",ExpressionUUID->"839c8a5b-e1bd-4735-b613-ae11b6742429"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Q1", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q2", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q3", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Q2", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q4", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q5", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Q3", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{"Q5", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}], ",",
RowBox[{
RowBox[{"-",
RowBox[{"Q1", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}], "-",
RowBox[{"Q4", "[",
RowBox[{"t", ",", "y1", ",", "y2"}], "]"}]}]}], "}"}]}],
"}"}]], "Output",
CellChangeTimes->{3.915078030104127*^9, 3.915078119677905*^9,
3.915079101172957*^9, 3.915088554949638*^9},
CellLabel->"Out[7]=",ExpressionUUID->"4d7486e3-7396-474b-ab48-4fd66eed0897"]
}, Open ]],
Cell[BoxData[{
RowBox[{
RowBox[{"LLIxiQ", "=",
RowBox[{
RowBox[{"DtConvE", "[",
RowBox[{"IxiQ", ",", "\"\<L\>\""}], "]"}], "//", "SimplifyProxy"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Vm", "=",
RowBox[{
RowBox[{"GetMaterialVelocity", "[", "]"}], "//", "SimplifyProxy"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"DVm", "=",
RowBox[{
RowBox[{"CompD", "[", "Vm", "]"}], "//", "SimplifyProxy"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Gm", "=",
RowBox[{
RowBox[{"ProjTan", "[", "DVm", "]"}], "//", "SimplifyProxy"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Sm", "=",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"Gm", "+",
RowBox[{"Transpose", "[", "Gm", "]"}]}], ")"}], "/", "2"}], "//",
"SimplifyProxy"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.915078166414899*^9, 3.915078273780748*^9}, {
3.915078316954118*^9, 3.915078329594247*^9}, {3.9150785619670267`*^9,
3.915078725754545*^9}, {3.915088695329957*^9, 3.915088713960581*^9}},
CellLabel->"In[20]:=",ExpressionUUID->"a5d5b85a-e3cc-4254-b462-2b159c879137"],
Cell[BoxData[
RowBox[{
RowBox[{"IdS", "=",
RowBox[{
RowBox[{"GetSurfaceIdentity", "[", "]"}], "//", "SimplifyProxy"}]}],
";"}]], "Input",
CellChangeTimes->{{3.915079604573807*^9, 3.915079621336783*^9}},
CellLabel->"In[13]:=",ExpressionUUID->"e83123b4-fba2-4d9e-827f-2af71ebb391e"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"DtJauE", "[",
RowBox[{"IdentityMatrix", "[", "3", "]"}], "]"}], "//",
"SimplifyProxy"}]], "Input",
CellLabel->"In[14]:=",ExpressionUUID->"c171ba6f-bf4e-43fe-82c5-0f472c5819c5"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{3.915079272500348*^9, 3.91508857333036*^9},
CellLabel->"Out[14]=",ExpressionUUID->"da9baa8a-25f8-49b6-97c8-e742d06f5ef2"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"LLIxiQ", "-",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "\[Xi]"}], "*",
RowBox[{"DtJauE", "[", "Q", "]"}]}], "+",
RowBox[{"Sm", ".", "IxiQ"}], "+",
RowBox[{"IxiQ", ".", "Sm"}]}], ")"}]}], "//", "SimplifyProxy"}]], "Input",\
CellChangeTimes->{{3.9150787434077272`*^9, 3.915078757908471*^9}, {
3.9150789413203783`*^9, 3.915078980121265*^9}, {3.915079346051629*^9,
3.915079346582505*^9}},
CellLabel->"In[15]:=",ExpressionUUID->"c556fadc-cb7d-4895-a46c-3fa20132cdc0"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{3.9150789964247932`*^9, 3.915079105510364*^9,
3.915079349768502*^9, 3.915088574288764*^9},
CellLabel->"Out[15]=",ExpressionUUID->"8e524a83-29a9-42a2-95da-16f7e7de671d"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"Tr", "[",
RowBox[{"DtJauE", "[", "Q", "]"}], "]"}], "//", "Simplify"}]], "Input",
CellChangeTimes->{{3.9150794459930563`*^9, 3.915079457354939*^9}},
CellLabel->"In[16]:=",ExpressionUUID->"3fc498b2-7ae8-4e65-b7f3-98e1889ffdb7"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{3.915079458181699*^9, 3.915088574329247*^9},
CellLabel->"Out[16]=",ExpressionUUID->"8351b6ec-cf86-46a5-b1a3-e4a76846dc5c"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Rac", "=",
RowBox[{
RowBox[{
RowBox[{"1", "/", "2"}], "*",
RowBox[{"Tr", "[", "LLIxiQ", "]"}]}], "//", "Simplify"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Rac", "-",
RowBox[{"InnerPE", "[",
RowBox[{"IxiQ", ",", "Sm"}], "]"}]}], "//",
"Simplify"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Rac", "-",
RowBox[{"InnerPE", "[",
RowBox[{
RowBox[{"IdS", "-",
RowBox[{"\[Xi]", "*",
RowBox[{"IdS", ".", "Q"}]}]}], ",", "DVm"}], "]"}]}], "//",
"Simplify"}]}], "Input",
CellChangeTimes->{{3.915079473491105*^9, 3.9150795907732353`*^9}, {
3.915079627773424*^9, 3.915079651001985*^9}},
CellLabel->"In[17]:=",ExpressionUUID->"2320fe78-84f8-4a33-969a-4e131430810f"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{3.9150795289962273`*^9, 3.915079652651115*^9},
CellLabel->"Out[45]=",ExpressionUUID->"4622518e-f213-4207-8a16-627fa959df65"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{3.9150795289962273`*^9, 3.915079652739798*^9},
CellLabel->"Out[46]=",ExpressionUUID->"ccf1374d-77bc-4da0-9ff2-946c4c15680b"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{
3.915000815668765*^9, 3.915001208341139*^9, 3.9150065026266193`*^9,
3.915076866743194*^9, 3.915078934525136*^9, 3.915079101029883*^9, {
3.915088553370406*^9, 3.915088574787006*^9}},
CellLabel->"Out[18]=",ExpressionUUID->"424adf13-57d6-4103-af20-dee4029fc2cb"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{
3.915000815668765*^9, 3.915001208341139*^9, 3.9150065026266193`*^9,
3.915076866743194*^9, 3.915078934525136*^9, 3.915079101029883*^9, {
3.915088553370406*^9, 3.9150885748727617`*^9}},
CellLabel->"Out[19]=",ExpressionUUID->"09235e00-39d8-4da0-a5e1-8f7d0674d2e5"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Ria", "=",
RowBox[{
RowBox[{"CompTDiv", "[", "Vm", "]"}], "//", "Simplify"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Rna", "=",
RowBox[{
RowBox[{
RowBox[{"-", "\[Xi]"}], "*",
RowBox[{"InnerPE", "[",
RowBox[{"Q", ",",
RowBox[{"ProjQ2E", "[", "Sm", "]"}]}], "]"}]}], "//", "Simplify"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Rac", "-",
RowBox[{"(",
RowBox[{"Ria", "+", "Rna"}], ")"}]}], "//", "Simplify"}]}], "Input",
CellChangeTimes->{{3.9150886789479103`*^9, 3.915088681769356*^9}, {
3.915088721088037*^9, 3.915088813669271*^9}, 3.91508890853478*^9},
CellLabel->"In[28]:=",ExpressionUUID->"55c1cbe7-420c-49a3-8660-8f85f9a6134d"],
Cell[BoxData["0"], "Output",
CellChangeTimes->{3.915088814801549*^9, 3.915088909328094*^9},
CellLabel->"Out[30]=",ExpressionUUID->"deca8be7-8e2b-41e0-a1df-a4fcebe5d34d"]
}, Open ]]
},
WindowSize->{1440, 774},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification:>0.9 Inherited,
FrontEndVersion->"13.2 for Linux x86 (64-bit) (January 30, 2023)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"328941ce-6db8-4a43-af4b-73af12e51e98"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 4147, 80, 209, "Input",ExpressionUUID->"7e42b61f-2302-4df4-bfab-65ccf677b957"],
Cell[4730, 104, 288, 4, 30, "Output",ExpressionUUID->"3a4b31e0-b243-4ba1-a6d8-35ec64ebb492"],
Cell[CellGroupData[{
Cell[5043, 112, 228, 6, 24, "Echo",ExpressionUUID->"a82544bc-c8dd-49c0-8cf3-b499a141301d"],
Cell[5274, 120, 228, 6, 24, "Echo",ExpressionUUID->"4e5ef113-da38-40b6-bcaa-fb70a555151e"],
Cell[5505, 128, 237, 6, 24, "Echo",ExpressionUUID->"837e0956-6de1-4fa0-a2d3-f59da11107f5"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[5791, 140, 817, 21, 46, "Input",ExpressionUUID->"839c8a5b-e1bd-4735-b613-ae11b6742429"],
Cell[6611, 163, 1216, 34, 30, "Output",ExpressionUUID->"4d7486e3-7396-474b-ab48-4fd66eed0897"]
}, Open ]],
Cell[7842, 200, 1154, 33, 103, "Input",ExpressionUUID->"a5d5b85a-e3cc-4254-b462-2b159c879137"],
Cell[8999, 235, 296, 7, 26, "Input",ExpressionUUID->"e83123b4-fba2-4d9e-827f-2af71ebb391e"],
Cell[CellGroupData[{
Cell[9320, 246, 219, 5, 28, "Input",ExpressionUUID->"c171ba6f-bf4e-43fe-82c5-0f472c5819c5"],
Cell[9542, 253, 402, 10, 30, "Output",ExpressionUUID->"da9baa8a-25f8-49b6-97c8-e742d06f5ef2"]
}, Open ]],
Cell[CellGroupData[{
Cell[9981, 268, 547, 14, 28, "Input",ExpressionUUID->"c556fadc-cb7d-4895-a46c-3fa20132cdc0"],
Cell[10531, 284, 452, 11, 30, "Output",ExpressionUUID->"8e524a83-29a9-42a2-95da-16f7e7de671d"]
}, Open ]],
Cell[CellGroupData[{
Cell[11020, 300, 267, 5, 26, "Input",ExpressionUUID->"3fc498b2-7ae8-4e65-b7f3-98e1889ffdb7"],
Cell[11290, 307, 171, 2, 30, "Output",ExpressionUUID->"8351b6ec-cf86-46a5-b1a3-e4a76846dc5c"]
}, Open ]],
Cell[CellGroupData[{
Cell[11498, 314, 777, 23, 67, "Input",ExpressionUUID->"2320fe78-84f8-4a33-969a-4e131430810f"],
Cell[12278, 339, 173, 2, 30, "Output",ExpressionUUID->"4622518e-f213-4207-8a16-627fa959df65"],
Cell[12454, 343, 173, 2, 30, "Output",ExpressionUUID->"ccf1374d-77bc-4da0-9ff2-946c4c15680b"],
Cell[12630, 347, 318, 5, 30, "Output",ExpressionUUID->"424adf13-57d6-4103-af20-dee4029fc2cb"],
Cell[12951, 354, 320, 5, 30, "Output",ExpressionUUID->"09235e00-39d8-4da0-a5e1-8f7d0674d2e5"]
}, Open ]],
Cell[CellGroupData[{
Cell[13308, 364, 757, 21, 67, "Input",ExpressionUUID->"55c1cbe7-420c-49a3-8660-8f85f9a6134d"],
Cell[14068, 387, 171, 2, 30, "Output",ExpressionUUID->"deca8be7-8e2b-41e0-a1df-a4fcebe5d34d"]
}, Open ]]
}
]
*)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment