14.2: Variable Rescaling
- Page ID
- 7851
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Variable rescaling of continuous field models comes with yet another bonus variable, i.e., space, which you can rescale to potentially eliminate more parameters from the model. In a 2-D or higher dimensional space, you can, theoretically, have two or more spatial variables to rescale independently. But the space is usually assumed to be isotropic (i.e., there is no difference among the directions of space) in most spatial models, so it may not be practically meaningful to use different rescaling factors for different spatial variables.
Anyway, here is an example. Let’s try to simplify the following spatial predator-prey model, formulated as a reaction-diffusion system in a two-dimensional space:
\[ \begin{align} \dfrac{\partial{r}}{\partial{t}} &=ar-brf+D_{r}\nabla^{2}r \nonumber \\[4pt] &= ar-brf+D_{r} \left(\dfrac{\partial^{2}{r}}{\partial{x^{2}}} +\dfrac{\partial^{2}{r}}{\partial{y^{2}}}\right) \label{(14.20)} \end{align} \]
\[ \begin{align} \dfrac{\partial{f}}{\partial{t}} &=-cf+drf +D_{f}\nabla^{2}f \nonumber \\[4pt] &=-cf+drf+D_{f} \left(\dfrac{\partial^{2}f}{\partial{x^{2}}} +\dfrac{\partial^{2}f}{\partial{y^{2}}}\right) \label{(14.21)} \end{align} \]
Here we use \(r\) for prey (rabbits) and \(f\) for predators (foxes), since \(x\) and \(y\) are already taken as the spatial coordinates. We can apply the following three rescaling rules to state variables \(r\) and \(f\), as well as time \(t\) and space \(x/y\):
\[ \begin{align*} r &\rightarrow \alpha{r'} \label{(14.22)} \\[4pt] f &\rightarrow \beta{f'} \label{(14.23)} \\[4pt] t &\rightarrow \gamma{t'} \label{(14.24)} \\[4pt] x, y &\rightarrow \delta x, \delta y \label{(14.25)} \end{align*} \]
With these replacements, the model equations (Equation \ref{(14.20)} and \ref{(14.21)}) can be rewritten as follows:
\[ \begin{align} \dfrac{\partial(\alpha{r'})}{\partial{\gamma{t'}}} &=a(\alpha{r'}) -b(\alpha{r'})(\beta{f'}) +D_{r} \left(\dfrac{\partial^{2}(\alpha{r'}) }{\partial{(\delta{x'})^{2}}} +\dfrac{\partial^{2}(\alpha{r'})}{\partial{(\delta{y'})^{2}}}\right) \label{(14.26)} \\[4pt] \dfrac{\partial(\beta{f'})}{\partial{(\gamma{t'})}} &=-c(\beta{f'}) +d(\alpha{r'})(\beta{f'}) +D_{f} \left(\dfrac{\partial^{2}(\beta{f'})}{\partial{(\delta{x'})^{2}}} +\dfrac{\partial^{2}(\alpha{r'})}{\partial{(\delta{y'})}}^{2}\right) \label{(14.27)} \end{align} \]
We can collect parameters and rescaling factors together, as follows:
\[ \begin{align} \dfrac{\partial{r'}}{\partial{t'}} &= a\gamma{r'} -b{\beta}{\gamma{r' f'}} +\dfrac{D_{r}{\gamma}}{\delta^{2}}(\dfrac{\partial^{2}{r'}}{\partial{x^{'2}}}+\dfrac{\partial^{2}{r'}}{\partial{y^{'2}}})\label{(14.28)} \\[4pt] \dfrac{\partial{f'}}{\partial{t'}} &=-c{\gamma}f'+ d{\alpha{\gamma{r' f'}}} +\dfrac{D_{f}\gamma}{\delta^{2}}(\dfrac{\partial^{2}f'}{\partial{x^{'2}}}+\dfrac{\partial^{2}f'}{\partial{y^{'2}}}) \end{align} \]
Then we can apply, e.g., the following rescaling choices
\[(\alpha, \beta, \gamma, \delta) = \left(\dfrac{a}{d},\dfrac{a}{b}, \dfrac{1}{a}, \sqrt{\dfrac{D_{r}}{a}}\right) \label{(14.30)} \]
to simplify the model equations into
\[\dfrac{\partial{r'}}{\partial{t'}} =r' -r' f' +\nabla^{2}r', \label{(14.31)} \]
\[\dfrac{\partial{f'}}{\partial{t'}} =-e{\gamma{f'}} +r' f' +D_{ratio}\nabla^{2}f' \label{(14.32)} \]
with \(e = c/a \) and \(D_{ratio} = D_{f}/D_{r}\). The original model had six parameters (\(a\), \(b\), \(c\), \(d\), \(D_r\), and \(D_f\)), but we were able to reduce them into just two parameters (\(e\) and \(D_{ratio}\)) with a little additional help from spatial rescaling factor \(δ\). This rescaling result also tells us some important information about what matters in this system: It is the ratio between the growth rate of the prey (\(a\)) and the decay rate of the predators (\(c\)) (i.e., \(e = c/a\)), as well as the ratio between their diffusion speeds (\(D_{ratio} = D_{f}/D_{r}\)), which essentially determines the dynamics of this system. Both of these new parameters make a lot of sense from an ecological viewpoint too.
Simplify the following Keller-Segel model by variable rescaling:
\[ \begin{align*} \dfrac{\partial{a}}{\partial{t}} &= \mu{\nabla}^{2}a -\chi{\nabla} \cdot(a\nabla{c}) \label{(14.33)} \\[4pt] \dfrac{\partial{c}}{\partial{t}} &= D\nabla^{2}c +fa -kc \label{(14.34)} \end{align*} \]