Derive $u(x,t)$ como solución al problema inicial/de valores límite.

Question

Dado $g : [0,\infty) \to \mathbb{R}$ con $g(0)=0$ deduzca la fórmula $$u(x,t)=\frac{x}{\sqrt{4\pi}}\int_0^t \frac 1{(t-s)^{3/2}}e^{-\frac{x^2}{4(t-s)}}g(s)\,ds$$ para una solución del problema inicial/de valor límite $$\begin{cases}u_t-u_{xx}=0 & \text{in }\mathbb{R}_+ \times (0,\infty) \\ \qquad \quad \, u=0 & \text{on } \mathbb{R}_+ \times \{t=0\} \\ \qquad \quad \, u= g & \text{on }\{x=0\} \times [0,\infty). \end{cases}$$ (Sugerencia $v(x,t):=u(x,t)-g(t)$ y ampliar $v$ a $\{x<0\}$ por reflexión impar.)

Este es el ejercicio 15 del capítulo 2 de PDE Evans, 2ª edición (páginas 87-88).

Siguiendo la pista y tomando la reflexión de impar, tenemos $$v(t)=\begin{cases}u(x,t)-g(t) & \text{if }x > 0 \\ -[u(-x,t)-g(t)] & \text{if }x<0 \end{cases}$$ y así obtenemos $$v_t=\begin{cases}u_t(x,t)-g'(t) & \text{if }x > 0 \\ -[u_t(-x,t)-g(t)] & \text{if }x<0 \end{cases} \quad \text{ and } \quad v(t)=\begin{cases} u_{xx}(x,t) & \text{if }x > 0 \\ -u_{xx}(-x,t) & \text{if }x<0. \end{cases}$$ Así, $\require{cancel} v(x,0)=\cancelto{0}{u(x,0)}-\cancelto{0}{g(0)}=0$ y $\require{cancel} v(0,t)=\cancelto{g'(t)}{u(0,t)}-g'(t)=0$ . Nuestra PDE es ahora $$\begin{cases}v_t - v_{xx} = \begin{cases} -g'(t) & \text{if } x > 0 \\ g'(t) & \text{if }x < 0 \end{cases} & \text{in } \mathbb{R}_+ \times (0,\infty) \\ \, \, \, v(x,0)=0 & \text{on } \mathbb{R}_+ \times \{t=0\} \\ \quad v(0,t) = 0& \text{on } \{x=0\} \times [0,\infty)\end{cases}$$

---

Ahora bien, la página 49 de PDE Evans, 2ª edición afirma que, para $x \in \mathbb{R}^n$ , $t > 0$ , $$u(x,t)=\int_0^t \frac 1{4\pi(t-s))^{n/2}} \int_{\mathbb{R}^n} e^{-\frac{|x-y|^2}{4(t-s)}} f(y,s) \, dy \, ds \tag{13}$$ es una solución del problema de valor inicial no homogéneo $$\begin{cases}u_t - \Delta u = f & \text{in } \mathbb{R}^n \times (0,\infty) \\ \qquad \quad u = 0 & \text{on } \mathbb{R}^n \times \{t=0\}. \tag{12} \end{cases}$$ (Las cifras $\text{(13)}$ y $\text{(12)}$ son de la página 49 del libro de texto. Además, estos $u$ 's para $\text{(13)}$ y $\text{(12)}$ son diferentes de los $u$ en el problema del ejercicio).

---

Aplicación de $\text{(13)}$ a nuestro problema de valor inicial/límite dado con $n=1$ obtenemos $$\begin{align} v(x,t)&= \int_0^t \frac 1{(4\pi(t-s))^{1/2}} \left[\int_{-\infty}^0 e^{-\frac{|x-y|^2}{4(t-s)}}g'(y,s)\,dy - \int_0^{\infty} e^{-\frac{|x-y|^2}{4(t-s)}}g'(y,s)\,dy \right]ds \end{align}$$

Por el lema de la página 46, $\int_{-\infty}^{\infty} \frac 1{(4\pi(t-s))^{1/2}} e^{-\frac{|x-y|^2}{4(t-s)}} \, dy=1$ . Así que tenemos $$\begin{align} g(t)&=\int_0^t g(s) \, ds \int_{-\infty}^{\infty} \frac 1{(4\pi(t-s))^{1/2}} e^{-\frac{|x-y|^2}{4(t-s)}} \, dy \\ &= \int_0^t \frac 1{(4\pi(t-s))^{1/2}} \left[ \int_{-\infty}^0 e^{-\frac{|x-y|^2}{4(t-s)}} \, dy \, ds + \int_0^{\infty} e^{-\frac{|x-y|^2}{4(t-s)}} \, dy \, ds \right]. \end{align}$$

Así, $$\begin{align} u(x,t)&= v(x,t)+g(t)\\ &= 2 \int_0^t \frac 1{(4\pi(t-s))^{1/2}} \int_{-\infty}^0 e^{-\frac{|x-y|^2}{4(t-s)}} \, dy \, g'(s) \, ds \end{align}$$

Pero parece que no puedo derivar en este punto la fórmula deseada de $$u(x,t)=\frac{x}{\sqrt{4\pi}}\int_0^t \frac 1{(t-s)^{3/2}}e^{-\frac{x^2}{4(t-s)}}g(s)\,ds.$$ Creo que en este punto habría que recurrir a la integración por partes. ¿Voy por buen camino?

Answer 1

Voy a publicar la respuesta de math24.files.wordpress.com/2013/02/evans_solutions-ch2.pdf. Todo el mérito es de Joe Benson, Denis Bashkirov, Minsu Kim, Helen Li y Alex Csar.

Definimos $$\begin{align*} v(x,t)\equiv\begin{cases} u(x,t)-g(t)&x>0\\ -u(-x,t)+g(t)&x\leq 0. \end{cases} \end{align*}$$ Así $$\begin{align*} v_t(x,t)= \begin{cases} u_t(x,t)-g'(t)&x>0\\ -u_t(-x,t)+g'(t)&x\leq 0 \end{cases} \end{align*}$$ y $$\begin{align*} v_{xx}(x,t)= \begin{cases} u_{xx}(x,t)&x>0\\ -u_{xx}(-x,t)&x\leq 0. \end{cases} \end{align*}$$ Por lo tanto, $$\begin{align*} \begin{cases} v_t(x,t)-v_{xx}(x,t)=\begin{cases}-g'(t)&x>0\\g'(t)&x\leq 0 \end{cases} \\ v(x,0)=0\\\\N-. v(0,t)=0. \end{casos} \end{align*}$$ Utilizando la fórmula de la solución del problema no homogéneo, obtenemos $$v(x,t)=\int_0^t\frac{1}{\sqrt{4\pi (t-s)}} \left(\int_{-\infty}^0e^{\frac{-(y-x)^2}{4(t-s)}}g'(s)\,dy-\int_0^\infty e^{\frac{-(y-x)^2}{4(t-s)}}g'(s)\,dy \right) \,ds.$$ La integral de la solución fundamental satisface $$\int_{-\infty}^\infty \frac{1}{\sqrt{4\pi(t-s)}}e^\frac{-(y-x)^2}{4(t-s)}\,dy=1,$$ y así, para $x>0$ , dejando $y-x\equiv-z$ conduce a $$\begin{align*} u(x,t)&=v(x,t)+g(t)\\ &=v(x,t)+\int_0^tg'(s)\,ds\int_{-\infty}^\infty \frac{1}{\sqrt{4\pi(t-s)}}e^\frac{-(y-x)^2}{4(t-s)}\,dy\\ &=2\int_0^t \frac{1}{\sqrt{4\pi}}(t-s)^{-\frac12}\int_{-\infty}^0e^\frac{-(y-x)^2}{4(t-s)}\,dy\, g'(s)\,ds\\ &=\int_0^t \frac{1}{\sqrt{\pi}}(t-s)^{-\frac12}\int_x^\infty e^\frac{-z^2}{4(t-s)}\,dzdg(s). \end{align*}$$ Integrando por partes, obtenemos $$\begin{align*} u(x,t)&=\frac{1}{\sqrt{\pi}}(t-s)^{-\frac12}\int_x^\infty e^\frac{-z^2}{4(t-s)}\,dz\,g(s)|_{s=0}^{s=t}\\ &\hspace{10mm}-\int_0^tg(s)\frac{1}{2\sqrt{\pi}}(t-s)^{-\frac32}\,ds\int_X^\infty e^\frac{-z^2}{4(t-s)}\,dz\\ &\hspace{20mm}-\int_0^tg(s)\frac{1}{\pi}(t-s)^{-\frac12}\,ds\int_x^\infty e^\frac{-z^2}{4(t-s)}\frac{-z^2}{4(t-s)^2}\,dz\\ &\equiv I_1-\int_0^tg(s)\frac{1}{2\sqrt{\pi}}(t-s)^{-\frac32}\,ds\int_X^\infty e^\frac{-z^2}{4(t-s)}\,dz\\ &\hspace{10mm}-\int_0^tg(s)\frac{1}{\pi}(t-s)^{-\frac12}\,ds\int_x^\infty \frac{-z}{2(t-s)}\,de^\frac{-z^2}{4(t-s)}\\ &=I_1-\int_0^tg(s)\frac{1}{2\sqrt{\pi}}(t-s)^{-\frac32}\,ds\int_X^\infty e^\frac{-z^2}{4(t-s)}\,dz\\ &\hspace{10mm}+\int_0^tg(s)\frac{1}{\sqrt{4\pi}}(t-s)^{-\frac32}\,ds(-z)e^\frac{-z^2}{4(t-s)}|_{z=x}^{z=\infty}\\ &\hspace{20mm}+\int_0^tg(s)\frac{1}{2\sqrt{\pi}}(t-s)^{-\frac32}\,ds\int_x^\infty e^\frac{-z^2}{4(t-s)}\,dz\\ &=I_1+\frac{x}{\sqrt{4\pi}}\int_0^t\frac{1}{(t-s)^\frac32}e^\frac{-x^2}{4(t-s)}g(s)\,ds. \end{align*}$$ Ahora nos centramos en $I_1$ y establece $w^2\equiv\frac{z^2}{4\epsilon}$ , $$\begin{align*} I_1&=\lim_{\epsilon\to 0^+}\frac{1}{\sqrt{\pi}}\epsilon^{-\frac12}\int_x^\infty e^\frac{-z^2}{4\epsilon}\,dz g(t-\epsilon)\\ &=g(t)\lim_{\epsilon\to 0^+}\frac{1}{\sqrt{\pi}}\int_{\frac{x^2}{4\epsilon}}^\infty 2e^{-w^2}\,dw=0. \end{align*}$$ Así demostramos $$u(x,t)=\frac{x}{\sqrt{4\pi}}\int_0^t\frac{1}{(t-s)^\frac32}e^\frac{-x^2}{4(t-s)}g(s)\,ds$$ para $x>0$ . A continuación, debemos demostrar que $$\begin{align*} \lim_{x\to 0^+}u(x,t)&=\lim_{x\to 0^+}\frac{x}{\sqrt{4\pi}}\int_{t-\delta}^t\frac{1}{(t-s)^\frac32}e^\frac{-x^2}{4(t-s)}g(s)\,ds\\ &\hspace{10mm}+\lim_{x\to 0^+}\frac{x}{\sqrt{4\pi}}\int_0^{t-\delta}\frac{1}{(t-s)^\frac32}e^\frac{-x^2}{4(t-s)}g(s)\,ds\\ &=g(t)\lim_{x\to 0^+}\frac{x}{\sqrt{4\pi}}\int_{t-\delta}^t\frac{1}{(t-s)^\frac32}e^\frac{-x^2}{4(t-s)}\,ds\\ &=g(t)\lim_{x\to 0^+}\frac{x}{\sqrt{4\pi}}\int_0^\delta \frac{1}{s^\frac32}e^\frac{-x^2}{4s\,ds}. \end{align*}$$ Para fijos $x$ , dejemos que $s\equiv\frac{x^2}{w^2}$ y obtener $$\begin{align*} \lim_{x\to 0^+} u(x,t)&=g(t)\lim_{x\to 0^+}\frac{x}{2\sqrt{\pi}}\int_{-\infty}^{x^2/\delta}\frac{w^2}{x^3}e^\frac{-w^2}{4}\frac{-2x^2}{w^3}\,dw\\ &=g(t)\lim_{x\to 0^+}\frac{1}{\sqrt{\pi}}\int_{x^2/\delta}^\infty e^\frac{-w^2}{4}\,dw\\ &=g(t)\frac{1}{\sqrt{\pi}}\int_0^\infty e^\frac{-w^2}{4}\,dw\\ &=g(t) \end{align*}$$ y hemos terminado.