Difference between revisions of "Advection in log coordinate"

From TORI
Jump to: navigation, search
 
m (Text replacement - "\$([^\$]+)\$" to "\\(\1\\)")
 
Line 1: Line 1:
 
[[Advection in log coordinate]] refers to the specific consideration of [[advection]], using the change of independent variable in the [[advection equation]]
 
[[Advection in log coordinate]] refers to the specific consideration of [[advection]], using the change of independent variable in the [[advection equation]]
   
$(1) ~ ~ ~ \dot F(r,t)= - V'(r) F(r,t) - V(r) F'(r,t)$
+
\((1) ~ ~ ~ \dot F(r,t)= - V'(r) F(r,t) - V(r) F'(r,t)\)
   
where $r$ has sense of the independent coordinate, for example, size of particles;
+
where \(r\) has sense of the independent coordinate, for example, size of particles;
   
$t$ has sense of time;
+
\(t\) has sense of time;
   
$V(r)$ determines the velocity, speed, rate of advection at coordinate $r$ and is assumed to be a given function,
+
\(V(r)\) determines the velocity, speed, rate of advection at coordinate \(r\) and is assumed to be a given function,
   
$F(r,t)$ is unknown function, that is supposed to be found as solution of equation (1), assuming some additional condition(s); for example $F(r,0)=F_0(r)$, where $F_0$ is given function.
+
\(F(r,t)\) is unknown function, that is supposed to be found as solution of equation (1), assuming some additional condition(s); for example \(F(r,0)=F_0(r)\), where \(F_0\) is given function.
   
 
Here and through [[TORI]], by default, the prime after the name of function indicates its derivative with respect to its first argument, and dot above the name of the function indicates its derivative with respect to the last argument.
 
Here and through [[TORI]], by default, the prime after the name of function indicates its derivative with respect to its first argument, and dot above the name of the function indicates its derivative with respect to the last argument.
   
This article considers the special case of transformation of the independent variable, namely, change from variable $r$ to variable $x=\ln(r)$ in such a way that $r=\exp(x)$.
+
This article considers the special case of transformation of the independent variable, namely, change from variable \(r\) to variable \(x=\ln(r)\) in such a way that \(r=\exp(x)\).
   
 
In order to write equation of advection in new coordinate, let
 
In order to write equation of advection in new coordinate, let
   
$(2) ~ ~ ~ \rho(x,t)=F(\mathrm e^x,t)~ \mathrm e^x$
+
\((2) ~ ~ ~ \rho(x,t)=F(\mathrm e^x,t)~ \mathrm e^x\)
   
The goal of this article is to write out the equation for function $\rho$, that follows from equation (1). It is desirable, however to have the resulting equation, that also can be interpreted as advection in the new coordinate.
+
The goal of this article is to write out the equation for function \(\rho\), that follows from equation (1). It is desirable, however to have the resulting equation, that also can be interpreted as advection in the new coordinate.
   
 
==Deduction==
 
==Deduction==
From equation (2), the distribution $F$ can be expressed through the new function $\rho$ as follows:
+
From equation (2), the distribution \(F\) can be expressed through the new function \(\rho\) as follows:
   
$(3) ~ ~ ~ \displaystyle F(r,t)=\frac{1}{r} \rho\Big(\ln(r),t\Big)$
+
\((3) ~ ~ ~ \displaystyle F(r,t)=\frac{1}{r} \rho\Big(\ln(r),t\Big)\)
   
 
Then, its derivative with respect to its first argument, that appears in the right hand side of equation (1), can be written as follows:
 
Then, its derivative with respect to its first argument, that appears in the right hand side of equation (1), can be written as follows:
   
$(4) ~ ~ ~ \displaystyle F'(r,t)= \frac{-1}{r^2} \rho\Big(\ln(r),t\Big) + \frac{1}{r}\rho\Big(\ln(r),t\Big) \frac{1}{r}$
+
\((4) ~ ~ ~ \displaystyle F'(r,t)= \frac{-1}{r^2} \rho\Big(\ln(r),t\Big) + \frac{1}{r}\rho\Big(\ln(r),t\Big) \frac{1}{r}\)
   
Change of variable $r \mapsto \exp(x)$ gives
+
Change of variable \(r \mapsto \exp(x)\) gives
   
$(5) ~ ~ ~ \displaystyle F'(\mathrm e^x,t)= - \mathrm e ^{-2x} \rho(x,t)+\mathrm e^{-2x} \rho'(x,t)$
+
\((5) ~ ~ ~ \displaystyle F'(\mathrm e^x,t)= - \mathrm e ^{-2x} \rho(x,t)+\mathrm e^{-2x} \rho'(x,t)\)
   
With this expression, from equation (1), the derivative of $\rho$ with respect to the last argument can be written as follows:
+
With this expression, from equation (1), the derivative of \(\rho\) with respect to the last argument can be written as follows:
   
$(6) ~ ~ ~ \displaystyle \dot \rho(x,t)= \dot F(\mathrm e^x,t) ~ \mathrm e^x$
+
\((6) ~ ~ ~ \displaystyle \dot \rho(x,t)= \dot F(\mathrm e^x,t) ~ \mathrm e^x\)
$= \Big(
+
\(= \Big(
 
- V'(\mathrm e^x)\, F(\mathrm e^x,t) -V(\mathrm e^x)\, F'(\mathrm e^x,t)
 
- V'(\mathrm e^x)\, F(\mathrm e^x,t) -V(\mathrm e^x)\, F'(\mathrm e^x,t)
\Big)~ \mathrm e^x$
+
\Big)~ \mathrm e^x\)
$=-V'(\mathrm e^x)\, F(\mathrm e^x,t)\,\mathrm e^x - V(\mathrm e^x)\,\mathrm e^x\, F'(\mathrm e^x,t)$
+
\(=-V'(\mathrm e^x)\, F(\mathrm e^x,t)\,\mathrm e^x - V(\mathrm e^x)\,\mathrm e^x\, F'(\mathrm e^x,t)\)
$= -V'(\mathrm e^x)\,\rho(x,t) - V(\mathrm e^x)\,\Big(- \mathrm e ^{-x} \rho(x,t)+\mathrm e^{-x} \rho'(x,t)\Big)$
+
\(= -V'(\mathrm e^x)\,\rho(x,t) - V(\mathrm e^x)\,\Big(- \mathrm e ^{-x} \rho(x,t)+\mathrm e^{-x} \rho'(x,t)\Big)\)
$= -\Big(V'(\mathrm e^x) - \mathrm e^{-x} V(\mathrm e^x)\Big)\, \rho(x)
+
\(= -\Big(V'(\mathrm e^x) - \mathrm e^{-x} V(\mathrm e^x)\Big)\, \rho(x)
-V(\mathrm e^x)\, \mathrm e^{-x} \rho'(x,t)$
+
-V(\mathrm e^x)\, \mathrm e^{-x} \rho'(x,t)\)
   
In order to interpret equation (6) as the [[advection equation]], we need to define a new function $W$, that have sense of velocity of drift; let
+
In order to interpret equation (6) as the [[advection equation]], we need to define a new function \(W\), that have sense of velocity of drift; let
   
$(7) ~ ~ ~ \displaystyle W(x)=V(\mathrm e^x) \, \mathrm e^{-x}$
+
\((7) ~ ~ ~ \displaystyle W(x)=V(\mathrm e^x) \, \mathrm e^{-x}\)
   
 
Then, its derivative
 
Then, its derivative
   
$(8) ~ ~ ~ \displaystyle W('x)=V'(\mathrm e^x) -V(\mathrm e^x)\, \mathrm e^{-x}$
+
\((8) ~ ~ ~ \displaystyle W('x)=V'(\mathrm e^x) -V(\mathrm e^x)\, \mathrm e^{-x}\)
   
and this is coefficient at the first term in the right hand side of equation (6). With function $W$, equation (6) can be rewritten as follows:
+
and this is coefficient at the first term in the right hand side of equation (6). With function \(W\), equation (6) can be rewritten as follows:
   
$(9) ~ ~ ~ \displaystyle
+
\((9) ~ ~ ~ \displaystyle
\dot \rho(x,t)= W'(x)\, \rho(x,t) - W(x)\, \rho'(x,t) $
+
\dot \rho(x,t)= W'(x)\, \rho(x,t) - W(x)\, \rho'(x,t) \)
   
This is also advection equation (without diffusion); and this article is loaded to [[TORI]] in order to simplify the search of the deduction of this equation. Equation (9) indicates, that $W(x)$ can be interpreted as the drift rate at coordinate $x$, while $\rho(x,t)$ has sense of density of distribution of some quantity along this coordinate.
+
This is also advection equation (without diffusion); and this article is loaded to [[TORI]] in order to simplify the search of the deduction of this equation. Equation (9) indicates, that \(W(x)\) can be interpreted as the drift rate at coordinate \(x\), while \(\rho(x,t)\) has sense of density of distribution of some quantity along this coordinate.
   
 
==Lognormal distribution==
 
==Lognormal distribution==
 
[[File:Lognot300.jpg|500px|thumb|Explicit plot of lognormal distribution by (14)]]
 
[[File:Lognot300.jpg|500px|thumb|Explicit plot of lognormal distribution by (14)]]
   
Change of variable $r \mapsto x=\ln(t)$ has sense in order to apply the Fourier methods for the advection with distribution similar to the lognormal distribution. This refers to the case, when the initial condition for $\rho$ id Gaussian, for example,
+
Change of variable \(r \mapsto x=\ln(t)\) has sense in order to apply the Fourier methods for the advection with distribution similar to the lognormal distribution. This refers to the case, when the initial condition for \(\rho\) id Gaussian, for example,
   
$(10) ~ ~ ~ \displaystyle \rho(x,0)= \frac{1}{\sqrt{\pi}} \exp(-x^2)$
+
\((10) ~ ~ ~ \displaystyle \rho(x,0)= \frac{1}{\sqrt{\pi}} \exp(-x^2)\)
   
This corresponds to the lognormal initial condition for the function $F$. This can be seen from the deduction below.
+
This corresponds to the lognormal initial condition for the function \(F\). This can be seen from the deduction below.
This deduction refers to $t\!=\!0$, so, the last argument in the deduction is omitted.
+
This deduction refers to \(t\!=\!0\), so, the last argument in the deduction is omitted.
 
Let
 
Let
   
$(11) ~ ~ ~ \rho(x)\, \mathrm dx = F(r)\, \mathrm d r= F(\mathrm e^x)\, \mathrm e^x\, \mathrm d x $
+
\((11) ~ ~ ~ \rho(x)\, \mathrm dx = F(r)\, \mathrm d r= F(\mathrm e^x)\, \mathrm e^x\, \mathrm d x \)
   
 
Then
 
Then
   
$(12) ~ ~ ~ \rho(x)=F(\mathrm e^x) \mathrm e^x$
+
\((12) ~ ~ ~ \rho(x)=F(\mathrm e^x) \mathrm e^x\)
   
$(13) ~ ~ ~ \rho\big(\ln(r)\big)=F(r) r$
+
\((13) ~ ~ ~ \rho\big(\ln(r)\big)=F(r) r\)
   
$(14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r} \rho\big(\ln(r)\Big)$
+
\((14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r} \rho\big(\ln(r)\Big)\)
   
 
For the Gaussian (10), one this gives the initial lognormal function
 
For the Gaussian (10), one this gives the initial lognormal function
   
$(14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r \sqrt{\pi}} \exp\Big(-\ln(r)^2\Big)$
+
\((14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r \sqrt{\pi}} \exp\Big(-\ln(r)^2\Big)\)
   
 
This function is shown in the Figure at right.
 
This function is shown in the Figure at right.

Latest revision as of 18:44, 30 July 2019

Advection in log coordinate refers to the specific consideration of advection, using the change of independent variable in the advection equation

\((1) ~ ~ ~ \dot F(r,t)= - V'(r) F(r,t) - V(r) F'(r,t)\)

where \(r\) has sense of the independent coordinate, for example, size of particles;

\(t\) has sense of time;

\(V(r)\) determines the velocity, speed, rate of advection at coordinate \(r\) and is assumed to be a given function,

\(F(r,t)\) is unknown function, that is supposed to be found as solution of equation (1), assuming some additional condition(s); for example \(F(r,0)=F_0(r)\), where \(F_0\) is given function.

Here and through TORI, by default, the prime after the name of function indicates its derivative with respect to its first argument, and dot above the name of the function indicates its derivative with respect to the last argument.

This article considers the special case of transformation of the independent variable, namely, change from variable \(r\) to variable \(x=\ln(r)\) in such a way that \(r=\exp(x)\).

In order to write equation of advection in new coordinate, let

\((2) ~ ~ ~ \rho(x,t)=F(\mathrm e^x,t)~ \mathrm e^x\)

The goal of this article is to write out the equation for function \(\rho\), that follows from equation (1). It is desirable, however to have the resulting equation, that also can be interpreted as advection in the new coordinate.

Deduction

From equation (2), the distribution \(F\) can be expressed through the new function \(\rho\) as follows:

\((3) ~ ~ ~ \displaystyle F(r,t)=\frac{1}{r} \rho\Big(\ln(r),t\Big)\)

Then, its derivative with respect to its first argument, that appears in the right hand side of equation (1), can be written as follows:

\((4) ~ ~ ~ \displaystyle F'(r,t)= \frac{-1}{r^2} \rho\Big(\ln(r),t\Big) + \frac{1}{r}\rho\Big(\ln(r),t\Big) \frac{1}{r}\)

Change of variable \(r \mapsto \exp(x)\) gives

\((5) ~ ~ ~ \displaystyle F'(\mathrm e^x,t)= - \mathrm e ^{-2x} \rho(x,t)+\mathrm e^{-2x} \rho'(x,t)\)

With this expression, from equation (1), the derivative of \(\rho\) with respect to the last argument can be written as follows:

\((6) ~ ~ ~ \displaystyle \dot \rho(x,t)= \dot F(\mathrm e^x,t) ~ \mathrm e^x\) \(= \Big( - V'(\mathrm e^x)\, F(\mathrm e^x,t) -V(\mathrm e^x)\, F'(\mathrm e^x,t) \Big)~ \mathrm e^x\) \(=-V'(\mathrm e^x)\, F(\mathrm e^x,t)\,\mathrm e^x - V(\mathrm e^x)\,\mathrm e^x\, F'(\mathrm e^x,t)\) \(= -V'(\mathrm e^x)\,\rho(x,t) - V(\mathrm e^x)\,\Big(- \mathrm e ^{-x} \rho(x,t)+\mathrm e^{-x} \rho'(x,t)\Big)\) \(= -\Big(V'(\mathrm e^x) - \mathrm e^{-x} V(\mathrm e^x)\Big)\, \rho(x) -V(\mathrm e^x)\, \mathrm e^{-x} \rho'(x,t)\)

In order to interpret equation (6) as the advection equation, we need to define a new function \(W\), that have sense of velocity of drift; let

\((7) ~ ~ ~ \displaystyle W(x)=V(\mathrm e^x) \, \mathrm e^{-x}\)

Then, its derivative

\((8) ~ ~ ~ \displaystyle W('x)=V'(\mathrm e^x) -V(\mathrm e^x)\, \mathrm e^{-x}\)

and this is coefficient at the first term in the right hand side of equation (6). With function \(W\), equation (6) can be rewritten as follows:

\((9) ~ ~ ~ \displaystyle \dot \rho(x,t)= W'(x)\, \rho(x,t) - W(x)\, \rho'(x,t) \)

This is also advection equation (without diffusion); and this article is loaded to TORI in order to simplify the search of the deduction of this equation. Equation (9) indicates, that \(W(x)\) can be interpreted as the drift rate at coordinate \(x\), while \(\rho(x,t)\) has sense of density of distribution of some quantity along this coordinate.

Lognormal distribution

Explicit plot of lognormal distribution by (14)

Change of variable \(r \mapsto x=\ln(t)\) has sense in order to apply the Fourier methods for the advection with distribution similar to the lognormal distribution. This refers to the case, when the initial condition for \(\rho\) id Gaussian, for example,

\((10) ~ ~ ~ \displaystyle \rho(x,0)= \frac{1}{\sqrt{\pi}} \exp(-x^2)\)

This corresponds to the lognormal initial condition for the function \(F\). This can be seen from the deduction below. This deduction refers to \(t\!=\!0\), so, the last argument in the deduction is omitted. Let

\((11) ~ ~ ~ \rho(x)\, \mathrm dx = F(r)\, \mathrm d r= F(\mathrm e^x)\, \mathrm e^x\, \mathrm d x \)

Then

\((12) ~ ~ ~ \rho(x)=F(\mathrm e^x) \mathrm e^x\)

\((13) ~ ~ ~ \rho\big(\ln(r)\big)=F(r) r\)

\((14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r} \rho\big(\ln(r)\Big)\)

For the Gaussian (10), one this gives the initial lognormal function

\((14) ~ ~ ~ \displaystyle F(r)=\frac{1}{r \sqrt{\pi}} \exp\Big(-\ln(r)^2\Big)\)

This function is shown in the Figure at right. Similar (a little bit more general) form of this lognormal distribution appears at Wikipedia, http://en.wikipedia.org/wiki/Log-normal_distribution

Discussion

Conclusion

Keywords

Advection, Derivative, Equation in partial derivatives, Exp, Log, Lognormal distribution

References

http://en.wikipedia.org/wiki/Log-normal_distribution

Deduction in this article is copied from Dima2014, page 115.