Skip to content

Commit

Permalink
changed pde->PDE and ode->ODE
Browse files Browse the repository at this point in the history
  • Loading branch information
victorballester7 committed Oct 26, 2023
1 parent cdeaa95 commit 528772a
Show file tree
Hide file tree
Showing 15 changed files with 271 additions and 116 deletions.
122 changes: 61 additions & 61 deletions Mathematics/3rd/Differential_equations/Differential_equations.tex

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
Expand Up @@ -885,7 +885,7 @@
Let $S\subseteq\RR^3$ be a surface and $\vf{X}$, $\vf{Y}$ be vector fields tangent to $S$ along a curve $\vf\alpha:I\rightarrow S$ of class $\mathcal{C}^\infty$ such that they are parallel. Then, $t\mapsto\langle \vf{X}(t),\vf{Y}(t)\rangle$ is constant. In particular, the norms $\|\vf{X}(t)\|$, $\|\vf{Y}(t)\|$ as well as the angle between $\vf{X}(t)$ and $\vf{Y}(t)$ are constant.
\end{proposition}
\begin{proposition}
Let $S\subseteq\RR^3$ be a surface, $(V,\vf\varphi(u,v))$ is a parametrization of $S$ and $\vf\alpha: I\rightarrow S$ be a parametrized curve of class $\mathcal{C}^\infty$ such that $\vf\alpha=\vf\varphi(u(t),v(t))$. Then, given $t_0\in I$ and $\vf{w}\in T_{\vf\alpha(t_0)}S$ there exists a unique parallel vector field $\vf{X}=a\vf\varphi_u+b\vf\varphi_v$ along $\vf\alpha$ such that $\vf{X}(t_0)=\vf{w}$. This vector field is called \emph{parallel transport} of the vector $\vf{w}$ along $\vf{\alpha}$, and it is defined on the entire interval $I$. It can be found by solving this system of odes:
Let $S\subseteq\RR^3$ be a surface, $(V,\vf\varphi(u,v))$ is a parametrization of $S$ and $\vf\alpha: I\rightarrow S$ be a parametrized curve of class $\mathcal{C}^\infty$ such that $\vf\alpha=\vf\varphi(u(t),v(t))$. Then, given $t_0\in I$ and $\vf{w}\in T_{\vf\alpha(t_0)}S$ there exists a unique parallel vector field $\vf{X}=a\vf\varphi_u+b\vf\varphi_v$ along $\vf\alpha$ such that $\vf{X}(t_0)=\vf{w}$. This vector field is called \emph{parallel transport} of the vector $\vf{w}$ along $\vf{\alpha}$, and it is defined on the entire interval $I$. It can be found by solving this system ofODEs:
$$\left\{
\begin{aligned}
a'+\Gamma_{11}^1au'+\Gamma_{12}^1av'+\Gamma_{21}^1bu'+\Gamma_{22}^1bv' & =0 \\
Expand Down Expand Up @@ -991,7 +991,7 @@
\end{lemma}
\begin{definition}
Let $U\subseteq\RR^n$ be an open set and $\vf{X}=\sum X^i\pdv{}{x^i}\in\mathcal{X}(U)$. We say that a parametrized curve $\vf{\gamma}:I\rightarrow\RR^n$ is an \emph{integral curve} of $\vf{X}$ if: $$\vf\gamma'(t)=\vf{X}(\vf\gamma(t))\qquad \forall t\in I$$
That is, the integral curve $\vf\gamma(t)=(x^1(t),\ldots,x^n(t))$ of $\vf{X}$ satisfies the following system of odes:
That is, the integral curve $\vf\gamma(t)=(x^1(t),\ldots,x^n(t))$ of $\vf{X}$ satisfies the following system ofODEs:
$$
\left\{
\begin{aligned}
Expand Down
6 changes: 3 additions & 3 deletions Mathematics/4th/Dynamical_systems/Dynamical_systems.tex
Original file line number Diff line number Diff line change
Expand Up @@ -36,7 +36,7 @@
\end{enumerate}
\end{theorem}
\begin{definition}
Let $f,g:\RR^2\rightarrow\RR$ be two functions and consider the system of odes:
Let $f,g:\RR^2\rightarrow\RR$ be two functions and consider the system ofODEs:
\begin{equation}\label{DS:plane}
\left\{
\begin{aligned}
Expand Down Expand Up @@ -645,14 +645,14 @@
A semiestable limit cycle $\Gamma_\mu$ of a family of rotated vector fields splits into two simple limit cycles, one stable and one unstable, as the parameter $\mu$ is varied in one sense and it disappears as $\mu$ is varied in the opposite sense.
\end{theorem}
\begin{theorem}[Melnikov's method]
Let $\vf{f}\in\mathcal{C}^1(\RR^2)$, $\vf{g}\in\mathcal{C}^1(\RR^2\times\RR^m)$ and $\varepsilon\simeq 0$. Consider the following ode:
Let $\vf{f}\in\mathcal{C}^1(\RR^2)$, $\vf{g}\in\mathcal{C}^1(\RR^2\times\RR^m)$ and $\varepsilon\simeq 0$. Consider the followingODE:
\begin{equation}\label{DS:melnikov}
\vf{x}'=\vf{f}(\vf{x})+\varepsilon\vf{g}(\vf{x},\vf{\mu})
\end{equation}
Suppose that for $\varepsilon =0$ the system has a one-parameter family of periodic orbits $\vf\gamma_h(t)$ of period $T_h$. Then for any simple zero $(\vf\mu_0,h_0)$ of the function $$M(\vf\mu, h)=\int_{0}^{T_h}\vf{f}(\vf\gamma_h(t))\times \vf{g}(\vf\gamma_h(t))\dd{t}$$ there exists a unique limit cycle $\vf\Gamma_\varepsilon$ for $\varepsilon\simeq 0$ such that $\displaystyle\lim_{\varepsilon\to 0}\vf\Gamma_\varepsilon=\vf\gamma_{h_0}$. On the other hand, if $M(\vf\mu_0,h_0)\ne 0$, for sufficiently small $\varepsilon$, the system of \mcref{DS:melnikov} with $\vf\mu=\vf\mu_0$ has no limit cycle in any sufficiently small neighborhood of $\vf\gamma_{h_0}$.
\end{theorem}
\begin{corollary}[Melnikov's method]
Let $H\in\mathcal{C}^2(\RR^2)$, $P,Q\in\mathcal{C}^1(\RR^2\times\RR^m)$ and $\varepsilon\simeq 0$. Consider the following system of odes:
Let $H\in\mathcal{C}^2(\RR^2)$, $P,Q\in\mathcal{C}^1(\RR^2\times\RR^m)$ and $\varepsilon\simeq 0$. Consider the following system ofODEs:
\begin{equation*}
\left\{
\begin{aligned}
Expand Down
6 changes: 3 additions & 3 deletions Mathematics/4th/Harmonic_analysis/Harmonic_analysis.tex
Original file line number Diff line number Diff line change
Expand Up @@ -166,7 +166,7 @@
Let $f(x)=\exp{-a x^2}$. Then, $\F f(\xi)=\sqrt{\frac{\pi}{a}}\exp{-\frac{{(\pi \xi)}^2}{a}}$ and moreover $\F^2f=f$. In particular if $a=\pi$, then $\F f=f$.
\end{lemma}
\begin{sproof}
$f$ satisfies the ode $y'=-2a x y$. Taking $\ \widehat{}\ $ on this expression and using \mcref{HA:diffFourierXf,HA:diffFourierTransf} we obtain that $\widehat{f}$ must satisfy the following ode:
$f$ satisfies theODE $y'=-2a x y$. Taking $\ \widehat{}\ $ on this expression and using \mcref{HA:diffFourierXf,HA:diffFourierTransf} we obtain that $\widehat{f}$ must satisfy the followingODE:
$$y'=-\frac{2\pi^2\xi}{a} y$$
with initial condition $y(0)=\int_{-\infty}^{+\infty}\exp{-a x^2}\dd{x}=\sqrt{\frac{\pi}{a}}$.
\end{sproof}
Expand Down Expand Up @@ -573,7 +573,7 @@
\end{theorem}
\subsubsection{Applications of the Fourier transform}
\begin{remark}
Probably the most important application of Fourier series is the resolution of pdes and it is a consequence of \mcref{HA:diffFourierTransf}, which reduces any order of a pde in the spatial variable to 1. The procedure is to compute the Fourier transform $\F$ of the pde, solve it, and then get back to the first function using the inverse transform.
Probably the most important application of Fourier series is the resolution of PDEs and it is a consequence of \mcref{HA:diffFourierTransf}, which reduces any order of a PDE in the spatial variable to 1. The procedure is to compute the Fourier transform $\F$ of the PDE, solve it, and then get back to the first function using the inverse transform.
\end{remark}
\begin{theorem}[Uncertainty principle]
Let $f\in L^2(\RR)$ be differentiable such that $x\abs{f}^2\in L^1(\RR)$ and $f'\in L^2(\RR)$. Then:
Expand Down Expand Up @@ -1266,7 +1266,7 @@
$$
\partial_t \widehat{E}+4\pi^2a^2\norm{\vf\xi}^2\widehat{E}=\delta_t
$$
because $\delta=\delta_{\vf{x}}\delta_t$. It can be seen that a solution of this ode is:
because $\delta=\delta_{\vf{x}}\delta_t$. It can be seen that a solution of thisODE is:
$$
\widehat{E}(t,\xi)=\indi{[0,\infty)}(t)\exp{-4\pi^2a^2\norm{\vf\xi}^2t}
$$
Expand Down
2 changes: 1 addition & 1 deletion Mathematics/4th/Numerical_calculus/Numerical_calculus.tex
Original file line number Diff line number Diff line change
Expand Up @@ -84,7 +84,7 @@
Moreover, we say that the algorithm has \emph{order of convergence} $p$ if $\norm{\vf{e}_n}=\O{h^p}$.
\end{definition}
\begin{remark}
Note that in a consistent method the difference equation for the method approaches the ode as the step size goes to zero, whereas in a convergent method is the solution to the difference equation that approaches the solution to the ode as the step size goes to zero.
Note that in a consistent method the difference equation for the method approaches theODE as the step size goes to zero, whereas in a convergent method is the solution to the difference equation that approaches the solution to theODE as the step size goes to zero.
\end{remark}
\begin{theorem}\label{NC:errorLipschitz}
Consider a consistent one-step explicit method such that its incremental function $\vf\phi$ is Lipschitz continuous (with constant $L$) with respect to $\vf{x}$. Then:
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
\subsection{Finite difference schemes}
\subsubsection{Introduction}
\begin{definition}
A linear system of $n$ first order of pdes for $\vf{u}(t,x)$ is a system of the form: $$\vf{A}(t,x)\vf{u}_t+\vf{B}(t,x)\vf{u}_x=\vf{C}(t,x)\vf{u}+\vf{D}(t,x)$$
A linear system of $n$ first order of PDEs for $\vf{u}(t,x)$ is a system of the form: $$\vf{A}(t,x)\vf{u}_t+\vf{B}(t,x)\vf{u}_x=\vf{C}(t,x)\vf{u}+\vf{D}(t,x)$$
for certain matrices $\vf{A},\vf{B},\vf{C}, \vf{D}\in \mathcal{M}_q(\RR)$. The system is called \emph{hyperbolic} if $\vf{A}^{-1}\vf{B}$ is diagonalizable.
\end{definition}
\begin{definition}
Expand Down Expand Up @@ -70,23 +70,23 @@
\end{definition}
\begin{definition}
Let $(G_j)$ be a sequence of grids such that the time and space steps $k_j,h_j>0$ of each one satisfy $\displaystyle \lim_{j\to\infty}k_j=\lim_{j\to\infty}h_j=0$.
We say that a finite difference scheme $v$ approximating a pde with initial condition $u_0(x)$ is \emph{unconditionally convergent} if for any solution $u(x,t)$ to the pde we have:
We say that a finite difference scheme $v$ approximating a PDE with initial condition $u_0(x)$ is \emph{unconditionally convergent} if for any solution $u(x,t)$ to the PDE we have:
\begin{itemize}
\item For all $x\in\domain u_0$ and all increasing sequence $(m_j)\in\NN$ such that $(\cdot,x_{m_j})\in G_j$ and $\displaystyle\lim_{j\to\infty} x_{m_j}=x$, we have $\displaystyle\lim_{j\to\infty} v_{m_j}^0=u_0(x)$.
\item For all $(t,x)\in\domain u$ and all increasing sequences $(m_j),(n_j)\in\NN$ such that $(t_{n_j},x_{m_j})\in G_j$ and $\displaystyle\lim_{j\to\infty} x_{m_j}=x$, $\displaystyle\lim_{j\to\infty} t_{n_j}=t$, we have $\displaystyle\lim_{j\to\infty} v_{m_j}^{n_j}=u(t,x)$.
\end{itemize}
The scheme is \emph{conditionally convergent} if $\forall j\in\NN$, $(k_j,h_j)\in\Lambda$, for some stability region $\Lambda$.
\end{definition}
\begin{definition}
Let $P$ be a partial differential operator and $\vf{f}$ be a function. Given the pde $P\vf{u}=\vf{f}$ and a finite difference scheme $P_{k,h}\vf{v}=R_{k,h}\vf{f}$ with $R_{k,h}\vf{1}=\vf{1}$, we say that the scheme is \emph{consistent} with the pde if for any smooth function $\vf\phi(t,x)$ we have: $$\lim_{k,h\to 0}R_{k,h}P\vf\phi-P_{k,h}\vf\phi=\vf{0}$$
Let $P$ be a partial differential operator and $\vf{f}$ be a function. Given the PDE $P\vf{u}=\vf{f}$ and a finite difference scheme $P_{k,h}\vf{v}=R_{k,h}\vf{f}$ with $R_{k,h}\vf{1}=\vf{1}$, we say that the scheme is \emph{consistent} with the PDE if for any smooth function $\vf\phi(t,x)$ we have: $$\lim_{k,h\to 0}R_{k,h}P\vf\phi-P_{k,h}\vf\phi=\vf{0}$$
where the convergence is pointwise at each point $(t,x)$ in the domain of solutions. We say that the consistency is of order $(p,q)$ in time and space if: $$\lim_{k,h\to 0}R_{k,h}P\vf\phi-P_{k,h}\vf\phi=\O{k^p}+\O{h^q}$$ The consistency is a \emph{conditional consistency} if the limit is for $(k,h)\in \Lambda$, for some stability region $\Lambda$. In that case, it makes sense to say that the consistency is of order $r$ in $k=\lambda(h)$ if:
$$\lim_{h\to 0}R_{\lambda(h),h}P\vf\phi-P_{\lambda(h),h}\vf\phi=\O{h^r}$$
\end{definition}
\begin{lemma}
The Lax-Friedrichs scheme is consistent if and only if $\displaystyle\lim_{h,k\to 0}\frac{h^2}{k}=0$.
\end{lemma}
\begin{remark}
The consistency is not enough to guarantee convergence. For example, consider the pde $u_t+au_x=0$, with $a>0$. The forward-time forward-space scheme is consistent with the pde, but it is not convergent if we take the initial condition $u_0(x)=\indi{\{x<0\}}$ on the domain $[-1,1]$. Indeed, looking at \mcref{NIPDE:upwind} we see that from some instant of time, the solution will be $0$ everywhere, which cannot be possible. In that case we should use the forward-time backward-space scheme, which is convergent. The usage of this latter method in these cases is called the \emph{upwind condition}.
The consistency is not enough to guarantee convergence. For example, consider the PDE $u_t+au_x=0$, with $a>0$. The forward-time forward-space scheme is consistent with the PDE, but it is not convergent if we take the initial condition $u_0(x)=\indi{\{x<0\}}$ on the domain $[-1,1]$. Indeed, looking at \mcref{NIPDE:upwind} we see that from some instant of time, the solution will be $0$ everywhere, which cannot be possible. In that case we should use the forward-time backward-space scheme, which is convergent. The usage of this latter method in these cases is called the \emph{upwind condition}.
\end{remark}
\begin{figure}[H]
\centering
Expand Down Expand Up @@ -204,7 +204,7 @@
\end{align*}
\end{proof}
\begin{proposition}
Consider the pde of \mcref{NIPDE:traffic} with $\lambda =k/h=\const$ Then:
Consider the PDE of \mcref{NIPDE:traffic} with $\lambda =k/h=\const$ Then:
\begin{itemize}
\item The FTFS scheme is stable if and only if $a\lambda\in [-1,0]$.
\item The FTBS scheme is stable if and only if $a\lambda\in [0,1]$.
Expand Down Expand Up @@ -397,7 +397,7 @@
It can also be shown that the consistency and convergence imply stability.
\end{remark}
\begin{theorem}
Consider a scheme of $J$ steps for a 1st-order-in-time linear pde of constant coefficients whose amplification factor is $g$. Let $\Phi(\theta, g)$ be the \emph{amplification polynomial}, that is the polynomial that satisfies $g$ of degree $J-1$. Then, the scheme is stable if and only if:
Consider a scheme of $J$ steps for a 1st-order-in-time linear PDE of constant coefficients whose amplification factor is $g$. Let $\Phi(\theta, g)$ be the \emph{amplification polynomial}, that is the polynomial that satisfies $g$ of degree $J-1$. Then, the scheme is stable if and only if:
\begin{itemize}
\item for any root $g_j(\theta)$ of $\Phi$ we have $\abs{g_j(\theta)}\leq 1$-
\item if $\exists \theta_0$ and $k$ such that $\abs{g_k(\theta_0)}=1$, then this root is simple.
Expand All @@ -417,9 +417,9 @@
$$
If $\abs{a\lambda}<1$, then $\abs{g_{\pm}}^2=1$ and the two roots are simple $\forall \theta\in\RR$. If $\abs{a\lambda}>1$ and $\theta=\frac{\pi}{2}$, then either $\abs{g_+}>1$ or $\abs{g_-}>1$ and the scheme is unstable. Finally, if $\abs{a\lambda}=1$ and $\theta=\frac{\pi}{2}$, then the scheme is unstable because there is a double root.
\end{proof}
\subsubsection{Second order pdes}
\subsubsection{Second order PDEs}
\begin{definition}
Consider a second order pde of the form:
Consider a second order PDE of the form:
\begin{equation}
A u_{tt}+2Bu_{tx} +Cu_{xx}+Du_t+Eu_x+F u=G
\end{equation}
Expand All @@ -431,7 +431,7 @@
u_x(f(s),g(s))=\psi(s)
\end{cases}
$$
which are tied to the \emph{compatibility condition} $h'=\phi f'+\psi g'$ that follows from the chain rule. The characteristic curves are the curves from which we cannot find the highest order derivatives of $u$ from the initial conditions and the pde. Differentiating $u_t(s)$ and $u_x(s)$ we get the system of equations for $u_{tt}$, $u_{tx}$ and $u_{xx}$:
which are tied to the \emph{compatibility condition} $h'=\phi f'+\psi g'$ that follows from the chain rule. The characteristic curves are the curves from which we cannot find the highest order derivatives of $u$ from the initial conditions and the PDE. Differentiating $u_t(s)$ and $u_x(s)$ we get the system of equations for $u_{tt}$, $u_{tx}$ and $u_{xx}$:
\begin{equation*}
\begin{cases}
A u_{tt}+2Bu_{tx} +Cu_{xx}=G-D\phi-E\psi -Fh \\
Expand All @@ -443,13 +443,13 @@
$$
A{\left(\dv{x}{t}\right)}^2-2B\dv{x}{t}+C=0
$$
The pde is called \emph{elliptic} if $AC-B^2>0$, \emph{hyperbolic} if $AC-B^2<0$ and \emph{parabolic} if $AC-B^2=0$.
The PDE is called \emph{elliptic} if $AC-B^2>0$, \emph{hyperbolic} if $AC-B^2<0$ and \emph{parabolic} if $AC-B^2=0$.
\end{definition}
\begin{definition}
Consider a finite difference scheme with $J$ steps for a 2n order homogeneous pde and $\Lambda$ be a stability region. We say that it is \emph{stable} is given $T>0$, there exists $C_T>0$ such that for any grid with $(k,h)\in \Lambda$ and for any initial values $\vf{v}_m^j$, $m\in\ZZ$, $j=0,\ldots,J-1$ we have $$\sum_{m\in\ZZ}\norm{\vf{v}_m^n}^2\leq (1+n^2)C_T\sum_{j=0}^{J-1}\sum_{m\in\ZZ}\norm{\vf{v}_m^j}^2$$ for all $n\in\NN$ such that $0\leq nk\leq T$.
Consider a finite difference scheme with $J$ steps for a 2n order homogeneous PDE and $\Lambda$ be a stability region. We say that it is \emph{stable} is given $T>0$, there exists $C_T>0$ such that for any grid with $(k,h)\in \Lambda$ and for any initial values $\vf{v}_m^j$, $m\in\ZZ$, $j=0,\ldots,J-1$ we have $$\sum_{m\in\ZZ}\norm{\vf{v}_m^n}^2\leq (1+n^2)C_T\sum_{j=0}^{J-1}\sum_{m\in\ZZ}\norm{\vf{v}_m^j}^2$$ for all $n\in\NN$ such that $0\leq nk\leq T$.
\end{definition}
\begin{theorem}
Consider a finite difference scheme with $J$ steps for a 2n order homogeneous pde whose amplification factor is $g$ and $\Phi(\theta, g)$ is the amplification polynomial. Then, the scheme is stable if and only if:
Consider a finite difference scheme with $J$ steps for a 2n order homogeneous PDE whose amplification factor is $g$ and $\Phi(\theta, g)$ is the amplification polynomial. Then, the scheme is stable if and only if:
\begin{itemize}
\item for any root $g_j(\theta)$ of $\Phi$ we have $\abs{g_j(\theta)}\leq 1$.
\item if $\exists \theta_0$ and $k$ such that $\abs{g_k(\theta_0)}=1$ then this root is at most double.
Expand Down Expand Up @@ -479,7 +479,7 @@
\end{proposition}
\subsubsection{Elliptic equations}
\begin{definition}
Let $Pu=f$ be an elliptic pde on $\Omega$. We define the following boundary conditions on $\Fr{\Omega}$:
Let $Pu=f$ be an elliptic PDE on $\Omega$. We define the following boundary conditions on $\Fr{\Omega}$:
\begin{enumerate}
\item \emph{Dirichlet}: $u=f$
\item \emph{Neumann}: $\pdv{u}{n}=g$
Expand Down
Loading

0 comments on commit 528772a

Please sign in to comment.