Documentation

Atlas.IntroductionToPartialDifferentialEquations.code.CM19.TransportBurgers

def TransportBurgers.SolvesTransport {n : } (X : (Fin (n + 1))Fin (n + 1)) (u : (Fin (n + 1))) :

The transport equation $\sum_\mu X^\mu \partial_\mu u = 0$ associated to a vector field $X$ on $\mathbb{R}^{n+1}$. A function $u$ solves the transport equation at every point $p$ iff the directional derivative of $u$ along $X(p)$ vanishes.

Instances For
    def TransportBurgers.IsIntegralCurve {n : } (X : (Fin (n + 1))Fin (n + 1)) (γ : Fin (n + 1)) :

    A curve $\gamma : \mathbb{R} \to \mathbb{R}^{n+1}$ is an integral curve of the vector field $X$ if $\gamma'(s) = X(\gamma(s))$ for all $s$.

    Instances For
      theorem TransportBurgers.transport_constant_along_chars {n : } (X : (Fin (n + 1))Fin (n + 1)) (u : (Fin (n + 1))) (γ : Fin (n + 1)) (hu : SolvesTransport X u) ( : IsIntegralCurve X γ) (s : ) (hu_diff : DifferentiableAt u (γ s)) (hγ_diff : DifferentiableAt γ s) :
      deriv (fun (s : ) => u (γ s)) s = 0

      Proposition 1.0.1 (Connection between transport equations and ODEs). If $u$ solves the transport equation $\sum_\mu X^\mu \partial_\mu u = 0$ and $\gamma$ is an integral curve of $X$, then $u$ is constant along $\gamma$, i.e. $\frac{d}{ds} u(\gamma(s)) = 0$.

      Burger's equation $u_t + u\,u_x = 0$ in $1 + 1$ dimensions. A function $u(t, x)$ solves Burger's equation iff the time and space partial derivatives satisfy this nonlinear conservation law at every spacetime point.

      Instances For
        theorem TransportBurgers.c1_implies_t_diff (u : ) (hC1 : ContDiff 1 (Function.uncurry u)) (t x : ) :
        DifferentiableAt (fun (t' : ) => u t' x) t

        Helper: if the uncurried form of $u : \mathbb{R} \to \mathbb{R} \to \mathbb{R}$ is $C^1$, then for fixed $x$ the slice $t' \mapsto u(t', x)$ is differentiable at $t$.

        theorem TransportBurgers.c1_implies_x_diff (u : ) (hC1 : ContDiff 1 (Function.uncurry u)) (t x : ) :
        DifferentiableAt (fun (x' : ) => u t x') x

        Helper: if the uncurried form of $u$ is $C^1$, then for fixed $t$ the spatial slice $x' \mapsto u(t, x')$ is differentiable at $x$.

        theorem TransportBurgers.c1_continuous_x (u : ) (hC1 : ContDiff 1 (Function.uncurry u)) (t : ) :
        Continuous fun (x : ) => u t x

        Helper: if $u$ is jointly $C^1$, then for fixed $t$ the spatial slice $x \mapsto u(t, x)$ is continuous.

        theorem TransportBurgers.deriv_t_eq_fderiv_10 (u : ) (hC1 : ContDiff 1 (Function.uncurry u)) (t x : ) :
        deriv (fun (t' : ) => u t' x) t = (fderiv (Function.uncurry u) (t, x)) (1, 0)

        The partial time-derivative of $u$ equals the joint Fréchet derivative of the uncurried $u$ paired against the unit tangent vector $(1, 0)$.

        theorem TransportBurgers.deriv_t_continuous_in_x (u : ) (hC1 : ContDiff 1 (Function.uncurry u)) (t : ) :
        Continuous fun (x : ) => deriv (fun (t' : ) => u t' x) t

        For $u$ jointly $C^1$ and fixed $t$, the map $x \mapsto u_t(t, x)$ is continuous in the spatial variable.

        theorem TransportBurgers.bounded_of_continuous_tendsto_zero (f : ) (hf_cont : Continuous f) (hf_top : Filter.Tendsto f Filter.atTop (nhds 0)) (hf_bot : Filter.Tendsto f Filter.atBot (nhds 0)) :
        ∃ (M : ), 0 M ∀ (x : ), |f x| M

        A continuous function $f : \mathbb{R} \to \mathbb{R}$ that decays to $0$ at both $\pm\infty$ is globally bounded: $\exists M \ge 0$ with $|f(x)| \le M$ for all $x$.

        A continuous, integrable function on $\mathbb{R}$ that decays to $0$ at $\pm\infty$ has an integrable square: $f^2 \in L^1(\mathbb{R})$.

        theorem TransportBurgers.c1_decay_implies_leibniz_hypotheses (u : ) (T : ) (_hC1 : ContDiff 1 (Function.uncurry u)) (_hdecay : tSet.Icc 0 T, Filter.Tendsto (fun (x : ) => u t x) Filter.atTop (nhds 0) Filter.Tendsto (fun (x : ) => u t x) Filter.atBot (nhds 0)) (hInt_u : tSet.Icc 0 T, MeasureTheory.Integrable (fun (x : ) => u t x) MeasureTheory.volume) (hLeibnizBound : t'Set.Icc 0 T, ∃ (s : Set ) (_ : s nhds t') (bound : ), (∀ᵐ (x : ), t''s, deriv (fun (t''' : ) => u t''' x ^ 2) t'' bound x) MeasureTheory.Integrable bound MeasureTheory.volume) (t' : ) (_ht' : t' Set.Icc 0 T) :
        MeasureTheory.Integrable (fun (x : ) => u t' x ^ 2) MeasureTheory.volume ∃ (s : Set ) (_ : s nhds t'), (∀ᶠ (t'' : ) in nhds t', MeasureTheory.AEStronglyMeasurable (fun (x : ) => u t'' x ^ 2) MeasureTheory.volume) MeasureTheory.AEStronglyMeasurable (fun (x : ) => deriv (fun (t'' : ) => u t'' x ^ 2) t') MeasureTheory.volume ∃ (bound : ), (∀ᵐ (x : ), t''s, deriv (fun (t''' : ) => u t''' x ^ 2) t'' bound x) MeasureTheory.Integrable bound MeasureTheory.volume ∀ᵐ (x : ), t''s, HasDerivAt (fun (t''' : ) => u t''' x ^ 2) (deriv (fun (t''' : ) => u t''' x ^ 2) t'') t''

        Packages the hypotheses needed to apply a Leibniz / differentiation-under-the-integral rule to $\int_{\mathbb{R}} u(t, x)^2 \, dx$: from $C^1$ regularity, spatial decay, spatial integrability, and a uniform local bound on $\partial_t (u^2)$, derive the exact ensemble of measurability, integrability, and pointwise differentiability hypotheses required.

        theorem TransportBurgers.hasDerivAt_integral (f : ) (t : ) (hf_int : MeasureTheory.Integrable (f t) MeasureTheory.volume) {s : Set } (hs : s nhds t) (hf_meas : ∀ᶠ (t' : ) in nhds t, MeasureTheory.AEStronglyMeasurable (f t') MeasureTheory.volume) (hf_deriv_meas : MeasureTheory.AEStronglyMeasurable (fun (x : ) => deriv (fun (t' : ) => f t' x) t) MeasureTheory.volume) {bound : } (h_bound : ∀ᵐ (x : ), t's, deriv (fun (t'' : ) => f t'' x) t' bound x) (bound_int : MeasureTheory.Integrable bound MeasureTheory.volume) (h_diff_on : ∀ᵐ (x : ), t's, HasDerivAt (fun (t'' : ) => f t'' x) (deriv (fun (t'' : ) => f t'' x) t') t') :
        HasDerivAt (fun (t' : ) => (x : ), f t' x) ( (x : ), deriv (fun (t' : ) => f t' x) t) t

        Differentiation under the integral sign (specialised wrapper): under the standard dominated-convergence hypotheses, the parameter integral $t \mapsto \int_{\mathbb{R}} f(t, x) \, dx$ has derivative $\int_{\mathbb{R}} \partial_t f(t, x) \, dx$ at $t$.

        theorem TransportBurgers.l2_norm_conservation_burgers (u : ) (hu : SolvesBurgers u) (T : ) (_hT : T 0) (hC1 : ContDiff 1 (Function.uncurry u)) (hu_decay : tSet.Icc 0 T, Filter.Tendsto (fun (x : ) => u t x) Filter.atTop (nhds 0) Filter.Tendsto (fun (x : ) => u t x) Filter.atBot (nhds 0)) (hInt_u : tSet.Icc 0 T, MeasureTheory.Integrable (fun (x : ) => u t x) MeasureTheory.volume) (hLeibnizBound : t'Set.Icc 0 T, ∃ (s : Set ) (_ : s nhds t') (bound : ), (∀ᵐ (x : ), t''s, deriv (fun (t''' : ) => u t''' x ^ 2) t'' bound x) MeasureTheory.Integrable bound MeasureTheory.volume) (t : ) (ht : t Set.Icc 0 T) :
        (x : ), u t x ^ 2 = (x : ), u 0 x ^ 2

        Proposition 2.0.1 (Burger's equation is a conservation law). Let $u(t, x)$ be a $C^1$ solution of Burger's equation on $[0, T] \times \mathbb{R}$ that decays to $0$ as $x \to \pm \infty$ uniformly for $t \in [0, T]$ (with suitable integrability of $u(t, \cdot)$ and a Leibniz-rule bound). Then the spatial $L^2$ norm is conserved: $$\int_{\mathbb{R}} u(t, x)^2 \, dx = \int_{\mathbb{R}} u(0, x)^2 \, dx \quad \text{for all } t \in [0, T].$$

        structure TransportBurgers.Characteristics.BurgersCharacteristic (u : ) (γ_t γ_x : ) :

        Definition 2.0.1 (Characteristic curves). A characteristic curve for Burger's equation with solution $u$ is a pair $(\gamma_t, \gamma_x) : \mathbb{R} \to \mathbb{R}^2$ satisfying the ODE system $$\frac{d}{ds} \gamma_t(s) = 1, \qquad \frac{d}{ds} \gamma_x(s) = u(\gamma_t(s), \gamma_x(s)).$$

        • time_param (s : ) : deriv γ_t s = 1
        • space_param (s : ) : deriv γ_x s = u (γ_t s) (γ_x s)
        Instances For
          theorem TransportBurgers.Characteristics.burgers_chars_straight_lines (u : ) (γ_t γ_x : ) (_hu : SolvesBurgers u) (hchar : BurgersCharacteristic u γ_t γ_x) (hconst : ∀ (s : ), u (γ_t s) (γ_x s) = u (γ_t 0) (γ_x 0)) (_ht_init : γ_t 0 = 0) (hx_diff : Differentiable γ_x) (s : ) :
          γ_x s = γ_x 0 + u (γ_t 0) (γ_x 0) * s

          Proposition 2.0.3 (Burger characteristics are straight lines). Given a characteristic $(\gamma_t, \gamma_x)$ for Burger's equation along which $u$ is constant, the spatial component is a straight line: $\gamma_x(s) = \gamma_x(0) + u(\gamma_t(0), \gamma_x(0)) \cdot s$.

          theorem TransportBurgers.Characteristics.burgers_char_zero_accel_time (u : ) (γ_t γ_x : ) (hchar : BurgersCharacteristic u γ_t γ_x) (s : ) :
          deriv (deriv γ_t) s = 0

          The time component of a Burger's characteristic has zero acceleration: $\gamma_t''(s) = 0$ for all $s$ (a direct consequence of $\gamma_t'(s) = 1$).

          theorem TransportBurgers.Characteristics.burgers_char_zero_accel_space (u : ) (γ_t γ_x : ) (hchar : BurgersCharacteristic u γ_t γ_x) (hconst : ∀ (s : ), deriv (fun (s : ) => u (γ_t s) (γ_x s)) s = 0) (hx_diff : Differentiable γ_x) (s : ) :
          deriv (deriv γ_x) s = 0

          The spatial component of a Burger's characteristic has zero acceleration: if $u$ is constant along the characteristic, then $\gamma_x''(s) = 0$.