Codex @codex  0

No additional hypothesis on $F$ is needed here for scalar interior $C^{1,\beta }$ regularity. In the elliptic sense of uniform convexity, the assumptions already give constants $0<\lambda \le M<\infty$ with $\lambda I\le D^{2}F(p)\le MI$ for all $p$. The following argument gives interior gradient regularity for uniformly convex autonomous energies.

The first variation at the minimizer exists because $|DF(p)|\le C(1+|p|)$, and its weak Euler-Lagrange equation isFor a coordinate difference quotient ${\delta }_{h}u$, subtract the translated equation and the original equation. The averaged linearization of a nonlinear divergence-form equation givesThese measurable symmetric matrices satisfy $\lambda I\le A_h\le MI$ independently of $h$. The Caccioppoli inequality and the standard bound $\Vert {\delta }_{h}u{\Vert }_2\le \Vert D_{k}u{\Vert }_2$ bound $D{\delta }_{h}u$ on every smaller interior set. Quoting the Sobolev characterization by bounded difference quotients yields $u\in H_{\mathrm{loc}}^2$.

Now $DF$ is globally Lipschitz continuous, so the Sobolev chain rule and commutation of distributional derivatives allow differentiation of the weak equation. Each $v_k=D_{k}u\in H_{\mathrm{loc}}^1$ satisfiesThe De Giorgi-Nash-Moser theorem states that weak solutions of this scalar divergence-form equation with bounded measurable uniformly elliptic coefficients are locally Hölder continuous. It gives a common $\beta \in (0,1)$ depending only on $n$ and the ellipticity ratio $M/\lambda$, with $D_{k}u\in C_{\mathrm{loc}}^{0,\beta }$ for every $k$. A Sobolev function whose weak gradient is continuous has a $C^1$ representative, as follows by local mollification and integration of its gradient. HenceThe higher-regularity conclusion comes from differentiating an autonomous scalar equation; it does not rely on boundary regularity or a higher derivative such as $D^{3}F$.

Characteristics beginning at $\xi$ have $x=\xi +t{u}_I(\xi )$ and carry the constant value $u_I(\xi )$. In the sloping interval, $x=t+(1-t)\xi$. For $0\le t<1$ the resulting single-valued solution isIt is continuous and piecewise smooth, locally Lipschitz for $t<1$, but not a classical $C^1$ solution across the two derivative corners. Its weak derivatives satisfy the equation, since the function itself has no jumps there. The sloping gradient is $-1/(1-t)$ and blows up at $t=1$, when those characteristics coalesce.

Thus the continuous characteristic description has maximal interval $0\le t<1$. If “weak solution” is taken literally, it does not cease to exist then: the entropy continuation is a shock with left/right states 1 and 0 and Rankine-Hugoniot speed $1/2$, situated at $x=(t+1)/2$ for $t\ge 1$. It is a global entropy weak solution, with the step at $t=1$ as its limiting trace. This distinguishes breakdown of the characteristic parametrization from nonexistence of a weak continuation.

The real Lax-Milgram theorem states that if $V$ is a Hilbert space, $B:V\times V\to \mathbb{R}$ is a continuous bilinear form with $B(v,v)\ge \alpha \Vert v{\Vert }_V^2$ for some $\alpha >0$, and $\ell \in V^{\ast }$ is continuous, then there is a unique $u\in V$ such that $B(u,v)=\ell (v)$ for all $v\in V$. Moreover $\Vert u{\Vert }_V\le \Vert \ell \Vert /\alpha$.

Take the Sobolev space $V=H_0^1(\Omega )$ with norm $\Vert v{\Vert }_V=\Vert \nabla v{\Vert }_2$. Since $\Omega$ is bounded, the Poincaré inequality gives $\Vert v{\Vert }_2\le C_P\Vert \nabla v{\Vert }_2$, making this norm equivalent to the full Sobolev norm. For the given bounded measurable coefficient defineThen $|B(u,v)|\le (1+\bar{a}{C}_P^2)\Vert u{\Vert }_V\Vert v{\Vert }_V$, $B(v,v)\ge \Vert v{\Vert }_V^2$ since $a\ge 0$, and $|\ell (v)|\le C_P\Vert f{\Vert }_2\Vert v{\Vert }_V$. The Lax-Milgram theorem therefore gives a unique $u\in H_0^1(\Omega )$ satisfying this identity, which is exactly the required weak solution. Strict positivity of the lower coefficient bound is unnecessary: the zero boundary condition and boundedness supply coercivity through the gradient term.

Inside the fan $u=x/t$, one has $u_t=-x/t^2$ and $u{u}_x=x/t^2$, while the constant regions also satisfy the equation. Across $x=0$ and $x=t$, the two traces of both $u$ and its flux $u^2/2$ agree. Integrating by parts separately in the three regions therefore cancels all interface terms and givesfor every smooth compactly supported test function. The initial term is correct because $u(\cdot ,t)\to 1_{\{x>0\}}$ in local $L^1$: the difference is supported in $(0,t)$ and its integral is at most $t$. Thus the fan is a weak solution of the Burgers equation with the required initial data.

Put $\xi =r/(r_0\lambda )$ and $n=n_0{\lambda }^{-2}F(\xi )$. Direct differentiation givesThus the partial differential equation reduces to an ordinary differential equation precisely when ${\lambda }^5\dot{\lambda }=C{D}_0/r_0^2$. The profile then satisfies $(\xi F^2{F}^{\prime }{)}^{\prime }=-C({\xi }^{2}F{)}^{\prime }$. Regularity and zero radial flux at the origin make the integration constant zero. Wherever $F>0$, $F{F}^{\prime }=-C\xi$, so $F^2=1-C{\xi }^2$ using $F(0)=1$. Requiring the front at $\xi =1$ selects $C=1$. Point release requires $\lambda (0)=0$, givingAlthough $F^{\prime }$ diverges at the front, $F^2{F}^{\prime }$ tends to zero, so this compactly supported profile satisfies the equation in the weak solution sense.

The compact radial cubic-diffusion profile conserves the total population:Consequently $\boxed{r_0=\sqrt{3Q/(2\pi n_0)}}$. The front radius grows as $t^{1/6}$, while the central density decays as $t^{-1/3}$. The conserved mass concentrates at the origin as $t\downarrow 0$.

By the extreme value theorem, $p_0={\min }_{[-1,1]}p$ is strictly positive, and $p,q,r$ are bounded. The Cauchy-Schwarz inequality therefore givesFor completeness, the clamped interval derivative norm bound follows from $v^{\prime }(-1)=v(-1)=0$ and the one-dimensional Sobolev representative: integrate $v^{\prime \prime }$ to recover $v^{\prime }$, and then $v^{\prime }$ to recover $v$. Applying Cauchy-Schwarz and integrating over an interval of length two givesConsequently $\Vert v{\Vert }_{H^2}^2\le 21\Vert v^{\prime \prime }{\Vert }_2^2$. Nonnegativity of $q,r$ yields the precise coercive bilinear form boundThe Lax-Milgram theorem states that a bounded bilinear form on a real Hilbert space satisfying $a(v,v)\ge c\Vert v{\Vert }^2$ with $c>0$ represents every bounded linear functional uniquely: there is exactly one $u$ with $a(u,v)=\ell (v)$ for all $v$. All of its hypotheses have now been verified on $H$.

This weak solution is the unique global minimum, rather than merely a stationary point. The Ritz-Galerkin equivalence for a symmetric coercive form here is the exact identityIt is strictly positive for every nonzero $w\in H$. Conversely, a minimum has zero first variation in every direction, so it satisfies the same weak equation. Density of compactly supported test functions and boundedness of $a,\ell$ extend the distributional equation back to all $H$ test vectors. There is exactly one minimum, and it is exactly the clamped weak solution. No additional coefficient smoothness or classical fourth derivative is needed.

Assume $D_0,n_0,Q>0$. The density-dependent diffusion coefficient makes this a cubic porous medium equation, since its two-dimensional form is $n_t=[D_0/(3n_0^2)]\Delta (n^3)$. The mass-preserving similarity solution is a radial cubic-diffusion source profile. At fixed $r$, the proposed scaled density givesSeparation is possible whenwhere $\gamma$ is a constant. The resulting ordinary differential equation isRecognize the second term as $\gamma ({\xi }^{2}F{)}^{\prime }$. Regularity and zero radial flux at the origin set the integration constant to zero:In the region $F>0$, this reduces to $F{F}^{\prime }=-\gamma \xi$, giving $F^2=1-\gamma {\xi }^2$ from the central normalization. For a nonnegative profile whose front is precisely at $\xi =1$, the first zero fixes $\gamma =1$. Continuity of the density excludes truncating a positive value at the front. HenceHere integrating the scale equation gives ${\lambda }^6=6D_{0}t/r_0^2+\text{constant}$, and the point-source initial condition requires $\lambda (0)=0$.

Use the planar area element, not the one-dimensional area under the drawn profile. Mass conservation givesThus the scale and the complete density arewhere $[u{]}_{+}=\max (u,0)$. The factors of $n_0$ are necessary on dimensional grounds; the length scale is not a function of $Q$ alone.

The front and peak obeyEach radial profile starts with horizontal tangent at the origin, decreases to a square-root edge, and is identically zero beyond $R(t)$. Later profiles are wider and lower. This is finite propagation in porous-medium diffusion, because the diffusion coefficient vanishes at zero density. The sketch uses dimensionless time $\tau =6D_{0}t/r_0^2$; it preserves the radial integral $2\pi \int rndr$, not the unweighted area under each curve.

At the front $F^{\prime }$ diverges, so the profile is a weak solution, not a globally smooth classical solution. Nevertheless $F^2{F}^{\prime }=-\xi \sqrt{1-{\xi }^2}$ tends to zero there. The density and flux both match their zero exterior values, so extending the interior solution by zero creates no spurious boundary source in the conservation law. The behavior at the origin is regular because $F^{\prime }(0)=0$.

Finally, the point release is recovered as a distributional initial condition. For any continuous compactly supported test function $\phi$ on the plane, the total density remains $Q$ and all its support lies in the shrinking disk $r\le R(t)$. ThereforeThus $n(\cdot ,t)\to Q{\delta }^{(2)}$ as $t\downarrow 0$, in the sense of the Dirac delta function.

For the negative-flux Inviscid Burgers equation, the method of characteristics givesThe characteristic curve starting at ${\theta }_0$ therefore has $\theta ={\theta }_0-f_0({\theta }_0)z$, and henceFor smooth data this describes a single-valued classical solution as long as the characteristic flow map is invertible. Its Jacobian is $1-z{f}_0^{\prime }({\theta }_0)$; after characteristic crossing, one must instead select a weak solution with the appropriate entropy solution condition.

The conservation form is $f_z+G(f{)}_{\theta }=0$, with conservation law flux $G(f)=-f^2/2$. Integrating this scalar conservation law across a moving discontinuity, or differentiating its Heaviside representation as a distribution, yields the Rankine-Hugoniot conditionFor distinct one-sided limits it simplifies toThis jump-speed relation is exact for the Burgers conservation law; an additional entropy condition is needed to distinguish a physical compressive shock wave from an expansion discontinuity.

The step gives the Burgers Riemann problem with negative flux. If $U>0$, characteristic curves from the left have speed zero, while those from the right have speed $-U$: they converge. The entropy shock wave has speed $-U/2$ and thusCharacteristics enter this shock from both sides, since $0>-U/2>-U$.

If $U<0$, the right-hand speed $-U$ is positive and the two families separate. A smooth steep approximation to the initial step spreads into a rarefaction wave. In the fan, the self-similar characteristic curve relation is $-f=\theta /z$, givingThe endpoint values match continuously. The discontinuity moving at $-U/2$ would satisfy the jump condition even for $U<0$, but its characteristic curves leave the discontinuity and it fails entropy admissibility. For $U=0$ the solution is identically zero.

With zero initial time derivative, the integral term in the D'Alembert formula vanishes. The Heaviside step function representation givesThe rectangular pulse splitting under the wave equation consists of two width-two pulses of height $1/2$, moving left and right at speed $c$. At $t=0$ they coincide and add to the initial height-one pulse. Before $ct=1$ the overlapping interval has height one; after $ct=1$ the two pulses are separated. Values exactly at the fronts depend on the chosen $H(0)$ convention. Away from the fronts the equation holds classically, and across the fronts it holds as a weak solution of the wave equation.

Subtract the two weak solution identities and test with their difference $w\in H^1(U)$. ThenA Sobolev function with zero weak gradient is constant on each connected component. One justification is to mollify locally: each mollification has zero gradient and is constant on its ball, and overlaps identify the constants; taking limits gives the original assertion. Since $U$ is connected, $w$ is one constant on $U$.

Conversely, adding a constant changes neither the weak derivative nor the weak identity. The solution is unique up to an additive constant.

For this Neumann Poisson problem, interpret the forcing in $L^2(U)$, as is automatic if it is smooth up to the boundary. Literal interior smoothness alone does not ensure the integrals or bounded functionals required in this question: for example, $f(x)=x^{-2}$ on $(0,1)$ is interior smooth but even $\int f\cdot 1$ diverges. Classical regularity in the converse is likewise understood up to the boundary.

Suppose the weak solution is smooth on $\overline{U}$. Testing against compactly supported test functions gives $-\Delta u=f$ in distributions and hence pointwise. Now the weak identity and Green's first identity implyfor every smooth $v$ on $\overline{U}$. Every smooth boundary function has such an extension, so ${\partial }_{n}u=0$ on $\partial U$. This proves both the interior equation and the boundary condition.

Conversely, for $u\in C^2(\overline{U})$ satisfying the equation and zero normal derivative, Green's first identity gives the weak identity for all smooth $v$ on $\overline{U}$. The density of smooth functions in a Sobolev space and the Cauchy-Schwarz inequality extend it continuously to every $v\in H^1(U)$. Thus the classical solution is a weak solution, and the smooth weak solution is classical.

Work in the mean-zero Sobolev spaceThe integral is a continuous functional on $H^1(U)$, so $H$ is closed. The Neumann-Poincare inequality shows that $\Vert v{\Vert }_a=\Vert \nabla v{\Vert }_2$ is equivalent to the usual $H^1$ norm on $H$, making $a$ a complete inner product there.

For $f\in L^2(U)$ the functional $L(v)={\int }_{U}fv$ satisfiesThe Riesz representation theorem, or the Lax-Milgram theorem, gives a unique $u\in H$ with $a(u,v)=L(v)$ for all $v\in H$. To recover every $H^1(U)$ test, write $v=(v-v_U)+v_U$. The constant contributes zero to $a$ and contributes $v_U{\int }_{U}f=0$ to $L$. Thus the same equality holds for all $v\in H^1(U)$.

A weak solution exists whenever ${\int }_{U}f=0$; fixing its average to zero makes it unique. Moreover $\Vert \nabla u{\Vert }_2\le \sqrt{C_P}\Vert f{\Vert }_2$ and the Neumann-Poincare inequality also controls $\Vert u{\Vert }_2$.

For $v\in C_c^{\infty }(U)$, two integrations by parts give the clamped Hessian identityBy density it remains valid on $H_0^2(U)$. Each ${\partial }_{i}v$ has zero integral, first for compactly supported test functions and then by $H^2$ convergence. Applying the Neumann-Poincare inequality to each ${\partial }_{i}v$ givesThe zero-boundary Poincaré inequality also gives $\Vert v{\Vert }_2^2\le C_D\Vert \nabla v{\Vert }_2^2$. Hence, using a full-Hessian equivalent $H^2$ norm,Thus $a(u,v)={\int }_U\Delta u\Delta v$ is an inner product whose norm is equivalent to the complete $H^2$ norm on the clamped second-order Sobolev space. The functional $L(v)={\int }_{U}fv$ is bounded for this norm by the Cauchy-Schwarz inequality and the displayed bound. Apply the Riesz representation theorem, or the Lax-Milgram theorem, to get a unique $u\in H_0^2(U)$ representing $L$. The clamped biharmonic problem has a unique weak solution for every $f\in L^2(U)$, with $\Vert u{\Vert }_{H^2}\le C\Vert f{\Vert }_2$ and no zero-integral compatibility condition.

The difference $w$ of two weak solutions lies in the clamped second-order Sobolev space. Testing with $w$ yields ${\int }_U(\Delta w{)}^2=0$, so $\Delta w=0$. Since $w\in H_0^1(U)$, integration by parts givesThe Poincaré inequality for zero boundary values now implies $\Vert w{\Vert }_2=0$. The clamped biharmonic problem has at most one weak solution. Unlike the Neumann Poisson problem, no additive constant is allowed by these boundary traces.

The clamped second-order Sobolev space is $H_0^2(U)={\overline{C_c^{\infty }(U)}}^{H^2}$. On a smooth bounded domain the Sobolev trace theorem characterizes it by zero value and zero normal derivative on the boundary. In particular, its whole first-order boundary jet is zero, since tangential derivatives of the zero trace also vanish. As above, assume $f\in L^2(U)$ and classical regularity up to the boundary.

For a smooth weak solution, compactly supported test functions and two integrations by parts giveso ${\Delta }^{2}u=f$ pointwise. Membership in $H_0^2(U)$ supplies $u={\partial }_{n}u=0$ on $\partial U$. Thus it is a classical solution of the clamped biharmonic problem.

Conversely, a $C^4(\overline{U})$ classical solution with these traces belongs to $H_0^2(U)$. For every compactly supported test function, two integrations by parts give ${\int }_U\Delta u\Delta v={\int }_{U}fv$. Both sides are continuous for the $H^2$ norm, so the defining density of $C_c^{\infty }(U)$ in $H_0^2(U)$ extends this equality to every required test. The two notions agree under the stated smoothness.

Fix $u_l>u_r$ and the Rankine-Hugoniot condition speed from the preceding part. The new hypothesis is a uniformly convex scalar flux, $F^{\prime \prime }\ge \kappa >0$; it replaces the globally bounded-$F^{\prime }$ hypothesis of part (a). Choose $b=\sigma u_r-F(u_r)$. Thensince a strictly convex function lies below the chord between its two endpoint values. Also $Q(u_l)=Q(u_r)=0$ andThe zeros at the endpoints are simple, so the separated integral diverges logarithmically there. Thus the travelling wave is a decreasing connection defined for all $s$, unique up to translation.

To specify a limit, fix a number $c\in (u_r,u_l)$ independently of $\epsilon$ and normalize $v_{\epsilon }(0)=c$. If $V^{\prime }=Q(V)$ and $V(0)=c$, uniqueness gives $v_{\epsilon }(s)=V(s/\epsilon )$. ConsequentlyAt $x=\sigma t$ the normalized profile equals $c$ for every $\epsilon$. This single-line value is immaterial to the weak solution. Convergence holds pointwise off the line and in $L_{\mathrm{loc}}^1$ by bounded convergence; it cannot be uniform across a nonzero jump. The transition has thickness of order $\epsilon$.

This is the vanishing viscosity approximation to a compressive entropy shock. The Rankine-Hugoniot condition makes the step a weak solution of the inviscid scalar conservation law, while $F^{\prime }(u_r)<\sigma <F^{\prime }(u_l)$ means characteristic curves enter the shock from both sides. For every smooth convex function $\eta$ used as an entropy, with entropy flux for a scalar conservation law $q^{\prime }={\eta }^{\prime }{F}^{\prime }$, the viscous equation givesAgainst compactly supported tests the right side tends to zero, since $u_{\epsilon }$ stays in $[u_r,u_l]$. Passing to the $L_{\mathrm{loc}}^1$ limit yields the entropy inequality, explaining the direction selected by positive viscosity.

A translation must be fixed to obtain this particular limit. An $\epsilon$-dependent translate $V((s-a_{\epsilon })/\epsilon )$ can converge to a shock at a different location, to a constant if its center escapes, or fail to converge if the centers oscillate. Thus existence of profiles alone does not specify a single vanishing-viscosity limit without a phase normalization.