Skip to content

Commit

Permalink
WIP backup for Sonoma upgrade
Browse files Browse the repository at this point in the history
  • Loading branch information
msleigh committed Jan 3, 2024
1 parent d408bf8 commit 74fe0fb
Show file tree
Hide file tree
Showing 3 changed files with 6 additions and 13 deletions.
1 change: 1 addition & 0 deletions .gitignore
Original file line number Diff line number Diff line change
Expand Up @@ -9,3 +9,4 @@
docs/html
docs/latex
docs/doxygen.log
.venv
2 changes: 1 addition & 1 deletion docs/Doxyfile
Original file line number Diff line number Diff line change
Expand Up @@ -1859,7 +1859,7 @@ LATEX_OUTPUT = latex
# the output language.
# This tag requires that the tag GENERATE_LATEX is set to YES.

LATEX_CMD_NAME = latex
LATEX_CMD_NAME = pdflatex

# The MAKEINDEX_CMD_NAME tag can be used to specify the command name to generate
# index for LaTeX.
Expand Down
16 changes: 4 additions & 12 deletions docs/theory.md
Original file line number Diff line number Diff line change
Expand Up @@ -147,35 +147,27 @@ $$

is rearranged to give:

$$
u_r^{n+1} \approx
\left[ \frac{1-\left(1-\alpha\right)\hat\beta c\Delta t\hat\sigma}{1+\alpha\hat\beta c\Delta t\hat\sigma} \right] u^n_r
+
\left[ \frac{\hat\beta\hat\sigma\Delta t}{1+\alpha\hat\beta c\Delta t\hat\sigma} \right] \phi^\lambda
$$
$$ u_r^{n+1} \approx \left[ \frac{1-\left(1-\alpha\right)\hat\beta c\Delta t\hat\sigma}{1+\alpha\hat\beta c\Delta t\hat\sigma} \right] u^n_r + \left[ \frac{\hat\beta\hat\sigma\Delta t}{1+\alpha\hat\beta c\Delta t\hat\sigma} \right] \phi^\lambda $$

which allows us to write:

$$ \alpha u^{n+1}_r + \left(1-\alpha\right) u^n_r =
\frac{\alpha\hat\beta\hat\sigma\Delta t}{1 + \alpha\hat\beta c\Delta t\hat\sigma} \phi^\lambda
+
\frac{\alpha\hat\beta\hat\sigma\Delta t}{1 + \alpha\hat\beta c\Delta t\hat\sigma} \phi^\lambda +
\frac{1}{1 + \alpha\hat\beta c\Delta t\hat\sigma} u^n_r
$$

Substituting this into the transport equation gives:

$$ \frac{1}{c} \frac{\partial I}{\partial t} + \mu \frac{\partial I}{\partial x} + \sigma I
= \frac{1}{2}\sigma
\left[ \frac{\alpha\hat\beta c\hat\sigma\Delta t}{1 + \alpha\hat\beta c\Delta t\hat\sigma} \phi^\lambda
+
\left[ \frac{\alpha\hat\beta c\hat\sigma\Delta t}{1 + \alpha\hat\beta c\Delta t\hat\sigma} \phi^\lambda +
\frac{c}{1 + \alpha\hat\beta c\Delta t\hat\sigma} u^n_r\right]
$$

which we can re-write as:

$$ \frac{1}{c} \frac{\partial I}{\partial t} + \mu \frac{\partial I}{\partial x} + \sigma I =
\frac{1}{2} \sigma_s \phi^\lambda
+
\frac{1}{2} \sigma_s \phi^\lambda +
\frac{1}{2} c\sigma_a u^n_r
$$

Expand Down

0 comments on commit 74fe0fb

Please sign in to comment.